The prolonged quiescence anomaly and the future trends of strong earthquakes in Yunnan region
-
摘要:
从1996年丽江MS7.0地震发生后至2023年8月，云南地区M≥6.7和M≥7.0地震出现了长达27.55 a的超长时间平静现象。本文从更大时空尺度上对该平静异常现象及其与青藏高原和川滇地块的强震活动之间的关系进行了分析和总结，并讨论了云南地区未来地震趋势。结果表明，1887年以来，云南地区5次M≥7.0地震平静期内其周边地震活动环境有很大差异性：第 Ⅰ ，第 Ⅱ 和第 Ⅴ平静期出现在青藏高原M≥7.0地震活跃背景下（即云南地区平静，但青藏高原活跃），第 Ⅲ和第Ⅳ平静期出现在青藏高原M≥7.0地震平静背景下（即云南地区平静，青藏高原也平静）。当前云南地区处于27.55 a超长时间平静期（即第Ⅴ平静期），其所处的周边地震活动环境与第Ⅰ 和第 Ⅱ平静期相似，据此推断，其后续地震趋势可能也与第Ⅰ 和第 Ⅱ活跃期相似，处于相对弱的活跃期。同时，未来较长时间青藏高原可能仍将处于M≥7.0地震活跃时段；2022年泸定MS6.8地震后，川滇地块可能会进入M≥6.7地震活跃时段。
Abstract:From the Lijiang MS7.0 earthquake in 1996 to August 22, 2023, a significant phenomenon of seismic quiescence has lasted for 27.55 years, which are the longest quiet period since 1900, both for M≥6.7 and M≥7.0 earthquakes in Yunnan Province. We further analyze and summarize the quiescence anomaly phenomenon of Yunnan’s strong earthquake and its relationship to the seismic status of the Qinghai-Xizang Plateau and Sichuan-Yunnan block on a larger spatio-temporal scale. Meanwhile, we discuss the future seismic trends in Yunnan region. Here are the conclusions.
1） The energy released by all M≥5.0 earthquakes per year in Yunnan region is equivalent to an M6.7 earthquake, which means an M6.7 represents the mean released energy per year, so we choose M6.7 as the lower limit of strong earthquake to study the seismicity in Yunnan region. The retrospective analysis based on historical earthquakes indicates that the characteristics of active-quiet alternation in Yunnan region, Sichuan-Yunnan block, and QinghaiXizang Plateau, which are closely related to tectonic settings of M≥6.7 or M≥7.0 earthquakes.
Since 1887, there are four active-quiet alternative periods for M≥6.7 earthquakes in Yunnan region, the active periods are 1913−1925, 1936−1955, 1970−1979, 1988−1996 with the mean duration of （11.96±5.35） years, while the quiet periods are 1887−1913, 1925−1936, 1955−1970, 1979−1988 with the mean duration （15.07±7.57） years. As for M≥7.0 earthquakes, the active periods are 1913−1925, 1941−1950, 1970−1976, and 1988−1996 with the mean duration （8.40±2.12） years, while the quiet are 1887−1913, 1925−1941, 1950−1970, 1976−1988 with the mean duration （18.58±5.87） years. The annual rate of occurrence is 0.42−0.83 times per year for M≥6.7 earthquakes and 0.27−0.63 times per year for the M≥7.0 earthquakes in Yunnan region.
Since 1887, there have been two active-quiet alternative periods for M≥6.7 earthquakes in the Sichuan-Yunnan block, and the active periods are 1887−1925, 1948−1996 with the mean duration （42.47±7.38） years, while the quiet periods are 1925−1948, 1996−2022 with the mean duration （24.89±2.40） years. The annual rate of occurrence is 0.16−0.25 times per year for M≥6.7 earthquakes in Sichuan-Yunnan block.
Since 1897, there have been three active-quiet alternative periods for M≥7.0 earthquakes in Qinghai-Xizang Plateau. The active periods are 1904−1956, 1970−1976, and 1988 to now, while the quiet periods are 1897−1904, 1956−1970, 1976−1988 with the mean duration （10.88±3.26） years. The annual occurrence rate is 0.88−1.36 times per year for M≥7.0 earthquakes in Qinghai-Xizang Plateau.
2） Since 1887, there have been significant differences in seismicity around Yunnan region for its five quiet periods of M≥6.7 earthquakes. The quiescence periods Ⅰ , Ⅱ , and Ⅴ appeared under the background of the seismic activity of M≥7.0 earthquakes in QinghaiXizang Plateau, namely, seismic quiescence appeared in Yunnan region while activity in Qinghai-Xizang Plateau. However the quiescence periods Ⅲ and Ⅳ appeared under the background of the quiescence of M≥7 earthquakes in Qinghai-Xizang Plateau, namely both quiescence in Yunnan region and Qinghai-Xizang Plateau. Currently, the period Ⅴ has lasted for 27.55 years in Yunnan region, and the surrounding seismic activity is similar to that in the quiet periods Ⅰ and Ⅱ . Therefore we deduce the following seismic trend may be in relatively weak active periods, similar to the periods Ⅰ and Ⅱ .
3） The comparison of seismicity between the Big Triangle area and Qinghai-Xizang Plateau shows that the active period （1904−1956） of M≥7.0 earthquakes in Qinghai-Xizang Plateau corresponds to the active period （1879−1957） of M≥8.0 earthquakes in Big Triangle area; from 1956 to 1988 there were two quiet periods in Qinghai-Xizang Plateau which correspond to the quiet period （1954−2001） in Big Triangle area, seismicity of M≥7.0 earthquakes in Qinghai-Xizang Plateau has been in an active period since 1988 while the seismicity of M≥8.0 earthquakes in Big Triangle area has come into a new active period since 2001. These seem to indicate that the active periods are consistent, as well as the quiet periods of M≥7.0 earthquakes in Qinghai-Xizang Plateau and M≥8.0 earthquakes in the Big Triangle area. The analogy of the current and the period （1840−1957） of M≥8.0 earthquakes in the Big Triangle area suggests that the current active period of M≥7.0 earthquakes and M≥8.0 earthquakes in the Big Triangle area of East Asian continent will last for decades in the future.
4） On September 5, 2022, an MS6.8 earthquake occurred in Luding, Sichuan Province. The earthquake was located on the Moxi fault, the southeast section of the Xianshuihe fault zone. The quiet period of M≥6.7 earthquakes that lasted for 26.59 years in the SichuanYunnan block since 1996 has been terminated by the Luding MS6.8 earthquake. The following seismicity may be a new active period of M≥6.7 earthquakes in the Sichuan-Yunnan block that analogizes the seismicity of quiet and active periods around 1948.
5） Since 1997, all earthquakes with M≥7.0 in Qinghai-Xizang Plateau have concentrated in the Bayan Har block except for several occurred in the Himalayan belt and Myanmar arc region. The retrospective analysis based on historical earthquakes indicates that the main region of seismic activity is alternative in the Sichuan-Yunnan block and Bayan Har block: Seismic activity appeared in the Bayan Har block while quiescence in the Sichuan-Yunnan block, by contrast, seismic activity appeared in the Sichuan-Yunnan block while less activity in Bayan Har block. Based on the inference, seismicity will be active in the Sichuan-Yunnan block and weak in the Bayan Har block after the Luding MS6.8 earthquake in 2022.
-
引言
高精度相对、绝对重力测量的潮汐改正精度与地球表面潮汐模型的精度密切相关。潮汐模型可利用适当的地球模型理论计算建立，也可利用实测重力时间序列解算建立。实测模型的建立主要取决于力源信号的分离和观测系统的误差改正（Habel，Meurers，2014），而大面积地建立高密度连续重力观测站面临着高成本的挑战（孙和平等，2002），因此，目前重力测量的精密潮汐改正大多采用理论潮汐模型。早在二十世纪初，Love （1909）用一个球对称、无自转的弹性地球体来描述固体潮，给出了全球潮汐因子的平均值为1.16；Molodensky （1961）建立起适合于龙格-库塔（Runge-Kutta）法求解的微分方程组后，将地球勒夫数的计算扩展到29阶，并利用均质、刚性圆球状的地球模型给出了由5个潮波组成的全球理论潮汐模型。此后，Wahr （1981a，b，c）采用无海洋和无气压影响的符合流体静力学原理的均质、刚性旋转椭球状地球模型首次提出潮汐因子可用常数项和纬度依赖项来描述，并根据不同的地球模型给出了3个全球理论潮汐模型，这已成为目前建立全球潮汐模型的主要模式。Dehant等在Wahr模式的基础上建立了流体静力弹性地球（hydrostatic elastic，缩写为He）潮汐模型和流体非静力非弹性地球（non-hydrostatic in-elastic，缩写为NHi）潮汐模型（Dehant，Ducarme，1987；Dehant et al，1999），两模型分别被命名为DDW-He和DDW-NHi模型。此后Mathews （2001）利用甚长基线干涉测量（very long baseline interferometry，缩写为VLBI）修正了初始参考地球模型（preliminary reference Earth model，缩写为PREM）的液核动力学扁率，计算了新的全球理论潮汐模型，该模型简称为M2001。该模型的建立模式与Wahr模式基本一致，只是为了描述潮汐模型的相位，常数项部分使用复数来表示。从Dehant模型开始，潮汐因子的精度从最早的±0.001提高至现今的±0.000 01。
随着全球潮汐模型精度的提高，潮汐观测精度也在提高，且两个模型参数值的相对误差逐步减小，已可利用实测资料对理论潮汐模型进行检验和评估。全球重力潮汐台网主要潮波的潮汐因子（Melchior，de Becker，1983）与Wahr （1981a，b，c）模型值的相对误差为0.060—0.015。我国的重力潮汐剖面（毛慧琴等，1989）参数与Wahr模型值的相对误差为0.030—0.006，我国1989—1993年的连续重力观测台网（魏望生，喻节林，1995）与DDW模型的潮汐因子最小相差约0.002，这得益于当时我国重力潮汐因子的观测精度从±0.005提高至±0.001左右。随着全球海潮负荷改正模型和气压改正的不断完善，超导重力仪格值的精密测定工作持续展开，早期的T型（大容量杜瓦瓶型）、CT型（紧凑型）和CD型（双球型）超导重力仪周日和半日波潮汐因子的观测精度可达到±0.000 4—±0.000 7，为建立以观测为主的全球实验重力潮汐模型（global experimental model for gravity tides，缩写为GEMGT）创造了条件（Xu et al，2004）。近一二十年里，孙和平等（2000）和韦进等（2012）利用现代观象台型超导重力仪的观测结果来检验DDW和M2001两模型，其结果显示，潮汐因子偏差仅为0.002 4—0.000 6，其中M2和O1波的偏差仅0.002 0—0.000 6。上述成果都仅利用了观测结果与模型的偏差来评价重力仪的潮汐观测能力，并未对全球潮汐模型的潮汐改正能力予以评估，且以观测为主的实验模型只用到了中国一个超导重力站（武汉）。该模型是否适合我国大面积相对重力观测数据的潮汐改正，在理论和实验模型中如何选取潮汐改正模型，这些问题依然没有解决。为研究全球理论潮汐模型对我国重力潮汐改正的精度，并指导我国大面积流动重力观测数据的处理，本文拟选取2016—2018年我国北部的10个gPhone连续重力观测站的数据，在评估其潮汐观测精度的基础上，利用观测的和理论的潮汐改正模型参数构建全球潮汐模型的评价指标，提出符合我国相对重力潮汐改正的最优模型，以期为后续建立我国潮汐改正模型提供初值，也为流动重力数据高精度潮汐改正提供改进方案。
1. 全球潮汐模型
1.1 全球潮汐模型简介
为评价全球潮汐模型在中国区域的适应性，本文首先收集了1961年以来的7个全球潮汐模型，列于表1。
表 1 全球潮汐模型的特征Table 1. Feature statistics of global body tide models来源 长周期波 周日波 半日波 1/3日波 地球模型 Molodensky （1961） - 4 1 - 圆球 Warh （1981a，b，c） 1 19 1 0 1066A，PREM-C，C2 Dehant et al （1999） 1 9 1 1 流体静力弹性地球模型 11 9 1 1 非流体非弹性地球模型 Mathews （2001） - 16 - - PREM 二十世纪中叶，Molodensky （1961）和Wahr （1981a）分别利用圆球、1066A、PREM-C、C2等地球模型反演了潮汐模型，并将潮汐因子表达为全球常数项和相对极其小的纬度依赖项之和的形式。Dehant等（1999）以及Mathews （2001）都使用VLBI技术观测地球自由核章动（free core nutation，缩写为FCN）的周期，并通过该周期改正液核动力学扁率，因而建立了更精确的重力潮汐理论模型。至此模型潮汐因子的精度可达到0.000 01。从所收集到的全球潮汐模型的潮波组成来看，Dehant模型虽然在周日波的波群数量上不及Mathews模型，但其覆盖频段最宽。
1.2 理论潮汐改正模型的正演
观测数据量决定了观测潮汐模型波群（以下简称为观测波群）的分波方式，地球模型决定了理论潮汐模型波群（以下简称为模型波群）。观测波群与模型波群之间会因频率范围的不同而导致两者进行比较时存在差异。通常一个观测波群的频段包含多个理论波群，所以以观测潮汐模型的分波方式为标准，将多个理论波群的潮汐振幅（Tamura，1987）加权平均作为潮汐改正模型的振幅。M2001模型没有长周期、半日波、1/3日波的波群，本文用DDW-NHi模型的相应波群代替。这样全球潮汐模型正演的潮汐改正模型不仅能在频域内计算重力残差矢量（Melchior，1994；Habel，Meurers，2014），还能计算正演时域内的潮汐改正值。
由于我国超导重力仪较少，弹簧重力仪数量居多，因此在数据处理和分析过程中，还大量使用全球潮汐因子平均值进行潮汐改正。本文对1961年以来的7个理论全球潮汐模型中振幅较大的5个潮波的常数项和纬度依赖项进行了对比（表2），结果表明，从Dehant模型开始，潮汐改正模型的精度提高了两个数量级，而且M2001还提供了潮汐改正模型的相位信息，这势必会为当前丰富的地面重力数据带来更精确的潮汐改正值。
表 2 7个全球潮汐模型潮汐因子的常数项和纬度依赖项Table 2. Constant term and latitude-dependent term of the gravimetric amplitude factor in seven global body tide models全球潮汐模型 O1波 P1波 K1波 Q1波 M2波 ${G_0} $ 
$G_ { \pm }'/10^{-5}$ 
${G_0} $ 
$G_ { \pm }'/10^{-5}$ 
${G_0} $ 
$G_ { \pm }'/10^{-5}$ 
${G_0} $ 
$G_ { \pm }'/10^{-5}$ 
${G_0} $ 
$G_ { \pm }'/10^{-5}$ 
Molodensky 圆球 1.16 0 1.153 0 1.137 0 1.137 0 1.16 0 1066A 1.152 −600 1.147 −600 1.132 −600 1.152 −600 1.16 −500 Warh PREM-C 1.152 −700 1.147 −700 1.132 −600 1.152 −700 1.16 −500 C2 1.151 −700 1.147 −600 1.132 −600 1.151 −700 1.159 −500 Dehant DDW-He 1.154 24 8 1.147 77 −13 1.147 77 −13 1.132 83 −60 1.160 30 7 DDW-NHi 1.154 24 8 1.149 15 −10 1.134 89 −57 1.154 03 9 1.161 72 10 Mathews M2001 1.154 07 5 1.148 91 5 1.136 1 5 1.154 1 5 - - 注： ${G_0} $ 为潮汐因子的全球常数项，
$G_ { \pm }'$ 为纬度依赖项系数，其中M2001模型的
${G_0} $ 为原模型常数项的模，
$G_ { \pm }'$ 为原模型的周期项。
测点周日波的潮汐因子可表示为
$\delta {\text{＝}} {G_0} {\text{＋}} G'_ { {\text{±} } }\frac{{\sqrt 3 }}{{2\sqrt 2 }}(7{\cos ^2}\theta {\text{－}} 3){\text{，}}$
(1) 半日波的潮汐因子可被表示为
$\delta {\text{＝}} {G_0} {\text{＋}} G_ {\text{±}} '\frac{{\sqrt 3 }}{2}(7{\cos ^2}\theta {\text{－}} 1){\text{，}}$
(2) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
2. 重力资料的精度评价和模型的检验
2.1 重力资料处理及其精度评价
2.1.1 数据预处理和潮汐分析
gPhone重力仪采用了零长弹簧技术，采样频率为1 Hz。在进行重力固体潮观测同时台站还提供了气压小时值数据。本研究对2016—2018年我国10台套该型重力仪的重力固体潮观测数据逐年进行降采样、格值系数校正、预处理、海潮负荷改正、联合气压潮汐分析等处理，建立了观测潮汐改正模型。本研究所用10个观测站的位置如图1所示，具体信息列于表3。
![]()
图 1 我国10个重力观测站空间分布示意图Figure 1. The spatial distribution of ten gravity observatories in China表 3 2016年重力潮汐观测站基本信息Table 3. Basic information of the gravity observatories in 2016序号 台站 省份 台站
类型观测环境
等级台站简介 中误差
/（10−8 m·s−2）1 高台 甘肃 地下室 优秀 位于黑河以北不足2 km处 ±0.637 2 格尔木 青海 山洞 优秀 位于昆仑布尔汉达山北麓山洞中 ±0.706 3 海拉尔 内蒙 山洞 优秀 位于海拉尔盆地与大兴安岭北坡海拉尔河与伊敏河交汇处 ±0.601 4 兰州 甘肃 山洞 优秀 位于兰州东盆地盐场堡以北、白塔山以东、黄河北岸Ⅲ级阶级地的后缘，
距黄河最近距离约1.5 km±0.740 5 牡丹江 黑龙江 山洞 良好 距市区约1.6 km，近年在距离台站500 m处建有成片住宅小区 ±0.718 6 沈阳 辽宁 山洞 良好 位于沈阳市东郊东陵区天柱山西南麓，南距沈抚公路北线约0.4 km ±0.747 7 泰安 山东 山洞 良好 位于泰山南麓正南，0.5 km处有小型水库 ±0.989 8 炭山 宁夏 山洞 良好 位于海原县炭山乡丘陵村，距海原县新区约20 km，台站供电系统时常故障 ±0.891 9 乌加河 内蒙 山洞 优秀 位于内蒙古自治区巴彦淖尔市乌拉特中旗乌加河镇正北1 km处 ±0.838 10 乌什 新疆 山洞 优秀 阿克苏地区乌什县西郊的七女坟旁 ±0.605 本文所选重力站大多分布在我国的中高纬度地区，观测场地与潮汐分析中误差存在一定相关性，由表3可见：观测环境等级为优秀的观测站，中误差为±0.60×10−8 m/s2；观测环境等级为良好的台站中误差虽稍大，也都小于±1×10−8 m/s2。10个观测站的中误差平均值为±0.74×10−8 m/s2，其中乌什重力站最小，约为±0.6×10−8 m/s2。
2.1.2 潮汐分析结果及比较
潮汐因子中误差是评价gPhone重力仪重力潮汐观测能力的一个重要指标。本文将2016—2018年期间的重力固体潮观测数据逐年进行调和分析，平均后统计分析振幅大于10×10−8 m/s2的5个潮波的潮汐因子δ和中误差σ，结果列于表4和表5。
表 4 连续重力观测站的潮汐因子δ和气压导纳值的统计结果Table 4. The result of gravimetric amplitude factor δ and barometric admittance of gravity observatories重力站 O1波 P1波 K1波 M2波 S2波 气压导纳/（10−8 m·s−2·hPa−1） δ σ/10−4 δ σ/10−4 δ σ/10−4 δ σ/10−4 δ σ/10−4 数值 中误差 格尔木 1.152 94 4.0 1.145 22 14.3 1.131 84 9.9 1.165 22 1.5 1.165 10 3.4 −0.269 0.015 兰州 1.157 16 4.7 1.148 77 9.0 1.136 50 3.4 1.158 77 2.4 1.160 37 5.4 −0.319 0.011 高台 1.156 37 3.8 1.150 97 6.7 1.137 49 2.5 1.161 62 1.2 1.163 21 2.7 −0.391 0.008 海拉尔 1.155 19 5.9 1.157 16 11.2 1.138 65 4.1 1.167 74 4.2 1.171 67 9.4 −0.371 0.010 牡丹江 1.153 84 4.7 1.151 92 8.4 1.137 44 3.2 1.163 77 3.4 1.167 07 7.5 −0.347 0.009 沈阳 1.153 70 7.3 1.150 00 13.1 1.137 17 4.9 1.159 45 2.8 1.159 77 6.2 −0.356 0.010 泰安 1.158 32 11.3 1.153 09 20.6 1.139 41 7.7 1.155 73 2.4 1.156 20 5.6 −0.288 0.018 炭山 1.154 74 7.8 1.145 53 16.9 1.133 36 6.4 1.161 86 2.4 1.161 45 5.9 −0.312 0.024 乌加河 1.155 54 3.1 1.150 56 5.6 1.137 52 2.1 1.157 09 1.3 1.159 06 3.1 −0.299 0.010 乌什 1.155 14 3.9 1.151 40 7.0 1.137 26 2.6 1.160 97 1.4 1.161 00 3.2 −0.430 0.010 注：δ为潮波的潮汐因子，σ为中误差。 表 5 乌什重力观测站的潮汐分析结果Table 5. The tidal analysis results of Wushi gravity observatory起始频率
/cpd结束频率
/cpd潮波名称 潮汐振幅
/（10−9 m·s−2）潮汐因子 相位滞后/° 数值 中误差 数值 中误差 0.501 37 0.911 39 Q1 58.938 1 1.154 66 0.001 80 −0.097 0.09 0.911 39 0.947 99 O1 307.830 6 1.155 14 0.000 39 0.020 0.019 0.947 99 0.981 85 M1 24.209 7 1.504 93 0.004 20 −0.427 0.160 0.981 86 0.998 63 P1 143.232 9 1.151 40 0.000 70 −0.060 0.035 0.998 63 1.001 37 S1 3.387 2 1.209 77 0.042 58 16.412 2.319 1.001 37 1.004 11 K1 432.929 1 1.137 26 0.000 26 0.007 0.014 1.004 11 1.006 85 ψ1 3.387 1 1.260 70 0.030 07 −0.923 1.368 1.006 85 1.023 62 φ1 6.164 9 1.169 38 0.016 14 −0.622 0.791 1.023 62 1.057 49 J1 24.208 8 1.153 00 0.004 61 −0.174 0.229 1.057 49 1.470 24 OO1 13.246 6 1.162 69 0.011 38 0.117 0.561 1.470 24 1.880 26 2N2 13.001 1 1.161 73 0.003 27 0.201 0.161 1.880 27 1.914 13 N2 81.411 6 1.159 33 0.000 71 −0.046 0.035 1.914 13 1.950 42 M2 425.214 3 1.160 97 0.000 14 0.013 0.007 1.950 42 1.984 28 L2 12.018 8 1.157 81 0.006 06 −0.428 0.300 1.984 28 2.002 74 S2 197.832 1 1.161 00 0.000 32 −0.289 0.020 2.002 74 2.451 94 K2 53.780 9 1.158 76 0.001 54 −0.137 0.076 2.451 94 7.000 00 M3M6 6.287 4 1.073 14 0.006 61 0.109 0.353 由表4可见：10套仪器中除泰安和海拉尔的部分波外，其它波的潮汐因子中误差均在±0.000 4—±0.000 7附近。潮汐观测精度虽不及OSG型超导重力仪的精度，但已经与二十世纪八九十年代的超导重力仪结果相当。只是gPhone重力仪测定的气压导纳平均值为（−0.319±0.028）×10−8 m·s−2·hPa−1，精度较超导重力仪结果低一个量级（表6）。此外，10个重力观测站中，乌什重力站的潮汐分析结果精度最高，其潮汐分析结果（表5）呈现潮波振幅大、精度高和振幅小、精度低的特征，其中振幅最大的M2波的中误差可达±0.000 14，振幅最小的ψ1波仅为±0.003 27。这表明，最优gPhone重力仪（乌什重力站）的潮汐观测精度虽只有OSG型超导重力仪观测精度的0.5倍，然已高于早期超导重力仪0.5个数量级（表6）。
表 6 各型重力仪潮汐观测精度和稳定度分析Table 6. Accuracy and stability analysis of tidal observatory for different gravimeters重力仪
类型数据时段 M2波潮汐因子
中误差中误差
/（10−8 m·s−2）导纳中误差
/（10−8 m·s−2·hPa−1）稳定度 来源 LRC 1973—1990 0.002 00—0.005 00 - - - Melchior和de Becher （1983） GEO 1971—1990 0.001 00—0.002 00 - - - GWR 1975—1988 0.000 08 - - - LRC 1983—1986 ≤0.005 00 - - 0.005 0 毛慧琴等（1989） GS 1989—1993 0.000 30—0.003 00 - - - 魏望生和喻节林（1995） DZW 0.000 80—0.004 00 - - - LRC 0.000 40—0.002 00 - - - GEO 0.001 00—0.005 00 - - - T，CT，CD 1989—2000 - ±0.05—±0.70 ±0.000 4—0.003 0 0.001 5 Ducarme et al （2002） ASK - - ±0.010 0 - LRC 2000—2001 0.000 51 ±0.20 ±0.022 0 - 孙和平等（2002） OSG 0.000 06 ±0.02 ±0.002 0 - T，CT，CD 1997—2001 0.000 70 - ±0.000 4—0.003 0 0.001 0 Xu et al （2004） OSG 2009—2010 0.000 04 ±1.10—±1.40 ±0.000 5 - 韦进等（2012） 2.1.3 潮汐观测的稳定性检验
重力仪的潮汐观测能力还可以通过潮汐稳定度（Ducarme et al，2002）来评估。本文将10套仪器2016—2018年各潮波的年潮汐因子与3年均值的差作为潮汐稳定度（下文简称为稳定度）来评估重力仪的潮汐观测能力。
图2统计了10个重力站4个潮波的稳定度，可见：除其中7个观测站的P1波（图2c）外，其它潮波的稳定度均小于0.001 5，这一结果已经远高于流动重力仪组网观测的重力潮汐剖面各站点的稳定度（≤0.005 0）（毛慧琴等，1989）；7个重力站4个潮波的稳定度均达到了早期超导重力仪的观测水平（≤0.001 5）。
![]()
图 2 2016—2018年主要潮波的潮汐因子稳定度变化（a） O1波；（b） M2波；（c） P1波；（d） K1波Figure 2. Variation of the stability of gravimetric amplitude factor of the main tidal waves from 2016 to 2018（a） O1 wave；（b） M2 wave；（c） P1 wave；（d） K1 wave表6给出了二十世纪七十年以来重力仪潮汐观测能力的指标统计表。可见：gPhone重力仪周日和半日波的潮汐观测精度优于二十世纪八九十年代的LCR （LaCoste&Romberg）、GEO （大地型）、DZW （中国地震局地震研究所）、GS （阿斯卡尼亚公司一款重力仪）、ASK （阿斯卡尼亚公司一款重力仪）等各型弹簧重力仪约0.5—1个数量级，其中精度高的与早期的T型（大容量杜瓦瓶型）、CT型（紧凑型）和CD型（双球型）的相当，但与现代OSG型超导重力仪相比要低0.5—1个数量级。气压负荷改正参数的精度方面，gPhone可以解算出精度为±0.009×10−8—±0.028×10−8 m·s−2·hPa−1的气压导纳值，较超导重力仪结果低一个数量级。
全球潮汐模型关注的频段是gPhone重力仪的优势频段（周日和半日波），且Xu等（2004）在利用全球超导重力仪建立全球实验潮汐模型时就指出，潮汐因子中误差优于0.000 7即可进行全球潮汐模型的建立工作。本文据此选取重力仪来评估全球潮汐模型的精度。
2.2 观测和全球潮汐模型的比较
上述10个重力站的潮汐观测精度均达到了建立实验潮汐模型的标准，因此利用2016—2018年各年份的观测模型和表2中所列的7个全球潮汐模型进行比较分析，建立两者间的相对误差、单一波群的均方根（root mean square，缩写为RMS）、所有波群的和方根（root sum square，缩写为RSS）和纬度依赖关系等指标进行全球潮汐模型的精度评定。
2.2.1 基本方法
gPhone重力仪的优势频段为周日和半日波，且潮波振幅越大、精度越高。本研究采用均方根（RMS）来衡量理论模型与重力网观测模型的某潮波偏差和方根RSS来综合评估全球潮汐模型的精度（李大炜等，2012），具体表达式如下：
$ {\rm{RMS}} {\text{＝}} {\left\{ {\frac{1}{K}\sum\limits_{n {\text{＝}} 1}^K {[{A_0^2}{{\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!}}{{\!\!\!\!{\text{（}}\!{\delta _0}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!\cos {\varphi _0} {\text{－}} {\delta _m}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!\cos {\varphi _m}\!{\text{）}}\!\!\!\!}^2} {\text{＋}} {A_0^2}{{\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!}}{{\!\!\!\!{\text{（}}\!{\delta _0}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!\sin {\varphi _0} {\text{－}} {\delta _m}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!\sin {\varphi _m}\!{\text{）}}\!\!\!\!}^2}]} } \right\}^{1/2}}{\text{，}} $
(3) $\begin{split} {\rm{RSS}} {\text{＝}}& \left\{ {\frac{1}{K}\sum\limits_{i {\text{＝}} 1}^I {\sum\limits_{n {\text{＝}} 1}^K {\left[ {{{\!\!\!\!{\text{（}}\!{\delta _{0{\text{，}}\!\!\!\!\!i}}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!{A_{0{\text{，}}\!\!\!\!\!i}}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!\cos {\varphi _{0{\text{，}}\!\!\!\!\!i}} {\text{－}} {\delta _{m{\text{，}}\!\!\!\!\!i}}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!{A_{0{\text{，}}\!\!\!\!\!i}}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!\cos {\varphi _{m{\text{，}}\!\!\!\!\!i}}\!{\text{）}}\!\!\!\!}^2} {\text{＋}} } \right.} } } \right.\\& {\left. {\left. {{{\!\!\!\!{\text{（}}\!{\delta _{0{\text{，}}\!\!\!\!\!i}}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!{A_{0{\text{，}}\!\!\!\!\!i}}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!\sin {\varphi _{0{\text{，}}\!\!\!\!\!i}} {\text{－}} {\delta _{m{\text{，}}\!\!\!\!\!i}}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!{A_{0{\text{，}}\!\!\!\!\!i}}\!\!\!\!{\text{（}}\!n\!{\text{）}}\!\!\!\!\sin {\varphi _{m{\text{，}}\!\!\!\!\!i}}\!{\text{）}}\!\!\!\!}^2}} \right]} \right\}^{1/2}}{\text{，}} \end{split}$
(4) 式中：δ0，A0和φ0分别为重力站的潮汐因子、理论振幅、理论相位（通常为0°）；δm和φm分别为模型在重力站点经过气压、海潮改正后的潮汐因子和相位，K为重力站的个数，n为重力站的序号，m表示观测结果，I为潮波分波数，本文为5 （M2，K1，O1，P1，Q1，共计5个）。理论潮汐振幅分别从Tamura （1987）的潮波表和郗钦文的潮波表（郗钦文，侯天航，1987）中选一。表7为利用7个理论模型计算的RMS和RSS。由表7可见：除Molodensky均质圆球模型外，考虑地球扁率的M2001和DDW-NHi模型的RSS更小，其中M2001结果与观测数据吻合得更好，这表明考虑地球扁率的模型更能反映我国潮汐的观测规律，采用郗钦文潮波表（郗钦文，侯天航，1987）计算的RSS优于采用Tamura潮汐表所得的RSS，由此可见DDW-NHi和M2001模型的精度高于其它模型。
表 7 7个全球潮汐模型计算所得重力残差的均方根RMS及和方根RSSTable 7. RMS and RSS of the residual gravity corrected by seven global body tide models潮汐模型 地球模型 RMS/（10−8 m·s−2） RSS/（10−8 m·s−2） M2波 K1波 O1波 P1波 Q1波 Tamura
潮波表郗钦文
潮波表Molodensky 均质圆球 0.143 9 0.168 7 0.173 4 0.064 8 0.047 8 0.292 8 0.291 2 Warh 1066A 0.190 9 0.179 8 0.169 3 0.068 9 0.039 9 0.322 1 0.322 0 PREM-C 0.190 9 0.179 8 0.169 4 0.069 3 0.040 0 0.322 3 0.320 6 C2 0.199 9 0.179 8 0.171 2 0.068 9 0.039 9 0.328 5 0.329 0 DDW DDW-He 0.143 6 0.227 5 0.167 1 0.066 6 0.053 7 0.328 0 0.328 1 DDW-NHi 0.144 5 0.171 3 0.167 1 0.065 2 0.040 2 0.289 8 0.289 6 Mathews M2001 0.144 5 0.168 3 0.167 2 0.065 2 0.040 2 0.288 1 0.287 9 两个全球潮汐模型与观测站潮汐因子的相对误差统计结果（表8）显示，DDW-NHi模型下各重力站的M2，O1，K1和P1波潮汐因子的相对误差均值分别为0.16%，0.17%，0.35%，0.28%，M2001模型下的各潮波潮汐因子的相对误差均值分别为0.15%，0.15%，0.37%，0.18%，这表明DDW-NHi和M2001模型下各重力站各潮波潮汐因子的相对误差为0.15%—0.37%，一致性高，其中高台重力站M2波潮汐因子的偏差在0.000 15附近，达到了早期超导重力仪的潮汐观测稳定度（Xu et al，2004；Habel，Meurers，2014）。从表8还发现，M2001模型的RMS比DDW-NHi更小，这与累计残差振幅结果（表7）一致。
表 8 两个全球潮汐模型与观测站的潮汐因子的相对误差Table 8. The relative errors of the gravimetric amplitude factors between DDW-NHi and M2001 models and the observatories台站 M2波潮汐因子相对误差 O1波潮汐因子相对误差 P1波潮汐因子相对误差 K1波潮汐因子相对误差 RMS/（10−8 m·s−2） DDW-NHi M2001 DDW-NHi M2001 DDW-NHi M2001 DDW-NHi M2001 DDW-NHi M2001 格尔木 0.251 3% 0.251 3% 0.178 2% 0.150 4% 0.253 3% 0.247 4% 0.259 5% 0.459 4% 0.120 6 0.141 0 高台 0.015 2% 0.015 2% 0.224 4% 0.253 1% 0.187 8% 0.189 1% 0.312 9% 0.094 1% 0.087 5 0.068 3 海拉尔 0.015 3% 0.495 8% 0.062 1% 0.093 9% 0.815 2% 0.800 9% 0.574 1% 0.290 9% 0.223 8 0.198 0 兰州 0.015 4% 0.211 4% 0.354 1% 0.381 8% 0.269 4% 0.275 8% 0.291 5% 0.093 8% 0.114 6 0.105 1 牡丹江 0.015 5% 0.200 8% 0.188 0% 0.157 6% 0.582 2% 0.575 4% 0.244 1% 0.008 6% 0.227 8 0.220 6 乌什 0.015 6% 0.236 3% 0.025 6% 0.003 9% 0.168 3% 0.165 9% 0.246 3% 0.011 9% 0.107 8 0.095 5 2.2.2 潮汐振幅的纬度依赖关系
在全球潮汐模型中，潮汐因子通常被表示为全球常数项与纬度依赖项之和，潮汐因子的值与重力站纬度有关。将10个重力站2016—2018年逐年的潮汐参数和3个RSS较小的全球潮汐模型的潮汐参数按纬度排列，如图3所示。
![]()
图 3 O1波，M2波和K1波的潮汐振幅（a）和潮汐因子（b）的纬度依赖关系Figure 3. The latitude dependence of the gravimetric amplitude （a） and gravimetric amplitude factors （b） of O1，M2 and K1 waves由图3a可见，潮汐振幅与纬度有明显的相关性，其中：M2波和O1波的理论振幅与观测振幅差异不大；DDW-NHi以及M2001模型的K1波一致性较好，而与Molodensky模型相差1×10−8—2×10−8 m/s2，所以观测振幅与DDW-NHi或M2001模型一致性更好。
从潮汐因子与纬度的相关性（图3b）可见，本文所选站点的纬度与潮汐因子的相关性不明显，观测值分居模型的两侧，观测值与模型值相差0.000 15—0.005 00。
综合RMS、RSS和潮汐振幅的纬度依赖关系表明，DDW-NHi和M2001模型的RSS值相对较小，处于0.28×10−8—0.29×10−8 m/s2之间（表7）。在潮汐振幅与纬度的相关性中，Molodensky模型的K1波在高纬度地区（35°N—50°N）明显偏离另外两个模型。因此本文初步认定DDW-NHi和M2001模型可被我国相对重力仪用于潮汐改正。
2.2.3 理论和观测潮汐改正模型的比较
我国流动重力观测站数目远多于连续重力观测站，流动重力数据采用的还是Molodensky均质圆形地球确定的潮汐改正模型。为评估该模型和本文筛选模型时域范围内的改正特征，依然将Molodensky模型纳入评价范围。因此本文利用3个理论的和1个观测的潮汐改正模型对观测精度最优的乌什重力站2016年的观测数据（图4a）进行时域改正，比较改正后的重力残差时间序列和功率谱密度，结果如图4所示。
![]()
图 4 潮汐模型潮汐改正后的重力残差时间序列及其功率谱密度（a） 乌什重力站固体潮预处理数据；（b） Molodensky模型潮汐改正的重力残差；（c） DDW-NHi模型潮汐改正的重力残差；（d） M2001模型潮汐改正的重力残差；（e） 观测模型潮汐改正的重力残差；（f） 重力残差的功率谱密度Figure 4. Residual gravity and power spectral density corrected by several tidal models（a） The preprocessing gravity data from Wushi observatory；（b） The residual gravity corrected by Molodensky model；（c） The residual gravity corrected by DDW-NHi model；（d） The residual gravity corrected by M2001 model；（e） The residual gravity corrected by observed model；（f） The power spectral density of the residual gravity由图4可见：观测模型的重力残差最小（图4e），其日均值中误差仅为±0.1×10−8—±0.5×10−8 m/s2，且功率谱密度最大能量不到−110 dB （1 dB＝10lg （10−8 m·s−2）2/Hz），几乎被淹没在残差噪声中（图4f）；3个理论模型中DDW-NHi的重力残差日均值中误差最小，仅为±0.4×10−8—±1.0×10−8 m/s2，优于Molodensky模型的中误差±0.5×10−8—±1.5×10−8 m/s2，也优于M2001模型的中误差±0.7×10−8—±1.4×10−8 m/s2。
对理论潮汐模型的潮汐改正精度进行评价时，RMS和RSS指标未考虑小振幅潮波。而利用它们在时域内进行潮汐改正，则可以分析这些小振幅潮波对潮汐改正的影响。理论模型中，DDW-NHi和M2001模型的周日波频段上O1波，P1波和K1波的残差能量相当（图5a），但与Molodensky模型的K1波存在较大差异，这与潮汐振幅与纬度的相关性结果（图3e）一致。另一方面，Molodensky模型的K1波功率谱密度大于其它两模型，应与该模型K1波的潮汐因子不准确有关。
![]()
图 5 周日波（a）和半日波（b）的重力残差功率谱密度频段特征Figure 5. The power spectral density of the residual gravity for diural （a） and semi-diurnal （b） waves由于M2001模型的M2波模型参数采用的是其它模型的，带入时域改正后存在−90 dB的重力残差信号（图5b），这可能是M2波模型参数相位为0°所致，因此这也是M2001模型的改正结果不及DDW-NHi模型的原因之一。此外，模型与观测的差异应该是模型中未考虑台站观测场地、观测环境等因素的影响所导致。
综上所述，DDW-NHi模型是本文统计的理论模型中最优者，而传统的Molodensky潮汐改正模型的改正精度仅为±0.5×10−8—±1.5×10−8 m/s2。
3. 讨论与结论
在全球的7个潮汐模型中，DDW-NHi模型的潮汐改正结果精度最高，能达到±0.4×10−8—±1.0×10−8 m/s2，而传统的Molodensky模型的潮汐改正结果只能达到±0.5×10−8—±1.5×10−8 m/s2。前文分析主要集中于不同地球模型的改正精度评估，然而除了地球模型之外，影响全球潮汐模型改正精度的因素还包括计算精度时所使用的海潮模型、仪器多年运行观测系统的格值系数变化等。为讨论本文结果的有效性，下面将进一步讨论上述两个因素对全球潮汐模型改正精度评估的影响。
全球潮汐模型仅含地球岩石圈及其内部圈层的影响，而实际观测数据中富含诸如海潮、气压等地球外部圈层的影响。近年来国际上推出多个高精度全球海洋潮汐模型和区域海潮模型，经过验潮站、全球导航卫星系统（global navigation satellite system，缩写为GNSS）、重力等手段的检验获知联合近海海潮模型方能获得更高精度的海潮改正。
3.1 区域海潮改正模型对利用全球模型潮汐改正的影响
本文利用重力残差矢量（Melchior，1994，Hábel，Meurers，2014）分潮波比较Nao.99b全球海潮模型、Nao.99b＋osu.chinasea.2010组合以及Nao.99b＋naoregional.1999组合这3种海潮模型对DDW-NHi模型精度的影响特征，结果如图6所示。由于DDW-NHi模型并不包含相位的影响，因此在进行重力固体潮频域改正时也不考虑相位。
![]()
图 6 海潮负荷改正前后观测的与模型的重力残差振幅比较图中P代表进行了气压改正；＋N代表进行了Nao.99b海潮负荷改正；＋NC代表进行了Nao.99b和osu.chinasea.2010海潮负荷改正；＋NN代表进行了Nao.99b和naoregional.1999海潮负荷改正Figure 6. Comparison of the residual gravity amplitude before and after corrected by ocean loading models and global body tide modelP represents the gravimetric amplitude factors corrected by atmospheric pressure；＋N represents those corrected by global ocean model of Nao.99b；＋NC represent the those corrected by the composite ocean model of Nao. 99b and osu.chinasea.2010；＋NN represent those corrected by the composite ocean model of Nao.99b and naoregional.1999未进行海潮改正时（图6中的红色线柱），各潮波重力残差振幅和各潮波累加值在0.1×10−8—2×10−8 m/s2之间，并呈现东部振幅大、西部小的特点；东部振幅最大为2×10−8 m/s2，西部最小为0.5×10−8 m/s2。海潮负荷改正后，所有潮波的重力残差振幅显著减小，分潮波的海潮负荷改正量最大为1.5×10−8 m/s2，最小只有0.3×10−8 m/s2，改正量为改正前重力残差振幅的70%—80%，这一结果与其它内陆台站的改正量（周江存等，2009）基本一致。由于海潮负荷改正量级已经远大于表7中单一潮波的重力残差振幅，因此海潮负荷改正是本文必不可少的工作。
由图6还可以看到，不同组合模型之间的重力残差振幅差异不足0.1×10−8 m/s2，该量级小于表7中单一潮波重力残差振幅，因此本文认为区域海潮模型的差异不会对全球潮汐模型的改正产生影响。
3.2 全球海潮模型对利用全球模型潮汐改正的影响讨论
为进一步讨论全球海潮模型对潮汐改正的影响，本文对海潮负荷时序改正前后的潮汐因子与DDW-NHi模型之间的差异进行对比，结果如图7所示。可见：所有东部台站的海潮负荷改正量均较大，西部台站均较小，这符合我国东部重力站受太平洋海潮影响强于西部的规律。从改正量级看，牡丹江重力站O1波的改正量最大，可达0.020 0，平均值也有0.003 0，该量级与观测系统的潮汐稳定度相当，高于潮汐因子的观测精度，说明海潮负荷改正后的潮汐因子大小均趋于DDW-NHi全球潮汐模型值，海潮负荷改正在研究非常必要。
![]()
图 7 海潮改正对潮汐因子的影响图中P代表进行了气压改正的潮汐因子；DUT10，EOT11a，FES2004，GOT4.7，HAM11a，Nao.99b和TPXO7.2分别代表进行了气压和对应全球海潮模型改正的潮汐因子；DDW-NHi代表了DDW-NHi理论潮汐因子。横坐标台站按经度由东向西排列Figure 7. Impaction of the gravimetric amplitude factors corrected by seven global ocean modelsP represents the gravimetric amplitude factors corrected by atmospheric pressure；DUT10，EOT11a，FES2004，GOT4.7，HAM11a，Nao.99b and TPXO7.2 represent those corrected by atmospheric pressure and corresponding global ocean tide model respectively；DDW-NHi represents those in DDW-NHi. The stations are arranged from east to west according to their longitude比较不同海潮负荷改正模型之间的差异可知，图7中的7个全球海潮负荷改正结果的差异较小，为0.000 5—0.002 0，小于改正前后的潮汐因子差异0.020 0—0.003 0和观测系统稳定度0.001 5，与潮汐因子的观测精度0.000 7相当。因此，在进行全球潮汐模型评估时须进行海潮负荷改正，而不同海潮负荷改正模型之间的差异对结果影响并不明显。
3.3 格值系数对模型检验的影响
相对重力仪格值系数偏差会直接影响观测潮汐改正模型的准确性（Francis et al，1998；Hábel，Meurers，2014），gPhone重力仪的格值系数的范围是0.999 9—1.019 0 （刘子维等，2011；张锐等，2011）。本研究选择青海格尔木重力站2016—2018年逐年观测数据讨论潮汐因子、残差中误差与格值系数的关系，结果如图8所示。
![]()
图 8 格尔木观测站潮汐因子和调和分析中误差随格值系数的变化（a） M2波；（b） O1波；（c） K1波；（d） P1波；（e） Q1波；（f） 中误差随格值系数的变化Figure 8. Varaition of gravimetric amplitude factor and RMS of the harmonic analysis with the scale factor for Golmud observatory（a） M2 wave；（b） O1 wave；（c） K1 wave；（d） P1 wave；（e） Q1 wave；（f） The RMS of the residual gravity以DDW-NHi模型为潮汐改正模型，利用重力残差（图8f）作为约束测定格值系数。从图8a-e可见：观测潮汐因子与理论仅相差0.001—0.002。如果不进行这样的标定（即格值系数为1），潮汐因子偏差可能大于0.005，该量级已远超当前gPhone重力仪的稳定度和潮汐观测精度。因此进行观测系统的格值系数改正是必不可少的。
3.4 结论
本文比较了7个常用的全球潮汐模型，并利用10个安装gPhone重力仪的重力站超过3年的观测资料对模型进行了检验。gPhone重力仪的周日和半日波潮汐观测能力虽然比现代OSG型超导重力仪低0.5—1.0个数量级，但是依然可达到早期超导重力仪的水平，这样的精度可以对全球潮汐模型精度进行评估。对潮汐观测和理论模型的均方根、和方根、潮汐因子相对误差、潮汐因子的纬度依赖关系和重力残差及其功率谱密度的评估结果显示DDW-NHi的全球潮汐模型能够提供最优的重力残差精度为±0.4×10−8—±1.0×10−8 m/s2，比大规模流动重力观测常用的Molodensky模型的精度±0.5×10−8—±1.5×10−8 m/s2高出1×10−8—2×10−8 m/s2。
在讨论海潮负荷和格值系数对潮汐模型精度评估影响时，仅进行全球潮汐模型的改正而未进行海潮和格值系数校正引起的潮汐因子偏差分别可达0.003 0—0.020 0和0.001—0.005，这一量级已经远大于gPhone重力仪的潮汐观测稳定度和潮汐因子的观测精度。不同海潮模型改正的影响量级（0.000 5—0.002 0）与潮汐因子的观测精度相当，因此可忽略其对模型精度评估的影响。
本文所选观测站只有10个，且多分布在我国北方中高纬度地区；无南方台站特别是沿海台站参与本次评估，难免会使得结果仅具有区域性规律。从观测资料来看，南方和沿海地区的潮湿使重力观测系统在长期连续观测过程中的故障率较高，使用现有分析方法和近期数据无法达到评估模型所需的精度。如何筛选此类观测数据、研制新的观测精度评估方法成为解决此类问题的一个重要内容。在进行潮汐改正时，气压负荷影响不容忽视，而并非所有重力站目前均具备同址同采样率的气压观测条件。为弥补站点辅助测量数据的缺失可探讨国际已经发布的全球气压模型的潮汐改正可行性。此外，本文还讨论了格值系数对模型评估的影响，认为需要研制出高效的连续重力观测数据标定方法，使得标定数据达到潮汐模型评估和模型建立所需要的精度。总之，如何科学产出符合全球潮汐模型评估精度的重力潮汐数据是本文工作特别是建立潮汐模型的根本所在。
感谢国家重力台网中心提供10套高精度gPhone重力仪超过3年的重力固体潮数据，感谢审稿专家提出的宝贵意见和建议。
-
![]()
图 1 青藏高原活动地块划分及历史强震震中分布图
（a） 青藏高原M≥7.0地震分布。Ⅰ ：拉萨地块；Ⅱ ：羌塘地块；Ⅲ ： 巴颜喀拉地块；Ⅳ：柴达木地块；Ⅴ：祁连地块；Ⅵ：川滇地块；Ⅶ：滇南地块；Ⅷ：滇西地块；地块边界数据来引自张培震等（2003）。（b） 云南地区M≥6.7地震分布. F1：安宁河—小江断裂带；F2：红河断裂带；F3：澜沧江断裂带；断裂数据引自张国民等（2005）
Figure 1. Distribution of active blocks and strong earthquakes in Qinghai-Xizang Plateau
（a） Distribution of M≥7.0 earthquakes in Qinghai-Xizang Plateau. Ⅰ ：Lhasa block；Ⅱ ：Qiangtang block；Ⅲ ：Bayan Har block；Ⅳ ：Qaidam block；Ⅴ ：Qilian block；Ⅵ ：Sichuan-Yunnan block；Ⅶ ：Southern Yunnan block；Ⅷ：Western Yunnan block. Block data are from Zhang et al （2003）. （b） Distribution of M≥6.7 earthquakes in Yunnan region. F1：Anninghe-Xiaojiang fault zone；F2：Honghe fault zone；F3：Lancangjiang fault zone；the fault data refer from Zhang et al （2005）
![]()
图 2 1887—2022年云南地区M≥6.7和M≥7.0地震M-t图
图中彩色条带中数字为持续时间，下同
Figure 2. M-t plot of M≥6.7 and M≥7.0 earthquakes in Yunnan region in the period of 1887−2022
The numbers in color bar represent the duration，the same below
![]()
图 3 1880—2022年青藏高原M≥7.0地震时序分布
（a） M-t图；（b） 时间间隔dt-t图，tavg为青藏高原M≥7.0平均间隔时间；1std为1倍均方差线；2std为2倍均方差线，下同
Figure 3. Distribution with time series for M≥7.0 earthquakes in Qinghai-Xizang Plateau in 1880−2022
（a） M-t；（b） dt-t， where tavg represent the mean time interval for M≥7.0 earthquake in Qinghai-Xizang Plateau，1std represent the 1 time mean square variance，2std represent the 2 times mean square variance，the same below
![]()
图 4 云南地区5个M≥7.0地震平静期内青藏高原M≥7.0地震震中分布图
Figure 4. Distribution of M≥7.0 earthquakes in Qinghai-Xizang Plateau during five quiescence periods of M≥7.0 earthquakes in Yunnan region
![]()
图 5 1880—2022年川滇地块M≥6.7地震时序分布图
（a） M-t图；（b） 时间间隔dt-t图
Figure 5. Distribution of M≥6.7 earthquake with time in Sichuan-Yunnan block （1880−2022）
（a） M-t；（b） Time interval dt-t
![]()
图 6 1880—2022年巴颜喀拉地块和川滇地块M≥6.7地震的分布
（a） 震中分布图；（b） 时序分布图M-t图，①，②，③，④为划分的4个时段
Figure 6. Distribution of M≥6.7 earthquakes in Bayan Har block and Sichuan-Yunnan block in 1880−2022
（a） Distribution of epicenters；（b） M-t plot，where ①，②，③ and ④ represent four divided periods
![]()
图 7 1800—2022年大三角区内M≥8.0地震震中分布图（a）、M-t图（b）和时间间隔dt-t图（c）
Figure 7. Distribution of M≥8.0 earthquakes （1800−2022） in Big Triangle area of East Asian continent （a），M-t （b），time interval dt-t （c）
表 1 云南地区M≥6.7和M≥7.0地震活跃-平静分期特征
Table 1 Characteristics of active-quiet periods of earthquakes with M≥6.7 and M≥7.0 in Yunnan region
序号 发震日期
年-月-日震中位置 震级 震中
位置间隔时间/a 活动分期 活动特征 北纬/° 东经/° M≥6.7 M≥7.0 M≥6.7 M≥7.0 0 1 887-12-16 23.70 102.50 7.0 石屏 第Ⅰ平静期 t＝26.01 a
n＝0
r＝0t＝26.9 a
n＝0
r＝01 1 913-12-21 24.15 102.45 7.0 峨山 26.01 t＝12.00 a t＝11.24 a 2 1 917-07-31 28.00 104.00 6.7 大关 3.61 3 1 923-06-22 22.75 98.75 7.3 缅甸 5.89 9.51 第Ⅰ活跃期 n＝5 n＝3 4 1 925-03-16 25.70 100.40 7.0 大理 1.73 1.73 r＝0.42 r＝0.42 5 1 925-12-22 21.00 101.50 6.8 老挝 0.77 第Ⅱ平静期 t＝10.35 a
n＝0
r＝0t＝16.17 a
n＝0
r＝06 1 936-04-27
1 936-05-1628.90
28.50103.60
103.606.7
6.7马边
马边10.35
0.057 1 941-05-16 23.60 99.40 7.0 耿马 5.00 16.17 t＝19.41 a t＝8.71 a 8 1 941-12-26 22.70 99.90 7.0 澜沧 0.61 0.61 第Ⅱ活跃期 n＝9 n＝3 9 1 942-02-01 23.10 100.30 6.7 思茅 0.10 r＝0.46 r＝0.34 10 1 950-02-03
1 950-02-0321.70
22.10100.10
99.907.0
6.7勐海
澜沧8.01
08.10
11 1 952-09-30 28.30 102.20 6.7 冕宁 2.66 12 1 955-09-23 26.60 101.80 6.7 永仁 2.98 第Ⅲ平静期 t＝14.28 a
n＝0
r＝0t＝19.92 a
n＝0
r＝013 1 970-01-05 24.20 102.68 7.8 通海 14.28 19.92 t＝9.20 a t＝6.40 a 14 1 971-04-28 23.00 101.10 6.7 思茅 1.31 n＝7 n＝4 15 1 974-05-11 28.20 104.10 7.1 大关 3.04 4.35 第Ⅲ活跃期 r＝0.76 r＝0.63 16 1 976-05-29
1 976-05-2924.50
24.6099.00
98.707.3/
7.4龙陵
龙陵2.05
02.05
017 1 976-11-07 27.60 101.10 6.7 盐源 0.44 1 8 1 979-03-15 23.20 101.10 6.8 宁洱 2.36 1 9 1 983-06-24 21.40 102.60 6.9 越南 4.27 第Ⅳ平静期 t＝9.64 a
n＝1
r＝0.10t＝12.22 a
n＝0
r＝020 1 988-11-06
1 988-11-0622.92
23.1699.79
99.557.4
7.2澜沧
耿马5.37
012.22
0t＝7.24 a t＝7.24 a 21 1 992-04-23
1 992-04-2322.30
22.6099.10
99.006.7
6.9缅甸
缅甸3.46
0第Ⅳ活跃期 n＝6 n＝4 22 1 995-07-12 22.00 99.30 7.3 缅甸 3.22 6.68 r＝0.83 r＝0.55 23 1 996-02-03 27.20 100.30 7.0 丽江 0.56 0.56 第Ⅴ平静期 t＝27.49 a
n＝0
r＝0t＝27.49 a
n＝0
r＝0注：表中t表示持续时间，单位为a；n表示地震次数；r表示年发生率，r＝n/t。 
下载: 导出CSV
表 2 青藏高原、川滇地块、云南地区在活跃期中M≥6.7或M≥7.0地震的年发生率
Table 2 Annual occurrence rate of M≥6.7 or M≥7.0 earthquakes during the active period in Qinghai-Xizang Plateau，Sichuan-Yunnan block and Yunnan region
空间尺度 地震年发生率 M≥6.7 M≥7.0 青藏高原 0.88—1.36 川滇地块 0.16—0.25 云南地区 0.42—0.83 0.27—0.63
下载: 导出CSV
-
邓起东,程绍平,马冀,杜鹏. 2014. 青藏高原地震活动特征及当前地震活动形势[J]. 地球物理学报,57(7):2025–2042. doi: 10.6038/cjg20140701 Deng Q D,Cheng S P,Ma J,Du P. 2014. Seismic activities and earthquake potential in the Tibetan Plateau[J]. Chinese Journal of Geophysics,57(7):2025–2042 (in Chinese).
傅征祥,刘杰,王晓青,郝平,吕梅梅. 2003. 中国大陆西部及邻区板内8级大地震活动过程和趋势预测初步研究[J]. 地学前缘,10(特刊):112–117. Fu Z X,Liu J,Wang X Q,Hao P,Lü M M. 2003. Research on process and tendency prediction of great earthquakes (MS≥8) in western China continent and its surroundings[J]. Earth Science Frontiers,10(Z1):112–117 (in Chinese).
傅征祥,吕晓健,郝平,邵志刚,陈丹. 2012. 东亚大陆大三角地震区地震活动性概要[M]. 北京:地震出版社:22−27. Fu Z X,Lü X J,Hao P,Shao Z G,Chen D. 2012. Summary of Seismic Activity in the Great Seismic Triangle Area of Eastern Asian Continent[M]. Beijing:Seismological Press:22−27 (in Chinese).
国家地震局震害防御司. 1995. 中国历史强震目录(公元前23世纪—公元1911年)[M]. 北京:地震出版社:1–514. Department of Earthquake Disaster Prevention,State Seismological Bureau. 1995. Catalogue of Strong Earthquakes in China (23rd Century BC−AD 1911)[M]. Beijing:Seismological Press:1–514 (in Chinese).
国家地震局震害防御司. 1999. 中国近代地震目录(公元1912—1990,MS≥4.7)[M]. 北京:地震出版社:1–617. Department of Earthquake Disaster Prevention,State Seismological Bureau. 1999. Catalogue of Strong Earthquakes in China (AD 1912−1990,MS≥4.7)[M]. Beijing:Seismological Press:1−617 (in Chinese).
皇甫岗,秦嘉政. 2006. 云南地区大震活动规律研究[J]. 地震地质,28(1):37–47. doi: 10.3969/j.issn.0253-4967.2006.01.004 Huangfu G,Qin J Z. 2006. Study of the seismicity of strong earthquakes in Yunnan area[J]. Seismology and Geology,28(1):37–47 (in Chinese).
皇甫岗,石绍先,苏有锦. 2000. 20世纪云南地震活动研究[J]. 地震研究,23(1):1–9. doi: 10.3969/j.issn.1000-0666.2000.01.001 Huangfu G,Shi S X,Su Y J. 2000. Study on seismicity in Yunnan in the 20th century[J]. Journal of Seismological Research,23(1):1–9 (in Chinese).
皇甫岗,苏有锦,张建国,邵德盛. 2013. 新世纪滇缅交界大震活动研究[M]. 昆明:云南科技出版社:131−149. Huangfu G,Su Y J,Zhang J G,Shao D S. 2013. Research on Large Earthquake Activities in Yunnan and Myanmar Border in the New Century[M]. Kunming:Yunnan Science and Technology Press:131−149 (in Chinese).
皇甫岗,苏有锦,张建国. 2014. 缅甸孟帕亚7.2级地震及其对云南强震活动趋势的意义[J]. 地震地质,36(3):598–608. doi: 10.3969/j.issn.0253-4967.2014.03.005 Huangfu G,Su Y J,Zhang J G. 2014. The Mong Hpayak M7.2 earthquake and the significance for the tendency of large earthquake in Yunnan[J]. Seismology and Geology,36(3):598–608 (in Chinese).
刘祖荫,苏有锦,秦嘉政,李中华,张俊伟. 2002. 20世纪云南地震活动[M]. 北京:地震出版社:10−19. Liu Z Y,Su Y J,Qin J Z,Li Z H,Zhang J W. 2002. Seismic Activity in Yunnan Region in the 20th Century[M]. Beijing:Seismological Press:10−19 (in Chinese).
秦嘉政,钱晓东,刘丽芳. 2009. 云南强震活动间隔时间特征研究[J]. 地震研究, 32 (3):221−227. Qin J Z,Qian X D,Liu L F. 2009. Study on interval time characteristics of strong earthquake activities in Yunnan region[J]. Journal of Seismological Research, 32 (3):221−227 (in Chinese).
宋治平,张国民,刘杰,尹继尧,薛艳,宋先月. 2011. 全球地震灾害信息目录(9999 BC—AD 2010)[M]. 北京:地震出版社:17−448. Song Z P,Zhang G M,Liu J,Yin J Y,Xue Y,Song X Y. 2011. Disaster Information Catalog of Global Earthquakes (9999 BC−AD 2010)[M]. Beijing:Seismological Press:17−448 (in Chinese).
苏有锦,秦嘉政. 2001. 川滇地区强地震活动与区域新构造运动的关系[J]. 中国地震,17(1):24–34. doi: 10.3969/j.issn.1001-4683.2001.01.004 Su Y J,Qin J Z. 2001. Strong earthquake activity and relation to regional neotectonic movement in Sichuan-Yunnan region[J]. Earthquake Research in China,17(1):24–34 (in Chinese).
苏有锦,李忠华,刘祖荫,蔡民军. 2001. 20世纪云南地区MS≥5.0级地震活动的基本特征[J]. 地震研究,24(1):1–9. doi: 10.3969/j.issn.1000-0666.2001.01.001 Su Y J,Li Z H,Liu Z Y,Cai M J. 2001. Basic characteristics of seismicity with MS≥5.0 in Yunnan in the 20th century[J]. Journal of Seismological Research,24(1):1–9 (in Chinese).
苏有锦,李忠华. 2011. 云南地区6级以上强震时间分布特征及其概率预测模型研究[J]. 地震研究,34(1):1–7. doi: 10.3969/j.issn.1000-0666.2011.01.001 Su Y J,Li Z H. 2011. Interval distribution and probability model of the strong earthquakes with MS≥6.0 in Yunnan[J]. Journal of Seismological Research,34(1):1–7 (in Chinese).
吴中海,赵根模,龙长兴,周春景,范桃园. 2014. 青藏高原东南缘现今大震活动特征及其趋势:活动构造体系角度的初步分析结果[J]. 地质学报, 88 (8):1401−1416. Wu Z H,Zhao G M,Long C X,Zhou C J,Fan T Y. 2014. The seismic hazard assessment around south-east area of Qinghai-Xizang Plateau:A preliminary results from active tectonics system analysis[J]. Acta Geologica Sinca, 88 (8):1401−1416 (in Chinese).
张国民,马宏生,王辉,王新岭. 2005. 中国大陆活动地块边界带与强震活动[J]. 地球物理学报,48(3):602–610. doi: 10.3321/j.issn:0001-5733.2005.03.018 Zhang G M,Ma H S,Wang H,Wang X L. 2005. Boundaries between active tectonic blocks and strong earthquakes in the China mainland[J]. Chinese Journal of Geophysics,48(3):602–610 (in Chinese). doi: 10.1002/cjg2.693
张培震,邓起东,张国民,马瑾,甘卫军,闵伟,毛凤英,王琪. 2003. 中国大陆的强震活动与活动地块[J]. 中国科学(D辑),33(S1):12–20. Zhang P Z,Deng Q D,Zhang G M,Ma J,Gan W J,Min W,Mao F Y,Wang Q. 2003. Active tectonic blocks and strong earthquakes in the continent of China[J]. Science in China:Series D,46(S2):13–24. doi: 10.1360/03dz0002
赵小艳,付虹,王强. 2008. 汶川8.0级地震后云南地震趋势研究[J]. 地震研究, 31 (增刊):453−457. Zhao X Y,Fu H,Wang Q. 2008. Earthquake tendency in Yunnan after the Wenchuan MS8.0 earthquake[J]. Journal of Seismological Research, 31 (S):453−457 (in Chinese).
An Y R,Wang,D,Ma Q,Xu Y R,Li Y,Zhang Y Y,Liu Z M,Huang C M,Su J R,Li J L,Li M X,Chen W K,Wan Z F,Kang D J,Wang B S. 2023. Preliminary report of the September 5,2022 MS6.8 Luding earthquake,Sichuan,China[J]. Earthq Res Adv,3(1):100184. doi: 10.1016/j.eqrea.2022.100184
Bilham R,Gaur V K,Molnar P. 2001. Himalayan seismic hazard[J]. Science,293(5534):1442–1444. doi: 10.1126/science.1062584
计量
- 文章访问数: 237
- HTML全文浏览量: 53
- PDF下载量: 58
