文章快速检索     高级检索
  大地测量与地球动力学  2017, Vol. 37 Issue (6): 633-637, 643  DOI: 10.14075/j.jgg.2017.06.016

引用本文  

黎哲君, 张毅, 杨光亮, 等. 安徽及周边地区岩石圈有效弹性厚度研究[J]. 大地测量与地球动力学, 2017, 37(6): 633-637, 643.
LI Zhejun, ZHANG Yi, YANG Guangliang, et al. Effective Elastic Thickness of Anhui and Its Surrounding Regions[J]. Journal of Geodesy and Geodynamics, 2017, 37(6): 633-637, 643.

项目来源

中国地震局地震研究所所长基金(IS201416141);安徽省地震科研基金青年项目(20150404)。

Foundation support

Director Fund of Institute of Seismology, CEA, No. IS201416141; Youth Project of Seismological Research Foundation of Anhui Province, No. 20150404.

第一作者简介

黎哲君,工程师,主要从事地球重力场观测和研究工作,E-mail:zhejunli@126.com

About the first author

LI Zhejun, engineer, majors in earth gravity field observation and research, E-mail:zhejunli@126.com.

文章历史

收稿日期:2016-05-08
安徽及周边地区岩石圈有效弹性厚度研究
黎哲君1     张毅1     杨光亮2     贾宇鹏3     肖伟鹏1     谈昕1     龙剑锋1     
1. 安徽省地震局,合肥市长江西路558号,230031;
2. 中国地震局地震研究所,武汉市洪山侧路40号,430071;
3. 中国地震局地球物理勘探中心,郑州市文化路75号,450002
摘要:基于EIGEN-6C4重力异常和ETOPO1地形数据,采用基于fan小波的相关性分析方法计算安徽及周边地区较高分辨率的各向同性岩石圈有效弹性厚度Te(effective elastic thickness),结合区域地质构造和地球动力学背景进行讨论。结果表明,研究区Te处于5~75 km之间,大别造山带和下扬子地区为Te低值区,郯庐断裂两侧Te值相对较低,河淮盆地和苏北盆地Te值较高, Te分布反映了区域特定地质构造特征和地球动力学背景;地震活动主要分布在10~50 km Te值范围内,集中在Te低值区和Te值陡变带,表明相关区域岩石圈抵抗形变能力较弱,易于发生地震。
关键词岩石圈有效弹性厚度安徽地震活动

开展岩石圈Te研究有利于清晰地认识与各构造块体形成和发展相关的消减作用和造山运动等地球动力学过程[1]。近10 a来,利用最新的计算方法,结合重力、地形数据,国内外众多学者开展了多项岩石圈Te研究工作。Kirby等[2]、郑勇等[3]、Mao等[4]和Chen等[5]采用重力和地形导纳或小波相关性分析等方法,计算澳大利亚和华北克拉通等多个地区的岩石圈Te,分析岩石圈Te与大地构造、地表热流、地震活动等的关系。

安徽及周边地区(图 1)处于华北地块、秦岭-大别造山带和扬子地块相互作用和影响的关键部位[6-7]。郯庐断裂带属岩石圈尺度的深大断裂构造系统,对中国东部区域构造、岩浆活动以及现代地震活动等具有控制作用,大别造山带附近区域具有造山带和盆地交汇、深大断裂横贯等复杂构造特征,较高程度的地球物理研究及频繁地震活动等一系列因素使得该区域成为开展岩石圈Te研究的理想场所。本文基于EIGEN-6C4[8]重力异常数据和ETOPO1[9]数字高程数据,采用荷载分离方法的小波版本计算安徽及周边地区(114°~120°E,29°~35°N)的岩石圈Te及对应荷载比,根据计算结果分析岩石圈Te及对应荷载比及其与区内各构造单元地质构造特征、地球动力学背景、地震活动之间的关系。

图 1 安徽及周边地区地质构造背景 Fig. 1 Geotectonic setting of Anhui and its surrounding regions
1 数据

ETOPO1由国际地球物理数据中心(national geophysical data center,NGDC)发布,网格分辨率为1′,包含全球范围内的陆地和海底地形数据[9],已有多项岩石圈Te研究[4-5, 10-11]将ETOPO1数据作为地形资料。为消除计算结果边界效应,参与计算的数据范围要远大于研究区域并涉及到海域。对于海底地形h,通过式h′=(ρc-ρw)h/ρc将其转换为等效高程h′,其中ρc=2 670 kg/m3为地壳平均密度,ρw= 1 030 kg/m3为海水密度。

根据地球重力场模型EIGEN-6C4计算1′×1′网格分辨率自由空气异常,截断阶数为2 190。采用柱体积分公式进行地形改正[12],在此基础上加上布格板改正得到布格重力异常(图 2)。将莫霍界面作为内部荷载作用参考面,莫霍界面深度、地壳密度和上地幔密度等岩石圈结构信息采用CRUST1.0模型[13]结果,地形、重力异常和地壳结构数据均采用基于WGS84参考椭球的兰伯特圆锥正形投影方法投影至平面,后重采样为10 km×10 km的网格数据。

图 2 安徽及周边地区布格重力异常 Fig. 2 Bouguer gravity anomaly of Anhui and its surrounding regions
2 计算方法

岩石圈Te体现岩石圈挠曲刚度,两者之间的关系如式(1)所示[1]。Kirby等[2]、Swain等[14]发展了一套基于小波变换的谱分析技术,能够同时解析空间和频率信息,克服多窗口谱方法缺陷,已在多项研究[2, 5, 10-11, 14]中得到应用。Kirby[15]开发的fan小波方法通过将二维Morlet小波进行叠加来获得fan小波,Nθ=int(Δθ/δθ)是构建fan小波的Morlet小波个数,其中δθ是Morlet小波方位角的采样间隔,本文取Δθ=180°,获得的各向同性fan小波用以反演各向同性Te。Morlet小波中间波数决定空间域和波数域分辨率,其值越大波数分辨率越高,空间分辨率越低。中间波数取值通常不小于2.668[10],采用多个中间波数(如2.668、3.081、3.773和5.336)进行计算发现,不同中间波数计算的Te分布结果非常类似,高值区和低值区分布一致,仅在空间分辨率上存在一定的差别,故本文仅采用中间波数为3.773的计算结果进行讨论。

(1)

式中,E为杨氏模量,ν为泊松比,D为岩石圈挠曲刚度,实测相关性表达式为:

(2)

式中,γ2为相关系数,k为波数,G为重力,H为地形,Re表示取实部, 星号表示共轭复数,尖括号表示在波数域内按角度求平均。确定重力异常和地形的相关性后,通过比较理论和实测谱参数求Te和荷载比,荷载比F与Forsyth[16]定义的f之间的关系为:

(3)

式中,F表示地下荷载与总荷载的比,完全为地表荷载时f=F=0,完全为地下荷载时f=∞、F=1。计算过程中用到的物理常量值见表 1

表 1 物理常量值 Tab. 1 Values of physical constants
3 计算结果及分析 3.1 岩石圈Te的空间分布

图 3为安徽及周边地区岩石圈Te计算结果。由图可见,Te值整体处于5~75 km之间,平均值为26.9 km,Te分布总体呈由南向北逐渐增大的趋势,表现出一定的分区性:大别造山带Te值处于10~30 km之间,靠近扬子断块和华北地块的东部地区Te值要略低于西部;下扬子地区Te值处于5~20 km之间,由南向北逐渐增大,至苏北盆地Te值上升至30~45 km左右;华北断块内Te值变化较大,河淮盆地Te值较高,最大超过70 km,向北和向南逐渐减小至20 km左右;苏鲁造山带Te值处于25~45 km之间,由靠近华北地块向靠近扬子地块一侧逐渐增大。

图 3 安徽及周边地区岩石圈Te Fig. 3 The Te of Anhui and its surrounding regions
3.2 岩石圈Te与区域地质构造

通常古老而冷的岩石圈(如克拉通和地盾)对应高Te值(> 60 km),抵抗形变能力较强;造山带地区通常对应中Te值(30~60 km);热而年轻的岩石圈构造者(如岩石圈拉张地区)Te值(< 30 km)通常较低,在多种构造过程中易产生变形[11]

河淮盆地和苏北盆地尽管经历了强烈的新生代弧后扩张[4],且受郯庐断裂带活动的影响,仍然表现为较高的Te值(50~70 km),除说明河淮盆地和苏北盆地岩石圈结构相对比较完整、强度较高外,也可能是Te计算结果受盆地特殊构造环境影响。大范围岩石圈Te研究表明,华南地区的江汉盆地和四川盆地等盆地区域均表现为局部高Te[4],可能是盆地区域长期沉积作用削弱了地形和布格重力异常之间的相关性,导致Te偏高。河淮盆地和苏北盆地同为新生代沉降区,从Airy均衡重力异常[12]来看,两者均表现出明显的负异常,表明相关区域可能存在导致质量亏损的较厚沉积层,造成Te偏高。郯庐断裂带对河淮盆地和苏北盆地两个Te局部高值区具有分割作用,且郯庐断裂带附近地区Te值普遍较低,最大不超过40 km,说明郯庐断裂带(深入岩石圈的深大断裂)破坏了板块几何结构完整性,降低了Te值。

大别-苏鲁造山带位于我国东部,扬子地块与华北地块间的陆-陆碰撞造山带,受郯庐断裂带影响,造山带被错开成大别、苏鲁两部分[6]。两者在岩石圈Te上存在较大差异,苏鲁造山带Te值处于25~45 km左右;大别造山带Te值明显低于苏鲁造山带,处于10~30 km之间。大别造山带较低Te值与特殊深部构造有关,多项宽角反射与折射地震探测[17]和天然地震波成像[6]等地球物理探测表明,大别造山带与华北地台及扬子地块的地壳结构明显不同。在大别造山带内部,莫霍界面也存在若干垂直错断。地震反射成果显示,扬子陆块岩石圈向北俯冲至少可以到达大别山下,俯冲深度超过莫霍面,进入岩石圈地幔并发生了断离。正是由于华北地块南向俯冲及扬子地块北向俯冲并在大别造山带下发生断离,俯冲和碰撞挤压造成大别造山带莫霍面垂向错断,加之郯庐断裂带影响,导致大别造山带岩石圈结构完整性受到破坏,从而表现出较低的Te值。

下扬子地区(不包括苏北盆地)的Te值较低,处于5~20 km之间。下扬子地区地表热流较高,下地壳和上地幔P波速度较低,同华南地块广阔的褶皱带和岩浆活动带一样,下扬子地区低Te值可能与长期大洋岩石圈俯冲消减作用有关,强度较弱的华南地块岩石圈在周边地块挤压之下产生变形,表现为广阔的褶皱带,且岩浆活动分布广泛[4]

3.3 岩石圈荷载

地下荷载通常由基性岩入侵、下地幔增生、热异常和成分变化等引起的深部物质密度横向变化造成,地表荷载通常由地形和大尺度地表密度变化引起[10]图 4为安徽及周边地区荷载比F分布,由图可见,F值处于0.4~0.9之间,东大别和下扬子地区靠近郯庐断裂带一侧F值较高,达0.8以上,河淮盆地和苏北盆地局部地区F值较低,处于0.4~0.6之间。对比图 3看出,FTe之间存在一定的负相关关系。在下扬子地区靠近郯庐断裂带的区域及大别造山带,F值较高,地下荷载起主导作用,说明岩石圈深部存在较大的物质密度横向变化,推测与扬子地块和华北地块向大别造山带下俯冲至地幔所导致的大别造山带地壳增厚和莫霍界面垂向错断有关,郯庐断裂带的切割作用对大别造山带和扬子地块接触部位的岩石圈完整性造成破坏,导致岩石圈内部质量分布不均匀,表现出较高的地下荷载。河淮盆地大部分地区、苏鲁造山带和苏北盆地局部地区F值均约为0.5,地表荷载和地下荷载所占比重基本相当。

图 4 安徽及周边地区荷载比 Fig. 4 Loading ratio of Anhui and its surrounding regions
3.4 岩石圈Te与地震活动

岩石圈Te空间分布和量级控制岩石圈在长期构造荷载下发生形变的位置、程度和方式,对地震活动有潜在的影响,岩石圈Te与地震分布存在一定联系[4, 10]

图 5为安徽及周边地区岩石圈Te与公元294年以来5.0级以上地震震中分布,共85次地震。从震中分布看出,大别造山带和华北地块地震活动均较多,扬子地块内地震活动多分布在苏北盆地南部边缘,苏鲁造山带和苏北盆地内部地震活动相对较少。从岩石圈Te空间分布和震中位置来看,地震活动主要集中在大别造山带特别是东大别、下扬子地区靠近郯庐断裂带一侧(Te低值区),以及河淮盆地、苏北盆地南部边缘(Te等值线较密集区),而在苏北盆地内部和苏鲁造山带南部(Te高值、等值线稀疏区)地震较少。为直观表达地震活动性与岩石圈Te之间的关系,将Te值与地震数量的对应关系以直方图形式表示,图 6(a)为公元294年以来5.0级以上地震与Te值的对应关系,图 6(b)为1970~2015年3.0级以上地震(来自中国地震台网)与Te值的对应关系。可以看出,5.0级以上地震主要集中分布在Te值10~40 km范围内,占地震总数的80%,3.0级以上地震主要集中在Te值10~50 km范围内,占地震总数的89%,发生地震次数最多的Te值范围均是10~20 km,震级越高,地震数量相对于Te值的集中性越好。

图 5 安徽及周边地区岩石圈Te与地震震中分布 Fig. 5 Distribution of Te and epicenters of Anhui and its surrounding regions

图 6 岩石圈Te与地震数量的对应关系 Fig. 6 Numbers of earthquakes correspond to different Te value ranges
4 结语

利用EIGEN-6C4重力异常和ETOPO1地形资料,采用基于fan小波的相关性分析方法计算得到安徽及周边地区岩石圈Te值和荷载比F,分析岩石圈Te和荷载比F分布特征,探讨其构造和地球动力学等含义。结果显示:

1) 安徽及周边地区岩石圈Te在各构造块体间存在明显差异。大别造山带和扬子地块内的下扬子地区、华北地块、河淮盆地、苏北盆地、苏鲁造山带等构造单元均表现出各自不同的Te分布特征。

2) 安徽及周边地区岩石圈Te和荷载比F主要体现区域地质构造特征和地球动力学背景。大别造山带较低的Te值和高F值主要反映大别造山带在俯冲和挤压作用之下,岩石圈结构完整性遭到破坏;下扬子地区的低Te值与华南地块长期以来的大洋岩石圈俯冲和消减作用有关;郯庐断裂带两侧较低的Te值与郯庐断裂带对两侧岩石圈结构造成破坏有关;河淮盆地和苏北盆地较高的Te值说明相关地区岩石圈强度较高,或沉积层较厚。

3) 研究区域内地震趋向于发生在岩石圈Te较小的大别造山带和下扬子地区靠近郯庐断裂带一侧,苏北盆地南部边缘和河淮盆地附近岩石圈Te水平梯度较大,地震活动较集中。绝大多数地震分布在Te值为10~50 km范围,10~20 km的Te值范围对应地震数量最多。较低的Te值说明岩石圈强度较弱,Te水平梯度较大的地方岩石圈强度发生显著变化,抵抗形变能力较弱而易发生地震。

致谢: 感谢澳大利亚科廷大学空间科学系Kirby J F博士提供岩石圈Te计算程序。

参考文献
[1]
Kirby J F, Swain C J. A Reassessment of Spectral Te Estimation in Continental Interiors: The Case of North America[J]. Journal of Geophysical Research, 2009, 114(B8): 4 918 (0)
[2]
Kirby J F, Swain C J. Mapping the Mechanical Anisotropy of the Lithosphere Using a 2D Wavelet Coherence, and Its Application to Australia[J]. Physics of the Earth and Planetary Interiors, 2005, 158(2-4): 122-138 (0)
[3]
郑勇, 李永东, 熊熊. 华北克拉通岩石圈有效弹性厚度及其各向异性[J]. 地球物理学报, 2012, 55(11): 3 576-3 590 (Zheng Yong, Li Yongdong, Xiong Xiong. Effective Lithospheric Thickness and Its Anisotropy in the North China Craton[J]. Chinese Journal of Geophysics, 2012, 55(11): 3 576-3 590) (0)
[4]
Mao X L, Wang Q, Liu S W, et al. Effective Elastic Thickness and Mechanical Anisotropy of South China and Surrounding Regions[J]. Tectonophysics, 2012, 550-553: 47-56 DOI:10.1016/j.tecto.2012.05.019 (0)
[5]
Chen B, Liu J, Chen C, et al. Elastic Thickness of the Himalayan-Tibetan Orogen Estimated from the Fan Wavelet Coherence Method, and Its Implications for Lithospheric Structure[J]. Earth and Planetary Science Letters, 2015, 409: 1-14 DOI:10.1016/j.epsl.2014.10.039 (0)
[6]
徐佩芬, 刘福田, 王清晨, 等. 大别-苏鲁碰撞造山带的地震层析成像研究——岩石圈三维速度结构[J]. 地球物理学报, 2000, 43(3): 377-386 (Xu Peifen, Liu Futian, Wang Qingchen, et al. Seismic Tomography beneath the Dabie-Sulu Collision Orogeny: 3-D Velocity Structures of Lithosphere[J]. Chinese Journal of Geophysics, 2000, 43(3): 377-386 DOI:10.3321/j.issn:0001-5733.2000.03.011) (0)
[7]
郭兴伟, 张训华, 温珍河, 等. 中国海陆及邻域大地构造格架图编制[J]. 地球物理学报, 2014, 57(12): 4 005-4 015 (Guo Xingwei, Zhang Xunhua, Wen Zhenhe, et al. Complication of the Tectonic Framework Map of China Seas and Land and Adjacent Regions[J]. Chinese Journal of Geophysics, 2014, 57(12): 4 005-4 015) (0)
[8]
Förste C, Bruinsma S L, Abrikosov O, et al. EIGEN-6C4: The Latest Combined Global Gravity Field Model Including GOCE Data up to Degree and Order 2 190 of GFZ Potsdam and GRGS Toulouse[EB/OL]. http://icgem.gfz-potsdam.de/ICGEM/documents/Foerste-et-al-EIGEN-6C4.pdf, 2016 (0)
[9]
Amante C, Eakins B W. ETOPO1 Arc-Minute Global Relief Model: Proce-Dures, Data Sources and Analysis[EB/OL]. http: //www.ngdc.noaa.gov/mgg/global/global.html, 2009 (0)
[10]
Alberto J D, Javier R, Marta P G, et al. Spatial Variations of Effective Elastic Thickness of the Lithosphere in Central America and Surrounding Regions[J]. Earth and Planetary Science Letters, 2014, 391: 55-66 DOI:10.1016/j.epsl.2014.01.042 (0)
[11]
Chen B, Liu J X, Kaban M K, et al. Elastic Thickness, Mechanical Anisotropy and Deformation of the Southeastern Tibetan Plateau[J]. Tectonophysics, 2014, 637: 45-56 DOI:10.1016/j.tecto.2014.09.007 (0)
[12]
黎哲君, 徐如刚, 张玮晶, 等. 华北断块区均衡重力异常、构造特征及地震活动[J]. 地震地质, 2015, 37(2): 496-509 (Li Zhejun, Xu Rugang, Zhang Weijing, et al. Isostatic Gravity Anomaly, Structural Features and Seismicity of North China Craton[J]. Seismology and Geology, 2015, 37(2): 496-509 DOI:10.3969/j.issn.0253-4967.2015.02.012) (0)
[13]
Laske G, Masters G, Ma Z, et al. Update on CRUST1.0 a 1-Degree Global Model of Earth's Crust[J]. Geophys Res Abstr, 2013, 15: 2658 (0)
[14]
Swain C J, Kirby J F. An Effective Elastic Thickness Map of Australia from Wavelet Transforms of Gravity and Topography Using Forsyth's Method[J]. Geophys Res Lett, 2006, 33(2) (0)
[15]
Kirby J F. Which Wavelet Best Reproduces the Fourier Power Spectrum?[J]. Computers & Geosciences, 2005, 31(7): 846-864 (0)
[16]
Forsyth D W. Subsurface Loading and Estimates of the Flexural Rigidity of Continental Lithosphere[J]. J Geophys Res, 1985, 90(B14): 12 623-12 632 DOI:10.1029/JB090iB14p12623 (0)
[17]
刘福田, 徐佩芬, 刘劲松, 等. 大陆深俯冲带的地壳速度结构-东大别造山带深地震宽角反射/折射研究[J]. 地球物理学报, 2003, 46(3): 366-372 (Liu Futian, Xu Peifen, Liu Jinsong, et al. The Crustal Velocity Structure of the Continental Deep Subduction Belt: Study on the Eastern Dabie Orogeny by Seismic Wide-Angle Reflection/Refraction[J]. Chinese Journal of Geophysics, 2003, 46(3): 366-372 DOI:10.3321/j.issn:0001-5733.2003.03.014) (0)
Effective Elastic Thickness of Anhui and Its Surrounding Regions
LI Zhejun1     ZHANG Yi1     YANG Guangliang2     JIA Yupeng3     XIAO Weipeng1     TAN Xin1     LONG Jianfeng1     
1. Earthquake Administration of Anhui Province, 558 West-Changjiang Road, Hefei 230031, China;
2. Institute of Seismology, CEA, 40 Hongshance Road, Wuhan 430071, China;
3. Geophysical Exploration Center, CEA, 75 Wenhua Road, Zhengzhou 450002, China
Abstract: The Te map of Anhui and surrounding regions is calculated from EIGEN-6C4 gravity anomalies and ETOPO1 data using the isotropic fan wavelet method. The tectonic structure and geodynamic features of this area are discussed according to the Te map. The results indicate that the Te of Anhui and its surrounding regions varies from 5 to 75 km. The Te values of the Dabie orogenic belt and lower Yangtze region are low, as are the Te values along Tan-Lu fault. High Te values appear at Hehuai basin and Subei basin. The Te distribution mainly reflects the particular characteristics of tectonic structure and geodynamic features of corresponding areas.Most of the seismic activities located in regions with Te values vary from 10 to 50 km; the spatial distribution of earthquakes indicates that earthquakes tend to occur in regions with low Te or steep Te gradients, further indicating that lithospheres of regions with low Te or steep Te gradients show higher seismicity rateand are not strong enough to resist deformation.
Key words: lithosphere; effective elastic thickness; Anhui; seismicity