Search

Highlights • We identify the largest active fault systems of the Alboran Basin. • Characterization of faults is key for accurate tsunamigenic potential estimations. • Alboran largest fault systems may generate Mw > 7 earthquakes. • These earthquakes have the potential to generate significant tsunami waves approaching the coast. Abstract The westernmost Mediterranean hosts part of the plate boundary between the European and African tectonic plates. Based on the scattered instrumental seismicity, this boundary has been traditionally interpreted as a wide zone of diffuse deformation. However, recent seismic images and seafloor mapping studies support that most of the plate convergence may be accommodated in a few tectonic structures, rather than in a broad region. Historical earthquakes with magnitudes Mw > 6 and historical tsunamis support that the low-to-moderate instrumental seismicity might also have led to underestimation of the seismogenic and tsunamigenic potential of the area. We evaluate the largest active faults of the westernmost Mediterranean: the reverse Alboran Ridge, and the strike-slip Carboneras, Yusuf and Al-Idrissi fault systems. For the first time, we use a dense grid of modern seismic data to characterize the entire dimensions of the main fault systems, accurately describe the geometry of these structures and estimate their seismic source parameters. Tsunami scenarios have been tested based on 3D-surfaces and seismic source parameters, using both uniform and heterogeneous slip distributions. The comparison of our results with previous studies, based on limited information on the fault geometry and kinematics, indicates that accurate fault geometries and heterogeneous slip distributions are needed to properly assess the seismic and tsunamigenic potential in this area. Based on fault scaling relations, the four fault systems have a large seismogenic potential, being able to generate earthquakes with Mw > 7. The reverse Alboran Ridge Fault System has

The Alboran Basin is located in the westernmost Mediterranean Sea. This basin was formed during the Miocene, and since the late Miocene, has been deformed due to the Iberia – Africa tectonic plates convergence, producing the contractive reorganization of some structures at the basin. Thus, the Alboran Basin is a seismically active area, which hosts the plate boundary between the European and African tectonic plates. This plate boundary has been traditionally considered a wide deformation zone, in which several small faults are accommodating the deformation. Based on a modern set of active seismic data, we were able for the first time to quantify the total slip accommodated by the most prominent tectonic structures of the area, late Miocene - early Pliocene in age. Our results show that the estimated total slip accommodated by the main fault systems may be similar (with error bounds) to the estimated plate convergence value since the Messinian time (~24 km). Thus, slip on that faults may have accommodated most of the Iberian – African plate convergence during the Plio-Quaternary, revealing that the contractive reorganization of the Alboran basin is focused on a few first-order structures that act as lithospheric boundaries, rather than widespread and diffuse along the entire basin. These results have implications not only for kinematic and geodynamic models, but also for seismic and tsunami hazard assessments. Using the most complete dataset until the date, we performed a revision of the geometry and characteristics of the main fault systems offshore. Based on this data, we perform a first appraisal of the seismogenic and tsunamigenic potential of the main fault systems offshore. Our simulations show that the seismogenic and tsunamigenic potential of the offshore structures of the Alboran Basin may be underestimated, and a further characterization of their associated hazard is needed

Highlights • We identify the largest active fault systems of the Alboran Basin. • Characterization of faults is key for accurate tsunamigenic potential estimations. • Alboran largest fault systems may generate Mw > 7 earthquakes. • These earthquakes have the potential to generate significant tsunami waves approaching the coast. Abstract The westernmost Mediterranean hosts part of the plate boundary between the European and African tectonic plates. Based on the scattered instrumental seismicity, this boundary has been traditionally interpreted as a wide zone of diffuse deformation. However, recent seismic images and seafloor mapping studies support that most of the plate convergence may be accommodated in a few tectonic structures, rather than in a broad region. Historical earthquakes with magnitudes Mw > 6 and historical tsunamis support that the low-to-moderate instrumental seismicity might also have led to underestimation of the seismogenic and tsunamigenic potential of the area. We evaluate the largest active faults of the westernmost Mediterranean: the reverse Alboran Ridge, and the strike-slip Carboneras, Yusuf and Al-Idrissi fault systems. For the first time, we use a dense grid of modern seismic data to characterize the entire dimensions of the main fault systems, accurately describe the geometry of these structures and estimate their seismic source parameters. Tsunami scenarios have been tested based on 3D-surfaces and seismic source parameters, using both uniform and heterogeneous slip distributions. The comparison of our results with previous studies, based on limited information on the fault geometry and kinematics, indicates that accurate fault geometries and heterogeneous slip distributions are needed to properly assess the seismic and tsunamigenic potential in this area. Based on fault scaling relations, the four fault systems have a large seismogenic potential, being able to generate earthquakes with Mw > 7. The reverse Alboran Ridge Fault System has

The westernmost Mediterranean hosts part of the plate boundary between the European and African tectonic plates. Based on the scattered instrumental seismicity, this boundary has been traditionally interpreted as a wide zone of diffuse deformation. However, recent seismic images and seafloor mapping studies support that most of the plate convergence may be accommodated in a few tectonic structures, rather than in a broad region. Historical earthquakes with magnitudes Mw > 6 and historical tsunamis support that the low-to-moderate instrumental seismicity might also have led to underestimation of the seismogenic and tsunamigenic potential of the area. We evaluate the largest active faults of the westernmost Mediterranean: the reverse Alboran Ridge, and the strike-slip Carboneras, Yusuf and Al-Idrissi fault systems. For the first time, we use a dense grid of modern seismic data to characterize the entire dimensions of the main fault systems, accurately describe the geometry of these structures and estimate their seismic source parameters. Tsunami scenarios have been tested based on 3D-surfaces and seismic source parameters, using both uniform and heterogeneous slip distributions. The comparison of our results with previous studies, based on limited information on the fault geometry and kinematics, indicates that accurate fault geometries and heterogeneous slip distributions are needed to properly assess the seismic and tsunamigenic potential in this area. Based on fault scaling relations, the four fault systems have a large seismogenic potential, being able to generate earthquakes with Mw > 7. The reverse Alboran Ridge Fault System has the largest tsunamigenic potential, being able to generate a tsunami wave amplitude greater than 3 m in front of the coasts of Southern Spain and Northern Africa

The SW Iberian margin is one of the most seismogenic and tsunamigenic areas in W-Europe, where large historical and instrumental destructive events occurred. To evaluate the sensitivity of the tsunami impact on the coast of SW Iberia and NW Morocco to the fault geometry and slip distribution for local earthquakes, we carried out a set of tsunami simulations considering some of the main known active crustal faults in the region: the Gorringe Bank (GBF), Marquês de Pombal (MPF), Horseshoe (HF), North Coral Patch (NCPF) and South Coral Patch (SCPF) thrust faults, and the Lineament South (LSF) strike-slip fault. We started by considering for all of them relatively simple planar faults featuring with uniform slip distribution; we then used a more complex 3D fault geometry for the faults constrained with a large 2D multichannel seismic dataset (MPF, HF, NCPF, and SCPF); and finally, we used various heterogeneous slip distributions for the HF. Our results show that using more complex 3D fault geometries and slip distributions, the peak wave height at the coastline can double compared to simpler tsunami source scenarios from planar fault geometries. Existing tsunami hazard models in the region use homogeneous slip distributions on planar faults as initial conditions for tsunami simulations and therefore underestimate tsunami hazard. Complex 3D fault geometries and non-uniform slip distribution should be considered in future tsunami hazard updates. The tsunami simulations also support the finding that submarine canyons attenuate the wave height reaching the coastline, while submarine ridges and shallow shelves have the opposite effect.

[EN] The SW Iberian margin is one of the most seismogenic and tsunamigenic areas in Western Europe, where large historical and instrumental destructive earthquakes and tsunamis have occurred. For this reason, we present a new study to evaluate the sensitivity of the tsunami impact on the coast of SW Iberia and NW Morocco to the fault geometry and slip distribution for local earthquakes. We carried out several tsunami simulations considering the detailed geometry of the main active faults in the region: the Gorringe Bank (GBF), Marquês de Pombal (MPF), Horseshoe (HF), North Coral Patch (NCPF) and South Coral Patch (SCPF) thrust faults, and the Lineament South (LSF) strike-slip fault. We started from relatively simple planar faults featuring uniform slip for all faults; then, based on a large dataset of 2D multi-channel seismic profiles, 3D complex sub-surface models of the fault planes (MPF, HF, NCPF, and SCPF) have been generated. Finally, we used several heterogeneous slip-distributions for the HF case. Our tsunami models suggest that using more complex fault geometries and heterogeneous slip distributions, the peak wave height (about the sea level) at the coastline can be doubled compared to simpler tsunami source scenarios generated by simple fault-plane geometries, [ES] El margen SO de la Península Ibérica es una de las zonas más sismogénicas y tsunamigénicas de Europa Occidental, donde han ocurrido grandes terremotos y tsunamis destructivos, tanto en periodo histórico como instrumental. Por este motivo, presentamos un nuevo estudio evaluando la sensibilidad de los escenarios de tsunamis a la geometría de los planos de falla y la distribución del deslizamiento aplicado a las costas del SO de Iberia y del NE de Marruecos. En este trabajo realizamos diversas simulaciones de tsunami considerando la geometría en detalle de las principales fallas activas de la región, incluyendo: las fallas inversas de Gorringe Bank (GBF), Marquês de Pombal (MPF), Horseshoe (HF), North Coral Patch (NCPF) y South Coral Patch (SCPF); y la falla transcurrente del Lineament South (LSF). Empezamos considerando planos de falla simples con deslizamiento uniforme para todas las fallas. A continuación, generamos modelos complejos en 3D de los planos de falla para MPF, HF, NCPF y SCPF en base a un gran conjunto de datos de perfiles de sísmica multicanal 2D. Finalmente, utilizamos varias distribuciones de deslizamiento heterogéneo para el caso de HF. Nuestros modelos de tsunami sugieren que utilizando geometrías de falla complejas y distribuciones de deslizamiento heterogéneas, la altura máxima de la ola (sobre el nivel del mar) en la costa se puede duplicar en comparación con escenarios simples, dónde el tsunami es generado por una falla con geometría planar

The SW Iberian margin is one of the most seismogenic and tsunamigenic areas in W-Europe, where large historical and instrumental destructive events occurred. To evaluate the sensitivity of the tsunami impact on the coast of SW Iberia and NW Morocco to the fault geometry and slip distribution for local earthquakes, we carried out a set of tsunami simulations considering some of the main known active crustal faults in the region: the Gorringe Bank (GBF), Marquês de Pombal (MPF), Horseshoe (HF), North Coral Patch (NCPF) and South Coral Patch (SCPF) thrust faults, and the Lineament South (LSF) strike-slip fault. We started by considering for all of them relatively simple planar faults featuring with uniform slip distribution; we then used a more complex 3D fault geometry for the faults constrained with a large 2D multichannel seismic dataset (MPF, HF, NCPF, and SCPF); and finally, we used various heterogeneous slip distributions for the HF. Our results show that using more complex 3D fault geometries and slip distributions, the peak wave height at the coastline can double compared to simpler tsunami source scenarios from planar fault geometries. Existing tsunami hazard models in the region use homogeneous slip distributions on planar faults as initial conditions for tsunami simulations and therefore underestimate tsunami hazard. Complex 3D fault geometries and non-uniform slip distribution should be considered in future tsunami hazard updates. The tsunami simulations also support the finding that submarine canyons attenuate the wave height reaching the coastline, while submarine ridges and shallow shelves have the opposite effect.

The SW Iberian margin is one of the most seismogenic and tsunamigenic areas in W-Europe, where large historical and instrumental destructive events occurred. To evaluate the sensitivity of the tsunami impact on the coast of SW Iberia and NW Morocco to the fault geometry and slip distribution for local earthquakes, we carried out a set of tsunami simulations considering some of the main known active crustal faults in the region: the Gorringe Bank (GBF), Marquês de Pombal (MPF), Horseshoe (HF), North Coral Patch (NCPF) and South Coral Patch (SCPF) thrust faults, and the Lineament South (LSF) strike-slip fault. We started by considering for all of them relatively simple planar faults featuring with uniform slip distribution; we then used a more complex 3D fault geometry for the faults constrained with a large 2D multichannel seismic dataset (MPF, HF, NCPF, and SCPF); and finally, we used various heterogeneous slip distributions for the HF. Our results show that using more complex 3D fault geometries and slip distributions, the peak wave height at the coastline can double compared to simpler tsunami source scenarios from planar fault geometries. Existing tsunami hazard models in the region use homogeneous slip distributions on planar faults as initial conditions for tsunami simulations and therefore underestimate tsunami hazard. Complex 3D fault geometries and non-uniform slip distribution should be considered in future tsunami hazard updates. The tsunami simulations also support the finding that submarine canyons attenuate the wave height reaching the coastline, while submarine ridges and shallow shelves have the opposite effect.

The SW Iberian margin is one of the most seismogenic and tsunamigenic areas in W-Europe, where large historical and instrumental destructive events occurred. To evaluate the sensitivity of the tsunami impact on the coast of SW Iberia and NW Morocco to the fault geometry and slip distribution for local earthquakes, we carried out a set of tsunami simulations considering some of the main known active crustal faults in the region: the Gorringe Bank (GBF), Marquês de Pombal (MPF), Horseshoe (HF), North Coral Patch (NCPF) and South Coral Patch (SCPF) thrust faults, and the Lineament South strike-slip fault. We started by considering for all of them relatively simple planar faults featuring with uniform slip distribution; we then used a more complex 3D fault geometry for the faults constrained with a large 2D multichannel seismic dataset (MPF, HF, NCPF, and SCPF); and finally, we used various heterogeneous slip distributions for the HF. Our results show that using more complex 3D fault geometries and slip distributions, the peak wave height at the coastline can double compared to simpler tsunami source scenarios from planar fault geometries. Existing tsunami hazard models in the region use homogeneous slip distributions on planar faults as initial conditions for tsunami simulations and therefore underestimate tsunami hazard. Complex 3D fault geometries and non-uniform slip distribution should be considered in future tsunami hazard updates. The tsunami simulations also support the finding that submarine canyons attenuate the wave height reaching the coastline, while submarine ridges and shallow shelves have the opposite effect

The complexity of coseismic slip distributions influences the tsunami hazard posed by local and, to a certain extent, distant tsunami sources. Large slip concentrated in shallow patches was observed in recent tsunamigenic earthquakes, possibly due to dynamic amplification near the free surface, variable frictional conditions or other factors. We propose a method for incorporating enhanced shallow slip for subduction earthquakes while preventing systematic slip excess at shallow depths over one or more seismic cycles. The method uses the classic k(-2) stochastic slip distributions, augmented by shallow slip amplification. It is necessary for deep events with lower slip to occur more often than shallow ones with amplified slip to balance the long-term cumulative slip. We evaluate the impact of this approach on tsunami hazard in the central and eastern Mediterranean Sea adopting a realistic 3D geometry for three subduction zones, by using it to model similar to 150,000 earthquakes with M-w from 6.0 to 9.0. We combine earthquake rates, depth-dependent slip distributions, tsunami modeling, and epistemic uncertainty through an ensemble modeling technique. We found that the mean hazard curves obtained with our method show enhanced probabilities for larger inundation heights as compared to the curves derived from depth-independent slip distributions. Our approach is completely general and can be applied to any subduction zone in the world.

Finite-fault models for the 2010 Mw 8.8 Maule, Chile earthquake indicate bilateral rupture with large-slip patches located north and south of the epicenter. Previous studies also show that this event features significant slip in the shallow part of the megathrust, which is revealed through correction of the forward tsunami modeling scheme used in tsunami inversions. The presence of shallow slip is consistent with the coseismic seafloor deformation measured off the Maule region adjacent to the trench and confirms that tsunami observations are particularly important for constraining far-offshore slip. Here, we benchmark the method of Optimal Time Alignment (OTA) of the tsunami waveforms in the joint inversion of tsunami (DART and tide-gauges) and geodetic (GPS, InSAR, land-leveling) observations for this event. We test the application of OTA to the tsunami Green’s functions used in a previous inversion. Through a suite of synthetic tests we show that if the bias in the forward model is comprised only of delays in the tsunami signals, the OTA can correct them precisely, independently of the sensors (DART or coastal tide-gauges) and, to the first-order, of the bathymetric model used. The same suite of experiments is repeated for the real case of the 2010 Maule earthquake where, despite the results of the synthetic tests, DARTs are shown to outperform tide-gauges. This gives an indication of the relative weights to be assigned when jointly inverting the two types of data. Moreover, we show that using OTA is preferable to subjectively correcting possible time mismatch of the tsunami waveforms. The results for the source model of the Maule earthquake show that using just the first-order modeling correction introduced by OTA confirms the bilateral rupture pattern around the epicenter, and, most importantly, shifts the inferred northern patch of slip to a shallower position consistent with the slip models obtained by applying more complex physics-based corrections to the tsu

Finite-fault models for the 2010 Mw 8.8 Maule, Chile earthquake indicate bilateral rupture with large-slip patches located north and south of the epicenter. Previous studies also show that this event features significant slip in the shallow part of the megathrust, which is revealed through correction of the forward tsunami modeling scheme used in tsunami inversions. The presence of shallow slip is consistent with the coseismic seafloor deformation measured off the Maule region adjacent to the trench and confirms that tsunami observations are particularly important for constraining far-offshore slip. Here, we benchmark the method of Optimal Time Alignment (OTA) of the tsunami waveforms in the joint inversion of tsunami (DART and tide-gauges) and geodetic (GPS, InSAR, land-leveling) observations for this event. We test the application of OTA to the tsunami Green’s functions used in a previous inversion. Through a suite of synthetic tests we show that if the bias in the forward model is comprised only of delays in the tsunami signals, the OTA can correct them precisely, independently of the sensors (DART or coastal tide-gauges) and, to the first-order, of the bathymetric model used. The same suite of experiments is repeated for the real case of the 2010 Maule earthquake where, despite the results of the synthetic tests, DARTs are shown to outperform tide-gauges. This gives an indication of the relative weights to be assigned when jointly inverting the two types of data. Moreover, we show that using OTA is preferable to subjectively correcting possible time mismatch of the tsunami waveforms. The results for the source model of the Maule earthquake show that using just the first-order modeling correction introduced by OTA confirms the bilateral rupture pattern around the epicenter, and, most importantly, shifts the inferred northern patch of slip to a shallower position consistent with the slip models obtained by applying more complex physics-based corrections to the tsu

The complexity of coseismic slip distributions influences the tsunami hazard posed by local and, to a certain extent, distant tsunami sources. Large slip concentrated in shallow patches was observed in recent tsunamigenic earthquakes, possibly due to dynamic amplification near the free surface, variable frictional conditions or other factors. We propose a method for incorporating enhanced shallow slip for subduction earthquakes while preventing systematic slip excess at shallow depths over one or more seismic cycles. The method uses the classic k(-2) stochastic slip distributions, augmented by shallow slip amplification. It is necessary for deep events with lower slip to occur more often than shallow ones with amplified slip to balance the long-term cumulative slip. We evaluate the impact of this approach on tsunami hazard in the central and eastern Mediterranean Sea adopting a realistic 3D geometry for three subduction zones, by using it to model similar to 150,000 earthquakes with M-w from 6.0 to 9.0. We combine earthquake rates, depth-dependent slip distributions, tsunami modeling, and epistemic uncertainty through an ensemble modeling technique. We found that the mean hazard curves obtained with our method show enhanced probabilities for larger inundation heights as compared to the curves derived from depth-independent slip distributions. Our approach is completely general and can be applied to any subduction zone in the world.

Finite-fault models for the 2010 Mw 8.8 Maule, Chile earthquake indicate bilateral rupture with large-slip patches located north and south of the epicenter. Previous studies also show that this event features significant slip in the shallow part of the megathrust, which is revealed through correction of the forward tsunami modeling scheme used in tsunami inversions. The presence of shallow slip is consistent with the coseismic seafloor deformation measured off the Maule region adjacent to the trench and confirms that tsunami observations are particularly important for constraining far-offshore slip. Here, we benchmark the method of Optimal Time Alignment (OTA) of the tsunami waveforms in the joint inversion of tsunami (DART and tide-gauges) and geodetic (GPS, InSAR, land-leveling) observations for this event. We test the application of OTA to the tsunami Green’s functions used in a previous inversion. Through a suite of synthetic tests we show that if the bias in the forward model is comprised only of delays in the tsunami signals, the OTA can correct them precisely, independently of the sensors (DART or coastal tide-gauges) and, to the first-order, of the bathymetric model used. The same suite of experiments is repeated for the real case of the 2010 Maule earthquake where, despite the results of the synthetic tests, DARTs are shown to outperform tide-gauges. This gives an indication of the relative weights to be assigned when jointly inverting the two types of data. Moreover, we show that using OTA is preferable to subjectively correcting possible time mismatch of the tsunami waveforms. The results for the source model of the Maule earthquake show that using just the first-order modeling correction introduced by OTA confirms the bilateral rupture pattern around the epicenter, and, most importantly, shifts the inferred northern patch of slip to a shallower position consistent with the slip models obtained by applying more complex physics-based corrections to the tsu

The complexity of coseismic slip distributions influences the tsunami hazard posed by local and, to a certain extent, distant tsunami sources. Large slip concentrated in shallow patches was observed in recent tsunamigenic earthquakes, possibly due to dynamic amplification near the free surface, variable frictional conditions or other factors. We propose a method for incorporating enhanced shallow slip for subduction earthquakes while preventing systematic slip excess at shallow depths over one or more seismic cycles. The method uses the classic k(-2) stochastic slip distributions, augmented by shallow slip amplification. It is necessary for deep events with lower slip to occur more often than shallow ones with amplified slip to balance the long-term cumulative slip. We evaluate the impact of this approach on tsunami hazard in the central and eastern Mediterranean Sea adopting a realistic 3D geometry for three subduction zones, by using it to model similar to 150,000 earthquakes with M-w from 6.0 to 9.0. We combine earthquake rates, depth-dependent slip distributions, tsunami modeling, and epistemic uncertainty through an ensemble modeling technique. We found that the mean hazard curves obtained with our method show enhanced probabilities for larger inundation heights as compared to the curves derived from depth-independent slip distributions. Our approach is completely general and can be applied to any subduction zone in the world.

Finite-fault models for the 2010 Mw 8.8 Maule, Chile earthquake indicate bilateral rupture with large-slip patches located north and south of the epicenter. Previous studies also show that this event features significant slip in the shallow part of the megathrust, which is revealed through correction of the forward tsunami modeling scheme used in tsunami inversions. The presence of shallow slip is consistent with the coseismic seafloor deformation measured off the Maule region adjacent to the trench and confirms that tsunami observations are particularly important for constraining far-offshore slip. Here, we benchmark the method of Optimal Time Alignment (OTA) of the tsunami waveforms in the joint inversion of tsunami (DART and tide-gauges) and geodetic (GPS, InSAR, land-leveling) observations for this event. We test the application of OTA to the tsunami Green’s functions used in a previous inversion. Through a suite of synthetic tests we show that if the bias in the forward model is comprised only of delays in the tsunami signals, the OTA can correct them precisely, independently of the sensors (DART or coastal tide-gauges) and, to the first-order, of the bathymetric model used. The same suite of experiments is repeated for the real case of the 2010 Maule earthquake where, despite the results of the synthetic tests, DARTs are shown to outperform tide-gauges. This gives an indication of the relative weights to be assigned when jointly inverting the two types of data. Moreover, we show that using OTA is preferable to subjectively correcting possible time mismatch of the tsunami waveforms. The results for the source model of the Maule earthquake show that using just the first-order modeling correction introduced by OTA confirms the bilateral rupture pattern around the epicenter, and, most importantly, shifts the inferred northern patch of slip to a shallower position consistent with the slip models obtained by applying more complex physics-based corrections to the tsu

The complexity of coseismic slip distributions influences the tsunami hazard posed by local and, to a certain extent, distant tsunami sources. Large slip concentrated in shallow patches was observed in recent tsunamigenic earthquakes, possibly due to dynamic amplification near the free surface, variable frictional conditions or other factors. We propose a method for incorporating enhanced shallow slip for subduction earthquakes while preventing systematic slip excess at shallow depths over one or more seismic cycles. The method uses the classic k(-2) stochastic slip distributions, augmented by shallow slip amplification. It is necessary for deep events with lower slip to occur more often than shallow ones with amplified slip to balance the long-term cumulative slip. We evaluate the impact of this approach on tsunami hazard in the central and eastern Mediterranean Sea adopting a realistic 3D geometry for three subduction zones, by using it to model similar to 150,000 earthquakes with M-w from 6.0 to 9.0. We combine earthquake rates, depth-dependent slip distributions, tsunami modeling, and epistemic uncertainty through an ensemble modeling technique. We found that the mean hazard curves obtained with our method show enhanced probabilities for larger inundation heights as compared to the curves derived from depth-independent slip distributions. Our approach is completely general and can be applied to any subduction zone in the world.

The Alboran Basin located at the western Mediterranean hosts the boundary between the European and African plates. Previous studies support that seismicity is constrained to crustal depths. Thus, recent data acquired in the Alboran Basin during the last 15 years are key to constrain the crustal thickness as well as to characterize the sedimentary and tectonic structure of the area. The outcome of these successive surveys allows us to characterize with detail the seismic parameters of the main active faults systems (i.e. Carboneras, Al-Idrissi, Yusuf and Alboran Ridge), and to reveal their seismogenic and tsunamigenic potential. To estimate the tsunamigenic potential for each fault, we follow a two-step workflow: First, we analyse the faults in the area and generate 3D models of the fault planes. Our dense grid of multichannel seismic reflection data, together with selected pre-stack depth-migrated sections and a wide-angle seismic refraction profile showing the velocity structure across the basin, allows us to analyse for the first time in the area, the realistic geometry of the selected faults. Second, we calculated the possible tsunami scenarios for each fault system, using our 3D fault models and their simplified rectangular plane version, all together with their associated seismic parameter estimations. Tsunami simulations were perform by using the Tsunami-HySEA non-linear shallow water GPU-optimised code developed by the EDANYA Group of the University of Malaga, Spain (de la Asunción et al., 2013, Comput. Fluids). Comparing the tsunami obtained with the 3D and the rectangular models, we systematically observe distinct wave propagation patterns with enhanced wave heights. Thus, our results demonstrate the importance of fault characterization to improve the definition of the tsunamigenic potential of the fault systems of the Alboran Basin, providing an example for applications in other areas. In a further analysis, we are testing the importance of the slip distri

The Gulf of Cadiz is located in the southwestern margin of Iberia, where the slow NE-SW trending convergence between the African - Eurasian plates is of 3.8-5.6 mm/yr. Despite hosting a moderate/intermediate seismic activity, large historical and instrumental destructive earthquakes, such as the 1755 Lisbon earthquake (Mw 8.5) and tsunami, and the 1969 Horseshoe earthquake (Mw 8) occurred. On the basis of turbidite paleo-seismologic studies, we found that great earthquakes of Mw ≥ 8 occur in the external part of the Gulf of Cadiz following a recurrence interval of 1800 years. This study focuses on the active structures located in this area and their potential tsunami hazard. These structures are: a WNW-ESE dextral strike-slip fault (Lineament South) and NW-verging thrust faults such as the Marquês de Pombal, the Horseshoe Fault, the North Coral Patch Ridge Fault, the South Coral Patch Ridge Fault and Gorringe Bank Fault. Using a scenario approach on the basis of geological data acquired in the Gulf of Cadiz during the last 20 years, we are able to determine the locations of highest tsunami hazard along the coasts of the Iberian Peninsula and North Africa. The workflow involves the following tasks: 1. Interpretation of all seismic sections (in time), defining the traces of the main active faults in the area; 2. Mapping the seafloor fault trace on using multibeam bathymetry and a GIS software; 3. Mesh of the fault surface and sub-surface information, and their respective horizons in order to generate a 3D model of the subsurface for each fault; 4. Conversion of the 3D model from time to depth on the basis of available velocity models (i.e. assigning a velocity value to the intervals between horizons); 5. Defining the specific attributes for each fault, such as Length, Width, Depth, Strike, Dip and Rake. 6. Determine the maximum magnitude (Max Mw) that single faults can generate and its maximum slip, which is mechanically compatible with the length and width previously

The Alboran Basin, located in the westernmost Mediterranean, hosts the boundary between the European and African tectonic plates. The deep structure of this basin is poorly constrained by regional tomographic models, and the tsunamigenic potential of individual faults is partly undefined. Modern data acquired in the area during the last 15 years (e.g., active source seismic data from WESTMED cruise in 2006, EVENT-DEEP cruise in 2010 and TOPOMED cruise in 2011; and passive seismic data from TOPO-Europe and TOPO-Iberia projects in 2009-2010) are fundamental to constrain the crustal thickness and to image the sedimentary and tectonic structure of this area. Results of these surveys allowed us to advance in the characterization of the main active faults in the Alboran Basin (i.e., Carboneras, Al-Idrissi, Yusuf and Alboran Ridge front faults), and to review their tsunamigenic potential. To estimate the possible tsunamigenic potential for each fault, we followed a two-step workflow. First, we analysed the faults in the area, selecting the potentially most tsunamigenic ones and describing them in detail. We created 3D models of the fault planes, using bathymetric and active and passive seismic data. The dense grid of multichannel seismic data, together with the pre-stack depth-migrated sections and a refraction profile showing the velocity structure across the basin, allowed us to consider for the first time in the area realistic geometry of the major faults. Second, we calculated the possible tsunami scenarios for each fault, using our 3D fault plane models and their simplified rectangular plane version, all with associated seismic parameter estimates. Tsunami simulations were done with Tsunami-HySEA non-linear shallow water GPU-optimised code developed by the EDANYA Group of the University of Malaga, Spain (de la Asunción et al., 2013, Comput. Fluids). Main preliminary results show that: 1) the Alboran Ridge front fault expresses the highest tsunamigenic potential in the

The southwestern margin of the Iberian Peninsula, which includes the Gulf of Cadiz, is characterized by a present -day active deformation mainly driven by the NW-SE trending convergence (3.8-5.6 mm.yr-1) between the Nubia and Eurasia plates. The SW Iberian margin is a seismogenic area characterized by low to moderate magnitudes (Mw ≤ 5.5). In addition, this area hosts some of the largest earthquakes occurred in Western Europe, such as the 1st of November 1755 Lisbon Earthquake and tsunami (Mw ≥ 8.5). The active fault structures can be classified in two main families: a) WNW-ESE trending dextral strike-slip faults, and b) NE-SW trending thrusts faults. To characterize the seismogenic and tsunamigenic potential of each fault system, we develop several tsunami models. The workflow involves the following tasks: 1. Interpretation of the seismic profiles (in time) defining the traces of main active faults; 2. Mapping the trace of the faults using multibeam bathymetry; 3. Mesh of the fault surface and their respective horizons to generate a 3D model of the subsurface for each fault; 4. Conversion of the 3D model from time-to-depth assigning a velocity value (i.e. from available velocity models of the area) to the interval between horizons; 5.Defining the specific attributes for each fault, such as Length, Width, Depth, Strike, Dip and Rake; 6. Determine the maximum magnitude and slip for each fault. The maximum magnitude should be compatible with the length and the width previously defined, so we use the Leonard (2014) scaling-law; 7. Finally, the tsunami simulations for each fault have been run using “Tsunami-HySEA” software. We run two simulations for each fault, the first one considering the fault as an inclined planar surface and the second simulations used the 3D mesh

Seismological, tsunami and geodetic observations have shown that subduction zones are complex systems where the properties of earthquake rupture vary with depth as a result of different pre-stress and frictional conditions. A wealth of earthquakes of different sizes and different source features (e.g. rupture duration) can be generated in subduction zones, including tsunami earthquakes, some of which can produce extreme tsunamigenic events. Here, we offer a geological perspective principally accounting for depth-dependent frictional conditions, while adopting a simplified distribution of on-fault tectonic pre-stress. We combine a lithology-controlled, depth-dependent experimental friction law with 2D elastodynamic rupture simulations for a Tohoku-like subduction zone cross-section. Subduction zone fault rocks are dominantly incohesive and clay-rich near the surface, transitioning to cohesive and more crystalline at depth. By randomly shifting along fault dip the location of the high shear stress regions (“asperities”), moderate to great thrust earthquakes and tsunami earthquakes are produced that are quite consistent with seismological, geodetic, and tsunami observations. As an effect of depth-dependent friction in our model, slip is confined to the high stress asperity at depth; near the surface rupture is impeded by the rock-clay transition constraining slip to the clay-rich layer. However, when the high stress asperity is located in the clay-to-crystalline rock transition, great thrust earthquakes can be generated similar to the Mw 9 Tohoku (2011) earthquake.

Seismological, tsunami and geodetic observations have shown that subduction zones are complex systems where the properties of earthquake rupture vary with depth as a result of different pre-stress and frictional conditions. A wealth of earthquakes of different sizes and different source features (e.g. rupture duration) can be generated in subduction zones, including tsunami earthquakes, some of which can produce extreme tsunamigenic events. Here, we offer a geological perspective principally accounting for depth-dependent frictional conditions, while adopting a simplified distribution of on-fault tectonic pre-stress. We combine a lithology-controlled, depth-dependent experimental friction law with 2D elastodynamic rupture simulations for a Tohoku-like subduction zone cross-section. Subduction zone fault rocks are dominantly incohesive and clay-rich near the surface, transitioning to cohesive and more crystalline at depth. By randomly shifting along fault dip the location of the high shear stress regions (“asperities”), moderate to great thrust earthquakes and tsunami earthquakes are produced that are quite consistent with seismological, geodetic, and tsunami observations. As an effect of depth-dependent friction in our model, slip is confined to the high stress asperity at depth; near the surface rupture is impeded by the rock-clay transition constraining slip to the clay-rich layer. However, when the high stress asperity is located in the clay-to-crystalline rock transition, great thrust earthquakes can be generated similar to the Mw 9 Tohoku (2011) earthquake.

Seismological, tsunami and geodetic observations have shown that subduction zones are complex systems where the properties of earthquake rupture vary with depth as a result of different pre-stress and frictional conditions. A wealth of earthquakes of different sizes and different source features (e.g. rupture duration) can be generated in subduction zones, including tsunami earthquakes, some of which can produce extreme tsunamigenic events. Here, we offer a geological perspective principally accounting for depth-dependent frictional conditions, while adopting a simplified distribution of on-fault tectonic pre-stress. We combine a lithology-controlled, depth-dependent experimental friction law with 2D elastodynamic rupture simulations for a Tohoku-like subduction zone cross-section. Subduction zone fault rocks are dominantly incohesive and clay-rich near the surface, transitioning to cohesive and more crystalline at depth. By randomly shifting along fault dip the location of the high shear stress regions (“asperities”), moderate to great thrust earthquakes and tsunami earthquakes are produced that are quite consistent with seismological, geodetic, and tsunami observations. As an effect of depth-dependent friction in our model, slip is confined to the high stress asperity at depth; near the surface rupture is impeded by the rock-clay transition constraining slip to the clay-rich layer. However, when the high stress asperity is located in the clay-to-crystalline rock transition, great thrust earthquakes can be generated similar to the Mw 9 Tohoku (2011) earthquake.

Subduction zones are the favorite places to generate tsunamigenic earthquakes, where friction between oceanic and continental plates causes the occurrence of a strong seismicity. The topics and the methodologies discussed in this thesis are focussed to the understanding of the rupture process of the seismic sources of great earthquakes that generate tsunamis. The tsunamigenesis is controlled by several kinematical characteristic of the parent earthquake, as the focal mechanism, the depth of the rupture, the slip distribution along the fault area and by the mechanical properties of the source zone. Each of these factors plays a fundamental role in the tsunami generation. Therefore, inferring the source parameters of tsunamigenic earthquakes is crucial to understand the generation of the consequent tsunami and so to mitigate the risk along the coasts. The typical way to proceed when we want to gather information regarding the source process is to have recourse to the inversion of geophysical data that are available. Tsunami data, moreover, are useful to constrain the portion of the fault area that extends offshore, generally close to the trench that, on the contrary, other kinds of data are not able to constrain. In this thesis I have discussed the rupture process of some recent tsunamigenic events, as inferred by means of an inverse method. I have presented the 2003 Tokachi-Oki (Japan) earthquake (Mw 8.1). In this study the slip distribution on the fault has been inferred by inverting tsunami waveform, GPS, and bottom-pressure data. The joint inversion of tsunami and geodetic data has revealed a much better constrain for the slip distribution on the fault rather than the separate inversions of single datasets. Then we have studied the earthquake occurred on 2007 in southern Sumatra (Mw 8.4). By inverting several tsunami waveforms, both in the near and in the far field, we have determined the slip distribution and the mean rupture velocity along the causative fault. S

Subduction zones are the favorite places to generate tsunamigenic earthquakes, where friction between oceanic and continental plates causes the occurrence of a strong seismicity. The topics and the methodologies discussed in this thesis are focussed to the understanding of the rupture process of the seismic sources of great earthquakes that generate tsunamis. The tsunamigenesis is controlled by several kinematical characteristic of the parent earthquake, as the focal mechanism, the depth of the rupture, the slip distribution along the fault area and by the mechanical properties of the source zone. Each of these factors plays a fundamental role in the tsunami generation. Therefore, inferring the source parameters of tsunamigenic earthquakes is crucial to understand the generation of the consequent tsunami and so to mitigate the risk along the coasts. The typical way to proceed when we want to gather information regarding the source process is to have recourse to the inversion of geophysical data that are available. Tsunami data, moreover, are useful to constrain the portion of the fault area that extends offshore, generally close to the trench that, on the contrary, other kinds of data are not able to constrain. In this thesis I have discussed the rupture process of some recent tsunamigenic events, as inferred by means of an inverse method. I have presented the 2003 Tokachi-Oki (Japan) earthquake (Mw 8.1). In this study the slip distribution on the fault has been inferred by inverting tsunami waveform, GPS, and bottom-pressure data. The joint inversion of tsunami and geodetic data has revealed a much better constrain for the slip distribution on the fault rather than the separate inversions of single datasets. Then we have studied the earthquake occurred on 2007 in southern Sumatra (Mw 8.4). By inverting several tsunami waveforms, both in the near and in the far field, we have determined the slip distribution and the mean rupture velocity along the causative fault. S