Search

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. [ABSTRACT FROM AUTHOR]

In past work, tsunami detection algorithms (TDAs) have been proposed, and successfully applied to offline tsunami detection, based on analyzing tsunami currents inverted from high-frequency (HF) radar Doppler spectra. With this method, however, the detection of small and short-lived tsunami currents in the most distant radar ranges is challenging due to conflicting requirements on the Doppler spectra integration time and resolution. To circumvent this issue, in Part I of this work, we proposed an alternative TDA, referred to as time correlation (TC) TDA, that does not require inverting currents, but instead detects changes in patterns of correlations of radar signal time series measured in pairs of cells located along the main directions of tsunami propagation (predicted by geometric optics theory); such correlations can be maximized when one signal is time-shifted by the pre-computed long wave propagation time. We initially validated the TC-TDA based on numerical simulations of idealized tsunamis in a simplified geometry. Here, we further develop, extend, and apply the TC algorithm to more realistic tsunami case studies. These are performed in the area West of Vancouver Island, BC, where Ocean Networks Canada recently deployed a HF radar (in Tofino, BC), to detect tsunamis from far- and near-field sources, up to a 110 km range. Two case studies are considered, both simulated using long wave models (1) a far-field seismic, and (2) a near-field landslide, tsunami. Pending the availability of radar data, a radar signal simulator is parameterized for the Tofino HF radar characteristics, in particular its signal-to-noise ratio with range, and combined with the simulated tsunami currents to produce realistic time series of backscattered radar signal from a dense grid of cells. Numerical experiments show that the arrival of a tsunami causes a clear change in radar signal correlation patterns, even at the most distant ranges beyond the continental shelf, thus making an early tsunami detection possible with the TC-TDA. Based on these results, we discuss how the new algorithm could be combined with standard methods proposed earlier, based on a Doppler analysis, to develop a new tsunami detection system based on HF radar data, that could increase warning time. This will be the object of future work, which will be based on actual, rather than simulated, radar data. [ABSTRACT FROM AUTHOR]

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 generation models are used, in which the bottom (slide) layer is depth-integrated and represented by either a dense Newtonian fluid or a granular flow, in which inter-granular stresses are governed by Coulomb friction (Savage and Hutter model). In both cases, the upper (water) layer flow is simulated with the non-hydrostatic 3D σ-layer model NHWAVE. Both models are validated by simulating laboratory experiments for SMFs made of glass beads moving down a steep plane slope. In those, we assess the convergence of results (i.e., SMF motion and surface wave generation) with model parameters and their sensitivity to slide parameters (i.e., viscosity, bottom friction, and initial submergence). The historical Currituck SMF is simulated with the viscous slide model, to estimate relevant parameters for simulating tsunami generation from a possible SMF sited near the Hudson River Canyon. Compared to a rigid slump, we find that deforming SMFs of various rheology, despite having a slightly larger initial acceleration, generate a smaller tsunami due to their spreading and thinning out during motion, which gradually makes them less tsunamigenic; the latter behavior is controlled by slide rheology. Coastal tsunami hazard is finally assessed by performing tsunami simulations with the Boussinesq long wave model FUNWAVE-TVD, initialized by SMF tsunami sources, in nested grids of increasing resolution. While initial tsunami elevations are very large (up to 25 m for the rigid slump), nearshore tsunami elevations are significantly reduced in all cases (to a maximum of 6.5 m). However, at most nearshore locations, surface elevations obtained assuming a rigid slump are up to a factor of 2 larger than those obtained for deforming slides. We conclude that modeling SMFs as rigid slumps provides a conservative estimate of coastal tsunami hazard while using a more realistic rheology, in general, reduces coastal tsunami impact. [ABSTRACT FROM AUTHOR]

This work is part of a tsunami inundation mapping activity carried out along the US East Coast since 2010, under the auspice of the National Tsunami Hazard Mitigation program (NTHMP). The US East Coast features two main estuaries with significant tidal forcing, which are bordered by numerous critical facilities (power plants, major harbors,...) as well as densely built low-level areas: Chesapeake Bay and the Hudson River Estuary (HRE). HRE is the object of this work, with specific focus on assessing tsunami hazard in Manhattan, the Hudson and East River areas. In the NTHMP work, inundation maps are computed as envelopes of maximum surface elevation along the coast and inland, by simulating the impact of selected probable maximum tsunamis (PMT) in the Atlantic ocean margin and basin. At present, such simulations assume a static reference level near shore equal to the local mean high water (MHW) level. Here, instead we simulate maximum inundation in the HRE resulting from dynamic interactions between the incident PMTs and a tide, which is calibrated to achieve MHW at its maximum level. To identify conditions leading to maximum tsunami inundation, each PMT is simulated for four different phases of the tide and results are compared to those obtained for a static reference level. We first separately simulate the tide and the three PMTs that were found to be most significant for the HRE. These are caused by: (1) a flank collapse of the Cumbre Vieja Volcano (CVV) in the Canary Islands (with a 80 km volume representing the most likely extreme scenario); (2) an M9 coseismic source in the Puerto Rico Trench (PRT); and (3) a large submarine mass failure (SMF) in the Hudson River canyon of parameters similar to the 165 km historical Currituck slide, which is used as a local proxy for the maximum possible SMF. Simulations are performed with the nonlinear and dispersive long wave model FUNWAVE-TVD, in a series of nested grids of increasing resolution towards the coast, by one-way coupling. Four levels of nested grids are used, from a 1 arc-min spherical coordinate grid in the deep ocean down to a 39-m Cartesian grid in the HRE. Bottom friction coefficients in the finer grids are calibrated for the tide to achieve the local spatially averaged MHW level at high tide in the HRE. Combined tsunami-tide simulations are then performed for four phases of the tide corresponding to each tsunami arriving at Sandy Hook (NJ): 1.5 h ahead, concurrent with, 1.5 h after, and 3 h after the local high tide. These simulations are forced along the offshore boundary of the third-level grid by linearly superposing time series of surface elevation and horizontal currents of the calibrated tide and each tsunami wave train; this is done in deep enough water for a linear superposition to be accurate. Combined tsunami-tide simulations are then performed with FUNWAVE-TVD in this and the finest nested grids. Results show that, for the 3 PMTs, depending on the tide phase, the dynamic simulations lead to no or to a slightly increased inundation in the HRE (by up to 0.15 m depending on location), and to larger currents than for the simulations over a static level; the CRT SMF proxy tsunami is the PMT leading to maximum inundation in the HRE. For all tide phases, nonlinear interactions between tide and tsunami currents modify the elevation, current, and celerity of tsunami wave trains, mostly in the shallower water areas of the HRE where bottom friction dominates, as compared to a linear superposition of wave elevations and currents. We note that, while dynamic simulations predict a slight increase in inundation, this increase may be on the same order as, or even less than sources of uncertainty in the modeling of tsunami sources, such as their initial water elevation, and in bottom friction and bathymetry used in tsunami grids. Nevertheless, results in this paper provide insight into the magnitude and spatial variability of tsunami propagation and impact in the complex inland waterways surrounding New York City, and of their modification by dynamic tidal effects. We conclude that changes in inundation resulting from the inclusion of a dynamic tide in the specific case of the HRE, although of scientific interest, are not significant for tsunami hazard assessment and that the standard approach of specifying a static reference level equal to MHW is conservative. However, in other estuaries with similarly complex bathymetry/topography and stronger tidal currents, a simplified static approach might not be appropriate. [ABSTRACT FROM AUTHOR]

Where coastal tsunami hazard is governed by near-field sources, such as submarine mass failures or meteo-tsunamis, tsunami propagation times may be too small for a detection based on deep or shallow water buoys. To offer sufficient warning time, it has been proposed to implement early warning systems relying on high-frequency (HF) radar remote sensing, that can provide a dense spatial coverage as far offshore as 200-300 km (e.g., for Diginext Ltd.'s Stradivarius radar). Shore-based HF radars have been used to measure nearshore currents (e.g., CODAR SeaSonde system; ), by inverting the Doppler spectral shifts, these cause on ocean waves at the Bragg frequency. Both modeling work and an analysis of radar data following the Tohoku 2011 tsunami, have shown that, given proper detection algorithms, such radars could be used to detect tsunami-induced currents and issue a warning. However, long wave physics is such that tsunami currents will only rise above noise and background currents (i.e., be at least 10-15 cm/s), and become detectable, in fairly shallow water which would limit the direct detection of tsunami currents by HF radar to nearshore areas, unless there is a very wide shallow shelf. Here, we use numerical simulations of both HF radar remote sensing and tsunami propagation to develop and validate a new type of tsunami detection algorithm that does not have these limitations. To simulate the radar backscattered signal, we develop a numerical model including second-order effects in both wind waves and radar signal, with the wave angular frequency being modulated by a time-varying surface current, combining tsunami and background currents. In each 'radar cell', the model represents wind waves with random phases and amplitudes extracted from a specified (wind speed dependent) energy density frequency spectrum, and includes effects of random environmental noise and background current; phases, noise, and background current are extracted from independent Gaussian distributions. The principle of the new algorithm is to compute correlations of HF radar signals measured/simulated in many pairs of distant 'cells' located along the same tsunami wave ray, shifted in time by the tsunami propagation time between these cell locations; both rays and travel time are easily obtained as a function of long wave phase speed and local bathymetry. It is expected that, in the presence of a tsunami current, correlations computed as a function of range and an additional time lag will show a narrow elevated peak near the zero time lag, whereas no pattern in correlation will be observed in the absence of a tsunami current; this is because surface waves and background current are uncorrelated between pair of cells, particularly when time-shifted by the long-wave propagation time. This change in correlation pattern can be used as a threshold for tsunami detection. To validate the algorithm, we first identify key features of tsunami propagation in the Western Mediterranean Basin, where Stradivarius is deployed, by way of direct numerical simulations with a long wave model. Then, for the purpose of validating the algorithm we only model HF radar detection for idealized tsunami wave trains and bathymetry, but verify that such idealized case studies capture well the salient tsunami wave physics. Results show that, in the presence of strong background currents, the proposed method still allows detecting a tsunami with currents as low as 0.05 m/s, whereas a standard direct inversion based on radar signal Doppler spectra fails to reproduce tsunami currents weaker than 0.15-0.2 m/s. Hence, the new algorithm allows detecting tsunami arrival in deeper water, beyond the shelf and further away from the coast, and providing an early warning. Because the standard detection of tsunami currents works well at short range, we envision that, in a field situation, the new algorithm could complement the standard approach of direct near-field detection by providing a warning that a tsunami is approaching, at larger range and in greater depth. This warning would then be confirmed at shorter range by a direct inversion of tsunami currents, from which the magnitude of the tsunami would also estimated. Hence, both algorithms would be complementary. In future work, the algorithm will be applied to actual tsunami case studies performed using a state-of-the-art long wave model, such as briefly presented here in the Mediterranean Basin. [ABSTRACT FROM AUTHOR]

Since the devastating earthquake of 2010 in Haiti, significant efforts have been devoted to estimating future seismic and tsunami hazard in Hispaniola. In 2013, following a workshop of experts, UNESCO commissioned an initial modeling study to assess tsunami hazard, essentially from seismic sources, along the North shore of Hispaniola (NSOH), which is shared by the Republic of Haiti (RH) and the Dominican Republic (DR). The scope of this study included detailed tsunami inundation mapping for two selected critical sites, Cap Haitien in RH and Puerto Plata in DR. Results of this effort are reported here, and, although still limited in scope, they are within the framework and contribute to the advancement of the UNESCO IOC Tsunami and other Coastal Hazards Warning System for the Caribbean and Adjacent Regions (CARIBE EWS; von Hillebrandt-Andrade in Science 341:966-968, ). In similar work done for critical areas of the US east coast (under the auspice of the US National Tsunami Hazard Mitigation Program), the authors have modeled the most extreme far-field tsunami sources in the Atlantic Ocean basin, including: (1) a hypothetical $$M_w$$ 9 seismic event in the Puerto Rico Trench (PRT); (2) a repeat of the historical 1755 $$M_w$$ 9 earthquake in the Azores convergence zone (LSB); and (3) a hypothetical extreme $$450\,\hbox {km}^3$$ flank collapse of the Cumbre Vieja Volcano (CVV) in the Canary Archipelago. Here, tsunami hazard assessment is performed along the NSOH for these three sources, plus two additional near-field coseismic tsunami sources: (1) a $$M_w$$ 8 earthquake in the western segments of the nearshore Septentrional fault (SF), as a repeat of the 1842 event; and (2) a $$M_w$$ 8.7 earthquake occurring in selected segments of the North Hispaniola Thrust Fault (NHTF). Initial tsunami elevations are modeled based on each source's parameters and propagated with FUNWAVE-TVD (a nonlinear and dispersive long-wave Boussinesq model) in a series of increasingly fine-resolution nested grids (from 1 arc-min to 205 m) using a one-way coupling methodology. For the two selected sites, coastal inundation is computed with TELEMAC (a Nonlinear Shallow Water wave model), in finer-resolution (12-30 m) unstructured nested grids. While for the EC, PRT is a far-field source, for RH and DR, this would be local source as some of the NSOH would be affected within 1 h or is within 200 km of the PRT. This is per definitions of UNESCO IOC. Regional goes from 200 to 1000 km and within 1 and 3 h, and distant is greater than 3 h and more than 1000 km. We find that among the far-field sources CVV causes the largest impact, with up to 20-m runup at the critical sites while PRT, which is a local source for the NSOH, only causes up to 4-m runup due to its directionality; PRT, however, has both a much shorter return period and would impact the NSOH within 30 min of the earthquake. Among near-field sources, the SF event, as could be expected from a strike-slip fault, only causes a small tsunami, but the NHTF event causes up to 12-m runup in the critical sites, with the tsunami arriving within minutes of the earthquake. Hence, the latter event can be considered as the 'Probable Maximum Tsunami' (PMT; following, e.g., the US Nuclear Regulatory Commission terminology) for the NSOH. Results of detailed coastal modeling for this PMT can be used to develop maps of vulnerability for the critical sites and prepare for mitigating measures and evacuation; a few examples of such maps are given in the paper. Although a number of earlier studies have dealt with each of the far-field tsunami sources, the modeling of their impact on the NSOH and that of the near-field sources, presented here as part of a comprehensive tsunami hazard assessment study, are novel. Future work should model additional coastal sites and may consider effects of tsunamis generated by near-field submarine mass failures. [ABSTRACT FROM AUTHOR]

The coastal states and territories of the United States (US) are vulnerable to devastating tsunamis from near-field or far-field coseismic and underwater/subaerial landslide sources. Following the catastrophic 2004 Indian Ocean tsunami, the National Tsunami Hazard Mitigation Program (NTHMP) accelerated the development of public safety products for the mitigation of these hazards. In response to this initiative, US coastal states and territories speeded up the process of developing/enhancing/adopting tsunami models that can be used for developing inundation maps and evacuation plans. One of NTHMP's requirements is that all operational and inundation-based numerical (O&I) models used for such purposes be properly validated against established standards to ensure the reliability of tsunami inundation maps as well as to achieve a basic level of consistency between parallel efforts. The validation of several O&I models was considered during a workshop held in 2011 at Texas A&M University (Galveston). This validation was performed based on the existing standard (OAR-PMEL-135), which provides a list of benchmark problems (BPs) covering various tsunami processes that models must meet to be deemed acceptable. Here, we summarize key approaches followed, results, and conclusions of the workshop. Eight distinct tsunami models were validated and cross-compared by using a subset of the BPs listed in the OAR-PMEL-135 standard. Of the several BPs available, only two based on laboratory experiments are detailed here for sake of brevity; since they are considered as sufficiently comprehensive. Average relative errors associated with expected parameters values such as maximum surface amplitude/runup are estimated. The level of agreement with the reference data, reasons for discrepancies between model results, and some of the limitations are discussed. In general, dispersive models were found to perform better than nondispersive models, but differences were relatively small, in part because the BPs mostly featured long waves, such as solitary waves. The largest error found (e.g., the laboratory experiment case of a solitary wave on a simple beach) was 10 % for non-breaking wave conditions and 12 % for breaking conditions; these errors are equal or smaller than the thresholds (10 % and 20 %, respectively) defined by the OAR-PMEL-135 for predicting the surface profile; hence, all models examined here are deemed acceptable for inundation mapping purposes. [ABSTRACT FROM AUTHOR]

The December 26, 2004 tsunami was perhaps the most devastating tsunami in recorded history, causing over 200,000 fatalities and widespread destruction in countries bordering the Indian Ocean. It was generated by the third largest earthquake on record (Mw=9.1–9.3) and was a truly global event, with significant wave action felt around the world. Many measurements of this event were made with seismometers, tide gauges, global positioning system stations, and a few satellite overpasses. There were numerous eyewitness observations and video digital recordings of coastal tsunami impact, as well as subsequent coastal field surveys of runup and flooding. A few ship-based expeditions also took place in the months following the event, to measure and map seafloor disturbances in the epicenter area. Based on these various data sets, recent seismic analysis estimates of rupture propagation speed, and other seismological and geological constraints, we develop a calibrated tsunami source, in terms of coseismic seafloor displacement and rupture timing along 1,200 km of the Andaman–Sunda trench. This source is used to build a numerical model of tsunami generation, propagation, and coastal flooding for the December 26, 2004 event. Frequency dispersion effects having been identified in the deep water tsunami wavetrain, we simulate tsunami propagation and coastal impact with a fully nonlinear and dispersive Boussinesq model (FUNWAVE). The tsunami source is specified in this model as a series of discrete, properly parameterized, dislocation source segments [Okada, 1985, Bull. Seismol. Soc. Am., 75(4), 1135–1154], triggered in a time sequence spanning about 1,200 s. ETOPO2’s bottom bathymetry and land topography are specified in the modeled ocean basin, supplemented by more accurate and denser data in selected coastal areas (e.g., Thailand). A 1 min grid is used for tsunami simulations over the Indian Ocean basin, which is fine enough to model tsunami generation and propagation to nearshore areas. Surface elevations simulated in the model agree well, in both amplitude and timing, with measurements at tide gauges, one satellite transect, and ranges of runup values. These results validate our tsunami source and simulations of the December 26, 2004 event and indicate these can be used to conduct more detailed case studies, for specific coastal areas. In fact, part of the development of our proposed source already benefitted from such regional simulations performed on a finer grid (15 s), as part of a Thailand case study, in which higher frequency waves could be modeled (Ioualalen et al. 2007, J. Geophys. Res., 122, C07024). Finally, by running a non-dispersive version of FUNWAVE, we estimate dispersive effects on maximum deep water elevations to be more than 20% in some areas. We believe that work such as this, in which we achieve a better understanding through modeling of the catastrophic December 26, 2004 event, will help the scientific community better predict and mitigate any such future disaster. This will be achieved through a combination of forecasting models with adequate warning systems, and proper education of the local populations. Such work must be urgently done in light of the certitude that large, potentially tsunamogenic, earthquakes occur along all similar megathrust faults, with a periodicity of a few centuries. [ABSTRACT FROM AUTHOR]

Large scale, three-dimensional, laboratory experiments are performed to study tsunami generation by rigid underwater landslides. The main purpose of these experiments is to both gain insight into landslide tsunami generation processes and provide data for subsequent validation of a three-dimensional numerical model. In each experiment a smooth and streamlined rigid body slides down a plane slope, starting from different initial submergence depths, and generates surface waves. Different conditions of wave nonlinearity and dispersion are generated by varying the model slide initial submergence depth. Surface elevations are measured with capacitance gauges. Runup is measured at the tank axis using a video camera. Landslide acceleration is measured with a microaccelerometer embedded within the model slide, and its time of passage is further recorded at three locations down the slope. The repeatability of experiments is very good. Landslide kinematics is inferred from these measurements and an analytic law of motion is derived, based on which the slide added mass and drag coefficients are computed. Characteristic distance and time of slide motion, as well as a characteristic tsunami wavelength, are parameters derived from these analyses. Measured wave elevations yield characteristic tsunami amplitudes, which are found to be well predicted by empirical equations derived in earlier work, based on two-dimensional numerical computations. The strongly dispersive nature and directionality of tsunamis generated by underwater landslides is confirmed by wave measurements at gauges. Measured coastal runup is analyzed and found to correlate well with initial slide submergence depth or characteristic tsunami amplitude. [ABSTRACT FROM AUTHOR]

Based on numerical simulations presented in Part I, we derive predictive empirical equations describing tsunami generation by submarine mass failure (SMF) that are only valid in the vicinity of the tsunami sources. We give equations for slides and slumps, along with some cautions about their appropriate use. We further discuss results obtained here and in Part I and their practical application to case studies. We show that initial acceleration is the primary parameter describing SMF center of mass motion during tsunami generation. We explain an apparent paradox, raised in Part I, in slump center of mass motion, whereby the distance traveled is proportional to shear strength along the failure plane. We stress that the usefulness of predictive equations depends on the quality of the parameters they rely on. Parameter ranges are discussed in the paper, and we propose a method to estimate slump motion and shear strength and discuss SMF thickness to length values, for case studies. We derive the analytical tools needed to characterize SMF tsunami sources in propagation models. Specifically, we quantify three-dimensional (3D) effects on tsunami characteristic amplitude, and we propose an analytical method to specify initial 3D tsunami elevations, shortly after tsunami generation, in long wave tsunami propagation models. This corresponds to treating SMF tsunami sources like coseismic displacement tsunami sources. We conduct four case studies of SMF tsunamis and show that our predictive equations can provide rapid rough estimates of overall tsunami observations that might be useful in crisis situations, when time is too short to run propagation models. Thus, for each case, we show that the characteristic tsunami amplitude is a reasonable predictor of maximum runup in actual 3D geometry. We refer to the latter observation as the correspondence principle, which we propose to apply for rapid tsunami hazard assessment, in combination with the predictive tsunami amplitude equations. [ABSTRACT FROM AUTHOR]

Numerical simulations are performed with a two-dimensional (2D) fully nonlinear potential flow (FNPF) model for tsunami generation by two idealized types of submarine mass failure (SMF): underwater slides and slumps. These simulations feature rigid or deforming SMFs with a Gaussian cross section, translating down a plane slope. In each case, the SMF center of mass motion is expressed as a function of geometric, hydrodynamic, and material parameters, following a simple wavemaker formalism, and prescribed as a boundary condition in the FNPF model. Tsunami amplitudes and runup are obtained from computed free surface elevations. Model results are experimentally validated for a rigid 2D slide. Sensitivity studies are performed to estimate the effects of SMF–shape, type, and initial submergence depth—on the generated tsunamis. A strong SMF deformation during motion is shown to significantly enhance tsunami generation, particularly in the far-field. Typical slumps are shown to generate smaller tsunamis than corresponding slides. Both tsunami amplitude and runup are shown to depend strongly on initial SMF submergence depth. For the selected SMF idealized geometry, this dependence is simply expressed by power laws. Other sensitivity analyses are presented in a companion paper, and results from numerical simulations are converted into empirical curve fits predicting characteristic tsunami amplitudes as functions of nondimensional governing parameters. It should be stressed that these empirical formulas are only valid in the vicinity of the tsunami sources and, because of the complexity of the problem, many simplifications were necessary. It is further shown in the companion paper how 2D results can be modified to account for three-dimensional tsunami generation and used for quickly estimating tsunami hazard or for performing simple case studies. [ABSTRACT FROM AUTHOR]

Abstract: This work reports on the application and experimental validation, for idealized geometries, of a multiple-fluid Navier–Stokes model of waves generated by rigid and deforming slides, with the purpose of improving predictive simulations of landslide tsunamis. In such simulations, the computational domain is divided into water, air, and slide regions, all treated as Newtonian fluids. For rigid slides, a penalty method allows for parts of the fluid domain to behave as a solid. With the latter method, the coupling between a rigid slide and water is implicitly computed (rather than specifying a known slide kinematics). Two different Volume of Fluid algorithms are tested for tracking interfaces between actual fluid regions. The simulated kinematics of a semi-elliptical block, moving down a water covered plane slope, is first compared to an earlier analytical solution. Results for the vertical fall of a rectangular block in water are then compared to earlier experimental results. Finally, more realistic simulations of two- and three-dimensional wedges sliding down an incline are compared to earlier experiments. Overall, in all cases, solid block velocities and free surface deformations are accurately reproduced in the model, provided that a sufficiently resolved discretization is used. The potential of the model is then illustrated on more complex scenarios involving waves caused by multi-block or deformable slides. [Copyright &y& Elsevier]