中国大陆地处印度-欧亚板块碰撞引起的复杂变形区内，地震活动频发.统计表明，全球约有33%的大陆地震发生在中国境内(Zhang et al.，2003).随着社会经济的不断发展和城市化进程的不断加剧，出现了越来越多的人口高密区域.因此，分析中国大陆的地震危险性具有非常重要的意义.

随着GPS(Global Positioning System)等空间大地测量技术的发展，研究中国大陆的变形特征已经有很多研究成果(Jiang and Ma.，2003; Calais et al.，2006; 郑文俊等，2009).但是，前人研究局限于对形变场模型反映的研究区域变形特征的分析，在前人的研究结果中，变形较大的区域存在应变积累，也即为潜在的震源区.然而，GPS数据得到的形变场模型包含瞬时变形、断层蠕滑等非震形变、由已发地震释放的永久形变和震间弹性变形.因此，用大地测量资料分析研究区域的地震危险性，需要扣除已发地震产生的永久形变及其它非震形变.

Holt等(1995)利用历史地震资料反演了中、东亚的速度场；任金卫等(1999，2002)和Holt等(2000)利用地震资料、GPS等大地测量资料和第四纪活断层资料联合反演了中国大陆及其邻区的应变率场和速度场.上述模型分别从运动学和动力学的角度实际客观地描述了中国大陆及其邻区现今变形特征，为进一步分析该区域的地震危险性提供了基础资料.

本文在任金卫等(2002)的研究基础上，详细搜集整理了1900-2013年间发生在中国大陆及其邻区的4级及其以上历史地震资料，并对其中47个地表破裂带较长(40 km)的大震做了分段.其次，利用历史地震资料、第四纪活断层资料和GPS等大地测量资料联合反演了研究区域的应变率场和速度场.最后，利用历史地震资料和反演的应变率场估算了研究区域的地震矩亏损，为分析研究区域的地震危险性提供参考资料.

首先搜集了1900-2013年间发生在中国大陆及其邻区的1349个4级及其以上历史地震资料(图 1)，震源球大小用矩震级(M_W)表示.其中，1976-2013年间的地震资料来自哈佛大学全球地震中心矩张量目录(CMT)；1964-1976年的历史地震资料通过总结前人的研究成果得到(Chen et al.，1977; Molnar et al.，1983，1984; Holt et al.，1991a)；1900-1964的地震资料主要来自于Molnar等(1984)和Anne等(1984)的研究成果.图 1中，1900-2013年间，中国东北、蒙古东部和华南地区地震较少.

图 1

Fig. 1

图 1 中国大陆及其邻区地震震源机制 Fig. 1 Focal mechanisms of Chinese mainland and its vicinity

任金卫等(2002)将中国大陆及其邻区划分成不规则的地质构造单元，并结合地震资料、GPS等大地测量资料和第四纪活断层资料，反演了研究区域现今应变率场和速度场.为了减小大震的影响范围和提高反演结果的空间分辨率，本文对网格进行加密，从原来的12×11个增加到96×88个不规则的四边形网格.利用改进的双三次样条方法(Haines et al.，1993; Beaven et al.，2001)反演了中国大陆及其邻区的应变率场(^mod，图 2)，并利用^mod估算了模型矩率(₀^mod)(Haines et al.，1993; Beaven et al.，2001)为

图 2

Fig. 2

图 2 利用地震资料、第四纪活断层资料和GPS资料联合反演的研究区域的应变率场(单位:10-9/a)，图例表示50×10-9/a的拉张应变率 Fig. 2 Strain rate model determined from earthquake data,Quaternary fault rates and GPS observations (units:10-9/a),Legend is 50×10-9/a of tensional strain rate

利用Kostrov求和公式(Molnar et al.，1984; Kostrov.，1974)(式2)和历史地震资料(图 1)估算各网格的应变率观测值(^obs)为

利用式1和^obs估算各网格的观测矩率(₀^obs)，则剩余矩率Δ₀=₀^mod-₀^obs.各网格的剩余应变率张量Δ=^mod-^obs(图 3a)，再利用Δ求地震震源参数(式4~11)(Shentu et al.，1998)，则可将Δ₀用震源球表示(图 3b).图 3(a，b)中，红色表示地震矩率亏损，即Δ₀>0；蓝色为地震矩率盈余，即Δ₀ < 0.断层面走向(θ)，倾角(δ)，滑动角(λ)表示为

其中，

图 3

Fig. 3

图 3 (a)剩余应变率,模型应变率(mod)和应变率观测值(obs)之差. 红色为Δ0>0,蓝色为Δ0<0,图例表示50×10-9/a的拉张应变率. (b)剩余地震矩率.其中,红色表示地震矩率亏损(Δ0>0); 蓝色表示地震矩率盈余(Δ0 < 0) Fig. 3 (a) Residual strain rate from the difference between predicted strain rate (mod) and strain rate observations (obs),in which red represents Δ0>0 and blue represents Δ0 < 0.Legend is 50×10-9/a of tensional strain rate. (b) Residual moment rate,in which red represents moment rate deficit (Δ0>0) and blue represents moment rate surplus (Δ0 < 0)

图 3中的应变率和地震矩率盈余量(蓝色)分布很不合理，主要表现在喜马拉雅逆冲断裂带，东昆仑断裂、海原断裂带和蒙古高原等有大震发生的地区.在喜马拉雅逆冲断裂带，1934年比哈尔-尼泊尔地震的发震断裂为近E-W的逆冲断裂，其地表破裂长度达到200 km以上(Singh and Harsh，1980)，然而，图 3所示的盈余量主要分布在震中和其北边的网格，并未按照地震产生的形变带分布.昆仑山口西地震的发震断裂长度达到420 km(Lin et al.，2002)，破裂长度至少横穿三个网格，破裂带东段所在的网格却存在较大的地震矩率和应变率亏损.上述不合理性存在的主要原因是估算剩余应变率和地震矩率时，将大震以点源的形式释放，使大震形变分布过于集中而并未按照其破裂带分布，从而使存在大震的网格应变率观测值出现峰值.

为了解决2.2中存在的问题，本文将历史地震分为两类：第一类包括所有小于7级和7级以上但地表破裂带较短(40 km，接近于本文网格尺度)的地震；第二类为47个7级及其以上且地表破裂带较长的地震.详细搜集整理了第二类地震的震源参数、地表破裂带和地表位移资料(表 1).根据资料的完整程度将其分成A，B，C三类：A类，震源参数完整且有地表位移观测资料；B类，震源参数相对完整，存在发震断裂和地表破裂带资料，但缺少地表位移观测数据；C类，缺少发震断裂和地表破裂带资料.

表 1(Table 1)

表 1 1900-2013年间中国大陆及其邻区7级及其以上地震资料总结 Table 1 Summary of materials about 47 earthquakes with magnitude greater than or equal to 7.0 in Chinese Mainland and its vicinity between 1900-2013 编号 发震位置 发震时间 破裂长度(深度) 走向(°) 发震断裂 震级 地震矩(×10 20 N.m) 位错(m) 类型 分段个数 主要参考资料 1 蒙古 1905/7/9 130 60 车车尔勒格断裂 8.4 2.2±0.6 2.3 A 2 I.Baljinnyam, et al ., 1993 2 蒙古 1905/7/23 375 80 22 90 140 0 杭爱断裂 8.7 21±0.49 1.6±0.5 0.2 8 3.1 2 A 8 I.Baljinnyam, et al ., 1993 3 俄罗斯 1957/6/27 24 26 100 100 穆亚/乌多坎断裂 7.9 1.4 3.3 2.2 A 2 Kurushin, R. A. 2008 4 蒙古 1957/12/4 260 105 95 115 博格多断裂 8.3 18.0 6.2 6.5 4.1 A 10 I.Baljinnyam, et al ., 1993Jin-Hyuck choi,et al .,2012 5 蒙古 1967/1/5 12(10) 15(10) 20(6) 2 20 310 图利特乌尔断裂 7.7 0.36 1.6 3.2 2.2 A 3 A.Bayasgalan, et al .,1999 6 富蕴新疆 1931/10/10 180 160 富蕴断裂带 7.9 9.7±2.4 6.3 A 5 Lin,et al ., 1998新疆地震局, 1985 7 阿图什新疆 1902/8/22 250 260 托特拱拜孜断裂 8.6 30.0 4.0 B 6 冯先岳等,1985郭增建等,1988 8 玛纳斯新疆 1906/12/22 120 270 霍尔果斯-玛纳斯-吐古鲁逆断裂 8.3 1.3±0.9 2.0 B 3 Molnar, P, et al .,2000王春镛等,2001 9 吉尔吉斯斯坦 1911/1/11 265(70) 265 克敏断裂带 8.7 7.5±2.9 2.3 B 7 D. Delvaux, et al ., 2001 10 吉尔吉斯斯坦 1946/11/2 80 120 塔拉索-费尔干纳断裂 7.6 1.1±0.3 2.0 B 3 Simpson, D. W. et al ., 1981 11 且末新疆 1924/7/3 46 80 皮什盖-塔勒看里克 7.2 0.15 1 B 5 冯先岳等,1985 12 民丰新疆 1924/7/11 46 80 皮什盖-塔勒看里克 7.2 0.15 1 B 5 冯先岳等,1985 13 乌恰新疆 1974/8/11 180(26) 32.8 - 7.3 3.7 4.3 C 8 杨懋源等1985 14 于田新疆 2008/3/21 12 12 21 194 167 200 阿尔金断裂 7.5 a 0.54 0.5 3.2 3.0 A 3 徐锡伟等,2011Elliot, J. R, et al .,2010 15 印度 1905/4/4 100 300 喜马拉雅逆冲断裂带 8.6 59.6 3.6 B 6 Nicholas,et al .,2000Wallace, et al .,2005 16 比哈尔-尼泊尔 1934/1/15 129 280 喜马拉雅逆冲断裂带 8.4 44.4 2.8 B 6 Molnar, P, et al .,1989D.D.Singh, et al .,1980 17 巴基斯坦 1935/5/30 209 33 奎塔断裂 7.5 1.75 3.3 B 4 D.D.Singh, et al .,1980 18 朗县西藏 1947/7/29 100 270 — 7.7 2.0 3.0 C 3 Molnar,P. et al ., 1984 19 察隅西藏 1950/8/15 100 250 170 260 — 8.7 20 60 8 C C 9 8 Rolando, et al ., 1989 Menahem, 1974 20 当雄西藏 1951/11/18 91(10) 145 147 135 崩错断裂 8.0 2.7 8 A 3 吴章明等,1989 21 西藏中部 1973/7/14 64 190 — 7 0.1 0.6 C 2 Molnar, P. et al .,1983 22 西藏南部 1975/1/19 70 0 — 7 0.16 0.8 C 2 Molnar, P. et al .,1983 23 玛尼西藏 1997/11/8 120 79 西昆仑断裂西段 7.9 a 2.23 6.0 A 4 单新建等,2002Gilles,et al ., 1999 24 昆仑山口西西藏 2001/11/14 420 94 东昆仑断裂 8.2 a 5.90 4.5 A 9 陈杰等,2003Xu, Xiwei,et al .,2006 25 海原宁夏 1920/12/20 225 115 海原断裂带 8.7 12 8.0 A 20 宁夏地震局, 1990时振梁,苏公望等, 1987 26 古浪甘肃 1927/5/23 110 75 33 290 280 342 皇城-双塔 冷龙岭断裂 中坝-磨咀子隐伏断裂 8.3 4.3 4.7 A 10 郭增建等,1988侯康明等,1998郑文俊等,2004 27 玉门甘肃 1932/12/25 116(15) 115 昌马断裂带 7.6 3.55 3.8 A 4 兰州地震研究所, 1992侯康明等,1990 28 阿兰湖东青海 1937/1/7 240 110 东昆仑断裂带 7.6 4.8 7 A 4 青海地震局, 1999李陈侠,2006 29 阿兰湖青海 1963/4/19 72 277 东昆仑断裂带 7.1 0.14 0.74 C 2 沈正康等,2003 30 达日青海 1947/3/17 150 135 日查-克授滩断裂 7.75 3.0 4 C 6 叶世强等,1993戴华光等,1983 31 道孚四川 1904/8/30 50 322 鲜水河断裂 7.0 0.13 0.78 B 2 Papadimitriou,et al .,2004Allen, et al ., 1991 32 炉霍-道孚四川 1923/3/24 60 116 鲜水河断裂带 7.25 0.62 3 B 2 Papadimitriou,et al .,2004Allen, et al ., 1991 33 茂汶四川 1933/8/25 100 14 岷江断裂 7.4 1.90 0.5 B 3 唐荣昌等,1983 34 炉霍四川 1973/2/6 105 116 鲜水河断裂 7.7 1.90 4 A 3 Allen, et al ., 1991Huilan,Zhou,et al .,1983 35 汶川四川 2008/5/12 240 75 6 222 225 320 龙门山断裂带 8.2 a 9.0 4 3.5 4.8 A 7 Ran Yongkang, et al ., 2009Xiwei, Xu,et al ., 2009 36 昭通云南 1974/5/11 10 23 45 200 巧家-莲峰断裂 7.1 0.20 2.8 A 2 朱成男等,1976刘万琴等,1979 37 通海云南 1970/1/4 56(9) 116 曲江断裂 7.7 0.87 2.1 A 3 王春镛等,1978中国震例,1988 38 澜沧云南 1988/11/6 50(26) 313 木戛断裂和大塘子断裂 7.4 a 0.46 2.0 B 4 俞维贤等,1991吴小平等,2007 39 缅甸 1908/12/12 120 25 实皆断裂 7.5 b 1.2 2.3 C 3 Le Dain Anne Ybonne, et al ., 1984 40 缅甸 1929/8/8 80 -10 实皆断裂 7.0 b 0.4 1.3 C 2 Le Dain Anne Ybonne, et al ., 1984 41 缅甸 1930/5/5 131(15) -10 实皆断裂 7.4 b 0.9 3.0 A 3 Le Dain Anne Ybonne, et al ., 1984 42 缅甸 1930/12/3 120(15) -10 实皆断裂 7.5 b 0.9 3.0 A 3 Hiroyuki Tsutsumi,et al ., 2009 43 缅甸 1931/1/27 180(15) 25 实皆断裂 7.7 b 1.5 2.6 B 5 Le Dain Anne Ybonne, et al ., 1984 44 缅甸 1946/9/12 155(15) -5 实皆断裂 7.5 b 1.2 2.3 B 6 Hurukawa Nobuo, et al ., 2011 45 缅甸 1946/9/12 180(15) -5 实皆断裂 7.8 b 1.2 2.3 B 5 Hurukawa Nobuo, et al ., 2011 46 海城辽宁 1975/2/4 52(12) 288 海城-岫岩断裂 7.3 0.3 2.8 A 3 John Cipar, 1979林邦慧等,1979 47 唐山河北 1976/7/28 24(16) 24(16) 35 5 唐山断裂带 7.8 a 2.77 3.54 A 2 Butler, 1997Huang,et al .,1997 注：a为面波震级，其震源参数来自于CMT；b震源参数来自于参考文献Anne(1984)；其它地震震源参数来自参考文献Molnar(1984).限于文章的篇幅，表1中的参考文献并未在文中列出，将有后文予以详述(李长军等,2014).

表 1 1900-2013年间中国大陆及其邻区7级及其以上地震资料总结 Table 1 Summary of materials about 47 earthquakes with magnitude greater than or equal to 7.0 in Chinese Mainland and its vicinity between 1900-2013

对A类地震，首先根据地表破裂带展布和本文网格之间的位置关系将破裂带分段，再根据地表位移分布估算各子破裂段的地震矩和矩震级.对缺少地表位移观测资料的B类地震，由于地震地表位移总体呈现中间大、两边小的分布特征，在此假定其服从三角形分布(Martin and Beroza，2000)，利用与A类地震相同的分段方法将其分段.C类地震缺少地表位移和破裂带资料，首先，利用经验公式(Wesnousky，2008)估算破裂带长度，假定震源位于破裂带中间，根据走向参数和震源位置确定地表破裂带位置，再利用B类地震的方法分段.地震分段时，由于各破裂子段长度已知，通过调整破裂矩形面宽度，使分段后子地震的地震矩之和与该地震的地震矩观测值保持一致，47个大震的分段结果见图 4.

图 4

Fig. 4

图 4 47个大震子地震的震源机制 Fig. 4 Focal mechanisms of sub earthquakes of 47 major events

将所有第一类地震和分段处理的第二类地震利用Kostrov求和公式估算各网格的应变率观测值(^obs)，再根据(2.2)的方法估算各网格的剩余应变率(Δ)和剩余地震矩率(Δ).为了使表述更为清晰，将Δ转化成与地震矩率相关的应变率Δ(式12)，如图 5a所示.其中，红色表示地震矩率亏损(即Δ₀>0)，蓝色表示地震矩率盈余(即Δ₀ < 0)，公式(11)为

图 5

Fig. 5

图 5 (a)中国大陆及其邻区的剩余应变率(扣除蠕动等非震形变);(b)由图5a的剩余应变率估算的地震矩亏损(红色)和盈余(蓝色) Fig. 5 (a)Residual strain rate of Chinese Mainland and its vicinity (Part of aseismic deformation is removed);(b)Moment deficit (red) and surplus (blue) evaluated from residual strain rate in figure 5a

图 5a中，应变率亏损较大的区域主要包括喜马拉雅断裂带、天山、阿尔金断裂带、青藏高原中南部、川滇地区、鄂尔多斯块体周缘和东南沿海等区域.在鄂尔多斯块体内部，塔里木盆地，华南块体和中国东北的亏损量相对较小，与该区域较低的构造活动特征(张培震等，2013)一致.图 5b为利用剩余应变率估算的研究区域113年间的地震矩亏损(红色)和盈余(蓝色)，其中的地震震源机制反映了相关区域的力学性质，震源球由矩震级大小确定.表 2为研究区域主要活动块体和断裂的地震矩亏损量及对应的矩震级.

表 2(Table 2)

表 2 中国大陆及其邻区主要活动地区的地震矩亏损 Table 2 Moment deficit of main active regions in Chinese Mainland and its vicinity 区域 WH CH EHS CST WK EKL NCB AT TS CD NLM SLM OB ESC BK M 0 10.4 19.4 3.4 8.9 3.7 0.3 2.5 5.1 5.1 9.7 0.4 0.3 2.4 0.9 2.7 M W 7.9 8.1 7.6 7.9 7.6 6.9 7.5 7.7 7.7 7.9 7.0 6.9 7.5 7.2 7.6 注： M 0单位为10 27 dyne.cm， M W为矩震级.AT, 阿尔金断裂带；BK, 贝加尔裂谷带；CD，川滇地区；CH，喜马拉雅断裂带中部；CST，青藏高原中南部；EHS，喜马拉雅东构造结；EKL，昆仑断裂带东段；ESC，东南沿海；NCB，柴达木盆地北缘；NLM，龙门山断裂北段；OB，鄂尔多斯块体周缘；SLM，龙门山断裂南段；TS，天山；WH，喜马拉雅断裂带西段；WK，青藏高原西北部及喀喇昆仑断裂带.

表 2 中国大陆及其邻区主要活动地区的地震矩亏损 Table 2 Moment deficit of main active regions in Chinese Mainland and its vicinity

在青藏高原地区，地震矩亏损较大的区域主要位于喜马拉雅断裂带、阿尔金断裂带、柴达木盆地北缘和青藏高原中南部.在低海拔的喜马拉雅逆冲推覆带和祁连山以逆冲型地震为主，而在高原内部以正断和走滑型地震为主，与该区域的构造活动特征一致(张培震等，2013).喜马拉雅逆冲断裂带由于受康格拉地震(1905)、比哈尔-尼泊尔地震(1934)和察隅地震(1950)的影响，该断裂带的东段和中段局部范围内存在较大的盈余，其它大部分地区存在较大的亏损.其中，西段的亏损量为1.03×10²⁸ dyne.cm，相当于矩震级为7.9级的地震；中段的亏损量为3.4×10²⁷ dyne.cm，对应的矩震级为7.6级；东段的喜马拉雅构造结存在3.4×10²⁷ dyne.cm的矩亏损，矩震级为7.6级.西藏中南部的拉萨块体、羌塘块体和巴颜喀拉块体内部均存在地震矩亏损，总亏损量为8.9×10²⁷ dyne.cm，对应的矩震级为7.9级，金戈等(2009)的研究认为拉萨附近存在地震活动的可能.在柴达木块体北缘，地震矩亏损量为2.5×10²⁷ dyne.cm(M_W 7.5).阿尔金断裂113年间的地震记录较少，因此该区域存在较大的地震矩亏损，亏损量达到了5.1×10²⁷ dyne.cm(M_W 7.7).西秦岭北缘断裂的亏损量为1.1×10²⁷ dyne.cm(M_W 7.3)，葛伟鹏等(2013)认为该区域现今的形变场是不同活动断裂差异性相对运动、区域内部逆冲挤压和块体旋转共同作用的结果，在西秦岭北缘断裂存在地震危险性.东昆仑断裂存在较大的地震矩盈余，但其东段存在2.5×10²⁶ dyne.cm(M_W 6.9)的能量亏损，李正芳等(2012)通过古地震资料的分析研究表明，东昆仑断裂的塔藏段和玛曲段未来百年发生大震的可能性较高.

南北地震带中南段构造活动比较强烈(张培震等，2008，2013)，其地震矩总亏损量为9.7×10²⁷ dyne.cm，主要分布于鲜水河断裂带南段、安宁河断裂带、则木河断裂带、红河断裂带、龙门山断裂北段和南段等区域.图 6反映了该区域地震矩亏损量的分布，其中，鲜水河断裂亏损量为0.72×10²⁷ dyne.cm，安宁河断裂带的亏损量为1.4×10²⁷ dyne.cm，则木河的亏损量为0.95×10²⁷ dyne.cm，红河断裂带的亏损量为5.2×10²⁷ dyne.cm，龙门山断裂带北、南段的亏损量分别为0.39×10²⁷ dyne.cm和0.3×10²⁷ dyne.cm，对应的矩震级分别为M_W 7.2、M_W 7.4、M_W 7.3、M_W 7.7、M_W 7.0和M_W 6.9.冉勇康等(2008)认为，安宁河断裂北端若还处于集群期，最后一期事件的离逝时间为527a，断裂已进入复发期，未来一段时间有发震的可能.徐锡伟等(2005)认为鲜水河断裂南端和安宁河断裂北端为地震危险区.易桂喜等(2008)认为安宁河断裂冕宁附近、则木河断裂西昌附近和鲜水河断裂中南段将是未来大震或强震的震源区.Wang等(2011)的研究表明，红河断裂带存在约7.4级地震矩亏损.易桂喜等(2013)估算了龙门山南段的视应力，结果表明天全至宝兴段应力水平相对较高，具有发生中-强震的条件.此外，在四川盆地南端，本文的结果中存在约6级左右的地震矩亏损，该亏损量可能跟芦山(M_W 6.6)地震相关.

图 6

Fig. 6

图 6 南北地震带中南段地震矩亏损 Fig. 6 Moment deficit in the middle and southern section of North-South seismic belt

鄂尔多斯块体南缘和东缘的亏损量分别为0.74×10²⁷ dyne.cm和0.34×10²⁷ dyne.cm，对应的矩震级分别为M_W 7.2和M_W 7.0，朱艾斓等(2010)的研究表明鄂尔多斯块体西南缘南东段和南缘自汶川地震以后的地震活动率提高显著而且在其东缘的山西断陷带处于应力增加区，存在发生大震的可能，Wang等(2011)认为南缘的秦岭-大别断裂带为华北和华南块体分界线，该区域的较大的地震矩亏损可以解释为分布于过渡带中的应变能.鄂尔多斯块体北缘和西缘的亏损量分别为0.6×10²⁷ dyne.cm(M_W 7.1)和0.76×10²⁷ dyne.cm(M_W 7.2)，高立新等(2012)的研究表明，鄂尔多斯块体北缘处于低b值高应力积累区域，其西北缘的地震危险趋势进一步加剧.此外，天山的地震矩亏损量为5.1×10²⁷ dyne.cm(M_W 7.7)，东南沿海的地震矩亏损量为0.9×10²⁷dyne.cm(M_W 7.2)，贝加尔裂谷带整体经历拉张构造和较大的剪切作用，其地震活动性较高(刁发启等，2009)，地震矩亏损量为2.7×10²⁷ dyne.cm(M_W 7.6).

本文估算了中国大陆及其邻区113年间的地震矩亏损，为进行研究区域的地震危险性分析和地震中长期预报提供了参考资料.然而，估算结果存在几方面的问题：首先，断裂蠕滑等非震形变极其复杂，因时因地而异(Gratier et al.，2011; 刘冠中等，2013)，本文只对其做了粗糙的估计；其次，本文使用的地震目录时间尺度较短，未能包括完整的地震周期，因此，估算结果存在局限性.

3.1    本文详细搜集整理了1900-2013年间发生在中国大陆及其邻区的地震资料，对其中的47个大震做了分段.然后，利用地震资料和GPS等资料联合反演了中国大陆及其邻区的运动学模型，在此基础上计算了1900年以来中国大陆及其邻区的地震矩亏损，并从中扣除了10%的非震形变.最后，结合前人研究成果讨论了中国大陆及其邻区的地震危险性.

3.2    结果表明：研究区域存在较大地震矩亏损的地区主要包括喜马拉雅断裂带、天山、阿尔金断裂带、柴达木盆地北缘、青藏高原中南部、鲜水河断裂带南段、安宁河断裂带、则木河断裂带、红河断裂带、龙门山断裂带北段和南段、鄂尔多斯块体周缘、东南沿海和贝加尔裂谷带等.前人研究表明上述部分区域未来存在发生强震的可能.然而，由于非震形变的复杂性和地震数据的不完备性，若要更准确的分析研究区域的地震危险性，就需要更准确的估算各个断裂带的非震形变，完善历史地震和古地震资料.