Finite-frequency Rayleigh wave tomography of the western Mediterranean: Mapping its lithospheric structure

We present a 3-D shear wave velocity model for the crust and upper mantle of the western Mediterranean from Rayleigh wave tomography. We analyzed the fundamental mode in the 20¿167 s period band (6.0¿50.0 mHz) from earthquakes recorded by a number of temporary and permanent seismograph arrays. Using the two-plane wave method, we obtained phase velocity dispersion curves that were inverted for an isotropic Vs model that extends from the southern Iberian Massif, across the Gibraltar Arc and the Atlas mountains to the Saharan Craton. The area of the western Mediterranean that we have studied has been the site of complex subduction, slab rollback, and simultaneous compression and extension during African-European convergence since the Oligocene. The shear velocity model shows high velocities beneath the Rif from 65 km depth and beneath the Granada Basin from ¿¿70 km depth that extend beneath the Alboran Domain to more than 250 km depth, which we interpret as a near-vertical slab dangling from beneath the western Alboran Sea. The slab appears to be attached to the crust beneath the Rif and possibly beneath the Granada Basin and Sierra Nevada where low shear velocities (3.8 km/s) are mapped to >55 km depth. The attached slab is pulling down the Gibraltar Arc crust, thickening it, and removing the continental margin lithospheric mantle beneath both Iberia and Morocco as it descends into the deeper mantle. Thin lithosphere is indicated by very low upper mantle velocities beneath the Alboran Sea, above and east of the dangling slab and beneath the Cenozoic volcanics.This research was funded by the U.S. National Science Foundation EAR-0808939. The deployment of the IberArray broadband seismic network is part of the CONSOLIDER CSD2006-00041 (Geosciences in Iberia: Integrated studies on Topography and 4-D Evolution) grant from the Spanish Ministry of Science and Innovation. Additional funding was provided by the Spanish ministry under grants CGL2010-17280 and by Generalitat de Catalunya under grant 2009 SGR 6.Peer Reviewe

Mapping the onshore/offshore crustal transition at the Westernmost Mediterranean from seismic profiling

Cembrowski, Marcel... et. al.-- European Geosciences Union General Assembly 2013, 7-12 April, Vienna, Austria.-- 1 pageThe evolution of theWesternMediterranean is strongly affected by the collision of the African and Eurasian plates. The plate boundary as seen from earthquakes is diffuse over a wide area extending north and south of Gibraltar strait. The Western end of the Mediterranean is delineated by the Gibraltar Arc System, comprising the arcuate Spanish Betic and Moroccan Rif Mountain Belts, together with the Alboran-Sea Basin in-between. The extension of the Alboran Basin which started from Late Eocene and which coexisted with the Africa -Europe conversion is still under debate and is one of the key points to constrain the evolution of the Western Mediterranean. This motivated our interest to map the still unknown crustal transition from the Betic-Rif chain into the Alboran Sea, taking advantage of the coincidence in time (October 2011) of two seismic experiments in the area, on land (Rifsis project) and at sea (Gassis-WestMed project). For this purpose we deployed several tens of seismic stations, both in Morocco and Spain, to record the air-gun shooting (every 50 m) of the Sarmiento de Gamboa Spanish vessel performing multichannel reflection profiles at the Alboran sea, and hence to extend these marine lines to wide-angle distances in-land. The airguns were calibrated for the near zero-offset marine reflection study and it turns out to be difficult to observe clear signals on the records in-land at offsets larger than about 70 km. The data has therefore been processed with a frequency-dependent lateral coherence filter to enhance coherent reflection/refraction signals through the frequency-dependent attenuation of incoherent noise and signals. This processing has permitted to track signals (seismic energy) up to more than 200 km on some profiles. Hence, a classical procedure of forward modeling (ray tracing approach) to fit the travel times of the identified wide-angle phases is now underway, taking advantage of the sedimentary/basement sequences inferred from the multichannel sections to constrain the upper part of the velocity-depth model. The first structural results delineate significant lateral variations in crustal depths, particularly at the Rif-Alboran transition. In our presentation we will show and discuss data processing examples which enabled signal detection to large offsets, the signal identification and their interpretation, and the different 2-D cross sections which image the crustal transition to the Alboran BasinPeer Reviewe

Seismic structure and composition of the southern central Iberian crust: The ALCUDIA wide angle seismic reflection transect

The authors thank the Associated Editor and the anonymous reviewers for their thoroughly valuable suggestions and comments that improved the manuscript. Seismic data were collected in 2012 with funding provided by the Spanish Ministry of Science and Innovation (grants: CGL2004-04623/BTE, CGL2007-63101/BTE, CGL2011-24101, CSD2006-00041). Instrumentation was provided by the IRIS-PASSCAL instrument center, Socorro, New Mexico, USA. The seismic data, including experiment geometry are stored in the IRIS-PASSCAL facilities and can be accessed through the IRIS-PASSCAL data management center. I.P. is funded by the Spanish Government and the Universidad de Salamanca with a Beatriz Galindo grant (BEGAL 18/00090). S.A. Ehsan is funded by the European Commission grant Marie Curie Actions (264517-TOPOMOD-FP7-PEOPLE-2010-ITN). We thank Instituto Geologico y Minero de Espana for providing the logistic help and an academic crew for data acquisition. GMT was used to prepare some of the figures shown in the paper.The nature of the crust beneath central Iberia was estimated by a wide-angle seismic reflection/refraction transect, ALCUDIA-WA, which sampled the southern half of the Variscan Central Iberian Zone, covered in the north by the Cenozoic Tajo Basin. The shot gathers recorded by vertical component sensors revealed well defined P- and S-wave phases. These arrivals were modeled by an iterative forward approach providing 2D crustal models showing variations in the velocity distribution with upper crustal P- and S-wave velocities increasing northwards. The lower crust P-wave velocities are homogeneous along the profile while the S-wave velocities slightly increase northwards. The Moho is placed at 32 km depth in the southern edge of the profile, deepening northward down to 35 km beneath the Tajo Basin. The Poisson's ratio, calculated from P- and S-wave velocities, varies along the profile at upper crustal depths. The highest values are located below the Mora and Pedroches batholiths. These resulting physical properties can serve to constrain the crustal composition by comparing them with laboratory measurements on rock samples. Our results suggest that the upper crust in the southern and central segments of the ALCUDIA profile is made up of low-grade metasedimentary rocks, while the northern segment is dominated by igneous rocks, in agreement with the surface geology. Separated by a sharp boundary located between 12 km (south) and 18 km (north) depth, the lower crust is more homogeneous and shows low Poisson’ ratios compatible with a rather felsic composition. However, outstanding lamination described in coincident vertical incidence data indicates some degree of intercalation with mafic components.Spanish Government CGL2004-04623/BTE CGL2007-63101/BTE CGL2011-24101 CSD2006-00041Spanish GovernmentEuropean Commission BEGAL 18/00090Universidad de Salamanca BEGAL 18/00090European Commission grant Marie Curie Actions 264517-TOPOMOD-FP7-PEOPLE-2010-IT

Subduction-driven recycling of continental margin lithosphere

Whereas subduction recycling of oceanic lithosphere is one of the central themes of plate tectonics, the recycling of continental lithosphere appears to be far more complicated and less well understood(1). Delamination and convective downwelling are two widely recognized processes invoked to explain the removal of lithospheric mantle under or adjacent to orogenic belts(2-5). Here we relate oceanic plate subduction to removal of adjacent continental lithosphere in certain plate tectonic settings. We have developed teleseismic body wave images from dense broadband seismic experiments that show higher than expected volumes of anomalously fast mantle associated with the subducted Atlantic slab under northeastern South America and the Alboran slab beneath the Gibraltar arc region(6,7); the anomalies are under, and are aligned with, the continental margins at depths greater than 200 kilometres. Rayleigh wave analysis(8,9) finds that the lithospheric mantle under the continental margins is significantly thinner than expected, and that thin lithosphere extends from the orogens adjacent to the subduction zones inland to the edges of nearby cratonic cores. Taking these data together, here we describe a process that can lead to the loss of continental lithosphere adjacent to a subduction zone. Subducting oceanic plates can viscously entrain and remove the bottom of the continental thermal boundary layer lithosphere from adjacent continental margins. This drives surface tectonics and pre-conditions the margins for further deformation by creating topography along the lithosphere-asthenosphere boundary. This can lead to development of secondary downwellings under the continental interior, probably under both South America and the Gibraltar arc(8,10), and to delamination of the entire lithospheric mantle, as around the Gibraltar arc(11). This process reconciles numerous, sometimes mutually exclusive, geodynamic models proposed to explain the complex oceanic-continental tectonics of these subduction zones(12-17)This research was supported by US National Science Foundation grants EAR 0003572, 0607801 and 0808939 (A.L.), EAR 0808931 (E.D.H.), EAR 0809023 and 1054638 (M.S.M.), the Venezuelan National Fund for Science, Technology and Innovation grant G-2002000478 and PDVSA-INTEVEP-FUNVISIS cooperative agreement 2004-141 (M.S.), the Spanish Ministry of Science and Innovation grants CSD2006-00041, CGL2009-09727 and CGL2010-15146 (J.G. and R.C.), and by an A. v. Humboldt Foundation Research Prize (A.L.).Peer Reviewe

The Crust beneath Morocco: From the surface topography to the upper mantle a 700 km long seismic section across Morocco.

The most characteristic topographic features of Morocco are the Atlas Mountains and the Rif Coordillera. These two orogenic belts are the response of different geodynamic processes acting at lithospheric scale caused by a complex plate interaction. Both are located within the diffuse plate boundary zone separating Africa and Europe. The boundary zone is characterized by a relatively broad zone of deformation that includes mountain chains in southern Iberia, the Betics and in Morocco, the Rif Cordillera, separated by the Alboran basin. The zone delineates an arcuate arc system known as the Gibraltar arc. The area is characterized by a relatively large amount of earthquake activity at various depths and with a broad spectra of focal mechanisms. Within the last decade a large international effort have been devoted to the area. The topic has fostered a strong collaborations between Spanish and international research teams form Europe and USA. Key multi-seismic projects have been developed that aim to constrain the structure, composition and tectonic scenario from south of the Atlas to the Betics, across the Rif cordillera and the Alboran basin. The multidisciplinary research program included: natural source (earthquakes) recording with temporal deployments of broad band (BB) instrumentation and, controlled source seismic acquisition experiments where, spatially dense recording of wide-angle seismic reflection shot gathers were acquired. The natural source experiments consisted on a transect from Merzouga across the Gibraltar Arc and into the Iberian Peninsula (until south of Toledo) and, a nearly regular grid of BB. The controlled source data-sets were able to constrain the crustal structure and provide seismic P-wave propagation velocity models from the coast across the Rif and the Atlas. From south to north the crust features a relatively moderate crustal root beneath the Middle Atlas which can reach 40 km clearly differing from the 35 km thickness value observed at both sides of this root. Travel time inversion results position the crustal root just south of the High Atlas defining a thrusted mantle wedge and, also a limited crustal imbrication is suggested in the Middle Atlas. The most surprising feature is a prominent and unexpected crustal root (over 50 km) located beneath the external Rif and identified by both the wide-angle data and receiver function studies. To the east of this feature the crust thins rapidly by 20 km across the Nekkor fault zone, suggested to be related to the sharp change in crustal thickness. On shore-offshore recording of marine shots reveal further complexities in the transition to the Alboran basin. The low values of the Bouguer gravity anomalies beneath the Rif Cordillera are consistent with the crustal models derived from the new seismic data. The detailed knowledge on the crustal structure achieved by this high resolution imaging geophysical techniques is an asset to evaluate the earthquake and potential tsunami hazard for the coasts of North Africa and western Europe.This work has been primarily funded by the Spanish MEC project CGL2007–63889. Additional funding was provided by projects CGL2010–15416, CSD2006-00041, and CGL2009–09727 (Spain), CGL2008–03474-E, 07- TOPO_EUROPE_FP-006 (ESF Eurocores) and EAR-0808939 (US, NSF).Peer Reviewe

Characterization of the shallow subsurface structure across the Carrascoy Fault System (SE Iberian Peninsula) using P-wave tomography and Multichannel Analysis of Surface Waves

Acknowledgement. The authors would like to acknowledge the project INTERGEO (CGL2013-47412-C2-1-P) GEO3BCNCSIC for the data access. Data are public access through SeisDARE (DeFelipe et al., 2021), dataset Martí et al.(2015). The Ministry of Education and Culture of the Republic of Indonesia is thanked for the main author's Ph.D. scholarship (D3.2/KD.02.01/2019). JA is funded by MICINN (IJC2018-026335-I). I.P. is funded by the Spanish Government and the Universidad de Salamanca (Beatriz Galindo grant BEGAL 18/00090). IDF is funded by a FEDER-Junta de Castilla y León Postdoctoral contract (SA0084P20). We thank the GIPP-GFZ, (Germany) and Lisbon University (Portugal) for the instrumentation provided. Generalitat de Catalunya (AGAUR) grant 2017SGR1022 (GREG); EU (H2020) 871121 (EPOSSP); and EIT-RawMaterias 17024 (SIT4ME). We sincerely thank Seismic Unix CWP (Center for the Wave Phenomena, Colorado School of Mines (Cohen and Stockwell, 2019)). We also thank all the people involved directly or indirectly in this work.The seismicity in the SE Iberian Peninsula is distributed parallel to the coast in a well-developed strike-slip fracture system known as the Eastern Betic Shear Zone (EBSZ). This work focuses on the characterization of the shallow subsurface structure of the Algezares-Casas Nuevas Fault, within the Carrascoy Fault System of the EBSZ. The Carrascoy Fault borders the Guadalentín Depression to the south, which is a densely populated area with extensive agricultural activity. Therefore, this faults system represents a seismic hazard with significant social and economic implications. We have constructed two velocity-depth models based on P-wave tomography and Multichannel Analysis of Surface Waves (MASW) acquired from seismic reflection data. The resulting velocity models have allowed us to interpret the first ~250m depth and have revealed: i) the thickness of the critical zone; ii) the geometry of the Algezares-Casas Nuevas Fault; iii) the depth of the Messinian/Tortonian contact and iv) the presence of blind thrusts and damage zones under the Guadalentín Depression. Our results have also helped us to estimate an apparent vertical slip rate of 0.66±0.06m/ky for the Algezares-Casas Nuevas Fault since 209.1±6.2ka. Our results provide a methodological and backflow protocol to study the shallow subsurface of active faults, complementing previous geological models based on paleoseismological trenches, and can be used to improve the seismic hazard assessment of tectonically active regions around the world

From the Surface Topography to the Upper Mantle, Seismic constraints on the Crustal structure Across Morocco

The most characteristic topographic features of Morocco are the Atlas Mountains and the Rif Coordillera. These two orogenic belts are the response of different geodynamic processes acting at lithospheric scale caused by a unique driver, the collision between two tectonic plates. Both are located within the diffuse plate boundary zone separating Africa and Europe. The boundary zone is characterized by a relatively broad zone of deformation that includes Mountain chains in southern Iberia, the Betics and the Rif cordillera in Morocco. The zone delineates the arcuate arc system of Gribraltar. Within the last decade a large international effort have been devoted to the area mostly leaded by Spanish groups with the collaboration of international research teams (including scientist form Europe and USA). Key multi-seismic projects have been developed that aim to constrain the structure, composition and tectonic scenario from south of the Atlas to the Betics, across the Rif cordillera and the Alboran basin. The multidisciplinary research program includes: natural source (earthquakes) recording with temporal deployments of broad band (BB) instrumentation and, controlled source seismic acquisition experiments where, spatially dense recording of wide-angle seismic reflection shot gathers were acquired. The passive experiments consisted on: a transect from Merzouga across the Gibraltar arc and into the Iberian peninsula (untill south of Toledo); a nearly regular grid of BB which was achieve by multiple deployments of a number of BB. The controlled source datasets were able to constrain the crustal structure and provide seismic P-wave propagation velocity models from the coast across the Rif and the Atlas. Travel-time inversion of the controlled source seismic data across the Atlas constraints a crustal root to the south of the High Atlas, and reveals mantle wedge. A limited crustal imbrication also appears in the Middle Atlas. The crustal thickness, does not exceeded 40 km in the root zone and is less than 35 km elsewhere. Approaching the Rif, the controlled source data reveals a prominent crustal root over 50 km depth which is located where the topography does not exceed 1400 m. These features indicate that complex structure and processes beneath the crust play a key role in supporting the particular geometry of the surface topography of this part of the western Mediterranean. On one hand the Atlas is being supported by the mantle, on the other the abrupt change in crustal thickness at the Nekkor fault and the deep Rif crustal root can be attributed to interaction of the subducting Alboran slab with the North African Neo-Tethys passive margin.Peer Reviewe

Characterization of the shallow subsurface structure across the Carrascoy Fault System (SE Iberian Peninsula) using P-wave tomography and Multichannel Analysis of Surface Waves

The seismicity in the SE Iberian Peninsula is distributed parallel to the coast in a well-developed strike-slip fracture system known as the Eastern Betic Shear Zone (EBSZ). This work focuses on the characterization of the shallow subsurface structure of the Algezares-Casas Nuevas Fault, within the Carrascoy Fault System of the EBSZ. The Carrascoy Fault borders the Guadalentín Depression to the south, which is a densely populated area with extensive agricultural activity. Therefore, this faults system represents a seismic hazard with significant social and economic implications. We have constructed two velocity-depth models based on P-wave tomography and Multichannel Analysis of Surface Waves (MASW) acquired from seismic reflection data. The resulting velocity models have allowed us to interpret the first ~250m depth and have revealed: i) the thickness of the critical zone; ii) the geometry of the Algezares-Casas Nuevas Fault; iii) the depth of the Messinian/Tortonian contact and iv) the presence of blind thrusts and damage zones under the Guadalentín Depression. Our results have also helped us to estimate an apparent vertical slip rate of 0.66±0.06m/ky for the Algezares-Casas Nuevas Fault since 209.1±6.2ka. Our results provide a methodological and backflow protocol to study the shallow subsurface of active faults, complementing previous geological models based on paleoseismological trenches, and can be used to improve the seismic hazard assessment of tectonically active regions around the world.The authors would like to acknowledge the project INTERGEO (CGL2013-47412-C2-1-P) GEO3BCNCSIC for the data access. Data are public access through SeisDARE (DeFelipe et al., 2021), dataset Martí et al. (2015). The Ministry of Education and Culture of the Republic of Indonesia is thanked for the main author’s Ph.D. scholarship (D3.2/KD.02.01/2019). JA is funded by MICINN (IJC2018-026335-I). I.P. is funded by the Spanish Government and the Universidad de Salamanca (Beatriz Galindo grant BEGAL 18/00090). IDF is funded by a FEDER-Junta de Castilla y León Postdoctoral contract (SA0084P20). We thank the GIPP-GFZ, (Germany) and Lisbon University (Portugal) for the instrumentation provided. Generalitat de Catalunya (AGAUR) grant 2017SGR1022 (GREG); EU (H2020) 871121 (EPOSSP); and EIT-RawMaterias 17024 (SIT4ME). WPeer reviewe

Crustal Imbrication in an Alpine Intraplate Mountain Range: A Wide-Angle Cross-Section Across the Spanish-Portuguese Central System

Intraplate ranges are topographic features that can occur far from plate boundaries, the expected position of orogens as described in the plate tectonics theory. To understand the lithospheric structure of intraplate ranges, we focused on the Spanish-Portuguese Central System (SPCS), the most outstanding topographic feature in the central Iberian Peninsula. The SPCS is an Alpine range that exhumes Precambrian-Paleozoic rocks and is located at >200 km from the northern border of the Iberian microplate. Here, we provide a P-wave velocity model based on wide-angle seismic reflection/refraction data of the central SPCS (Gredos sector). Our results show: (a) a layered lithosphere characterized by three major interfaces: Conrad, Mohorovicic, and Hales discontinuities, (b) an asymmetry of the crust-mantle boundary under the SPCS, (c) the extent of the Variscan batholith forming the main outcrops of Gredos, and (d) the thinning of the lower crust toward the south. This model suggests that the exhumation of the SPCS basement was driven by a south-vergent thick-skinned thrust system, developed in the southern part of the SPCS and that promoted crustal imbrication and a Mohorovicic discontinuity's offset under the SPCS. Thus, the deformation mechanisms of the crust seem to be controlled by the presence of the late- to post-Variscan granitoids that assimilated the Variscan mid-crustal detachment creating a new rheological boundary. This tectonic structure allowed the formation of Alpine crustal-scale thrust systems that eased coupled deformation of the upper and lower crust, leading to limited underthrusting of both crustal layers.This study has been funded by the Ministry of Science, Innovation and Competitiveness through the Project CIMDEF (CGL2014-56548-P). IP is funded by the Spanish Government and the University of Salamanca (Beatriz Galindo grant BEGAL 18/00090). JA is funded by grant IJC2018-036074-I, funded by MCIN/AEI/10.13039/501100011033. DMP and FGL are also funded by grants CGL2015-71692 (MINECO/ FEDER) and PID2020-118822GB-I00 (MCIN/AEI/10.13039/501100011033).Peer reviewe

Crustal Imbrication in an Alpine Intraplate Mountain Range: A Wide-Angle Cross-Section Across the Spanish-Portuguese Central System

Intraplate ranges are topographic features that can occur far from plate boundaries, the expected position of orogens as described in the plate tectonics theory. To understand the lithospheric structure of intraplate ranges, we focused on the Spanish-Portuguese Central System (SPCS), the most outstanding topographic feature in the central Iberian Peninsula. The SPCS is an Alpine range that exhumes Precambrian-Paleozoic rocks and is located at >200 km from the northern border of the Iberian microplate. Here, we provide a P-wave velocity model based on wide-angle seismic reflection/refraction data of the central SPCS (Gredos sector). Our results show: (a) a layered lithosphere characterized by three major interfaces: Conrad, Mohorovicic, and Hales discontinuities, (b) an asymmetry of the crust-mantle boundary under the SPCS, (c) the extent of the Variscan batholith forming the main outcrops of Gredos, and (d) the thinning of the lower crust toward the south. This model suggests that the exhumation of the SPCS basement was driven by a south-vergent thick-skinned thrust system, developed in the southern part of the SPCS and that promoted crustal imbrication and a Mohorovicic discontinuity's offset under the SPCS. Thus, the deformation mechanisms of the crust seem to be controlled by the presence of the late- to post-Variscan granitoids that assimilated the Variscan mid-crustal detachment creating a new rheological boundary. This tectonic structure allowed the formation of Alpine crustal-scale thrust systems that eased coupled deformation of the upper and lower crust, leading to limited underthrusting of both crustal layers.Ministry of Science, Innovation and Competitiveness through the Project CIMDEF (CGL2014-56548-P)Spanish Government and the University of Salamanca (Beatriz Galindo grant BEGAL 18/00090)Grant IJC2018-036074-I, funded by MCIN/AEI/10.13039/501100011033Grants CGL2015-71692 (MINECO/ FEDER) and PID2020-118822GB-I00 (MCIN/AEI/10.13039501100011033