日本东部内陆地区2011年大地震前后的应力状态与诱发地震

通过反演震源机制数据,得到了东日本大地震前后日本东部内陆地区的应力场。地震前,在东北部(Tohoku),σ1轴呈EW方向,但在关东-东部(Kanto-Chubu)地区则呈NW-SE方向。地震后,在Tohoku北部和靠近磐城市(Iwaki City)的Tohoku东南部,应力场发生了变化,主应力取向大致等同于与地震相关的静应力变化取向。这说明地震前这些地区的差应力值小于1MPa。在Tohoku中部,虽然地震后加载的应力取向几乎逆转,但其应力场没有发生变化,说明那里的差应力值大大超过1MPa。在Kanto-Chubu,地震后加载了与背景应力方向近乎相同的应力,且应力场没有发生预期的变化。这可能就是Kanto-Chubu地区高地震活动水平的诱因。     

汶川大地震前后GPS观测的精密单历元解得到的震区地壳运动

利用BERNESE软件的单历元观测值双差定位功能,获取了接近2008年5月12日汶川Ms8.0大地震震中5个GPS观测站(PIXI,CHDU,JYAN,NEIJ和LUZH站)震前的短临位移过程及震后短时间的位移.通过拟合和对比分析,讨论了单历元解的精度.计算结果表明,以YANC站为固定站,LUZH站震前5月9日N(南北...     

Strong motion displacement waveforms using 10‐Hz precise point positioning GPS: an assessment based on free oscillation experiments

The recently developed precise point positioning (PPP) technique permits to compute instantaneous coordinates of a GPS station relatively to distant reference stations and waveforms of ground displacements during strong motions at the 1 Hz level. This is another application of GPS, different from the computation of static coseismic movements or of accurate monitoring of dynamic displacements of structures using a static receiver and a nearby moving receiver (DGPS). Recently, earthquake ground displacement waveforms using 10‐Hz GPS data have also been calculated, but no independent evidence to assess their quality exists. To overcome this problem, we evaluated the output of 10‐Hz PPP results on the basis of supervised learning experiments. Semistatic and dynamic displacements (damped harmonic oscillations) of known characteristics of the order of a few centimeter were produced and were recorded by GPS, an accelerometer, and a robotic total station. Time series of instantaneous displacements were analyzed using different PPP techniques and were compared with reference (true) values derived from DGPS and the other sensors. Our analysis revealed that the PPP‐derived coordinates are contaminated by long‐period noise but they can display the details of semistatic displacements, while their short‐period component describes well the pattern of waveforms and spectra (at least up to 4 Hz) of dynamic displacements, with up to 20 mm accuracy for isolated points. These results indicate that 10‐Hz PPP‐GPS is useful for earthquake engineering and can safely be used to reconstruct waveforms of deflections of the ground and of various points on structures during strong motions. Copyright © 2014 John Wiley & Sons, Ltd.     

利用日本GPS网探测2011年Tohoku海啸引发的电离层扰动

海平面的海啸波会产生大气重力波进而引发电离层扰动.本文利用日本GPS总电子含量数据来探测2011年3月11日Tohoku海啸引发的电离层扰动.观测结果表明,在日本上空的电离层中存在两种重力波信号,分别由海平面的海啸波以及地震破裂过程产生.地震产生的电离层重力波分布在震中周围（包括海洋上空以及远离海洋的区域）,而海啸引发的电离层重力波主要分布在海洋上空.地震产生的电离层重力波具有不同的水平速度,包括约210 m·s^-1以及170 m·s^-1,其频率为1.5 mHz;而海啸引发的电离层重力波水平速度快于前者,约为280 m·s^-1,其频率为1.0 mHz.此外,海啸引发电离层重力波与海平面上的海啸波有相似的水平速度、方向、运行时间、波形以及频率等传播特征.本文的研究将电离层中的海啸信号与地震信号区分开来,进一步确认电离层对海啸波的敏感性.

     

基于GPS同震位移场约束反演2008年5.12汶川大地震破裂空间分布

2008年5月12日发生在四川汶川的大地震造成映秀—北川断裂和灌县—江油断裂同时破裂,分别形成了240多公里和70多公里的地表破裂带.本文以GPS观测获得的同震位移场为约束,反演地震破裂的空间分布.反演结果显示映秀—北川主破裂带倾向北西,沿破裂带的走向从南到北倾角逐渐变大,破裂断层的平均宽度在10~18 km左右.破裂断层的错动在南段以逆冲为主,在北段走滑分量逐步加大,右旋走滑成为断层破裂的主要特征.断层破裂最大段落错动量分别达到了7.8 m和7.4 m,恰好对应这次地震中地表破坏最为严重的映秀和北川地区.本次地震释放地震矩6.70×1020N·m,相应矩震级Mw=7.9.     

用多种数据构建2008年汶川特大地震同震位移场

本文主要以GPS、精密水准观测和卫星SAR遥感图像分析2008年汶川特大地震同震位移特征.GPS数据包括:(1)四川盆地和川西高原地区各类国家等级GPS网点复测;(2)沿破裂带国家天文大地网GPS复测.前者推算的同震位移测定精度优于2 cm,后者6~8 cm.SAR遥感资料包括:(1)ALOS 卫星升轨相位干涉图像,精度优于8 cm;(2)ALOS和ENVISAT卫星影像合成的三维位移图,精度优于0.5 m.同震位移场显示,断层下盘(四川盆地)变形总体呈扇形集中指向震中,断层上盘(龙门山)变形总体上呈逆时针旋转态势,最大的实测水平位移5.5 m.汶川、理县、茂县等地测站位移指向破裂带方向,而平武、青川等地测站逐渐转变为平行,乃至远离破裂带方向,与汶川地震逆冲兼走滑的破裂特征一致.断层上盘大幅隆升,下盘靠近断层的区域以下沉为主,远场表现为幅度很小的隆升,垂直升降区域间,有一条与龙泉山断裂带平行的升降过渡带,调节龙泉断层的应力状态.用实测变形场检验多个地震波破裂模型表明,近场(距离断层50 km) 模型形变准确度可达40~50 cm, 远场精度优于5 cm.     

用多种数据构建2008年汶川特大地震同震位移场

本文主要以GPS、精密水准观测和卫星SAR遥感图像分析2008年汶川特大地震同震位移特征.GPS数据包括:(1)四川盆地和川西高原地区各类国家等级GPS网点复测;(2)沿破裂带国家天文大地网GPS复测.前者推算的同震位移测定精度优于2 cm,后者6~8 cm.SAR遥感资料包括:(1)ALOS 卫星升轨相位干涉图像,精度优于8 cm;(2)ALOS和ENVISAT卫星影像合成的三维位移图,精度优于0.5 m.同震位移场显示,断层下盘(四川盆地)变形总体呈扇形集中指向震中,断层上盘(龙门山)变形总体上呈逆时针旋转态势,最大的实测水平位移5.5 m.汶川、理县、茂县等地测站位移指向破裂带方向,而平武、青川等地测站逐渐转变为平行,乃至远离破裂带方向,与汶川地震逆冲兼走滑的破裂特征一致.断层上盘大幅隆升,下盘靠近断层的区域以下沉为主,远场表现为幅度很小的隆升,垂直升降区域间,有一条与龙泉山断裂带平行的升降过渡带,调节龙泉断层的应力状态.用实测变形场检验多个地震波破裂模型表明,近场(距离断层50 km) 模型形变准确度可达40~50 cm, 远场精度优于5 cm.     

2011年日本MW9.0地震引发的海啸对地震背景噪声的影响

基于北美沿岸和内陆地震台站的连续地震波形记录,并结合沿岸台站附近布设的DART系统记录的海底压力数据以及预测潮汐数据,利用时频分析和极化分析方法对2011年3月11日日本东北部海域MW9.0大地震所激发的海啸对地震背景噪声所产生的影响予以深入分析.结果显示:海啸对高频噪声(1.3—1.5 Hz)以及短周期双频微地动噪声...     

用钻孔应变数据探测2011年日本Tohoku大地震激发的低频地球环型自由振荡及其谱峰分裂

地震激发的地球自由振荡中, 环型自由振荡(特别是低频环型自由振荡)的精确探测相比于球型自由振荡更加困难. 环型自由振荡对SH波速度结构更加敏感, 能更好地监测走滑断层运动, 因此精确估算环型自由振荡的谱峰参数可为地球内部大尺度SH波速度结构和大地震的震源机制提供更强约束. 本研究利用2011年日本Tohoku地震后四分量钻孔应变仪的观测资料, 提取了该地震激发的1.3 mHz以下低频环型和球型自由振荡信号, 与STS-1宽频带地震仪的结果相比较, 发现四分量钻孔应变仪对1 mHz以下的环型自由振荡具有更强的探测能力. 本研究还采用频率域自回归方法精确估计了₀T₂—₀T₆的简并频率, 结果与PREM模型理论值相差非常小, 平均相对误差仅为2.52×10^-4.此外还用钻孔应变仪观测资料探测到了₀T₂谱峰分裂现象. 研究结果表明, 钻孔应变仪能有效探测到大地震激发的低频地球环型自由振荡信号, 可为地幔大尺度各向异性结构研究提供高质量的基础数据.

     

2008年汶川大地震中北川地区极重震害的物理机制研究

2008年汶川特大地震中北川地区受灾特别严重,而该地区远离汶川主震震中超过100多千米.尽管汶川地震发生已近10年,但这种不正常现象一直困扰着地学工作者,至今没有给出合理的解答.为此,本文利用有限单元方法模拟汶川地震的主要发震断裂带（即：映秀—北川断层）的自发破裂动力学过程.模型中,几乎平直的映秀—北川断裂带在高川地区发生拐折,形成高川右弯.模拟结果显示,破裂沿着映秀—北川断裂带向东北方向前行,当破裂前端到达高川右弯时,破裂形态发生剧烈变化,破裂传播速度由亚剪切破裂突然转化为超剪切破裂.与断层阶区促进超剪切破裂明显不同的是,这种超剪切破裂形态的转化不需要时间停顿.由于超剪切地震破裂的产生,破裂辐射的地震波发生相长干涉,形成马赫波,地震动被大大放大.计算给出的强地面运动峰值加速度在北川地区不仅数值高,而且高值区分布范围广,造成了北川地区的震害特别严重.模拟结果还表明,若高川右弯不存在（即映秀—北川断层面为平直）,则北川地区的断层不会形成超剪切破裂;如果高川右弯附近断层的不连续程度较大,则破裂会被终止;这些情况下,北川地区的震害都没有那么严重.所以,高川右弯的几何形状对于超剪切破裂的形成起着决定性作用,也是造成远离震中的北川地区震害特别严重的主要原因.因此,分析发震断层几何对于研究震源动力学过程及震害评估等有着非常重要的科学意义和实际应用价值.

     

2003年日本北海道8.0级地震次声波特征研究

利用在北京昌平安装的次声三点阵, 记录了2003年9月26日日本北海道地震的前兆次声波和震后次声波。 实际记录的P-t曲线及波速波向图, 经快速傅立叶变换得到其三维动态频谱。 结果表明: ① 强震前2～7天能测到振幅很强、 方向可测的地震前兆次声波。 其频谱特点是: 振幅由弱(10～15 Pa)逐渐加强到(50～80 Pa)或更大; 先为短周期波(10～30 min), 紧跟着长周期波(30～50 min), 然后又出现更长的周期波, 最后长短周期的波一起出现连成一片(2～62 min); ② 地震前兆次声波的产生和记录较容易, 而当地次声波和震中次声波记录和识别比较困难; ③ 若两地的三点阵波向相交, 可预测震中位置。 故地震前兆次声波的测量研究, 有可能发展成为临震预报中一种有效的新方法     

Stress field change around the Mount Fuji volcano magma system caused by the Tohoku megathrust earthquake, Japan

Crustal deformation by the M _w 9.0 megathrust Tohoku earthquake causes the extension over a wide region of the Japanese mainland. In addition, a triggered M _w 5.9 East Shizuoka earthquake on March 15 occurred beneath the south flank, just above the magma system of Mount Fuji. To access whether these earthquakes might trigger the eruption, we calculated the stress and pressure changes below Mount Fuji. Among the three plausible mechanisms of earthquake–volcano interactions, we calculate the static stress change around volcano using finite element method, based on the seismic fault models of Tohoku and East Shizuoka earthquakes. Both Japanese mainland and Mount Fuji region are modeled by seismic tomography result, and the topographic effect is also included. The differential stress given to Mount Fuji magma reservoir, which is assumed to be located to be in the hypocentral area of deep long period earthquakes at the depth of 15 km, is estimated to be the order of about 0.001–0.01 and 0.1–1 MPa at the boundary region between magma reservoir and surrounding medium. This pressure change is about 0.2 % of the lithostatic pressure (367.5 MPa at 15 km depth), but is enough to trigger an eruptions in case the magma is ready to erupt. For Mount Fuji, there is no evidence so far that these earthquakes and crustal deformations did reactivate the volcano, considering the seismicity of deep long period earthquakes.     

Site effects by generalized inversion technique using strong motion recordings of the 2008 Wenchuan earthquake

The generalized inversion of S-wave amplitude spectra from the free-field strong motion recordings of the China National Strong Motion Observation Network System (NSMONS) are used to evaluate the site effects in the Wenchuan area. In this regard, a total of 602 recordings from 96 aftershocks of the Wenchuan earthquake with magnitudes of M3.7-M6.5 were selected as a dataset. These recordings were obtained from 28 stations at a hypocenter distance ranging from 30 km to 150 km. The inversion results have been verified as reliable by comparing the site response at station 62WUD using the Generalized Inversion Technique (GIT) and the Standard Spectral Ratio method (SSR). For all 28 stations, the site predominant frequency F _p and the average site amplification in different frequency bands of 1.0–5.0 Hz, 5.0–10.0 Hz and 1.0–10.0 Hz have been calculated based on the inversion results. Compared with the results from the horizontal-to-vertical spectral ratio (HVSR) method, it shows that the HVSR method can reasonably estimate the site predominant frequency but underestimates the site amplification. The linear fitting between the average site amplification for each frequency band and the V _s20 (the average uppermost-20 m shear wave velocity) shows good correlation. A distance measurement called the asperity distance D _Aspt is proposed to reasonably characterize the source-to-site distance for large earthquakes. Finally, the inversed site response is used to identify the soil nonlinearity in the main shock and aftershocks of Wenchuan earthquake. In ten of the 28 stations analyzed in the main shock, the soil behaved nonlinearly, where the ground motion level is apparently beyond a threshold of PGA > 300 cm/s² or PGV > 20 cm/s, and only one station coded 51SFB has evidence of soil nonlinear behavior in the aftershocks.     

Paleoseismic events and the 1596 Keicho-Fushimi large earthquake produced by a slip on the Gosukebashi fault at the eastern Rokko Mountains,Japan

Field surveys and trench excavation investigations revealed that there were at least four large seismic events produced by slips on the Gosukebashi fault in the Holocene in the southeastern Rokko Mountains of Japan. The characteristics of deformed topographies and three-dimensionally excavated exposures show that this fault is a right-lateral strike–slip fault having an average slip rate of 1.0 mm/year, with a reverse displacement component. The principle indicators of past faulting events are: (i) termination of secondary faults; (ii) sedimentary deposits related to faulting; and (iii) injection veins of fault gouge related to seismic faulting in the fractured zone. Radiocarbon dates indicate that the events occurred pre-1660 BC , 1660 BC –220 AD , from ～ 30–220 to 600 AD and 15th century AD . The youngest event is probably associated with the large 1596 AD Keicho-Fushimi earthquake which occurred in the area around Kyoto and Kobe Cities. The second younger event is probably correlated with the 416 AD Yamato earthquake, which is the oldest historic earthquake in Japanese historic records. The results of trench surveys show that the horizontal displacement produced by an individual event is ～ 1.5 m, and the recurrence of seismic event intervals is ～ 1200 years in the Gosukebashi fault.     

Empirical fragility analysis of building damage caused by the 2011 Great East Japan tsunami in Ishinomaki city using ordinal regression,and influence of key geographical features

Tsunamis are disastrous events typically causing loss of life, and extreme damage to the built environment, as shown by the recent disaster that struck the East coast of Japan in 2011. In order to quantitatively estimate damage in tsunami prone areas, some studies used a probabilistic approach and derived fragility functions. However, the models chosen do not provide a statistically sound representation of the data. This study applies advanced statistical methods in order to address these limitations. The area of study is the city of Ishinomaki in Japan, the worst affected area during the 2011 event and for which an extensive amount of detailed building damage data has been collected. Ishinomaki city displays a variety of geographical environments that would have significantly affected tsunami flow characteristics, namely a plain, a narrow coast backed up by high topography (terrain), and a river. The fragility analysis assesses the relative structural vulnerability between these areas, and reveals that the buildings surrounding the river were less likely to be damaged. The damage probabilities for the terrain area (with relatively higher flow depths and velocities) were lower or similar to the plain, which confirms the beneficial role of coastal protection. The model diagnostics show tsunami flow depth alone is a poor predictor of tsunami damage for reinforced concrete and steel structures, and for all structures other variables are influential and need to be taken into account in order to improve fragility estimations. In particular, evidence shows debris impact contributed to at least a significant amount of non-structural damage.     

2011年3月11日日本地震海啸越洋传播及对中国影响的数值分析

北京时间2011年3月11日13时46分(05:46 UTC)日本东北部近海(38.3°N,142.4°E)发生M_w9.0级特大地震,此次地震的强度为日本近1200a来最强.随后环太平洋的数十个国家和地区的验潮站和海啸监测浮标均监测到了强震引发的越洋海啸,海啸奔袭23 h到达南美洲的智利沿岸;此次海啸除了对近场的日本东北部沿岸地区造成了巨大灾害,还对太平洋东岸的部分国家和地区造成了一定程度的影响.地震发生4 h后海啸波到达我国台湾东部沿海,6~8 h海啸波到达我国大陆东南沿海,受此影响我国发布了第一份海啸蓝色警报.本文利用海啸数值模型对此次地震海啸的产生、越洋传播过程进行了数值模拟,给出了海啸波能量在我国近海及泛太平洋区域分布特征;同时重点模拟分析了海啸波在日本及中国近海传播的波动特征,模拟结果与观测数据吻合良好.最后通过对数值模拟结果的分析,阐述了此次海啸对中国的影响,给出了潜在的日本地震海啸对中国的风险估计.     

2011年3月11日日本地震海啸越洋传播及对中国影响的数值分析

北京时间2011年3月11日13时46分(05:46 UTC)日本东北部近海(38.3°N,142.4°E)发生Mw9.0级特大地震,此次地震的强度为日本近1200a来最强.随后环太平洋的数十个国家和地区的验潮站和海啸监测浮标均监测到了强震引发的越洋海啸,海啸奔袭23 h到达南美洲的智利沿岸;此次海啸除了对近场的日本东北部沿岸地区造成了巨大灾害,还对太平洋东岸的部分国家和地区造成了一定程度的影响.地震发生4 h后海啸波到达我国台湾东部沿海,6~8 h海啸波到达我国大陆东南沿海,受此影响我国发布了第一份海啸蓝色警报.本文利用海啸数值模型对此次地震海啸的产生、越洋传播过程进行了数值模拟,给出了海啸波能量在我国近海及泛太平洋区域分布特征;同时重点模拟分析了海啸波在日本及中国近海传播的波动特征,模拟结果与观测数据吻合良好.最后通过对数值模拟结果的分析,阐述了此次海啸对中国的影响,给出了潜在的日本地震海啸对中国的风险估计.     

Automatic recognition of damaged town buildings caused by earthquake using remote sensing information: Taking the 2001 Bhuj, India, earthquake and the 1976 Tangshan, China, earthquake as examples

Introduction As we well know, the hazard of earthquake is very wide especially in cities. The conventionalmethods to investigate the damage are difficult to meet the requirements in applications. In recentyears, with the rapid development of remote sensing, especially the successful launch and applica-tion of high-resolution commercial remote sensing satellite, it has become possible to recognize andextract damage information by using remote sensing. The researchers at home and abroad hav…     

Analysis of landslide damage caused by the 2008 Wenchuan earthquake using strong motion data: A case study in the Beichuan county town

The county town of Beichuan county, China, experienced catastrophic destruction due to landslides induce by the 2008 Wenchuan earthquake. In consideration of the special location of the county town, this paper selected the landslides induced in the town as representative of large-scale near-rupture landslides, and quantitatively analyzed why the landslide damage was so destructive in the town by using strong motion data obtained from the Wenchuan earthquake in the Longmenshan area. Three methods were employed to estimate the landslide damage using strong motion data. (1) Peak ground accelerations (PGAs) on the hanging wall were used to evaluate the PGAs on the landslide sites in the town. The evaluated average PGAs were all greater than 1 g, indicating that the ground motion intensity was very strong during the earthquake. (2) Acceleration time histories, from another station with similar geological conditions to the town, were used to evaluate the critical acceleration changing range, and the estimated values showed the geological conditions were very susceptible to earthquakes. (3) Acceleration time histories, from two stations on the hanging and foot walls of the rupture, and near the town, were used to calculate the Newmark displacements, and all the evaluated displacements indicated that landslides were very likely. The results show that the slopes, susceptible to earthquakes in the Beichuan county town, were easily triggered under such strong ground-motion intensity and developed into large-scale catastrophic events.     

利用VP型垂直摆倾斜仪观测数据检测2011日本M_w9.0级地震激发的低频地球自由振荡

甚长周期VP型垂直摆倾斜仪对某些地球物理信号有特殊的敏感性,除了在观测小地震、慢地震等方面有优势之外,还可有效应用于低频地球自由振荡信号的探测.利用安置在武汉大学珞珈山的我国自行研制的VP型垂直摆倾斜仪在2011日本Mw9.0级大地震之后不同长度的观测数据,联合EEMD方法、自回归估计(AR)方法和bootstrap法,本文不仅检测到该地震激发的频率小于4.7 mHz的零阶球型振荡(0S2至0S38)和环型振荡(0T4至0T35)几乎所有振型以及15个谐频振型,还检测到5个其频率低于1mHz的低阶球型多线态(0S2、2S1、0S3、0S4和1S2)的部分或全部谱峰分裂现象,并给出了所有检测结果的精度评估.此外,本文分析了某些球型和环型振荡之间的耦合效应,结果表明耦合效应将显著影响地球自由振荡信号的相关参数.