Tomographic Imaging of Lg and Sn Propagation in the Middle East

?—?Observations based on relatively limited data recorded by sparsely distributed stations have indicated that regional seismic phase propagation (Lg and Sn) is very complex in the Middle East. Accurate characterization of regional seismic wave propagation in this region necessitates the use of a large number of seismic stations. We have compiled a large data set of regional and local seismograms recorded in the Middle East. This data set comprises approximately four years of data from national short-period networks in Turkey and Syria, data from temporary broadband arrays in Saudi Arabia and the Caspian Sea region, and data from GSN, MEDNET, and GEOFON stations in the Middle East. We have used this data set to decipher the character and pattern of regional seismic wave propagation. We have mapped zones of blockage as well as inefficient and efficient propagation for Lg, Pg, and Sn throughout the Middle East. Two tomographic techniques have been developed in order to objectively determine regions of lithospheric attenuation in the Middle East.¶We observe evidence of major increase in Lg attenuation, relative to Pg, across the Bitlis suture and the Zagros fold and thrust belt, corresponding to the boundary between the Arabian and Eurasian plates. We also observe a zone of inefficient Sn propagation along the Dead Sea fault system which coincides with low Pn velocities along most of the Dead Sea fault system and with previous observations of poor Sn propagation in western Jordan. Our observations indicate that in the northern portion of the Arabian plate (south of the Bitlis suture) there is also a zone of inefficient Sn propagation that would not have been predicted from prior measurements of relatively low Pn velocities. Mapped high attenuation of Sn correlates well with regions of Cenozoic and Holocene basaltic volcanism. These regions of uppermost mantle shear-wave attenuation most probably have anomously hot and possibly thin lithosphere.     

Lg Coda Q and its Relation to the Structure and Evolution of Continents: A Global Perspective

—Tomographic maps of Lg coda Q (Q ^c _Lg) variation are now available for nearly the entire African, Eurasian, South American, and Australian continents, as well as for the United States. Q ^c _Lg at 1 Hz (Q ₀) varies from less than 200 to more than 1000 and Q ^c _Lg frequency dependence (<eta>) varies between 0.0 and nearly 1.0. Q ₀ appears to increase in proportion to the length of time that has elapsed since the most recent major episode of tectonic or orogenic activity in any region. A plot of Q ₀ versus time since that activity indicates that a single Q ₀-time relation approximates most mean Q ₀ values. Those that deviate most from the trend lay in Australia, the Arabian Peninsula, and the East African rift. The increase in Q ₀ with time may be due to a continual increase in crustal shear wave Q (Q ) caused by the loss of crustal fluids and reduction of crustal permeability following tectonic or orogenic activity. Extrapolated values of Q ^c _Lg at 5 Hz (using Q ₀ and <eta> values measured at 1 Hz and assuming that <eta> is constant in all regions between 1 and 5 Hz) show a similar percentage-wise increase with times that has elapsed since the most recent activity. Other factors that can reduce Q ₀ in continental regions include thick accumulations of sediment (especially sandstone and shale of Mesozoic age and younger), severe velocity gradients at the crust-mantle transition and, possibly, lateral variations in the depth, thickness, and severity of those gradients. Severe and large increases of Q in the mid-crust of some regions can cause relatively large values of <eta>, even if the frequency dependence of Q is small.     

Lg Wave Attenuation in the Isparta Angle and Anatolian Plateau (Turkey)

We estimate Lg wave attenuation using local and regional seismic phases in the Isparta Angle and the Anatolian Plateau (Turkey). The Isparta Angle (IA) is a tectonically active zone forming the boundary between the African Plate and the Anatolian Plateau, and is currently undergoing N–S extensional deformation. The Anatolian Plateau contains many intra-continental faults including the North Anatolian Fault Zone and the East Anatolian Fault Zone as well as the Menderes Massif. A large waveform data set was compiled from a variety of local and regional seismic networks including 121 digital seismic stations (broad-band and short period) between 1999 and 2008 spanning the IA, the Anatolian Plateau and Azerbaijan. The data set was used to determine the nature of Lg wave propagation and characterize the nature of seismic attenuation within the crust of these regions. Lg waveforms were used to calculate the frequency-dependent Lg-Q _o and Lg- $ \eta $ . A wide range of Lg-Q _o values was obtained between ~52 ± 6 and 524 ± 227. Low Lg-Q _o values (~90–155) are calculated towards the north of IA, Iskenderun Gulf and its vicinity, Bingöl-Karl?ova, Izmit and its vicinity. Lg-Q _o values are especially low (<90) along the Menderes Massif and the Aksehir-Simav Fault Zones. This may be due to intrinsic attenuation of Lg associated with the partially molten crust and young volcanism. The high Lg-Q _o values (~350) are probably caused by the crust not being subject to large amounts of extensional deformation like the Antalya Gulf and apparently being thick enough to support Lg propagation. Relatively higher values along the border of this subduction zone and plate boundary might be related to the Taurus Mountain belts and Bitlis-Zagros Suture Zone. The lateral frequency dependency Lg- $ \eta $ is also consistent with high tectonic activity in this region.     

三维椭球形盆地对球面波的散射

针对三维沉积盆地对球面波的散射问题,发展一种快速宽频间接边界元方法(IBEM).利用ANSYS建立求解模型,基于Intel-Fortran编译器编译相应的计算程序,对基岩半空间三维半椭球形盆地对球面波的散射进行了数值分析,着重探讨入射波频率、波源埋深、波源与不规则地形(沉积盆地)距离等参数对地震动特性的影响规律,计算方...     

Efficiency of Lg Propagation from SmS Dynamic Ray Tracing in Three-dimensionally Varying Crustal Waveguides

The effects of variations in crustal and basin thicknesses on the efficiency of Lg propagation in Central Asia are investigated by comparing the predictions of ray theoretical modeling of Lg as multiple, critically reflected, SmS rays with observations at the KNET, CNET, and ILPA arrays. Ray modeling predicts that strong spatial gradients in crustal thickness within 10 to 50 km of either the source or receiver can affect Lg efficiency over paths exceeding 500 km. Variations in the efficiency of Lg propagation to KNET, including instances of strong variations across the aperture of the array, are consistent with predicted blockage of SmS by high slopes of Moho topography associated with the Hindu Kush and Pamir mountains. At CNET and ILPA the blocking effects of Moho slopes associated with the Caucausus and Zagros mountain belts and Black and Caspian Seas are accentuated by the effects of thick basins. Thick basins (greater than 2000 m) are found to have large effects on the efficiency of Lg propagation from shallow (1 km) crustal sources, blocking either a portion or all of the Lg phase along paths crossing these basins. The basin effects are strongly sensitive to source depth, nearly vanishing for sources 15 km and deeper.     

Simulation of Medium Wave Propagation in the Magnetosphere

Geomagnetism and Aeronomy - A model of the medium of the ionosphere and magnetosphere, including the distributions of the concentrations and temperatures, collision frequencies, and magnetic field...     

SH Wave Propagation in Viscoelastic Media

When comparing solutions for the propagation of SH waves in plane parallel layered elastic and viscoelastic (anelastic) media, one of the first things that becomes apparent is that in the elastic case the location of the saddle points required to obtain a high frequency approximation are located on the real p axis. This is true of the branch points also. In a viscoelastic medium this is not typical. The saddle point corresponding to an arrival lies in the first quadrant of the complex p-plane as do the branch points. Additionally, in the elastic case the saddle point and branch points lie on a straight line drawn through the origin (the positive real axis in the complex p-plane), while in the viscoelastic case this is generally not the case and the saddle point and branch points lie in such a manner as to indicate the degree of their complex values.In this paper simple SH reflected and transmitted particle displacement arrivals due to a point torque source at the surface in a viscoelastic medium composed of a layer over a half space will be considered. The path of steepest descent defining the saddle point in the first quadrant will be parameterized in terms of a real variable and the high frequency solutions and intermediate analytic results obtained will be used to formulate more specific constraints and observations regarding saddle point location relative to branch point locations in the complex p-plane.As saddle point determination for an arrival is, in general, the solution of a non-linear equation in two unknowns (the real and imaginary parts of the complex saddle point p ₀), which must be solved numerically, the use of analytical methods for investigating this problem type is somewhat limited.Numerical experimentation using well documented solution methods, such as Newton's method, was undertaken and some observations were made. Although fairly basic, they did provide for the design of algorithms for the computation of synthetic traces that displayed more efficient convergence and accuracy than those previously employed. This was the primary motivation for this work and the results from the SH problem may be used with minimal modifications to address the more complicated subject of coupled P-SV wave propagation in viscoelastic media.Another reason for revisiting a problem that has received some attention in the literature was to approach it in a fairly comprehensive manner so that a number of specific observations may be made regarding the location of the saddle point in the complex p-plane and to incorporate these into computer software. These have been found to result in more efficient algorithms for the SH wave propagation and a significant enhancement of the comparable software in the P-SV problem.     

Sedimentary Basins: Tectonic Recorders of Sediment Discharge from Drainage Catchments

The sedimentary basin is to geology what the drainage system is to geomorphology, with sediment flux being a unifying link between the two features and disciplines. The global sediment discharge from contents to oceans contributes to major mass recycling associated with plate subduction. Major plate reorganization is accompanied by small fluctuations in continental volumes due to recycling of sediment stored along the sink-like passive margins. In continental sedimentary basins the SEdiment Preservation RAtio (SEPRA) depends upon the tectonic ‘drawdown’ of sediment flux preserved in the subsurface until basin inversion. Major unsteadiness in sediment discharge arises from teconic or climatic causes, the latter prominent in the Quaternary to Holocene record. Contrasting trends in sediment discharge reflect both the nature of climate and vegetation change and source rock lithology. © 1997 by John Wiley & Sons, Ltd.     

A Numerical Experiment on the Horizontal to Vertical Spectral Ratio in Flat Sedimentary Basins

—?In this paper we study the seismic response of flat sedimentary basins and carry out numerical experiments to determine the extent to which we could go using the Horizontal to Vertical Spectral Ratio (HVSR) for a given site. The HVSR has been used by many researchers to characterize local conditions in terms of the dynamic response of the soil, and one of its variants, that proposed by NAKAMURA (1989) in which records of microtremors are used, is one of the most applied in recent years. We study the response of different configurations under incident waves coming from an explosive source using the Indirect Boundary Element Method (IBEM), and we investigate two cases: low- and high-velocity contrast. We compute the seismic response using the HVSR technique at various locations in the free surface of the basins, and compare it with the response calculated with the horizontal Sediment to Bedrock Spectral Ratio (SBSR) and with the Horizontal Component (HC) of the transfer function for the displacement at the same locations. The comparison shows that, in general, HVSR cannot provide the predominant period of a site due to the fact that this technique cannot predict accurately the Spectral amplification levels. On the other hand, the HVSR provides an erroneous response in the sedimentary basins which have a low-impedance contrast, with respect to bedrock, and with shape ratios like the one studied here, whereas it can reasonably well predict the fundamental local frequency when there is a high-impedance contrast, except in the center of the basin.     

An inversion of site response and Lg attenuation using Lg waveform

Based on spectral ratio method, a joint inversion method was used to obtain parameters of Lg wave attenuation and site response. The inversion method allows simple and direct (two-parameter) determination of Lg wave attenua- tion and site response from sparse spectral data, which are not affected by radiation pattern factor and different response of same instrument after geometrical spreading. The method was used successfully for estimating site re- sponse of stations of Zhejiang Seismic Network and measuring Lg wave attenuation. The study is based on 20 earth- quakes occurred in northeast of Taiwan with magnitude MS5.0~6.7 and 960 seismic wave records from 16 stations in Zhejiang area from 2002 to 2005. The parameters of site response and Lg attenuation were calculated with a fre- quency interval of 0.2 Hz in the range of 0.5 Hz to 10 Hz. Lg wave attenuation coefficient corresponding to U-D, E-W and N-S components are γ ( f )=0.00175 f 0.43485, γ ( f )=0.00145f 0.48467 and γ ( f )=0.0021f 0.41241, respectively. It is found that the site response is component-independent. It is also found that the site response of QIY station is significant above the frequency of 1.5 Hz, and that the site response of NIB station is low for most frequency     

应力诱导各向异性介质中地震波的传播

主要讨论了应力变化如何影响各向异性介质中波速度的问题。推导了一般各向异性介质在初始应力下的Christoffel方程,得到介质中3种波的相速度和初始应力的关系表达式;通过实验数据验证了单轴应力能够诱导各向异性,当施加单轴应力时,速度在沿应力的方向增加最大,在垂直应力的方向增加最小,实验结果与理论推导一致;用Christoffel方程的数值解模拟在3种对称情况下的弹性各向异性介质中初始应力对波速度的影响。数值结果表明:初始应力对各向异性介质中波传播速度的影响,随着各向异性强度的增加而增大,而且速度越慢,影响越大。     

Mapping Crustal Heterogeneity Using Lg Propagation Efficiency Throughout the Middle East, Mediterranean, Southern Europe and Northern Africa

—?In this paper we describe a technique for mapping the lateral variation of Lg characteristics such as Lg blockage, efficient Lg propagation, and regions of very high attenuation in the Middle East, North Africa, Europe and the Mediterranean regions. Lg is used in a variety of seismological applications from magnitude estimation to identification of nuclear explosions for monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty (CTBT). These applications can give significantly biased results if the Lg phase is reduced or blocked by discontinuous structure or thin crust. Mapping these structures using quantitative techniques for determining Lg amplitude attenuation can break down when the phase is below background noise. In such cases Lg blockage and inefficient propagation zones are often mapped out by hand. With our approach, we attempt to visually simplify this information by imaging crustal structure anomalies that significantly diminish the amplitude of Lg. The visualization of such anomalies is achieved by defining a grid of cells that covers the entire region of interest. We trace Lg rays for each event/station pair, which is simply the great circle path, and attribute to each cell a value equal to the maximum value of the Lg/P-coda amplitude ratio for all paths traversing that particular cell. The resulting map, from this empirical approach, is easily interpreted in terms of crustal structure and can successfully image small blockage features often missed by analysis of raypaths alone. This map can then be used to screen out events with blocked Lg prior to performing Q tomography, and to avoid using Lg-based methods of event identification for the CTBT in regions where they cannot work.¶For this study we applied our technique to one of the most tectonically complex regions on the earth. Nearly 9000 earthquake/station raypaths, traversing the vast region comprised of the Middle East, Mediterranean, Southern Europe and Northern Africa, have been analyzed. We measured the amplitude of Lg relative to the P-coda and mapped the lateral variation of Lg propagation efficiency. With the relatively dense coverage provided by the numerous crossing paths we are able to map out the pattern of crustal heterogeneity that gives rise to the observed character of Lg propagation. We observe that the propagation characteristics of Lg within the region of interest are very complicated but are readily correlated with the different tectonic environments within the region. For example, clear strong Lg arrivals are observed for paths crossing the stable continental interiors of Northern Africa and the Arabian Shield. In contrast, weakened to absent Lg is observed for paths crossing much of the Middle East, and Lg is absent for paths traversing the Mediterranean. Regions that block Lg transmission within the Middle East are very localized and include the Caspian Sea, the Iranian Plateau and the Red Sea. Resolution is variable throughout the region and strongly depends on the distribution of seismicity and recording stations. Lg propagation is best resolved within the Middle East where regions of crustal heterogeneity on the order of 100?km are imaged (e.g., South Caspian Sea and Red Sea). Crustal heterogeneity is resolvable but is poorest in seismically quiescent Northern Africa.     

Lg magnitude scale in the mainland of China

The data of short-period seismograms had been collected widely in the mainland area of China not including Xinjiang and Tibet. The physical quantities of Lg wave are determined respectively in the five subregions. The group velocities of priminary arrival and maximum amplitude of Lg wave are equal to 3.54±0.02 km/s and 3.30±0.05 km/s, respectively. The periods of Lg waves are between 0.2s to 1.2s, averaging 0.7s. The γ-values of Lg waves in the five subregions are equal to 0.0034±0.0001 km⁻¹ for East, 0.0031±0.0004 km⁻¹ for Southwest, 0.0027±0.0004 km⁻¹ for Northeast, 0.0022±0.0001 km⁻¹ for South, and 0.0021±0.0002 km⁻¹ for Northwest subreqion, respectively. The average γ-value for the five subregions, γ=0.0027±0.0006 km⁻¹. The relations among the amplitude ratioH/Z, the station correctionD _z andD _h of amplitudes, and among them and station site condition are discussed. The subregional magnitude calibration functions ofm _Lg had been established according to each regional γ-value. From these, the unified magnitude calibration function of Chinese mainland not including Xinjiang and Tibet was given by

Wave Propagation Modeling and Amplitude-Variation-with-Offset Response in a Fractured Coalbed

Coalbed methane can be detected employing the amplitudevariation-with-offset technique. However, there are two issues in applying this technique to a coalbed: strong azimuthal anisotropy resulting from high-density fractures, and the seismic response being composed of many or several individual reflections within the coalbed. To overcome these difficulties, we present an exact solution for reflections in extensive dilatancy anisotropy media. First, we build a three-layer model and simulate the wave propagation in this model. Then we derive an exact P- and converted S-wave reflection coefficient equation based on boundary conditions. Finally, substituting given model parameters into the exact equation, we obtain the variation in the reflection coefficient with incidence angle. The results show that the fracture factors, wavelet frequency and thickness of the coalbed have different effects on the reflection coefficient. Furthermore, we create a synthetic seismogram by forward calculation, and the result fits well with results of the exact equation.     

中国大陆Lg波衰减成像

本文利用2009—2019年间中国大陆固定台站记录的宽频带地震波形数据, 采用反双台法, 获得了约36万个台站对的Lg波Q值, 并反演得到了中国大陆Lg波的衰减模型.结果显示, 中国大陆Lg波衰减的Q值结构存在显著的横向非均一性.青藏高原、渤海湾盆地以及祁连山造山带等构造活动块体展现出明显的地壳低Q值强衰减特征, 而华南块体、鄂尔多斯盆地以及塔里木盆地等相对稳定的地区则呈现地壳高Q值弱衰减特征.Lg波衰减与构造活动强度、松散沉积物、地下流体或部分熔融以及地质块体结构等因素密切相关.渤海湾盆地的强衰减主要由上地壳巨厚的松散沉积物和较高的地壳温度导致.青藏高原和祁连山造山带的地壳强衰减主要归因于现今活跃的构造活动、复杂的断裂系统和地质结构以及富含地下流体甚至是部分熔融等因素.华南造山带和大兴安岭地区的地壳弱衰减特征可能与侵入岩冷却形成高强度的结晶基底以及古老造山带地壳构造趋于稳定相关.

     

爆炸波在地铁车站中的传播规律研究

应用LS-DYNA软件建立了北京比较典型的2层地铁车站有限元模型,通过数值模拟方法分析地铁车站发生内部爆炸时,爆炸波在地铁车站内部的传播规律。研究表明：①在0.3m/kg~（1/3）〈Z〈1.6m/kg~（1/3）范围内,车站内部超压峰值可用Henrych经验公式估算;②爆炸波在地铁车站的封闭空间不断出现反射、绕射现象,使得结构构件遭受多次冲击作用;③楼梯口起到了一定的泄爆作用。本文研究结果可供地铁车站的抗爆设计与安全分析参考。     

SH波在表面多层介质中传播的精确模拟

针对地震横波在地表低速层内的振幅放大效应问题，提出了一种模拟SH波在地表层状介质中传播的递推算法，并用它模拟了新西兰Alfredton盆地A10场址的SH波地震动响应特性。这个方法适用于具线性吸收性质的粘弹性介质。由于方法不受介质层厚薄制约，层厚可以无限薄化，实践上可以用许多薄层逼近的办法来模拟纵向上任意变化的连续介质。通过求取不同频率不同波数平面简谐波解并按实际问题的加权迭加可求解具特定波形和传播方向组合的任意SH波场。此方法在计算上具有解析解特有的精确性，稳定性和方便性，特别适用于模拟薄层介层，次波长现象及需要进行大量而又精确模拟计算的情形。     

川滇地区Lg波Q值层析成像研究

本研究联合使用川滇地区36个台站(云南区域数字台网中22个台站,四川区域数字地震台网中14个台站)记录的644个帆3.5～6.0地震共7468条波形数据,开展川滇地区Lg波Q值层析成像研究.研究中以2.8～3.6km/s的传播速度窗选取Lg波波段,采用平移窗方法计算Lg波垂直向观测位移谱.我们首先利用Lg波观测谱反演得到了研究区的平均QLg(f),通过棋盘测试方法确定将研究区划分为0.5°×0.5°的网格,用平均值QLg-1作为每个网格的初始值,分别反演了0.5Hz,1Hz,2Hz,5Hz等多个频率的QLg(f)分布.