摘要
地震层析成像方法是研究地球内部速度结构最有效的手段之一,它能获取地球内部不同深度不同尺度的三维速度图像,为人们深入研究地球内部精细结构提供可靠的信息。华北地区不仅是我国的经济政治中心,也是中国大陆地震活动最为强烈的地区之一,对该地区地下精细结构进行深入研究具有重要意义。
为进一步深入研究华北地区地下介质结构,中国地震局地球物理研究所于2006年在华北地区开展了大规模流动地震台阵观测。本文利用华北地震科学台阵和首都圈台网记录的近震和远震P波到时资料,采用近震和远震联合层析成像方法得到了华北地区(111°E-120°E,37°N-42°N)深至700km的高分辨率(水平方向约36km)的地壳上地幔三维P波速度结构,对研究区地壳上地幔横向不均匀性、强震发生的深部构造环境和区内重要的断裂带的深部速度结构特征等进行了研究和讨论。
本文研究获得的主要认识如下:
(1)华北地区不同构造块体速度结构存在明显差异。华北盆地下方在60-120km深度存在明显低速异常,可能与软流圈有关;在150-350km深度存在明显高速异常,这可能意味着华北盆地岩石圈遭受破坏后仍有古老的克拉通岩石圈残留体存在于上地幔中。渤海盆地下方地幔转换带存在一个较大的低速体,该低速体位于太平洋俯冲板块前缘的西侧,可能与太平洋俯冲板块的相变脱水等有关。山西裂陷盆地一带下方在60-350km存在明显的低速异常,可能与该区上地幔高温低速物质上涌有关。燕山隆起下方60-150km深度表现为连续的高速异常,我们推测是由该地区较大的岩石圈厚度引起的。
(2)研究区三维P波速度结构与强震分布对应关系的研究表明多数强震(M≧6.0)往往发生在高、低速异常过渡带上,且震源下方多表现为P波低速异常;强震的发生可能与地壳内存在的流体、壳幔过渡带的不均匀性和上地幔热物质作用等动力学背景密切相关。
(3)区内重要的断裂带(太行山山前断裂和张家口-蓬莱断裂带)处在研究区不同块体的交接带上,太行山山前断裂带在地壳中位于高低速异常的分界线附近,至壳幔过渡带这种界线变得很模糊;张家口-蓬莱断裂带在地壳中速度结构较复杂,在60-90km深度位于高低速异常的分界线附近,即位于燕山隆起区较厚的岩石圈块体边缘,表明其深部结构主要受不同块体岩石圈厚度变化的控制。
Seismic tomography has become a powerful tool for studying the three-dimensional crust and mantle structure.It can determine the different scales and the different depths of the three-dimensional velocity image,which provides reliable information for people studying the fine structure of the Earth's interior.North China is our country's economic and political center,and the earthquake activity of this area is one of the strongest area in China.Therefore,in-depth studying the fine structure of this area has great significance.
In order to in-depth study the fine structure of the North China,Institute of Geophysics China Earthquake Administration launched a large-scale movement of seismic array observations in North China in 2006.We determined a high-resolution (horizontal resolution is about 36km) 3D P-wave velocity model of the crust and upper mantle down to the depth of 700km under North China by jointly inverting local events and teleseismic events recorded by North China Seismic Array and Capital Seismic Network.We detailed study the lateral heterogeneity of the crust and upper mantle,the deep tectonic environment of strong earthquake's occurrence and the deep structure of important fault in North China.
The main results are as follows:
(1) Different structure blocks in North China have the existence of significant differences in velocity structure.Under North China basin,a pronounced low-velocity anomaly exists at depth of 60-120km,which may relate with shallow depth of asthenosphere,and a pronounced high-velocity anomaly exists at depth of 150-350km,suggesting that old lithosphere residues still exists in the upper mantle through large-scale lithosphere thinning.A large pronounced low-velocity anomaly exists in the mantle transition zone under Bohai sea basin and locates on the west side of the front subducting Pacific slab,which may result from phase transition and dehydration of the subducting Pacific slab.At the depth of 60-350km,a pronounced low-velocity anomaly exists under Shanxi rift basin,which may result form the upwelling of high temperature and low-velocity mantle material under the studying zone.At the depth of 60-150km,a pronounced high-velocity anomaly exists under Yanshan uplift,which may relate with the thick lithosphere.
(2)Through studying the relationship between the 3D P-wave velocity structure and strong earthquake distribution,we found that most strong earthquake(M≧6.0) often occurred between high-velocity and low-velocity anomaly transition zone.Low P-wave velocity anomaly exists under the hypocenter.The occurrence of strong earthquakes has the nearly relationship with the existence of fluid in the crust, heterogeneity in the crust-mantle transition zone and the upwelling mantle material.
(3)The Taihangshan front fault and the Zhangjiakou-Penglai fault become flimsy tectonic boundaries.The Taihangshan front fault becomes the boundary between high-velocity and low-velocity anomaly in the crust,but the boundary becomes blurred in the crust-mantle transitional zone.The velocity structure of Zhangjiakou-Penglai fault is very complex.At the depth of 60-90km,the Zhangjiakou-Penglai fault also becomes the boundary between high-velocity and low-velocity anomaly,and this fault locates on the edge of the Yanshan uplift which has thick lithosphere,suggesting that the deep structure of this fault may be controlled by listhosphere thickness of the different block.
引文
陈国英,宋仲和,安昌强,等.1991.华北地区地壳上地幔结构.地球物理学报,34(2):172-181
蔡学林,朱介寿,曹家敏,等.2006.中国及邻近陆海地区软流圈三维结构及其与岩石圈的相互作用.中幽地质,33(4):804-815
丁志峰.1999.近震层析成像的理论及应用.[博士学位论文].北京:中国地震局地球物理研究所
邓晋福,魏文博,等.2007.中国华北地区岩石圈三维结构及演化.北京:地质出版社,31-32
邓晋福,苏尚国,刘翠,等.2006.关于华北克拉通燕山期岩石圈减薄的机制与过程的讨论:是拆沉,还是热侵蚀和化学交代?地学前缘,13(2):105-119
邓晋福,赵海玲,莫宣学等著.1996.中国大陆根-柱构造-大陆动力学的钥匙.第一版.北京:地质出版社,26-27
邓启东,等,2004.中国活动断层分布图(1:400万),北京:科学出版社
段永红,张先康,放盛明.2002.华北地区上部地壳结构的三维有限差分层析成像.地球物理学报,45(3):362-369
房立华,吴建平.2009.华北地区基于噪声的瑞利面波群速度层析成像.地球物理学报,52(3):663-671
国家地震局.1989.中国岩石圈动力学地图集.北京:中国地图出版社
高战武,徐杰,宋长青,等.2001.张家口-蓬莱断裂带的分段特征.华北地震科学,19(1):35-54
何正勤,丁志峰,叶太兰,等.2002.中国大陆及其邻域的瑞利波群速度分布图与地壳上地幔速度结构.地震学报,24(3):252-259
何正勤,叶太兰,丁志峰.2009.华北东北部的面波相速度层析成像研究.地球物理学报,52(5):1233-1242
洪学海,朱介寿,曹家敏,等.2003.中国大陆地壳上地幔S波品质因子三维层析成像.地球物理学报,46(5):642-651
侯贵廷,叶良新,杜庆娥.1999.渤张断裂带的构造机制及其地质意义.地质科学,34(3):375-380
嘉世旭,张先康.2005a.华北不同构造块体地壳结构及其对比研究.地球物理学报,48(3):611-620
嘉世旭,齐诚,王夫运,等.2005b.首都圈地壳网格化三维结构.地球物理学报,48(6): 1316-1324
金安蜀,刘福田,孙永智.1980.北京地区地壳和上地幔的三维P波速度结构.地球物理学报,23(2):172-182
刘启元,李顺成,沈扬,等.1997.延怀盆地及其邻近地区地壳上地幔速度结构的宽频带地震台阵研究.地球物理学报,40(6):763-772
刘启元,王峻,陈九辉,等.2007.1976年唐山大地震的孕震环境:密集地震台阵观测得到的结果.地学前缘,14(6):205-213
路凤香,郑建平,邵济安等.2006.华北东北中生代晚期-新生代软流圈上涌与岩石圈减薄.地学前缘,13(2):86-92
刘福田,曲克信,吴华,等.1986.华北地区地震层面成像.地球物理学报,29(5):442-449
刘福田,曲克信,吴华,等.1989.中国大陆及其邻近地区的地震层析成像.地球物理学报,32(3):281-291
刘国栋,史书林,王宝钧.1984.华北地区壳内高导层及其与地壳构造活动性的关系.中国科学B辑,9:839-848
刘建华,刘福田,阎晓蔚,等.2004.华北地区Lg尾波衰减研究--Lg尾波Q_0地震成像.地球物理学报,47(6):1044-1052
刘建华,胥颐,郝天珧,等.2005.中国边缘海及邻区地壳和上地幔速度结构特征.地球物理学进展,20(3):594-599
李铁明,沈正康,徐杰,等.2007.华北地区Ms≧6.5级地震震源断层参数的研究.地球物理学进展,22(1):95-103
李乐,陈棋福,陈顒.2007.首都圈地震活动构造成因的小震精定位分析.地球物理学进展,22(1):24-34
李建华.1983.燕山南缘北西西向昌平-丰南隐伏活断层的初步研究.地震地质,5(4):61-66
李强,王椿镛,刘瑞丰,等.1999.应用层析成像技术研究华北地壳速度结构.地震地磁观测与研究,20(5):88-97
李志伟,胥颐,郝天珧,等.2006.环渤海地区的地震层析成像与地壳上地幔结构.地球物理学报,49(3):797-804
赖晓玲,张先康,孙译.2006.张北地震区壳幔边界复杂性特征及其构造意义.地震学报,28(3):230-237
赖晓玲,张先康,成双喜,等.2004.张家口-渤海断裂带西缘壳幔过渡带研究.地球物理学报,47(5):798-804
马杏垣,刘和普,王维囊,等.1983.中国东部中、新生代裂陷作用和伸展构造.地质学报,1:22-32
马杏垣.1987.中国岩石圈动力学纲要[M].北京:地质出版社,76
秦馨菱,Pedersen,L.B.,赵玉林,等.1991.唐山地震区地壳电性结构及MT探索潜在震源的可能性.地震学报,13(3):354-363
齐诚,赵大鹏,陈顒,等.2006.首都圈地区地壳P波和S波三维速度结构及其与大地震的关系.地球物理学报,49(3):805-815
邱瑞照,邓晋福,周肃,等.2004.华北地区三类岩石圈的壳幔岩石学结构与化学结构及其地球动力学意义.岩石矿物学杂志,23(2):127-140
邱瑞照,李廷栋,邓晋福,等.2006.中国大地构造单元新格局--从岩石圈角度的思考.中国地质,33(2):401-410
孙若昧,刘福田.1995.京津唐地区地壳结构与强震的发生-I.P波速度结构.地球物理学报,38(5):599-607
宋儒,姜振泉,王飞跃.2004.汾渭地堑地裂缝的控热现象及其机制分析.中国煤田地质,16(2):26-33
宋仲和,陈国英,安昌强,等.1993.中国大陆及其海域地壳上地幔三维速度结构.中国科学(B),23(2):180-188
滕吉文,曾融生,闫雅芬,等.2002.东亚大陆及周边海域Moho界面深度分布和基本构造格局.中国科学D辑,32(2):89-100
田有,赵大鹏,孙若昧,等.2007.1992年美国加州兰德斯地震--地壳结构不均匀性对地震发生的影响.地球物理学报,50(5):1488-1496
田晓峰,张先康,嘉世旭,等.2006.华北北部不同构造单元地壳速度模型-均匀分层介质反射波走时反演的结果.华北地震科学,24(3):1-5
汪素云,许忠淮,裴顺平.华北地区上地幔顶部Pn波速度结构及其构造含义.中国科学D辑,33(增刊):91-98
王未来,吴建平,房立华.2009.唐海-商都地震台阵剖面下方的地壳上地幔S波速度结构研究,地球物理学报,52(1):81-89
王若柏,顾国华,徐杰,等.2004.张家口-渤海地震构造带的地壳形变研究.地震地质,26(4):586-595
王帅军,张先康,方盛明,等.2008.渤海湾西北缘及其邻近地区地壳结构与构造特征.地球物理学报,51(5):1451-1458
王帅军,张先康,刘宝峰,等.2005.张渤地震构造带中西段及邻区深部构造探测.大地测量与地球动力学,25(3):110-115
汪一鹏,邓启东,朱世龙.1989.华北地区岩石圈动力学特征.马杏垣主编.中国岩石动力学地图集.北京:中国地图出版社,59
徐杰,宋长青,楚全芝.1998.张家口-蓬莱断裂带地震构造特征的初步探讨.地震地质,20(2):146-154
徐杰,高战武,宋长青等.2000.太行山山前断裂带的构造特征.地震地质,22(2):111-122
徐锡伟,吴卫民,张先康,等.2002.首都圈地区地壳最新构造变动与地震.北京:科学出版社,2-8
徐锡伟,邓起东.1988.晋北张性区盆岭构造及其形成的力学机制.中国地震,4(2):19-27
徐佩芳,刘福田,王清晨,等.2000.大别-苏鲁碰撞造山带的地震层析成像研究-岩石圈三维速度结构.地球物理学报,43(3):377-385
徐常芳.2003.中国大陆壳内与上地幔高导层成因及唐山地震机理研究.地学前缘,10(特刊):101-111
邢作云,赵斌,涂美义,等.2005.汾渭裂谷系与造山带耦合关系及其形成机制研究.地学前缘,12(2):247-262
邢作云,邢集善,赵斌.2006.华北地区深部构造特征.地质科技情报,25(6):17-36
邢集善,姚典群,黎明,等.1989.试从地球物理资料论山西地堑系的构造特征.山西地质,4(2):95-109
邢集善,叶志光,孙振国,等.1991.山西板内构造及其演化特征初探.山西地质,6(1):3-1
于湘伟,陈运泰,王培德.2003.京津唐地区中上地壳三维P波速度结构.地震学报,25(1):1-14
杨文采,李幼铭,等著.1993.应用地震层析成像.地质出版社
朱露培,曾融生,刘福田.1990.京津唐张地区地壳上地幔三维P波速度结构.地球物理学报,33(3):267-277
张先康,王椿镛,刘国栋,等.1996.延庆-怀来地区地壳细结构-利用深地震反射剖面.地球物理学报,39(3):356-364
张先康,祝志平,张成科,等.1998.张家口-渤海地震带及其两侧地壳上地幔构造与速度结构研究.见:国家地震局地质研究所编.活动断裂研究(6).北京:地震出版社.1-16
朱介寿,曹家敏,李显贵,等.1997.中国及其邻区地球三维结构初始模型的建立.地球物理学报,40(5):627-650
朱介寿,宣瑞卿,刘魁,等.2005.用瑞利面波研究东亚及西太平洋地壳上地幔三维结构.物探化探计算技术,27(3):286-193
朱介寿,曹家敏,蔡学林,等.2002.东亚及西太平洋边缘海高分辨率面波层析成像.地球物理学报,45(5):646-664
朱介寿,曹家敏,蔡学林,等.2004.欧亚大陆及西太平洋边缘海岩石圈结构.地球科学进展,19(3):387-391
朱介寿,曹家敏,蔡学林,等.2003.中国及邻近陆域海域地球内部三维结构及动力学研究.地球科学进展,18(4):497-503
张岭,刘劲松,郝天珧,等.2007.渤海湾盆地及其邻域地区地壳与上地幔层析成像,中国科学D辑,37(11):1444-1455
朱日祥.2007.地球内部结构探测研究-以华北克拉通为例.地球物理进展,22(4):1090-1100
Aki,K.,Lee,W.,1976,Determination of three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes,1.A homogeneous initial model.J.Geophys.Res.81:4381-4399
Aki,K.,Christoffersson,A.,Husebye,E.,1977,Determination of the three-dimensional seismic structure of the lithosphere.J.Geophys.Res.,82:277-296
Crosson R S.1976a,Crustal structure modeling of earthquake data,partl,Simultameous least squares estimation of hypocenter and velocity parameters.J.Geophys.Res.,81(17):3036-3046
Crosson R S.1976b,Crustal structure modeling of earthquake data,part2,Velocity structures of the Puget Sound region,Washington.J.Geophys.Res.,81(17):3047-3053
Carlson,R.,Pearson,D.,James,D.,et al.,2005,Physical,chemical,and chronological characteristics of continental mantle.Rev.Geophys.,43,dol:10.1029/2004RG000156
Chen,L.,Tao,W.,Zhao,L.,et al.,2008,Distinct lateral variation of lithospheric thickness in the Northeastern North China Craton.Earth Planet.Sci.Lett.,267:56-68
Dueker K.G.and Sheehan A.F.,1997,Mantle discontinuity structure of converted P to S waves across the Yellowstone hotspot track.J.Geophys.Res.,102(B4):8313-8328
Dziewonski,A.M.,and D.L.,Anderson,1984,Seismic tomography of the Earth's interior.Scientiffic American.72:483-494
Dziewonski,A.M.,Gilbert,F.,1976,The effect of small,aspherical perturbations on travel times and a re-examination of the corrections for ellipticity,Geophys.J.R.astr.Soc.,44:7-17
Eberhart-Phillips,D.,Christensen,D.H.,Brocher,T.M.et al.,2003,Imaging the transition from
Aleutian subduction to Yakutat collision in central Alaska, with local earthquakes and active source data, Eos Trans. AGU, 84(46), Fall Meet. Suppl., Abstr. S21C-01
Eberhart-Phillips, D., Christensen, D. H., Brocher T. M. et al., 2006, Imaging the transition from Aleutian subduction to Yakutat collision in central Alaska, with local earthquakes and active source data, J.Geop. Res., submitted
Grand, S. P., van der Hilst, R. D. and Widiyantoro, S., 1997, Global seismic tomography: a snapshot of convection in the Earth, GSA today, 7:1-7
Griffin, W., Zhang, A., O'reilly,S., et al. 1998. Phanerozoic evolution of the lithosphere beneth the Sino-Korean craton, in mnatle Dynamics and plate Interactions in East Asia, Geodyn. Ser., 27,edited by Flower, M. et al., AGU, Washington, D. C. 107-126
Geiger, I., 1912, probability method for the determination of earthquake epicenters from the arrival time only. Bull St. Louis Univ,8 (1):60-71
Gregory, P., et al., 2006, VP and VS structure of the Yellowstone hot spot from teleseismic tomography: Evidence for an upper mantle plume. J. Geophys. Res., 111, B04303, doi: 10.1029/2005JB003867
Huang, J., Zhao, D., 2006, High-resolution mantle tomography of China and surrounding regions.J. Geophys. Res. 111, B09305. doi:10.1029/2005JB004066
Humphreys, E. R., Clayton, R. W., 1988, Adaptation of back projection tomography to seismic travel time problems, J. Geophys. Res., 93:1073-1085
Huang, J., Zhao, D. 2009, Seismic imaging of the crust and upper mantle under Beijing and surrounding regions. Phys. Earth Planet. Inter.,173(3-4):330-348
Huang, J., Zhao, D. 2006, High-resolution mantle tomography of China and surrounding regions. J. Geophys. Res.,111,B0935.doi:10.1029/2005JB004066
Horiuchi, S., Ishii, H., Takagi, A., 1982a, Two-dimensional depth structure of the crust beneath the Tohoku district, the northeastern Japan arc, Part1. Method and Conrad discontinuity. J. Phys.Earth, 30:47-69
Horiuchi, S., Yamamoto, A.,Ueki, A., et al.,1982b, Two-dimensional depth structure of the crust beneath the Tohoku district, the northeastern Japan arc, Part II, Moho discontinuity and p wave velocity. J. Phys. Earth, 30:71-86
Kennett, B. L. N., Engdahl, E. R., 1991, Traveltimes for global earthquake location and phase identification, Geophys. J. Int., 105:429-465
Kissling, E., 1988, Geotomography with local earthquake data. Rev. Geophys., 26:659-698
Lei, J. Xie, F. Lan, C, et al., 2008, Seismic image under the Beijing region inferred from P and PmP data. Phys. Earth Planet. Inter., 162(3-4): 165-185
Lei, J., Zhao, D., 2005, P-wave tomography and origin of the Changbai intraplate volcano in Northeast Asia, Tectonophysics, 397:281-295
Langston, C., 1979, Structure under Mout Rainier, Washington, inferred form teleseismic body waves. J. Geophys. Res., 84(B9):4749-4762
Lippitsch, R., Kissling, E. and Ansorge, J., 2003, Upper mantle structure beneath the Alpine orogen from high-resolution teleseismic tomography, J. Geophys. Res., 108(B8), 2376, dor : 10.1029/2002JB002016
Owens, T.J., Zandt, G, and Taylor, S.R., 1984, Seismic evidence for an ancient rift beneath the
Cumberland Plateau, Tennessee: A detailed analysis of broadband teleseismic P waveforms. J.Geophys. Res., 89:7783-7795
Ozalaybey, S., M. K. Savage, A. F., Sheehan, et al.,1997, Shear wave velocity structure in the Northern Basin and Reange Province from the combined analysis of receiver function and surface waves. Bull. Seism. Soc. Am., 87:183-199
Pavlis, G., Booker, J. R.,1980, The mixed discrete-continuous inverse problem: Application of thes imultaneousd eterminationo f earthquakeh ypocenters and velocitys tructure. J .Geophys.Res.,85(B9):4801-4810
Paige, C. C., Saunders, M. A., 1982, LSQR: Sparse linear equations and least squares problems, AMC Transactions. Math., 8:195-209
Spencer, C., Gubbins, D., 1980, Travel-time inversion for simultaneous earthquake location and velocity structure determination in laterally varying media. Geophys. J. Roy. Astr. Soc, 63(1) :95-116
Thurber, C., Ellsworth, W., 1980, Rapid solution of ray tracing problems in heterogeneous media,Bull. Seism. Soc. Amer., 70:1137-1148
Thurber, C., 1983, Earthquake locations and threedimensional crustal structure in the Coyote Lake area, central California, J. Geophys. Res., 88:8226-8236
Thurber, C., 1993, Local earthquake tomography: velocities and Vp/Vs - theory, in Seismic
Tomography: Theory and Practice, edited by Iyer, H. M. and Hirahara, K., Chapman and Hall,London,563-583
Thurber, C., Eberhart-Phillips, D., 1999, Local earthquake tomography with flexible gridding,Computers and Geosciences, 25:809-818
Tian, Y., Zhao, D., Sun, R.,et al.,2009. Seismic imaging of the crust and upper mantle beneath the North China Craton. Phys. Earth Planet. Inter., 172(3-4):169-182
Um, J., Thurber, C, 1987, A fast algorithm for two-point seismic ray tracing. Bull. Seismol. Soc.Am., 77(3):972-986
VanDecar, J., Crosson, R..,1990, Determination of teleseismic relative phase arrivel times using multi-channel cross-correlation and least squares. Bull. Seismol. Soc. Am., 80(1): 150-169
Wu, F., Lin, J., Wilde, S., et al.,2005, Nature and significance of the Early Cretaceous giant igneous event in eastern China. Earth Planet. Sci. Lett., 233:103-119
Walck, M.C., Clayton, R.W., 1987, P wave velocity variations in the Coso Regio, California,derived form local earthquake travel times, J. Geophys. Res., 92,393-405
Waldhauser, F., Lippitsch, R., Kissling, E. et al., 2002, High-resolution teleseismic tomography of upper mantle structure using an a priori 3D crustal model. Geophys. J. Int., 150:141-403
Xu, Y. 2001, Thermo-tectonic destruction of the Archaean lithospheric keel beneath the Sino-Korean Craton in China: evidence, timing and mechanism. Phys. Chem. Earth, 26:747-757
Yuan X, Kind R, Mechie J, et al. 1997, Lithospheric and upper mantle structure of southern Tibet form a seismological passive source experiment. J. Geophys. Res.,102(B12):27491-27500
Zhao, D., Horiuchi, S., Hasegawa, A., 1990, 3-D seismic velocity structure of the crust and uppermost mantle in the northeastern Japan arc. Tectonophysics, 181:135-149
Zhao, D., Hasegawa, A., Horiuchi, S., 1992a, Tomographic imaging of P-and S-wave velocity structure beneath northeastern Japan. J. Geophys. Res.,97:19909-19928
Zhao, D., Hasegawa, A., Kanamori, H., 1994, Deep structure of Japan subduction zone as derived from local, regional, and teleseismic events. J. Geophys. Res., 99:22313-22329
Zhao, D., Negishi, H., 1998, The 1995 Kobe earthquake: seismic image of the source zone and its implications for the rupture nucleation. J. Geophys. Res.,103: 9967-9986
Zhao, D., Matsuzawa, T., Hasegawa, A., 1997a, Morphology of the subducting slab boundary in the northeastern Japan arc. Phys. Earth Planet. Int., 102:89-104
Zhao, D., Ochi, F., Hasegawa, A., Yamamoto, A., 2000a, Evidence for the location and cause of large crustal earthquakes in Japan. J. Geophys. Res., 105:13579-13594
Zhao, D., 2001b, Seismic structure and origin of hotspots and mantle plumes,Earth Planet. Sci. Lett., 92:251-265
Zhao, D., Mishra, O and Sanda, R., 2002, Influence of fluids and magma on earthquakes:seismological evidence. Phys. Earth Planet. Inter., 132:249-267
Zhao, D., 2004, Global tomographic images of mantle plumes and subducting slabs: insight into deep Earth dynamics. Phys. Earth Planet. Inter., 146:3-34
Zhao, D., Todo, S. and Lei, J., 2005, Local earthquake reflection tomography of the Landers aftershock area, Earth Planet. Sci. Lett., 235:623-631
Zhao, D., Lei, J., Inoue, T. et al., 2006, Deep structure and origin of the Baikal rift zone, Earth Planet. Sci. Lett., 243:681-691
Zhao, D., Wang, Z., Umino, N., Hasegawa, A., 2007, Tomographic imaging outside a seismic network: application to the Northeast Japan arc. Bull. Seismol. Soc. Am. 97:1121-1132