L Buontempo

G H R Bokelmann
email: bokelmann@gm.univ-montp2.fr

G Barruol
email: barruol@gm.univ-montp2.fr

J Morales
email: morales@iag.ugr.es

Seismic anisotropy beneath southern Iberia from SKS splitting

Seismic anisotropy of the south Iberian upper mantle is investigated using shear-wave splitting of SKS phases. We analyzed teleseismic events recorded by sixteen permanent broadband stations installed on the southern Iberian Peninsula and in northern Africa, and we determined fast polarization directions ϕ, and delay times δt between fast and slow components. The area of investigation extends across two important geological structures in the Variscan Iberian Peninsula: the Variscan Iberian Massif in its center, and the Gibraltar arc in the Southeast, that represents the most westerly Alpine belt in the western Mediterranean. Shear-wave splitting measurements from stations in the Betic domain show homogeneous ENE-WSW fast directions nearly parallel to the trend of the mountain belt, and smooth spatial variations. Stations in the North, toward the southern part of the Variscan Iberian Massif show homogeneous fast directions however trending NS to NE-SW, different from those recorded in the Betic. These observations may reflect a post-Hercynian (Variscan) deformation of the Ossa-Morena zone, related to the main stages in the tectonic evolution of this part, namely transpressional stage, transtensional stage and shortening episode, or a deformation related to the posterior Alpine orogeny. Along the Gibraltar arc, we observe a smoothly varying ϕ trend changing from ENE-WSW in the Eastern Betics to NS in the area of Gibraltar and Ceuta, following more or less the general trend of the mountain belt around the Alboran Sea, and the coastline. Since a similar rotation is also visible in results from Pn anisotropy, this suggests that the anisotropy is vertically coherent starting from just below the Moho. Comparing the anisotropy pattern expected from various geodynamic models with the observed SKS splitting suggests that the anisotropy is best explained by a model of slab rollback, rather than by delamination models.

Introduction

During the last decades, seismological techniques have been developed to analyze the deep deformation of the Earth, and particularly the relation between plate tectonics and mantle flow as well as the internal deformation of the plate (e.g., [START_REF] Silver | Implications for continental structure and evolution from seismic anisotropy[END_REF]. The most popular technique uses the birefringence in anisotropic media that induces the splitting of shear waves and allows one to obtain quantitative information about seismic anisotropy in the upper mantle below a seismic station. Seismic anisotropy in the upper mantle is dominated by elastic anisotropy of rockforming minerals which develop lattice-preferred orientations (LPO) or fabrics in response to strain (e.g. [START_REF] Nicolas | Formation of anisotropy in upper mantle peridotites-A review[END_REF][START_REF] Mainprice | The seismic anisotropy of the Earth's mantle: from single crystal to polycrystal[END_REF]. The primary constituent of the upper mantle is olivine which is more anisotropic than most other minerals. Olivine deforms by dislocation-creep at upper mantle conditions [START_REF] Nicolas | Formation of anisotropy in upper mantle peridotites-A review[END_REF], and it generally aligns the main crystallographic axes with respect to the structural directions (lineation, pole of the foliation), therefore producing a large-scale anisotropy that is detectable by seismic waves. The anisotropy is thus correlated with the strain in the upper mantle, and it therefore allows constraining deformation at that depth within the Earth remotely from the surface.

The most useful seismological phase for studying shear-wave splitting is the SKS phase. It travels through the mantle as an S wave, through the liquid outer core as a P wave after a S to P conversion at the core-mantle boundary, and is converted back into a radially polarized S wave when passing back in the mantle. On its way to the Earth's surface, this SKS phase is split into two mutually orthogonally polarized shear waves as it propagates through an anisotropic medium. From three-component seismic records, two parameters can be measured to quantify the anisotropy: 1) the difference in arrival time (δt) between the two split waves that depends on thickness and the intrinsic anisotropy of the medium, and 2) the polarization plane azimuth ϕ of the faster split shear wave that is related to the orientation of the anisotropic structure. Shear-wave splitting measurements can therefore characterize the orientation and to a lesser degree the depth extent of the mantle strain field (see review by [START_REF] Savage | Seismic anisotropy and mantle deformation: what have we learned from shear wave splitting?[END_REF].

Determining the vertical location of the anisotropic layer(s) is always challenging since the splitting of SKS phases may occur somewhere between the core-mantle boundary and the station.

Other phases such as local S (e.g. [START_REF] Bowman | Shear-wave splitting in the upper-mantle wedge above the Tonga subduction zone[END_REF], may help to constrain the anisotropy distribution beneath a station, but they require the presence of nearby seismogenic zones. Petrophysical analyses of mantle minerals and rocks [START_REF] Mainprice | The seismic anisotropy of the Earth's mantle: from single crystal to polycrystal[END_REF][START_REF] Mainprice | Pressure sensitivity of olivine slip systems: implications for the interpretation of seismic anisotropy of the Earth's upper mantle[END_REF] suggest that anisotropy lies in the uppermost 400 km of the Earth. Recent theoretical studies using finite-frequency sensitivity kernels provide similar conclusions since SKS splitting are shown to be more sensitive to the upper mantle above the transition zone than to the lower mantle [START_REF] Sieminski | Finite-frequency sensitivity of body waves to anisotropy based upon adjoint methods[END_REF]. If one assumes that the crust may contribute only a few tenths of a second to the total SKS splitting [START_REF] Barruol | A quantitative evaluation of the contribution of crustal rocks to the shear wave splitting of teleseismic SKS waves[END_REF][START_REF] Godfrey | Anisotropy of schists: contribution of crustal anisotropy to active source seismic experiments and shear wave splitting observations[END_REF], these arguments suggest that the anisotropy lies primarily in the subcrustal upper mantle. It is however not clear whether the anisotropy lies within the lithosphere and is related to a 'frozen' deformation [START_REF] Silver | Implications for continental structure and evolution from seismic anisotropy[END_REF], or within the asthenosphere and is thus related to present-day plate motion [START_REF] Vinnik | Global patterns of azimuthal anisotropy and deformations in the continental mantle[END_REF]. This question can be addressed by SKS splitting only indirectly since these observations constrain only weakly the depth of anisotropy within the upper mantle. On the other hand, this question can in principle be resolved combining SKS observations with surface wave anisotropy and Pn observations. We will in fact use Pn in this study to compare with SKS splitting results.

The area of interest of this paper, in southern Iberia and around the Gibraltar arc, represents a key area of geodynamic interest for the Africa-Eurasia collision and for the various stages of opening and closure of the Mediterranean basins. Various geodynamic models and mechanisms have been proposed to explain the large-scale structures and timing constraints, based on regional-scale recycling of lithosphere into the mantle, such as active continental subduction [START_REF] Morales | Active continental subduction beneath the Betic Cordillera and the Alborán Sea[END_REF], delamination (Seber et al., 1996a) or active oceanic subduction [START_REF] Lonergan | Origin of the Betic-Rif mountain belt[END_REF][START_REF] Gutscher | Evidence for active subduction beneath Gibraltar[END_REF].

In this study, we present new shear-wave splitting measurements obtained at permanent stations located in the Betic and Rif Cordillera and in the southern part of the Variscan Iberian Massif. We compare these new results with previous ones from Iberia [START_REF] Diaz | Anisotropy beneath the Iberian Peninsula: the contribution of the ILIHA-NARS broad-band experiment[END_REF][START_REF] Schmid | Delay times and shear wave splitting in the Mediterranean region[END_REF], as well as from Pn phase anisotropy (Calvert et al., 2000b). In the light of tomographic models for this region we propose that all these anisotropy observations may help to constrain and to better understand the geodynamic setting of this complex area.

Geodynamic setting

The Betic and Rif belts, connected through the Gibraltar arc, represent the most westerly Alpine range in southern Europe (Fig. 1). Together with the Alboran basin, these belts constitute one of the most complex and most debated geological regions in the western Mediterranean. A major characteristic of the Alboran basin is that it is an extensional regime [START_REF] Comas | The origin and tectonic history of the Alborán Basin: insights from Leg 161 results[END_REF] within a general collision regime, induced by the convergence of Africa and Eurasia from the Cretaceous to Paleogene [START_REF] Dewey | Kinematics of the western Mediterranean, in: Conference on Alpine tectonics[END_REF]. The geological records and geophysical observations shed some light on the evolution of the region. Formally, the Gibraltar arc is divided into three pre-Neogene crustal domains: the first consists of Mesozoic and Tertiary sedimentary rocks deformed during the Neogene as a thrust-and-fold belt mainly directed to the NW and WNW (e.g., [START_REF] Garcia-Dueñas | Les unités allochtones de la zone Subbétique,dans la transversale de Grenade (Cordillères Bétiques, Espagne)[END_REF][START_REF] Garcia-Hernandez | Mesozoic paleogeographic evolution of the External Zones of the Betic Cordillera[END_REF], which covered the South-Iberian and Maghrebian continental paleomargins and constitute the External Zones of the Betics and Rif. The second corresponds to the sediments deposited on oceanic troughs or thin continental crust [START_REF] Luján | The Flysch Trough thrust imbricate (Betic Cordillera):a key element of the Gibraltar Arc orogenic wedge[END_REF] and belongs to the so-called Flysch Trough that is well-represented in the Rif and Tell regions. The last one, the Alboran crustal domain, represents the internal zones of the orogen, composed mainly of Paleozoic to Mesozoic rocks affected by low to very high-grade metamorphism that were deformed mainly during the late Cretaceous to Paleogene and thrust onto the south Iberian and north African margins during the Early Miocene [START_REF] Balanyá | El Cabalgamiento cortical de Gibraltar y la tectónica de Béticas y Rif, presentado en II Congreso Geológico de España[END_REF]. Three tectono-metamorphic complexes are recognized in this zone, the Nevado-Filabride, the Alpujarride and the Malaguide in the Betics. The Alboran domain shows N-S continuity below the sea and constitutes the metamorphic basement of the Alboran basin (Platt et al., 1996). That basin represents a post-orogenic back arc basin [START_REF] Lonergan | Origin of the Betic-Rif mountain belt[END_REF], similar to the Calabrian arc [START_REF] Faccenna | Lateral slab deformation and the origin of the western Mediterranean arcs[END_REF][START_REF] Rosenbaum | Neogene and Quaternary rollback evolution of the Tyrrhenian Sea, the Apennines and the Sicilian Maghrebides[END_REF], formed during the westward migration of the orogenic wedge. The main extension occurred during the Early Miocene [START_REF] Comas | The origin and tectonic history of the Alborán Basin: insights from Leg 161 results[END_REF], coeval with shortening in the orogen. The extension was associated with subsidence during the Miocene and accompanied by tholeitic and calcoalkaline magmatism [START_REF] Turner | Magmatism associated with orogenic collapse of the Betic-Alboran Domain[END_REF][START_REF] Duggen | Post-collisional transition from subduction-to intraplate-type magmatism in the westernmost Mediterranean: evidence for continental-edge delamination of subcontinental lithosphere[END_REF]. Since the Upper Miocene-Pliocene until present, a contractive episode with a roughly NW-SE direction, although varying from the Upper Miocene to the Pliocene [START_REF] Ott D'estevou | Evolution structurale de la zone bétique orientale (Espagne) du Tortonian a l'Holocene[END_REF][START_REF] De Larouzière | The Betic segment of the lithospheric Trans-Alboran shear zone during the Late Miocene[END_REF][START_REF] Galindo-Zaldívar | Crustal structure of the central sector of the Betic Cordillera (SE Spain)[END_REF] manifests itself by strike-slip tectonics, mainly in the eastern Betics and Alboran basin. This contractive episode was accompanied by magmatism of mainly shoshonitic, lamproitic and alkaline basalts [START_REF] Lonergan | Origin of the Betic-Rif mountain belt[END_REF][START_REF] Zeck | Post-collisional volcanism in a sinking slab setting-crustal anatectic origin of pyroxene-andesite magma, Caldear Volcanic Group, Neogene Alborán volcanic Province, southeastern Spain[END_REF][START_REF] Duggen | Post-collisional transition from subduction-to intraplate-type magmatism in the westernmost Mediterranean: evidence for continental-edge delamination of subcontinental lithosphere[END_REF].

At present, the kinematics of the Gibraltar arc as defined by GPS and seismic moment tensor observations shows different deformation styles that are characterized by NNW-directed thrust faulting near the Algerian coast and the SW of Iberia and an E-W extensional regime in the Alboran Sea [START_REF] Stich | Moment tensor solutions for small and moderate earthquakes in the Ibero-Maghreb region[END_REF][START_REF] Stich | Kinematics of the Iberia-Maghreb plate contact from seismic moment tensors and GPS observations[END_REF]. A crucial role is played by the so-called Trans-Alboran shear zone, a major structure formed by NE-SW sinistral strike-slip faults, and NW-SE conjugates in the partitioning of the strain in the region [START_REF] Fernández-Ibáñez | Present-day stress field in the Gibraltar Arc (western Mediterranean)[END_REF] that allow coeval shortening and extensional structures in a oblique convergent plate margin between Africa and Europe.

Several models have been proposed to explain the formation of the Gibraltar arc. These models may be compared by considering the style of recycling of lithospheric material back into the mantle that may occur in the form of lithospheric delamination (e.g. Seber et al., 1996a;[START_REF] Platt | Thermal evolution, rate of exhumation, and tectonic significance of metamorphic rocks from the floor of the Alboran extensional basin, western Mediterranean[END_REF]Calvert et al., 2000b), convective removal [START_REF] Platt | Extensional collapse of thickened continental lithosphere: a working hypothesis for the Alboran Sea and Gibraltar Arc[END_REF], or subduction [START_REF] Lonergan | Origin of the Betic-Rif mountain belt[END_REF][START_REF] Gutscher | Evidence for active subduction beneath Gibraltar[END_REF]. [START_REF] Platt | Extensional collapse of thickened continental lithosphere: a working hypothesis for the Alboran Sea and Gibraltar Arc[END_REF] proposed that a lithospheric thickening as a result of convergence created a collisional ridge of thick crust underlain by a thick root of cold lithospheric mantle, during late Eocene to Oligocene time. Such a gravitationally unstable root could have been removed by convection and replaced by an asthenospheric mantle during the late Oligocene, resulting in the uplift and extension of the region and in the exhumation of metamorphic rocks and hightemperature peridotites from the base of the crust. The continuous convergence between Africa and Eurasia is then proposed to be accommodated by shortening in the external zones of the belts (Subbetic and External Rif).

The model proposed by Seber et al. (1996a) is also based on delamination, and is somewhat similar to the model of [START_REF] Platt | Extensional collapse of thickened continental lithosphere: a working hypothesis for the Alboran Sea and Gibraltar Arc[END_REF]. The difference is the style of removal of the thickened lithosphere. Seber et al. (1996a) proposed indeed that, as a consequence of compression, a thickened lithosphere was underlying the (thickened) crust, and that gravitationally unstable lithospheric root was progressively replaced by hot asthenosphere by convective removal from the crust toward W and NW. A consequence of the resulting asthenospheric inflow is the heating of the crust inducing an uplift and the extension in the Alboran area.

The slab rollback model of [START_REF] Royden | Evolution of retreating subduction boundaries formed during continental collision[END_REF] and [START_REF] Lonergan | Origin of the Betic-Rif mountain belt[END_REF] is very different from the above models, since it focuses primarily on subduction of oceanic lithosphere under the Alboran Sea. As the rate of subduction exceeds the rate of convergence, the subduction zone began migrating to the west while the crust of the former collision ridge (Internal Zone) was broken up and dispersed, and the crust behind this region was thinned by extension. In the northern and southern portions of the subduction zone, the westward movement slowed down as Iberian and African lithosphere was encountered and the emplacement of the Internal Zone was accommodated by shortening in the continental margins. The central portion of the subducting trench continued its westward migration and the continued subduction-zone rollback induced extension beneath the Alboran Sea. On the other hand, several other models consider different geometries and dip direction for the subducting slab: the above models invoke eastward subduction of oceanic lithosphere rolling back to the west [START_REF] Lonergan | Origin of the Betic-Rif mountain belt[END_REF][START_REF] Gutscher | Evidence for active subduction beneath Gibraltar[END_REF]; others consider an extinct subduction oriented northwards (e.g., [START_REF] Zeck | Betic-Rif orogeny:subduction of Mesozoic Tethys lithosphere under eastward drifting Iberia, slab detachment shortly before 22 Ma, and subsequent uplift and extensional tectonics[END_REF]. Recent analyses of body-wave dispersion at Ceuta [START_REF] Bokelmann | Mantle structure under Gibraltar constrained by dispersion of body waves[END_REF] indicate very different dispersion characteristics of events arriving from the West from those arriving from the East, which suggests the presence of an oceanic slab beneath the Alboran Sea.

Data and method

In this study, we used seismic data recorded by 13 permanent broadband stations of the Instituto Andaluz de Geofisica (IAG), and by 3 broadband seismic stations of the ROA-UCM-GEOFON (Real Observatorio de la Armada in San Fernando-Universidad Computence Madrid-Geofon). The station locations are reported in Fig. 1 and Table 1.

Locations and centroid times of teleseismic events used to extract the data from the continuous record were taken from the Harvard catalogue, http://www.seismologyharvard.edu/CMTsearch.html and arrival times of theoretical SKS phases at the stations were calculated using the theoretical Earth model IASP91 [START_REF] Kennett | Seismic traveltime table[END_REF] with the IASP91TTIM software [START_REF] Buland | The computation of seismic travel times[END_REF][START_REF] Kennett | Travel times for global earthquake location and phase identification[END_REF]. From the teleseismic events covering the period 2001-2006, we selected events occurring at distances larger than 85°a nd of magnitude (Mw) larger than 5.7. The selected events are shown in Fig. 2, with a global projection centered on southern Spain and preserving the azimuths. For the distance range between 130 and 165°w e also selected SKKS phases. A careful visual inspection of the data allowed us to keep about 207 seismograms from 67 teleseismic events.

For each selected event, we calculated the two splitting parameters, that is, the azimuth of the fast axis ϕ and the delay time δt between the fast and slow component of the shear waves. We performed these measurements using SplitLab [START_REF] Wüstefeld | Null detection in shear-wave splitting measurements[END_REF][START_REF] Wüstefeld | SplitLab: a shear-wave splitting environment in Matlab[END_REF], which utilizes three different techniques simultaneously: the rotation-correlation method (e.g. [START_REF] Bowman | Shear-wave splitting in the upper-mantle wedge above the Tonga subduction zone[END_REF], that is maximizing the cross-correlation between the radial and transverse component of the SKS phase, the minimum energy method [START_REF] Silver | Shear wave splitting and subcontinental mantle deformation[END_REF], that is minimizing the energy on the transverse component, and the minimum eigenvalue method [START_REF] Silver | Shear wave splitting and subcontinental mantle deformation[END_REF]. Fig. 3 shows an example of the application of these techniques for an event that arrives from the North at station ACBG. Under good measurement conditions, i.e. good signal-to-noise ratio and a favorable backazimuth, these methodologies give similar splitting parameters. The comparison between the different methods is helpful for characterizing cases of small δt where seismograms do not show clear SKS splitting, the so-called Null measurements, and for From left to right: 1) fast and slow components (dashed and continuous lines, respectively), corrected for the calculated splitting delay time; 2) corrected radial and transverse components (dashed and continuous lines, respectively), note that the energy on the transverse component is well removed after anisotropy correction; 3) the particle motion in the horizontal plane (dashed) becomes linear after the correction for anisotropy (solid); 4) contour plot for the maximum value of correlation coefficient and for the energy on transverse component as function of delay time and fast polarization angle. The shaded area marks the 95% confidence interval. (bottom) Numerical values that summarize the results for the three methodologies are given at the bottom, with their respective error ranges. discriminating them from cases of existent but weak splitting. Null measurements occur when the SKS phase is not split. This may happen when the medium is isotropic or when the backazimuth of the incoming SKS wave is parallel to either the slow or the fast directions in the anisotropic layer. In the case of a simple anisotropic medium, Null and non-Null measurements have to be consistent with each other. For simplicity in the discussion and pictures we present only measurements obtained with the minimum energy method. Table 1 gives resulting splitting parameters (ϕ and δt) for each seismogram.

We qualified the splitting results as "good", "fair", and "poor" as proposed by [START_REF] Barruol | Seismic anisotropy in the eastern US: deep structure of a complex continental plate[END_REF], i.e., depending on the quality of the seismograms, on the signal-to-noise ratio of the initial phase, on the amount of energy on the transverse component, on the correlation of the two split waveforms, and on the elliptical particle motion before anisotropy correction and its linearization after correction. Waveforms with good or fair quality generally provide similar results with and without filtering, indicating that results do not overly depend on the filter parameters. We generally choose to not filter the data unless required to remove high frequencies and/or long-period noise.

We will initially assume that the anisotropy is characterized by layer (s) with horizontal symmetry axes. Dipping axes of symmetry (or some orthorhombic symmetry systems) can cause systematic variations of splitting parameters as functions of backazimuth. However the small number of observations at each station precludes testing such cases.

Fig. 6. Backazimuthal distribution of SKS splitting for stations SESP, JAND, HORN, ARAC, and ANER that are not easily explained by a single anisotropic layer. The panels for ANER also show the prediction for the best-fitting two-layer model. See text for explanation.

Results of SKS splitting

The individual splitting measurements (see Table 2, given in the Electronic supplement) are plotted in Fig. 4 at the various stations to evaluate the coherence of the results. Lines show the fast anisotropy directions ϕ, obtained from the minimum energy method, and their lengths are proportional to the splitting delay time δt. "Good" data are shown by thick dark lines, and "fair" data by thin lines; "poor" data are not plotted. Most stations in the Eastern Betics, CART, VELZ, ACBG, GORA, ACLR and ASCB show a homogeneous ϕ trending roughly N70°E and high delay times of around 1.5 s. Given the weak spreading of fast direction measurements ϕ at stations in that region, and the spatial coherence, the splitting observed at those stations is likely originating within a single anisotropic layer that is well-constrained by a suitable number of measurements. We record a similar trend for the station MELI in Northern Africa, where the fast direction is N68°E with a delay time of roughly 0.95 s. Interestingly, ϕ is nearly North-South at stations SFS and CEUT, i.e., parallel to the trend of the Gibraltar arc. At each of these two stations unfortunately, measurement conditions were relatively difficult, and there are only two measurements available at CEUT from 3 yr of data analyzed, (N12°E, δt = 1.44 s and N19°W δt = 2.68 s), and two non-Null measurements at SFS (N19°E, δt = 1.80 s and N21°W δt = 0.4 s). [START_REF] Diaz | Anisotropy beneath the Iberian Peninsula: the contribution of the ILIHA-NARS broad-band experiment[END_REF] used data from temporary stations and found a fast direction oriented N80˚E at a station in the southern part of the Variscan Iberian Massif. In the Betics, they observed a NNE-SSW trend in the eastern part that changes to E-W in the central part, near Granada, and again NNE-SSW near the coast of Malaga interpreted in terms of their particular lithospheric geodynamic settings. [START_REF] Schmid | Delay times and shear wave splitting in the Mediterranean region[END_REF] studied the station CART, in the eastern part of Betic, and MELI in the north of Africa, finding fast directions in the Betic trend, oriented N70°E and N56°E respectively.

At a number of stations, there is a clear backazimuthal dependence of the splitting parameters, like at stations ARAC, HORN, JAND, SESP, and ANER. Such a backazimuthal dependence of splitting parameters may be explained either by invoking vertical variation of anisotropy, e.g. by the presence of two (or more) anisotropic layers, or by lateral variations of anisotropy. First, we address the case of lateral heterogeneity. In order to better visualize the lateral sampling of the anisotropy measurements and the regional-scale stability of the anisotropy parameters beneath the stations, we choose to project the individual splitting measurements of good and fair quality along their incoming rays to a depth of 150 km (Fig. 5). Although the splitting of the SKS waves is not well-located vertically, this representation is helpful since it plots the splitting parameters closer to the location where they were acquired. The measurements with well-constrained anisotropy ("good") are presented as thick dark lines while measurements of lower quality ("fair") are shown by thin lines. In addition, Fig. 5 presents the Null measurements by open circles. In the eastern and internal Betics, this figure clearly indicates smooth spatial variation of the anisotropic structures beneath most stations (with the exception of ANER, see discussion below). Such a pattern was less clear in Fig. 4. It suggests that the variation is not caused by noise in the observations. It is apparently due to real anisotropic variations at depth. This suggests that the upper mantle beneath the Internal Betics is dominated by a single anisotropic layer, with smoothly varying fast directions. Station ESTP located in the Subbetic domain shows a homogeneous distribution of the fast directions that led us to conclude that this part is characterized by a single layer. We will address the nature of this anisotropy and its smooth spatial variation in detail below.

The splitting parameters observed at stations ARAC, HORN and JAND near the transition to the Variscan Iberian Massif, SESP in the Prebetic units, and ANER in the Internal Betics, are not fully explained by lateral variation of splitting. These 5 stations are characterized by an apparent backazimuthal variation of the splitting parameters as shown in Fig. 4. The projection at depth of the individual splitting measurements (Fig. 5) suggests strong lateral variations at small scalelength. Such a rapid variation is not easily explained by lateral variation since the Fresnel zone of SKS phases in the upper mantle that quantifies the width of the zone over which the SKS 'average' the structure, is 100 to 200 km wide for SKS waves with 10 s period. For these 5 stations, we show the detailed ϕ and δt backazimuthal variation in Fig. 6. A single anisotropic layer would produce a constant fast direction ϕ, independent of backazimuth, as well as a constant splitting delay δt. With the possible exception of station ARAC, such a simple single layer does not explain the observations at these stations. However, the stations are characterized by varying amounts of data, partly due to the different running periods of the stations. Data at station JAND were available only from July 2006 and therefore provide only few SKS measurements. Nevertheless, the station exhibits a general NE trend except for the event on day 237 (year 2006, see Table 2 Electronic supplement) with backazimuth of 234°that shows a different ϕ. The small number of data cannot constrain the results for that station. For ANER, there is a nearly North-South trend for such events that arrive from ENE, but data with NNE backazimuth show NE-SW fast directions. The stations HORN, SESP and ARAC exhibit a good number of data which will allow the discussion of a possible two-layer case that we will present below.

Discussion

A striking feature in the anisotropy pattern observed in this study is that the Eastern Betics seems to be well-described by a single and rather homogeneous anisotropic layer, whereas a more complex anisotropic model likely involving lateral and vertical variations of anisotropy seems to be present in the upper mantle beneath stations in the North, and perhaps the West and the South. The Eastern Betics show indeed a smooth lateral variation of anisotropy, whereas stations in the Western and Southern part of this region indicate much stronger lateral variations in the splitting parameters. There are also strong lateral variations associated with the transition from the Betic to Variscan Iberia, in addition to two-layer anisotropy as we will show below, as is required at some of the stations at least.

Eastern Betics

Fast directions are oriented ENE-WSW in the Eastern Betics, and there is a smooth spatial variation clearly visible in Fig. 5 that is more or less parallel to the mountain belt. The anisotropy in that region is easily explained by a single anisotropic layer. The average splitting delay time in that region is about 1.5 s, which is consistent with a layer that is 150 km thick assuming the typically used value of 4% anisotropy in the upper mantle [START_REF] Mainprice | Interpretation of SKS-waves using samples from the subcontinental lithosphere[END_REF][START_REF] Mainprice | The seismic anisotropy of the Earth's mantle: from single crystal to polycrystal[END_REF]. Tomographic models for the region, (e.g. [START_REF] Blanco | The P-wave velocity structure of the mantle below the Iberian Peninsula: evidence for subducted lithosphere below southern Spain[END_REF], show a pronounced low-velocity anomaly in the upper mantle under the Betics between 30 and 100 km depth that is also apparent and even more pronounced in the Pn tomographic model of Calvert et al. (2000a). This area is also the locus of a pronounced high attenuation anomaly of Sn waves (Calvert et al., 2000b). This may suggest that continental mantle material under the Betics has been replaced by asthenospheric material. Geochemical studies favor such a hypothesis by suggesting an edge delamination that could have occurred under the Betics [START_REF] Duggen | Post-collisional transition from subduction-to intraplate-type magmatism in the westernmost Mediterranean: evidence for continental-edge delamination of subcontinental lithosphere[END_REF]. If the lithosphere has indeed been removed, this should suggest that the deformation recorded by the seismic anisotropy under the Betics is associated rather with hot and thus relatively low-viscosity material. At larger depths, however, velocities are again relatively high in tomographic models (e.g., [START_REF] Blanco | The P-wave velocity structure of the mantle below the Iberian Peninsula: evidence for subducted lithosphere below southern Spain[END_REF] down to the transition zone, including the zone of very deep seismicity around 600 km depth [START_REF] Buforn | Source mechanism of intermediate and deep earthquakes in southern Spain[END_REF].

The ENE-WSW fast directions observed in the Eastern Betics are nearly parallel to the absolute plate motion direction, which is 239° and 2 cm/yr for the motion model Nuvel1A relative to the hotspot reference system HS3 [START_REF] Gripp | Young tracks of hotspots and current plate velocities[END_REF]. Although the Eurasian plate motion vector is slow and therefore still matter of debate, a simple deformation model related to plate motion of the Iberian/Eurasian plate over the deeper mantle might thus explain the anisotropy. In this context, the gradual rotation of fast direction along the coast might be explained reasonably by invoking a deviation of mantle flow around the Iberian lithosphere, following similar suggestions of [START_REF] Bormann | Upper mantle anisotropy beneath Central Europe from SKS wave splitting: effects of absolute plate motion and lithosphere asthenosphere boundary topography?[END_REF] for central Europe and of [START_REF] Barruol | Seismic anisotropy in the eastern US: deep structure of a complex continental plate[END_REF] and [START_REF] Fouch | Shear wave splitting, continental keels, and patterns of mantle flow[END_REF] for North America. This might also explain the minor difference of 10°between the average fast direction in the Eastern Betics and the absolute plate motion direction.

On the other hand, the observed fast directions in the Eastern Betics also show an interesting correlation with crustal features which is not expected in the case of plate motion induced asthenospheric anisotropy. The convergence direction between the African and the Iberian/Eurasian plate is NNW-SSE, and the current style of deformation in the crust is predominately left-lateral strike-slip faults with a NE-SW trend (e.g. Alhama de Murcia, Palomares and Carboneras faults). [START_REF] Vauchez | Mountain building: strike parallel motion and mantle anisotropy[END_REF] suggested that in many collisional belts the dominant motion is strike-slip parallel to the main trend of the belt. In this case, both crustal and mantle fabrics may reflect these movements, and the entire lithosphere develops a fabric and becomes anisotropic. They attribute this anisotropy to a dominant mantle flow parallel to the mountain belt during orogeny. Such a structure is well illustrated in the Ronda peridotite massif where pervasive structures and fabric in the lithospheric and asthenospheric mantle are well preserved [START_REF] Vauchez | Seismic properties of an asthenospherized lithospheric mantle: constraints from lattice preferred orientations in peridotite from the Ronda Massif[END_REF]. The present-day extension direction inferred by [START_REF] Stich | Kinematics of the Iberia-Maghreb plate contact from seismic moment tensors and GPS observations[END_REF], using moment tensors of crustal earthquakes, is N240°E, which is indeed parallel to the inferred fast polarization directions. The anisotropy is thus also consistent with a notion of vertically coherent deformation of crust and mantle as initially proposed by Silver andChan (1988, 1991), except that we are dealing with hot deformed mantle in this region, and thus probably with current deformation, rather than fossil deformation remaining from the creation of continental lithosphere. For the close-by Gulf of Cadiz region, [START_REF] Stich | Crust-mantle coupling in the Gulf of Cadiz (SW-Iberia)[END_REF] have shown that the stresses imposed by the Africa-Iberia plate convergence have the same geometry in the crust and the mantle. The convergence thus imposes a similar lateral boundary condition over a considerable depth range. Crust and mantle in the Betics may thus both be constrained to deform in a similar fashion. The observed anisotropy under the Eastern Betics could be thus caused by a combination of this kind of vertically coherent deformation, and absolute-motion-related deformation, both producing a seismic anisotropy such as the observed one. The relative importance of each is not constrained though. Interestingly, a model of slab rollback is also consistent with such coherence between crustal and upper mantle deformation. The retreat to the WSW of the Alboran slab may have induced a large transcurrent deformation in the mantle beneath the eastern Betics but also in its crust.

Variscan belt

The northernmost stations of our study area are situated on or near the Variscan Iberian Massif. The fast anisotropy directions we observe at these stations are trending more NE-SW than those recorded by stations in the Betics that trend more ENE-WSW. These directions do not appear to correspond to orientations of geological structures in the Variscan basement, since the latter trend more or less NW-SE in the Northwest portion of Iberia and East-West in eastern Spain. While we do not obtain those "Variscan" directions for our northernmost stations that are near the southern end of the Variscan belt, they are in fact observed further North, in the center of Variscan Iberia [START_REF] Silver | Implications for continental structure and evolution from seismic anisotropy[END_REF][START_REF] Schmid | Delay times and shear wave splitting in the Mediterranean region[END_REF], where they are roughly E-W. Three of our northernmost stations show a wide range of fast orientations (Figs. 5 and6). Different anisotropy parameters are obtained for measurements arriving from the North and South. This is not easily explained by lateral heterogeneity, unless all these stations are located above a geological suture zone so that events arriving from the North and the South could experience a different medium. Since these stations are indeed near the Betics deformation front, we wish to further test both hypotheses for the northernmost stations, namely the 'lateral heterogeneity hypothesis' and the 'two-layer anisotropy hypothesis'. Fig. 5 illustrates the 'lateral variation hypothesis', where splitting results were projected along the incoming rays to 150 km depth. In this view, fast directions that show up to the North of the stations, within the Iberian Variscan Massif, are generally associated with NE-SW fast direction, while results generally show more East-West trending fast directions, for southern backazimuth i.e., close to the general Betics trend. This is particularly clear at station SESP, and it may indicate that the events arriving from the south are experiencing the same anisotropic structure as the stations located further south in the Internal Betics. This would suggest that SESP may lie close to an important boundary between lithospheric blocks. Among the northern stations of our study area, ARAC shows little backazimuthal dependence of the anisotropy parameters: NE-SW fast directions are visible both to the North and to the South of the station.

On the other hand, we may assume a two-layer model for each of these four northern stations, to try to explain the observed variation of the splitting parameters apparent in Fig. 6. Constraining the four parameters of a two-layer model requires a relatively good backazimuthal coverage (see [START_REF] Walker | Seismic anisotropy in the asthenosphere beneath the Eifel region[END_REF], which we do not have. We therefore decided to apply a similar approach as [START_REF] Fontaine | Upper-mantle flow beneath French Polynesia from shear wave splitting[END_REF] for the SKS splitting obtained at oceanic island stations, by considering the area as homogeneous enough to group the individual splitting measurements together into a single "virtual" station. We processed measurements obtained from the four stations together, since they have similar backazimuthal coverages. Different from the other three stations, SESP is located within the Betic Cordillera, in the so-called Pre-Betic. Nevertheless splitting parameters are more like those of the stations at the limit, e.g. HORN, rather than the other station on the Betics.

If there were several anisotropic layers present at depth, mantle deformation would not necessarily mimic crustal deformation. A tectonic decoupling somewhere within the lithosphere might explain why anisotropic directions in the area are different from the trend of the mountain belt. Seismic reflection profiles in the south-western portion of the Iberian Massif show that the Moho is discontinuous [START_REF] Simancas | Crustal structure of the transpressional Variscan orogen of SW Iberia: SW Iberia deep seismic reflection profile (IBERSEIS)[END_REF], particularly under the Ossa-Morena and Central Iberian Massif. In fact, the Moho is probably a decoupling zone between the mantle and the crust in the area.

Following the approach defined by [START_REF] Silver | The interpretation of shear-wave splitting parameters in the presence of two anisotropic layers[END_REF], we therefore tested two-layer models by varying ϕ and δt in the upper and lower layer by increments of 2°and 0.2 s, respectively, and compared the expected backazimuthal variations with our observations. We use the statistical technique presented in [START_REF] Walker | Shear-wave splitting around the Eifel hotspot: evidence for a mantle upwelling[END_REF] to judge the significance of the variance reduction over the best-fitting one-layer model. Testing all possible two-layer models (Fig. 7), we find that the best-fitting models belong to three groups of models, 1) those having NE-SW fast directions in both layers, 2) those having NE-SW for the upper layer and NW-SE for the lower layer, and 3) vice-versa. While the best model formally falls into the second category, we cannot rule out the other two cases, including the one of a shallower layer with NW-SE fast azimuths, that would indeed correspond to the trend of the Variscan basement in its southern portion. This ambiguity is in part due to the backazimuthal distribution of the events being rather limited, but more fundamentally since a second layer with exactly perpendicular orientation to the first one only varies the splitting delay but does not produce a backazimuthal variation of fast direction and splitting delay. We note that the splitting delay associated with the lower layer is quite small (0.2 s). While the additional variance reduction (39%) associated with this weak second layer is slightly statistically significant, we choose to not attribute much significance to this. The backazimuthal variation that is present in the data may likely represent a complexity that is perhaps better characterized by an anisotropy that varies laterally.

Gibraltar arc: constraints on geodynamic models

Fast directions in Southernmost Spain, close to Gibraltar show rather distinct orientations from those in the Internal or External Betics: ϕ trends more or less North-South at the two stations CEUT and SFS. Station ANER near Malaga provides evidence for higher complexity and perhaps multi-layer anisotropy. We have therefore performed a two-layer modeling as detailed above but consider only the results obtained at this station. The thousand best two-layer models obtained for ANER give a best-fitting model with fast direction ϕ for the upper layer of N28°E and for the lower layer of N68°E (Fig. 8), and splitting delays 1.0 and 0.6 s. This result suggests that this station lies in a region of transition between the eastern Betics where ϕ trend ENE (i.e., close to the ANER lower layer ϕ) and the Gibraltar arc where ϕ trend close to NS (i.e., close to the ANER upper layer ϕ). Interestingly, the close NS trend of the upper layer seems also to coincide with the N-S line of intermediate-depth earthquakes [START_REF] Morales | Active continental subduction beneath the Betic Cordillera and the Alborán Sea[END_REF]Calvert et al., 2000b). Nevertheless, we note that for the fast directions in the two layers, the maxima of the histograms in Fig. 8 are only 20°a part. Even though the statistical improvement is slightly significant, we defer a deeper discussion of a possible two-layer anisotropy until data from more stations around ANER become available. We consider both Pn and SKS anisotropy measurements, which provide different and complementary depth sampling (e.g., [START_REF] Pera | Petrophysical properties of the upper mantle beneath the Torre Alfina area (Northern Apennines, Central Italy)[END_REF]. SKS waves travel almost vertically and are not able to localize the precise depth range of the anisotropy, except that a layer of more than 100 km thickness is required to explain the observed splitting delay times [START_REF] Mainprice | Interpretation of SKS-waves using samples from the subcontinental lithosphere[END_REF][START_REF] Ben Ismaïl | An olivine fabric database: an overview of upper mantle fabrics and seismic anisotropy[END_REF], and that it should be in the upper mantle. On the other hand, Pn travels horizontally just below the Moho, and thus constrains anisotropy at that depth level. Pn models for the area have been produced by Calvert et al. (2000b) and by [START_REF] Serrano | Seismic anisotropy and velocity structure beneath the southern half of the Iberian peninsula[END_REF] which give rather similar results in the area of interest. In each case, arrival times of regional events have been used to determine isotropic and anisotropic velocity variations. The agreement of anisotropic fast directions from Pn and SKS are in fact remarkable on the Spanish side of the Alboran arc (Fig. 9). Not only are these fast directions parallel in the eastern Betics, but they seem to rotate together with the largescale geological structures of the Alboran arc, and parallel with the coast, in a similar fashion for the shallowest and somewhat deeper mantle. The model may perhaps suggest that this parallelism is not present on the African side of the Gibraltar arc where fast Pn trend NS to NNW-SSW and fast split SKS trend NE-SW. The model is less wellconstrained in that area, but if the feature is real, this may suggest a complete decoupling between the sub-Moho and the deeper anisotropy. However, the similarity between Pn and SKS fast directions may reflect different tectonic processes at depth, particularly in the vicinity of a subduction zone. Indeed, if there is subduction in the area, as has often been argued (e.g., [START_REF] Gutscher | Evidence for active subduction beneath Gibraltar[END_REF], Pn may be affected by the anisotropy above the slab, whereas SKS might be sensitive to the anisotropy below, above, and possibly within the slab. Fast split directions parallel to the slab is a rather common feature in subduction zones, both above and below the slab. This has been found for instance in New Zealand (Marson-Pidgeon and [START_REF] Marson-Pidgeon | Shear-wave splitting variations across an array in the southern North Island, New Zealand[END_REF] and in Northern Italy [START_REF] Mele | Compressional velocity structure and anisotropy in the uppermost mantle beneath Italy and surrounding regions[END_REF]. Trench-parallel fast directions, may be explained either by invoking trench-parallel flow (e.g., [START_REF] Russo | Trench-parallel flow beneath the Nazca Plate from seismic anisotropy[END_REF], or by invoking a hydrated mantle wedge above the slab [START_REF] Jung | Water-induced fabric transitions in olivine[END_REF]. The degree of hydration of the mantle wedge and of magnitude of the applied stress above a subduction zone appears to have a strong influence on the style of seismic anisotropy (e.g., [START_REF] Jung | Water-induced fabric transitions in olivine[END_REF][START_REF] Kaminski | Timescales for the evolution of seismic anisotropy in mantle flow[END_REF]: Slab-parallel seismic fast axes from SKS can in principle be explained either by slab-parallel flow in an anhydrous mantle wedge or by slab-normal compression in a hydrous mantle wedge. Such a pattern of anisotropy may show apparently erratic behaviour in some places; this may perhaps help to explain the complexity that we find at station ANER.

Deep structure and dynamics beneath the Alboran area has been widely debated throughout the last decades and a wide variety of geodynamic models has been proposed for the area. These models can be categorized by the way material is recycled into the mantle. If mantle lithosphere has been removed as has been proposed in the convection-removal hypothesis [START_REF] Platt | Extensional collapse of thickened continental lithosphere: a working hypothesis for the Alboran Sea and Gibraltar Arc[END_REF]) and the delamination model (e.g., Seber et al., 1996a), this would require the removed material to be replaced by asthenospheric material. This replacement would correspond to a flow that is directed radially inwards. An idealized convective-removal model would in fact imply that such a flow would be acting radially inwards from all directions. One should not expect to observe a radial flow from all directions though, since any asymmetry or heterogeneity would easily lead to an asymmetric flow, more like that to be expected in a delamination model. However, a radial flow would necessarily exist somewhere, to replace the material that has been recycled into the mantle. This flow would produce anisotropic fast directions that are oriented radially. Fig. 5 shows however, that the flow is oriented more or less tangentially to the Alboran Sea region. The required radial flow is thus not observed. This suggests that convective-removal and delamination models are unlikely to explain the anisotropy observations. On the other hand, subduction to the Northwest, as proposed by [START_REF] Zeck | Betic-Rif orogeny:subduction of Mesozoic Tethys lithosphere under eastward drifting Iberia, slab detachment shortly before 22 Ma, and subsequent uplift and extensional tectonics[END_REF], would require a rollback to the South to produce trenchparallel flow. Focal mechanisms of intermediate-depth seismicity are consistent with this, but hypocenters do not show a dip toward the North. In fact, they do not show a clear Benioff zone at all, but a line of seismicity dipping steeply to the South. For that same area, teleseismic and regional tomographic studies indicate a pronounced high-velocity anomaly located beneath the Alboran Sea (e.g. [START_REF] Blanco | The P-wave velocity structure of the mantle below the Iberian Peninsula: evidence for subducted lithosphere below southern Spain[END_REF]Calvert et al., 2000b;Seber et al., 1996a). Calvert et al. (2000b) produced a tomographic model that showed a robust high-velocity anomaly beneath the west Alboran Sea and west and central Betics between 100 and 150 km. Assuming that these positive velocity anomalies indicate a region of cooler mantle, the most likely interpretation is that it indicates the presence of lithospheric mantle in view of the coherent continuation of the anomaly at depth. The positive anomaly extends from lithospheric depths beneath the Strait of Gibraltar and southern Spain to depths of about 350-400 km beneath the Alboran Sea, and it is interpreted as a lithospheric body that has descended into the upper mantle. More recently, waveform studies [START_REF] Bokelmann | Mantle structure under Gibraltar constrained by dispersion of body waves[END_REF] have indicated dispersed body-wave arrivals propagating through the mantle under the Alboran Sea, that favor the presence of a subducted slab under the Alboran Sea, still containing a continuous low-velocity crustal layer. [START_REF] Wortel | Subduction and slab detachment in the Mediterranean-Carpathian region[END_REF] also proposed a slab structure isolated in the mantle below the Betic-Alboran region at 200 km depth and that seems to be connected with deep seismicity at depths of about 640 km. [START_REF] Lonergan | Origin of the Betic-Rif mountain belt[END_REF] proposed that there was a subduction zone beneath the Alboran Sea that is now extinct. They suggested that the ancient subducting slab in the Western Mediterranean was split into two fragments, the eastern one of which has continued to roll back toward the south-east, generating the Tyrrhenian Sea and forming the present-day Calabrian arc. The other slab fragment rolled back to the west, generating the Alboran Sea and the Betic-Rif orocline. Our SKS splitting observations, together with the previously published Pn results, show a spectacular rotation following the curvature of the roll back and of the mountain chain. If the SKS fast directions are to be explained by a rotation similar to that shown by the Pn observations, this would be explained most easily by a toroidal mantle flow under the subducting lithosphere that shows rollback toward the west. There is an interesting similarity with the Calabrian arc in Southern Italy, for which shear-wave splitting results have been presented by [START_REF] Civello | Toroidal mantle flow around the Calabrian slab (Italy) from SKS splitting[END_REF] and more recently by [START_REF] Baccheschi | Seismic anisotropy reveals focused mantle flow around the Calabrian slab (Southern Italy)[END_REF]. The asthenospheric toroidal flow induced by the Tyrrhenian slab roll back is presently beneath Calabria [START_REF] Civello | Toroidal mantle flow around the Calabrian slab (Italy) from SKS splitting[END_REF][START_REF] Baccheschi | Seismic anisotropy reveals focused mantle flow around the Calabrian slab (Southern Italy)[END_REF] but may have left some imprint in the mantle all along its Neogene travel path from southern France [START_REF] Barruol | Mapping Upper mantle anisotropy beneath SE France by SKS splitting indicates a Neogene asthenospheric flow induced by the Apenninic slab rollback and deflected by the deep Alpine roots[END_REF] to Sicily [START_REF] Lucente | Seismic anisotropy reveals the long route of the slab through the western-central Mediterranean mantle[END_REF]. Seismic anisotropy beneath the Betics is likely related to a similar phenomenon, and thus to the large-scale dynamics of ocean opening in the Southwestern Mediterranean. Calvert et al., 2000b). See text for explanations.

Conclusions

We have studied upper mantle seismic anisotropy in Southern Spain and around the Gibraltar arc. Our observations suggest a consistent splitting direction parallel to the mountain belt in the Internal Betics, a more northerly direction at stations located further north in the transition with the Variscan Iberia, and a rotation of fast split directions toward NS azimuths around the Gibraltar arc.

The Eastern Betics anisotropy observations are well explained by a single layer of anisotropy with a fast anisotropic direction trending ENE, with smooth and small variations in space. In the transition region between the Variscan and the Betics, the anisotropy shows backazimuthal variations of splitting parameters. Models of twolayer anisotropy and lateral heterogeneous anisotropy have been tested, and a strong lateral change in anisotropy is likely for the transition region toward the Variscan Iberian Massif. Close to Gibraltar, fast directions are oriented close to N-S at CEUT and SFS and this suggests a gradual rotation around the arc. Combining Pn and SKS anisotropy, geological features and tomographic images allow us to discuss the various families of geodynamic models. Our SKS measurements are difficult to reconcile with geodynamic models based on delamination or convective removal. The observations instead favor models of subduction and rollback of the western Mediterranean slab.

Fig. 1 .

 1 Fig. 1. Tectonic sketch of the Mediterranean area showing the Neogene basins and Alpine chains (top, redrawn from Comas et al., 1999) and geological setting of South of Iberia and North of Africa (bottom, redrawn from Lonergan and White, 1997; Platt el al., 2003). Main faults are shown in the External and Internal Zones and in the Alboran Domain (CVF: Crevillente fault, AMF: Alhama de Murcia fault, PF: Palomares fault, CFZ: Carboneras fault, NF: Nekor fault). Locations of seismological stations used in this study are also shown by triangles.

Fig. 2 .

 2 Fig.2. Spatial distribution of earthquakes used in this study in a projection that preserves backazimuths, centered on the Alboran. The 85°(minimum distance for SKS splitting measurement) and 130°of epicentral distance are also shown.

Fig. 3 .

 3 Fig.3. Example of data processing using SplitLab software(Wüstefeld et al., 2008) at ACBG station. (top) The upper panel shows the event information and the initial radial and transverse components. Dashed lines are the predicted phase arrival times from IASP91; the shaded area is the selected window for the shear-wave splitting processing. (middle) Figures for inspecting the performance are shown with four plots for the rotation-correlation method in panel a) and four plots for the minimum energy method in panel b). From left to right: 1) fast and slow components (dashed and continuous lines, respectively), corrected for the calculated splitting delay time; 2) corrected radial and transverse components (dashed and continuous lines, respectively), note that the energy on the transverse component is well removed after anisotropy correction; 3) the particle motion in the horizontal plane (dashed) becomes linear after the correction for anisotropy (solid); 4) contour plot for the maximum value of correlation coefficient and for the energy on transverse component as function of delay time and fast polarization angle. The shaded area marks the 95% confidence interval. (bottom) Numerical values that summarize the results for the three methodologies are given at the bottom, with their respective error ranges.

Fig. 4 .

 4 Fig. 4. Splitting results shown at the seismic stations used in the present study. Each measurement is characterized by a fast azimuth ϕ, and the length of the segment is proportional to the delay time δt. Good (thick dark lines) and fair (thin lines) are represented.

Fig. 5 .

 5 Fig. 5. Individual splitting measurements projected at 150 km depth along the ray paths. Good measurements are shown by thick lines, fair measurements by thin lines. Null measurements (events with no apparent splitting) are shown by circles. White arrows are average stress extension directions calculated by Stich et al. (2006), the dark arrow shows absolute plate motion direction for the model HS3-Nuvel-lA (Gripp and Gordon, 2002).

Fig. 7 .

 7 Fig. 7. Statistical analysis of two-layer modelling for the joint data set of 4 stations in the north (ARAC + HORN + JAND + SESP): parameters for the 1000 best-fitting models are shown as histograms for a) fast direction in upper layer; b) delay time in the upper layer c) fast direction in the lower layer; d) delay time in the lower layer. The arrows show the optimum values for each model, ϕ = N52°E and ϕ = N136°E for a) and c) respectively, and δt = 1.2 s and δt = 0.2 s for b) and d) respectively.

Fig. 8 .

 8 Fig. 8. Statistical analysis of two-layer modelling for station ANER: parameters for the 1000 best-fitting models are shown as histograms for a) fast direction in upper layer; b) delay time in the upper layer; c) fast direction in the lower layer; d) delay time in the lower layer. The arrows show the optimum values for each model, ϕ = N28°E and ϕ = N68°E for a) and c) respectively, and δt = 1.0 s and δt = 0.6 s for b) and d) respectively.

Fig. 9 .

 9 Fig. 9. Comparison of SKS fast directions (weighted-mean values from this study) with Pn-fast directions (redrawn fromCalvert et al., 2000b). See text for explanations.

Acknowledgements

We warmly thank Andreas Wüstefeld for his useful help in the use of SplitLab software, and Alain Vauchez, Daniel Stich, Carlos Garrido and Francisco Vidal for fruitful discussions. We also thank Martha Savage and an anonymous reviewer for their comments and corrections. SplitLab is available at http://www.gm.univ-montp2.fr/ splitting/ together with its user guide, the related publications and the SKS World splitting database. Figures are created using the GMT software of Wessel and Smith (1995). This study was possible thanks to the effort of collecting high-quality broadband seismic data of the IAG-UGR and ROA-UCM-Geofon networks. We received support through Spanish project CGL2005-04541-C03-01-BTE, European Commission project NEAREST (GoCE-037110) and Junta de Andalucia RNM104 and through the project Topo-Iberia CSD2006.

Appendix A. Supplementary data

Supplementary data associated with this article can be found, in the online version, at doi:10.1016/j.epsl.2008.06.024.