从震后形变探讨青藏高原下地壳黏滞系数.pdfVIP

中国科学 D 辑 :地球科学 2008 年 第 38 卷 第 10 期 : 1250 ~ 1257 SCIENCE IN CHINA PRESS 从震后形变探讨青藏高原下地壳黏滞系数 ①②* ② ② 张晁军 , 曹建玲 , 石耀霖 ① 中国地震局地震台网中心, 北京 100045; ② 中国科学院研究生院, 计算地球动力学实验室, 北京 100049 * E-mail: zcj72@ 收稿日期: 2008-01-02; 接受日期: 2008-07-25 国家重点基础研究发展计划(编号: 2004CB418405)、国家自然科学基金(批准号:、中国科学院知识创新工程重要方向项目(编号: KZCX2-YW-123)和地震基本科研业务项目(批准号:05)资助 摘要 青藏高原下地壳黏滞系数究竟是多少, 已成为深入定量研究中的突出问题, 它的数值量 关键词 级将极大地影响定量模拟的结果. 为了获得青藏高原下地壳的黏滞系数, 从3 条途径对该参数进 流变模型 行了计算: 一是基于对青藏高原深部温度状态改进了的估计, 用流变定律和 GPS 求得的应变速 应变速率 率对高原北部下地壳黏滞系数重新进行了估算, 获得的昆仑山地区中地壳等效黏滞系数为 流变定律 20 22 19 21 震后形变 10 ~10 Pa ·s 量级, 下地壳等效黏滞系数在10 ~10 Pa ·s 之间. 二是用3 种流变模型对2001 年 黏滞系数 昆仑山MS8.1 级地震震后跨断层GPS 站点记录到的震后变形进行了模拟, 得到的下地壳黏滞系 数为 1017 Pa ·s 量级. 三是用黏弹性模型对炉霍地震后的跨断层形变曲线进行了拟合, 得到的下 地壳黏滞系数为1019 Pa ·s 量级. 前人研究等效黏滞系数时忽视了等效黏滞系数与应变速率存在 的非线性关系, 本研究结果协调地解释了实验室实验、大地震后较短时期的变形和大地震后较长 时间变形下, 其等效黏滞系数存在差异的问题. 青藏高原下地壳黏滞系数究竟是多少, 已成为 黏滞系数为 1018 Pa ·s量级, 在四川盆地边缘得到的 深入定量研究中突出的问题, 它的数值量级将极大 21 Pa ·s量级, 拟合盆地平坦地区得到 黏滞系数为 10 地影响定量模拟的结果. 许多学者用不同的方法对 16 [4] 的黏滞系数为 10 Pa ·s量级. Hilley 等 用地震复发 这一问题进行了探讨. 例如, Royden 等[1]利用解析方 周期方法得到的青藏高原下地壳黏滞系数在 法研究青藏高原抬升高度及地壳厚度随时间变化问 18 21 [5] 10 ~10 Pa ·s 之间. Ryder 等 用 Maxwell 流体模 题时, 认为青藏高原地区地壳非常柔软, 其黏滞系数 型、标准线弹性固体模型, 借助于 1997 年 11 月 8 日 范围很宽, 介于 1012~1021 Pa ·s之间, 当应变速率为 Mw7.6 级玛尼地震后 4 a非连续的 InSar 观测数据资 3×10− 15/s 时, 得到的黏滞系数为 1018 Pa