Search

Abstract: Tsunamis are frequent events in modern marine environments but evidence of their passing in subtidal ancient epicontinental sea deposits is elusive. It has been suggested that this is due to mis-identification or poor preservation potential. Herein a numerical modelling approach is used to show that tsunami propagation in one large ancient epicontinental sea was hindered by the damping effect of shallow bathymetries and reflection, refraction and diffraction from emergent landmasses. The Imperial College Ocean Model (ICOM) is used for this study and is first validated against data from the Sumatra–Andaman Tsunami of December 2004. A palaeobathymetric dataset is then presented for the Hettangian (Lower Jurassic) Laurasian Seaway with idealised tsunami sources situated on the continental shelf and within the adjacent oceanic basin. Results show that tsunamis forced from within ocean basins adjacent to the epicontinental sea are rapidly attenuated over the continental slope and fail to propagate great distances onto the shelf. Similarly, the sedimentological effect of tsunamis forced from within the epicontinental sea is also restricted. It is concluded that tsunami deposits in ancient epicontinental seas are most likely to occur in relative proximity to the source region and this must contribute to their scarcity in the geological record. [Copyright &y& Elsevier]

Tsunami generated by submarine slides are arguably an under-considered risk in comparison to earthquake-generated tsunami. Numerical simulations of submarine slide-generated waves can be used to identify the important factors in determining wave characteristics. Here we use Fluidity, an open source finite element code, to simulate waves generated by deformable submarine slides. Fluidity uses flexible unstructured meshes combined with adaptivity which alters the mesh topology and resolution based on the simulation state, focussing or reducing resolution, when and where it is required. Fluidity also allows a number of different numerical approaches to be taken to simulate submarine slide deformation, free-surface representation, and wave generation within the same numerical framework. In this work we use a multi-material approach, considering either two materials (slide and water with a free surface) or three materials (slide, water and air), as well as a sediment model (sediment, water and free surface) approach. In all cases the slide is treated as a viscous fluid. Our results are shown to be consistent with laboratory experiments using a deformable submarine slide, and demonstrate good agreement when compared with other numerical models. The three different approaches for simulating submarine slide dynamics and tsunami wave generation produce similar waveforms and slide deformation geometries. However, each has its own merits depending on the application. Mesh adaptivity is shown to be able to reduce the computational cost without compromising the accuracy of results. [ABSTRACT FROM AUTHOR]

The ∼8.15 ka Storegga submarine slide was a large (∼3000 km 3 ), tsunamigenic slide off the coast of Norway. The resulting tsunami had run-up heights of around 10–20 m on the Norwegian coast, over 12 m in Shetland, 3–6 m on the Scottish mainland coast and reached as far as Greenland. Accurate numerical simulations of Storegga require high spatial resolution near the coasts, particularly near tsunami run-up observations, and also in the slide region. However, as the computational domain must span the whole of the Norwegian-Greenland sea, employing uniformly high spatial resolution is computationally prohibitive. To overcome this problem, we present a multiscale numerical model of the Storegga slide-generated tsunami where spatial resolution varies from 500 m to 50 km across the entire Norwegian-Greenland sea domain to optimally resolve the slide region, important coastlines and bathymetric changes. We compare results from our multiscale model to previous results using constant-resolution models and show that accounting for changes in bathymetry since 8.15 ka, neglected in previous numerical studies of the Storegga slide-tsunami, improves the agreement between the model and inferred run-up heights in specific locations, especially in the Shetlands, where maximum run-up height increased from 8 m (modern bathymetry) to 13 m (palaeobathymetry). By tracking the Storegga tsunami as far south as the southern North sea, we also found that wave heights were high enough to inundate Doggerland, an island in the southern North Sea prior to sea level rise over the last 8 ka. [ABSTRACT FROM AUTHOR]

The development of reliable, sophisticated hydro-morphodynamic models is essential for protecting the coastal environment against hazards such as flooding and erosion. There exists a high degree of uncertainty associated with the application of these models, in part due to incomplete knowledge of various physical, empirical and numerical closure related parameters in both the hydrodynamic and morphodynamic solvers. This uncertainty can be addressed through the application of adjoint methods. These have the notable advantage that the number and/or dimension of the uncertain parameters has almost no effect on the computational cost associated with calculating the model sensitivities. Here, we develop the first freely available and fully flexible adjoint hydro-morphodynamic model framework. This flexibility is achieved through using the pyadjoint library, which allows us to assess the uncertainty of any parameter with respect to any model functional, without further code implementation. The model is developed within the coastal ocean model Thetis constructed using the finite element code-generation library Firedrake. We present examples of how this framework can perform sensitivity analysis, inversion and calibration for a range of uncertain parameters based on the final bedlevel. These results are verified using so-called dual-twin experiments, where the 'correct' parameter value is used in the generation of synthetic model test data, but is unknown to the model in subsequent testing. Moreover, we show that inversion and calibration with experimental data using our framework produces physically sensible optimum parameters and that these parameters always lead to more accurate results. In particular, we demonstrate how our adjoint framework can be applied to a tsunami-like event to invert for the tsunami wave from sediment deposits. • Development of first free and fully flexible adjoint hydro-morphodynamic framework. • Verification of adjoint capabilities using dual-twin experiments. • Application of adjoint methods to analyse sensitivity to spatially-varying parameters. • Use of adjoint methods to calibrate for several uncertain parameters in an experiment. • Inversion of a tsunami-like wave using final sediment deposits from experimental data. [ABSTRACT FROM AUTHOR]

Large tsunamis can be generated by submarine slides, but these events are rare on human timescales and challenging to observe. Experiments and numerical modelling offer methods to understand the mechanisms by which they generate waves and what the potential hazard might be. However, to fully capture the complex waveform generated by a submarine slide, the slide dynamics must also be accurately modelled. It is computationally difficult to model both a three-dimensional submarine slide whilst simultaneously simulating oceanic-scale tsunamis. Past studies have either coupled localised models of the slide generation to oceanic-scale tsunami simulations or simplified the slide dynamics. Here, we present a new methodology of model coupling that generates the wave in the ocean-scale model via boundary-condition coupling of a two-dimensional dynamic slide simulation. We verify our coupling methodology by comparing model results to a previous simulation of a tsunamigenic slide in the Gulf of Mexico. We then examine the effect of slide deformation on the risk posed by hypothetical submarine slides around the UK. We show the deformable submarine slide simulations produce larger waves than the solid slide simulations due to the details of acceleration and velocity of the slide, although lateral spreading is not modelled. This work offers a new methodology for simulating oceanic-scale tsunamis caused by submarine slides using the output of a two-dimensional, multi-material simulation as input into a three-dimensional ocean model. This facilitates future exploration of the tsunami risk posed by tsunamigenic submarine slides that affect coastlines not normally prone to tsunamis. • Numerical modelling framework simulates tsunami waves generated by deformable submarine slides. • Predicted waves agree with numerical model for waves generated by submarine slide in Gulf of Mexico. • Used methodology to simulate waves from hypothetical submarine slides around the UK. • We show how slide speed, acceleration, deformation and direction affect the wave. [ABSTRACT FROM AUTHOR]