利用相对反投影方法对2010年Mw 8.8级智利地震破裂过程成像

基于美国流动台网的三分量远震P波资料,利用相对反投影方法,对2010年2月27日Mw 8.8智利地震的破裂过程成像。根据美国地质调查局(USGS)给出的震源参数及源区地质资料,设定一个小倾角断层面。将远震P波反投影至该断层面上,获取震源破裂过程。结果表明,智利地震破裂由3个子事件组成,破裂长度至少为513 km,宽度至少为100 km。破裂从震中开始向北方和南方两个方向双侧扩展,但北侧的破裂明显强于南侧。北侧的破裂长度约为340 km,破裂持续时间约125 s,破裂速度约为2.87 km/s;而南侧破裂长度约为173 km,破裂持续时间约99 s,破裂速度约为1.84 km/s。另外,北侧两个子事件的能量释放高峰分别出现在破裂开始之后16 s和79 s,与初始破裂点的距离分别为75 km和230 km。南侧子事件的能量释放高峰出现时刻为80 s,与初始破裂点的距离为145 km。智利地震在破裂过程中可能受到了破裂区域的障碍体影响。     

Extent, duration and speed of the 2004 Sumatra-Andaman earthquake imaged by the Hi-Net array

The disastrous Sumatra-Andaman earthquake of 26 December 2004 was one of the largest ever recorded. The damage potential of such earthquakes depends on the extent and magnitude of fault slip. The first reliable moment magnitude estimate of 9.0 was obtained several hours after the Sumatra-Andaman earthquake, but more recent, longer-period, normal-mode analyses have indicated that it had a moment magnitude of 9.3, about 2.5 times larger. Here we introduce a method for directly imaging earthquake rupture that uses the first-arriving compressional wave and is potentially able to produce detailed images within 30 min of rupture initiation. We used the Hi-Net seismic array in Japan as an antenna to map the progression of slip by monitoring the direction of high-frequency radiation. We find that the rupture spread over the entire 1,300-km-long aftershock zone by propagating northward at roughly 2.8 km s(-1) for approximately 8 minutes. Comparisons with the aftershock areas of other great earthquakes indicate that the Sumatra-Andaman earthquake did indeed have a moment magnitude of approximately 9.3. Its rupture, in both duration and extent, is the longest ever recorded.     

Slip model for the 2011 M_w 9.0 Sendai （Japan） earthquake and its M_w 7.9 aftershock derived from GPS data

Using GPS-measured coseismic and post-seismic displacements for the 8 h following the M w 9.0 Sendai earthquake of March 11, 2011, coseismic and post-seismic fault slip models were developed based on a layered crustal model. The geodetic moment magnitude of the main shock was measured as approximately M w 8.98. The slip exhibits clear reverse characteristics, with a maximum near the hypocenter, and a magnitude of about 23.3 m. Some strike-slip behavior may occur on the two sides of the peak rupture zone. Almost 90% of the seismic moments released by the main shock occurred at depths less than 40 km. The energy released by the fault slip in the 8 h following the main shock is approximately equal to an earthquake of M w 8.13. With a maximum of ~1.5 m, the post-seismic slip was concentrated in the southwestern part of the coseismic rupture fault, which agrees well with the location and behavior of the M w 7.9 aftershock. This implies that the post-seismic deformation in the 8 h after the main shock was mainly induced by the M w 7.9 aftershock. In addition, a post-seismic slip of 0.2-0.4 m was observed at the down-dip extension of the coseismic rupture, which may have been caused by the effect of after-slip during this period.     

Tracking the rupture of the Mw = 9.3 Sumatra earthquake over 1,150 km at teleseismic distance

On 26 December 2004, a moment magnitude Mw = 9.3 earthquake occurred along Northern Sumatra, the Nicobar and Andaman islands, resulting in a devastating tsunami in the Indian Ocean region. The rapid and accurate estimation of the rupture length and direction of such tsunami-generating earthquakes is crucial for constraining both tsunami wave-height models as well as the seismic moment of the events. Compressional seismic waves generated at the hypocentre of the Sumatra earthquake arrived after about 12 min at the broadband seismic stations of the German Regional Seismic Network (GRSN), located approximately 9,000 km from the event. Here we present a modification of a standard array-seismological approach and show that it is possible to track the propagating rupture front of the Sumatra earthquake over a total rupture length of 1,150 km. We estimate the average rupture speed to be 2.3-2.7 km s(-1) and the total duration of rupture to be at least 430 s, and probably between 480 and 500 s.     

Static slip model of the M w 9.0 Tohoku (Japan) earthquake: Results from joint inversion of terrestrial GPS data and seafloor GPS/acoustic data

Based on co-seismic displacements recorded by terrestrial GPS stations and seafloor GPS/acoustic stations, the static slip model of the 2011 Mw 9.0 Tohoku earthquake was determined by inverting the data using a layered earth model. According to a priori information, the rupture surface was modeled with a geometry that is close to the actual rupture, in which the fault dip angle increases with depth and the fault strike varies with the trend of the trench. As shown by the results inferred from the joint inversion, the "geodetic" moment is 3.68 × 10 22 Nm, corresponding to Mw 9.01, and the maximum slip is positioned at a depth of 13.5 km with a slip magnitude of 45.8 m. Rupture asperities with slip exceeding 10 m are mainly distributed from 39.6 to 36.97°N, over a length of almost 240 km along the trench. The slip was mostly concentrated at depths shallower than 40 km, up-dip of the hypocenter. "Checkerboard" tests reveal that a joint inversion of multiple datasets can resolve the slip distribution better than an inversion with terrestrial GPS data only-especially when aiming to resolve slip at shallow depths. Thus, the joint inversion results obtained by this work may provide a more reliable slip model than the results of other studies that are only derived from terrestrial GPS data or seismic waveform data.     

Focal mechanism analysis and parameter estimation of Zhangbei-Shangyi earthquake from differential SAR interferometry

On 10 January, 1998 an earthquake of Ms=6.2occurred in the Zhangbei-Shangyi region of North China.The surface seismic deformation was measured in the previous study using the 3 pass ERS-1/2 SAR differential interferometric technology (D-INSAR). In this note the focal mechanism of Zhangbei-Shangyi earthquake is estimated from the D-INSAR measurement of surface deformation based on a standard elastic dislocation model for seismic displacement. The inversion procedure is an iterative, linear least-squares algorithm. Through the relation between the focal parameters and displacement in the line of sight direction measured in the radar interferogram, the optimum focal parameter set is derived. The results show that the seismic fault of Zhangbei-Shangyi earthquake is a thrust fault dipping SW with a large right-lateral displacement component.The strike and dip are 95° and 30° respectively on a fault patch of 12 km long by 14 km wide. Its hypocenter is located at N40°58', E114°21', and 7.5 km in depth. The estimated slip vector is 0.728 m with a rake of 105.95°, the trend of slip is NW13.26°, and M0is 2.69×1018 N @ m.     

汶川地震滑坡详细编录及其与全球其他地震滑坡事件对比

 汶川地震滑坡编录结果表明,2008年5月12日汶川地震触发了197481处滑坡,面积约为1160km²,大致分布在一个面积约1.1x10⁵km²的区域内,是当前汶川地震滑坡最全面详细的编录成果,也是有记录以来单次地震事件触发滑坡数量最多的记录。其中196007处滑坡(占滑坡总数量的99.25%),面积为1150.622km²(占总滑坡面积的99.19%),分布在一个面积为44031km²、以映秀—北川地表破裂为近似中心的近椭圆形区域内。分析这196007处滑坡面积与滑坡数量的幂律关系,结果表明,面积在1x10⁴—1x10⁶m2之间的滑坡,滑坡面积与滑坡数量累加之间的幂律关系式可表示为lgN(A)=-2.0745A+13,判定系数R²=0.9931,表明该滑坡面积区间内滑坡面积与滑坡数量之间具有良好的幂律关系。对比结果表明,汶川地震滑坡影响区面积略高于全球其他地震事件滑坡影响区面积的上包络线。汶川地震滑坡编录成果为后续汶川地震滑坡相关科学研究提供了重要的基础数据。     

Rupture process of the Wenchuan <Emphasis Type="Italic">M</Emphasis>8.0 earthquake of Sichuan China: the segmentation feature

Moment tensor solution, rupture process and rupture characteristics of the great Wenchuan M8.0 earthquake are studied by using 39 long-period P and SH waveforms with evenly azimuth coverage of stations. Our results reveal that the Wenchuan M8.0 event consisted of 5 sub-events of Mw≥7.3 occurring succesively in time and space. Rupture started with a Mw7.3 introductory strike-slip faulting in the first 12 s, then within 12?40 s, two sub-events with Mw7.6 and Mw7.4 occurred within 80 km northeast from the init...     

The M S7.1 Yushu earthquake surface rupture and large historical earthquakes on the Garzê-Yushu Fault

As revealed by field investigations, the co-seismic surface rupture zone of the 2010 M_S7.1 Yushu earthquake, Qinghai is a char-acteristic sinistral strike-slip feature consisting of three distinct sinistral primary ruptures, with an overall strike of 310°–320° and a total length of 31 km. In addition, an approximately 2-km-long en-echelon tensile fissure zone was found east of Longbao Town; if this site is taken as the north end of the rupture zone, then the rupture had a total length of ~51 km. The surface rupture zone is composed of a series of fissures arranged in an en-echelon or alternating relationship between compressive bulges and tensile fissures, with a measured maximum horizontal displacement of 1.8 m. The surface rupture zone extends along the mapped Garzê-Yushu Fault, which implicates it as the seismogenic fault for this earthquake. Historically, a few earthquakes with a magnitude of about 7 have occurred along the fault, and additionally traces of paleoearthquakes are evident that characterize the short-period recurrence interval of large earthquakes here. Similar to the seismogenic process of the 2008 Wenchuan earthquake, the Yushu earthquake is also due to the stress accumulation and release on the block boundaries resulting from the eastward expansion of Qinghai-Tibet Plateau. However, in contrast with the Wenchuan earthquake, the Yushu earthquake had a sinistral strike-slip mechanism resulting from the uneven eastward extrusion of the Baryan Har and Sichuan-Yunnan fault blocks.     

Evidence for fault weakness and fluid flow within an active low-angle normal fault.

Determining the composition and physical properties of shallow-dipping, active normal faults (dips < 35 degrees with respect to the horizontal) is important for understanding how such faults slip under low resolved shear stress and accommodate significant extension of the crust and lithosphere. Seismic reflection images and earthquake source parameters show that a magnitude 6.2 earthquake occurred at about 5 km depth on or close to a normal fault with a dip of 25-30 degrees located ahead of a propagating spreading centre in the Woodlark basin. Here we present results from a genetic algorithm inversion of seismic reflection data, which shows that the fault at 4-5 km depth contains a 33-m-thick layer with seismic velocities of about 4.3 km s(-1), which we interpret to be composed of serpentinite fault gouge. Isolated zones exhibit velocities as low as approximately 1.7 km s(-1) with high porosities, which we suggest are maintained by high fluid pressures. We propose that hydrothermal fluid flow, possibly driven by a deep magmatic heat source, and high extensional stresses ahead of the ridge tip have created conditions for fault weakness and strain localization on the low-angle normal fault.     

Predicting the endpoints of earthquake ruptures

The active fault traces on which earthquakes occur are generally not continuous, and are commonly composed of segments that are separated by discontinuities that appear as steps in map-view. Stress concentrations resulting from slip at such discontinuities may slow or stop rupture propagation and hence play a controlling role in limiting the length of earthquake rupture. Here I examine the mapped surface rupture traces of 22 historical strike-slip earthquakes with rupture lengths ranging between 10 and 420 km. I show that about two-thirds of the endpoints of strike-slip earthquake ruptures are associated with fault steps or the termini of active fault traces, and that there exists a limiting dimension of fault step (3-4 km) above which earthquake ruptures do not propagate and below which rupture propagation ceases only about 40 per cent of the time. The results are of practical importance to seismic hazard analysis where effort is spent attempting to place limits on the probable length of future earthquakes on mapped active faults. Physical insight to the dynamics of the earthquake rupture process is further gained with the observation that the limiting dimension appears to be largely independent of the earthquake rupture length. It follows that the magnitude of stress changes and the volume affected by those stress changes at the driving edge of laterally propagating ruptures are largely similar and invariable during the rupture process regardless of the distance an event has propagated or will propagate.     

汶川地震序列震源机制解研究

利用2008年5月12日至8月7日包括国家地震台网、 地方区域地震台网和流动地震台网在内的1023个地震台站及IRIS台站的数字地震波形记录资料, 读取P波初动极性, 并采用改进的格点尝试法对汶川地震序列震源机制解进行求解, 给出121个4.0级以上质量可靠的震源机制解。根据求解结果发现: 汶川Ms 8.0级主震P波初动解为纯逆冲型; 汶川余震序列震源机制解绝大部分是逆冲型或走滑型地震, 并具有分段分布特征; 其P轴方位有两个优势分布方向: 北西西?南东东向和北东东?南西西向。     

Source rupture process of Lushan M S7.0 earthquake, Sichuan, China and its tectonic implications

The source rupture process of the M _S7.0 Lushan earthquake was here evaluated using 40 long-period P waveforms with even azimuth coverage of stations. Results reveal that the rupture process of the Lushan M _S7.0 event to be simpler than that of the Wenchuan earthquake and also showed significant differences between the two rupture processes. The whole rupture process lasted 36 s and most of the moment was released within the first 13 s. The total released moment is 1.9×10¹⁹N m with M _W=6.8. Rupture propagated upwards and bilaterally to both sides from the initial point, resulting in a large slip region of 40 km×30 km, with the maximum slip of 1.8 m, located above the initial point. No surface displacement was estimated around the epicenter, but displacement was observed about 20 km NE and SW directions of the epicenter. Both showed slips of less than 40 cm. The rupture suddenly stopped at 20 km NE of the initial point. This was consistent with the aftershock activity. This phenomenon indicates the existence of significant variation of the medium or tectonic structure, which may prevent the propagation of the rupture and aftershock activity. The earthquake risk of the left segment of Qianshan fault is worthy of attention.     

地震监测预报的误区分析与工作建议

四川省汶川县巨震之后,人们对地震前兆问题广为讨论,由此引发诸多思考。基于地球系统复杂变化过程中的突变现象,分析了地震前兆信息的基本问题;指出目前社会公众及国际地震学界一般存在以下3个误区,即认为前兆与地震应一一对应,陆区地震预报主要依赖于台网监测,地震前兆分析方法比较单一;进而提出中国地震监测预报工作应重视卫星对地观测数据的快速处理与分析,建立地震数据共享机制和地震信息汇合机构,期望中国地震监测预报理论和技术尽快实现创新发展。     

North China sub-craton lithospheric structure elucidated through coal mine blasting

We present a super-range seismic observation along the >1300-km-long profile passing through the Yinchuan basin and the Ordos block from the blasting point towards the southeast triggered by a large-(dynamite) scale coal blast in the Ningxia Hui Autonomous Region’s Helan Mountain (yielded a profile). The seismic wave information from the uppermost mantle reflecting different depths was obtained by the China continental seismic survey. Pn refracted waves from the uppermost mantle were effectively traced up t...     

Synthetic seismograms of ground motion near earthquake fault using simulated Green’s function method

Seismograms near source fault were synthesized using the hybrid empirical Green＇s function method where the discretely simulated seismic waveforms are used for Green＇s functions instead of the observed waveforms of small earthquakes. The Green＇s function seismic waveforms for small earthquake were calculated by solving wave equation using the pseudo-spectral method with the staggered grid real FFT strategy under a detailed 2-D velocity structure in Kobe region. Magnitude and seismic moment of simulated Green＇s function waveforms were firstly determined by using the relationship between fault length and corner frequency of source spectrum. The simulated Green＇s function waveforms were employed to synthesize seismograms of strong ground motion near the earthquake fault. The synthetic seismograms of the target earthquake were performed based on the model with multiple source rupture processes. The results suggest that synthesized seismograms coincide well with observed seismic waveforms of the 1995 Hyogo-ken Nanbu earthquake. The simulated Green＇s function method is very useful for prediction of the strong ground motion in region without observed seismic waveforms. The present technique spreads application field of the empirical Green＇s function method.     

联合slipBERI与D-InSAR技术的青海省玛多县地震形变场提取与模拟

针对2021年5月22日青海玛多县Ms7.4地震震区形变信息、形变特征以及滑动断裂特性的提取与模拟等问题,获取了玛多县地震区的Sentinel-1A影像,采用双轨差分干涉法,并优化各项参数,提取出同震形变场,利用slipBERI(slip from Bayesian Regularized Inversion)方法对断层的几何参数及形变场进行反演和模拟。结果表明：玛多地震同震形变场的形状近似于一个椭圆,断层整体呈西北-东南走向,其上部为沉降区,下部为隆升区,最大LOS(line of sight,视线向)形变分别为0.65m和0.81m。地震形变场的运动主要以东西方向的水平运动为主,并伴有明显的左旋走滑,断层上下方相对视线向运动可达1.50m,表明此次地震的地表破裂有明显的错位移动。通过分析形变信息和地表破裂特征,可以判断该破裂带位于巴颜喀拉块体,为昆仑山口-江口断裂,反演结果与观测结果相符,这表明观测结果较可靠。     

Fine three-dimensional P-wave velocity structure beneath the capital region and deep environment for the nucleation of strong earthquakes

A detailed 3-D P-wave velocity model of the crust and uppermost mantle under the capitol region is determined with a spatial resolution of 25 km in the horizontal direction and 4-17 km in depth. We used 48750 precise P-wave arrival time data from 2973 events of local crustal earthquakes, controlled seismic explosions and quarry blasts. These events were recorded by 123 seismic stations. The data are analyzed by using a 3-D seismic tomography method. Our tomographic model provides new information on the geological structure and complex seismotectonics of this region. Different patterns of velocity structures show up in the North China Basin, the Taihangshan and the Yanshan Mountainous areas. The velocity images of the upper crust reflect well the surface geological, topographic and lithological features. In the North China Basin, the depression and uplift areas are imaged as slow and fast velocity belts, respectively, which are oriented in NE-SW direction. The trend of velocity anomalies is the same as that of major structure and tectonics. Paleozoic strata and Pre-Cambrian basement rocks outcrop widely in the Taihangshan and Yanshan uplift areas, which exhibit strong and broad high-velocity anomalies in our tomographic images, while the Quaternary intermountain basins show up as small low-velocity anomalies. Most of large earthquakes, such as the 1976 Tangshan earthquake (M 7.8) and the 1679 Sanhe earthquake (M 8.0), generally occurred in high-velocity areas in the upper to middle crust. However, in the lower crust to the uppermost mantle under the source zones of the large earthquakes, low-velocity and high-conductivity anomalies exist, which are considered to be associated with fluids, just like the 1995 Kobe earthquake (M 7.2) and the 2001 Indian Bhuj earthquake (M 7.8). The fluids in the lower crust may cause the weakening of the seismogenic layer in the upper and middle crust and thus contribute to the initiation of the large crustal earthquakes.     

Upper mantle SH velocity structure beneath Qiangtang Terrane by modeling triplicated phases

We constrain SH wave velocity structure for the upper mantle beneath western Qiangtang Terrane by comparing regional distance seismic triplicated waveforms with synthetic seismograms, based on an intermediate event （-220 km） recorded by the INDEPTH-Ⅲ seismic array. The ATIP model reveals a low-velocity anomaly with up to -4% variation at the depth of 190-270 km and a relatively small velocity gradient above the depth of 410 km in the upper mantle, which is in agreement with previous results. In combination with other geological studies, we suggest that the depth of top asthenosphere is 190 km and no large-scale lithosphere thinning occurs in western Qiangtang Terrane, besides, Qiangtang Terrane has the same kind of upper mantle structure as the stable Eurasia.     

Preliminary results pertaining to coseismic displacement and preseismic strain accumulation of the Lushan M S7.0 earthquake, as reflected by GPS surveying

This paper presents the coseismic displacement and preseismic deformation fields of the Lushan M _S7.0 earthquake that occurred on April 20, 2013. The results are based on GPS observations along the Longmenshan fault and within its vicinity. The coseismic displacement and preseismic GPS results indicate that in the strain release of this earthquake, the thrust rupture is dominant and the laevorotation movement is secondary. Furthermore, we infer that any possible the rupture does not reach the earth’s surface, and the seismogenic fault is most likely one fault to the east of the Guanxian-Anxian fault. Some detailed results are obtainable. (1) The southern segment of the Longmenshan fault is locked preceding the Lushan earthquake. After the Wenchuan earthquake, the strain accumulation rate in the southeast direction accelerates in the epicenter of the Lushan earthquake, and the angle between the principal compressional strain and the seismogenic fault indicates that a sinistral deformation background in the direction of the seismogenic fault precedes the Lushan earthquake. Therefore, it is evident that the Wenchuan M _S8.0 earthquake accelerated the pregnancy of the Lushan earthquake. (2) The coseismic displacements reflected by GPS data are mainly located in a region that is 230 km (NW direction) × 100 km (SW direction), and coseismic displacements larger than 10 mm lie predominantly in a 100-km region (NW direction). (3) On a large scale, the coseismic displacement shows thrust characteristics, but the associated values are remarkably small in the near field (within 70 km) of the earthquake fault. Meanwhile, the thrust movement in this 70-km region does not correspond with the attenuation characteristics of the strain release, indicating that the rupture of this earthquake does not reach the earth’s surface. (4) The laevorotation movements are remarkable in the 50-km region, which is located in the hanging wall that is close to the earthquake fault, and the corresponding values in this case correlate with the attenuation characteristics of the strain release.