3-D surface wave group velocity distribution in Central-Southern Africa

Group velocities estimated from fundamental mode Love and Rayleigh waves are used in a tomography process in central-southern Africa. The waves were generated by eighteen earthquakes, which occurred along the East African Rift and recorded at BOSA, LBTB and SLR seismic stations in southern Africa. The group velocities from Love and Rayleigh waves were isolated using the Multiple Filter Technique (MFT) at the period range of 10 to 50 seconds. The tomography method developed by Ditmar and Yanovskaya (1987) and Yanovskaya and Ditmar (1990), was applied to calculate the lateral distribution of surface wave group velocities in central-southern Africa. The results of the tomographic inversion were plotted as distribution maps. In addition to the maps, I also produced two velocity cross-sections across the area of study. The velocity distribution maps show the regional tectonic units, though with poor resolution. The azimuthal bias of the surface wave paths is reflected in the distribution of the group velocities. The Moho depth appears to correlate with velocities at a period of about 30 s. A low velocity feature observed beneath the Zimbabwe craton implies a thickening upper asthenosphere and lithospheric thinning beneath the Zimbabwe craton. Also estimated was a shear wave velocity model beneath the Zimbabwe craton.     

Structure of the upper mantle in Asia from phase and group velocities of Rayleigh waves

Dispersion curves of phase velocities of Rayleigh waves are determined by the method of frequency-time analysis in a range of periods of 10–200 s from data of 43 interstation traces in Central Asia. Because the joint use of phase and group velocities significantly decreases the uncertainty in the determination of S wave velocity structures, the same traces were used for calculating group velocities from tomographic reconstructions obtained in [Yanovskaya and Kozhevnikov, 2003, 2006] and determining average velocity structures along these traces. The velocity structures were calculated by the Monte Carlo and linear inversion methods, which gave consistent results. Using velocity values obtained at fixed depths by the 2-D tomography method, lateral variations in velocities at these depths were estimated, which allowed us to construct smoothed vertical velocity structures at some points in the region. The resulting structures were used as initial approximations for constructing local velocity structures solely from previously obtained local dispersion curves of group velocities in the area (32°–56°N, 80°–120°E). Based on these structures, we mapped the lateral distribution of velocity variations at upper mantle depths of 75–400 km and along three vertical profiles. The inferred velocity variations are in good agreement with data on the tectonics of the region.     

利用地震噪声准实时监测短周期面波波速变化

依据由噪声信号提取面波格林函数的原理,利用福建省地震台网25个宽频带台站2007年7月2日—8月29日的脉动观测资料,得到了瑞雷波群速度在福建地区的分布,并将该结果作为下一步相对变化动态成像的背景速度分布．分布结果表明,福建地区周期为3—5s的瑞雷波群速度大致在2.9—3.1km/s之间,平均速度为3.0km/s,瑞雷波群速度分布呈北高南低的现象,这与福建地区北部多山、南部多平原盆地的地理环境有很好的吻合．而且,该群速度分布图在漳州盆地地区表现出一个非常明显的低速,这主要是因为受到盆地沉积层的影响．通过滑动窗（窗长为20d,步长1d）技术得到了观测区内周期约为3-5s的瑞雷波波速分布变化的连续图像;再利用扣除背景影响的技术,得到了2007年8月14日—2008年7月1日福建地区瑞雷波波速的相对变化时空动态图像．通过分析相对变化时空动态图像与该时间范围内发生的地震的对应关系,表明福建地区瑞雷波波速在多次网内中等强度的地震或震群（ML＞3.0）中均表现出震前波速升高,震后下降恢复的变化趋势．初步分析认为,这可能与震前整个地区受到的应力增大震后应力得到释放所导致的介质变化有关．      

Probabilistic inversion of Rayleigh-wave dispersion data: an application to Mt. Etna, Italy

We present a methodology for determining the elastic properties of the shallow crust from inversion of surface wave dispersion characteristics through a fully nonlinear procedure. Using volcanic tremor data recorded by a small-aperture seismic array on Mount Etna, we measured the surface waves dispersion curves with the multiple signal classification technique. The large number of measurements allows the determination of an a priori probability density function without the need of making any assumption about the uncertainties on the observations. Using this information, we successively conducted the inversion of phase velocities using a probabilistic approach. Using a wave-number integration method, we calculated the predicted dispersion function for thousands of 1-D models through a systematic grid search investigation of shear-wave velocities in individual layers. We joined this set of theoretical dispersion curves to the experimental probability density function (PDF), thus obtaining the desired structural model in terms of an a posteriori PDF of model parameters. This process allowed the representation of the objective function, showing the non-uniqueness of the solutions and providing a quantitative view of the uncertainties associated with the estimation of each parameter. We then compared the solution with the surface wave group velocities derived from diffuse noise Green’s functions calculated at pairs of widely spaced (~5–10 km) stations. In their gross features, results from the two different approaches are comparable, and are in turn consistent with the models presented in several earlier studies.     

Crustal Structure of the Korean Peninsula Using Surface Wave Dispersion and Numerical Modeling

Surface wave dispersion is studied to obtain the 1-D average velocity structure of the crust in the Korean Peninsula by inverting group- and phase-velocities jointly. Group velocities of short-period Rayleigh and Love waves are obtained from cross-correlations of seismic noise. Multiple-filter analysis is used to extract the group velocities at periods between 0.5 and 20 s. Phase velocities of Rayleigh waves in 10- and 50-s periods are obtained by applying the two-station method to teleseismic data. Dispersion curves of all group and phase velocities are jointly inverted for the 1-D average model of the Korean Peninsula. The resultant model from surface wave analysis can be used as an initial model for numerical modeling of observations of North Korean events for a velocity model appropriated to the Korean Peninsula. The iterative process is focused especially on the surface sedimentary layer in the numerical modeling. The final model, modified by numerical modeling from the initial model, indicates that the crust shear wave velocity increases with depth from 2.16 km/s for a 2-km-thick surface sedimentary layer to 3.79 km/s at a Moho depth of 33 km, and the upper mantle has a velocity of 4.70 km/s.     

Velocity structure of the upper mantle of the East European platform according to seismic noise data

Cross-correlation functions of noise are constructed on 119 interstation paths from seismic noise records at stations of Eastern Europe. Dispersion curves of the group velocity of Rayleigh waves obtained from the cross-correlation functions are used for constructing the three-dimensional distribution of the velocity of transverse waves on the East European platform and in adjacent regions by methods of surface-wave tomography. The mean velocity in the crust is minimum in the region of the Caspian depression and Black Sea basin (<3.3 km/s) and maximum in the Baltic shield area (>3.7 km/s). The upper mantle beneath the Baltic and Ukrainian shields is characterized by increased velocity and the absence of the asthenospheric layer. Reduced velocities are noted in the upper mantle of the Black Sea basin. A low-velocity anomaly in the shape of a vertical column is revealed at depths of 200–300 km in the central part of the Dnieper-Donets aulacogen, which confirms the existence of a paleorift in this region.     

影响面波勘探精度的因素探讨

分析了影响面波勘探精度的主要因素,针对计算面波相速度的精度问题,讨论了常用测量台间以相速度方法的优缺点及其适用范围,并提出通过测量台间格林函数进行相位校正,以得到两台之间的精解相位。以合成的高频理论地震图作为记录信号,并在理论地震图上叠加了２０％的随机噪声,利用这些方法分析理论信号,数字实验表明：在相干频率范围内,格林函数方法能够精确测量台间相速度,而互相关法和窄带通滤波互相关法所得结果较为离散。     

Surface Wave Group Velocity Tomography Imaging from Ambient Noise for Fujian Province and Its Adjacent Areas

Two-month continuous waveforms of 108 broadband seismic stations in Fujian Province and its adjacent areas are used to compute noise cross-correlation function (NCF). The signal quality of NCF is improved via the application of time-frequency phase weighted stacking. The Rayleigh and Love waves group velocities between 1s-20s are measured on the symmetrical component of the NCF with the multiple filter method. More than 5,000 Rayleigh wave dispersion curves and about 4,000 Love wave dispersion curves are obtained and used to invert for group velocity maps. This data set provides about 50km resolution that is demonstrated with checkerboard tests. Considering the off great circle effect in inhomogeneous medium, the ray path is traced based on the travel time field computed with a finite difference method. The inverted group velocity maps show good correlation with the geological features in the upper and middle crust. The Fuzhou basin and Zhangzhou basin showed low velocity on the short period group velocity maps. On the long period group velocity maps, the low velocity anomaly in the high heat flow region near Zhangzhou and clear velocity contrast across the Zhenghe-Dapu faults, which suggests that the Zhenghe-Dapu fault might be a deep fault.     

基于背景噪声经验格林函数的地震准确定位精度分析——以2008年甘肃武都地震为例

准确的震源位置能够为抗震减灾工作和地球内部结构研究提供关键信息.在震中附近观测台网密集且方位角覆盖良好的情况下,通过拟合多个台站P波、S波的观测到时,能够得到准确的震源位置.而在台网稀疏的地区,由于缺少可靠的三维速度结构模型,往往造成较大的地震定位误差.近年来发展了基于背景噪声经验格林函数（EGFs）对地震波形进行校正的重定位方法,能够有效地压制路径上复杂速度结构体的影响,提高了地震定位精度,为稀疏台网情形下地震准确定位研究提供了一个新思路.本文选取由InSAR观测到准确位置的2008年甘肃武都M_S5.5地震作为测试案例,对稀疏台网下基于背景噪声格林函数地震准确震中测定方法进行了定量评估.利用震中附近多个流动台分别作为参考台,提取其与固定台站之间的背景噪声格林函数（EGFs）,然后使用噪声Rayleigh面波格林函数对地震波形进行校正,重定位武都地震震中,并与真实震中位置进行对比.结果表明：使用距离震中30 km以内的参考台,利用10~30 s频段的面波走时信息,噪声叠加时长为一个月,定位精度在5 km以内;当噪声叠加时间一年以上,重定位精度优于1 km.本文进一步针对缺少近台以及台网更加稀疏的情况进行了测试,发现使用10个左右固定台进行重定位,基于较高质量的噪声EGFs频散数据,定位结果精度可达3 km,从而给出了该方法的高精度定位所需的固定台网及参考台站的观测指标体系.

     

Rayleigh wave phase velocities of South China block and its adjacent areas

Using records of continuous seismic waveforms from 609 broadband seismic stations in the South China Block and its adjacent areas in 2010–2012, empirical Green's functions of surface waves were obtained from cross-correlation functions of ambient noise data between these stations. High quality phase velocity dispersion curves of Rayleigh waves were obtained using time-frequency analysis. These interstation dispersion curves were then inverted to build Rayleigh wave phase velocity maps at periods of 6–50 s. The results of phase velocity maps indicate that phase velocities at 6–10 s periods are correlated with the geological features in the upper crust. Major basins and small-scale grabens and basins display slow velocity anomalies; while most of the orogenic belts and the fold belts display high velocity anomalies. With the gravity gradient zone along Taihang Mountain to Wuling Mountain as the boundary for the phase velocity maps at period of 20–30 s, the western area mainly displays low velocity anomalies, while the eastern side shows high velocity anomalies. Phase velocities in the eastern South China Block south to the Qinling-Dabie orogenic belt is higher than that in the eastern North China Block to the north, which is possibly due to the differences of tectonic mechanisms between the North China Craton and the South China Block. The phase velocities at periods of40–50 s are possibly related to the lateral variations of the velocity structure in the lower crust and upper mantle: The low-velocity anomalies in the eastern part of the Tibetan Plateau are caused by the thick crust; while the Sichuan Basin and the southern part of the Ordos Basin display distinct high-velocity anomalies, reflecting the stable features of the lithosphere in these blocks. The lateral variation pattern of phase velocities in the southern part of the South China Block is not consistent with the surface trace of the block boundary in the eastern Yunnan Province and its vicinities. The phase velocities in the Sichuan Basin are overall slow at short periods and gradually increase with period from the central part to the edge of the basin, indicating the features of shallower basement in the center and overall stable lithospheric mantle of the basin. The middle and upper crust of the southern Ordos Basin in the North China Block is heterogeneous, while in lower crust and the uppermost mantle the phase velocities mainly exhibit high anomalies. High-velocity anomalies are widespread at the middle of the Qinling-Dabie orogenic belt, as well as the areas in southeastern Guangxi with Caledonian granite explosion, but its detailed mechanism is still unclear.     

Surface wave tomography method based on data from distant earthquakes

The standard method of surface wave tomography assumes the use of data on surface wave velocities measured on paths wholly located inside the region being investigated and oriented in different directions. If the data on group velocities obtained in the regional network of stations are used for this purpose, then the earthquake sources also should be situated within the limits of this region. If the region is small, then this requirement restricts the range of periods of surface waves and, correspondingly, the depth of research. The use of data from distant earthquakes enlarges the range of periods. However, it does not enable one to use the usual tomography method because of an inappropriate configuration of the system of paths: outside of the region in which the stations are located, the paths do not intersect, but, at the same time, the paths from a particular earthquake source toward the different stations of networks are very close over a large distance, and the mean correction to the velocity in the path sections outside of the network of stations can be accepted as identical. This assumption forms the basis of the proposed modification of the method of surface wave tomography, in which group velocities measured at stations of the local network from distant earthquakes are used. The lateral variations in the velocity are determined within the limits of the network, and the mean corrections to the velocity are determined on the paths from different earthquake sources in the sections outside of the network. As earlier, the condition of smoothness is imposed on the distribution of velocity variations, and the condition of smallness of the sum of their squares is imposed on the values of the corrections to the velocity. Testing of this method based on model examples and on real data obtained in the network of stations in Tibet showed its advantage over the standard tomography method.     

Surface wave tomography of the crust and upper mantle of Chinese mainland and its neighboring region

Introduction The three-dimensional S wave velocity of Chinese mainland and its neighboring region in-verted by surface wave dispersion data plays an important role in studying the lateral variation of lithosphere and geodynamic process, and understanding the forming and evolution of Chinese mainland and the relationship between shallow and deep structures. The three-dimensional veloc-ity structures of China and its major tectonic blocks were respectively studied by SONG, et al (1993), ZHOU…     

Seismic crust structure beneath the Aegean region in southwest Turkey from radial anisotropic inversion of Rayleigh and Love surface waves

The Turkish plate is covered by hundreds of accelerometer and broadband seismic stations with less than 50 km inter-station distance providing high-quality earthquake recordings within the last decade. We utilize part of these stations to extract the fundamental mode Rayleigh and Love surface wave phase and group velocity data in the period range 5–20 s to determine the crust structure beneath the Aegean region in southwest Turkey. The observed surface wave signals are interpreted using both single-station and two-station techniques. A tomographic inversion technique is employed to obtain the two-dimensional group velocity maps from the single-station group velocities. One-dimensional velocity–depth profiles under each two-dimensional mesh point, which are jointly interpreted to acquire the three-dimensional image of the shear-wave velocities underneath the study area, are attained by utilizing the least-squares inversion technique, which is repeated for both Rayleigh and Love surface waves. The isotropic crust structure cannot jointly invert the observed Rayleigh and Love surface waves where the radial anisotropic crust better describes the observed surface wave data. The intrusive magmatic activity related to the northward subducting African plate under the Turkish plate results the crust structure deformations, which we think, causing the observed radial anisotropy throughout complex pattern of dykes and sills. The magma flow resulting in the mineral alignment within dykes and sills contributes to the observed anisotropy. Due to the existence of dykes, the radial anisotropy in the upper crust is generally negative, i.e., vertically polarized S-waves (Vsv) are faster than horizontally polarized S-waves (Vsh). Due to the existence of sills, the radial anisotropy in the middle-to-lower crust is generally positive, i.e., horizontally polarized S-waves (Vsh) are faster than vertically polarized S-waves (Vsv). Similar radial anisotropic results to those of the single-station analyses are obtained by the two-station analyses utilizing the cross-correlograms. The widespread volcanic and plutonic rocks in the region are consistent with the current seismic interpretations of the crustal deformations.     

Joint inversion of teleseismic P waveforms and surface-wave group velocities from ambient seismic noise in the Bohemian Massif

Joint inversion of teleseismic P-waveforms and local group velocities of surface waves retrieved from ambient seismic noise has been performed to model velocity structure of the crust and uppermost mantle of the Bohemian Massif. We analysed P-waveforms of 381 teleseismic earthquakes recorded at 54 broadband seismic stations located on the territory of the Czech Republic and in its close surroundings. Group velocities of Rayleigh and Love surface waves were obtained by cross-correlating long-term recordings of seismic noise. The basis for waveform inversion is the well-known methodology of P-to-S receiver functions constructed from converted phases. Due to instabilities in direct inversion of receiver functions caused by the necessity of applying deconvolution, we propose an alternative formulation to fit observed and calculated radial components of P waveforms. The joint inversion is transformed into a search for the minimum of the cost function defined as a weighted sum of waveform and group velocity misfits. With the use of the robust stochastic optimizer (Differential Evolution Algorithm), neither derivatives nor a starting model are needed. The task was solved for 1D layered isotropic models of the crust and the uppermost mantle. We have performed a sequence of inversions with models containing one, two, three and four layers above a half-space. By using statistical criteria (F-test) we were able to select the simplest velocity models satisfying data and representing local geological structures. Complex crustal models are typical for stations located close to boundaries of major tectonic units. The relatively low average P to S wave-velocity ratio is in agreement with the generally accepted view that the BM crust is predominantly felsic.     

<Emphasis Type="Italic">S</Emphasis>-wave velocity model for several regions of the Kamchatka Peninsula from the cross correlations of ambient seismic noise

The data from the seismic networks of the Kamchatka Branch of the Geophysical Survey of the Russian Academy of Sciences are used for calculating the cross correlations of seismic noise for the stationary digital stations over 2013 and for radio telemetric stations (RTS) in the region of the Klyuchevskoy volcano over the period from January 1, 2009 to May 31, 2013. Four hundred and two correlations overall are calculated. The fundamental-mode group velocities of the Rayleigh waves are calculated in the periods ranging from 5 to 50 s. The calculations for the region of the Klyuchevskaya group of volcanoes are based on the RTS data and cover the periods from 2 to 8 s. The two-dimensional (2D) maps of group velocity distributions in different periods are constructed with the use of the algorithm of surface wave tomography (Barmin, 2001). The velocity sections for the selected Kamchatka regions are reconstructed by the dispersion curve inversion technique (Mordret, 2014). For each region, the structure of the Earth’s crust and upper mantle down to a depth of 50 km was obtained.     

THREE-DIMENSIONAL S-WAVE VELOCITY DISTRIBUTION BASED ON AMBIENT NOISE ANALYSIS IN EASTERN NORTH

We apply ambient noise tomography to continuous vertical component broadband seismic data between January 1, 2010 and December 31, 2011from the regional networks of 190 stations deployed by China Earthquake Administration in Hebei, Shanxi and Inner Mengolia. Ambient noise cross-correlations were performed to produce the Green's functions of each station-pair. Firstly, we used the multiple-filter analysis method to extract surface wave group and phase velocity dispersion curves from inter-station paths at periods from 7 to 40s. Then the study area was discretized into a 0.2°×0.2° grid to obtain the group and phase velocity distributions using O'ccam inversion method. After that, three dimensional (3-D) S-wave velocity structures from the surface down to 50km are inverted from group and phase velocities dispersion results. the results of S wave velocity distribution maps generally demonstrate good correlations with surface geological and tectonic features, and they also clearly revealed the lateral velocity variation in the crust. In the mid-upper crust, the basins are clearly resolved with low S wave velocity due to its thick sedimentary layer, and the Taihang and Yanshan uplifts show relative higher S wave velocity distribution. With the increase of depth (>30km), the S wave velocity distribution presents a contrary characteristic compared to that of the shallow layer, and the S wave velocity beneath the Taihang and Yanshan uplifts are much lower than basin areas, which is possibly correlated with the thickness of the crust. 3-D S wave velocity shows a low-velocity zone at~10~20km depth observed beneath the Tanshan-Hejian-Xintai-Cixian belt and Bohai Bay. the low-velocity zone at~20~30km depth beneath the Datong area may be associated with the thermal material in the crust-mantle. Our S wave velocity distribution maps clearly show that Taihang Mountains is not only the boundary of topography and tectonic zone, but also the transition zone of high and low S wave velocity.     

Anisotropic models of the upper mantle

Long period Rayleigh wave and Love wave dispersion data, particularly for oceanic areas, have not been simultaneously satisfied by an isotropic structure. In this paper available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. We assume that the surface wave data represents an azimuthal average of actual velocities. Thus, we can treat the mantle as transversely isotropic. The resulting models for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-velocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studies. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SH > SV in the vicinity of the low-velocity zone. Neat the East Pacific Rise, however, SV > SH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of ～50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures. A new model for average Earth is obtained which includes a thin LID. This model extends the fit of a PREM, type model to shorter period surface waves.     

Upper mantle structure beneath central Western Europe from data on phase velocities of surface waves

The deep structure of the upper mantle is determined from data on phase velocities of Love and Rayleigh waves measured by a differential method on traces between two stations in central Western Europe. One-dimensional velocity structures are first constructed from data of each pair of stations, after which two-dimensional distributions of SH and SV velocities are calculated by the method of two-dimensional tomography from S wave velocities at fixed depths. The results are presented in the form of 2-D vertical structures of the average S wave velocity (S = (SV + SH)/2) constructed along profiles crossing the region in directions of the best resolution. The main structural features are a higher velocity zone at depths of 60–80 km in the area (48°–50°N, 9°–11°E) and a lower velocity zone in the western part of the region at depths of 100–150 km, probably extending farther beyond the studied area.     

南海瑞雷面波群速度层析成像及其地球动力学意义

南海处于欧亚板块、 菲律宾海板块、 太平洋板块和印度-澳大利亚板块的交汇处, 其地质和构造作用十分复杂．通过面波群速度成像, 给出了南海及邻区的三维横波速度分布并分析了其地球动力学意义．南海西部和南部新布设的地震台站使得利用单台法时路径覆盖比过去更好. 特别是在华南地区, 新的台站分布能够弥补该地区地震少且台站少造成的射线密度不够的缺点． 首先运用多重滤波法得到南海周边48个台站周期为14——130 s范围内的基阶瑞雷波频散曲线图; 接着通过子空间反演得到整个区域在不同周期时的群速度分布; 最后通过阻尼最小二乘反演得到不同深度切片上的横波速度分布及不同纵剖面上的横波速度分布． 结果显示: ① 海盆速度较高, 且速度分布很好地勾勒出海盆的轮廓. 浅层较高的横波速度说明海盆都具有洋壳性质, 而深部较高的横波速度则可能对应扩张中心生成洋壳后残留的高速物质． 不同海盆速度上的差异与它们的热流值和年龄大小一致．海盆下的高速异常在60 km以下消失, 且在一定深度范围内由低速区替代． 在低速区下200 km深度, 在南海海盆观测到一条NE-SW走向的高速异常, 可能与古俯冲带有关. ② 环南海出现明显的高速区, 对应俯冲带特征, 且这些高速区速度差异明显且有间断, 说明俯冲带的非均质性和俯冲角度的差异. ③ 在环南海高速区内侧（向南海侧）观测到不连续的低速区． 在浅层, 这些低速区反映了沉积层和地壳的厚度特征. 在地幔, 这些低速区可能对应于古太平洋俯冲带的地幔楔或者也可能反映了南海海盆停止扩张后残留的地幔熔融物质. ④ 南海海盆岩石圈的厚度为60——85 km．      

华北克拉通中部地区背景噪声成像

利用华北地区的流动与固定地震观测台阵在2007年1月~12月期间的垂直分量记录,采用背景噪声层析成像方法获得了华北克拉通中部地区6~40 s的相速度和群速度分布图像.成像结果显示研究区地壳以及上地幔的速度结构存在很大的横向不均匀性.8~16 s的相速度以及8~20 s的群速度分布图像清晰地勾绘出华北地区盆地和隆起之间的边界.华北盆地以及山西断陷带呈明显的低速异常,表明相应区域存在较厚的沉积层并且在中、上地壳内存在低速带;而太行山隆起以及鄂尔多斯高原则呈现高速异常.相速度16~25 s以及群速度20~30 s的图像显示华北盆地由低速异常变为高速异常,表明华北盆地的地壳较薄.长周期(30~35 s)相速度图像表明华北盆地的上地幔顶部存在低速异常.华北盆地较薄的地壳和低速异常(明显低于全球平均水平),可能与该地区岩石圈减薄,软流圈物质上涌有关.8~35 s的群速度和相速度图像都显示,大同火山区是一个低速异常区,可能是由目前仍在活动的岩浆的加热而引起地壳升温造成的.