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.     

Extensive-Dilatancy Anisotropy (Eda) Inferred From Observations of Crustal Shear Waves Generated By A Refraction Experiment In Northern Scandinavia

Summary. Particle-motion plots of shear waves have been studied for the section FG of the FENNOLORA seismic experiment. Shear-wave splitting is observed on some records and the polarization of the first arriving shear waves show two peaks at about N35°W and N65°E. These results can be interpreted as being due to crack-induced anisotropy with the crack direction dominated by a (dominant) horizontal stress around N35°W. This is consistent with in situ stress measurements and focal mechanism studies in Scandinavia. the results show that seismic refraction experiments may be useful in providing evidence of crack-induced anisotropy in the stable continental crust.     

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.     

利用地震面波和重力资料联合反演地壳—上地幔三维密度结构的方法探讨

本文探讨了利用地震面波和重力资料联合反演地壳－地幔三维密度结构的反演问题。首先建立了地震面波和重力资料的观测方程,然后应用广义线性反线理论给出了反演问题的解。模型试验结果表明,与单纯面波反演的结果相比,联合反演的解在分析率和方差两个方面都有改善,当加入扰动重力数据时,在０－３００ｋｍ深度范围内联合反演的结果明显好于单纯面波反演的结果。     

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.     

Response of GPS occultation signals to atmospheric gravity waves and retrieval of gravity wave parameters

We show that the amplitude of the Global Positioning System (GPS) signals in the radio occultation (RO) experiments is an indicator of the activity of the gravity waves (GW) in the atmosphere. The amplitude of the GPS RO signals is more sensitive to the atmospheric wave structures than is the phase. Early investigations used only the phase of the GPS occultation signals for statistical investigation of the GW activity in the height interval 10–40 km on a global scale. In this study, we use the polarization equations and Hilbert transform to find the 1-D GW radio image in the atmosphere by analyzing the amplitude of the RO signal. The radio image, also called the GW portrait, consists of the phase and amplitude of the GW as functions of height. We demonstrate the potential of this method using the amplitude data from GPS/Meteorology (GPS/MET) and satellite mission Challenge Mini-satellite Payload (CHAMP) RO events. The GW activity is nonuniformly distributed with the main contribution associated with the tropopause and the secondary maximums related to the GW breaking regions. Using our method we find the vertical profiles of the horizontal wind perturbations and its vertical gradient associated with the GW influence. The estimated values of the horizontal wind perturbations are in fairly good agreement with radiosonde data. The horizontal wind perturbations v(h) are ±1 to ±5 m s with vertical gradients dv/dh ±0.5 to ±15 m s km at height 10–40 km. The height dependence of the GW vertical wavelength was inferred through the differentiation of the GW phase. Analysis of this dependence using the dispersion relationship for the GW gives the estimation of the projection of the horizontal background wind velocity on the direction of the GW propagation. For the event considered, the magnitude of this projection changes between 1.5 and 10 m s at heights of 10–40 km. We conclude that the amplitude of the GPS occultation signals contain important information about the wave processes in the atmosphere on a global scale.     

Condensed hyperelements method of non-vertical consistent boundaries for wave propagation analysis in irregular media

The study of wave propagation in finite/infinite media has many applications in geotechnical and structural earthquake engineering and has been a focus of research for the past few decades. This paper presents an analysis of 2D anti- plane problems （Love waves） and 2D in-plane problems （Rayleigh waves） in the frequency domain in media consisting of a near-field irregular and a far-field regular part. The near field part may contain structures and its boundaries with the far-field can be of any shape. In this study, the irregular boundaries of the near-field are treated as consistent boundaries, extending the concept of Lysmer＇s vertical consistent boundaries. The presented technique is called the Condensed Hyperelements Method （CHM）. In this method, the irregular boundary is limited to a vertical boundary at each end that is a consistent boundary at the far-field side. Between the two ends, the medium is discretized with hyperelements. Using static condensation, the stiffness matrix of the far-field is derived for the nodes on the irregular boundary. Examples of the application of the CHM illustrate its excellent accuracy and efficiency.