基于融合重力异常的大别造山带东段莫霍面反演研究

利用实测布格异常和EGM2008重力异常融合结果,采用频率域位场反演方法计算大别造山带东段莫霍面三维空间分布,结合区域地质构造和地震活动等资料讨论大别造山带东段莫霍面分布特征及其构造含义等。研究结果显示:①莫霍面空间分布体现了块体构造差异,大别造山带莫霍面最深,最大深度达42km,显示东大别造山带存在明显山根,扬子地块深度次之,华北地块最浅;郯庐断裂带及其东侧区域存在莫霍面上隆带;②大别造山带北侧和南侧莫霍面陡变带分别位于青山—晓天断裂附近和襄樊—广济断裂以北,体现华北地块和扬子地块向大别造山带之下俯冲的构造特征,指示了深部构造缝合带位置;莫霍面深度特征表明郯庐断裂带区域构造边界带属性明显,其切割深度至少达到壳幔边界,大别造山带商城—麻城断裂两侧经历了差异隆升;③研究区域绝大多数地震发生在莫霍面以上深度,莫霍面深度陡变带、上隆带及不同莫霍面深度特征区的转换带为区域地震活动提供了深部构造条件。     

安徽霍山地震区深部电性结构和发震构造特征

霍山地震区位于大别造山带北缘华北板块与扬子板块接触带上,是大别造山带及周边地震活动最频繁、最集中的地区.83个大地电磁测点组成的大地电磁三维阵列覆盖了整个霍山地震区.用多重网格法、印模迭代重构法和非线性共轭梯度法对阵列数据进行三维带地形反演,获得了地震区深部三维电性结构.电性结构显示,北大别、北淮阳区的中上地壳为电阻率1000Ωm以上的高阻区,中下地壳为电阻率数十欧姆米的相对低阻区;六安盆地电阻率整体较低,中地壳存在显著的电阻率为几欧姆米的壳内高导层.北西向的晓天—磨子潭断裂分隔了北大别高阻层和北淮阳高阻层,在浅部向NE倾,深部向SW倾;北东向的落儿岭—土地岭断裂切穿北大别上地壳高阻层.小震双差定位结果表明,地震主要发生在NE向延伸的落儿岭—土地岭断裂附近的北大别、北淮阳中上地壳的高阻区,并集中于NW向的晓天—磨子潭断裂运动所造成的构造薄弱带中;2014年M S4.3霍山地震震源深度较深,位于北大别高阻区内部的电性梯度较大的区域.综合上述结果我们认为,霍山地震区的主要发震断裂为落儿岭—土地岭断裂,断裂的运动变形充分利用了晓天—磨子潭断裂早先活动所形成的构造薄弱带,断裂下方壳源高导体中的流体沿断层传播使断层强度弱化,使得这些薄弱带区易于发生小地震.由于北大别、北淮阳构造区显著高阻层的存在,我们认为霍山地震区存在发生6级以上中强震的深部孕震环境.     

江苏及邻区地壳上地幔结构研究

根据江苏及邻区大地构造单元和人工地震研究程度的不同,划分为下扬子地区、郯庐断裂带、大别山东段、苏鲁地块等地区,介绍了利用人工地震、天然地震层析成像及其它地球物理探测等方法对这一地区的深部构造研究结果。其中下扬子地区莫霍面的深度为28~33 km,横向不均匀,由西向东逐渐减薄;大别—苏鲁地区地表超高压变质带地壳波速特征主要表现在上地壳顶部的高速层,其厚度一般小于10 km;苏鲁地块地壳厚度约在32~33 km之间,明显高于周边地区;大别山地区沿造山带方向莫霍面变化平缓,地壳厚度33 km左右;横穿造山带方向起伏变化较大,莫霍面最深达41 km左右。天然层析成像研究范围大,分辨率较低。大别苏鲁地区人工地震测深工作较深入,但调查范围较小,主要在浅层,认识相对局部,有待于面上的大范围深部结构图像研究。     

郯庐断裂带南段大地电磁主轴方位统计分析及其意义

本文利用多测点-多频点阻抗张量统计成像技术统计分析了三条横跨郯庐断裂带南段的大地电磁剖面上的电性主轴分布特征,从电性主轴的角度描绘了大别郯庐构造区不同构造的特征及构造间的几何关系.统计结果显示:郯庐断裂带南段在断裂带浅部的电性主轴方位及二维有效因子中有明显的显示,同时郯庐断裂带南段可能未影响断裂深部介质;扬子板块及大别块体东南部地区浅部和深部的主轴方位明显不一致,表明浅部和深部之间存在结构不整合,主轴方位表明不同深度的介质受不同方向构造力的作用;大别造山带南部及以南地区深部介质中EW向电性主轴可能是板块汇聚时期扬子板块深部韧性层受NS向挤压应力作用的结果,郯庐断裂带以东地区深部韧性层中NWW向的电性主轴则是扬子板块沿郯庐带向NNE方向运动形成的;北大别块体深浅一致的电性主轴分布特征符合中生代穹隆抬升、增厚的演化过程.对该地区电性主轴的统计结果为认识该地区的深部结构及地下介质的运动极性提供了重要的依据,也为认识大别郯庐构造区的演化过程提供了新的证据.     

郯庐断裂北带地区主要构造单元壳幔结构特征与地震活动性.

以东北地区的地学断面及其相关资料为基础,分析了郯庐断裂北带地区的地壳-上地幔结构特征与不同地壳构造单元地震活动性之间的关系．结果表明:郯庐断裂北带两侧的壳幔结构存在明显差异,断裂带本身的发育与壳幔结构变异带关系密切;断裂的活动特点及分段明显受到地壳介质力学特性影响;郯庐断裂北带地区软流层的起伏与地貌形态呈明显的镜像关系,前者还制约着不同地壳构造单元的构造活动性及变形方式;郯庐断裂北带地区南、北两端地震较为活跃地段均与壳内低速 高导层发育地段相对应的事实,证明了刚性地块中壳幔结构较为复杂的活动断裂段是构成中强地震的主要控震构造．      

郯庐断裂带南段重力异常及不同深度的横向构造特征

为研究郯庐断裂带南段及周边构造的深部特征、空间展布、交切关系,文中利用小波多尺度分析方法对该地区的布格重力场进行场源分离,剖析不同深度下的地壳横向构造;同时采用Parker变密度模型对莫霍面深度进行反演分析。研究表明,郯庐断裂带南段表现为NNE走向的大型重力区域场梯度条带,切割深度达岩石圈地幔,其两侧密度结构及构造特征差异明显。沉积层及上地壳密度结构复杂,郯庐带东支2条断裂形成线性异常纵贯区域;而西支2条断裂线性异常较弱,截切EW向的重力异常体断续延展。中、下地壳密度结构简单,断裂带形成宽缓的低异常条带,反映了白垩纪—古近纪伸展环境下造成的地堑式构造。西支2条主干断裂纵穿合肥市,沿肥东凹陷西缘向S延伸,由于舒城以南的高密度圈闭体遮挡,尖灭于舒城县一带;而肥中、六安—合肥、肥西—韩摆渡EW向的深大断裂交切于西支,未延伸至东支。该区近代小震多发生于断裂构造所对应的重力高、低异常转化带之间,尤其是断裂之间的交会处与郯庐构造带内部,郯庐断裂带南段为历史强震空区段,考虑到断裂带莫霍面陡变及与多条深大断裂交切等深部环境,中强震危险性不容忽视。     

大别─苏鲁碰撞造山带的地震层析成像研究──岩石圈三维速度结构

对包含大别－苏鲁碰撞造山带在内的东经 １１２°－１２４°,北纬２８°－３９°区域进行地 震层析成像研究,重建其地壳及上部地幔的三维速度图像．结果表明：造山带岩石圈速度横 向不均匀性显著;大别造山带以商城－麻城断裂为界,东侧的大别地块与西侧红安地块在地 壳速度上是两个不同的速度块体;中地壳 １５-２５ｋｍ深度范围内存在低速带,与伸展滑脱构造 有关;南、北大别构造单元之下,莫霍面下凹,地壳内发育了速度为６．５－６．６ｋｍ／ｓ、向北倾斜的 相对高速体,与超高压变质岩体相对应;在大别－苏鲁造山带下方的上部地幔中存在向北倾 斜的板片状高速体,结合已有地质、地球化学证据推测,它是三叠纪俯冲的扬子地块的残留 体;俯冲板片在深部发生了断离．本文利用地震层析成像方法揭示的造山带岩石圈速度结构 细节,对研究与地表地质有关的地球动力学无疑是十分重要的．     

落儿岭-土地岭断裂几何结构及晚第四纪活动特征初探

落儿岭-土地岭断裂是东大别地区重要的发震构造,然而限于自然地理条件等因素,前人对其的研究并不充分。本文以发生多次中强地震的落儿岭-土地岭断裂为研究对象,在高精度卫星影像解译的基础上,通过详细的野外地质地貌调查,尤其是对典型断层剖面进行分析,研究落儿岭-土地岭断裂的断错地质地貌特征、几何结构及活动特征。通过野外调查并结合其他资料分析认为落儿岭-土地岭断裂为发育于大别造山带内部的一条走向NE、向NW陡倾的断裂带。依据地质地貌特征及地震活动性,断裂可以划分为杨树沟-黑石渡段和黑石渡-横塘岗乡段两个几何段落。断裂最新活动时代为中更新世晚期-晚更新世早期,断裂最新活动继承了中生代以来的运动方式,以兼具右旋走滑的拉张正断为主。     

苏鲁造山带及邻区深部电性结构研究

苏鲁造山带及毗邻华北地块,是中国东部地学研究的热点区域,研究其深部结构可以为讨论苏鲁超高压变质带以及华北克拉通演化提供重要的证据．对横切苏鲁造山带获得的一条大地电磁测深剖面资料进行了解释,剖面沿SE129°,西起华北地块,跨郯庐断裂带、苏鲁超高压变质带、苏鲁高压变质带,止于扬子地块．大地电磁反演解释采用了二维非线性共扼梯度法,用TE和TM联合模式得到了关于测区地下150km以上的电性剖面．该电性剖面在横向上,沿剖面自西而东,划分出了7个电性分区,电性边界带与郯庐断裂带、海州一泗阳断裂以及嘉山一响水断裂等重要的边界断裂具有很好的对应关系;纵向上,划分出6个电性构造单元．发现了在华北地块与扬子地块内存在着（壳内的）高导区,而苏鲁造山带下部没有发现高导区,这一点与大别造山带存在较大的差异．发现了在50—90km之间层位,存在较连续的相对低阻带,推测为上地幔顶部的软弱带,在该低阻带下部分别对应华北地块上地幔浅部相对高阻区、苏鲁造山带上地幔浅部相对低阻区以及扬子地块上地幔浅部相对高阻区．从整个二维电性结构模型来看,在苏鲁造山带及邻区上地幔浅部不存在异常低的电阻率,这表明现今已不存在与岩石圈减薄有关的热软流圈物质．     

大别─苏鲁碰撞造山带的地震层析成像研究──岩石圈三维速度结构

对包含大别－苏鲁碰撞造山带在内的东经 １１２°－１２４°,北纬２８°－３９°区域进行地 震层析成像研究,重建其地壳及上部地幔的三维速度图像．结果表明：造山带岩石圈速度横 向不均匀性显著;大别造山带以商城－麻城断裂为界,东侧的大别地块与西侧红安地块在地 壳速度上是两个不同的速度块体;中地壳 １５—２５ｋｍ深度范围内存在低速带,与伸展滑脱构造 有关;南、北大别构造单元之下,莫霍面下凹,地壳内发育了速度为６．５－６．６ｋｍ／ｓ、向北倾斜的 相对高速体,与超高压变质岩体相对应;在大别－苏鲁造山带下方的上部地幔中存在向北倾 斜的板片状高速体,结合已有地质、地球化学证据推测,它是三叠纪俯冲的扬子地块的残留 体;俯冲板片在深部发生了断离．本文利用地震层析成像方法揭示的造山带岩石圈速度结构 细节,对研究与地表地质有关的地球动力学无疑是十分重要的．     

沂沭断裂带重力场及地壳结构特征

沂沭断裂带为郯庐断裂带山东段,新构造运动显著,是华北地区的强震活动带之一。文中收集了该地区的布格重力数据,利用小波多尺度分析方法对重力场进行有效分离,研究区域地壳结构特征及断裂空间展布,并应用Parker变密度模型对区域莫霍面进行反演分析,得到以下几点结论:1)重力区域场显示,沂沭断裂带形成了NNE走向的大型重力梯度带,分隔了鲁西、鲁东地块,成为区域内重要的地球物理分界线。2)重力局部场显示,中上地壳结构复杂,沂沭带内部呈现两堑一垒的重力异常格局,5条主干断裂形成线性梯度带分布于东、西地堑内,鲁西块体的多条NW向活动断裂交切于沂沭断裂带,多数断裂只交切于西地堑,而蒙山山前断裂和苍尼断裂横穿沂沭断裂带;下地壳结构相对简单,发生明显的褶曲构造,表现出大规模高、低密度异常相间排列的典型特征。3)区域莫霍面形态东高西低,沂沭断裂带形成了莫霍面陡变带,造成了东西分异格局,潍坊东—莒县—临沂一线出现莫霍面上隆区,具有强震发生的深部孕震环境。4)区域内地震多发于高、低重力异常转化带之间,特别是活动断裂对应的重力梯度条带之上,地震的发生与断裂活动有着密切的关系,沂沭断裂带地震活动性最强,且东地堑强于西地堑。     

大别山北缘深部结构的高精度重磁电震解析

大别山北缘位于大别造山带与华北陆块会聚地带,其深部精细结构还存在一些争议问题.根据近年来在研究区域内采集的地球物理数据,通过OMEGA、OASIS和FUGRO-LCT等多个重磁电震软件处理,获得了深部信息丰富的多种地球物理属性图像.本文从中选取了2条平行的电法和地震剖面,辅以高精度重磁数据,揭示了大别山北缘深部地质结构总体呈现华北陆块南向俯冲、北淮阳构造带向北逆冲的特征,并可清楚地识别大别山北缘造山带、肥南山前坳陷带和肥北斜坡带.其特征分别如下:①大别山北缘造山带,地震反射杂乱,电性高阻大于2000Ωm,以磨子潭-晓天断裂为界分为北大别杂岩带和北淮阳构造带,主要由大别群、卢镇关杂岩和佛子岭群组成,沿舒城-信阳断裂逆冲于中、下侏罗统之上,浅部发育北倾的金寨-龙门冲滑覆断裂.②肥南山前坳陷带,上部为中、新生界,呈现中低阻特征,地震反射界面清楚,断面波明显;下部为华北型前中生界沉积岩(上部可能残存古生界),电性特征为低阻(5~50 Ωm),地震反射较连续,内部反射弱.以六安断裂为界分为舒城凹陷和肥中断裂带,其北侧边界为肥中断裂.③肥北斜坡带,上部主要为侏罗系,呈现中低阻、水平层状地震反射特征;下部主要为霍邱群,电性为高阻,地震反射杂乱;中间为华北型前中生界沉积岩,向北逐渐减薄、直至缺失.上述认识将为华北陆块南向俯冲、南北板块碰撞缝合线和油气资源勘查等研究提供深部地质约束.     

太行山重力梯级带的密度结构及其地质解释

斜贯中国东部地区的大兴安岭 -太行山 -武陵山重力梯级带是一条地壳深部构造变异带。选择与其中段———太行山重力梯级带相垂直的阜平 -定州剖面 ,利用人工地震测深与布格重力异常资料 ,进行了地壳 -上地幔密度反演分析。结果表明 ,该剖面的地壳 -上地幔密度分布具有明显的不均匀性。总体特征是剖面西部密度较低 ,东部密度较高 ,地壳厚度的横向变化基本上决定了布格重力异常梯级带的形成。结合华北地区的中新生代构造演化 ,认为太行山重力梯级带是燕山期以后形成的 ,它与板块之间的相互作用及上地幔热物质的运动相关     

青藏高原东南缘多尺度重力场及构造含义

本文基于EGM2008重力场模型研究了青藏高原东南缘均衡重力异常和多尺度的布格重力异常特征,以鲁甸和景谷地震为例,认识其深部构造环境和动力学过程,为该区域的构造运动和地震孕育环境研究提供依据.结果表明,研究区布格重力异常和均衡重力异常与地质构造格局相关性较好,川滇地块剧烈的区域布格重力异常和非均衡状态与其强烈的地壳变形...     

TRANSVERSE STRUCTURES FEATURES OF DIFFERENT DEPTHS DERIVED FROM BOUGUER GRAVITY ANOMALIES IN THE SOUTHERN SEGMENT OF TAN-LU FAULT ZONE

To research the faults distribution and deep structures in the southern segment of Tan-Lu fault zone(TLFZ) and its adjacent area, this paper collects the Bouguer gravity data and makes separation by the multi-scale wavelet analysis method to analyze the crustal transverse structure of different depths. Meanwhile Moho interface is inversed by Parker variable density model. Research indicates that the southern segment of TLFZ behaves as a NNE-directed large-scale regional field gravity gradient zone, which separates the west North China-Dabie orogen block and the east Yangtze block, cutting the whole crust and lithosphere mantle. There are quite differences of density structures and tectonic features between both sides of this gradient belt. The sedimentary and upper crustal density structure is complex. The two east branches of TLFZ behave as linear gravity anomalous belt throughout the region, whereas the two west branches of TLFZ continue to extend after truncating the EW-trending gravity anomaly body. The lower crustal density structure is relatively simple. TLFZ behaves as a broad and gentle low abnormal belt, which reflects the Cretaceous-Paleogene extension environment caused graben structure. The two west branches of TLFZ, running through Hefei city, extend southward along the west margin of Feidong depression and pinch out in Shucheng area due to the high density trap occlusions in the south of Shucheng. The Feizhong Fault, Liu'an-Hefei Fault, and Feixi-Hanbaidu Fault intersect the two west branch faults of TLFZ without extending to the east. Recent epicenters are mainly located in conversion zones between the high-density and the low-density anomaly, especially in TLFZ and the junction of the faults, where earthquakes frequently occurred in the upper and middle crust. As strong earthquakes rarely occur in the southern segment of TLFZ, considering its deep feature of abrupt change of the Moho and intersections with many EW-trending faults, the hazard of strong earthquake cannot be ignored.     

Wavelet analysis and interpretation of gravity data in Sichuan-Yunnan region, China

Introduction Sichuan-Yunnan region, located in the east margin of Qinghai-Xizang (Tibetean) Plateau, is a transitional zone between the rapidly upheaving Tibetean Plateau and relatively steady Yangtze Platform. Under the pressure exerted by the northward movement of Indian Plate, Sichuan-Yunnan region has been undergone strong deformation and regmagenesis, becoming one of the regions with the most intensive seismicity in the world. The research on the tectonics and seismicity there is alw…     

THE SEISMOGENIC ENVIRONMENT ANALYSIS OF LUDIAN MS6.5 EARTHQUAKE USING GRAVITY DATA

The main rupture of Ludian M_S6.5 earthquake is directed to the northwest, which occurred in the east of Xianshuihe-Xiaojiang fault zone. The epicenter is in the transitional zone of the Sichuan-Yunnan block and the South China block, where there are many slip and nappe structures. Some controversy still remains on the earthquake tectonic environment. So, Bouguer gravity anomalies calculated by EGM2008 were broken down into 1-5 ranks using the way of Discrete Wavelet Transform(DWT), then we get the lateral heterogeneity in different depths of the crust. The distribution characteristics of Bouguer gravity anomaly are analyzed using measured gravity profile data. We also get its normalized full gradient(NFG)picture, and study the differences between different depths in crust. The results show that: (1)the characteristic of Buoguer gravity anomaly in southwest to northeast is high-low-high between the Lianfeng Fault(LFF)and Zhaotong-Ludian Fault(ZLF). The mainshock and aftershocks are distributed in the middle of the low-value zone, which means that the east moving materials of Qinghai-Tibet plateau broke through the southern section of Lianfeng Fault(LFF), moving along the Baogunao-Xiaohe zone(low-value belt)to the southeast, stopped by the Zhaotong-Ludian Fault(ZLF), and then earthquake occurred.(2)The third-order discrete wavelet transform(DWT)details show that: there is a good consistency between the negative gravity anomaly in upper crust and the distribution of major faults, which reflects that the rupture caused by the movements of the faults in crust has reduced gravity anomaly. There is a NW-trending negative anomaly belt near the epicenter, which may has some relationship to the southward development of the Daliangshan Fault(DLSF). So we speculate that the southward movement of Daliangshan Fault is the main direct force source of Ludian earthquake.(3)In the picture of the fourth-order DWT details, there is an obvious positive gravity anomaly under the epicenter of Ludian earthquake, which confirms the presence of a high-density body in the middle crust. While the fifth-order DWT details show that: A positive anomaly belt is below the epicenter too, which may be caused by mantle material intruding to the lower crust. Tensile force in crust caused by mantle uplift and extrusion-torsion force caused by Indian plate push are the main force source in the tensile and strike slip movement of the Ludian earthquake.(4)The normalized total gradient of Bouguer gravity anomalies of Huili-Ludian-Zhaotong profile shows that: there is obvious ‘deformation’ in the Xiaojiang fault zone which dips to the east and controls the local crust movement. There is a local ‘constant body’ at the bottom of the epicenter. The stable constant body in density has limiting effects to the earthquake rupture, which is the reason that the earthquake rupture' scale in strike and in depth are limited.(5)The ability of earthquake preparation in Zhaotong-Ludian Fault is lower than the Xianshuihe-Xiaojiang fault zone, and the maximum earthquake capacity in this area should be around magnitude 7.     

龙门山断裂带地壳密度结构

研究龙门山及邻区地壳密度结构对于认识该地区地震活动性具有重要意义.根据龙门山及邻区( 100°～105°E,28°～33°N)的布格重力异常资料,选取了跨越龙门山断裂带的6条重力测线,在深地震测深资料约束下,使用Geosoft软件分别反演出了龙门山地区地下的沉积层、康拉德界面和莫霍面的深度分布.研究结果表明:龙门山断裂带两侧的地壳结构明显不同,西面高原地区沉积层较薄,大部分为基岩出露;而东边盆地沉积层明显较厚,多在6km以上.莫霍面和康拉德面在两侧均相对平缓,康拉德面从东部的大约24km增加到青藏高原山区的35km左右;莫霍面深度从东部盆地的大约42km增加到西部青藏高原的67km左右.龙门山断裂带整体表现为一条近SN向的陡变重力梯度带,并在其地壳内各界面均发生错断,莫霍面和康拉德面错断距离分别达6 ～ 7km和3～ 5km.该区地壳的这种陡变和不均匀性是导致地震活动性强烈的主要原因之一.     

Deep tectonic setting of the 2008 Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 earthquake in southwestern China

Teleseismic P-wave receiver functions at 20 broadband seismic stations in the Longmenshan fault zone (LMFZ) and its vicinity were extracted, and the crustal thickness and the P- and S-wave velocity ratio were calculated by use of the H-k stacking algorithm. With the results as constraints, the S-wave velocity structures beneath each station were determined by the inversion of receiver functions. The crustal structure of the Rear-range zone is similar to that of the Songpan-Garze Block, whereas the velocity structure of the Fore-range zone resembles that of Sichuan Basin, implying that the Central Principal Fault of LMFZ is the boundary between the eastern Tibetan Plateau and the Yangtze Block. Lower velocity zone exists in lower crust of the Songpan-Garze Block and the central-southern segment of the Rear-range zone, which facilitates the detachment of the material in upper and middle crust. Joint analysis of the receiver functions and the Bouguer gravity anomalies supports the thesis on the detachment-thrust mode of the LMFZ. A double-detachment pattern is suggested to the tectonic setting in the Songpan-Garze Block. The upper detachment occurs at the depth of 10-15 km, and represents a high-temperature ductile shear zone. There is a lower detachment at the depth of about 30 km, below which the lower crust flow exists in the eastern Tibetan Plateau. Interpretation of the Bouguer gravity anomalies indicates that the Sichuan Basin is of higher density in upper and middle crust in comparison with that of the Songpan-Garze Block. The LMFZ with higher density is the result from the thrusting of the Songpan-Garze Block over the Sichuan Basin. In the lower crust, higher P velocity and higher density in the Sichuan Basin are related to more rigid material, while lower S velocity and lower density in the Songpan-Garze Block are related to the softened and weakened material. The higher density block beneath the Sichuan Basin obstructs the eastward flow of lower crustal material from the Tibetan Plateau, which is driven by the compression of northward movement of Indian Plate. The eastward movement of upper and middle crustal material is also obstructed by the rigid Yangtze Block, resulting in the stress concentrated and accumulated along the LMFZ. When the stress releases sharply, the Wenchuan M _s8.0 earthquake occurs. Supported by the National Natural Science Foundation of China (Grant Nos. 40334041, 40774037) and Joint Foundation of Earthquake Science (Grant No. 1040062)