2011年东日本大地震:位移延伸至海沟轴

2011年3月11日日本地震(Mw9.0)后发生的大海啸被认为是由延伸至日本海沟处一条俯冲带浅部的断层破裂所致.这是不同的地震和大地测量反演程序所揭示的结果[1-2];然而,精确的同震位移逆倾斜向上极限值仍然没有被明确界定.在此,我们介绍了地震前和地震后穿越破裂带内的海沟重复进行的多波束水深测量结果,以估算其逆倾斜向上极限值,并量化海底位移.     

2011年日本东北近海地震震源上方的位移

2011年3月11日,一次大的板间地震(MW9.0)发生在日本东北宫城县附近的板块交界处。根据余震分布推测的震源区为约500km长、200km宽的近海区(Japan Meteo-rological Agency,2011)。     

日本东北9.0级地震的同震与震后滑动

大部分强震都发生在海沟,那里是海洋板块向大陆板块俯冲的地方.大量矩震级MW9.0以上的地震发生在若干区域,包括智利,阿拉斯加,堪察加半岛和苏门答腊岛等.位于太平洋板块俯冲鄂霍茨克板块的日本海沟,历史记载上没有发生过MW9.0地震,除了至今震级还有争议的公元869年Jogan大地震[1](可能超过MW9.0).然而,根据...     

2011年日本MW9.0地震相关研究综述

2011年3月11日, 日本海沟发生的9级地震造成重大人员伤亡, 受到社会普遍关注, 本文基于此次日本9级地震相关研究结果, 尝试从不同侧面分析此次地震的观测、 现象和认识, 主要包括如下几点: ① 此次地震发生在太平洋板块西北边界上日本海沟俯冲带上, 同震破裂可能存在深浅两个位错集中区, 较深的位错集中区位错量相对较小, 但历史上7级地震多发; 而较浅的位错集中区位错较大, 但历史上强震活动相对较弱; ② 基于GPS观测资料为约束的相关断层运动研究结果表明, 日本海沟断层运动背景以大范围稳定闭锁为主(闭锁区空间尺度与同震破裂尺度相当), 自2003年日本北海道8级地震后日本海沟地区断层运动开始出现扰动, 2008年以后有几次7级左右地震震后余滑分布明显比主震位错量要大, 之后分别于2008年和2011年观测到显著慢滑移事件, 最后分别于2011年3月9日和3月11日发生7级前震和9级主震, 震前日本海沟俯冲带断层运动变化过程比较清楚; ③ 可能是由于监测的原因, 传统上的前兆观测并未出现显著异常, 其震前异常主要为: 部分地震活动参数表明强震震源区震前应力状态相对较高、 区域地表运动速率的短期异常等; ④ 对于震源区物理性质的分析引起了更多的科学问题, 例如, 震源区介质物性是否与周边存在显著差异、 断层摩擦性质是否决定了发震能力和破裂过程、 震前断层运动是否存在预滑、 震前深部流体是否影响到震源区断层运动等。 他山之石可以攻玉, 希望本文对地震预测预报基础研究工作能起到抛砖引玉的作用。     

2011年日本9级大地震的同震和震后滑移

多数大地震都是沿大洋板块向大陆板块俯冲的海沟发生.现在已经清楚矩震级Mw达到9级的巨大地震只发生在几个地区,包括智利、阿拉斯加、堪察加半岛和苏门答腊.日本海沟是太平洋板块向鄂霍次克海板块俯冲的地区,并没有历史记载表明这里曾发生过矩震级9级的大地震,唯一可能的例外是公元869年贞观(Jogan)大地震[1],其震级到现在...     

Slip-Weakening Models of the 2011 Tohoku-Oki Earthquake and Constraints on Stress Drop and Fracture Energy

We present 2D dynamic rupture models of the 2011 Tohoku-Oki earthquake based on linear slip-weakening friction. We use different types of available observations to constrain our model parameters. The distribution of stress drop is determined by the final slip distribution from slip inversions. As three groups of along-dip slip distribution are suggested by different slip inversions, we present three slip-weakening models. In each model, we assume uniform critical slip distance eastward from the hypocenter, but several asperities with smaller critical slip distance westward from the hypocenter. The values of critical slip distance are constrained by the ratio of deep to shallow high-frequency slip-rate power inferred from back projection source imaging. Our slip-weakening models are consistent with the final slip, slip rate, rupture velocity and high-frequency power ratio inferred for this earthquake. The average static stress drop calculated from the models is in the range of 4.5–7 MPa, though large spatial variations of static stress drop exist. To prevent high-frequency radiation in the region eastward from the hypocenter, the fracture energy needed there is in the order of 10 MJ/m², and the average up-dip rupture speed cannot exceed 2 km/s. The radiation efficiency calculated from our models is higher than that inferred from seismic data, suggesting the role of additional dissipation processes. We find that the structure of the subduction wedge contributes significantly to the up-dip rupture propagation and the resulting large slip at shallow depth.     

Review of the Research of the 2011 MW9.0 Tohoku-Oki Earthquake

On March 11, 2011, a M_W9.0 earthquake occurred in the Japan Trench, causing tremendous casualties, and attracting extensive concern. Based on the results of related research, this paper analyzes the observations, phenomena and understandings of the earthquake from varied aspects, and obtains four main conclusions. (1) The earthquake, occurring in the subduction zone in the Japan Trench located in the northwest boundary of the pacific plate has two zones of concentrated coseismic slip at different depths, and the slip in the deep zone is relatively small. Though there have been many M7.0 historical earthquakes, slips in the shallow zone are large, but there have been few historical strong earthquakes. (2) Constrained by GPS data, the study of fault movement shows that fault movement in the Japan Trench has a background of widely distributed stability and locking (the locking zone is equivalent that of coseismic rupture zone). Perturbation occurred after the 2008 M8.0 Hokkaido earthquake, several M7.0 events had after slips larger than the coseismic slip, and two obvious slow slip events were recorded in 2008 and 2011. Eventually, the March 9, 2011 M7.0 foreshock and the March 11, 2011 M_W9.0 mainshock occurred. The pre-earthquake changing of the fault movement in the Japan Trench is quite clear. (3) Traditional precursory observation show no obvious anomaly, possibly due to monitoring reason. Anomaly before earthquake consists of high stress state in focal zone reflected by some seismic activity parameters, short period anomaly in regional ground motion, etc. (4) The analysis of physical property in focal zone aroused more scientific issues, for example, is there obvious difference between physical property in focal zone and its vicinity? Does frictional property of fault determine seismogenic ability and rupture process? Whether pre-earthquake fault movement include pre-slips? Could deep fluid affect fault movement in focal zone? Experience is the best teacher, and authors hope this paper could be a modest spur to induce others in basic research in earthquake forecast and prediction.     

2011年日本东北地区近海大地震最大位移区的应力状态

2011年MW9.0的日本东北地区近海大地震在日本海沟附近产生了最大超过50m的同震滑动,这可能导致该区域进入了应力近乎完全卸载的状态。为检验上述假设,我们利用综合大洋钻探计划在震后1年左右时间得到的钻孔资料,以及对震前应力状态的推断结果,确定了陆楔前缘的原地应力状态。根据孔壁崩落法对水平应力方向与应力大小的估算结果,以及在破裂向海沟轴方向传播过程中同震位移增大的现象,可推断原地水平应力在地震过程中变小了。这种应力变化显示出板块界面的前缘发生了相对自由的滑动,该结论与地震引起了断层弱化和发生了近乎全部的应力降等观测结果是吻合的。     

Seafloor geodetic network establishment and key technologies

Seafloor geodetic network construction involves the development of geodetic station shelter, network configuration design, location selection and layout, surveying strategy, observation model establishment and optimization, data processing strategy and so on. This paper tries to present main technological problems involved in the seafloor geodetic network construction, and seek the technically feasible solutions. Basic conceptions of developing seafloor geodetic station shelters for shallow sea and deep-sea are described respectively. The overall criteria of seafloor geodetic network construction for submarine navigation and those of network design for crustal motion monitoring are both proposed. In order to enhance application performances of the seafloor geodetic network, the seafloor network configuration should prefer a symmetrical network structure. The sea surface tracking line measurements for determining the seafloor geodetic station position should also adopt an approximately symmetrical configuration, and we recommend circle tracking line observations combined with cross-shaped line(or double cross-shape line) observations for the seafloor positioning mode. As to the offset correction between the Global Navigation Satellite System antenna phase center and the acoustic transducer, it is recommended to combine the calibration through external measurements and model parameter estimation. Besides, it is suggested to correct the sound speed error with a combination of observation value correction and parameterized model correction, and to mainly use the model correction to reduce the influence of acoustic ray error on the seafloor positioning. Following the proposed basic designs, experiments are performed in shallow sea area and deep-sea area respectively. Based on the developed seafloor geodetic shelter and sufficient verification in the shallow sea experiment, a long-term seafloor geodetic station in the deep-sea area of 3000 m depth was established for the first time, and the preliminary positioning result shows that the internal precision of this station is better than 5 cm.     

2011年东日本大地震对中国东北地区地壳变形的持续影响

实际监测与理论研究均表明, 2011年东日本大地震对中国东北地区地壳变形产生了显著的影响。 然而, 对震后影响的持续时间却存在着不同的看法, 有的学者依据所掌握的资料认为即使受影响最大的东北地区其持续时间截止于2013年上半年等。 针对这一问题, 本文对中国地壳运动观测网络2016年前的多期复测资料进行了处理与分析, 结果表明: 这种影响不仅存在, 而且仍是目前东北地区时段性地壳变形的主体成分, 仅是随着时间的推移呈现逐渐减弱的态势。 因此, 推测震后影响至少在今后数年内仍将存在。     

Joint analysis of b-value and apparent stress before the 2011 MW9.0 Tohoku-Oki,Japan earthquake

Detecting tempo-spatial changes of crust stress associated with major earthquakes has implications for understanding earthquake seismogenic processes. We conducted a joint analysis of b-value and apparent stress in the source region before the March 11, 2011 M_W9.0 Tohoku-Oki, Japan earthquake. Earthquakes that occurred between January 1, 2000 and March 8, 2011 were used to estimate b-values, while source parameters of events with magnitudes of Ms5.0–6.9 between January 1, 1997 and March 8, 2011 were used to calculate the apparent stresses. Our results show that the average b-value decreased steadily from 1.26 in 2003 to 0.99 before the Tohoku-Oki mainshock. This b-value decrease coincided with an increase in the apparent stress from 0.65 MPa to 1.64 MPa. Our results reveal a clear negative correlation between the decrease in b-value and increase in apparent stress, which lasted for approximately eight years prior to the 2011 mainshock. Additionally, spatial pattern results of the relative change in b-value show that the area associated with drastic b-value decreases (25% or greater) was concentrated near the 2011 mainshock epicenter. The joint analysis of b-value and apparent stress provides a promising method for detecting anomalies that could serve as potential indicators of large earthquakes.     

利用 Yabuki＆Matsu'ura 反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震同震滑动分布

Yabuki & Matsu'ura反演方法是利用ABIC最佳模型参数选取方法和平滑的滑动分布作为约束条件,由形变观测数据计算发震断层滑动分布.本文基于日本列岛同震GPS观测数据和发震断层曲面构造模型,利用Yabuki&Matsu'ura反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震的发震断层同震滑动分布.反演结果表明,断层面上的最大滑动量为35 m,较大滑动分布在浅于30 km的震源中心上部,最大破裂集中在20 km深度的地方,其地震矩约为3.63×1022N·m,对应的矩震级为Mw9.0.模拟结果显示Yabuki&Matsu'ura反演方法更适用于倾角低于40°的断层模型反演.最后,本文基于上述方法获得的发震断层滑动模型,利用地球体位错理论正演计算该地震在中国及其邻区产生的远场形变,正演计算结果基本可以解释由中国GPS陆态网络观测到的同震形变.     

On the stochastic relations between the points of a geodetic network

Summary Three important characteristics of the stochastic relation between the position of one point of a geodetic network and the positions of the remaining points are investigated: a) Change of the ellipse of confidence of the position of a point of the network to the ellipse of confidence of the position of a point of the network, relative to any other point of the network. (This characteristic is called the change of the global ellipse of confidence to the relative ellipse). b) The regression between the positional vector of the investigated point and the positional vectors of some group of other points of the network and c) the so-called total correlation between any two points of the network. The latter characteristic is a question of determining a characteristic direction for each pair of investigated points, which would have the property that the correlation coefficient between the random displacement of the first point in its characteristic direction and the random displacement of the second point in its characteristic direction is maximum.     

利用同震和震后位移数据联合反演2011年日本MW9.0地震同震断层滑动

基于黏弹性球体位错理论,联合陆地和海底同震GPS数据以及日本本岛330个陆地GPS站点5～10年的震后数据,反演了日本Mw9.0地震的断层滑动模型,提升了断层滑动分布在细节上的合理性.首先,基于日本本岛330个陆地GPS站点震前2年和震后10年的连续观测数据,获取了日本Mw9.0地震震后5～10年的年平均位移,该时段的...     

超导重力仪检测2011年日本东北Mw9.0级地震前的重力扰动信号

对大地震前的扰动现象的研究有助于认识地震孕育的动力学过程,震前重力扰动已成为关注的热点之一. 对大地震前的扰动现象的研究有助于认识地震孕育的动力学过程,震前重力扰动已成为关注的热点之一. 为检验日本M_w9.0级地震是否存在震前扰动现象,本研究利用全球超导重力仪记录到的地震前后7天内20组秒采样数据进行分析. 经潮汐、大气改正等处理去除仪器的漂移及残余潮汐效应,得到非潮汐重力变化曲线.结果表明大部分振幅大于30×10^-8 m·s^-2的曲线反映了全球M_w≥6级地震引起的高频波动信号,其中11组数据在3月9日M_w7.3级前震之前出现了扰动现象.震前扰动可分解为三个频段,其中,低于0.1 Hz和高于0.18 Hz的分量分别反映了地震波动信号及非构造信息,中间频段(0.118~0.18 Hz)信号能够较大程度地压制地震波动信号、并同时保留异常扰动信息.它的振幅在3月7日10时之前基本保持约1×10^-8 m·s^-2,之后开始逐渐增大,到3月9日7.3级前震前后达到最大,此后振荡衰减,振幅保持约(5~10)×10^-8 m·s^-2,直至主震发生.中间频段信号的变化特征与主震前的应力迁移过程以及实验记录到的地震成核过程有许多相似之处;不过,震前重力异常是否与主震前的应力加速积累有关,仍待进一步研究.     

2011年M_W9.0东日本大地震动力学破裂过程的数值模拟

力学上,地震可以看作在应力场作用下由于断层带介质的突然损伤或软化导致的断层带失稳事件.本文基于这个地震动力学模型,利用一种可以模拟断层大位错的有限元方法,研究了2011年MW9.0东日本大地震(Tohoku-Oki)的动力学破裂过程.比较了无障碍体和具有不同刚度障碍体的断层带模型产生的断层位移、位错和应力降.主要结果表明,障碍体的存在并不明显地改变障碍体区域的初始构造应力场.对有障碍体情形,准静态结果显示断层上盘最大逆冲位移和最大剪切位错分别为51m和58m,均发生在海底表面海沟处,与无障碍体的结果(最大剪切位错约55m)相比差别不大;下盘最大倾向位移(-10m)并不与上盘最大值出现在同一位置,而是在障碍体处.障碍体处剪应力降(约11 MPa)大于周围非障碍体区域.障碍体处正应力降的最大值约为3 MPa.模拟结果似乎不支持海山是导致本次地震异乎寻常大位错的原因,而倾向于断层带剪切刚度在地震过程中极度损伤或软化.     

利用Yabuki & Matsu'ura反演方法计算 2011年日本东北地区太平洋海域Mw9.0级地震同震滑动分布

Yabuki & Matsu'ura反演方法是利用ABIC最佳模型参数选取方法和平滑的滑动分布作为约束条件,由形变观测数据计算发震断层滑动分布.本文基于日本列岛同震GPS观测数据和发震断层曲面构造模型,利用Yabuki & Matsu'ura反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震的发震断层同震滑动分布.反演结果表明,断层面上的最大滑动量为35 m,较大滑动分布在浅于30 km的震源中心上部,最大破裂集中在20 km深度的地方.其地震矩约为3.63×1022N·m,对应的矩震级为Mw9.0.模拟结果显示Yabuki & Matsu'ura反演方法更适用于倾角低于40°的断层模型反演.最后,本文基于上述方法获得的发震断层滑动模型,利用地球体位错理论正演计算该地震在中国及其邻区产生的远场形变,正演计算结果基本可以解释由中国GPS陆态网络观测到的同震形变.     

反馈网络对地震计性能的影响

通过对地震计数学模型的详细分析，阐明了不同形式的反馈网络在地震计中的不同作用，结果表明：①机械部分对反馈式地震计频率特性的影响在一定条件下可以忽略。②地震计的类型和通频带由反馈网络决定。③反馈扩大地震计的动态范围，降低了地震计的分辨率。④反馈系统对地震计的信噪比没有影响。