引言

近年来，地震界及遥感界学者利用不同卫星红外数据，研究了地震前地球表面温度场动态分布情况及其与地震活动的关系(Mogi，1984； Gorny et al，1988; Tramutoli et al，2001; Andrew et al，2002; Ouzounov，Freund，2004; Arun et al，2008)，结果均显示震前存在热异常． 这些利用卫星获取温度的研究实际反映了星下具有一定厚度的大气平均热状况变化，仅局限于震前水平方向范围内的热异常分析，而大气在不同高度上热属性差异明显，对大气在垂直方向上的异常分布是否与构造活动密切相关以及大气热异常是否源自本地区下垫面直接加热则鲜有讨论． 本文重点关注地震热异常垂直结构和热源，将有助于解决地震遥感热捕捉技术面临的如何识别非地震构造运动“热”的技术难题．

另外，地震是地球内部构造应力的一种表现，尽管其主要受内部构造应力影响，但地球不是一个孤立的天体，它的运动必然受到宏观天体运动的影响，天体潮汐引力就是主要的外部因素之一(Heaton，1975； McNutt，Beavan，1981; Kilston，Knopoff，1983)，其在遥感中的角色值得关注．

本文利用天体引潮力附加构造应力变化(changing of the additive tectonics stress caused by celestial tide-generating force，简写为CATSCTF)模型计算2005年11月26日江西九江M_S5.7地震的附加构造应力(Ma et al，2007)，探讨天体引力与本次地震构造环境的关系，并以天体引力变化周期为时间背景指导，利用美国国家环境预报中心(NCEP)提供的温度数据，研究九江地震热异常的时空演变特征，结果清晰地显示了断裂构造在应力作用下的活动(Wu et al，2006)．

1.   九江M_S5.7地震构造环境及区域附加构造应力分析

2005年11月26日江西省九江市(29.7°N，115.7°E)发生强烈地震． 在大地构造部位上，震中位于北部秦岭—大别造山带与南部扬子准地台两大地质单元交汇地带(图 1)，该区域地质运动主要受此二者控制(李传友等，2008). 地质构造方面，如图 1所示，主要有震中以北NE向庐江—广济断裂(F₁，郯庐断裂带南段)与NW向襄樊—广济断裂(F₂)相交汇； 震中以南存在NE向丁家山—郎君山—桂林桥—武宁断裂(F₃)、 九江—靖安断裂(F₄)、 湖口—新干断裂(F₅)等． 尽管地震活动相对微弱，但考虑到地震区域位于地震危险转换带上，仍然认为该区具备发生中强地震的构造背景．

图  1  2005年九江M_S5.7地震构造图F₁: 庐江—广济断裂; F₂: 襄樊—广济断裂; F₃: 丁家山—郎君山—桂林桥—武宁断裂；F₄: 九江—靖安断裂; F₅: 湖口—新干断裂

Figure  1.  Map showing seismogenic structures in Jiujiang M_S5.7 earthquake area F₁: Lujiang-Guangji fault; F₂: Xiangfan-Guangji fault; F₃: Dingjiashan-Langjunshan-Guilinqiao-Wuning fault; F₄: Jiujiang-Jing’an fault; F₅: Hukou-Xin’gan fault

下载: 全尺寸图片 幻灯片

为分析天体引潮力附加构造应力对发震断层的作用类型，需计算震中位置日、 月对地球的天体引潮力沿发震构造的主压、 主张应力轴上产生的附加应力P与T，探索其大小、 方向及发震时间与发震区域的关系(Ma et al，2008)，为温度异常背景的选取提供时间指导．

本文根据美国地质调查局(USGS)提供的震源机制解，采用Ma等(2007)的天体引潮力调制模型算法计算了2005年11月14日—12月4日九江M_S5.7地震过程中的附加构造应力值(图 2)． 此间共经历4个周期(标记为A，B，C和D)，其中A代表 11月1—15日，B代表 11月16—30日，C代表 12月1—15日，D代表 12月16—30日． 从图 2可以看出： ① P和T变化方向与P轴和T轴同向； ② 该次地震的天体引潮力附加构造应力变化作用属于“增压促滑型”； ③ 该地震发生在天体引力附加构造应力变化达到相对比较高值段的尾部，而不是发生在它刚刚到达峰值的时段．

图  2  2005年九江M_S5.7地震天体引潮力附加构造应力P和T的变化

Figure  2.  The change of additive tectonic stresses P and T caused by astro-tidal-triggering for the Jiujiang M_S5.7 earthquake in 2005

下载: 全尺寸图片 幻灯片

依据天体引力附加构造应力模型，附加构造应力增加了断层面的正压力，也增加了断层滑动力，从而诱发地震． 同时，地震发生在P和T到达高值段的尾部，表明附加构造应力的作用是连续的． 瞬间变化并不能立即诱发地震，只有构造应力达到岩层的破裂强度时才可能发生地震． 在其它相似的3个周期A，C和D内，由于没有出现类似B周期中的温度变化过程而无地震发生，说明天体引力不能依靠自身诱发地震，只有在活动断层的地应力达到破裂临界点时才能诱发地震．

2.   九江M_S5.7地震大气增温异常时空演变特征

一些研究人员曾利用美国国家海洋和大气管理局(National Oceanic and Atmospheric Administration，简写为NOAA)气象卫星上携带的高级甚高分辨率辐射仪(Advanced very high resolution radiometer，简写为AVHRR)提供的数据研究温度变化过程(Ouzounov，Freund，2004; Arun et al，2008)，但是红外不能穿透云层，当震源上空存在云系时，就不能有效记录红外热异常.Tramutoli等(2005)虽然通过计算热红外异常估计(Robust estimator of thermal infrared anomalies，简写为RETIRA)指数可消除云的影响，但这种方法需要大量的统计数据，因而通常不能够发现震前微弱的温度变化. 而利用NCEP提供的50年以上全球再分析数据资料可以较好地解决这个问题(Kalnay et al，1996)，能够满足震前热异常时空监测研究．

为了减少地形地貌、 地物类型和气象等非震因素干扰，获取由构造活动造成的温度增加信息，需要利用NCEP数据与背景温度相减的方法，获得发震前后该地区异常增温变化图像． 针对九江M_S5.7地震，本文依据天体引力附加构造应力周期(如图 2中所示，天体引力附加构造应力值经历由小—大—小的过程)，采用2005年11月16日温度值为正常背景值(天体引潮力附加构造应力变化值在最低点)，将11月21—29日的温度与之逐日相减，获得同时次、 同区域范围内近地表及其上空不同气压层面上的动态增温变化序列图像(图 3)，以此作为该次地震及其强余震序列临震增温异常分析的客观依据．

图  3  九江M_S5.7地震前后大气增温异常时空动态演变特征(世界时18:00)(图中F₁—F₅同图 1)

Figure  3.  Spatio-temporal dynamic evolution characteristics of abnormal temperature increase at 18:00 (UTC) before and after Jiujiang M_S5.7 earthquake in 2005 The faults F₁—F₅ have the same meaning as those in Fig. 1

下载: 全尺寸图片 幻灯片

如图 3所示，在九江M_S5.7地震过程中热异常增温主要表现为以下两方面:

1)在时间序列上具有持续性. 以近地表 950 hPa气压层为例，结合主余震时序分析，此次地震异常增温演变过程为起始增温(22日)→加强增温(23—25日)→主震(26日8时)→加剧(26日18时)→余震→增温高峰(27日)→余震→衰减(28日)→平静(29日)，这与马未宇等(2006)仅选取1000 hPa等压高度面上的温度异常变化进行分析所描述的现象基本一致． 该过程与遥感-岩石力学试验中岩石受力破裂的红外热像变化过程也相似(Wu et al，2006)，表明震前热异常可能是强构造运动的一种热表现．

2)在空间分布尺度上呈现出如下规律： ① 温度异常随大气高度变化，越接近地表热异常增温幅度越大，增温区与震中对应关系越明显． 随着大气抬升气体扩散，异常区域面积扩大，幅增减小，但仍与地表分布相似，直到逆温层顶部热异常逐渐消失，说明异常热源来自地面； ② 热异常垂直结构与活动构造分布相吻合． 随着地震的临近，增温异常条带呈SW--NE向，且主要沿F₁，F₃，F₄断裂带(图 1)展布，表明该区活动断裂块在发生挤压活动，且构造运动是温度异常变化的主控原因． 如果该热异常是由气象增温引起，那么气温的影响强度在空间上具有大范围、 逐渐过渡等特点； 在时间上则呈数天或数小时连续变化(屈春燕等，2007)． 然而地震热异常在空间域上受活动构造控制，主要集中在活动断裂带上． 因此，由地震热活动引起的地表温度变化在空间分布上表现出孤立性、 大幅度增温等特点； 在时间上，地震热异常的出现具有突发性，即在一定时段内断续突现． 图 3所示时空动态演变特征正好符合地震热异常的表现特征，因此我们认为该热异常增温现象与地震活动有关．

为了进一步验证上述方法的有效性，对2003年10月16日发生在我国云南大姚的M_S6.1地震进行了研究，发现天体引力附加构造应力值变化以及温度变化过程与九江M_S5.7地震基本一致(图 4，图 5)．

图  4  2003年大姚M_S6.1地震天体引潮力附加构造应力P和T的变化

Figure  4.  The change of additive tectonic stressesP and T caused by astro-tidal-triggering for the Dayao M_S6.1 earthquake in 2003

下载: 全尺寸图片 幻灯片

图  5  2003年大姚M_S6.1地震前大气增温异常时空动态演变特征图(世界时18：00)

Figure  5.  Spatio-temporal dynamic evolution characteristics of abnormal temperature increase at 18:00 (UTC) before Dayao M_S6.1 earthquake in 2003

下载: 全尺寸图片 幻灯片

3.   讨论与结论

天体潮汐引力的变化与构造应力的变化具有一定的关系，其对地震的触发和诱发作用并不是突然的，而是一个连续加速的过程． 本文依据天体引潮力附加构造应力变化值的变化周期，利用NCEP资料获得的温度异常时空动态演变图像清楚地反映了九江M_S5.7地震前的热异常变化过程．

1)震前热异常在时间上具有连续性，即增温起始→增温→加剧→高峰→衰减→平静的演变过程，可能正是体现了活动地震构造的初始微破裂→扩大微破裂→地震构造大破裂(发震)→牵动各向活动构造破裂→构造调整(余震)→趋近稳定发震过程的热像反映(马未宇等，2006)． 当然，更深层次的机制分析还需借助遥感-岩石实验的深入开展进行． 此外，还应对更长时间的震例图像进行研究，尤其需要关注出现温度异常而未发生地震的特殊情况．

2)热异常垂直结构与活动构造具有一致性，其主要集中在活动断裂带上，在空间域上受控于活动构造，并且随大气高度而变化，与近地表分布保持相似性． 当气压高度变化至800 hPa时，温度异常开始消失，这也可能意味着不同强度的地震与其影响造成的最大增温异常高度存在着某种关联．

另外，为了更准确地提取和识别震前热异常信息，在以后震例研究中还应顾及不同地形单元条件下的高程及其下垫面属性差异． NCEP 增温异常实质是构造应力急剧增加的表现，通过增温异常的时空动态演变图像，可较直观地显示地震构造活动过程． 因此，以NCEP增温异常图像为主导，以构造附加应力变化为诱导的地震前兆捕捉技术，可为研究短临地震预测提供一种很有前景的新途径．

作者向对本文提出重要修改意见的审稿专家表示衷心的感谢！