Crustal and upper mantle structure of stable continental regions in North America and northern Europe

From an analysis of many seismic profiles across the stable continental regions of North America and northern Europe, the crustal and upper mantle velocity structure is determined. Analysis procedures include ray theory calculations and synthetic seismograms computed using reflectivity techniques. TheP wave velocity structure beneath the Canadian Shield is virtually identical to that beneath the Baltic Shield to a depth of at least 800 km. Two major layers with a total thickness of about 42 km characterize the crust of these shield regions. Features of the upper mantle of these region include velocity discontinuities at depths of about 74 km, 330 km, 430 km and 700 km. A 13 km thickP wave low velocity channel beginning at a depth of about 94 km is also present.A number of problems associated with record section interpretation are identified and a generalized approach to seismic profile analysis using many record sections is described. TheS wave velocity structure beneath the Canadian Shield is derived from constrained surface wave data. The thickness of the lithosphere beneath the Canadian and Baltic Shields is determined to be 95–100 km. The continental plate thickness may be the same as the lithospheric thickness, although available data do not exclude the possibility of the continental plate being thicker than the lithosphere.     

Modeling crustal structure of the south-eastern Fennoscandia

We present new seismic velocity models of the crust and uppermost mantle along two refraction and wide-angle reflection profiles in the southern Fennoscandia: the Pribalt and 1-EB profiles. Some new results obtained along the Coast and the Baltic Sea profiles are also presented. The intercept time method and ray tracing are used for the modeling. The study shows that the lateral variations are small in the velocity structure of?? the crust up to the depth of 20?C25 km. The most significant lateral variations are observed in the Moho discontinuity topography and in the seismic velocities in the lower crust. In Paleoproterosoic Svekofennian domain, besides the well-known Moho depression in southern Finland, another Moho depression is revealed in the region from the Gotland Island to the Gulf of Riga. We suggest that this depression can correspond to the unknown crustal unit (we call it the Gotland-Riga belt). The Moho depth increases from the average of 40?C45 km to 55 km in this belt. The Moho depression is filled by the matter with velocities of 6.8?C7.1 km/s. Deep faults inclined to the north and strong variations of the mantle velocities are typical for the uppermost mantle of the Gotland-Riga belt.     

台湾海峡大容量气枪震源海陆联测初探

本文利用在我国台湾海峡采用大容量气枪震源开展海陆联测获得的广角地震测线HX9, 采用二维射线追踪法反演得到了HX9剖面的地壳二维速度结构和地壳界面形态, 初步探明了福建—台湾海峡海陆过渡带的深部构造． 结果表明: HX9剖面的地壳内存在两个速度间断面, 即C界面和莫霍面, 其中: C界面为上、 下地壳的分界面, 是一个小的速度不连续面, 速度变化值达0.08—0.16 km/s; 而地壳底部的莫霍面则有较大的速度反差, 变化值达1.02—1.29 km/s, 莫霍面上、 下的速度分别为6.75—6.97 km/s和8.00—8.07 km/s． 沿剖面的地壳界面形态总体起伏不大, 陆域上、 下地壳的厚度和界面变化趋势均相似, 从陆域到海域呈微倾斜变化趋势, 表现为减薄陆壳的特征． 莫霍面陆域埋深约为31.6 km, 向福建东南沿海逐渐减薄至27.4 km左右．      

The complexity of the crust from refraction/wide-angle reflection data

ThePmP wide-angle reflected signal from the Moho shows a wide variation in its characteristics when cornpared from region to region. If the earth's crust is simple and homogeneous and the Moho is a sharp discontinuity, supercriticalPmP wavelets are large and isolated in the 100 km to 300 km distance range when compared to the preceding signals on the seismic traces. If the crust itself has numerous short reflectors and the Moho transition zone is thick or badly disrupted from past tectonic disturbances, thenPmP is often poorly defined and difficult to identify in the coda of earlier arriving signals. These signals are poorly correlated because the reflectors themselves are in general not only discontinuous but tend to be distributed at various depths within the crust. The effect of the vertical velocity gradient in the crust is to make the reflective field for downgoing waves much greater than for upgoing waves. The large variations in reflection coefficients wiht angle of incidence has the effect of making the coda generated by a reflective machanism distance dependent. The reflected signals are also contaminated with scattered signals from smaller scale heterogeneities which may be distributed more uniformly within the crust.In conventional seismic inversion methods applied to crustal refraction experiments, emphasis is usually placed on obtaining a velocity model of the subsurface structure. These models are relatively simple and are limited in the amount of complexity which can be uniquely inferred from the data. The main problem arises because of the difficulties in the identification of the origin of the signals which are themselves often incoherent from trace to trace.In this paper it is shown how the conventional record section may be complemented with a normal moveout corrected intensity section which emphasizes areas of large signal complexities. Data from this intensity section is then used as input to obtain a quantitative measurement of a complexity parameter. These measurements may be used to infer or compare differences in crustal heterogeneity from one region to another. The discussion is illustrated with both numerical modeling data as well as data from recent crustal experiments which were conducted over the Canadian Shield.     

The crust-mantle transition and the Moho beneath the Vogtland/West Bohemian region in the light of different seismic methods

The structure of the crust and the crust-mantle boundary in the Vogtland/West Bohemian region have been a target of several seismic measurements for the last 25 years, beginning with the steep-angle reflection seismic studies (DEKORP-4/KTB, MVE-90, 9HR), the refraction and wide-angle experiments (GRANU’95, CELEBRATION 2000, SUDETES 2003), and followed by passive seismic studies (receiver functions, teleseismic tomography). The steep-angle reflection studies imaged a highly reflective lower crust (4 to 6 km thick) with the Moho interpreted in a depth between 30 and 32 km and a thinner crust beneath the Eger Rift. The refraction and wide-angle reflection seismic studies (CELEBRATION 2000) revealed strong wide-angle reflections in a depth of 26–28 km interpreted as the top of the lower crust. Long coda of these reflections indicates strong reflectivity in the lower crustal layer, a phenomenon frequently observed in the Caledonian and Variscan areas. The receiver function studies detected one strong conversion from the base of the crust interpreted as the Moho discontinuity at a depth between 27 and 37 km (average at about 31 km). The discrepancies in the Moho depth determination could be partly attributed to different background of the methods and their resolution, but could not fully explain them. So that new receivers function modelling was provided. It revealed that, instead of a first-order Moho discontinuity, the observations can be explained with a lower crustal layer or a crust-mantle transition zone with a maximum thickness of 5 km. The consequent synthetic ray-tracing modelling resulted in the model with the top of the lower crust at 28 km, where highly reflective lower crustal layer can obscure the Moho reflection at a depth of 32–33 km.     

Deep sounding with industrial power lines in conjunction with MTS measurements

This study is devoted to the development of a method for deep tensor electromagnetic inductive sounding of the Baltic Shield lithosphere with the use of mutually orthogonal industrial power lines alternately connected to a generator of up to 100 kW in power. The first results of the sounding are obtained at separations of up to 510 km from the source. The electrical conductivity of the Baltic Shield lithosphere in the crust-mantle transition zone (depths of 20–50 km) is numerically estimated from soundings with natural and controlled sources. The perspectives of future investigations are associated with synchronous observations on a 2-D network of measuring stations. This will enable the study of the relation of deep geoelectric parameters of the lithosphere to the geological structure of the crust in terms of 3-D models.     

Crustal structure of Northwest Zagros (Kermanshah) and Central Iran (Yazd and Isfahan) using teleseismic Ps converted phases

Receiver functions are widely employed to detect P-to-S converted waves and are especially useful to image seismic discontinuities in the crust. In this study we used the P receiver function technique to investigate the velocity structure of the crust beneath the Northwest Zagros and Central Iran and map out the lateral variation of the Moho boundary within this area. Our dataset includes teleseismic data (M _b ≥ 5.5, epicentral distance from 30° to 95°) recorded at 12 three-component short-period stations of Kermanshah, Isfahan and Yazd telemetry seismic networks. Our results obtained from P receiver functions indicate clear Ps conversions at the Moho boundary. The Moho depths were firstly estimated from the delay time of the Moho converted phase relative to the direct P wave beneath each network. Then, we used the P receiver function inversion to find the properties of the Moho discontinuity such as depth and velocity contrast. Our results obtained from PRF are in good agreement with those obtained from the P receiver function modeling. We found an average Moho depth of about 42 km beneath the Northwest Zagros increasing toward the Sanandaj-Sirjan Metamorphic Zone and reaches 51 km, where two crusts (Zagros and Central Iran) are assumed to be superposed. The Moho depth decreases toward the Urmieh-Dokhtar Cenozoic volcanic belt and reaches 43 km beneath this area. We found a relatively flat Moho beneath the Central Iran where, the average crustal thickness is about 42 km. Our P receiver function modeling revealed a shear wave velocity of 3.6 km/s in the crust of Northwest Zagros and Central Iran increasing to 4.5 km/s beneath the Moho boundary. The average shear wave velocity in the crust of UDMA as SSZ is 3.6 km/s, which reaches to 4.0 km/s while in SSZ increases to 4.3 km/s beneath the Moho.     

南北地震带南段地壳厚度重震联合最优化反演

重力反演方法是研究地壳结构和物性界面起伏的有效地球物理手段之一.本文收集了南北地震带南段67个已有的固定台站接收函数反演的Moho面深度结果,并使用基于EGM2008重力异常模型计算的布格重力异常,验证了本文提出的重震联合密度界面反演方法的有效性.利用接收函数对台站下方Moho面深度估计作为先验约束,定义了一类评价函数,通过对重力反演算法中尺度因子,平移因子和稳定性因子的最优选择,最小化重力反演结果与接收函数模型之间的差异.结果表明,本文提出的方法,可以有效地同化不同地球物理方法获得的反演模型,且通过重震联合反演可以改进由于对空间分布不均匀的接收函数结果插值可能而引起的误差.本文还通过引入Crust1.0的Moho面深度为初值,同时考虑地壳密度的横向不均匀分布,通过模型之间的联合反演有效改善了地球物理反演模型间的不一致性问题.本文反演得到的最优化Moho面深度模型与已知67个台站位置接收函数模型之间的标准差约1.9km,小于Crust1.0与接收函数结果模型之间标准差为3.73km的统计结果.本文研究结果对于同化重震反演结果、精化地壳密度界面模型,都具有十分重要的参考意义.     

大容量气枪震源陆地反射地震探测——以长江中下游铜陵地区为例

为了实验大容量气枪震源陆地水体流动激发反射地震探测效果,在长江中下游安徽省铜陵段,采用气枪船长江航道流动激发、沿江岸布设反射地震仪器接收的非纵弯线工作方式,得到了反映测线经过地区地壳深部结构和构造特征的反射地震数据。原始资料信噪比较低,但部分资料不同部位仍可辨认出来自地壳及莫霍面反射波组。就传播距离而言,地震波传播的水平距离最大可达21km,垂直深度可达30km以上。在数据处理中,根据原始资料特点,针对性采用了非纵弯线面元定义、三维层析静校正、叠前多域去噪及组合反褶积技术,最终得到的叠加时间剖面上具有丰富的壳内反射波组。结果显示,测线经过地区的地壳结构为双层结构,总厚度为30.0～36.0km。上地壳呈现隆坳相间的反射特征,下地壳存在多组叠层状弧型反射波组,莫霍面反射特征清晰,由2～3个反射同相轴组成,呈现SW端向NE段抬升的形态。剖面经过地区存在一个切穿下地壳和莫霍面的深部断裂,应该是长江深断裂的反映。研究结果充分说明,大容量气枪震源可应用于陆地流动水体地壳精细结构的深地震反射探测。     

Combined Gravimetric–Seismic Crustal Model for Antarctica

The latest seismic data and improved information about the subglacial bedrock relief are used in this study to estimate the sediment and crustal thickness under the Antarctic continent. Since large parts of Antarctica are not yet covered by seismic surveys, the gravity and crustal structure models are used to interpolate the Moho information where seismic data are missing. The gravity information is also extended offshore to detect the Moho under continental margins and neighboring oceanic crust. The processing strategy involves the solution to the Vening Meinesz-Moritz’s inverse problem of isostasy constrained on seismic data. A comparison of our new results with existing studies indicates a substantial improvement in the sediment and crustal models. The seismic data analysis shows significant sediment accumulations in Antarctica, with broad sedimentary basins. According to our result, the maximum sediment thickness in Antarctica is about 15 km under Filchner-Ronne Ice Shelf. The Moho relief closely resembles major geological and tectonic features. A rather thick continental crust of East Antarctic Craton is separated from a complex geological/tectonic structure of West Antarctica by the Transantarctic Mountains. The average Moho depth of 34.1 km under the Antarctic continent slightly differs from previous estimates. A maximum Moho deepening of 58.2 km under the Gamburtsev Subglacial Mountains in East Antarctica confirmed the presence of deep and compact orogenic roots. Another large Moho depth in East Antarctica is detected under Dronning Maud Land with two orogenic roots under Wohlthat Massif (48–50 km) and the Kottas Mountains (48–50 km) that are separated by a relatively thin crust along Jutulstraumen Rift. The Moho depth under central parts of the Transantarctic Mountains reaches 46 km. The maximum Moho deepening (34–38 km) in West Antarctica is under the Antarctic Peninsula. The Moho depth minima in East Antarctica are found under the Lambert Trench (24–28 km), while in West Antarctica the Moho depth minima are along the West Antarctic Rift System under the Bentley depression (20–22 km) and Ross Sea Ice Shelf (16–24 km). The gravimetric result confirmed a maximum extension of the Antarctic continental margins under the Ross Sea Embayment and the Weddell Sea Embayment with an extremely thin continental crust (10–20 km).     

云南思茅—中甸地震剖面的地壳结构

云南思茅—中甸宽角反射/折射地震剖面切割松潘—甘孜、扬子和华南三个构造单元的部分区域. 我们利用初至波和壳内反射波走时层析成像获得地壳纵波速度结构. 在获得新的地壳速度结构模型基础上，利用地震散射成像思想和低叠加次数的叠前深度偏移方法重建了研究区的地壳、上地幔反射结构. 综合分析研究区地壳P波速度模型和壳内地震反射剖面发现：沿测线从北至南地壳厚度从约50 km减薄至35 km左右，地壳厚度的减薄量主要体现在下地壳，剖面北段下地壳厚度约为30 km，剖面南段下地壳厚度仅为15 km左右；上地幔顶部局部位置P波速度值偏低，一般为76～78 km/s，反映出云南地区是典型的构造活动区的特点.剖面沿线地壳内地震反射发育，其中莫霍强反射出现在景云桥下方；在景云桥弧形断裂带8～10 km深处出现宽约50 km的强反射带.     

DEEP STRUCTURE OF NORTHERN HENAN PROVINCE AND ADJACENT AREAS DERIVED FROM GRAVITY AND SEISMIC SOUNDING DATA IN RELATION TO DISTRIBUTION OF EARTHQUAKES

We conduct the wave field separation of the gravity field for northern Henan Province and adjacent areas by the wavelet multi-scale decomposition method, and obtain multi-order gravity wavelet details and regional gravity field information. Then the Parker density surface inversion is used to invert the Moho interface. Based on the analysis of wavelet details in different orders and results of three seismic sounding profiles available in this area, we attempt to reveal the deep crustal structure of the study area. Research results show that the crustal structure is dominated by uneven density distribution accompanied by uplifts and depressions in the region with obvious heterogeneities of the density in horizontal and vertical directions. The gravity field characteristics in the middle-upper crust correspond to the surface topography, the lower crust is dominated by the large-scale high-low gravity anomalies, and several major depression basins show the characteristics of low velocity and low density. At the same time, the depth of the Moho interface changes greatly, which forms the block structure pattern of the regional crustal thickness. Among these features, the area with relatively large variations of the Moho is located in the transition zone of the basin to the Taihang Mountains, or exactly the Moho mutation belt. The Moho interface of the basin area as a whole is dominated by the uplift intertwined with local variations, of which the least and largest depths are 31km and 37km, respectively. Due to the gravity isostasy, the crustal thickness is larger(about 41km)in the northwest of the Taihang Mountains, with less average crustal density. In the study area, earthquakes tend to occur around the transition zone with density changes where the Moho is locally convex. The seismogenic mechanism may be associated with upwelling of upper mantle materials, low-velocity and low-density structures in the middle-lower crust and connection of deep large faults. Moreover, the deep large faults play a controlling role in the distribution of regional earthquakes.     

台湾-阿尔泰地学断面阿尔金-龙门山剖面的地壳纵波速度结构

通过分析阿尔金—龙门山地学断面的地震资料,建立了该剖面的地壳纵波速度结构。研究结果表明,阿尔金北侧的塔里木盆地地区莫霍面为50km,而在其南侧的祁连地块莫霍面突然加深至73km,在柴达木盆地莫霍面又抬升至58km左右,然后,在松潘甘孜地块莫霍面降至70km,并呈现为台阶状向龙门山方向抬升到60km左右,最低速层,而在其南部地区则没有低速层出现,推测低速层为地壳中部的局部熔融物质,阿尔金—龙门山剖面上的两个莫霍面坳陷区分别与祁连地块和松潘—甘孜地块上的两个莫霍面坳陷区相对应,指示出这个两个地块具有较深的山根,青藏高原北部的巨厚地壳很可能是由于中生代以来发生的印度板块与亚洲板块碰撞时受到来自东西及南北方向的挤压,使地壳缩短所致。     

南海北部地壳密度结构:基于约束三维重力反演

地壳结构的揭示是研究陆缘伸展机制的基础.尽管在南海北部陆缘已开展了大量地壳尺度的二维地震探测,但目前还存在许多覆盖空白的地区,这些调查所得到的结果无法提供地壳结构的区域视图.为了揭示南海北部的地壳结构,本研究将水深和沉积物厚度信息作为约束条件,对空间重力异常进行区域三维反演,并利用地震研究的结果来衡量反演结果的质量.沿地震测线的密度切片显示,重力反演结果与地震研究结果具有良好的一致性.本研究通过选择两个与地震研究结果最匹配的密度异常等值面分别作为康拉德面和莫霍面,获取了莫霍面深度和上、下地壳的厚度.根据假设的相应初始地壳厚度,本研究进一步计算了全地壳、上地壳和下地壳的拉张因子.通过与已发表的居里面深度比较,发现南海北部大部分地区的居里面深度均位于莫霍面之下,其中西沙海槽的拉张因子βw大于3.5并且缺乏岩浆活动,是地幔橄榄岩蛇纹石化的潜在区域.南海北部陆缘的地壳拉张因子显示其经历了伸展方向为128°和160°的两组张裂运动,分别对应于神弧运动和珠琼运动一幕及二幕,在张裂过程中应力场发生了顺时针旋转.此外,上、下地壳的拉张因子表明北部陆缘普遍存在正向和反向差异伸展,陆架区域表现为反向差异伸展,洋陆过渡带为正向差异伸展,推测这种正向和反向差异伸展可能是由下地壳流动导致的,由地壳厚度差异引起的横向梯度力、软流圈浮力和沉积物负载共同驱动.     

The crustal structure under Sanjiang and its dynamic implications: Revealed by seismic reflection/refraction profile between Zhefang and Binchuan,Yunnan

The Sanjiang area in southwest China is considered as a tectonic intersection belt between the Tethys-Alps and the western Pacific, and has endured three-phase evolution processes: Proto-Tethys,Paleo-Tethys and Meso-Tethys[1―4]. In this area, its tectonics and struc- ture are extremely complicated, and intensively extru-sive deformation and faults are widely developed[1―3]. For that, the area is considered as the ideal na- ture-laboratory to study the evolution of Paleo-Tethys and also …     

香港地区海陆地震联测及深部地壳结构研究

为了探明南海北部海陆过渡带的深部地壳结构,我们在香港外海域进行了一次海陆地震联洲的实验,利用固定地震台网远距离接收海上气枪信号,接收距离远达200多km,并利用此次实验的测线1剖面模拟得到了海陆过渡带的深部地壳速度结构.速度结构模型表明:研究区海陆过渡带的地壳结构非均匀性较明显,由陆至海沉积层有一个突然增厚的特点;莫霍面深度约为26～29 km,上地壳P波速度约为5.5～6.4 km/s,下地壳P波速度为6.5～6.9 km/s.在担杆列岛往海方向有一个低速破碎带,其上地壳P波速度为5.2～6.1 km/s,下地壳P波速度为6.2～6.4 km/s,结合野外地质调查的结果,推测它可能为滨海断裂带.在担杆列岛往陆方向香港和深圳之间的研究区域,莫霍面有较大起伏,可能与此处发育的海丰断裂有关.     

用接收函数方法研究汶川和芦山地震之间未破裂段的地壳结构

芦山与汶川地震之间存在约40 km的地震空区.震源区和地震空区的深部构造背景的研究对深入了解中强地震的深部孕育环境及地震空区的地震活动性具有重要科学意义.利用本小组布设的15个临时观测地震台以及21个芦山科考台站和21个四川省地震局固定台站记录的远震数据,用H-K叠加方法得到各个台站的地壳厚度和平均泊松比,并构建了接收函数共转换点(CCP)偏移叠加图像以及反演得到台站下方的S波速度模型.我们的结果揭示了震源区和地震空区地壳结构特征差异:(1)汶川震源区的地壳平均泊松比为~0.28;芦山震源区为~0.29;而地震空区处于泊松比变化剧烈的区域;(2)汶川地震与芦山地震的震源区以西下方的Moho面呈现深度上的突变(这与前人的研究成果基本一致),分别从~44 km突变到~59 km,~40 km突变到~50 km,而地震空区地壳平均厚度呈现渐变性变化;(3)地震空区Moho面下凹且具有低速的上地壳.综合一维S波速度结构和H-k以及CCP的初步结果,这可能显示汶川地震的发震断裂在深部方向上向西倾斜并形成切割整个地壳的大型断裂;芦山地震则可能是由于上、下地壳解耦引起的;而地震空区处于两种地震形成机制控制区域的过渡带中.     

East-west crustal structure and "down-bowing" Moho under the northern Tibet revealed by wide-angle seismic profile

We recognized 6 sets of reflecting P- and S-wave events from Moho and other interfaces within the crust, respectively, with the wide-angle seismic data acquired from 510 km-long Selincuo-Ya'anduo profile in the northern Tibet, fitted the observed events with forward modeling, and interpreted crustal structure of P- and S-wave velocities and Poisson's ratio under the profile. The results demonstrate that the crustal structure between Yarlungzangbo and Bangong-Nujiang sutures changes abruptly, and the crust is the thickest at the middle part of the profile with thickness of 80 km or more. The "down-bowing" Moho is the striking feature for the crustal variation along the west-east direction. The Moho uplifts with steps, and the uplifting rate westward is greater than that eastward. The heterogeneity of P- and S-wave velocities exists both vertically and horizontally, and one lower velocity layer (LVL) exists with the depth range of 27-34 km and the thickness range of 5-7 km. For the upper crust, Poisson's ratio is the lowest at the middle part of the profile; for the lower crust, the Poisson's ratio at the east segment is lower than that at west segment, which means that the crustal rigidity for the upper crust is different from the lower crust, and the lower crust under the east segment of the profile is more ductile. We infer that the substance in the lower crust endured eastward flow along with the collision between Eurasian and Indian plates, and the "down-bowing" Moho is attributable to the multi-phase E-W tectonic processes.     

用转换函数方法研究喜马拉雅地区速度结构

利用流动数字地震台网提供的三分量地震波形记录,应用转换函数及快速模拟退火算法对喜马拉雅山脉地区46个地震站下的地壳横波速度结构进行了反演,为进一步揭示青藏高原喜马拉雅山脉地区的动力学演化过程提供了可靠的地球物理证据.根据本文结果可清晰看到,喜马拉雅山脉地区作为当今地壳活动最活跃的地区,物质交换非常活跃,地下结构远远未达到平衡,地壳速度有很大差异,在板块边界处莫霍界面速度间断不是非常明显,自喜马拉雅南坡向高原腹地,地壳厚度大致从55 km增加到80 km;沿经度方向,莫霍面也有一定的起伏.通过研究得到另外一个证据是,在喜马拉雅的主中央逆冲断裂,由大陆碰撞产生的主要构造,其深度可能要大于80 km.     

Crustal Thickness and Velocity Ratio beneath National Seismic Stations in the Sichuan-Yunnan Area Based on Receiver Function Analysis

By using the teleseismic receiver function method, this paper analyzes the crustal thickness and v_P/v_S ratios beneath the 4 National seismic stations (KMI, TNC, CD2 and PZH) in the Sichuan-Yunnan area. This study gives the variance of Moho depths and velocity ratios of the 4 stations in different directions. The results show that the Moho depth beneath the Kunming station is around 50km, and the velocity ratio varies between 1.62 and 1.69. The thickness of crust and the velocity ratio do not change much with the direction. The crust beneath Tengchong station shows clear directivity, being 40.7km thick in the northeast and 49.7km thick in the southeast. The difference of the v_P/v_S values is remarkable between the two directions, reaching 0.2. The Chengdu station also has shallow Moho, about 40km, but is 8km deeper in the northeast and southwest and the velocity ratio has a change of 0.13 between the two directions. The crust beneath the Panzhihua station is stable. In all directions, the Moho depth is around 60km and the v_P/v_S ratio doesn't change significantly.