Receiver Functions for Three-layer Media

We have extended the H-k stacking method of receiver functions applicable to a three-layer model, which is useful in studying detailed crustal structures. We have demonstrated its application with two sample sites in Taiwan, making use of travel times of converted phases from the direct P waves of teleseismic events as the P reflects and refracts at different discontinuities in the crust. This three-layer extension allows a closer examination on the crust, as well as the relationship between crustal layers and the associated V_p/V_s ratios. Data were processed using the multiple-taper correlation technique to obtain the radial receiver functions (RRFs). The relative time delays of the converted phases measured from the RRFs were used to estimate the depths of the crustal discontinuities. Results not only yield the depths to principal crustal layers, including the Moho, but also give the corresponding ranges of V_p/V_s ratio which are comparable with findings from other tomographic studies.     

Receiver Functions from Autoregressive Deconvolution

Summary Receiver functions can be estimated by minimizing the square errors of Wiener filter in time-domain or spectrum division in frequency domain. To avoid the direct calculation of auto-correlation and cross-correlation coefficients in Toeplitz equation or of auto-spectrum and cross-spectrum in spectrum division equation as well as empirically choosing a damping parameter, autoregressive deconvolution is presented to isolate receiver function from three-component teleseismic P waveforms. The vertical component of teleseismic P waveform is modeled by an autoregressive model, which can be forward and backward, predicted respectively. The optimum length of the autoregressive model is determined by the Akaike criterion. By minimizing the square errors of forward and backward predicting filters, autoregressive filter coefficients can be recursively solved, and receiver function is also estimated in the similar procedure. Both synthetic and real data tests show that autoregressive deconvolution is an effective method to isolate receiver function from teleseismic P waveforms in time-domain.     

平面剪切断层系相互作用的研究

本文利用数值模拟方法对 2种典型地震断层系进行了研究。首次引入利用所有裂纹面上裂纹张开位移 ( COD)数据计算应力强度因子和应力场的方法 ,提高了计算精度。本文对此方法进行了验证 ,利用边界元算法获得 COD数据 ,将全裂纹面上的 COD数据代入公式 ,拟合出应力强度因子和应力场。此方法后处理少、可靠性好、计算时间也相应较短 ,且易于操作实现 ,适用范围较广。在上述工作的基础上 ,本文利用此算法讨论了 2种典型的平面双等长裂纹系问题。第一种是平面共线剪切双裂纹 ,主要讨论了这种非单一裂纹在外力作用下的相互作用 ,定量给出了裂…     

Imaging the Moho and Subducted Oceanic Crust at the Isthmus of Tehuantepec, Mexico, from Receiver Functions

Using teleseismic data recorded along a transect, which we call VEOX (for Veracruz-Oaxaca seismic line), of 46 broadband stations installed across the Isthmus of Tehuantepec in southern Mexico, we obtained receiver functions and stacked them to study the Moho topography and back projected them to visualize the subducted slab geometry beneath the isthmus. We observed a back-azimuth dependent Moho thickness across the transect, particularly beneath the Los Tuxtlas Volcanic Field. Also, we observed the Cocos plate which subducts with an angle of 26° between 140 and 310?km from the trench. Comparison with regional seismicity indicates that it occurs below the oceanic crust.     

The Fluid Dynamics of Solid Mechanical Shear Zones

Shear zones in outcrops and core drillings on active faults commonly reveal two scales of localization, with centimeter to tens of meters thick deformation zones embedding much narrower zones of mm-scale to cm-scale. The narrow zones are often attributed to some form of fast instability such as earthquakes or slow slip events. Surprisingly, the double localisation phenomenon seem to be independent of the mode of failure, as it is observed in brittle cataclastic fault zones as well as ductile mylonitic shear zones. In both, a very thin layer of chemically altered, ultra fine grained ultracataclasite or ultramylonite is noted. We present an extension to the classical solid mechanical theory where both length scales emerge as part of the same evolutionary process of shearing the host rock. We highlight the important role of any type of solid-fluid phase transitions that govern the second degree localisation process in the core of the shear zone. In both brittle and ductile shear zones, chemistry stops the localisation process caused by a multiphysics feedback loop leading to an unstable slip. The microstructural evolutionary processes govern the time-scale of the transition between slow background shear and fast, intermittent instabilities in the fault zone core. The fast cataclastic fragmentation processes are limiting the rates of forming the ultracataclasites in the brittle domain, while the slow dynamic recrystallisation prolongs the transition to ultramylonites into a slow slip instability in the ductile realm.     

远震接收函数的各向异性研究

用接收函数研究地壳各向异性,首先要提取高质量的接收函数,然后计算接收函数Ps震相的分裂或拟合Ps震相随后方位角的变化特性来提取地壳的各向异性参数。在提取接收函数的方法中,频率域方法简单,执行效率高;时间域方法精度高、稳定性好。与先前的方法相比,接收函数的各向异性研究针对整个地壳,对于更准确地分析地壳结构、了解壳幔耦合关系、岩石圈形变及动力学特征具有重要的意义。作为研究各向异性的方法之一,接收函数的各向异性研究仍然在发展之中,会有更广阔的发展前景。     

利用接收函数方法研究四川地区地壳结构

采用接收函数反演和共转换点(CCP)偏移叠加成像方法, 利用四川数字地震台网宽频带的52个区域固定地震台站和布设的两条52个宽频带流动地震观测台站的远震地震波形数据资料, 对四川地区地壳结构进行研究。 结果表明, 四川地区的Moho面深度在青藏高原和四川盆地差异明显, 在川西高原地区地壳厚度为52~68 km, 在川滇地块地壳厚度为50~60 km, 在中地壳内存在不连续的低速层分布; 而在四川盆地地壳厚度为38~45 km, 地壳内没有低速层存在。 Moho面深度从川西高原的60多公里至四川盆地的约40 km, 在二者的交界处龙门山断裂带下面, 存在厚度约30 km左右宽的下降过渡带, 说明其下的Moho面可能受断层影响, 结构比较复杂; 在高原地区的上地壳界面和下地壳上界面比四川盆地的相应界面深; 高原地区在中地壳的上部有不连续的低速层分布, 在松潘—甘孜地块的上地壳下部存在向南东运动的脆性推覆体, 在羌塘—理塘地块的上地壳下部存在向南东和南运动的脆性物质流动。     

Velocity Model of Western Volgo-Uralia from Receiver Functions

The results of studying the deep structure of the Earth’s crust and upper mantle in the central part of the Russian platform from receiver functions are presented. The records of teleseismic waves by the Monakovo small-aperture seismic array in the region of the northwestern slope of the Tokmovskii Arch of the Volga–Kama anteclise are used. The modification of the P-receiver function method (Vinnik, 1977) suggested in (Sanina et al., 2014) for analyzing the receiver functions in the regions with a complexly structured upper part of the section and the presence of a thick sedimentary cover is applied. The method is based on separating the high- and low-frequency components of the seismic record and successive reconstruction of the V-s velocity section in the upper part of the crust, which is performed first and, next, the entire deep section of the crust and the mantle down to a depth of ~300 km. The positions of the seismic conversion boundaries in the crust and upper mantle beneath the Monakovo array are determined. The upper mantle velocity section constructed based on the observations at the Mikhnevo array (Sanina et al., 2014) is compared with the world data on the ancient Precambrian platform.     

接收函数的H-k方法及其应用

接收函数的H-k方法是研究地球内部结构的重要方法,由于接收函数不依赖于震源参数,可以更为精确地确定观测台站下方的地壳结构;同时H-k方法简单明了,可以利用深度-波速比扫描同时获取台站下方的地壳厚度及地壳内部的波速比或泊松比。本文利用H-k方法,对桂东北地区的宽频带地震流动台网观测数据的地震数据进行应用,结果表明该区的地壳厚度变化大（28~34 km）,而且地壳内部的波速比（1.801~1.897）或泊松比（0.277~0.308）较高,与桂东北地区的基性、超基性岩大片出露的事实相吻合。      

利用接收函数方法研究腾冲地区S波速度结构

腾冲地区邻近印度板块与欧亚板块碰撞、俯冲的边界,地质环境和构造背景十分复杂,是我国地震、火山活动比较活跃的地区之一.本文采用最大熵谱反褶积方法提取腾冲地区1.0°×0.8°范围内5个流动数字地震台站的宽频带远震接收函数,反演得到台站下方0～100km深度范围的S波速度结构,分析讨论了该地区的深部构造特征.结果表明:1)腾冲地区地质结构存在明显的横向非均匀性;2)盈江断裂两侧莫霍面深度有较大差异;3)腾冲和高黎贡山之间是地壳厚度和S波速度变化的高梯度带;4)盈江断裂东南、新生破裂带以西附近地区存在明显的低速层;5)盈江断裂和新生破裂带都可能对火山区的熔融体具有阻隔作用.     

Definition of Exclusion Zones Using Seismic Data

—?In verifying compliance with the Comprehensive Nuclear-Test-Ban Treaty (CTBT), there is a motivation to be effective, efficient and economical and to prevent abuse of the right to conduct an On-site Inspection (OSI) in the territory of a challenged State Party. In particular, it is in the interest of a State Party to avoid irrelevant search in specific areas. In this study we propose several techniques to determine `exclusion zones', which are defined as areas where an event could not have possibly occurred. All techniques are based on simple ideas of arrival time differences between seismic stations and thus are less prone to modeling errors compared to standard event location methods. The techniques proposed are: angular sector exclusion based on a tripartite micro array, half-space exclusion based on a station pair, and closed area exclusion based on circumferential networks.     

Rapid Growth of the Shear Strain in Weakened Zones of the Lithosphere

Izvestiya, Physics of the Solid Earth - A weakened zone in the lithosphere plunging into the mantle can lead to an earthquake after the application of a shear stress only in the case if the...     

Crustal Structure Along Sunda-Banda Arc Transition Zone from Teleseismic Receiver Functions

We analyzed receiver function of teleseismic events recorded at twelve Indonesian-GEOFON (IA-GE) broadband stations using nonlinear Neighbourhood Algorithm (NA) inversion and H-k stacking methods to estimate crustal thickness, V _p /V _s ratios and S-wave velocity structure along Sunda-Banda arc transition zone. We observed crustal thickness of 34–37 km in Timor Island, which is consistent with the previous works. The thick crust (> 30 km) is also found beneath Sumba and Flores Islands, which might be related to the arc-continent collision causing the thickened crust. In Timor and Sumba Islands, we observed high V _p /V _s ratio (> 1.84) with low velocity zone that might be associated with the presence of mafic and ultramafic materials and fluid filled fracture zone. The high V _p /V _s ratio observed at Sumbawa and Flores volcanic Islands might be an indication of partial melt related to the upwelling of hot asthenosphere material through the subducted slab.     

利用接收函数计算地壳各向异性的可靠性分析及倾斜界面的影响*

使用接收函数研究壳幔速度间断面和速度结构已是常用的技术,但介质各向异性或倾斜的莫霍(Moho)界面,都会造成接收函数波形的复杂性。本文利用远震P波接收函数计算台站下方的地壳各向异性,通过信噪比测试和谐波分析两种方法来验证各向异性结果的可靠性。通过甘肃地震台网的两个台站记录,讨论各向异性和倾斜界面对接收函数的影响,结果显示,台站BYT(白银)下方具有各向异性,而台站WYT(渭源)下方由于可能存在倾斜界面,得到的各向异性结果则有待进一步分析。为了更好地认识倾斜界面对各向异性计算结果的影响,采用合成理论地震图,计算接收函数,然后利用合成接收函数进行各向异性分析。结果表明,计算得到的快波方向不会受到倾斜界面的影响,但是时间延迟会受到影响。     

利用接收函数研究河北及邻区地壳厚度与泊松比分布特征

本文收集了河北及邻区174个宽频带地震仪2007年1月至2011年12月记录到的全球范围内体波震级6.0以上、共488个远震波形资料。利用H-k叠加搜索方法得到河北及邻区143个台站下方接收函数、地壳厚度和泊松比值分布特征。结果表明,河北及邻区地壳厚度具有明显的分块特征并与区域构造背景有较好的相关性,泊松比值在0.25上下变化,与全球的泊松比估计值相当。河北及邻区地壳厚度以太行山山脉为界呈现出东薄西厚的特征并与地表地形起伏呈正相关;山西地震带南北地壳厚度存在明显差异,大同、宁武及安泽盆地地区泊松比值超过0.3;张—渤地震带的地壳厚度自东向西逐渐增厚,泊松比值则表现为相对较强的横向变化。华北平原带的薄地壳和低泊松比值特征可能与华北克拉通拆沉作用相关联。     

利用接收函数研究程海断裂带莫霍面深度和地壳内S波速度结构特征

利用程海断裂带附近27个数字地震台站远震波形资料,提取每一个台站的接收函数,计算出各台站莫霍面深度同时利用时间域线性反演方法,获得了各个台站下方的横波速度。结果显示:程海断裂带莫霍面深度从南部42km增至北部的54km,南部和北部莫霍面深度有明显的不同。从程海断裂带下方不同深度S波速度剖面可以看出,宾川及其北东部地区中下地壳存在明显的低速层,此低速层可能与还没有固结的热物质有关。而永胜南部地区,地壳中S波速度垂直变化剧烈,低速异常高速异常交替丛生,这可能是此区地震频发的主要原因。同时,本文对宽频带地震仪和短周期地震仪得到的接收函数进行了初步的对比分析。     

Stress Correlation Function Evolution in Lattice Solid Elasto-dynamic Models of Shear and Fracture Zones and Earthquake Prediction

--It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.