Characteristics of azimuthal anisotropy in SE Tibetan plateau and its relationship with the background of block structure

The southeastern Tibetan Plateau, which includes the Tibetan Plateau, Yangtze Block, and Cathaysia Block, is geotectonically situated in the compound part of the Tethys-Himalayan tectonic domain and the Peninsular Pacific tectonic domain. It is one of the critical regions for studying the scientific problems of plateau material lateral escape, lithosphere deformation, geotectonic properties of blocks, and deep dynamics. In this study, we use ambient noise data recorded by 401 broadband stations to obtain high-resolution short-period (T = 4–32 s) Rayleigh wave azimuthally anisotropic phase velocity maps. These could provide fresh clues for an in-depth understanding of the crust-mantle velocity structure, deformation mechanism, and geotectonic evolution in the southeastern Tibetan Plateau. Within the Simao block, the strikes of the faults and the orientations of the principal compressive stress of the stress field both generally coincide with the fast-wave polarization direction (FPD). The FPD near the Lancangjiang fault zone in the west is in the NE-SW direction, near the Wuliangshan fault zone in the center is near the NS direction, and near the Red River fault zone in the east is the NW-SE direction. We estimate that the compressive stress in the southwest direction of the Tibetan Plateau material has a controlling effect on the crustal deformation of the Simao block, which is likewise blocked by the Lincang granite belt, resulting in strong tectonic deformation. The FPD of the crust in the middle Red River fault zone is NS direction, significantly different from the fault strike. Combining with the seismic activity and GPS results, the depth of 8 km below the surface of the middle Red River fault is completely locked, and we conclude that the anisotropy of the upper crust of the middle part of the Red River fault zone is related to the action of the regional tectonic stress field. Taking into account geochemical and thermochemical results, we speculate that the complex tectonic stress at the junction of the blocks leads to prominent regional characteristics of the FPDs of azimuthal anisotropy in the crust, suggesting that the Shizong-Mile fault zone may be the western boundary between the Yangtze block and the Cathaysia block.


Introduction
The southeastern Tibetan Plateau (TP) is the Frontier zone and the deep material escape zone for the NE collision and extrusion of the Indian plate and the Eurasian continent ( Figure 1A). There are the Tengchong block (TCB), Baoshan block (BSB), Simao block (SMB), Chuanxibei block (CXB), Dianzhong block (DZB) and other tectonic units in the area. Since the Cenozoic, it has experienced intense tectonic deformation and frequent seismic activities (Su and Qin, 2001), making it an ideal place to study the interaction among blocks and deep tectonic deformation.
Seismic anisotropy is one of the effective methods to study the deformation of the crust and upper mantle. The upper crustal anisotropy is related to the shape-preferred orientation of microstructures such as lamellae, joints, and fractures in the rock. Furthermore, the anisotropy of the mid-to-lower crust and upper mantle is related to the lattice-preferred orientation of anisotropic minerals such as amphibole, mica, and peridotite. We study seismic anisotropy to understand the relationship between medium and anisotropy, anisotropy and strain, and the relationship between strain and crustal movement. These help us infer the tectonic activities that might have operated in the past and those working at present. Some scholars have studied the SKS-wave splitting measurement in the southeastern TP and believe that, unlike the vertically coherent deformation of the crust and mantle inside the TP, the crust-mantle deformation in the southeastern TP is decoupled (Silver and Chan, 1991;Flesch et al., 2005;Lev et al., 2006;Chang et al., 2015;Huang et al., 2015). However, Wang et al. (2008) believed that the vertical coherent deformation model could also explain SKS-wave splitting measurements outside the TP. The anisotropy obtained by SKS-wave splitting measurement is the cumulative effect of the anisotropy on the path from the coremantle boundary to the surface, which is difficult to determine the depth of the anisotropic layer. Chen et al. (2013) and Cai et al. (2016) speculated that the crust and upper mantle deformation in the southern study area decoupled. However, Pms anisotropy can only reveal the average anisotropy characteristic of the entire crust. Due to the dispersive properties of surface waves, surface waves of different periods are sensitive to structure information at different depths. Therefore, compared with SKS wave splitting and Pms anisotropy, surface wave azimuthal anisotropy can provide better constraints on the depth variation of anisotropy.
By using the seismic surface wave and ambient noise data, some researchers illustrate the azimuthal anisotropy of the TP and its vicinities Lu et al., 2014;Wang et al., 2015;Liu et al., 2019;Liang et al., 2020;Zhu et al., 2021). These studies reveal the anisotropy on a large scale, i.e., mid-to-lower crust and upper mantle depth range, which provides essential evidence for the internal   (Zhao et al., 2015) of the SE Tibetan Plateau relative to the stable Eurasia. The rose-red band indicates the Lincang (LC) granite belt. The bottom-right inset highlights Main block boundary lines. (B) Distributions of the stations used in our tomography study. Red, green, blue and yellow triangles show the location of the ChinArray seismic network (ChinArray, 2006), Xiaojiang Array (Wu et al., 2013), Chuanxi Array  and permanent seismic stations in the southeastern Tibetan Plateau (Zheng et al., 2010) structure and deformation dynamic mechanism of the crust and upper mantle in the southeastern TP. However, there is still a debate about the interactions between the sub-blocks and their natural properties. For example, the SMB, located between the Red River Fault zone (RRF) and the Lancangjiang Fault zone (LCJF), is spread from NW to SE (Figure 1). It is adjacent to the Dianzhong block in the north and the Baoshan block in the west. Significant tectonic deformation has occurred since Cenozoic. Geological studies propose a variety of different structural deformation modes. It mainly includes the pushing out effect during lateral rotational extrusion of crustal blocks in the southeastern TP (Tanaka et al., 2008;Kondo et al., 2012;Tong et al., 2013); the extrusion effect of the thrust fault in Lincang (LC) granite foundation from west to east , and the blocking effect of Lincang granite foundation on the lateral clockwise rotational extrusion movement of SMB (Xu et al., 2017). The eastern boundary between the YZB and the CYB can be made with certainty (Zhang et al., 2005;Zhang et al., 2015;Guo et al., 2019;Ma et al., 2022), but the western border is more controversial (Zhu et al., 2019). The edge of the west is difficult to discern because of the lack of outcrops of the combined geological landmarks (ophiolites, island arc igneous rocks), the thick cover in the later period, and the multi-stage structure and magmatic intrusion. The one point of views is the Shizong-Mile fault (SZMLF, Figure 1B) (Dong et al., 2002;Guo et al., 2009).
In our study, we collect the continuous waveform data recording by dense seismic arrays of the permanent China National-backbone stations and mobile broadband stations (Figure1B). We adopt ambient noise tomography to acquire high-resolution Rayleigh wave phase velocity azimuthal anisotropy in this area. It provides new evidence for the tectonic evolution and dynamic process of the crust and upper mantle in the southeastern TP.

Data acquisition and processing
In this study, we collected the continuous waveform data recorded by three temporal seismic arrays: 195 stations of the ChinArray (X1) between October 2011 and October 2012, 114 stations of the Chuanxi array (CX) between October 2006 and July 2009 , and 50 stations of the Xiaojiang array (XJ) between January 2009 and December 2010 (Wu et al., 2013). In addition, we also use data recorded by 42 permanent stations of the Sichuan and Yunnan regional seismic networks for 1 year in 2009. There are 401 seismic stations in total in our study region. The average interstation space is about 30 km. The dense station distribution enhances the azimuthal coverage and path density, especially in DZB, compared with previous studies (e.g., Lu et al., 2014;Wang et al., 2015). Figure 1B shows the station distribution. Most stations are equipped with Guralp CMG-3ESPS seismometer (60s-50 Hz) and Reftek 130 digitizer.
The data processing technique for ambient noise tomography is currently in a relatively advanced stage. It mainly includes single station data preprocessing, empirical Green's function extraction, phase velocity dispersion measurement, phase velocity tomography, and surface wave azimuthal anisotropy inversion. Empirical Green's functions (EGFs) are derived from the cross-correlation functions between each station pair using the vertical component of continuous waveform following the method of Bensen et al. (2007) and Fang et al. (2009). Phase velocity dispersions from EGFs are measured using the far-field representation of the surface-wave Green's function and image transformation technique following Yao et al. (2006). Figure 2 shows example of Z-Z cross-correlation functions (CCFs) between stations randomly selected from XJ array. The CCFs are filtered at 5-40 s period band. Figure 3 shows an example of the dispersion measurement of the phase velocity.

Surface wave tomography inversion
In the case of weak anisotropic medium and ignoring 4φ terms, Rayleigh wave phase velocity c(ω, M, φ) at an arbitrary point M for each angular frequency ω and azimuth φ can be expressed as (Smith and Dahlen, 1973): Where c 0 (ω) is the reference phase velocity, α 0 and α i (i 1, 2) are the isotropic phase velocity perturbation and the azimuthal anisotropy coefficients, respectively. During the inversion, these parameters are determined using the continuous regionalization method of Montagner (1986) and the generalized inversion algorithm of Tarantola and Valette (1982) for each period. The

FIGURE 2
Example of Z-Z CCFs between stations randomly selected from XJ array. The CCFs are filtered at 5-40 s period band.
The inversion for α i (i 0, 1, 2) is controlled by three parameters: the standard error of phase velocity measurements, δd, the a priori parameter error, δp (which constrains the anomaly amplitude), and the correlation length Lc (which constrains the smoothness of the resulting model). In our inversion, following Yao et al. (2010), we set δd to 2% for all measurements. For a given α 0 , δp is set to be twice that of the standard deviation (in percent) of all observed phase velocities at each period with a minimum value of 0.15 km/s. For α 1 and α 2 , δp is set to be 1.5% of the average phase velocity at each period. In this study, we focus on short and intermediate period surface wave. It may result in artificial velocity anomaly if the correlation length Lc is too small. So we set the correlation length Liso=max (30 km, C 0 *T/2) for the isotropic term, where C 0 is the average phase velocity at a certain period, T. The correlation length for the azimuthally anisotropic parameters is set to be 2*L iso at the corresponding period.

Resolution tests
The knowledge of the resolution is vital in the interpretation of the tomography results. Artificial values may be introduced in inversion due to insufficient data and poor ray path coverage. Following Simons (2002) and Yao et al. (2010), we performed checkerboard tests to assess the resolution of our model. We calculated the synthetic data from input models with a 1°×1°a nisotropic pattern along with the fast azimuthal axis-oriented NE and NW, with 5% isotropic anomalies at each period ( Figure 4A). We added the random error with a maximum of 1% to the synthetic data. We performed the inversion with the same parameters as for the actual data. The checkerboard tests indicate that for the azimuthal anisotropy parameters, the 1°×1°patterns are well recovered throughout our study area ( Figures 4B-D). The azimuthal anisotropy is not well recovered in the margin due to relatively poor azimuthal path coverage. Generally, checkerboard tests provide qualitative information about spatial resolution. However, we must remember that the checkerboard test differs from the inversion of actual data due to its inherent limitations, as pointed out by previous studies (e.g., Simons, 2002). For example, the limits of ray approximation, the horizontal smoothing imposed in the inversion, and the trade-off between lateral heterogeneities and anisotropy (Pandey et al., 2015). Frontiers in Earth Science frontiersin.org Tarantola and Valette (1982) pointed out that calculating the posterior covariance operator of the model parameters after inversion can estimate the posterior error and resolution of the inversion parameters. If the posterior error is significantly smaller than the prior error, the inversion results are reliable, and the resolution is high. Conversely, when the posterior error is close to the prior error, it indicates that the data cannot resolve the parameter. Figure 5 shows the posterior error distribution of the anisotropic intensity after inversion. Its amplitude is much smaller than the prior error of the anisotropic intensity and the anisotropic intensity in the study area. Taking 10s as an example, in the central, the posterior errors of the anisotropic intensity Ac are all less than 0.3%. Therefore, the inversion results are reliable for regions with anisotropy strength greater than 2%.

Results
The Rayleigh surface wave phase velocity and azimuthal anisotropy are obtained by inversion according to the mixed path dispersion curve. The different periods reflect the structural differences in different depth ranges. Compared with the S-wave velocity inverted by a single node, the reliability of the phase velocity at different periods is higher. We select six periods for discussion according to the phase velocity characteristics at each period. Figures 6A-F show isotropic phase velocity and azimuthal anisotropy maps at six periods. Moreover, the background is the absolute value of the isotropic phase velocity. The length and azimuth of the short black bars indicate the magnitude of the azimuthal anisotropy and the FPD, respectively. At short periods (4-10 s), the phase velocity mainly reflects the characteristics of the upper crustal structure. The phase velocity distribution is consistent with the surface geological structure. The SMB exhibits low velocity, associated with thick sedimentary layers in this area, consistent with the S-wave velocity results of Liu et al. (2019). There are high-velocity anomalies in the vicinity of Panzhihua (PZH) and the east of the Shizong-Mile fault (SZMLF) zone. An N-S-orientated high-velocity belt along the Lvzhijiang (LZJF) and Yimen faults (YMF) through the Red River Fault (RRF). The FPD of surface wave azimuthal anisotropy is nearly N-S direction in the Dianzhong block and SMB, which is consistent with the strike of the Xiaojiang fault (XJF). In the high-velocity body of PZH, the magnitude of azimuthal anisotropy is weak. The FPD to the east side of the XJF zone is NE. The FPD to the east of the SZMLF is NW-SE direction as a whole, which is significantly different from its adjacent areas.
At intermediate periods (14-20 s), the phase velocity mainly reflects the mid-crustal structure, and its sensitivity range is about 10-25 km. The area near PZH and east of the SZMLF still exhibits high-velocity anomalies. Still, the area of high-velocity anomalies near the PZH area has begun to decrease. After 20 s period, highvelocity anomalies began to appear on the west side of the SMB. The azimuthal anisotropic characteristics at intermediate periods are similar to the short periods. The FPDs to the northeastern SZMLF changed from NW-SE to N-S direction. Furthermore, the FPDs to the south of the eastern RRF zone are nearly E-W direction, and the magnitude of azimuthal anisotropy is relatively strong.
At long periods (>20 s), the phase velocity is sensitive to the structures in the mid-to-lower crust. For example, the SZMLF is the boundary between the YZB and CYB. The high-velocity anomaly on the east side is parallel to the strike of the SZMLF. As the period increases, the area of high-velocity anomalies near the PZH area continues to decrease. After 30 s period, this highvelocity anomaly almost disappeared. The azimuthal anisotropic characteristics seem to be those at 12-20 s periods. The FPDs in the DZB and SMB are still in a close N-S direction. After 24 s, the FPDs in the northern CYB changed from NW-SE and nearly NS to NE-SW. The FPDs to the south of the eastern RRF zone are still near the E-W direction.

Crustal deformation mechanism of simao block and its surrounding areas
The SMB has frequent seismic activities, dominantly about magnitude 6.0, showing the characteristics of stronger and less large earthquakes. The present crustal deformation and intense earthquake activities in this area are closely related to the collision of the Indian plate and the Eurasian plate and the lateral extrusion of the TP. Therefore, understanding the velocity structure and anisotropy characteristics will help us further understand the dynamic relationship between the present crustal deformation process and the material extrusion in the southeastern TP.
The velocity anomaly in the SMB is heterogeneous. At 14-32 s periods, it is mainly sensitive to the S-wave velocity structure in the depth range of 20-45 km (Figure 7). The area near the Wuliangshan fault zone (WLSF) continues to exhibit low-velocity anomalies, and the other areas show relatively high-velocity anomalies relative to short periods. Regarding azimuthal anisotropy ( Figures 8A, B), the FPDs in the SMB show noticeable regional variance. In the central, the FPDs are generally similar to the fault strikes and the orientation of the principal compressive stress (Sheng et al., 2022). Among them, the FPD near the LCJF in the west is NE-SW direction, the FPD near the WLSF zone in the middle is near the N-S direction, and the FPS near the RRF zone in the east is NW-SE direction. These phenomena are consistent with the phase velocity azimuthal anisotropy of Wang et al. (2015), the group velocity azimuthal anisotropy of Lu et al. (2014), and Pms anisotropy Sun et al., 2015;Cai et al., 2016) (Figure 9B).
Previous studies have shown that the RRF experienced a leftlateral strike-slip movement at 32-17 Ma and transformed into a right-lateral strike-slip movement at about 5 Ma (Gilley et al., 2003). In contrast, the XJF zone started to experience an initial left-lateral strike-slip movement from the Middle Miocene to the Early Pliocene (about 17-5 Ma) (Roger et al., 1995). At this time, the Chuandian diamond block (including the CXB and DZB) began to rotate and extrude clockwise along the XJF. During this process, the Chuandian

Frontiers in Earth Science
frontiersin.org diamond block was bound to have a robust southward extrusion on the SMB. Since the Oligocene (37-21 Ma), the Lincang granite in the southern SMB has experienced a gradual uplift from south to north (Shi et al., 2006). Furthermore, with the obstruction of the Lincang granite belt, the crustal structure of the central SMB was deformed strongly, forming a northwardly protruding honeycomb-shaped structure. Therefore, when subjected to compressive stress in the southwest direction of the Chuandian diamond block, the crust of the eastern SMB underwent shear deformation, forming crustal anisotropy parallel to the strike of the RRF. Meanwhile, the crust of the western SMB was blocked by the Lincang granite belt and strongly deformed, forming crustal anisotropy parallel to the strikes of the Lincang granite belt and the LCJF. Likewise, the central is in the transition zone of differential deformation on both sides, forming near N-S trending crustal anisotropy, which is almost consistent with the extensive development of folds and thrust structures in the Mesozoic-Cenozoic strata in the block (Guan et al., 2006).
The SMB has undergone 50°-70°clockwise rotational deformation under the combined action of the Chuandian diamond block since the Pliocene (Sato et al., 2007). The GPS velocity field shows that, relative to the South China block, the NNE pushing caused by the collision of the Indian plate and the Eurasian plate (Wang and Shen, 2020) and the gravitational potential energy generated by the plateau uplift resulted in the eastward extrusion of the TP . After encountering the obstruction from the stable South China block, the TP material moved southeast and southward, causing the Chuandian diamond block to rotate clockwise around the eastern Himalayan tectonic Variation of absolute Rayleigh phase velocities (color image) and azimuthal anisotropy (short black bars) in the crust and uppermost mantle beneath southeastern Tibetan Plateau at different periods (A-F). The color bars show the color scale of phase velocities (km/s). The short brown bars are the P-axis directions of the principal stress based on the focal mechanism solution (Ekström et al., 2012). The abbreviations of geological units are defined in the caption of Figure 1.

Frontiers in Earth Science
frontiersin.org knot. Therefore, we infer that the compressive stress in the southwest direction of the Chuandian diamond block has a controlling effect on the crustal deformation of the SMB, and the west of the block was blocked by the Lincang granite belt and experienced intense tectonic deformation ( Figure 10). In the middle RRF and its vicinity, the anisotropic FPD at 4-24 s periods (corresponding to the crustal depth range, Figure 7) is near the N-S direction (Figure 5), showing a large angle with the fault strike. The azimuthal anisotropy of the surface wave group velocity, Pms anisotropy, and S-wave azimuthal anisotropy ( Figure 9B) also show almost the same characteristics (Lu et al., 2014;Cai et al., 2016;Liu et al., 2019). The RRF zone, as a critical plate boundary cutting the Moho surface , did not cause the FPD of the anisotropy in the crust to be parallel to the strike of the fault. The tectonic stress field obtained from the focal mechanism (Xu, 2001;Wu et al., 2004;Xu et al., 2016) and the GPS velocity field (Jin et al., 2019) show that the maximum principal stress in the vicinity of the middle RRF zone is in the near N-S direction, which is consistent with the FPD of azimuthal anisotropy ( Figure 9A). The anisotropy of the upper crust is generally considered to be related to the shape preferred orientation of microstructures (Rabbel and Mooney, 1996;Crampin and Peacock, 2008). Secondly, fault zones may also affect it, but the influence is limited (Gao et al., 2011). The GPS velocity field ( Figure 1A) and seismic activity suggest that the middle RRF zone is completely locked from the surface to 8 km depth (Zhao et al., 2015;Wang et al., 2022). We speculate that the anisotropic material in the upper crust of the middle RRF zone may not have a dominant arrangement along the fault strike but is related to the action of the regional tectonic stress field (Figure 10).

Constraints on the western boundary between Yangtze Block and Cathaysia Block
At 4-32 s periods, the Rayleigh surface wave phase velocities in the eastern and western sides of the SZMLF have noticeable velocity differences ( Figures 8A, B, D), which correspond well with the S-wave and P-wave velocity variations (Wu et al., 2013;Yang et al., 2020). Unlike the velocities on the west side of the fault zone with the periods changing, the east side of the fault zone continues to show relatively high phase velocity anomalies. The apparent velocity differences may reflect the inhomogeneity of the crustal structures and lithologies on both sides. The receiver function showed the crustal thickness on both sides of the SZMLF differs up to 10 km   (Figure 8D). The geochemical study of the volcanic rocks in the northern SZMLF zone shows that it is mainly alkaline basalt, characterized by low TiO 2 and high Al 2 O 3 , which is different from the Emeishan continental overflow basalt with high TiO 2 and low Al 2 O 3 characteristics (Dong et al., 2002). The Bouguer gravity anomaly and the equilibrium gravity anomaly along both sides of the SZMLF zone are high in the north and low in the south, which shows that the SZMLF is the block boundary separating the YZB and CYB with different crustal structures. In addition, the thermochemical states of the upper mantle in South China based on multi-observation probability inversion (Yang et al., 2021) show that in the lithospheric mantle, the average density in this region also takes the SZMLF zone as the western boundary to separate the YZB and CYB. Coupled with the complex tectonic stress at the junction of the blocks leading to prominent regional characteristics of the FPDs, we infer that the SZMLF zone may be the western boundary between the YZB and CYB in the crust (Figure 10).
The SCB located at the northwest of the YZB. In the northeastern study area, the crustal anisotropy of the southern SCB at 4-20 s periods is dominated by NE-SW and NEE-SWW trends ( Figure 8C), which is consistent with the NE-SW trend of the surface eastern Sichuan fold belt Xiong et al., 2016). The shear wave splitting of local earthquakes (Tai et al., 2015), Pms anisotropy (Sun et al., 2012), and S-wave azimuthal anisotropy  also show the same crustal anisotropy characteristics. Since the Cenozoic, the eastward compressive stress of the TP may have been transmitted to the eastern Sichuan fold belt through the hard SCB, leading to an uplift of the east Sichuan fold belt (1-2 mm/ a) in this area (Shi et al., 2016;Yuan et al., 2018). We believe that the azimuthal anisotropy of the upper crust is related to the high-angle stratigraphic directional arrangement caused by the eastern Sichuan fold belt deformation, which can be considered as the response evidence for the eastward extrusion of the TP. At 22-34 s periods ( Figures 9B, C), the FPD of the azimuthal anisotropy is mainly in the NW-SE direction, and the anisotropy amplitude is greater than 2%, basically consistent with the direction and magnitude of the SKS anisotropy . We deduce that the lithosphere has undergone significant deformation.

Crustal structure characteristics and deformation mechanism of the dianzhong block
At the period of 4-30 s, the area near PZH exhibits high phase velocity anomaly ( Figure 8C), which is in good agreement with the

FIGURE 7
Sensitivity kernels of Rayleigh wave phase velocity at different periods.

Frontiers in Earth Science
frontiersin.org characteristics of high Poisson's ratio, high wave velocity, high resistivity, high density, positive magnetic anomaly, and low terrestrial heat flow Wang et al., 2017;Wu et al., 2013;Cheng et al., 2017;Shen et al., 2015;Xu et al., 2015;Teng et al., 2019), all distributed in the inner core of the Emeishan large igneous rock province. Moreover, the FPDs in the vicinity of PZH have the characteristics of the NS direction at 4-30 s periods. Another notable feature is a large area of low-velocity anomalies near the Xiaojinhe fault (XJHF) zone in the north, which is consistent with the S-wave velocity obtained by Liu et al. (2014). The phase velocities are sensitive to the S-wave velocities of the midto-lower crust at 20-34 s periods (Figure 7). Within this period, the FPD near the northern XJHF zone gradually changed from NS to NE-SW direction, roughly parallel to the fault zone. Geochemical and regional geological studies (Xu and Zhong, 2001) indicated that during the period of 263-251Ma in the Late Permian, a famous tectonic activity occurred in the western YZB, the magmatic activity of the Emeishan great igneous rock province (Xu and Chung, 2001), whose dynamic mechanism may be related to mantle plume activity. Both sources of intracrustal magmatic intrusions and eruptive basalts come from the deep mantle. According to the dome structure and basalt thickness, the active center of the mantle plume is located near PZH (Shen et al., 2007;Xu and Zhong, 2001). We speculate that the high-velocity anomaly in the crust may be formed by the cooling and crystallization of a large number of high-density basic and ultramafic rock intrusions under the PZH area during the uplifting process of the early mantle plume (Wu et al., 2013;Yang et al., 2014;Zheng et al., 2016).
The XJHF, as the western boundary of the Yangtze block, obliquely cuts the Chuandian diamond block into two sub-blocks in CXB and DZB (Xiang et al., 2002;Xu et al., 2003). The topographic differences between the two sides of the XJHF fault Absolute isotropic phase velocities across the four profiles shown in Figure 8D (black lines). The red bars indicate that the fast-wave polarization direction (FPD) of azimuthal anisotropy shows a regional variation. the FPD near the LCJF in the west is NE-SW direction, the FPD near the WLSF zone in the middle is near the N-S direction, and the FPS near the RRF zone in the east is NW-SE direction. The depth corresponding to each period is calculated in terms of the sensitive kernel in Figure 7. Topography is depicted above each profile as the black fonts and the red fonts above it, marking the location of major faults along each profile. The red bars denote the tendency of anisotropic direction. The abbreviations for fault names (red) are the same as in Geological investigations have shown that differential movements on both sides of the fault zone are apparent. The horizontal slip rate values of the sub-block in CXB toward SE are 2 mm/a higher than those of the DZB. Since the Late Quaternary, the average differential upward and downward movement rates have reached 1.0-1.3 mm/a (He et al., 1993;Xu et al., 2003). Furthermore, there are also considerable differences in vertical movements observed by level (Liang et al., 2013). What causes the significant differential motion between the two sub-blocks within the Chandian diamond block? Seismic tomography results show  that there is a pronounced lower crustal flow in the CXB, and the flow direction is consistent with the extrusion direction of the block to the south. We suggest that the southward extrusion of the TP material is blocked by the hard intracrustal masses with PZH as the core, resulting in a rapid uplift of the northern topography.

FIGURE 9
Comparison of the previous seismic anisotropy results and our azimuthal anisotropy at different periods. (A) Azimuthal anisotropy of Rayleigh wave at period 8 s vs. the Local Seismic Anisotropy (LSA) (Shi et al., 2009Gao et al., 2012); (B) Azimuthal anisotropy of Rayleigh wave at period 22 s vs. the Pms phase splitting of receiver function (Sun et al., 2012;Chen et al., 2013;Cai et al., 2016); (C) Azimuthal anisotropy of Rayleigh wave at period 32 s vs. the SKS/SKKS phases splitting (Chang et al., 2015); (D) Distribution of Moho depth .The abbreviations of geological units are defined in the caption of Figure 1.

Frontiers in Earth Science
frontiersin.org

Conclusion
We present an azimuthally anisotropic phase velocity model of the crust and uppermost mantle beneath the southeastern TP from surface wave dispersion. A summary of our major findings and the structural features revealed by the azimuthally anisotropic model is as follows.
1. The azimuthal anisotropy and phase velocity in the SMB show noticeable regional variance. Subject to compressive stress in the southwest direction of the Chuandian diamond block, the crust of the eastern SMB underwent shear deformation, forming crustal anisotropy parallel to the strike of the RRF. Meanwhile, the crust of the western Simao block was blocked by the Lincang granite belt and strongly deformed, forming crustal anisotropy parallel to the strikes of the Lincang granite belt and the LCJF. Likewise, the central is in the transition zone of differential deformation on both sides, forming near N-S trending crustal anisotropy, which is almost consistent with the extensive development of folds and thrust structures in the Mesozoic-Cenozoic strata in the block. We infer that the compressive stress in the southwest direction of the Chuandian diamond block has a controlling effect on the crustal deformation of the SMB, and the west was blocked by the Lincang granite belt and experienced strong tectonic deformation.
2. In the middle RRF and its vicinity, the anisotropic FPD at 4-24 s periods is near N-S direction, showing a large angle with the fault strike. The anisotropic material in the upper crust of the middle RRF zone may not have a dominant arrangement along the fault strike but is related to the action of the regional tectonic stress field. 3. At 4-32 s periods, the Rayleigh surface wave phase velocities in the eastern and western sides of the SZMLF have noticeable velocity differences, which correspond well with the S-wave and P-wave velocity variations. Coupled with the complex tectonic stress at the junction of the blocks leading to prominent regional characteristics of the FPDs, we infer that the SZMLF zone may be the western boundary between the Yangtze and Cathaysia blocks in the crust. 4. Azimuthal anisotropy and Rayleigh surface-wave phase velocity show a pronounced lower crustal flow in the Chuanxibei subblock, and the flow direction is consistent with the extrusion direction of the block to the south. The southward extrusion of the Tibetan Plateau material is blocked by the hard intracrustal masses with Panzhihua as the core, resulting in a rapid uplift of the northern topography.

Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: the China Seismic Array Data Management Center at the Institute of Geophysics, China Earthquake Administration.

Author contributions
JL and LF designed the method and programmed the codes; all authors participated in writing, revising, and approving the final manuscript.

Funding
This work was supported by grants 42204070, U2039204 and U2139205 from the National Natural Science Foundation of China, 2021M702988 from the China Postdoctoral Science Foundation, and 2020-S-39 from guiding science and technology projects in Sanming City.

FIGURE 10
A cartoon summarizing the deformation pattern of the crust in the southeastern Tibetan Plateau. We estimate that the compressive stress in the southwest direction of the Tibetan Plateau material has a controlling effect on the crustal deformation of the Simao block, which likewise blocked by the Lincang granite belt, resulting in strong tectonic deformation. The anisotropy of the upper crust of the middle Red River fault zone is related to the regional tectonic stress field. And the Shizong-Mile fault zone may be the western boundary between the Yangtze block and the Cathaysia block. The red ellipses denote the FPD of the Rayleigh-wave phase velocity in the crust. The black arrow denotes the movement direction of the Dianzhong block relative to the South China block . The black lines denotes the faults.

Frontiers in Earth Science
frontiersin.org Publisher's note All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.