Spatiotemporal variations in the statistical parameters of the shear wave splitting from local events beneath Eastern Hokkaido

Statistical analysis is carried out of the distribution of the fast Swave $(\overset{\lower0.5em\hbox{$(\overset{\lower0.5em\hbox{, polarization azimuths, time delays between the split shear waves $(\delta \overset{\lower0.5em\hbox{$(\delta \overset{\lower0.5em\hbox{ and the dominant frequencies of P waves $(\overset{\lower0.5em\hbox{$(\overset{\lower0.5em\hbox{ from weak seismic events that hit Eastern Hokkaido during the period of 1998–2003, which also includes strong events with M > 6 and the strong Tokachi-oki earthquake (M8.0) of September 26, 2003. The analysis revealed oscillatory pattern of variations in wave parameters along Eastern Hokkaido. In the rigid crustal blocks beneath the Tokachi and Kushiro Plains, stability is observed in the fast azimuth orientation and unstable behavior of $\delta \overset{\lower0.5em\hbox{$\delta \overset{\lower0.5em\hbox{ and $\overset{\lower0.5em\hbox{$\overset{\lower0.5em\hbox{ in the interval of periods 0.5–2.5 years. In the weakened zones beneath the Hidaka Range and the Nemuro Peninsula, periodical alternation of $\overset{\lower0.5em\hbox{$\overset{\lower0.5em\hbox{ along orthogonal directions is revealed. It is found that the parameters of the waves coming from the seismic events originating in different depth intervals experience anomalous changes before strong earthquakes. These anomalies reflect depth changes in the type of deformation of the medium. Active deformation is observed beneath the Hidaka Range. Dynamic changes in beneath the stations on either side of the range (ERM, MYR), their antiphase pattern, and the discrepancies in the behavior of the parameters from the events that originate in different depths indicate that these regions act as active buffers in the redistribution of stresses in the Kuril-Japan arc-arc junction zone. The onset of instability that develops before the Tokachi-oki earthquake is marked by gradual shortening of the period of oscillations in $\overset{\lower0.5em\hbox{$\overset{\lower0.5em\hbox{ and transition to the north and east azimuths as well as by the decreasing of $\overset{\lower0.5em\hbox{$\overset{\lower0.5em\hbox{ from the intermediate-depth events starting from middle 2002, and for deep events, starting from 2000. The acceleration of the process begins in the middle of 2003 and is recognized by a noticeable increase in and a decrease in $\overset{\lower0.5em\hbox{$\overset{\lower0.5em\hbox{ beneath almost all stations.     

Time variations in parameters of split S waves from weak local earthquakes beneath eastern Hokkaido

Parameters of seismic waves from clusters of local weak earthquakes that occurred at the upper boundary of the seismofocal zone in depth intervals of 40–60 and 70–90 km along eastern Hokkaido are investigated for the period 1998–2003, including the strong (M = 8.0) Tokachi-Oki earthquake of September 26, 2003. Analysis of data indicates that the distribution of anisotropic properties along Hokkaido is inhomogeneous and parameters of split waves (the azimuth of the fast S wave and the time delay between split S waves) are sensitive to variations in the stress-strain state of the medium. Unstable behavior of split wave parameters and increased values of the ratio V _P /V _S for clusters of events in the areas of the Hidaka Mountains and Nemuro Peninsula imply that the medium is in a mechanically weakened state (in the regime of intense dilatant deformation). On the contrary, the regions beneath the Tokachi and Kushiro plains are more rigid and are characterized by lower V _P /V _S values and comparatively stable behavior of wave parameters. Anomalous parameters of split waves from events of clusters in areas of different stations correlate with each other and are related to occurrence times of large earthquakes around Hokkaido, which may point to a redistribution of stresses and strains and fluid migration in the subduction zone.     

The distribution of split shear-wave parameters in the backarc region of the South Kuril subduction zone

This paper reports a study in the distribution of parameters of split shear waves excited by deep earthquakes in the Sakhalin and Hokkaido area, with the orientation of the axes of symmetry in the earth being estimated assuming a viscoelastic anisotropic model for the mantle. A mantle flow along NW 310° ± 20° has been identified beneath the Japan Sea. The flow is inclined at an angle of 20–30° relative to the horizontal plane and is consistent with the motion of the Pacific plate. Beneath the southern Sea of Okhotsk the fast a-axis [100] of olivine is oriented NE 30 ± 15°, nearly parallel to the trend of the Kuril arc, while the c-axis [001] is inclined at ～35° relative to the horizontal plane. Bearing in mind the increased heat flow in the region, we assume the development of a mantle flow along the olivine c-axis at NW 300° ± 20° and with the shear plane [010] in the conditions of partial mantle melting (the B-type LPO). The lowest anisotropy (1–2.5%) was identified beneath Sakhalin and the greatest (3–5%) beneath the Japan Sea. An increasing degree of anisotropy is also noted to occur with an increasing depth of focus (down to 350 km).     

Anisotropy and temporal variations in the fast azimuth of the fast shear wave beneath southern Kamchatka during 1993–2002

Parameters of split shear waves from local earthquakes in the area of the PET IRIS station (town of Petropavlovsk-Kamchatski) were measured over the period 1993–2002 for the study of anisotropic properties of rocks in the subduction zone and variations in the fast azimuth of the fast shear wave (?). The dominating fast shear wave polarization directions were oriented in 1993–2002 along N90°E ± 20° in agreement with the direction of the Pacific plate motion. The normalized shear wave delay times δt _SS increase to a depth of 150 km. The values of δt _SS are largest (up to 20 ms/km) for earthquakes at depths of 50–60 and 90–150 km and smallest (up to 6 ms/km) for earthquakes at depths greater than 200 km. The fast azimuths for events with H < 80 km are described in terms of a horizontal transversely isotropic (HTI) model of the medium, with the axis oriented northward. Temporal variations in the fast azimuths with an amplitude of up to 90° and a predominant period of about 400–600 days are observed for events at depths of 80–120 km. The anisotropy of rocks is described by effective models of the orthorhombic and HTI symmetries. The predominant fast shear wave fast azimuths from events at depths of 120–310 km vary with time: the polarization axis was oriented to the north in 1993–1995, to the north and east in 1996–1998, to the east in 1999–2000, and to the northeast and southeast in 2001–2002. The anisotropy of rocks can be described in terms of the HTI model with the symmetry axis subparallel to the focal zone dip.     

The velocity structure of the upper mantle of europe from the ambient noise surface wave tomography

The method for surface wave tomography based on the records of ambient seismic noise (Ambient Noise Tomography, ANT) is applied to the data from the East European and West European stations. In order to reduce the effects of the earthquakes at long periods, the cross correlation functions were calculated for the time interval of 2001–2003, when distinct clusters of the earthquakes were absent. Using the local dispersion curves in the range of 10–100 s, we reconstructed the vertical velocity sections at the nodes of the 3° × 3° grid. On the basis on these curves, we calculated the horizontal distributions of S-velocity variations in the upper mantle in the depth interval of 75–275 km and the vertical velocity sections along the profiles across the Vrancea zone and the region of the Baltic and Ukrainian shields. The velocity distribution in the Vrancea zone confirms the subduction of the ancient oceanic plate from the east westwards and the detachment of its bottom part, as hypothesized by some authors. Beneath the Baltic Shield lithosphere, there is a low-velocity zone, which can be interpreted as the asthenospheric layer. It is noted that the velocity distributions beneath the Baltic and Ukrainian shields are similar, which probably points to the genetic relationship between these two structures.     

中国大陆及邻区SKS波分裂研究

SKS波分裂测量是研究大陆地幔的形变特征、探索大陆动力学和演化过程的重要工具. 本文对中国大陆及邻区地震台站的SKS波分裂现象进行了研究. 选用中国数字化宽频带地震台网(CB台网)和美国IRIS数据中心提供的三分量宽频带数字化地震资料，使用SC(Silver and Chan,1991)方法，得到了中国大陆及周边地区80多个台站下方上地幔各向异性参数，即快波偏振方向φ和快慢波到时差δt. 快波偏振方位在相同地块有一定的优势排列方向，大多数台站快波偏振方向都能与过去或现今大规模的构造运动得到很好的符合. 整个研究区域所得到的分裂延迟时间在0.4～2.4s之间，平均为1.2s. 根据SKS波测量得到的分裂参数，分析了该研究区域各向异性介质的特性，从而探索与岩石弹性各向异性相关的地球内部动力学过程.     

Collision orogeny at arc-arc junctions in the Japanese Islands

Abstract In the Japanese Islands, collision tectonics are operating at arc-arc junctions in three regions: Hokkaido, Central Japan and Kyushu. Hokkaido is situated at the junction of the Kuril and Northeast Japan Arcs. The Kuril fore arc sliver collides with the Northeast Japan Arc, and the lower crust of the Kuril Arc thrusts upon the fore arc of the Northeast Japan Arc in Hokkaido. Outcrops of the lower crust are observed in the Hidaka Mountains in the fore arc of the junction area. Central Japan is in the juncture area among the Northeast Japan, Izu-Bonin, and Southwest Japan Arcs. The Izu-Bonin arc is colliding against the Honshu mainland, which has been bent by the collision. Kyushu is a juvenile collision area between the Southwest Japan and Ryukyu Arcs. The fore arc of the Southwest Japan Arc is starting to underthrust beneath the Kyushu islands along the Bungo Strait, where shallow seismicity within the crust is active in terms of the collision. Collision tectonics are observed at most of the arc-arc junctions in the circum-Pacific orogenic belts and may be an important process contributing to the relatively rapid growth of new continental crust in subduction zones.     

长白山-镜泊湖火山区地壳结构接收函数研究

利用71个远震的波形资料,用接收函数方法提取了布设在长白山—镜泊湖火山区的34个宽频带流动数字地震台站的接收函数,通过对接收函数反演,获得了台站下方的S波速度结构.研究结果表明,沈阳—敦化一线莫霍面深度32～33km,向西地壳厚度加厚,到长春附近地壳厚度约为36km.在天池火山口莫霍面深度为达38km,而镜泊湖火山口森林的莫霍面深度约为39km.总体看研究区的地壳厚度是南浅北深.长白山天池火山口附近地下10km左右有一明显的低速层存在;镜泊湖火山口森林附近30km也可能有低速体存在;研究发现莫霍面上S波速度梯度在火山口附近和远离火山口有明显区别.在火山口附近其莫霍面的S波速度梯度比非火山口地区的S波速度梯度明显小,说明火山口下与一般的地壳莫霍面结构有差别.研究发现沈阳—敦化一线两侧的莫霍面深度有较大变化,其位置与地表的敦化—密山断裂基本一致,说明敦化—密山断裂是研究区的一条非常重要的地质构造带.     

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.     

Palaeomagnetism of the Early Carboniferous Deer Lake Group,western Newfoundland: no evidence for mid-Carboniferous displacement of “Acadia”

Early Carboniferous (Viséan to possibly earliest Namurian) sedimentary rocks of the Deer Lake Group of western Newfoundland rest unconformably on Grenvillian basement rocks of the Canadian Shield which form the western border of the Early Palaeozoic Appalachian orogen. In addition to magnetically soft magnetizations directed along the present field, three families of magnetization directions are found. Two of them (referred to as N (north) and S (south)) are antiparallel and prefolding, and were probably acquired during the depositional process. N and S are roughly of equal frequency. They have a mean direction irrespective of sign of 0.7°, ?35.2°, k = 40, α₉₅ = 8.9°, and a palaeopole 21.5°N, 121.8°E (10.3°, 6.0°) corresponding to a palaeolatitude of 20 ± 6°S. This agrees with the palaeolatitude (17 ± 5°S) determined from the somewhat older Early Carboniferous (Tournaisian) Terenceville Formation of the Avalon Platform on the eastern side of the Appalachian orogen in Newfoundland. The third magnetization, referred to as H (roughly horizontal), has a mean direction 156.8°, ?13.3°, k = 37, α₉₅ = 10.1°, and a palaeopole 45.4°N, 140.3°E (10.3°, 5.3°) corresponding to a palaeolatitude of 7 ± 4°N; we interpret this to be an early Kiaman (latest Carboniferous to early Permian) overprint probably acquired chemically. The palaeolatitude determined from the H overprint agrees with that determined from Early Carboniferous rocks of cratonic North America west of the Appalachians. Therefore, we argue, Early Carboniferous palaeofield determinations for cratonic North America have been strongly biased by unremoved Kiaman overprints. Because of this, and because of the good agreement between Early Carboniferous palaeolatitudes obtained from opposite margins of the Appalachian orogen, we suggest that there is, at present, no palaeomagnetic evidence for the previously proposed 1500 km displacement from the south of an eastern portion of the Appalachians (“Acadia”) relative to cratonic North America during the Carboniferous.     

Space-time variations in the parameters of split S waves excited by local earthquakes in southern Kamchatka

An analysis of the distribution of split S-wave parameters due to small local earthquakes occurring down to a depth of 200 km beneath the PET IRIS station (in the town of Petropavlovsk-Kamchatskii) for the period 1993–2002 revealed an inhomogeneous distribution of anisotropic properties and time-dependent changes in the state of stress and strain in the sinking slab. Three time intervals have been identified from the behavior of the split-wave parameters: 1993–1995, 1996–1998, and 1999–2002. The most orderly orientation of fast-wave polarization azimuths at all depths and lower stresses in the earth was observed in 1996–1998. Regions of relatively stable and unstable parameter behavior consistent with increased and decreased P and S velocities have been identified. The central block down to a depth of 120 km and the lower part of the Benioff zone are classified as stable regions dominated by E-SE fast-wave azimuths. Unstable behavior of fast-wave azimuths was noticed for depths of 60 to 90 km and around the top of the Benioff zone. The highest values of S-wave delay time and regions where the delay times were migrating mostly occur at contacts between more rigid and weakened blocks.     

Evolution history of the Hidaka-oki (offshore Hidaka) basin in the southern central Hokkaido, as revealed by seismic interpretation, and related tectonic events in an adjacent collision zone

Off the southern coast of Hokkaido the Hidaka-oki (offshore Hidaka) basin has developed on the western flank of a collision suture under the influence of long-standing compressional plate motion and provoked tectonic stresses around the northwestern Pacific rim throughout the late Cenozoic. The basin forming history of the Japan arc and Kuril arc collision zone is described on the basis of seismic reflection data interpretation. We identify two stages of basin formation: the older (late Oligocene-Miocene) faulted en echelon graben (pull-apart basin) and younger (Plio-Pleistocene) regional downwarping. Paleoenvironmental changes recorded within the fore-arc sediments indicate that the older basin filled up by the late Miocene. We inferred the volumes of the distinctive basins from the depth-conversion of seismic data, which suggest episodic uplifts and massive erosion of the Hidaka Mountains in the middle-late Miocene and the Plio-Pleistocene. Estimated sediment supply rates into the basins have a similar level for the both stages. Cause of an episodic uplift in the older stage is attributed to the delayed opening of the Japan Sea. The eastern Eurasian margin underwent N-S right-lateral faulting at 25 Ma as a result of rifting of the Kuril back-arc basin. Formation of the Japan Sea back-arc basin since the early Miocene (ca. 20 Ma) caused eastward motion of the western Hokkaido block and transpressive regime along the pre-existing N-S shear deformation zone.     

S-wave velocity structure in Tangshan earthquake region and its adjacent areas from joint inversion of receiver functions and surface wave dispersion*

Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas (39°N–41.5°N, 115.5°E–119.5°E), we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion. We obtained some detailed information about the Tangshan earthquake region and its adjacent areas, including sedimentary thickness, Moho depth, and crustal and upper mantle S-wave velocity. Meanwhile, we also obtained the v_P/v_S structure along two sections across the Tangshan region. The results show that: (1) the Moho depth ranges from 30 km to 38 km, and it becomes shallower from Yanshan uplift area to North China basin; (2) the thickness of sedimentary layer ranges from 0 km to 3 km, and it thickens from Yanshan uplift region to North China basin; (3) the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure, while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust. Compared with the upper crust, the heterogeneity of the middle and lower crust is more obvious; (4) the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust, and the low v_S and high v_P/v_S beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.     

2012年江苏高邮、宝应交界MS4.9地震发震断层参数测定

基于一维单侧有限移动震源模式,根据地震波传播过程中的多普勒效应,分别利用P波和S波拐角频率的方位变化,反演2012年7月20日江苏高邮、宝应交界MS4.9地震的发震断层面参数。P波和S波拐角频率的反演结果一致显示:本次地震的断层面破裂方向为232°左右,破裂面呈NE-SW向;地震马赫数v/c为0.2左右,平均破裂速度小于S波速度,破裂长度较短,为0.2~0.3km左右。破裂面方位与震源机制解、宏观烈度调查和余震精定位的研究结果具有一致性,结合震区周边的地质构造背景,分析认为滁河断裂很可能是高邮、宝应交界MS4.9地震的发震构造。     

青藏高原中部INDEPTH-Ⅲ剖面上地幔远震S波层析成像

印度板块向欧亚俯冲前缘位于班公—怒江缝合带附近,但是印度岩石圈地幔的俯冲形态和形变过程仍然缺乏共识,在不同地区使用不同方法获得的结果之间存在明显差异.本文使用青藏高原中部INDEPTH-Ⅲ剖面远震S波波形数据,提取走时信息,通过层析成像方法获得剖面下方S波速度扰动图像.结果显示:在班公—怒江缝合带下方100至300km深度范围内存在一个高角度(约65°)北倾的S波高速体,推测可能是回退的印度岩石圈板片或/和小规模对流引起的岩石圈拆沉后残留的印度大陆岩石圈板片.     

Crustal thickness,VP/VS ratio,and shear wave velocity structures beneath Myanmar and their tectonic implications

The crustal structure in Myanmar can provide valuable information for the eastern margin of the ongoing Indo-Eurasian collision system. We successively performed H–k stacking of the receiver function and joint inversion of the receiver function and surface wave dispersion to invert the crustal thickness (H), shear wave velocity (V_S), and the V_P/V_S ratio (k) beneath nine permanent seismic stations in Myanmar. H was found to increase from 26 ?km in the south and east of the study area to 51 ?km in the north and west, and the V_P/V_S ratio was complex and high. Striking differences in the crust were observed for different tectonic areas. In the Indo-Burma Range, the thick crust (H ?～ ?51 ?km) and lower velocities may be related to the accretionary wedge from the Indian Plate. In the Central Myanmar Basin, the thin crust (H ?= ?26.9–35.5 ?km) and complex V_P/V_S ratio and V_S suggest extensional tectonics. In the Eastern Shan Plateau, the relatively thick crust and normal V_P/V_S ratio are consistent with its location along the western edge of the rigid Sunda Block.     

蒙古中南部地区地壳各向异性及其动力学意义

利用蒙古中南部地区布设的69套宽频带数字地震仪2011年8月—2013年7月记录的远震事件,使用时间域反褶积方法提取接收函数,并挑选高质量Pms震相,通过改进的剪切波分裂方法对研究区地壳各向异性参数进行了研究,最终获取了1473对各向异性参数.经过统计分析,有48个台站可以归纳出两个方向的各向异性,11台站得到单个方向的各向异性,而剩余10个台站各向异性方向比较发散.结果显示,各向异性在蒙古中南部地壳中呈不均匀分布,有54个台站得到了NE-SW向各向异性,快波偏振方向平均值为N58°E±16°,与最大水平主应力σHmax方向和区域内主要断层走向一致,说明这部分地壳各向异性的主要成因存在于上地壳,可能与流体填充的微裂隙有关.而NW-SE向各向异性在53个台站被观测到,各向异性方向变化范围平均N132°E±16°,与研究区大部分SKS分裂快波方向具有较好的一致性,说明下地壳成岩矿物晶体定向排列是各向异性的主要成因.研究区地壳各向异性的分层特征总体上支持岩石圈受到NE-SW向挤压的动力学模型.     

湖北地区上地幔各向异性及其动力学意义

SKS波分裂测量是研究大陆地幔的形变特征、探索大陆动力学和演化过程的重要工具.本文选用湖北数字地震台网(HB台网)和中国数字化宽频带地震台网(CB台网)提供的三分量宽频带数字化地震资料,使用旋转相关法和切向能量最小法,计算得到了湖北地区20个地震台站下方上地幔各向异性参数.结果表明,快波偏振方向较为复杂:大别造山带内,...     

Image enhancement by multi‐parameter characterization of common image gathers

We present an innovative approach for seismic image enhancement using multi‐parameter angle‐domain characterization of common image gathers. A special subsurface angle‐domain imaging system is used to generate the multi‐parameter common image gathers in a summation‐free image space. The imaged data associated with each common image gathers depth point contain direction‐dependent opening‐angle image contributions from all the available incident and scattered wave‐pairs at this point. Each direction‐dependent opening‐angle data can be differently weighted according to its coherency measure. Once the optimal migration velocity is used, it is assumed that in the actual specular direction, the coherency measure (semblance) along reflection events, from all available opening angles and opening azimuths, is larger than that along non‐specular directions. The computed direction‐dependent semblance attribute is designed to operate as an imaging filter which enhances specular migration contributions and suppresses all others in the final migration image. The ability to analyse the structural properties of the image points by the multi‐parameter common image gather allows us to better handle cases of complicated wave propagation and to improve the image quality at poorly illuminated regions or near complex structures. The proposed method and some of its practical benefits are demonstrated through detailed analysis of synthetic and real data examples.     

Parallel PSM/FDM Hybrid Simulation of Ground Motions from the 1999 Chi-Chi,Taiwan, Earthquake

-- A new technique for the parallel computing of 3-D seismic wave propagation simulation is developed by hybridizing the Fourier pseudospectral method (PSM) and the finite-difference method (FDM). This PSM/FDM hybrid offers a good speed-up rate using a large number of processors. To show the feasibility of the hybrid, a numerical 3-D simulation of strong ground motion was conducted for the 1999 Chi-Chi, Taiwan earthquake (M_w 7.6). Comparisons between the simulation results and observed waveforms from a dense strong ground motion network in Taiwan clearly demonstrate that the variation of the subsurface structure and the complex fault slip distribution greatly affect the damage during the Chi-Chi earthquake. The directivity effect of the fault rupture produced large S-wave pulses along the direction of the rupture propagation. Slips in the shallow part of the fault generate significant surface waves in Coastal Plain along the western coast. A large velocity gradient in the upper crust can propagate seismic waves to longer distances with minimum attenuation. The S waves and surface waves were finally amplified further by the site effect of low-velocity sediments in basins, and caused the significant disasters.