Crustal structure of Iraq from receiver functions and surface wave dispersion: implications for understanding the deformation history of the Arabian–Eurasian collision

We report the crustal structure for two locations in Iraq estimated by joint inversion of P -wave receiver functions (RFs) and surface (Rayleigh) wave group velocity dispersion. RFs were computed from teleseismic recordings at two temporary broad-band seismic stations located in Mosul (MSL) in the Zagros Fold Belt and Baghdad (BHD) in the Mesopotamian Foredeep. Group velocity dispersion curves at the sites were derived from continental-scale tomography. The inversion results show that the crustal thicknesses are 39 km at MSL and 43 km at BHD. We observe a strong Ps _Moho at BHD consistent with a sharp Moho discontinuity. However, at MSL we observe a weak Ps _Moho suggesting a transitional Moho where crustal thickening is likely to be occurring in the deep crust. Both sites reveal low velocity surface layers consistent with sedimentary thickness of about 3 km at station MSL and 7 km at BHD and agreeing well with the previous reports. Ignoring the sediments, the crystalline crustal velocities and thicknesses are remarkably similar at both stations. The similarity of crustal structure suggests that the crust of the northeastern proto-Arabian Platform was uniform before subsidence and deposition of the sediments in the Cenozoic. If crystalline crustal structure is uniform across the northern Arabian Platform then crustal thickness variations in the Zagros Fold Belt and Thrust Zone should reveal the history of deformation and crustal shortening in the Arabian–Eurasian collision zone and not reflect pre-existing crustal thickness variations in the Arabian Plate.     

On crustal corrections in surface wave tomography

Mantle models from surface waves rely on good crustal corrections. We investigated how far ray theoretical and finite frequency approximations can predict crustal corrections for fundamental mode surface waves. Using a spectral element method, we calculated synthetic seismograms in transversely isotropic PREM and in the 3-D crustal model Crust2.0 on top of PREM, and measured the corresponding time-shifts as a function of period. We then applied phase corrections to the PREM seismograms using ray theory and finite frequency theory with exact local phase velocity perturbations from Crust2.0 and looked at the residual time-shifts. After crustal corrections, residuals fall within the uncertainty of measured phase velocities for periods longer than 60 and 80 s for Rayleigh and Love waves, respectively. Rayleigh and Love waves are affected in a highly non-linear way by the crustal type. Oceanic crust affects Love waves stronger, while Rayleigh waves change most in continental crust. As a consequence, we find that the imperfect crustal corrections could have a large impact on our inferences of radial anisotropy. If we want to map anisotropy correctly, we should invert simultaneously for mantle and crust. The latter can only be achieved by using perturbation theory from a good 3-D starting model, or implementing full non-linearity from a 1-D starting model.     

Seismic motion in urban sites consisting of blocks in welded contact with a soft layer overlying a hard half-space

We address the problem of the response to a seismic wave of an urban site consisting of   N _b   blocks overlying a soft layer underlain by a hard substratum. The results of a theoretical analysis, appealing to a space–frequency mode-matching (MM) technique, are compared to those obtained by a space–time finite-element (FE) technique. The two methods are shown to give rise to the same prediction of the seismic response for   N _b = 1  , 2 and 40 blocks. The mechanism of the interaction between blocks and the ground, as well as that of the mutual interaction between blocks, are studied. It is shown, in the first part of this paper, that the presence of a small number of blocks modifies the seismic disturbance in a manner which evokes qualitatively, but not quantitatively, what was observed during the 1985 Michoacan earthquake in Mexico City. Anomalous earthquake response at a much greater level, in terms of duration, peak and cumulative amplitude of motion, is shown, by a theoretical and numerical analysis in the second part of this paper, to be induced by the presence of a large (≥10) number of identical equi-spaced blocks that are present in certain districts of many cities.     

Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis – II. Crustal and upper-mantle structure

We determine the 3-D shear wave speed variations in the crust and upper mantle in the southeastern borderland of the Tibetan Plateau, SW China, with data from 25 temporary broad-band stations and one permanent station. Interstation Rayleigh wave (phase velocity) dispersion curves were obtained at periods from 10 to 50 s from empirical Green's function (EGF) derived from (ambient noise) interferometry and from 20 to 150 s from traditional two-station (TS) analysis. Here, we use these measurements to construct phase velocity maps (from 10 to 150 s, using the average interstation dispersion from the EGF and TS methods between 20 and 50 s) and estimate from them (with the Neighbourhood Algorithm) the 3-D wave speed variations and their uncertainty. The crust structure, parametrized in three layers, can be well resolved with a horizontal resolution about of 100 km or less. Because of the possible effect of mechanically weak layers on regional deformation, of particular interest is the existence and geometry of low (shear) velocity layers (LVLs). In some regions prominent LVLs occur in the middle crust, in others they may appear in the lower crust. In some cases the lateral transition of shear wave speed coincides with major fault zones. The spatial variation in strength and depth of crustal LVLs suggests that the 3-D geometry of weak layers is complex and that unhindered crustal flow over large regions may not occur. Consideration of such complexity may be the key to a better understanding of relative block motion and patterns of seismicity.     

Surface wave tomography for azimuthal anisotropy in a strongly reduced parameter space

Large scale seismic anisotropy in the Earth's mantle is likely dynamically supported by the mantle's deformation; therefore, tomographic imaging of 3-D anisotropic mantle seismic velocity structure is an important tool to understand the dynamics of the mantle. While many previous studies have focused on special cases of symmetry of the elastic properties, it would be desirable for evaluation of dynamic models to allow more general axis orientation. In this study, we derive 3-D finite-frequency surface wave sensitivity kernels based on the Born approximation using a general expression for a hexagonal medium with an arbitrarily oriented symmetry axis. This results in kernels for two isotropic elastic coefficients, three coefficients that define the strength of anisotropy, and two angles that define the symmetry axis. The particular parametrization is chosen to allow for a physically meaningful method for reducing the number of parameters considered in an inversion, while allowing for straightforward integration with existing approaches for modelling body wave splitting intensity measurements. Example kernels calculated with this method reveal physical interpretations of how surface waveforms are affected by 3-D velocity perturbations, while also demonstrating the non-linearity of the problem as a function of symmetry axis orientation. The expressions are numerically validated using the spectral element method. While challenges remain in determining the best inversion scheme to appropriately handle the non-linearity, the approach derived here has great promise in allowing large scale models with resolution of both the strength and orientation of anisotropy.     

Seismic vertical array analysis for phase decomposition

We propose a vertical array analysis method that decomposes complex seismograms into body and surface wave time histories by using a velocity structure at the vertical array site. We assume that the vertical array records are the sum of vertically incident plane P and S waves, and laterally incident Love and Rayleigh waves. Each phase at the surface is related to that at a certain depth by the transfer function in the frequency domain; the transfer function is obtained by Haskell's matrix method, assuming a 1-D velocity structure. Decomposed P , S and surface waves at the surface are estimated from the vertical array records and the transfer functions by using a least-squares method in the frequency domain; their time histories are obtained by the inverse Fourier transform. We carried out numerical tests of this method based on synthetic vertical array records consisting of vertically incident plane P and S waves and laterally incident plane Love and Rayleigh waves. Perfect results of the decomposed P , S , Love and Rayleigh waves were obtained for synthetic records without noise. A test of the synthetic records in which a small amount of white noise was added yielded a reasonable result for the decomposed P , S and surface waves. We applied this method to real vertical array records from the Ashigara valley, a moderate-sized sedimentary valley. The array records from two earthquakes occurring at depths of 123 and 148 km near the array (epicentral distance of about 31 km) exhibited long-duration later phases. The analysis showed that duration of the decomposed S waves was a few seconds and that the decomposed surface waves appeared a few seconds after the direct S -wave arrival and had very long duration. This result indicated that the long-duration later phases were generated not by multireflected S waves, but by basin-induced surface waves.     

1961-2010年青海省人体舒适度指数时空分布特征

根据青海省50个气象观测站1961-2010年逐日气象资料, 对人体舒适度指数时空分布特征及其影响因子的权重进行了统计分析. 结果表明: 青海省人体舒适度主要为寒冷、冷、凉、凉爽和舒服等级, 整体呈现冷凉特征, 各区各等级年均日数分布差异较大. 青海省最不舒适的月份是1月份, 其次是12月; 最舒适的月份是7月, 其次是8月、6月. 1961-2010 年青海省人体舒适度指数上升趋势极显著, 与青藏高原气候变化趋势一致. 各区人体舒适度均呈显著上升趋势, 其中, 柴达木地区上升最明显; 春、秋季东部农业区和柴达木地区相似, 环青海湖地区和青南牧区相似, 各区秋季上升趋势均高于春季; 青南牧区各季人体舒适度年际波动幅度较大, 冬季环青海湖地区和青南牧区人体舒适度差异不大. 青海省人体舒适度季节差异明显, 夏季最高, 冬季最低, 春季略高于秋季, 近50 a来四季人体舒适度均呈显著上升趋势, 冬季波动幅度较大.温度(湿度、风速)与人体舒适度指数存在极显著的正(负)相关关系, 温度是影响人体舒适度指数的最主要因子, 风速和湿度主要通过温度影响人体舒适度指数, 且风速的负影响略大于湿度. 东部农业区和柴达木地区风速的负影响较大, 环青海湖地区和青南牧区相对湿度的负影响较大.