# Highlights on mantle deformation beneath the Western Alps with seismic anisotropy using CIFALPS2 data

Silvia Pondrelli, Simone Salimbeni, Judith M. Confal, Marco G Malusà, Anne Paul, Stephane Guillot, Stefano Solarino, Elena Eva, Coralie Aubert, Liang Zhao

► To cite this version:

Silvia Pondrelli, Simone Salimbeni, Judith M. Confal, Marco G Malusà, Anne Paul, et al. Highlights on mantle deformation beneath the Western Alps with seismic anisotropy using CIFALPS2 data. Solid Earth, In press. hal-04618966

HAL Id: hal-04618966

https://hal.science/hal-04618966

Submitted on 20 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Highlights on mantle deformation beneath the Western Alps with seismic anisotropy using CIFALPS2 data

Silvia Pondrelli, Simone Salimbeni, Judith M. Confal, Marco G Malusà, Anne Paul, Stephane Guillot, Stefano Solarino, Elena Eva, Coralie Aubert, Liang Zhao

► To cite this version:

Silvia Pondrelli, Simone Salimbeni, Judith M. Confal, Marco G Malusà, Anne Paul, et al. Highlights on mantle deformation beneath the Western Alps with seismic anisotropy using CIFALPS2 data. Solid Earth, 2024, 15 (7), 14p. 10.5194/egusphere-2024-468. http://hal-04618966

HAL Id: hal-04618966

https://hal.science/hal-04618966v1

Submitted on 20 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Highlights on mantle deformation beneath the Western Alps with seismic anisotropy using CIFALPS2 data

Silvia Pondrelli<fnref n="1" />, Simone Salimbeni<fnref n="1" />, Judith M. Confal<fnref n="1" />, Marco G. Malusà<fnref n="2" />, Anne Paul<fnref n="3" />, Stephane Guillot<fnref n="3" />, Stefano Solarino<fnref n="4" />, Elena Eva<fnref n="4" />, Coralie Aubert<fnref n="3" />, Liang Zhao<fnref n="5" />

*Correspondence to:* Silvia Pondrelli (silvia.pondrelli@ingv.it)

**Abstract.** There are still open questions about the deep structure beneath the Western Alps. Seismic velocity tomographies show the European slab subducting beneath the Adria plate, but all these images did not clarify completely the possible presence of tears, slab windows, or detachments. Seismic anisotropy, considered as an indicator of mantle deformation and studied using data recorded by dense networks, allows a better understanding of mantle flows in terms of location and orientation at depth. Using the large amount of shear wave splitting and splitting intensity measurements available in the Western Alps, collected through the CIFALPS2 temporary seismic network, together with already available data, some new patterns can be highlighted and gaps left by previous studies can be filled. Instead of the typical seismic anisotropy pattern parallel to the entire arc of the Western Alps, this study supports the presence of a differential contribution along the belt, only partly related to the European slab steepening. A nearly NS anisotropy pattern beneath the external Western Alps, a direction that cuts the morphological features of the belt, is clearly found with the new CIFALPS2 measurements. It is however confirmed that the asthenospheric flow from Central France towards the Tyrrhenian Sea, is turning around the southern tip of the European slab.

## 1 Introduction

Seismic anisotropy has become a convincing study tool, mainly in the areas where recent and past geodynamic evolution have left their marks in the mantle deformation and its patterns (e.g. Long and Silver, 2009; Long, 2013 for a complete review). Several methods exist to measure seismic anisotropy. Depending on the used seismic phase, signal, and frequency, it is possible to measure seismic anisotropy at different depths and relate it to different parts of the Earth's structure, last but not least with seismic anisotropy tomography (e.g. Zhao et al., 2023). Local seismicity and surface waves are used to measure crustal seismic anisotropy, usually attributed to fractures or/and the state of stress of crustal depths (i.e. Crampin and

<footnote n="1">¹ Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna, Bologna, Italy</footnote>

<footnote n="2">² Department of Earth and Environmental Sciences, University of Milano-Bicocca, Piazza della Scienza 4, 20126 Milan, Italy</footnote>

<footnote n="3">³ Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, UGE, ISTerre, Grenoble, France</footnote>

<footnote n="4">⁴ Istituto Nazionale di Geofisica e Vulcanologia, ONT, Genova, Italy</footnote>

<footnote n="5">⁵ State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China</footnote>
# Highlights on mantle deformation beneath the Western Alps with seismic anisotropy using CIFALPS2 data

Silvia Pondrelli<fnref n="1" />, Simone Salimbeni<fnref n="1" />, Judith M. Confal<fnref n="1" />, Marco G. Malusà<fnref n="2" />, Anne Paul<fnref n="3" />, Stephane Guillot<fnref n="3" />, Stefano Solarino<fnref n="4" />, Elena Eva<fnref n="4" />, Coralie Aubert<fnref n="3" />, Liang Zhao<fnref n="5" />

*Correspondence to:* Silvia Pondrelli (silvia.pondrelli@ingv.it)

**Abstract.** There are still open questions about the deep structure beneath the Western Alps. Seismic velocity tomographies show the European slab subducting beneath the Adria plate, but all these images did not clarify completely the possible presence of tears, slab windows, or detachments. Seismic anisotropy, considered as an indicator of mantle deformation and studied using data recorded by dense networks, allows a better understanding of mantle flows in terms of location and orientation at depth. Using the large amount of shear wave splitting and splitting intensity measurements available in the Western Alps, collected through the CIFALPS2 temporary seismic network, together with already available data, some new patterns can be highlighted and gaps left by previous studies can be filled. Instead of the typical seismic anisotropy pattern parallel to the entire arc of the Western Alps, this study supports the presence of a differential contribution along the belt, only partly related to the European slab steepening. A nearly NS anisotropy pattern beneath the external Western Alps, a direction that cuts the morphological features of the belt, is clearly found with the new CIFALPS2 measurements. It is however confirmed that the asthenospheric flow from Central France towards the Tyrrhenian Sea, is turning around the southern tip of the European slab.

## 1 Introduction

Seismic anisotropy has become a convincing study tool, mainly in the areas where recent and past geodynamic evolution have left their marks in the mantle deformation and its patterns (e.g. Long and Silver, 2009; Long, 2013 for a complete review). Several methods exist to measure seismic anisotropy. Depending on the used seismic phase, signal, and frequency, it is possible to measure seismic anisotropy at different depths and relate it to different parts of the Earth's structure, last but not least with seismic anisotropy tomography (e.g. Zhao et al., 2023). Local seismicity and surface waves are used to measure crustal seismic anisotropy, usually attributed to fractures or/and the state of stress of crustal depths (i.e. Crampin and

<footnote n="1">¹ Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna, Bologna, Italy</footnote>

<footnote n="2">² Department of Earth and Environmental Sciences, University of Milano-Bicocca, Piazza della Scienza 4, 20126 Milan, Italy</footnote>

<footnote n="3">³ Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, UGE, ISTerre, Grenoble, France</footnote>

<footnote n="4">⁴ Istituto Nazionale di Geofisica e Vulcanologia, ONT, Genova, Italy</footnote>

<footnote n="5">⁵ State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China</footnote>
Peacock, 2008; Okaya et al., 2018). Using seismic signals that travel deeper, it is possible to sample the deformation at lithospheric-asthenospheric depths. For instance, Pn phases record the seismic anisotropy immediately below the Moho in the lithospheric mantle (i.e. Diaz et al., 2013). On the other hand, using core phases (SKS, SKKS, etc) that record information on the receiver branch of their path, we can sample seismic anisotropy that is thought to be mainly concentrated in the upper mantle.

<figure><img src="image_2.png" /><figcaption>Figure 1 - a) Map of the study region, focusing on the Western Alps. In red are indicated the CIFALPS2 stations, while in blue are permanent and previous temporary stations (i.e. CIFALPS and AlpArray). FPF = Frontal Pennine Fault, BG = Bresse Graben, PP = Po Plain, LM = Ligurian Mountains; b) the red square is the study area reported in a); c) map of all seismic events used in this study, with the star centered in the study region.</figcaption></figure>

The azimuth of the fast velocity direction and the delay time, the two parameters that commonly result from shear wave splitting analysis of core phases, are interpreted respectively as the direction assumed by olivine crystals, the principal mineral component of the upper mantle when mantle undergoes deformation, and the amount of anisotropy crossed by a seismic ray.

In the Alps these kinds of measurements improved immensely with recent temporary experiments such as AlpArray (Hetényi et al., 2018) or CIFALPS and CIFALPS2 in the western sector of the chain (Zhao et al., 2015; 2016; 2018), which complemented the permanent seismic networks operating in the region.
Peacock, 2008; Okaya et al., 2018). Using seismic signals that travel deeper, it is possible to sample the deformation at lithospheric-asthenospheric depths. For instance, Pn phases record the seismic anisotropy immediately below the Moho in the lithospheric mantle (i.e. Diaz et al., 2013). On the other hand, using core phases (SKS, SKKS, etc) that record information on the receiver branch of their path, we can sample seismic anisotropy that is thought to be mainly concentrated in the upper mantle.

<figure><img src="image_2.png" /><figcaption>Figure 1 - a) Map of the study region, focusing on the Western Alps. In red are indicated the CIFALPS2 stations, while in blue are permanent and previous temporary stations (i.e. CIFALPS and AlpArray). FPF = Frontal Pennine Fault, BG = Bresse Graben, PP = Po Plain, LM = Ligurian Mountains; b) the red square is the study area reported in a); c) map of all seismic events used in this study, with the star centered in the study region.</figcaption></figure>

The azimuth of the fast velocity direction and the delay time, the two parameters that commonly result from shear wave splitting analysis of core phases, are interpreted respectively as the direction assumed by olivine crystals, the principal mineral component of the upper mantle when mantle undergoes deformation, and the amount of anisotropy crossed by a seismic ray.

In the Alps these kinds of measurements improved immensely with recent temporary experiments such as AlpArray (Hetényi et al., 2018) or CIFALPS and CIFALPS2 in the western sector of the chain (Zhao et al., 2015; 2016; 2018), which complemented the permanent seismic networks operating in the region.
The European Alps originated in the late Cretaceous from the oblique subduction of the Alpine Tethys under the Adria microplate. The subduction evolved in a continental collision during the late Cenozoic (e.g., Handy et al., 2010, 2013 and references therein). In the Western Alps, the tectonic lineament that worked as the suture accommodating the shortening between the two plates is the Frontal Pennine Fault (FPF, Fig. 1). Even though the geological history of this belt is one of the best studied and well known in the world, the geodynamic evolution of the European slab, in terms of position and possible presence of slab break off, is still poorly understood.

All travel time tomographic studies identified at mantle depth the presence of seismic velocity heterogeneities interpreted as the European slab subducting beneath the Adria plate (e.g. Piromallo and Morelli, 2003; Lippitsch et al., 2003; Kissling et al., 2006; Giacomuzzi et al., 2011; Paffrath et al., 2021; Rappisi et al., 2022). However, the existence of possible slab detachments, windows or tears is debated, for instance beneath the Western Alps (e.g., Zhao et al., 2016). The first CIFALPS experiment (see Malusà et al., 2021 and references therein) clarified several points, starting from the first seismic evidence of subducted European continental lithosphere beneath the Adria lithosphere (Zhao et al., 2015), to a tomographic model with a continuous slab beneath this region (Zhao et al., 2016). In addition, recent seismic anisotropy analyses of the Western to the Central Alps shed additional light on potential discontinuities of the slabs, thanks to the possible mapping of mantle flows that would occur through them (Petrescu et al., 2020; Salimbeni et al., 2018).

The additional contribution of CIFALPS2, a temporary experiment deployed for 14 months from 2018 to 2019 (Zhao et al., 2018), on mantle seismic anisotropy mapping and interpretation, was expected to fill a gap in the northwestern part of the Alpine arc (red dots in Figure 1). Receiver function and ambient-noise tomography studies have underlined the north-south differences in the lithospheric structure along the belt strike (Paul et al., 2022). Therefore, there is a need for measuring additional seismic anisotropy in the mantle from CIFALPS2 data and to compare them to previous results.

In this study, we present the results of the analysis of data recorded by CIFALPS2, describing them in an integrated view with previous shear wave splitting measurements (SWS) to identify new features in the mantle and draw hypotheses on their origin.

## 2 Data and Methods

Data used for the analysis are the recordings at CIFALPS2 stations (Figure 1; Zhao et al., 2018; doi: 10.15778/RESIF.XT2018) of teleseismic earthquakes with a magnitude M>6.0, that occurred between June 2018 and December 2019 and are located at a distance interval from the network between 88° and 120°, typical to guarantee well isolated SKS phases in the waveforms. 80 to 150 events for each of the 56 temporary stations have been analyzed (Figure 1).

The entire SWS analysis has been conducted using the code SplitRacer (Reiss and Rümpker, 2017), based on the Silver and Chan (1991) method and thus on the minimization of the energy on the transverse component. Different filters have been
The European Alps originated in the late Cretaceous from the oblique subduction of the Alpine Tethys under the Adria microplate. The subduction evolved in a continental collision during the late Cenozoic (e.g., Handy et al., 2010, 2013 and references therein). In the Western Alps, the tectonic lineament that worked as the suture accommodating the shortening between the two plates is the Frontal Pennine Fault (FPF, Fig. 1). Even though the geological history of this belt is one of the best studied and well known in the world, the geodynamic evolution of the European slab, in terms of position and possible presence of slab break off, is still poorly understood.

All travel time tomographic studies identified at mantle depth the presence of seismic velocity heterogeneities interpreted as the European slab subducting beneath the Adria plate (e.g. Piromallo and Morelli, 2003; Lippitsch et al., 2003; Kissling et al., 2006; Giacomuzzi et al., 2011; Paffrath et al., 2021; Rappisi et al., 2022). However, the existence of possible slab detachments, windows or tears is debated, for instance beneath the Western Alps (e.g., Zhao et al., 2016). The first CIFALPS experiment (see Malusà et al., 2021 and references therein) clarified several points, starting from the first seismic evidence of subducted European continental lithosphere beneath the Adria lithosphere (Zhao et al., 2015), to a tomographic model with a continuous slab beneath this region (Zhao et al., 2016). In addition, recent seismic anisotropy analyses of the Western to the Central Alps shed additional light on potential discontinuities of the slabs, thanks to the possible mapping of mantle flows that would occur through them (Petrescu et al., 2020; Salimbeni et al., 2018).

The additional contribution of CIFALPS2, a temporary experiment deployed for 14 months from 2018 to 2019 (Zhao et al., 2018), on mantle seismic anisotropy mapping and interpretation, was expected to fill a gap in the northwestern part of the Alpine arc (red dots in Figure 1). Receiver function and ambient-noise tomography studies have underlined the north-south differences in the lithospheric structure along the belt strike (Paul et al., 2022). Therefore, there is a need for measuring additional seismic anisotropy in the mantle from CIFALPS2 data and to compare them to previous results.

In this study, we present the results of the analysis of data recorded by CIFALPS2, describing them in an integrated view with previous shear wave splitting measurements (SWS) to identify new features in the mantle and draw hypotheses on their origin.

## 2 Data and Methods

Data used for the analysis are the recordings at CIFALPS2 stations (Figure 1; Zhao et al., 2018; doi: 10.15778/RESIF.XT2018) of teleseismic earthquakes with a magnitude M>6.0, that occurred between June 2018 and December 2019 and are located at a distance interval from the network between 88° and 120°, typical to guarantee well isolated SKS phases in the waveforms. 80 to 150 events for each of the 56 temporary stations have been analyzed (Figure 1).

The entire SWS analysis has been conducted using the code SplitRacer (Reiss and Rümpker, 2017), based on the Silver and Chan (1991) method and thus on the minimization of the energy on the transverse component. Different filters have been
applied according to the amount of noise at the various sites. For most of them, located in the Alps or Ligurian mountains, a bandpass filter between 7 and 20 s worked well, while for instance sites in the Po Plain needed different choices, i.e. 5-30 s. The signal to noise ratio (SNR) was also used to avoid noisy waveforms; initially, the threshold was 3, but where the amount of events to be analyzed was scarce we decreased it down to 1.5, again mainly for sites located in or close to the Po Plain. It is worth noticing that SWS analysis recovers fast-velocity directions, assuming a single layer of horizontal anisotropy. Moreover, the depth at which this anisotropy is located is difficult to define, but it is classically assumed that most measured anisotropy is in the upper mantle (Savage, 1999). Thus, it is common to visualize any lateral variation by plotting results at the piercing point of the incident ray at 150 km depth (Figure 2a).

<figure><img src="image_3.png" /><figcaption>Figure 2 - a) Map of single SWS measurements for the study region, plotted at the location of the piercing point of the ray at 150 km depth. In red and orange are good and fair measurements from CIFALPS2 stations (red circles) respectively. In the background, in light blue, previous SWS measurements and stations; b) histogram of back azimuths of events used in the analysis.</figcaption></figure>

To improve the resolution of the data to be discussed and interpreted, splitting intensity (SI) measurements have been performed on the same CIFALPS2 recordings used for SWS measures (Table S2 in Supplementary Material). Splitting intensity is measured by projecting the transverse component on the radial component derivative; it is related to the variations of the amplitude of the transverse component with the back azimuth (Chevrot, 2000; Monteiller and Chevrot, 2010). A routine based on Kong et al. (2015) and Confal et al. (2023) was used to calculate splitting intensity values from waveforms with a cut-off of 15 s before and 30 s after the supposed SKS arrival. A dominant period of 12 s is used for the Wiener filtering. To filter out low quality waveforms in this automatic process a cross-correlation coefficient of |0.7| and
applied according to the amount of noise at the various sites. For most of them, located in the Alps or Ligurian mountains, a bandpass filter between 7 and 20 s worked well, while for instance sites in the Po Plain needed different choices, i.e. 5-30 s. The signal to noise ratio (SNR) was also used to avoid noisy waveforms; initially, the threshold was 3, but where the amount of events to be analyzed was scarce we decreased it down to 1.5, again mainly for sites located in or close to the Po Plain. It is worth noticing that SWS analysis recovers fast-velocity directions, assuming a single layer of horizontal anisotropy. Moreover, the depth at which this anisotropy is located is difficult to define, but it is classically assumed that most measured anisotropy is in the upper mantle (Savage, 1999). Thus, it is common to visualize any lateral variation by plotting results at the piercing point of the incident ray at 150 km depth (Figure 2a).

<figure><img src="image_3.png" /><figcaption>Figure 2 - a) Map of single SWS measurements for the study region, plotted at the location of the piercing point of the ray at 150 km depth. In red and orange are good and fair measurements from CIFALPS2 stations (red circles) respectively. In the background, in light blue, previous SWS measurements and stations; b) histogram of back azimuths of events used in the analysis.</figcaption></figure>

To improve the resolution of the data to be discussed and interpreted, splitting intensity (SI) measurements have been performed on the same CIFALPS2 recordings used for SWS measures (Table S2 in Supplementary Material). Splitting intensity is measured by projecting the transverse component on the radial component derivative; it is related to the variations of the amplitude of the transverse component with the back azimuth (Chevrot, 2000; Monteiller and Chevrot, 2010). A routine based on Kong et al. (2015) and Confal et al. (2023) was used to calculate splitting intensity values from waveforms with a cut-off of 15 s before and 30 s after the supposed SKS arrival. A dominant period of 12 s is used for the Wiener filtering. To filter out low quality waveforms in this automatic process a cross-correlation coefficient of |0.7| and
splitting intensities values and error threshold of |2.0| and 0.5 respectively were used (Baccheschi et al., under revision). With splitting measurements from at least four different 10° bins of back azimuths, the classical evaluation of the anisotropy parameters, i.e. the azimuth of the fast direction phi and the delay time dt, is obtained by fitting a sinusoid to the back azimuthal dependence of splitting intensity values (Chevrot, 2000). In particular, the sinusoid amplitude and phase give dt and phi respectively (see Figure S1 in Supplementary Material as an example).

### 3 Shear wave splitting results

From SWS analysis we obtained more than 400 pairs of splitting parameters (phi and dt) if we consider together good (170) and fair (241) results (Figure 2a, good in red, fair in orange; all results are available at https://osf.io/nqxk4, Pondrelli et al., 2023). The quality assignment is given following the SplitRacer criteria (Reiss and Rümpker, 2017), considering the visibility of the phase, the ellipticity of the initial particle motion and its linearity in the final stage, and the errors associated with phi and dt values. In Figure S2 of the Supplementary Material some measurement examples are shown. In addition, nearly 600 null measurements have been obtained (Figure S3 in Supplementary Material), where a null is considered when no split appears in the signal (i.e. no energy in the transversal component). This is due either to the absence of anisotropy or to the initial polarization being parallel to the fast or slow anisotropic direction. A high percentage of good and fair results were obtained for events with a NE back azimuth, so it should be taken into account that this direction is oversampled (Figure 2b).

This new dataset fills the region between the north-western external Alps and the Ligurian Sea. The anisotropy directions of no-nulls and nulls mostly agree with previous measurements (Figures 2 and 3, Figure S3 in Suppl. Material). Along the part of the transect crossing the Alps (transect AA' in Figure 3), NE-SW direction dominates in the internal part of the belt, between the FPF and the boundary of the Po Plain (see average values, red dots in Figure 3). In the western part of the AA' transect, measurements are more scattered, with a coexistence of NE-SW and NS to NNE-SSW directions. In the outer part of the Alpine belt and in the Bresse Graben, the prevailing directions are NS to NNE-SSW. Anisotropy in this region is weaker, but fast velocity directions remain constant toward the NW end of the transect, confirmed also by null measurements (Fig.S3 in Supplementary Material). These two patterns, one NE-SW parallel to the Alps strike and the other nearly NS, and their location along the transect are well visible in Figure 3, mainly following average values (reddish dots).
splitting intensities values and error threshold of |2.0| and 0.5 respectively were used (Baccheschi et al., under revision). With splitting measurements from at least four different 10° bins of back azimuths, the classical evaluation of the anisotropy parameters, i.e. the azimuth of the fast direction phi and the delay time dt, is obtained by fitting a sinusoid to the back azimuthal dependence of splitting intensity values (Chevrot, 2000). In particular, the sinusoid amplitude and phase give dt and phi respectively (see Figure S1 in Supplementary Material as an example).

### 3 Shear wave splitting results

From SWS analysis we obtained more than 400 pairs of splitting parameters (phi and dt) if we consider together good (170) and fair (241) results (Figure 2a, good in red, fair in orange; all results are available at https://osf.io/nqxk4, Pondrelli et al., 2023). The quality assignment is given following the SplitRacer criteria (Reiss and Rümpker, 2017), considering the visibility of the phase, the ellipticity of the initial particle motion and its linearity in the final stage, and the errors associated with phi and dt values. In Figure S2 of the Supplementary Material some measurement examples are shown. In addition, nearly 600 null measurements have been obtained (Figure S3 in Supplementary Material), where a null is considered when no split appears in the signal (i.e. no energy in the transversal component). This is due either to the absence of anisotropy or to the initial polarization being parallel to the fast or slow anisotropic direction. A high percentage of good and fair results were obtained for events with a NE back azimuth, so it should be taken into account that this direction is oversampled (Figure 2b).

This new dataset fills the region between the north-western external Alps and the Ligurian Sea. The anisotropy directions of no-nulls and nulls mostly agree with previous measurements (Figures 2 and 3, Figure S3 in Suppl. Material). Along the part of the transect crossing the Alps (transect AA' in Figure 3), NE-SW direction dominates in the internal part of the belt, between the FPF and the boundary of the Po Plain (see average values, red dots in Figure 3). In the western part of the AA' transect, measurements are more scattered, with a coexistence of NE-SW and NS to NNE-SSW directions. In the outer part of the Alpine belt and in the Bresse Graben, the prevailing directions are NS to NNE-SSW. Anisotropy in this region is weaker, but fast velocity directions remain constant toward the NW end of the transect, confirmed also by null measurements (Fig.S3 in Supplementary Material). These two patterns, one NE-SW parallel to the Alps strike and the other nearly NS, and their location along the transect are well visible in Figure 3, mainly following average values (reddish dots).
<figure><img src="image_4.png" /><figcaption>Figure 3 - Distribution of splitting parameters along the two sections of the CIFALPS2 region. For each section topography (upper), fast axes direction (middle) and delay time (lower) distributions along a swath box of 30 km are displayed. In the fast axes and delay time panels, every single measurement is represented by a square; the results of CIFALPS2 are blueish and color coded in agreement with the back azimuth of the events analyzed, while the results of previous works are represented by smaller empty squares. In the fast axes panel, circles represent the average values calculated using a basic circular arithmetic mean inside the swath box with 30-km-step increment; they are coloured in agreement with the spreading distribution around the mean value (R=0 distribution completely scattered, R=1 distribution completely aligned with the mean direction). In the delay time panels, red dots are the average value and its error, calculated with arithmetic mean and standard deviation. FPF = Frontal Penninic Fault; BG = Bresse Graben; PP = Po Plain.</figcaption></figure>

In the Po Plain, measurements appear scattered, similar to previous studies (e.g. Salimbeni et al., 2018; Petrescu et al., 2020; Figures 2 and 3). In the transition between the belt and the plain, the typical NE-SW Alpine direction prevails. Even if measurements are scarce, this direction appears also in the center of the Po Plain, together with a few NW-SE directions
<figure><img src="image_4.png" /><figcaption>Figure 3 - Distribution of splitting parameters along the two sections of the CIFALPS2 region. For each section topography (upper), fast axes direction (middle) and delay time (lower) distributions along a swath box of 30 km are displayed. In the fast axes and delay time panels, every single measurement is represented by a square; the results of CIFALPS2 are blueish and color coded in agreement with the back azimuth of the events analyzed, while the results of previous works are represented by smaller empty squares. In the fast axes panel, circles represent the average values calculated using a basic circular arithmetic mean inside the swath box with 30-km-step increment; they are coloured in agreement with the spreading distribution around the mean value (R=0 distribution completely scattered, R=1 distribution completely aligned with the mean direction). In the delay time panels, red dots are the average value and its error, calculated with arithmetic mean and standard deviation. FPF = Frontal Penninic Fault; BG = Bresse Graben; PP = Po Plain.</figcaption></figure>

In the Po Plain, measurements appear scattered, similar to previous studies (e.g. Salimbeni et al., 2018; Petrescu et al., 2020; Figures 2 and 3). In the transition between the belt and the plain, the typical NE-SW Alpine direction prevails. Even if measurements are scarce, this direction appears also in the center of the Po Plain, together with a few NW-SE directions
located southeast, close to the Ligurian Alps. However section BB' in Figure 3 shows that anisotropy directions are diverse, resulting in a strong dispersion of average values computed along the section. In the Ligurian part of the transect, several directions are detected, ENE-SSW and NW-SE (Apenninic) on the eastern side with a few weaker (lower dt) NNE-SSW measurements in the western side.

In general, we do not find any particular pattern in the delay time measurements. Average values computed along the sections (Figure 3) are mostly constant, around 1.5 s with a large range in single measurement values.

<figure><img src="image_5.png" /><figcaption>Figure 4 - Maps of average SWS measurements (red, yellow and blue) and anisotropy parameters obtained using SI measurements (green). In red average values for CIFALPS2 stations obtained with more than 3 measurements, in yellow averages obtained with less than 3 values; in light blue, average measurements from previous works. Dots represent stations.</figcaption></figure>

In Figure 4, average SWS values for each station are mapped together with anisotropy values obtained by splitting intensity (SI) measurements (Table S1 in Supplementary Material). The first observation is that average SWS and SI data are very similar, mainly at sites where the single measurement results are more homogeneous, while differences are present where
located southeast, close to the Ligurian Alps. However section BB' in Figure 3 shows that anisotropy directions are diverse, resulting in a strong dispersion of average values computed along the section. In the Ligurian part of the transect, several directions are detected, ENE-SSW and NW-SE (Apenninic) on the eastern side with a few weaker (lower dt) NNE-SSW measurements in the western side.

In general, we do not find any particular pattern in the delay time measurements. Average values computed along the sections (Figure 3) are mostly constant, around 1.5 s with a large range in single measurement values.

<figure><img src="image_5.png" /><figcaption>Figure 4 - Maps of average SWS measurements (red, yellow and blue) and anisotropy parameters obtained using SI measurements (green). In red average values for CIFALPS2 stations obtained with more than 3 measurements, in yellow averages obtained with less than 3 values; in light blue, average measurements from previous works. Dots represent stations.</figcaption></figure>

In Figure 4, average SWS values for each station are mapped together with anisotropy values obtained by splitting intensity (SI) measurements (Table S1 in Supplementary Material). The first observation is that average SWS and SI data are very similar, mainly at sites where the single measurement results are more homogeneous, while differences are present where
back azimuthal variations in fast direction are more evident, for instance in the southernmost part of the transect, in the Ligurian Mountains.

In general, it is clear that CIFALPS2 results are coherent with the average distribution of the anisotropy from previous measurements. The main Alpine pattern that follows the belt arc, here NE-SW, is represented in most of the averaged values for sites located in the transition between the Po Plain and the FPF. From the FPF to the NW endpoint of the transect, the main direction is close to NS, a direction that does not find an agreement with the orogen trend, still NE-SW.

In order to extract as much information as possible, we split the dataset into groups according to the main morphological features: a) External zone to the west of the FPF; b) Internal zone between the FPF and the western boundary of the Po plain, c) the Po Plain, and d) the Ligurian Alps. By plotting a rose diagram for all good and fair measurements obtained for the stations of each group, taking into account the influence of back azimuths, we get an overview of the main features and lateral changes of the anisotropy detected through these SWS measurements.

<figure><img src="image_6.png" /><figcaption>Figure 5 - Normalised rose diagrams produced for groups of stations along the profile, a) External zone with respect to FPF; b) Internal zone with respect to FPF; c) Po Plain; d) the Ligurian Alps. Grey rose diagrams: all data; Red: events with SW (a,b) or W (c,d) back azimuths, blue: events with NE (a,b) or E (c,d) back azimuths.</figcaption></figure>

In Figure 5, the rose diagrams show evident differences along the CIFALPS2 transect, somewhere underlined by different patterns with respect to back azimuths. Only in group a), in the External zone, the N-S dominant direction remains, regardless of the amount of measurements or separation into opposite back azimuths (red and blue rose diagrams). For all
back azimuthal variations in fast direction are more evident, for instance in the southernmost part of the transect, in the Ligurian Mountains.

In general, it is clear that CIFALPS2 results are coherent with the average distribution of the anisotropy from previous measurements. The main Alpine pattern that follows the belt arc, here NE-SW, is represented in most of the averaged values for sites located in the transition between the Po Plain and the FPF. From the FPF to the NW endpoint of the transect, the main direction is close to NS, a direction that does not find an agreement with the orogen trend, still NE-SW.

In order to extract as much information as possible, we split the dataset into groups according to the main morphological features: a) External zone to the west of the FPF; b) Internal zone between the FPF and the western boundary of the Po plain, c) the Po Plain, and d) the Ligurian Alps. By plotting a rose diagram for all good and fair measurements obtained for the stations of each group, taking into account the influence of back azimuths, we get an overview of the main features and lateral changes of the anisotropy detected through these SWS measurements.

<figure><img src="image_6.png" /><figcaption>Figure 5 - Normalised rose diagrams produced for groups of stations along the profile, a) External zone with respect to FPF; b) Internal zone with respect to FPF; c) Po Plain; d) the Ligurian Alps. Grey rose diagrams: all data; Red: events with SW (a,b) or W (c,d) back azimuths, blue: events with NE (a,b) or E (c,d) back azimuths.</figcaption></figure>

In Figure 5, the rose diagrams show evident differences along the CIFALPS2 transect, somewhere underlined by different patterns with respect to back azimuths. Only in group a), in the External zone, the N-S dominant direction remains, regardless of the amount of measurements or separation into opposite back azimuths (red and blue rose diagrams). For all
other groups, there is a clear difference depending on the back azimuth. In group b) the entire dataset shows a main trend in a NE-SW direction, in agreement with the Alps strike in this region; the same direction is observed for events coming from NE (blue rose diagram). On the contrary, measurements obtained from events with a SW back azimuth are dominated by a nearly NS direction, quite similar to that shown by group a). The entire dataset for the Po Plain (c) and the subgroup with east back azimuths show a clear bimodal distribution, with both the Alpine NE-SW and the Apenninic NW-SE typical directions; west back azimuth events show instead N to NNW directions, very similar to those obtained in previous work at the first CIFALPS transect (e.g. Salimbeni et al., 2018) and again similar to group a) and SW back azimuth measurements of group b). In the Ligurian Alps, in the group d), opposite back azimuths show different results. It is worth noticing that in the group d) the ENE-WSW direction is dominant both in the all data plot and in the east back azimuth plot, while the west back azimuth plot shows a wide dispersion with a NW-SE direction prevailing.

## 4 Discussion and hypotheses

CIFALPS2 SWS measurements clearly provide new information since they cover large areas where no seismic station is currently deployed or has been operating in the past. They allow us to fill data gaps and draw conclusions on the seismic anisotropy pattern and mantle deformation beneath the Western Alps.

A large-scale summary, prepared using the average values to avoid the scatter of single SWS measurements, is shown in Figure 6. ENE-WSW fast velocity directions, parallel to the strike of the Alps, are present from the Central to the Western Alps, in the transition between Po Plain and the belt, and terminate where the European slab and the belt are bending (double headed arrows in the inset sketch of Figure 6). At this point anisotropy directions do not strictly follow the chain strike, as they rather cut the main tectonic and morphological features. In the same zone (highlighted by the green circle A, Figure 6) converges a mantle flow that strikes from NE-SW to NS, in a coherent direction that however cuts the arcuate shape of the belt (yellow lines in the northern part of study region, Figure 6); this pattern is different from those described in previous anisotropy studies of the Western Alps (Barruol et al., 2004; Lucente et al., 2006; Barruol et al., 2011; Salimbeni et al., 2018) where anisotropy has always been described as rotating with the belt direction. The discrepancy between the direction of tectonic lines and the nearly NS fast velocity is in agreement with a deep source of the anisotropy, in the mantle. Link and Rumpker (2023) in a recent analysis found that a similar nearly NS pattern would be mostly located in a shallower part of the mantle, and in part also present in a lower layer. They consider it as a shallow asthenospheric contribution.

On the contrary, the NW-SE asthenospheric flow identified from SE France toward the Ligurian Sea, which culminates in a flow around the southern tip of the European slab, is in agreement with previous studies. This flow apparently originates beneath Central France, in the Massif Central region (Barruol et al., 2004), and it seems that the NS mantle flow merges with it (left green circle, Figure 6) and then, flowing around the southern tip of the European slab, moves to the Tyrrhenian Sea.

The overlap with the teleseismic travel time tomographic image at 150 km depth by Zhao et al. (2016) indicates that all these fast velocity directions correspond to mantle deformation below the European slab. The mantle close to the slab has a NE-
other groups, there is a clear difference depending on the back azimuth. In group b) the entire dataset shows a main trend in a NE-SW direction, in agreement with the Alps strike in this region; the same direction is observed for events coming from NE (blue rose diagram). On the contrary, measurements obtained from events with a SW back azimuth are dominated by a nearly NS direction, quite similar to that shown by group a). The entire dataset for the Po Plain (c) and the subgroup with east back azimuths show a clear bimodal distribution, with both the Alpine NE-SW and the Apenninic NW-SE typical directions; west back azimuth events show instead N to NNW directions, very similar to those obtained in previous work at the first CIFALPS transect (e.g. Salimbeni et al., 2018) and again similar to group a) and SW back azimuth measurements of group b). In the Ligurian Alps, in the group d), opposite back azimuths show different results. It is worth noticing that in the group d) the ENE-WSW direction is dominant both in the all data plot and in the east back azimuth plot, while the west back azimuth plot shows a wide dispersion with a NW-SE direction prevailing.

## 4 Discussion and hypotheses

CIFALPS2 SWS measurements clearly provide new information since they cover large areas where no seismic station is currently deployed or has been operating in the past. They allow us to fill data gaps and draw conclusions on the seismic anisotropy pattern and mantle deformation beneath the Western Alps.

A large-scale summary, prepared using the average values to avoid the scatter of single SWS measurements, is shown in Figure 6. ENE-WSW fast velocity directions, parallel to the strike of the Alps, are present from the Central to the Western Alps, in the transition between Po Plain and the belt, and terminate where the European slab and the belt are bending (double headed arrows in the inset sketch of Figure 6). At this point anisotropy directions do not strictly follow the chain strike, as they rather cut the main tectonic and morphological features. In the same zone (highlighted by the green circle A, Figure 6) converges a mantle flow that strikes from NE-SW to NS, in a coherent direction that however cuts the arcuate shape of the belt (yellow lines in the northern part of study region, Figure 6); this pattern is different from those described in previous anisotropy studies of the Western Alps (Barruol et al., 2004; Lucente et al., 2006; Barruol et al., 2011; Salimbeni et al., 2018) where anisotropy has always been described as rotating with the belt direction. The discrepancy between the direction of tectonic lines and the nearly NS fast velocity is in agreement with a deep source of the anisotropy, in the mantle. Link and Rumpker (2023) in a recent analysis found that a similar nearly NS pattern would be mostly located in a shallower part of the mantle, and in part also present in a lower layer. They consider it as a shallow asthenospheric contribution.

On the contrary, the NW-SE asthenospheric flow identified from SE France toward the Ligurian Sea, which culminates in a flow around the southern tip of the European slab, is in agreement with previous studies. This flow apparently originates beneath Central France, in the Massif Central region (Barruol et al., 2004), and it seems that the NS mantle flow merges with it (left green circle, Figure 6) and then, flowing around the southern tip of the European slab, moves to the Tyrrhenian Sea.

The overlap with the teleseismic travel time tomographic image at 150 km depth by Zhao et al. (2016) indicates that all these fast velocity directions correspond to mantle deformation below the European slab. The mantle close to the slab has a NE-
SW seismic anisotropy direction because it is deformed by the slab steepening. This feature is visible only in a narrow stripe, in the transition between Western and Central Alps, probably because in this place the slab steepening was favored by the plate motion direction (e.g. Adria anticlockwise rotation, Figure 6). Moving north, in the outer part of the belt, the anisotropy directions are those of an undeformed mantle, substantially parallel to the APM direction (orange arrows in Figure 6). Moving (south)westward, the deflection from NE-SW to NS may be related to the mantle being squeezed by the retreating European slab and moved toward the south where the retreat process was weakening or probably already ended. Indeed, in the Western Alps, the arcuate shape of the trench and slab, together with the (upper) Adria plate rotating in anticlockwise direction with a rotation pole more or less located in the western Po Plain (e.g. Serpelloni et al., 2016; Le Breton et al., 2017), may have been less favorable to the retreat of the slab. Such differences between the northern and southern Western Alps, in particular in the European Moho geometry, have been also identified by comparing CIFALPS and CIFALPS2 receiver function sections (Paul et al., 2022). The European Moho is strongly dipping down to ~75 km depth along CIFALPS section, while its dip is much weaker in CIFALPS2 section. Moreover, the absence of slab retreat in the Western Alps has been described also by Malusà et al. (2015), studying the mechanisms for exhumation of (U)HP terranes along the Cenozoic Adria-Europe plate boundary.

<figure><img src="image_7.png" /></figure>

Figure 6 - Left: Zhao et al. (2016) tomography (dVp in %) - layer at 150 km depth - overlapped with average SWS measurements: results from CIFALPS2 data in light grey and from previous studies in black. Orange arrows on top represent the absolute European plate motion from GSRM v2.1 model (calculated from Plate Motion Calculator | Software | GAGE). Yellow traces simulate average mantle deformation directions and green areas A and B highlight points where directions converge. Right: a sketch of this part of alpine subduction, where, in light yellow is the European plate, the double-headed arrows represent the only anisotropy of the study region related to European slab retreat (whose direction is represented by the red arrows), thin black lines are the anisotropy parallel to the APM, the thick black arrow represents the SE France asthenospheric mantle flow, the purple circular arrow represents the Adria plate rotation.
SW seismic anisotropy direction because it is deformed by the slab steepening. This feature is visible only in a narrow stripe, in the transition between Western and Central Alps, probably because in this place the slab steepening was favored by the plate motion direction (e.g. Adria anticlockwise rotation, Figure 6). Moving north, in the outer part of the belt, the anisotropy directions are those of an undeformed mantle, substantially parallel to the APM direction (orange arrows in Figure 6). Moving (south)westward, the deflection from NE-SW to NS may be related to the mantle being squeezed by the retreating European slab and moved toward the south where the retreat process was weakening or probably already ended. Indeed, in the Western Alps, the arcuate shape of the trench and slab, together with the (upper) Adria plate rotating in anticlockwise direction with a rotation pole more or less located in the western Po Plain (e.g. Serpelloni et al., 2016; Le Breton et al., 2017), may have been less favorable to the retreat of the slab. Such differences between the northern and southern Western Alps, in particular in the European Moho geometry, have been also identified by comparing CIFALPS and CIFALPS2 receiver function sections (Paul et al., 2022). The European Moho is strongly dipping down to ~75 km depth along CIFALPS section, while its dip is much weaker in CIFALPS2 section. Moreover, the absence of slab retreat in the Western Alps has been described also by Malusà et al. (2015), studying the mechanisms for exhumation of (U)HP terranes along the Cenozoic Adria-Europe plate boundary.

<figure><img src="image_7.png" /></figure>

Figure 6 - Left: Zhao et al. (2016) tomography (dVp in %) - layer at 150 km depth - overlapped with average SWS measurements: results from CIFALPS2 data in light grey and from previous studies in black. Orange arrows on top represent the absolute European plate motion from GSRM v2.1 model (calculated from Plate Motion Calculator | Software | GAGE). Yellow traces simulate average mantle deformation directions and green areas A and B highlight points where directions converge. Right: a sketch of this part of alpine subduction, where, in light yellow is the European plate, the double-headed arrows represent the only anisotropy of the study region related to European slab retreat (whose direction is represented by the red arrows), thin black lines are the anisotropy parallel to the APM, the thick black arrow represents the SE France asthenospheric mantle flow, the purple circular arrow represents the Adria plate rotation.
The mantle flow from beneath the Central Alps converges with the asthenospheric flow coming from SE France where Western Alps subduction, precisely a continental subduction (Malusà et al., 2021; Paul et al., 2022) runs out. In the region, the European slab is almost vertical (Zhao et al., 2016) and is no more affected by slab retreat. Substantially a differential behaviour of slab movement along the chain is at the base of the anisotropy pattern variation from Central to the entire Western Alps.

In the Po Plain, the image is patchy and complex, with very few good quality measurements and some really low values of dt. This does not necessarily mean that the mantle is isotropic, but that anisotropy may be multilayered or fast velocity directions are vertical.

It is worth noting that beneath the Po basin, the European slab is nearly vertical (e.g., Zhao et al. 2016; Paffrath et al., 2021) and the space for the mantle above is really narrow. With such a complex mantle structure, it is not surprising that a unique and significant pattern of anisotropy cannot be identified. In the Ligurian Alps, our measurements show different orientations east and west of the transect with, in general, a prevailing EW direction (Figure 5d), but a minor NNE-SSW set of measurements from the western back azimuths. These results are certainly intriguing but not sufficient to support any new hypothesis.

## 5 Conclusion

New data collected by the CIFALPS2 project clearly fill the gap that has forced the interpolation between more sparse information in the past. The pattern of seismic anisotropy shown here, is not entirely parallel to the belt, as it is in the Central Alps and in the transition to the Western Alps, up to the point this portion of belt is arcuate. Indeed, in the central part of the study region, a nearly NS pattern coming from central Europe cuts all principal morphologic features of the belt. It converges with the part of the mantle deformed by the retreating slab in the point where the retreat is less favored or ended. The arcuate shape of the belt and of the slab, added to the Adria plate rotation with respect to Eurasia around a pole here particularly close to the boundary, reduce the effectiveness of a slab retreat. The NS mantle flow is then interpreted as the European mantle moving south to merge with the large asthenospheric flow that from beneath Central France moves toward the Tyrrhenian Sea turning around the southern tip of the European slab.

## Data availability

All shear wave splitting measurements from this work have been included and are available in the Italian shear wave splitting collection https://osf.io/nqxk4 (Pondrelli et al., 2023).
The mantle flow from beneath the Central Alps converges with the asthenospheric flow coming from SE France where Western Alps subduction, precisely a continental subduction (Malusà et al., 2021; Paul et al., 2022) runs out. In the region, the European slab is almost vertical (Zhao et al., 2016) and is no more affected by slab retreat. Substantially a differential behaviour of slab movement along the chain is at the base of the anisotropy pattern variation from Central to the entire Western Alps.

In the Po Plain, the image is patchy and complex, with very few good quality measurements and some really low values of dt. This does not necessarily mean that the mantle is isotropic, but that anisotropy may be multilayered or fast velocity directions are vertical.

It is worth noting that beneath the Po basin, the European slab is nearly vertical (e.g., Zhao et al. 2016; Paffrath et al., 2021) and the space for the mantle above is really narrow. With such a complex mantle structure, it is not surprising that a unique and significant pattern of anisotropy cannot be identified. In the Ligurian Alps, our measurements show different orientations east and west of the transect with, in general, a prevailing EW direction (Figure 5d), but a minor NNE-SSW set of measurements from the western back azimuths. These results are certainly intriguing but not sufficient to support any new hypothesis.

# 5 Conclusion

New data collected by the CIFALPS2 project clearly fill the gap that has forced the interpolation between more sparse information in the past. The pattern of seismic anisotropy shown here, is not entirely parallel to the belt, as it is in the Central Alps and in the transition to the Western Alps, up to the point this portion of belt is arcuate. Indeed, in the central part of the study region, a nearly NS pattern coming from central Europe cuts all principal morphologic features of the belt. It converges with the part of the mantle deformed by the retreating slab in the point where the retreat is less favored or ended. The arcuate shape of the belt and of the slab, added to the Adria plate rotation with respect to Eurasia around a pole here particularly close to the boundary, reduce the effectiveness of a slab retreat. The NS mantle flow is then interpreted as the European mantle moving south to merge with the large asthenospheric flow that from beneath Central France moves toward the Tyrrhenian Sea turning around the southern tip of the European slab.

# Data availability

All shear wave splitting measurements from this work have been included and are available in the Italian shear wave splitting collection https://osf.io/nqxk4 (Pondrelli et al., 2023).
**Author contribution**

PS, SS and CJM made measurements and analyses. PS prepared the manuscript with the contribution of all co-authors. All authors designed the CIFALPS2 experiment and carried it out.

**Competing interests**

The authors declare that they have no conflict of interest.

**Acknowledgements**

This research is funded by the National Natural Science Foundation of China and by NEWTON (NEw Window inTO Earth's iNterior), ERC StG funded project (grant ID:758199).

**References**

Barruol, G., Bonnin, M., Pedersen, H., Bokelmann, G., Tiberi, C.: Belt-parallel mantle flow beneath a halted continental collision: the Western Alps, Earth Planet. Sci. Lett. 302, 429–438, 2011.

Barruol, G., Deschamps, A., Coutant, O.: 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, Tectonophys., 394, 125–138, 2004.

Barruol, G. and M. Granet: A Tertiary asthenospheric flow beneath the southern French Massif Central related to the west Mediterranean extension evidenced by upper mantle seismic anisotropy, Earth Plan. Sci. Lett., 202, 31–47, 10.1016/S0012-821X(02)00752-5, 2002.

Chevrot, S.: Multichannel analysis of shear wave splitting, J. Geophys. Res. Solid Earth, 105(B9), 21579-21590, 2000.

Confal, J. M., Baccheschi, P., Pondrelli, S., Karakostas, F., VanderBeek, B. P., Huang, Z., Faccenda, M.: Reproducing Complex Anisotropy Patterns at Subduction Zones from Splitting Intensity Analysis and Anisotropy Tomography, Geophys. J. Int., 235, 1725–1735, https://doi.org/10.1093/gji/ggad329, 2023.

Crampin, S. and Peacock S.: A review of the current understanding of seismic shear-wave splitting in the Earth’s crust and common fallacies in interpretation, Wave Mot., 45 (2008), pp. 675-722, 10.1016/j.wavemoti.2008.01.003, 2008.

Díaz, J., Gil, A., and Gallart, J.: Uppermost mantle seismic velocity and anisotropy in the Euro-Mediterranean region from Pn and Sn tomography, Geophys. J. Int., 192(1), 310-325, 2013.

Monteiller, V., and Chevrot, S.: How to make robust splitting measurements for single-station analysis and three-dimensional imaging of seismic anisotropy, Geophys. J. Int., 182(1), 311-328, 2010.

Giacomuzzi, G., Chiarabba, C., and De Gori, P.: Linking the Alps and Apennines subduction systems: New constraints revealed by high-resolution teleseismic tomography, Earth and Plan. Science Lett., 301(3–4), 531–543, 2011.
**Author contribution**

PS, SS and CJM made measurements and analyses. PS prepared the manuscript with the contribution of all co-authors. All authors designed the CIFALPS2 experiment and carried it out.

**Competing interests**

The authors declare that they have no conflict of interest.

**Acknowledgements**

This research is funded by the National Natural Science Foundation of China and by NEWTON (NEw Window inTO Earth's iNterior), ERC StG funded project (grant ID:758199).

**References**

Barruol, G., Bonnin, M., Pedersen, H., Bokelmann, G., Tiberi, C.: Belt-parallel mantle flow beneath a halted continental collision: the Western Alps, Earth Planet. Sci. Lett. 302, 429–438, 2011.

Barruol, G., Deschamps, A., Coutant, O.: 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, Tectonophys., 394, 125–138, 2004.

Barruol, G. and M. Granet: A Tertiary asthenospheric flow beneath the southern French Massif Central related to the west Mediterranean extension evidenced by upper mantle seismic anisotropy, Earth Plan. Sci. Lett., 202, 31–47, 10.1016/S0012-821X(02)00752-5, 2002.

Chevrot, S.: Multichannel analysis of shear wave splitting, J. Geophys. Res. Solid Earth, 105(B9), 21579-21590, 2000.

Confal, J. M., Baccheschi, P., Pondrelli, S., Karakostas, F., VanderBeek, B. P., Huang, Z., Faccenda, M.: Reproducing Complex Anisotropy Patterns at Subduction Zones from Splitting Intensity Analysis and Anisotropy Tomography, Geophys. J. Int., 235, 1725–1735, https://doi.org/10.1093/gji/ggad329, 2023.

Crampin, S. and Peacock S.: A review of the current understanding of seismic shear-wave splitting in the Earth’s crust and common fallacies in interpretation, Wave Mot., 45 (2008), pp. 675-722, 10.1016/j.wavemoti.2008.01.003, 2008.

Díaz, J., Gil, A., and Gallart, J.: Uppermost mantle seismic velocity and anisotropy in the Euro-Mediterranean region from Pn and Sn tomography, Geophys. J. Int., 192(1), 310-325, 2013.

Monteiller, V., and Chevrot, S.: How to make robust splitting measurements for single-station analysis and three-dimensional imaging of seismic anisotropy, Geophys. J. Int., 182(1), 311-328, 2010.

Giacomuzzi, G., Chiarabba, C., and De Gori, P.: Linking the Alps and Apennines subduction systems: New constraints revealed by high-resolution teleseismic tomography, Earth and Plan. Science Lett., 301(3–4), 531–543, 2011.
Hetényi G., I. Molinari, J. Clinton et al.: The AlpArray Seismic Network: a large-scale European experiment to image the Alpine orogeny, Surv. Geophysics, 39, 1009-1033. doi:10.1007/s10712-018-9472-4, 2018.

Kissling, E., Schmid, S. M., Lippitsch, R., Ansorge, J., and Fügenschuh, B.: Lithosphere structure and tectonic evolution of the Alpine arc: New evidence from high-resolution teleseismic tomography, Geological Society, London, Memoirs, 32, 129–145, 2006.

Kong, F., Gao, S. S., and Liu, K. H.: Applicability of the multiple-event stacking technique for shear-wave splitting analysis, Bull. Seism. Soc. Am., 105(6), 3156-3166, 2015.

Le Breton, E., Handy, M. R., Molli, G., and Ustaszewski, K.: Post-20 Ma motion of the Adriatic plate: New constraints from surrounding Orogens and implications for crust-mantle decoupling, Tectonics, 36, 3135–3154. doi:10.1002/2016TC004443, 2017.

Lippitsch, R., Kissling, E., and Ansorge, J.: Upper mantle structure beneath the Alpine orogen from high-resolution teleseismic tomography, J. Geophys. Res.: Solid Earth, 108(B8). https://doi.org/10.1029/2002JB002016, 2003.

Long, M.: Constraints on subduction geodynamics from seismic anisotropy. Rev. Geophys., 51, 76–112, https://doi.org/10.1002/rog.20008, 2013.

Long, M.D., and Silver, P. G.: Shear Wave Splitting and Mantle Anisotropy: Measurements, Interpretations, and New Directions, Surv Geophys 30, 407–461, https://doi.org/10.1007/s10712-009-9075-1, 2009.

Malusà, M. G., C. Faccenna, S. L. Baldwin, P. G. Fitzgerald, F. Rossetti, M. L. Balestrieri, M. Danišík, A. Ellero, G. Ottria, and C. Piromallo: Contrasting styles of (U)HP rock exhumation along the Cenozoic Adria-Europe plate boundary (Western Alps, Calabria, Corsica), Geochem. Geophys. Geosyst., 16, 1786–1824, doi:10.1002/2015GC005767, 2015.

Malusà, M. G., Guillot, S., Zhao, L., Paul, A., Solarino, S., Dumont, T., Schwartz, S., Aubert, C., Baccheschi, P., Eva, E., Lu, Y., Lyu, C., Pondrelli, S., Salimbeni, S., Sun, W., & Yuan, H.: The deep structure of the Alps based on the CIFALPS seismic experiment: A synthesis, Geochem. Geophys. Geosyst., 22(3), e2020GC009466, 2021.

Okaya, D., S.S. Vel, W.J. Song, S.E. Johnson: Modification of crustal seismic anisotropy by geological structures (“structural geometric anisotropy”), Geosphere, 15, pp. 146-170, 10.1130/GES01655, 2018.

Paffrath, M., Friederich, W., Schmid, S. M., Handy, M. R., and AlpArray and AlpArray-Swath D Working Group: Imaging structure and geometry of slabs in the greater Alpine area–A P-wave travel-time tomography using AlpArray Seismic Network data, Solid Earth, 12, 2671-2702, doi: 10.5194/se-12-2671-2021, 2021.

Paul, A., M.G. Malusà, S. Solarino, S. Salimbeni, E. Eva, A. Nouibat, S. Pondrelli, C. Aubert, T. Dumont, S. Guillot, S. Schwartz, L. Zhao: Along-strike variations in the fossil subduction zone of the Western Alps revealed by the CIFALPS seismic experiments and their implications for exhumation of (ultra-) high-pressure rocks, Earth Plan. Science Lett., 598, 117843, https://doi.org/10.1016/j.epsl.2022.117843, 2022.

Piromallo, C., and Morelli, A.: P wave tomography of the mantle under the Alpine-Mediterranean area, J. Geophys. Res.: Solid Earth, 108(B2). https://doi.org/10.1029/2002JB001757, 2003.

Pondrelli S., S. Salimbeni, P. Baccheschi, J. M. Confal, L. Margheriti (2023). Peeking inside the mantle structure beneath the Italian region through SKS shear wave splitting anisotropy: a review, Ann. Geophys., 66, 2, doi:1044.01/ag-8872.
Hetényi G., I. Molinari, J. Clinton et al.: The AlpArray Seismic Network: a large-scale European experiment to image the Alpine orogeny, Surv. Geophysics, 39, 1009-1033. doi:10.1007/s10712-018-9472-4, 2018.

Kissling, E., Schmid, S. M., Lippitsch, R., Ansorge, J., and Fügenschuh, B.: Lithosphere structure and tectonic evolution of the Alpine arc: New evidence from high-resolution teleseismic tomography, Geological Society, London, Memoirs, 32, 129–145, 2006.

Kong, F., Gao, S. S., and Liu, K. H.: Applicability of the multiple-event stacking technique for shear-wave splitting analysis, Bull. Seism. Soc. Am., 105(6), 3156-3166, 2015.

Le Breton, E., Handy, M. R., Molli, G., and Ustaszewski, K.: Post-20 Ma motion of the Adriatic plate: New constraints from surrounding Orogens and implications for crust-mantle decoupling, Tectonics, 36, 3135–3154. doi:10.1002/2016TC004443, 2017.

Lippitsch, R., Kissling, E., and Ansorge, J.: Upper mantle structure beneath the Alpine orogen from high-resolution teleseismic tomography, J. Geophys. Res.: Solid Earth, 108(B8). https://doi.org/10.1029/2002JB002016, 2003.

Long, M.: Constraints on subduction geodynamics from seismic anisotropy. Rev. Geophys., 51, 76–112, https://doi.org/10.1002/rog.20008, 2013.

Long, M.D., and Silver, P. G.: Shear Wave Splitting and Mantle Anisotropy: Measurements, Interpretations, and New Directions, Surv Geophys 30, 407–461, https://doi.org/10.1007/s10712-009-9075-1, 2009.

Malusà, M. G., C. Faccenna, S. L. Baldwin, P. G. Fitzgerald, F. Rossetti, M. L. Balestrieri, M. Danišík, A. Ellero, G. Ottria, and C. Piromallo: Contrasting styles of (U)HP rock exhumation along the Cenozoic Adria-Europe plate boundary (Western Alps, Calabria, Corsica), Geochem. Geophys. Geosyst., 16, 1786–1824, doi:10.1002/2015GC005767, 2015.

Malusà, M. G., Guillot, S., Zhao, L., Paul, A., Solarino, S., Dumont, T., Schwartz, S., Aubert, C., Baccheschi, P., Eva, E., Lu, Y., Lyu, C., Pondrelli, S., Salimbeni, S., Sun, W., & Yuan, H.: The deep structure of the Alps based on the CIFALPS seismic experiment: A synthesis, Geochem. Geophys. Geosyst., 22(3), e2020GC009466, 2021.

Okaya, D., S.S. Vel, W.J. Song, S.E. Johnson: Modification of crustal seismic anisotropy by geological structures (“structural geometric anisotropy”), Geosphere, 15, pp. 146-170, 10.1130/GES01655, 2018.

Paffrath, M., Friederich, W., Schmid, S. M., Handy, M. R., and AlpArray and AlpArray-Swath D Working Group: Imaging structure and geometry of slabs in the greater Alpine area–A P-wave travel-time tomography using AlpArray Seismic Network data, Solid Earth, 12, 2671-2702, doi: 10.5194/se-12-2671-2021, 2021.

Paul, A., M.G. Malusà, S. Solarino, S. Salimbeni, E. Eva, A. Nouibat, S. Pondrelli, C. Aubert, T. Dumont, S. Guillot, S. Schwartz, L. Zhao: Along-strike variations in the fossil subduction zone of the Western Alps revealed by the CIFALPS seismic experiments and their implications for exhumation of (ultra-) high-pressure rocks, Earth Plan. Science Lett., 598, 117843, https://doi.org/10.1016/j.epsl.2022.117843, 2022.

Piromallo, C., and Morelli, A.: P wave tomography of the mantle under the Alpine-Mediterranean area, J. Geophys. Res.: Solid Earth, 108(B2). https://doi.org/10.1029/2002JB001757, 2003.

Pondrelli S., S. Salimbeni, P. Baccheschi, J. M. Confal, L. Margheriti (2023). Peeking inside the mantle structure beneath the Italian region through SKS shear wave splitting anisotropy: a review, Ann. Geophys., 66, 2, doi:1044.01/ag-8872.
Rappisi, F., VanderBeek, B. P., Faccenda, M., Morelli, A., Molinari, I.: Slab geometry and upper mantle flow patterns in the Central Mediterranean from 3D anisotropic P-wave tomography, J. Geophys. Res.: Solid Earth, 127, e2021JB023488, doi: 10.1029/2021JB023488, 2022.

Reiss, M. C. and G. Rümpker: SplitRacer: MATLAB Code and GUI for Semiautomated Analysis and Interpretation of Teleseismic Shear-Wave Splitting, Seism. Res. Lett., 88, 2A, 392-409, doi: https://doi.org/10.1785/0220160191, 2017.

Salimbeni, S., M. G. Malusà, L. Zhao, S. Guillot, S. Pondrelli, L. Margheriti, A. Paul, S. Solarino, C. Aubert, T. Dumont, S. Schwartz, Q. Wang, X. Xu, T. Zheng and R. Zhu: Active and fossil mantle flows in the western Alpine region unravelled by seismic anisotropy analysis and high-resolution P wave tomography, Tectonophys., 731-732, 35-47, doi:10.1016/j.tecto.2018.03.002, 2018.

Serpelloni, E., G. Vannucci, L. Anderlini, and R.A. Bennett: Kinematics, seismotectonics and seismic potential of the eastern sector of the European Alps from GPS and seismic deformation data, Tectonophys., https://doi.org/10.1016/j.tecto.2016.09.026, 2016.

Silver, P. G. and W. W. Chan: Shear wave splitting and subcontinental mantle deformation, J. Geophys. Res.: Solid Earth, 96, B10, 16429-16454, 1991.

Zhao, L., Paul, A., Guillot, S., Solarino, S., Malusà, M.G., Zheng, T.Y., Aubert, C., Salimbeni, S., Dumont, T., Schwartz, S., Zhu, R.X., Wang, Q.C.: First seismic evidence for continental subduction beneath the Western Alps, Geology 43, 815-818, 2015.

Zhao, L., Paul, A., Malusà, M. G., Xu, X., Zheng, T., Solarino, S., Guillot, S. Schwartz, S., Dumont, T., Salimbeni, S., Aubert, C., Pondrelli, S., Wang, Q., and Zhu, R.: Continuity of the Alpine slab unraveled by high-resolution P wave tomography, J. Geophys. Res., 121, 8720-8737, 2016.

Zhao, L., Paul, A., Solarino, S., RESIF: Seismic network YP: CIFALPS temporary experiment (China-Italy-France Alps seismic transect) [Data set]. RESIF - Réseau Sismologique et géodésique Français. https://doi.org/10.15778/RESIF.YP2012, 2016.

Zhao, L., Paul, A., Solarino, S., RESIF: Seismic network XT: CIFALPS-2 temporary experiment (China-Italy-France Alps seismic transect #2 [Data set]. RESIF - Réseau Sismologique et Géodésique Franc,ais. https://doi.org/10.15778/RESIF.XT2018, 2018.

Zhao D., Liu X., Wang Z. and Gou T.: Seismic Anisotropy Tomography and Mantle Dynamics, Surv. Geophys., 44, 947-982, https://doi.org/10.1007/s10712-022-09764-71, 2023.
Rappisi, F., VanderBeek, B. P., Faccenda, M., Morelli, A., Molinari, I.: Slab geometry and upper mantle flow patterns in the Central Mediterranean from 3D anisotropic P-wave tomography, J. Geophys. Res.: Solid Earth, 127, e2021JB023488. doi: 10.1029/2021JB023488, 2022.

Reiss, M. C. and G. Rümpker: SplitRacer: MATLAB Code and GUI for Semiautomated Analysis and Interpretation of Teleseismic Shear-Wave Splitting, Seism. Res. Lett., 88, 2A, 392-409, doi: https://doi.org/10.1785/0220160191, 2017.

Salimbeni, S., M. G. Malusà, L. Zhao, S. Guillot, S. Pondrelli, L. Margheriti, A. Paul, S. Solarino, C. Aubert, T. Dumont, S. Schwartz, Q. Wang, X. Xu, T. Zheng and R. Zhu: Active and fossil mantle flows in the western Alpine region unravelled by seismic anisotropy analysis and high-resolution P wave tomography, Tectonophys., 731-732, 35-47, doi:10.1016/j.tecto.2018.03.002, 2018.

Serpelloni, E., G. Vannucci, L. Anderlini, and R.A. Bennett: Kinematics, seismotectonics and seismic potential of the eastern sector of the European Alps from GPS and seismic deformation data, Tectonophys., https://doi.org/10.1016/j.tecto.2016.09.026, 2016.

Silver, P. G. and W. W. Chan: Shear wave splitting and subcontinental mantle deformation, J. Geophys. Res.: Solid Earth, 96, B10, 16429-16454, 1991.

Zhao, L., Paul, A., Guillot, S., Solarino, S., Malusà, M.G., Zheng, T.Y., Aubert, C., Salimbeni, S., Dumont, T., Schwartz, S., Zhu, R.X., Wang, Q.C.: First seismic evidence for continental subduction beneath the Western Alps, Geology 43, 815-818, 2015.

Zhao, L., Paul, A., Malusà, M. G., Xu, X., Zheng, T., Solarino, S., Guillot, S. Schwartz, S., Dumont, T., Salimbeni, S., Aubert, C., Pondrelli, S., Wang, Q., and Zhu, R.: Continuity of the Alpine slab unraveled by high-resolution P wave tomography, J. Geophys. Res., 121, 8720-8737, 2016.

Zhao, L., Paul, A., Solarino, S., RESIF: Seismic network YP: CIFALPS temporary experiment (China-Italy-France Alps seismic transect) [Data set]. RESIF - Réseau Sismologique et géodésique Français. https://doi.org/10.15778/RESIF.YP2012, 2016.

Zhao, L., Paul, A., Solarino, S., RESIF: Seismic network XT: CIFALPS-2 temporary experiment (China-Italy-France Alps seismic transect #2 [Data set]. RESIF - Réseau Sismologique et Géodésique Franc,ais. https://doi.org/10.15778/RESIF. XT2018, 2018.

Zhao D., Liu X., Wang Z. and Gou T.: Seismic Anisotropy Tomography and Mantle Dynamics, Surv. Geophys., 44, 947-982, https://doi.org/10.1007/s10712-022-09764-71, 2023.