Search

The Ligurian coast, located at the French-Italian border, is densely populated as well as a touristic area. It is also a location where earthquakes and underwater landslides are recurrent. The nature of the local tsunamigenesis is therefore a legitimate question, because no tsunami warning system can resolve tsunami arrival times of a few minutes, which is the case for the area. As far as the seismicity of the area is concerned, the frequent recurrent earthquakes are generally of moderate magnitude: most of them are lower than M (w) 5. However, the relatively large M (w) 6.9 earthquake (Larroque et al., in Geophys J Int, 2012. doi: 10.1111/j.1365-246X.2012.05498.x) that occurred on the February 23, 1887, offshore of Imperia (Italian Riviera) is quite emblematic. This unusual event for the region merits a complete study: the quantification of its rupture mechanism is essential (1) to understand the regional active deformation, but also (2) to evaluate its tsunamigenesis potential by deriving relevant rupture scenarios obtained from our knowledge of the event; for that purpose the event is extensively described here. The first point has been the subject of quite a few studies based on the seismotectonics of the area. The last documented approach has been completed by Larroque et al. (Geophys J Int, 2012. doi: 10.1111/j.1365-246X.2012.05498.x) who proposed a rupture scenario involving a reverse faulting along a north dipping fault and favoring a M (w) 6.9 magnitude. In the present paper (1) we study the accuracy of their solutions in relation to the computational grid spacing and the dispersive/nondispersive parameterization, (2) based on an uncertainty on the recorded wave amplitude of the Genoa tide gauge they used, we propose a M (w) 6.7 earthquake magnitude solution for the event (the kinematics is unchanged), co-existing with the M (w) 6.9, (3) we evaluate the tsunami coastal impact of the 1887 event, and (4) we test a range of possible ruptures that local faults may undergo in order to propose a synoptic mapping of the tsunami threat in the area. The spatial distribution of the maximum wave height (MWH) is provided with a tentative identification of the processes that are responsible for it. This latter issue is imperative in order to make our mapping as generic as possible in the framework of our deterministic approach (based on realistic scenarios and not on ensemble statistics). The predictions suggest that the wave impact is mostly local, considering the relatively moderate size of the rupture planes. Although the present-day seismicity in this region is moderate, stronger earthquakes (M > 6.5) have occurred in the past. The studied scenarios show that for such events specific localities along the French-Italian Riviera may experience very significant MWH related to the shallow focal depth tested for such scenarios. We may reasonably conclude that the tsunami threat is relatively significant and uniform at the Italian side of the Riviera (from Ventimiglia to Imperia), while it is more localized (sporadic) at the French side from Antibes to Menton with, however, higher local level of inundation, e.g., Nice city center.

Following the 2007, August 15th, M-w 8.0, Pisco earthquake in central Peru, SLADEN et al. (J Geophys Res 115: B02405, 2010) have derived several slip models of this event. They inverted teleseismic data together with geodetic (InSAR) measurements to look for the co-seismic slip distribution on the fault plane, considering those data sets separately or jointly. But how close to the real slip distribution are those inverted slip models? To answer this crucial question, the authors generated some tsunami records based on their slip models and compared them to DART buoys, tsunami records, and available runup data. Such an approach requires a robust and accurate tsunami model (non-linear, dispersive, accurate bathymetry and topography, etc.) otherwise the differences between the data and the model may be attributed to the slip models themselves, though they arise from an incomplete tsunami simulation. The accuracy of a numerical tsunami simulation strongly depends, among others, on two important constraints: (i) A fine computational grid (and thus the bathymetry and topography data sets used) which is not always available, unfortunately, and (ii) a realistic tsunami propagation model including dispersion. Here, we extend Sladen's work using newly available data, namely a tide gauge record at Callao (Lima harbor) and the Chilean DART buoy record, while considering a complete set of runup data along with a more realistic tsunami numerical that accounts for dispersion, and also considering a fine-resolution computational grid, which is essential. Through these accurate numerical simulations we infer that the InSAR-based model is in better agreement with the tsunami data, studying the case of the Pisco earthquake indicating that geodetic data seems essential to recover the final co-seismic slip distribution on the rupture plane. Slip models based on teleseismic data are unable to describe the observed tsunami, suggesting that a significant amount of co-seismic slip may have been aseismic. Finally, we compute the runup distribution along the central part of the Peruvian coast to better understand the wave amplification/attenuation processes of the tsunami generated by the Pisco earthquake.