青藏高原下地壳流动方式的数值模拟研究

摘要/Abstract

摘要： 青藏高原存在柔性下地壳流动被越来越多的学者接受, 但是关于下地壳流动方式及速度存在争议。地表运动有GPS等直接测量, 上地幔运动有S波分裂间接反映, 下地壳运动目前没有直接观测手段, 使得开展数值分析非常重要。本文利用三维球壳黏弹性有限元模型研究了青藏高原下地壳柔性流动方式和流动速度。本文通过对地表GPS观测资料的拟合与不同数值模型的对比分析, 认为青藏高原柔性下地壳东向流动遇到四川盆地的抵阻, 下地壳物质可能仅在高原东南方向存在物质外溢通道, 而在高原东北方向不存在类似的物质通道; 下地壳的流动速度比地表运动速率每年快几毫米至十几毫米, 对应的黏滞系数为10¹⁸~10¹⁹ Pa·s。

关键词: 青藏高原, 下地壳流动, 三维黏弹性有限元模型, 数值模拟

Abstract: The idea of ductile flow beneath the Tibetan Plateau is widely accepted. However, there is a dispute on the way and speed of the lower crustal flow. Surface motion has direct measurements such as GPS, and mantle movement has indirect methods such as S-wave splitting. As to the lower crustal movement, there is no direct observation, so that numerical analysis is very important. In this paper, we constructed a 3-D spherical viscoelastic finite element model of the Tibetan Plateau to simulate the crustal flow pattern. By fitting GPS data and comparing different models, we considered that the eastward flow of the ductile lower crust in the Tibetan Plateau has been resisted by the Sichuan Basin, a channel exists in the south-eastern corner of the plateau, and similar channels of substances do not exist on the northeast corner of the plateau. The lower crust moves several milimeter faster than the upper crust, the corresponding viscosity of the lower crust was between 10¹⁸~10¹⁹ Pa·s.

Key words: Lower crustal flow, Finite element modeling, Viscos-elastic, Tibetan Plateau

中图分类号:

P315.7

曹建玲, 王辉, 张晶. 青藏高原下地壳流动方式的数值模拟研究[J]. 地震, 2013, 33(4): 55-63.

CAO Jian-ling, WANG Hui, ZHANG Jing. Numerical Modeling of Lower Crustal Flow beneath the Tibetan Plateau[J]. EARTHQUAKE, 2013, 33(4): 55-63.

参考文献

[1] Yang Y, Liu M. Crustal thickening and lateral extrusion during the Indo-Asian collision: a 3D viscous flow model[J]. Tectonophysics, 2009, 465: 128-135.
[2] Wang C Y, Han W B, Wu J P, et al. Crustal structure beneath the eastern margin of the Tibetan Plateau and its tectonic implications[J]. Journal of Geophysical Research, 2007, 112: B07307, doi:10.1029/2005JB003873.
[3] 雷建设, 赵大鹏, 苏金蓉, 等. 龙门山断裂带精细结构与汶川地震发震机理[J]. 地球物理学报, 2009, 52(5): 339-345.
[4] Xu Y, Li Z W, Huang R Q, et al. Seismic structure of the Longmen Shan region from S-wave tomography and its relationship with the Wenchuan M_S8.0 earthquake on 12 May 2008, southwestern China[J]. Geophysical Research Letters, 2010, 37: L02304, doi:10.1029/2009GL 041835.
[5] Royden L H, Burchfiel B C, King R W, et al. Surface Deformation and Lower Crustal Flow in Eastern Tibet[J]. Science, 1997, 276(2): 788-790.
[6] Beaumont C, Jamieson R A, Nguyen M H, et al. Crustal channel flows: 1. Numerical models with applications to the tectonics of the Himalayan-Titetan orogen[J]. Journal of Geophysical Research, 2004, 109: B06406, doi:10.1029/2003JB002809.
[7] Copley A, McKenzie D. Models of crustal flow in the India-Asia collision zone[J]. Geophysical Journal International, 2007, 169(2): 683-698.
[8] 熊熊, 许厚泽, 滕吉文. 青藏高原物质东流的岩石层力学背景探讨[J]. 地壳形变与地震, 2001, 21(2): 2-6.
[9] 汪一鹏, 沈军, 王琪, 等. 川滇块体的侧向挤出问题[J]. 地学前缘, 2003, 10(特刊): 188-192.
[10] An M J, Shi Y L. Lithospheric thickness of the Chinese continent[J]. Physics of The Earth and Planetary Interiors, 2006, 159(3-4): 257-266.
[11] 安美建, 石耀霖. 中国大陆地壳和上地幔三维温度场[J]. 中国科学(D辑), 2007, 37(6): 736-745.
[12] Bai D, Unsworth M J, Meju M A, et al. Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging[J]. Nature Geosci, 2010, 3(5): 358-362.
[13] Shen F, Royden L H, Burchfiel B C. Large-scale crustal deformation of the Tibetan Plateau[J]. Journal of Geophysical Research, 2001, 106(4): 6793-6816.
[14] Royden L. Coupling and decoupling of crust and mantle in convergent orogens: Implications for strain partitioning in the crust[J]. Journal of Geophysical Research, 1996, 101(8): 17679 -17705.
[15] 曹建玲, 石耀霖, 张怀, 等. 青藏高原GPS位移绕喜马拉雅东构造结顺时针旋转成因的数值模拟[J]. 科学通报, 2009, 54(2): 224-234.
[16] 朱守彪, 石耀霖. 用遗传有限单元法反演川滇下地壳流动对上地壳的拖曳作用[J]. 地球物理学报, 2004, 47(2): 232-239.
[17] Royden L, Burchfiel B C, van der Hilst R D. The geological evolution of the Tibetan Plateau[J]. Science, 2008, 321: 1054-1058.
[18] Clark M K, Bush J W M, Royden L H. Dynamic topography produced by lower crustal flow against rheological heterogeneities bordering the Tibetan Plateau[J]. Geophysical Journal International, 2005, 162: 575-590.
[19] 张中杰, 滕吉文, 李英康, 等. 藏南地壳速度结构与地壳物质东西向“逃逸”以佩枯错-普莫雍错宽角反射剖面为例[J]. 中国科学(D辑), 2002, 32(10): 793-798.
[20] 石耀霖, 曹建玲. 中国大陆岩石圈等效粘滞系数的计算与讨论[J]. 地学前缘, 2008, 15(3): 82-95.
[21] Gan W J, Zhang P Z, Shen Z K, et al. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements[J]. Journal of Geophysical Research 2007, 112: B08416, doi:10.1029/2005JB004120.
[22] Liu M, Yang Y. Extensional collapse of the Tibetan Plateau : Results of three-dimensional finite element modeling[J]. Journal of Geophysical Research, 2003, 108(8): 2361-2380.
[23] 张晁军, 曹建玲, 石耀霖. 从震后变形探讨青藏高原下地壳黏滞系数[J]. 中国科学(D), 2008, 38(10): 1250-1257.
[24] Rippe D, Unsworth M J. Quantifying crustal flow in Tibet with magnetotelluric data[J]. Physics of the Earth and Planetary Interiors, 2010, 179: 107-121.
[25] Yao H, Beghein C, van der Hilst R. Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis: II-Crustal and upper-mantle structure[J]. Geophys J Int, 2008, 173: 205-219.
[26] Wang H, Cao J L, Zhang H, et al. Numerical simulation of the influence of lower-crustal flow on the deformation of the Sichuan-Yunnan region[J]. Acta Seismological Sinica, 2007, 20: 617-627.
[27] 安美建, 石耀霖. 中国大陆岩石圈厚度分布研究[J]. 地学前缘, 2006, 13(3): 23-30.
[28] 韦伟, 孙若昧, 石耀霖. 青藏高原东南缘地震层析成像及汶川地震成因[J]. 中国科学(D辑), 2010, 40(7): 831-839.
[29] 刘明军, 李松林, 方盛明, 等. 利用地震波速研究青藏高原东北缘地壳组成及其动力学[J]. 地球物理学报, 2008, 51(2): 412-430.
[30] 崔笃信, 王庆良, 胡亚轩, 等. 青藏高原东北缘岩石圈变形及其机理[J]. 地球物理学报, 2009, 52(6): 1490-1499.