Teleseismic Receiver Functions Research Articles

Crustal silica content of East China: A seismological perspective and its significance

East China has experienced multiple periods of tectonic movements, which have contributed to the composition and rheological properties of its present crust. Estimating the composition of the crust is crucial for understanding the tectonic processes. Based on the teleseismic receiver functions with data from the National Seismic Network of China, we applied the H-κ-c method to obtain the crustal bulk VP/VS ratio and to constrain the SiO2 content in the crust. We estimated the SiO2 content to range from 50.89 wt% to 73.51 wt%, with an average value of 65.87 wt%, indicating a predominant felsic composition of the East China's crust. Our study suggests that the North-South Gravity Lineament (NSGL) is an approximate delimitator of the felsic and mafic crust in East China, hinting at a widespread deficiency of mafic lower crust in the east of the NSGL. The mafic crust is extensively distributed in the Taihang orogenic region (TSR) and WuLingshan gravity gradient belts (WLG), particularly in the Datong volcanic area, which manifests the mantle materials intraplating. The scatteredly distributed mafic crust at the east of the NSGL is mainly concentrated in the southeast coast of China and in the intersection region of the Tanlu Fault (TLF) and Sulu region. In Sulu region, the TLF may primarily provide a channel for the thermal intrusion from the underlying mantle lithosphere, which has increased the mafic content. The Pacific/Philippine Sea plate subduction has triggered a significant amount of crust-mantle material exchange below southeast China that resulted in a high degree of mafic crustal composition.

Inversion of the Deep Structure of the Northern Margin of the Ordos Block Using Teleseismic Receiver Functions and Gravity Data

Inversion of the Deep Structure of the Northern Margin of the Ordos Block Using Teleseismic Receiver Functions and Gravity Data

INSIGHTS INTO THE MOHO DEPTH VARIATIONS AND CRUSTAL CHARACTERISTICS IN THE GUELMA–CONSTANTINE BASIN, NORTHEASTERN ALGERIA: A SHORT-PERIOD SEISMIC-RECEIVER-FUNCTION PERSPECTIVE

Teleseismic receiver functions (RF) were extracted from data collected at eight short-period, three-component seismic recording stations over the Guelma–Constantine Basin, northeastern Algeria, to improve the understanding of crustal structure and geodynamic processes. The H-κ stacking method was used to determine the Moho depths and average vP/vS ratios at each station. Careful linear inversion of RF was performed to determine the most appropriate average shear-wave and P-wave velocity profiles at each site. Both methods have yielded highly congruent results, with Moho depths showing robust correlations with previous seismological and geophysical studies. The previously observed pattern of the increasing Moho depth from north to south in the Tell Atlas has been confirmed. Furthermore, the identified transitional nature of the Moho in the Constantine Basin is consistent with a recent study. In addition, we identify a low-velocity zone (LVZ) at approximately 20 km depth within the southern Guelma Basin, confirming previous results in the Constantine Basin and suggesting an eastward elongation of the LVZ, at least into the southern periphery of the Guelma Basin. Examination of data from the northern tip of the Hammam Debbagh–Roknia NW–SE fault, the western boundary of the Guelma pull-apart basin, revealed a shallow Moho depth (22 km), less than the basin average depth of 25 km. The LVZ observed in the lower crust (12 km) suggests the presence of partial melts, consistent with gravimetric and chemical analyses of hydrothermal sources in the area. The extensional tectonic activity along this boundary, coupled with the low-viscosity zone and low average vP/vS ratio, is potentially associated with delamination processes. The effectiveness of our approach underscores its potential as a viable alternative or complementary method for investigating variations in the Moho depth.

Crustal structures and anisotropy in southeastern China: Continental collision and intraplate earthquakes

This study investigates the role of crustal structures and anisotropy in southeastern China regarding intraplate earthquakes within the context of concentrated and perturbed tectonic stresses. The region is characterized by a complex multiphase tectonic history,including the ongoing collision between the Philippine Plate and the South China Sea margin. Utilizing the teleseismic receiver function method, we offer a comprehensive depiction of the crustal structures and layered anisotropy in Guangdong and Fujian, two earthquake-prone areas. Our results reveal that the crustal anisotropy is stratified into two layers, with a slow axis perpendicular to the orogenic belt in the upper-middle crust, whereas the fast axis aligns perpendicularly at lower crust levels. The lateral heterogeneity and anisotropy of the crustal structure emphasize both the boundary delineation and differential movement between continental blocks. The observed increase in density and intensity of intraplate earthquakes across this boundary supports our conclusion regarding a connection between decoupled crustal deformation and seismicity. Specifically, the foliation of metamorphic rocks during the Mesozoic collision generated orogen-parallel fabrics in the upper-middle crust, while ductile shearing resulted in orogen-perpendicular fast directions at lower crust levels. In contrast, mismatch evolution between continental blocks reactivated by present collision provides stress conditions conducive to earthquake occurrence in southeastern China. Our findings provide valuable insights into mechanisms behind intraplate earthquakes as well as tectonic evolution within continental collision zones.Plain Language SummaryThe mechanism of intraplate earthquakes is a topic of great interest in the study of continental collision and regional tectonic evolution. The controlling factors behind these earthquakes are highly complex and uncertain. In this study, we aim to investigate the role of crustal structures and anisotropy in southeastern China within the context of intraplate seismic activity. This region, which is currently experiencing a collision between the Philippine Plate and the South China Sea Margin, has a multiphase tectonic history that includes remnants of paleo-Pacific subduction. By utilizing the receiver function method, we provide a detailed picture of crustal structures and layered anisotropy in Guangdong and Fujian, two earthquake-prone regions in southeastern China. Our findings reveal that anisotropy in the crust can be stratified into two layers: an upper-middle crust and a lower crust with a slow axis and a fast axis, both perpendicular to the orogeny, respectively. The lateral heterogeneity and anisotropy of the crustal structure highlight the boundary and differential movement between the continental blocks. The enhanced density and intensity of the intraplate earthquakes across this boundary provide evidence for a connection between decoupled crustal deformation and seismicity. Specifically, our observations suggest that (1) the Mesozoic continental collision induced metamorphic foliation in the upper-middle crust, (2) ductile shearing is characteristic of the lower crust, and (3) ongoing collision has led to mismatch evolution between continental blocks, which creates stress conditions favorable for earthquake occurrence.

Passive Source Reverse Time Migration Based on the Spectral Element Method

AbstractIncreasing deployment of dense arrays has facilitated detailed structure imaging for tectonic investigation, hazard assessment and resource exploration. Strong velocity heterogeneity and topographic changes have to be considered during passive source imaging. However, it is quite challenging for ray‐based methods, such as Kirchhoff migration or the widely used teleseismic receiver function, to handle these problems. In this study, we propose a 3‐D passive source reverse time migration strategy based on the spectral element method. It is realized by decomposing the time reversal full elastic wavefield into amplitude‐preserved vector P and S wavefields by solving the corresponding weak‐form solutions, followed by a dot‐product imaging condition to get images for the subsurface structures. It enables us to use regional 3‐D migration velocity models and take topographic variations into account, helping us to locate reflectors at more accurate positions than traditional 1‐D model‐based methods, like teleseismic receiver functions. Two synthetic tests are used to demonstrate the advantages of the proposed method to handle topographic variations and complex velocity heterogeneities. Furthermore, applications to the Laramie array data using both teleseismic P and S waves enable us to identify several south‐dipping structures beneath the Laramie basin in southeast Wyoming, which are interpreted as the Cheyenne Belt suture zone and agree with, and improve upon previous geological interpretations.

THE CRUSTAL THICKNESS AND Vp/Vs RATIO IN MONGOLIA

Crustal thickness beneath Mongolia was determined using the teleseismic receiver function H-k stacking method, which employed waveform data from broadband seismic stations in Mongolia. The findings reveal that the crust of Mongolia thickens from east to west, with the thickness reaching 63 km in west Mongolia and 35.9 km in east Mongolia. In the central part of Mongolia, the average thickness of the crust is approximately 45 km. The characteristics of P-to-S converted phases and negative amplitudes in stacked receiver functions, may be attributed to regional tectonostratigraphic terranes in Mongolia. Subsequently, a general map of the crustal thickness of Mongolia was generated by combining the results of this study with those of other seismological studies and the CRUST1.0 model of the Earth’s crust.

Influence of upwelling mantle magmas on cratonic crust implied from Vp/Vs beneath South America platform

SUMMARY The crust of the South American platform recorded imprints of dynamic processes related with the opening of the Central and South Atlantic but has not been well measured. Crustal structure can be retrieved from teleseismic receiver functions using H–κ stacking, but nearly parallel stripes of high stacking values existing in stacking images for seismic stations in sedimentary area cause difficulties in identifying solutions. We show that some seemingly spurious stripes that do not point to any layer solution are helpful in the identification of the solution position. With the aid of the auxiliary stripes, we retrieved thicknesses and Vp/Vs of sedimentary and crystalline crust for 65 permanent stations of the Brazilian Seismographic Network and six new portable seismic stations in Brazil and Uruguay. The resulted sedimentary thickness and Vp/Vs exhibit a good correlation with the Phanerozoic sediments in the South American basins. The crust of Paraná–Etendeka Large Igneous Province (LIP) had been expected to be more mafic since it had ever been penetrated by mantle magma in the Cretaceous related to the south Atlantic opening. However, we found very low Vp/Vs (1.67) in the crystalline crust beneath the LIP, implying a more felsic crust and that no significant mafic intruding/underplating has occurred in the region. The more felsic crust may be formed in a special evolution early than the magmatic event, or during the magmatic event by releasing crustal volatiles. The resulted sedimentary thickness and Vp/Vs ratios exhibit a good correlation with the Phanerozoic sediments in the South American basins, which implies that Triassic–Jurassic and Cretaceous magmatism did not cause significant metamorphism in sediments formed before the magmatic events.

A translithospheric magmatic system revealed beneath Changbaishan volcano

Abstract Changbaishan volcano (CBV), located on the border between China and North Korea, has undergone violent eruptions since the Oligocene, making it one of the most captivating and explosive volcanoes on Earth. However, the lack of precise characterization regarding the magmatic system makes it difficult to decipher its eruption risk and mechanism, despite its significant size and past devastating effects. In this study, we employed newly developed teleseismic receiver function techniques, including Ps and Sp waves, to construct a lithospheric structure model beneath the CBV region. The results show a thick crust (~37 km) and a weak, thin lithosphere under the CBV, with a low-amplitude width of 80 km at the Moho depth and 200 km at the depth of the lithosphere-asthenosphere boundary. These features depict the seismological response of a translithospheric magmatic system beneath the CBV, where hot upwelling material rises through the lithospheric mantle, underplates at the base of the crust, and forms the magma chamber(s) at shallow depth. Such magmatic system features can be taken as a unifying paradigm for large volcanic regions worldwide.

Injection-induced basement seismicity beneath the Raton Basin: constraints from refined fault architectures and basin structure.

Intraplate earthquakes induced by anthropogenic fluid injection present unexpected seismic risk to previously quiescent or low seismicity-rate regions. Despite many studies of induced seismicity, there are relatively few with detailed openly accessible constraints on the interaction between seismic sources and subsurface structures. In this study of the Raton Basin, we refine source observations from a dense nodal array and constrain basin structure using teleseismic receiver functions. The cross-correlation-based relocated hypocentres and a new set of focal mechanisms light up active fault segments and show clear spatiotemporal patterns. The geometric complexity of reactivated fault clusters appears greatest near higher rate injection wells. Simpler normal fault structure is found farther from injection wells and near abrupt structural transitions suggested by receiver functions. While less induced seismicity in the crystalline basement is expected when injection is >1 km from the top of the basement (like Raton), our receiver function analysis identified a basin thickness ~3 km beneath the nodal array and lateral variations in sedimentary structures. Our results explain potential fluid connectivity between the injection depths focused at ~1-1.5 km below the surface and basement fault activity that begins at ~3 km and reaches peak activity at ~4-8 km depths. This article is part of the theme issue 'Induced seismicity in coupled subsurface systems'.

A Seismic Nodal Deployment to Understand Magmatic Structure in the Vicinity of the Pāhala Earthquake Swarm

Abstract In Summer–Fall 2022, 80 three-component SmartSolo IGU-BD3C-5 nodal seismometers were deployed surrounding the Pāhala seismic swarm on the Island of Hawai‘i, with the goal of improving seismicity catalogs and seismic velocity images of the crust and upper mantle in this region. The Pāhala swarm, located south of Mauna Loa and Kīlauea, has been the site of a multiyear sustained swarm of seismicity at the depths of ∼25–40 km, with order of magnitude increases in rate in 2015, and then again in 2019. This seismicity is possibly related to the input of magma from the mantle plume below, which may then be subsequently transported to volcanic edifices. However, these processes remain enigmatic, in part due to a lack of precise earthquake locations and seismic velocity models in this region. Here, we provide an overview of the deployment, an assessment of the quality of the collected data, and discuss the viability of the data set for local earthquake relocation, tomography, and teleseismic receiver functions. Through comparisons with proximal permanent broadband and short-period instruments, we find that the nodes produce high-quality data, particularly at periods shorter than 5 s, although we find, document, and correct discrepancies with the gain and polarities of the instruments. We successfully record signals from teleseismic earthquakes, even at periods longer than 5 s (the corner of the flat response of the nodes). We also record local earthquakes, including details related to source characteristics. This indicates that the data are likely to prove useful for investigations using both local and teleseismic earthquake signals to better understand the connections between the deep and shallow magmatic systems of Hawai‘i. Although this deployment provides a snapshot in time, its success may provide a useful benchmark for future studies as the volcanic systems of Hawai‘i continue to evolve in the future.

A Displaced Lower Mantle Source of the Hainan Plume in South China Revealed by Receiver Function Imaging of the CEArray

AbstractWe analyzed 49,592 teleseismic receiver functions (RFs) recorded by 278 CEArray stations to image the mantle transition zone (MTZ) beneath the South China Block to understand the origins of deep velocity anomalies and their potential links to subduction and intraplate volcanism. We employed a fast‐marching method and a high‐resolution 3‐D velocity model (FWEA18) derived from full waveform inversion in computing P‐to‐S conversion times to better image the 410‐ and 660‐km discontinuities. Our results indicate that the common‐conversion‐point stacking of RFs using 3‐D conversion times yielded better migration images of the two discontinuities. The images revealed a slightly depressed 410‐km with a few small uplifted patches, and showed that the 660‐km beneath the western Yangtze Craton is depressed by 10–25 km, which is likely caused by the stagnant Paleo‐Pacific slab. The 660‐km beneath the southern Cathaysia Block has a 5–15 km high plateau with a topographic low at its central part. The lateral dimension of the topographic low is ∼150 km and is located beneath the central Pearl River Mount Basin near Hong Kong. We speculate that the topographic low occurs within the Hainan plume with a temperature excess of ∼300–400 K and is caused by the garnet phase transition. The displaced deep plume enters the MTZ and spreads nearly horizontally at the base. The plume evolves into two channels with a minor one toward the northeast and a major one toward the southwest, which keep moving upward to the 410‐km. The southwest channel is likely the source that feeds the Hainan volcanoes.

Crustal thickness variation of Kalimantan and Sulawesi Region from Teleseismic receiver function

Crustal thickness variation of Kalimantan and Sulawesi Region from Teleseismic receiver function

A lower crust shear zone facilitates delamination and continental subduction under the Apennines

Physical properties and structure of the lithosphere are the first step to constrain the evolution of mountain belts. Here we show detailed shear wave velocity profiles of the lithosphere in the Apennines that clarify a controversial aspect of continental subduction: the intricate mechanism of crust delamination from the downgoing plate. From the analysis of complete and dense teleseismic Receiver Function data set, we find that the delamination of the continental lithosphere is favored by the development of a low seismic shear wave velocity zone in the middle-lower crust. We observe a double Moho below the external portions of the present mountain range, suggesting the progressive formation of the shallow interface. The delamination edge is located in the forearc, far eastward than expected, implying that the re-equilibration of the thermal unbalance, generated by the mantle substitution, may last 10-7 Myr.

Crustal Imaging with Noisy Teleseismic Receiver Functions Using Sparse Radon Transforms

ABSTRACT The receiver function (RF) is a widely used crustal imaging technique. In principle, it assumes relatively noise-free traces that can be used to target receiver-side structures following source deconvolution. In practice, however, mode conversions and reflections may be severely degraded by noisy conditions, hampering robust estimation of crustal parameters. In this study, we use a sparsity-promoting Radon transform to decompose the observed RF traces into their wavefield contributions, that is, direct conversions, multiples, and incoherent noise. By applying a crustal mask on the Radon-transformed RF, we obtain noise-free RF traces with only Moho conversions and reflections. We demonstrate, using a synthetic experiment and a real-data example from the Sierra Nevada, that our approach can effectively denoise the RFs and extract the underlying Moho signals. This greatly improves the robustness of crustal structure recovery as exemplified by subsequent H−κ stacking. We further demonstrate, using a station sitting on loose sediments in the Upper Mississippi embayment, that a combination of our approach and frequency-domain filtering can significantly improve crustal imaging in reverberant settings. In the presence of complex crustal structures, for example, dipping Moho, intracrustal layers, and crustal anisotropy, we recommend caution when applying our proposed approach due to the difficulty of interpreting a possibly more complicated Radon image. We expect that our technique will enable high-resolution crustal imaging and inspire more applications of Radon transforms in seismic signal processing.

Characterization of anisotropy in basin-scale subsurface using teleseismic receiver function analysis

Teleseismic receiver function (RF) analysis offers a passive-source analogue to detect impedance boundaries using the converted body waves generated by earthquakes. While the technique traditionally has targeted deep earth structures such as the Moho and transition zones, there is growing interest in assessing its applicability in basin-scale seismic characterization, ultimately aimed for onshore commercial integration as a cost-effective complement to existing active-source seismic surveys. Specifically, the conventional broadband seismometers used in global observational seismic experiments are not only logistically simple in terms of data acquisition, but they also record ground motions in three mutually orthogonal time series, enabling effective detection of shear waves and directional variations of observed signals. Here, we perform teleseismic RF analysis to detect shear-wave anisotropy and related symmetry axes orientations in a basin setting, using open-source seismic data recorded at 55 closely spaced seismic stations in the LaBarge Passive Seismic Experiment deployed in Wyoming between November 2008 and June 2009. We find that the strengths and geometry of the observed anisotropy are variable along the array. Significantly, not only can anisotropy effectively delineate subsurface fault interfaces, it can also substantiate and reveal additional interpretable signals that are otherwise disregarded. The estimated fast axes orientations compare favorably with the complex fracture systems documented in the region. Finally, we show that P-to-S amplitude variations with P incidence are systematic and modelable using existing computational tools, offering an opportunity to develop an analysis technique similar to amplitude variation with offset with the products of RF analysis.

Crustal properties investigation in the Meratus Complex Zone, Southeast Borneo, using H-κ stacking of receiver function analysis

We study crustal properties using teleseismic receiver functions extracted from two broadband seismic stations around the Meratus Complex region, Southeast Borneo. The two seismic stations, BBKI and KBKI, are distributed in the Meratus Mountains and Pulau Laut. We then employ the H-κ stacking technique to measure the crustal thickness (H) and the crustal average Vp/Vs ratio (κ). The method works quite well for receiver functions recorded on seismic station KBKI. The thickness of the crust and Vp/Vs ratio underneath this station are around 38 km and 1.64, respectively. On the contrary, the standard H-κ method fails in estimating the crustal properties beneath station BBKI due to significant sedimentary contamination on the receiver functions. Thus, we apply the sequential H-κ stacking analysis by removing the sedimentary effect and increasing the Moho Ps amplitude. The result shows the thickness of the crust is around 29 km, although the Vp/Vs is not well determined. The research findings from both stations are comparable to those derived from the previous receiver function study using inversion modelling. The crustal properties around this area may reflect the continental material part of the Sunda continental block.

Moho Depth Estimation and The Presence of Subducting Slab in West Sumatra, Jambi, and Riau Islands Regions Using Teleseismic Receiver Function Method: A Preliminary Result

The teleseismic receiver function method was used to investigate the crustal structure in West Sumatra, Jambi, and Riau Islands regions. These regions have a high level of seismicity because it is in an active tectonic zone where the oblique subduction between the Indo-Australian and Eurasian plates occurs. This study aims to determine the depth of Moho discontinuity and the presence of subducting slab beneath 6 BMKG seismic stations that form a line perpendicular to the trench and cover 3 different geological zones. PPSI and RASM are in the forearc zone, SKSM is in the volcanic arc zone, and MBBI, JMBI, and DSRI are in the back-arc zone. This study used teleseismic earthquake record data with epicentral distance between 30°-90° from the receiver and magnitude 6 or more. The iterative time-domain deconvolution and receiver function migration techniques were applied to estimate the depth of Moho discontinuity and the presence of subducting slab. The depth of Moho discontinuity in West Sumatra, Jambi, and Riau Islands regions ranges from 21-39 km and generally deeper in the volcanic arc zone possibly due to the isostatic effect. Moho is at a depth of 21-29 km in the forearc zone, 36-39 km in the volcanic arc zone, and 30-36 km in the back-arc zone. Then the subducting slab was observed at a depth of 20 km under PPSI station to 200 km under MBBI station.

A Preliminary Results: Study of Crustal Structure in Cimandiri Fault Based on Teleseismic Receiver Function Analysis

West Java is one of the regions in Indonesia that has a complex geological structure caused by plate subduction in the south of Java Island. This resulted the formation structure of the main faults in West Java, including the Cimandiri Fault. The high level of seismicity on the Cimandiri Fault is interesting to be further researched to understand the structure of the Earth crustal, such as variations depth of the Moho layer and the determination of local velocity models of P and S waves. This study used receiver function analysis from teleseismic earthquake data with an epicentral distance of 30°-90° from the station location and magnitude ≥ 6. This study used 3 seismic stations belonging to the Meteorological, Climatological, and Geophysical Agency (BMKG) close to the Cimandiri Fault which are Sukabumi station (SKJI), Warung Kiara (WSJM), and Sukaraja (SSJM). The depth of the Moho layer at SKJI station is 34-36 km, WSJM station is 28-30 km, and SSJM station is 29-31 km which is located close to the segment from the Cimandiri Fault. A negative phase or decrease in speed was observed at the SKJI and WSJM stations which is estimated to come from the Cimandiri Fault structure which is at a depth of 17-19 km. Furthermore, the P wave velocity model at the three stations is at an interval of 7.6-8.6 km/s and the S wave velocity model is at an interval of 4.2-4.8 km/s.

Crustal Thickness and Vp/Vs Ratio Determination beneath Tangkuban Parahu Volcano, West Java, Indonesia using Receiver Function Analysis

Tangkuban Parahu, an active stratovolcano located approximately 15 km north of Bandung City in Indonesia, is the focal point of this study. The research aims to examine seismic structures within the crust beneath the Tangkuban Parahu Complex. To accomplish this, a teleseismic receiver function method was employed, utilizing broadband seismic stations along a seismological profile within the Tangkuban Parahu Volcano area. Waveform data were collected during the period from October 24, 2021, to March 31, 2022, for analysis. Various seismic wave phases, including P, Ps, PpPs, and PsPs, were identified in the recorded waveforms, within a time window extending from -2 seconds before the onset of the P-wave to 20 seconds after its onset. The analysis and modelling were conducted using the H-κ stacking receiver function technique. The findings of the H-κ stacking receiver function revealed notable variations in crustal structures, encompassing crustal thickness and the average crustal composition, as indicated by the Vp/Vs ratio, within the Tangkuban Parahu Volcano area. Preliminary results suggest that the crust beneath Tangkuban Parahu, estimated based on data from single seismic station, possesses a thickness of approximately 32.3 km, with a corresponding Vp/Vs ratio of around 1.746.

Orogen-parallel discontinuity of the Apennines subduction zone in Southern Italy as seen from mantle wedge seismic structure

We investigate the seismic structure of the mantle wedge of the Apennines subduction zone (Central Mediterranean) using teleseismic receiver function (RF). We inverted RF for both isotropic and anisotropic properties of the mantle wedge, from below the overriding Moho to the “plate boundary”, i.e. the interface that separate the slab from the mantle wedge. Given the distribution of the seismic network, we are able to map out the change in the elastic properties at the transition between southern apennines and the Calabrian arc, given by the change in the subduction style (i.e from the subduction of continental materials to oceanic plate). We found that the anisotropy in the mantle wedge is similar between all seismic stations, generally highly anisotropic (> 10%), with a direction of the symmetry axis that rotates clockwise from North to South, following the Calabrian arc geometry and likely indicating the mantle flow driven by the slab retreat. The elastic properties of the subducted crust are more heterogeneous. To the North, the subducted crust shows a highly anisotropic (> 10%) behavior, and it occurs at larger depth (around 70 km depth), where to the South anisotropy is less intense (around 7%) and the subducted crust is shallower (around 60 km depth). These results point out a change in the subduction style that can be given by either a change in the metamorphic phase (more evolved blueschist facies stage to the North, initial greenschist facies stage to the South) or a different origin for the subducted materials (continental to the North and oceanic to the South). The differences in the anisotropic behavior of the subducted crust are reflected in the topography of the plate boundary, which becomes shallower from North to South, suggesting the existence of either a step in the slab topography or a more gentle ramp.