Statistical analysis of microwave radiation anomaly before earthquake: A case study of Kamchatka Peninsula
-
摘要: 岩石加载实验证实受载岩体能够产生特定频段的微波辐射能量增强，卫星遥感数据分析也表明，个别震例强震前可能存在微波辐射异常，这一重要发现亟待普遍性检验与统计特征研究. 本文收集了俄罗斯堪察加半岛2003—2011年Aqua AMSR-E 的微波辐射亮温数据，并对地震和火山活动信息进行了统计，进而利用改进后的异常识别方法及构建的背景场对9年内多次强震的震前短临异常的位置和时间信息进行了检测与统计. 结果表明：① 勘察加半岛平均亮温变化存在明显的季节特征，1—4月和11—12月为平稳期、5—10月为剧变期，平稳期的亮温曲线基本保持平稳，而剧变期呈抛物线型，自4月起持续增温至7月，随后持续降温至11月；② 受地形地貌影响，半岛南部和东南部亮温较高，中西部亮温较低；③ 火山点亮温小于地震点和非地震点，火山活动对亮温影响有限；④ 使用四分位数法判断辐射异常检测指数，发现17次地震中有13次可以识别到震前异常，其中MS≥6.0地震均可检测到异常，且异常发生在地震前1个月以内.Abstract: It was confirmed that microwave energy is emitted during the compression of rock and the energy can be detected by some special frequency bands. The analysis of satellite remote sensing data found that there may be abnormal microwave radiation before several earthquakes. This important discovery is in urgent need of universal inspection and statistical analysis.In this study, the data of Aqua AMSR-E microwave radiation brightness temperature of nine years (2003−2011) in Kamchatka Peninsula were collected, and the information of earthquake and volcanic activity in this period was counted. The time and location information of microwave radiation anomalies before all the strong earthquakes over the past nine years were studied by using the improved method of anomaly identification and the construction of the background field. The results show that: ① The average brightness temperature change of Kamchatka Peninsula has obvious seasonal characteristics, which can be divided into two periods, from January to April and from November to December are stationary periods and from May to October is fluctuation period. The brightness temperature curves of the stationary periods are basically steady and that of the fluctuation period is parabolic, for it shows continual warming from April to July and then continual cooling until November. ② Due to the effect of topography, the brightness temperature is higher in the southern and southeastern peninsula, while that of the central and western part is low. ③ The brightness temperature of volcanoes is less than that of the seismic and non-seismic points, which indicates that the volcanic activity has a limited influence on the brightness temperature. ④ Using the four percentile method to determine the radiation anomaly detection index, it was found that there are 13 times of 17 earthquakes can be identified anomalies before earthquakes, and allMS≥6.0 earthquakes can be detected anomalies, which occurred in one month before the earthquake.
-
Keywords:
- before earthquake /
- microwave radiation /
- Kamchatka Peninsula /
- AMSR-E /
- stati-stical analysis
-
引言
自20世纪90年代起，中国地震局地球物理研究所、中国地震局地质研究所、中国科学院遥感应用研究所和中国矿业大学等单位基于受载岩石破裂过程中的热红外辐射效应实验开展了岩石加载微波辐射观测实验研究. 实验结果表明，受载岩体能够产生特定频段的微波辐射能量增强，其岩体破裂辐射规律具有明显的前兆特征，能够反映岩体的变形破裂过程（邓明德等，1995；耿乃光等，1995；房宗绯等，2000；王恩元等，2011）. 同时，日本研究人员也通过岩石单轴压缩破裂实验观测到了300 MHz，2 GHz，22 GHz频段的微波辐射异常（Maki et al，2006 ）. 这是继岩石破裂红外辐射异常现象之后，又一新的岩石电磁辐射发现，为地震岩石破裂的微波遥感监测提供了理论依据.
美国地球观测系统（Earth observing system，简写为EOS） Aqua卫星搭载有先进的微波扫描辐射计（the advanced microwave scanning radiometer for the Earth observing system，简写为AMSR-E），可获取双极化12通道（6.925，10.65，18.7，23.8，36.5，89 GHz）的微波辐射亮温数据（王永前等，2011）. 近年来，国内外研究人员使用AMSR-E的18.7 GHz微波亮温数据检测到了地震孕育过程中岩石挤压破裂产生的微波能量辐射异常. 例如：Maeda和Takano （2010）使用AMSR-E的18.7 GHz双极化微波辐射数据，在数据重建之后提出一种基于相邻像元亮温差计算异常指数的算法，且此前已使用该算法在2008年汶川MS7.9地震后1天在震中附近检测到异常增强信号（Maeda，Takano，2009）；陈昊和金亚秋（2010）在该算法基础上进一步定义了辐射异常指数（radiation anomaly index，简写为RAI）和极化指数，并对2010年4月14日玉树MS7.1地震孕育过程中产生的岩石破裂辐射异常进行初步检测，同时引入23.8 GHz频段数据对地表温度和大气水汽在检测过程的影响进行评估. 上述研究均表明个别震例在震前可能存在微波辐射异常，但这一结论并未得到普遍性检验. 虽可检测到微波辐射异常增强，但并未对异常出现的时间和位置进行统计分析. 鉴于此，本文使用改进后的异常检测方法对俄罗斯堪察加半岛17次MS≥5.0地震进行微波辐射异常识别与统计分析，拟在证实强震前存在微波辐射异常这一普遍现象的基础上，进一步探索并试图总结微波辐射异常时空特征的基本规律.
1. 研究区域概况
堪察加半岛位于俄罗斯的远东地区、欧亚大陆的尽头，西濒鄂霍茨克海，东临太平洋和白令海（图1）. 整个半岛面积约为47.23万km2，从东北向西南延伸超过1 200 km. 在地质构造上，堪察加半岛属于新生代阿尔卑斯褶皱带，处于太平洋、欧亚和北美三大板块交汇处，主要受太平洋板块沿NNW方向向鄂霍茨克板块下方的俯冲作用. 板块之间的俯冲和挤压在堪察加地区形成一条典型的板块边界地震带，该地震带基本平行于堪察加半岛的中南部（孟国杰等，2009）. 堪察加半岛俯冲地区分为南、北两部分，南部大约位于50°N—55°N，且在55°N与阿留申海沟交汇，俯冲方位角为55°，而北部的俯冲角由55°降为约35° （Gorbatov et al，1997 ；Kozhurin et al，2006 ）.
本文选取堪察加半岛为研究区域主要基于两个方面：一方面搭载了AMSR-E扫描仪的Aqua卫星在高纬度地区覆盖度较好，而在低纬度地区数据较少；另一方面，堪察加半岛的地震和火山活动频繁，是理想的统计分析试验区域. 根据USGS （2011）的统计数据，堪察加半岛自2003年1月至2011年9月共发生MS≥5.0地震17次，其中MS≥6.0地震4次（图1）.
![]()
图 1 2003—2011年堪察加半岛MS≥5.0地震分布Figure 1. Distribution of MS≥5.0 earthquakes in Kamchatka Peninsula from 2003 to 20112. AMSR-E数据与研究方法
2.1 AMSR-E微波数据
AMSR-E是搭载在美国地球观测系统Aqua卫星上的被动微波辐射计，具有6频率双极化探测的特性，过境时间为当地时间13:30和1:30，目的在于提供全球实时的地球物理学和大气参数数据以支持全球气候变化和水汽循环的科学研究和监测工作（Kawanishi et al，2003 ; Shibata et al，2003 ；王永前等，2011）. 表1列出了AMSR-E辐射计各类设计参数和运行模式，包括频率、瞬间视场和扫描入射角等. 通过美国国家冰雪数据中心（National Snow and Ice Data Center，2011）对2003年1月至2011年9月约9年的AMSR-E L3级数据进行搜索并下载.Maeda和Takano （2008） 的研究表明，AMSR-E的18.7，23.8 GHz通道均接近于岩石破裂所释放微波能量的22 GHz频段，但是由于23.8 GHz频段处于大气水汽吸收峰，而大气水汽是影响微波亮温的重要因素之一，因此本文选取18.7 GHz频段作为主要检测频段.
表 1 AMSR-E辐射计的主要性能参数Table 1. Main parameters and characteristics of AMSR-E中心频率/GHz 带宽/MHz 动态范围/K 瞬间视场/km 主波束效率 入射角/° 天底偏角/° 扫描间距/km 扫描方式 6.925 350 2.7—340 43×75 95.1% 55.0 47.5 10.1×9.0 以40 r/min的
速度呈圆锥
形扫描10.650 100 2.7—340 29×51 94.8% 55.0 47.5 10.1×9.0 18.70 200 2.7—340 16×27 95.8% 55.0 47.5 10.1×9.0 23.80 400 2.7—340 18×32 94.8% 55.0 47.5 10.1×9.0 36.50 1 000 2.7—340 8.2×14 93.9% 55.0 47.5 10.1×9.0 89.00（A） 3 000 2.7—340 3.7×6.5 94.5% 55.0 47.5 10.1×9.0 89.00（B） 3 000 2.7—340 3.5×5.9 93.7% 54.5 47.0 4.5×4.4 2.2 异常识别方法
Maeda和Takano （2009）等在证实了岩石破裂可产生微波辐射的基础上，对AMSR-E数据进行重建和重采样，并提出了一种基于相邻点微波亮温差异的异常分析方法（Maeda，Takano，2010），该方法在摩洛哥和汶川等地震实例的应用中取得了较好的异常分析结果. 本文在该方法的基础上对其进行改进，并利用统计学中的四分位数法构建了适用于堪察加半岛地区的背景场.
空间自相关性，即同一变量的不同观测值之间因观测点在空间上相邻而形成的相关性，在影像上表现为相邻像素的微波亮温值之间存在局部空间自相关性. 基于该原理，本文对日本研究人员相邻两点间微波亮温差异的异常分析方法进行改进，参考遥感锐化中的拉普拉斯算子，将相邻两像素点之间的比较改进为将像元F与周围4个像元P进行综合比较，以便更加有效地提取由于地震岩石破裂所释放的微波辐射增强.
在对数据进行去除背景、辐射校正和数据频段选取等处理后，针对每次地震，选取以震中像元F为中心的5×5像元矩阵为研究区域，并对25个像元进行1—25编码，13号点则表示震中（图2），然后以震前60天震后10天为研究时段，对2003—2011年统计得到的震例逐次进行分析.
![]()
图 2 异常识别方法的研究区域选取与计算改进Figure 2. Research area selection and computing improvement of abnormal recognition methodThis 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.
$\Delta T_{ B18V(18H)} = T_{{{B}}F} - \frac{1}{4}\left( {T_{{BU}}{{ + }}T_{{BD}}{{ + }}T_{{BL}}{{ + }}T_{{BR}}} \right){\text{，}}$
(1) 式中：TB为微波亮温值；下标U，D，L，R分别为目标点上、下、左、右4个像元点.
辐射检测S指数的计算公式为
${S\!_f} = \left\{\!\!\!\!\!\! {\begin{array}{*{20}{c}}{\sqrt {\Delta T_{{ B}f{ V}}^2 + \Delta T_{{ B}f{ H}}^2} }\, ,\\\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\! {0, }\end{array}\begin{array}{*{20}{c}}{\quad\quad \Delta T_{{ B}f{ V}} > 0,\;\;\Delta T_{{ B}f{ H}} > 0} \\\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!{\text{其它}}, \end{array}} \right.$
(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.
统计学上，一组数据由小到大重新排列之后，将其均分为4等份，其3个分割点上的数值即为四分位数，依次分别为下四分位数Q1、中位数M和上四分位数Q2，通过这3个四分位数进一步计算得到四分位距IQR
${ IQR} = {Q_2} - {Q_1}$
(3) 和上下边界. 通过计算这些统计量，生成一个箱体（图3），箱体内包含大部分正常数据，而箱体上下边界之外的数据则被认为是异常数据.
在得到上一步计算的S指数的时间序列之后，为了提取S指数的异常变化，利用四分位法构建了辐射检测S指数的背景场，具体过程为以目标地震年以外的其它年份同时间同目标点像元的S指数构成数组，利用四分位数法，分别计算该数组的下四分位数Q1、中位数M、上四分位数Q2和四分位距IQR. 由于所检测的微波异常属于正异常，因此结合实际情况只计算S指数的上限L1，即
${L_1} = M + 1.5{ IQR}.$
(4) 最后根据逐像元计算得到的S指数和构建的背景场对研究区域各像元的S指数时间序列分布进行筛选，最终对微波异常的位置和时间进行判定.
3. 结果分析
3.1 亮温月均值变化及空间分布
理解堪察加半岛微波亮温的月变化规律及其空间分布特征可以为震前异常识别提供基础认知. 选取2003—2011年频率为18.7 GHz的微波亮温数据，对堪察加半岛的月均值进行统计. 图4给出了堪察加半岛平均微波亮温的月变化曲线，可以看出，其大体可分为两个阶段：1—4月和11—12月为平稳期，微波亮温均值稳定地维持在较低水平，即250 K以下；4—10月为剧变期，其中4—7月急剧增温，8—10月急剧降温，微波亮温均值曲线呈抛物线型，随着夏秋季节地表温度的变化，微波亮温均值在7月和8月达到最大，约为263 K.
![]()
图 4 2003—2011年堪察加半岛微波亮温月均值Figure 4. Monthly mean of microwave brightness temperature in Kamchatka Peninsula from 2003 to 2011图5给出了堪察加半岛不同地形微波亮温的差异，从图中可以看出，中部的堪察加河和谷地以及西部的平原地区微波亮温较低，而火山与地震频发的南部和东南部区域微波亮温较高.
![]()
图 5 堪察加半岛2003—2011年亮温月均值影像Figure 5. The image of monthly mean brightness temperature of Kamchatka Peninsula from 2003 to 20113.2 火山点、地震点与非震点亮温对比
由于火山和地震带均位于堪察加半岛的东南区域，因此对地震微波异常进行识别分析时，火山的辐射是不可忽视的. 根据美国国家海洋气象局（National Oceanic and Atmospheric Administration，2011）中的全球火山监测数据对2003—2011年间勘察加半岛的火山进行统计，结果列于表2. 堪察加半岛共有7座火山在上述时间范围内存在活动迹象，火山类型大部分为成层火山，网站记录的活动次数最高达到16次.
表 2 2003—2011年堪察加半岛剧烈活动火山统计Table 2. The statistics of active volcanoes in Kamchatka Peninsula from 2003 to 2011火山名称 东经/° 北纬/° 海拔/m 火山类型 上次显著喷发时间 活动次数 Koryaksky 158.7 53.3 3 430 成层火山 2009-09 4 Karymsky 159.4 54.0 1 513 成层火山 2009-08 13 Bezymianny 160.6 56.0 2 882 成层火山 2009-11 12 Gorely 158.0 52.6 1 799 破火山口 2010-07 1 Kizimen 160.3 55.1 2 334 成层火山 2011-10 3 Klyuchevskoy 160.6 56.1 4 754 成层火山 2010-06 9 Sheveluch 161.4 56.7 3 283 成层火山 2011-04 16 结合收集的AMSR-E微波亮温数据和表2所列7座火山在2003—2011年间的具体活动，筛选出Koryaksky，Karymsky，Bezymianny，Klyu-chevskoy，Sheveluch这5座火山作为研究对象，用V1—V5进行标记. 选取与上述5座火山地理位置最近的地震点作为参考对象，分别对应于2009年12月12日、2008年3月1日、2004年6月10日、2009年4月23日和2010年5月2日这5次地震的震中像元，用S1—S5表示. 非震点依据微波亮温月均值的空间分布差异选取（图5），即中部山脉亮温相差较大，东西两侧的亮温较为接近. 考虑到半岛亮温以中部山脉为界呈两侧对称，选择与地震点同纬度、与中部山脉大致对称的像元点作为非震点，用N1—N5表示. 图6给出了具体火山点、地震点和非震点的选取位置.
图7所示为火山点、地震点及非震点2003—2011年亮温平均值的逐日变化. 最初猜想，火山点和地震点的亮温均应大于非震点. 但研究结果显示，火山点亮温均小于地震点和非震点，地震点与非震点的微波亮温相互交错. 此结果虽然与猜想不同，但其实更符合半岛地表的实际情况. 虽然研究期间火山活动较为频繁，但活跃程度并不高，除相对较大的火山喷发以外，其它火山因常年被积雪覆盖且海拔较高，所以平均温度较低. 通过地震点和非震点的亮温比较可以发现，4—8月，地震点的亮温略大于非震点，与这些月份地震活动更频繁有关.
![]()
图 7 火山点、地震点和非震点2003—2011年微波亮温TB日均值Figure 7. Daily mean of the microwave brightness temperature TB of the volcanoes，seismic and non-seismic points from 2003 to 20113.3 震前微波异常识别与统计
利用2.2节所述的统计方法，对堪察加半岛2003—2011年17次中强地震进行统计分析，针对每次震例，选取特定研究区域和研究时段，计算研究区域内各像元的S指数值并构建非地震年份S指数的背景场，识别微波异常的时空特征. 本文以4次MS≥6.0地震为例说明异常识别与统计分析过程，其中L1和M分别为由非震年份同期同位置计算的S指数构成的数组按四分位法计算得到的上限和中位数.
图8为2003年6月16日MS6.9地震经异常提取后得到的点位时序分布图，可以看出，9号点在5月18日、5月25日、6月3日和6月10日，21号点在5月18日和6月19日（震后3天），辐射检测S指数均显著增长且均大于设定背景场的阈值上限L1，表现为不同程度的异常变化，其中9号点的S指数在6月10日变化最为明显.
![]()
图 8 2003年6月16日MS6.9地震的S异常指数变化在9号点（a）和21号点（b）的统计结果Figure 8. The statistical results of S index variation at point 9 （a） and point 21 （b） for the MS6.9 earthquake on June 16，2003图9为2004年6月10日MS6.9地震经异常提取后得到的点位时序分布图，可见9号点在5月25日和5月27日异常表现相对较弱. 与2003年6月16日MS6.9地震的对比显示，这两次地震的震中位置为同一个像元，即两次地震的研究区域是一致的，故认为2004年6月10日地震的检测结果可能受到2003年6月16日MS6.9地震的影响，9点号在6月8日和21号点在6月5日S指数的突变增长同样可判定为强异常变化. 对比相似的异常检测结果，推测这两次地震的孕震机理是相同的.
![]()
图 9 2004年6月10日MS6.9地震的S异常指数变化在9号点（a）和21号点（b）的统计结果Figure 9. The statistical results of S index variation at point 9 （a） and point 21 （b） for the MS6.9 earthquake on June 10，2004图10为2006年5月22日MS6.2地震经异常提取后得到的点位时序分布图，其中1号点5月1日的S指数异常增高并达到整个时间序列的最大值；5号点在4月24日、5月4日和5月10日表现为弱异常； 11号点和12号点的S指数也在某些时刻出现峰值，但数值较低，判断为可能的弱异常.
![]()
图 10 2006年5月22日MS6.2地震的S异常指数变化在1号点（a）、5号点（b）、11号点（c）和12号点（d）的统计结果Figure 10. The statistical results of S index variation for the MS6.2 earthquake on May 22，2006 at the point 1 （a），point 5 （b），point 11 （c） and point 12 （d）图11为2007年5月30日MS6.4地震经异常提取后得到的点位时序分布图，从S指数时间序列和背景场的统计分析可以看出，3号点在4月29日，9号点在5月21日和5月31日（震后1天）检测到S指数的异常变化.
![]()
图 11 2007年5月30日MS6.4地震在3号点（a）和9号点（b） S异常指数变化的统计结果Figure 11. The statistical results of S index variation at point 3 （a） and point 9 （b） for the MS6.4 earthquake on May 30，2007对所有17次地震中探测到S指数的地震进行统计，结果列于表3. 17次地震中有13次地震在震前1个月内均检测到不同程度的微波异常，4次地震未检测到异常，MS≥6.0地震均检测到异常，MS5.0—6.0地震中有9次检测到S指数异常，4次未发现. 分析4次未检测到异常的地震案例，认为可能的原因为：① 地震震中虽在陆地但却靠近海岸，使得研究区域内的数据完整性不足或有效值不够；② 部分地震如2003年2月6日MS5.0地震的研究区域或研究时间无法满足预先的设定，用较小的研究区域和时间进行替代而无法检测到异常.
表 3 17次地震S指数的检测结果Table 3. S index detection results of the seventeen earthquakes发震时刻 MS 检测异常点编号 异常出现的时间/d 年-月-日 震前 震后 2003-02-06 5.0 无 无 无 2003-06-16 6.9 9，21 4，13，22，29 3 2004-06-10 6.9 9，21 2，5，14，16 无 2005-06-12 5.0 2，3 11，30 无 2006-05-22 6.2 1，5，11，12 3，12，20，21，24 无 2006-08-17 5.7 11，12 10，25 无 2007-05-30 6.4 3，9 1，9 2 2007-11-25 5.2 1，2，4 15，20，23 无 2008-03-01 5.4 无 无 无 2009-04-23 5.2 2，18，19 3，12，14，30 2 2009-06-06 5.2 无 无 无 2009-11-27 5.0 3，18 1 4 2009-12-12 5.1 7，16 16 无 2010-05-02 5.4 无 无 无 2010-10-19 5.0 1，6，22 3，16，23 无 2010-11-27 5.1 5，21 6，7，8，10 无 2010-12-15 5.0 5，17，18 4，14，16，21 无 4. 讨论与结论
在发现卫星遥感可探测到震前微波辐射异常这一现象之后，研究人员多倾向于选择致灾影响较大的单一震例进行分析，而缺少针对典型地震区域的岩石破裂微波辐射异常普遍性检验与统计特征研究. 本文通过对堪察加半岛地震的统计分析，揭示了半岛全年微波亮温具有明显的季节特征及空间差异性的特点，同时对火山点的微波辐射影响进行了评估. 利用改进后的异常识别方法与背景场的构建对2003—2011年共17次中强地震的震前短临异常的位置和时间信息进行统计，发现17次地震中有13次可以识别出震前异常，MS≥6.0地震均可检测到这些异常，且异常大多发生在地震前1个月内. 这一统计分析充分验证了震前微波辐射异常这一现象的普遍性.
本文的研究仍有很多不足，一方面在于新的统计研究方法虽能够有效地消除地形地貌及气候因素的影响，但对土壤湿度和大气水汽等参数的影响无法剔除；另一方面，所提取的震前短临异常的位置和时间信息中，缺少堪察加半岛的地形构造和地震破裂带等相关资料，因此统计结果需要其它数据的佐证. 本文对于地震前卫星遥感微波辐射异常的统计分析还处于基础阶段，对于堪察加半岛地震频发的特点，认为相邻地震的微波异常之间可能存在相关性，对异常信息之间的相关性进行统计分析可能是未来研究的重点，研究结果亦有助于更好地理解微波异常产生的机制.
-
![]()
图 1 2003—2011年堪察加半岛MS≥5.0地震分布
Figure 1. Distribution of MS≥5.0 earthquakes in Kamchatka Peninsula from 2003 to 2011
![]()
图 2 异常识别方法的研究区域选取与计算改进
Figure 2. Research area selection and computing improvement of abnormal recognition method
![]()
图 4 2003—2011年堪察加半岛微波亮温月均值
Figure 4. Monthly mean of microwave brightness temperature in Kamchatka Peninsula from 2003 to 2011
![]()
图 5 堪察加半岛2003—2011年亮温月均值影像
Figure 5. The image of monthly mean brightness temperature of Kamchatka Peninsula from 2003 to 2011
![]()
图 7 火山点、地震点和非震点2003—2011年微波亮温TB日均值
Figure 7. Daily mean of the microwave brightness temperature TB of the volcanoes，seismic and non-seismic points from 2003 to 2011
![]()
图 8 2003年6月16日MS6.9地震的S异常指数变化在9号点（a）和21号点（b）的统计结果
Figure 8. The statistical results of S index variation at point 9 （a） and point 21 （b） for the MS6.9 earthquake on June 16，2003
![]()
图 9 2004年6月10日MS6.9地震的S异常指数变化在9号点（a）和21号点（b）的统计结果
Figure 9. The statistical results of S index variation at point 9 （a） and point 21 （b） for the MS6.9 earthquake on June 10，2004
![]()
图 10 2006年5月22日MS6.2地震的S异常指数变化在1号点（a）、5号点（b）、11号点（c）和12号点（d）的统计结果
Figure 10. The statistical results of S index variation for the MS6.2 earthquake on May 22，2006 at the point 1 （a），point 5 （b），point 11 （c） and point 12 （d）
![]()
图 11 2007年5月30日MS6.4地震在3号点（a）和9号点（b） S异常指数变化的统计结果
Figure 11. The statistical results of S index variation at point 3 （a） and point 9 （b） for the MS6.4 earthquake on May 30，2007
表 1 AMSR-E辐射计的主要性能参数
Table 1 Main parameters and characteristics of AMSR-E
中心频率/GHz 带宽/MHz 动态范围/K 瞬间视场/km 主波束效率 入射角/° 天底偏角/° 扫描间距/km 扫描方式 6.925 350 2.7—340 43×75 95.1% 55.0 47.5 10.1×9.0 以40 r/min的
速度呈圆锥
形扫描10.650 100 2.7—340 29×51 94.8% 55.0 47.5 10.1×9.0 18.70 200 2.7—340 16×27 95.8% 55.0 47.5 10.1×9.0 23.80 400 2.7—340 18×32 94.8% 55.0 47.5 10.1×9.0 36.50 1 000 2.7—340 8.2×14 93.9% 55.0 47.5 10.1×9.0 89.00（A） 3 000 2.7—340 3.7×6.5 94.5% 55.0 47.5 10.1×9.0 89.00（B） 3 000 2.7—340 3.5×5.9 93.7% 54.5 47.0 4.5×4.4 
下载: 导出CSV
表 2 2003—2011年堪察加半岛剧烈活动火山统计
Table 2 The statistics of active volcanoes in Kamchatka Peninsula from 2003 to 2011
火山名称 东经/° 北纬/° 海拔/m 火山类型 上次显著喷发时间 活动次数 Koryaksky 158.7 53.3 3 430 成层火山 2009-09 4 Karymsky 159.4 54.0 1 513 成层火山 2009-08 13 Bezymianny 160.6 56.0 2 882 成层火山 2009-11 12 Gorely 158.0 52.6 1 799 破火山口 2010-07 1 Kizimen 160.3 55.1 2 334 成层火山 2011-10 3 Klyuchevskoy 160.6 56.1 4 754 成层火山 2010-06 9 Sheveluch 161.4 56.7 3 283 成层火山 2011-04 16
下载: 导出CSV
表 3 17次地震S指数的检测结果
Table 3 S index detection results of the seventeen earthquakes
发震时刻 MS 检测异常点编号 异常出现的时间/d 年-月-日 震前 震后 2003-02-06 5.0 无 无 无 2003-06-16 6.9 9，21 4，13，22，29 3 2004-06-10 6.9 9，21 2，5，14，16 无 2005-06-12 5.0 2，3 11，30 无 2006-05-22 6.2 1，5，11，12 3，12，20，21，24 无 2006-08-17 5.7 11，12 10，25 无 2007-05-30 6.4 3，9 1，9 2 2007-11-25 5.2 1，2，4 15，20，23 无 2008-03-01 5.4 无 无 无 2009-04-23 5.2 2，18，19 3，12，14，30 2 2009-06-06 5.2 无 无 无 2009-11-27 5.0 3，18 1 4 2009-12-12 5.1 7，16 16 无 2010-05-02 5.4 无 无 无 2010-10-19 5.0 1，6，22 3，16，23 无 2010-11-27 5.1 5，21 6，7，8，10 无 2010-12-15 5.0 5，17，18 4，14，16，21 无
下载: 导出CSV
-
陈昊, 金亚秋. 2010. 星载微波辐射计对玉树地震岩石破裂辐射异常的初步检测[J]. 遥感技术与应用, 25（6）: 860-866. Chen H, Jin Y Q. 2010. A preliminary detection of anomalous radiation of rock failures related with Yushu earthquake by using satellite-borne microwave radiometers[J]. Remote Sensing Technology and Application, 25（6）: 860-866（in Chinese）.
邓明德, 樊正芳, 崔承禹, 籍全权, 耿乃光. 1995. 无源微波遥感用于地震预报的实验研究[J]. 红外与毫米波学报, 14（6）: 401-406. Deng M D, Fan Z F, Cui C Y, Ji Q Q, Geng N G. 1995. The experimental study for earthquake prediction by passive microwave remote sensing[J]. Journal of Infrared and Millimeter Waves, 14（6）: 401-406（in Chinese）.
房宗绯, 邓明德, 钱家栋, 尹京苑, 耿乃光, 刘晓红, 樊正芳, 荆惠连. 2000. 无源微波遥感用于地震预测及物理机理研究[J]. 地球物理学报, 43（4）: 464-470. Fang Z F, Deng M D, Qian J D, Yin J Y, Geng N G, Liu X H, Fan Z F, Jing H L. 2000. Application of passive microwave remote sensing technology to the seismic prediction and its physical mechanism[J]. Chinese Journal of Geophysics, 43（4）: 464-470（in Chinese）.
耿乃光, 樊正芳, 籍全权, 崔承禹, 邓明德. 1995. 微波遥感技术在岩石力学中的应用[J]. 地震学报, 17（4）: 482-486. Geng N G, Fan Z F, Ji Q Q, Cui C Y, Deng M D. 1995. Application of microwave remote sensing technology in the rock mechanics[J]. Acta Seismologica Sinica, 17（4）: 482-486（in Chinese）.（未找到本条文献英文信息, 请核对）
孟国杰, 申旭辉, Smirnov V, Rogozhin E A, 伍吉仓. 2009. 堪察加地区现今地壳运动与变形特征研究[J]. 地球物理学报, 52（3）: 720-731. Meng G J, Shen X H, Smirnov V, Rogozhin E A, Wu J C. 2009. Research on characteristics of present‐day crustal motion and deformation in Kamchatka area[J]. Chinese Journal of Geophysics, 52（3）: 720-731（in Chinese）.
王恩元, 王云刚, 李忠辉, 刘晓斐, 赵恩来. 2011. 受载煤体变形破裂微波辐射前兆规律的实验研究[J]. 地球物理学报, 54（9）: 2429-2436. Wang E Y, Wang Y G, Li Z H, Liu X F, Zhao E L. 2011. Experimental study on the microwave radiation precursor laws of loading coal in deformation and fracture process[J]. Chinese Journal of Geophysics, 54（9）: 2429-2436（in Chinese）.
王永前, 施建成, 蒋灵梅, 杜今阳, 田帮森. 2011. AMSR-E低频亮温数据空间分辨率提高以及不同波段亮温数据组合应用时分辨率匹配的算法[J]. 中国科学: 地球科学, 41（2）: 253-264. Wang Y Q, Shi J C, Jiang L M, Du J Y, Tian B S. 2011. The development of an algorithm to enhance and match the resolution of satellite measurements from AMSR-E[J]. Science China Earth Sciences, 54（3）: 410-419.
Gorbatov A, Kostoglodov V, Suárez G, Gordeev E. 1997. Seismicity and structure of the Kamchatka Subduction Zone[J]. J Geophys Res, 1021（B8）: 17883-17898.
Kawanishi T, Sezai T, Ito Y, Imaoka K, Takeshima T, Ishido Y, Shibata A, Miura M, Inahata H, Spencer R W. 2003. The Advanced Microwave Scanning Radiometer for the Earth Observing System（AMSR-E）, NASDA's contribution to the EOS for global energy and water cycle studies[J]. IEEE Trans Geosci Remote Sens, 41（2）: 184-194.
Kozhurin A, Acocella V, Kyle P R, Lagmay F M, Melekestsev I V, Ponomareva V, Rust D, Tibaldi A, Tunesi A, Corazzato C, Rovida A, Sakharov A, Tengonciang A, Uy H. 2006. Trenching studies of active faults in Kamchatka, eastern Russia: palaeoseismic, tectonic and hazard implications[J]. Tectonophysics, 417（3/4）: 285-304.
Maeda T, Takano T. 2008. Discrimination of local and faint changes from satellite-borne microwave-radiometer data[J]. IEEE Trans Geosci Remote Sens, 46（9）: 2684-2691.
Maeda T, Takano T. 2009 Detection of Earthquake-related Rock Failures Associated with the 2008 Sichuan Earthquake（WSANE 2009（Workshop for Space, Aeronautical and Navigational Electronics））[J]. Technical Report of Ieice Sane, 109: 229-232.
Maeda T, Takano T. 2010. Detection algorithm of earthquake-related rock failures from satellite-borne microwave radiometer data[J]. IEEE Trans Geosci Remote Sens, 48（4）: 1768-1776.
Maki K I, Takano T, Soma E, Ishii K, Yoshida S, Nakatani M. 2006. An experimental study of microwave emissions from compression failure of rocks[J]. J Seismol Soc Japan, 58(4): 375-384.
National Oceanic and Atmospheric Administration. 2011. Volcano location database search[EB/OL]. [2016-07-12]. http://volcano.si.edu/search_volcano.cfm.
National Snow and Ice Data Center. 2011 . Search data of AMSR-E/Aqua daily L3 surface soil moisture[EB/OL]. [2016-06-10].https://nsidc.org/data/search/#keywords=AMSRE%252FAqua+Daily+L3+Surface+Soil+Moisture/sortKeys=score,,desc/facetFilters=%257B%257D/pageNumber=1/itemsPerPage=25.
Shibata A, Imaoka K, Koike T. 2003. AMSR/AMSR-E level 2 and 3 algorithm developments and data validation plans of NASDA[J]. IEEE Trans Geosci Remote Sens, 41（2）: 195-203.
USGS. 2011. Search earthquake catalog[EB/OL]. [2016-05-16].https://earthquake.usgs.gov/earthquakes/search/.
-
期刊类型引用(3)
1. 丁逸凡 ，吴立新 ，齐源 ，毛文飞 . 2022年6月阿富汗帕克提卡M6.0级地震前微波亮温异常的多态性及其归因辨析. 遥感学报. 2024(10): 2621-2631 . 
百度学术
2. 刘善军，纪美仪，宋丽美，魏恋欢. 基于改进两步差法的玛多M_S7.4地震微波异常研究. 地震学报. 2023(02): 328-340 . 本站查看
3. 吴立新，齐源，毛文飞，刘善军，丁逸凡，荆凤，申旭辉. 多波段多极化被动微波遥感地震应用研究进展与前沿方向探索. 测绘学报. 2022(07): 1356-1371 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 1697
- HTML全文浏览量: 907
- PDF下载量: 41
- 被引次数: 6
