利用地震前地球表面温度异常动态变化预测地震备受关注(Mogi，1984)，震前地球表面温度异常动态变化研究结果表明：震前热异常与地震存在关联(Tronin，2002；Saraf 等，2008；Tramutoli 等，2001；Ouzounov和Freund，2004)。前苏联学者Gorny等人(1988)发现震前出现大面积卫星热红外(10.5—12.5 μm)异常波动，开辟了利用热遥感技术大范围、快速、准实时获取地震信息的新途径。RTS(Robust Satellite Technique)技术被引入该研究领域中(Tramutoli，2007；Lisi 等，2010)。由于该方法采用统计算法，需要多年的统计数据，因此在观测不连续地区的地震中无法开展。Yang和Guo(2010)采用了一种新的相减方法，它利用今天的18:00 (世界时)温度减去昨天的18:00 (世界时)温度，依此类推，获得不同时间序列下一系列的温度变化图像，由于气温变化是非常复杂的，不同背景温度的选择所获得结果差异也很大。Zhang等人(2010)利用卫星遥感热辐射数据通过小波变换方法获得亮温数据，但是选择不同的标准函数傅里叶分裂窗会造成结果的差异。另外，由于热红外不能穿透云层，导致有云区域无法准确获得云下温度，研究中为了获得热异常特征，都尝试剔除有云天数据来解决问题。但云是大气水分相态变化的产物，包含大量的热能交换，去云运算可能造成数据观测人为干扰影响。上述研究集中在分析震前在水平范围内的热异常变化，而大气温度在空间垂直高度上的分布是不均匀的，研究表明将复杂的下垫面环境区分识别将对地震热异常识别具有重要意义(Qin 等，2012, 2013)。由构造运动引起的地震热异常，其热力来源是下垫面构造运动引起的热量交换，热动态演化应呈现出与构造紧密关联的非均匀加热，以及自下而上、从地表开始增温，并随大气运动抬升扩散，在一定高度的高空逐渐消亡的过程。

此外，地震是地球内部构造应力的一种表现，尽管其主要受内部构造应力决定，但地球不是一个孤立的天体，它的运动必然受到宏观天体运动的影响，天体潮汐引力就是主要的外部因素之一(McNutt和Beavan，1981；Heaton，1975)，其在地震热变化中的角色值得关注。

2001年11月14日在昆仑山口西青海和新疆交界处发生了Ms8.1级地震。昆仑山西口地震主破裂带主要位于东昆仑断裂南麓冲洪积台地后缘的地貌陡变带和断层谷地内，对青藏高原成因、演化和现今活动性的认识具有独特的地位(单新建 等，2005；任金卫和王敏，2005)。本文利用天体引潮力变化模型计算此次地震的引潮力，探讨引力与地震的关系，并依据引力变化周期为时间背景指导，利用NCEP温度数据，研究2001年中国昆仑山Ms8.1地震热异常时空演变特征。

引潮力是周期性连续变化的，其触诱发地震多发生在周期振幅相位最大值的临近时刻(Ma等，2012)。引潮力计算过程如下：

根据卡尔文的计算方法，任意天体对地球内部任意一点P产生的引潮力位 ${W_i}(P)$ 为

式中， ${P_n}\left({\cos {Z_{\rm{m}}}} \right)$ 为 $\left({\cos {Z_{\rm{m}}}} \right)$ 的勒让德多项式。 ${Z_{\rm{m}}}$ 为星体的天顶距。 $M$ 是月球、地球的质量。 $k$ 万有引力常数， $r$ 震中与地心距离， ${r_{\rm{m}}}$ 月心与地心间距离。

对于月亮在目前精度条件取 $n = 2$ 和 $n = 3$ 则分别有：

同理对于太阳取 $n = 2$ 有

对于地球整体则有：

式中， ${\delta _{\rm{S}}}$ 和 ${\delta _{\rm{m}}}$ 是日、月赤尾 $\varphi $ 震中纬度。

利用以上公式，计算昆仑山Ms8.1地震前后天体引潮力随时间变化曲线(图1)。横坐标为时间序列，纵坐标为天体引潮力，单位为Gal。由图1可以看出，天体引潮力相位变化具有明显的周期性变化，经历低谷、高峰、低谷的循环，地震发生时引潮力值所处最大振幅相位附近。反映引潮力对地应力达到临界点的活动构造发生地震具有触、诱发的作用，但是地应力的强度变化情况如何？为此我们进一步分析潮汐周期背景下的NCEP大气温度变化。

图 1 昆仑山Ms8.1地震天体引潮力时序变化曲线 Figure 1 Variations of tidal force of celestial body before and after Kunlun Ms8.1 earthquakes

为了减少地形地貌、地物类型和气象等非震因素干扰，获取由构造活动造成的温度增加，需要利用NCEP数据与背景温度相减的方法，获得发震前后该地区异常增温变化图像。同时，为了尽可能降低白天太阳辐射对大气加热的影响，主要采用夜间(世界时18:00)的大气温度数据进行处理。NCEP再分析资料数据是由美国国家环境预测中心和大气研究中心网上准实时发布，包含50年以上数据(Kalnay 等，1996)，该资料集是对来自地面、船舶、无线电探空、探空气球、飞机及卫星等的观测资料进行质量控制和同化处理后得到的，空间分辨率为1°×1°，时间间隔为6 h(秦凯 等，2011)，能够满足震前热异常时空监测研究。针对昆仑山地震，本文依据天体引力附加构造应力周期(图1)，采用2001年11月9日温度值为正常背景值(潮汐值在最低点)，将11月10到17日的温度与之逐日相减，获得同时次、相同区域范围内近地表及其上空不同气压层面上的动态增温变化序列图像(图2)，作为该次地震序列临震增温异常分析的客观依据。

图 2 NCEP增温异常时空动态变化特征图(世界时18:00) Figure 2 Temperature evolution of NCEP before and after Kunlun earthquakes (at 18:00 UTC)

如图2所示，在此次地震过程中热异常增温主要表现为，一方面，在时间序列上具有持续性，以近地表600 Pa气压层为例，2001年11月10日出现异常增温现象；11月11日在地震中心以南，持续出现异常增温，较背景温度高出7 K；11月12日异常增温面积、增幅明显减小，仅较背景温度高出4 K；11月13日在西南方向再次出现较大范围的增温，增温幅度较小，震中较背景只增加3 K，表明孕震体中变形的范围正在变大；14日增温异常区向震中逐渐迁移，增幅达到7 K以上，并覆盖震中；15日主震发生后，地形完整性被破坏，形成诸多的塌陷区域或裂隙，成为热通道，释放大量热，在NCEP图像上反映出增温面积进一步扩大，16日增温幅度达到最高10 K以上，并伴有大量的余震发生。17日增温快速衰减，逐渐恢复正常。此次地震异常增温演变时序过程为：起始增温(10日)、加强增温(11日)、增温衰减(12日)、再增温(13—14日)、主震、加剧(15日)、余震、增温高峰(16日)、快速衰减(17日)的演变过程。该过程与遥感—岩石力学试验中岩石受力破裂的红外热像变化过程相似(吴立新 等，2004)，表明震前热异常可能是强构造运动的一种热表现。

另一方面，在空间分布尺度上呈现出如下规律：(1)温度异常随大气高度变化，越接近地表热异常增温幅度越大，增温区与震中对应关系越明显。随着大气抬升气体扩散，异常区域面积扩大，幅增减小，但仍与地表分布相似，直到逆温层顶部热异常逐渐消失，说明异常热源来自地面；(2)热异常垂直结构与活动构造分布相吻合。表明该区活动断裂块在发生挤压活动，且构造运动是温度异常变化的主控原因。因为如果该热异常是由气象增温引起，那么气温的影响强度在空间上具有大范围性、逐渐过渡性等特点；在时间上则呈现数天或数小时连续变化(屈春燕 等，2006)。然而地震热异常在空间域上受活动构造控制，它主要集中在活动断裂带上。因此，由地震热活动引起的地表温度变化在空间分布上则表现出孤立性、大幅度增温等特点；在时间上，地震热异常的出现具有突发性，即在一定时段内断续突现(刘军 等，2014)。图2所示时空动态演变特征图恰好符合地震热异常表现特点，因此认为该热异常增温现象与地震活动有关。为了检验方法的可信度，增加更长时间周期的数据研究，分析处理了与发震时刻引潮力相位变化相似的震前与震后两个时段的温度图像，结果均未出现与发震时段相似的连续大气增温现象。

另外，本次地震震前11月12日大气温度的衰减与2013年4月20日芦山Ms7.0地震震前4月19日地面红外热异常的衰减表现相似(马未宇 等，2014)，这种临震前的衰减可能表征的岩石应力闭锁期，能否为地震的最终到来提供指示，值得探索。

应用大气分层技术，分析昆仑山Ms8.1地震前后大气温度垂直分布动态演化，结果不仅与岩石在应力加载下发生形变—破裂的过程中热释放过程相吻合，而且符合地面对大气加热导致近地表大气升温、抬升、扩散、消亡的大气热动力学基本特性，反映出下垫面构造运动是本次温度异常变化的主控原因，而非区域外部热量输入造成增温。该方法也展示出利用地震前后大气热异常垂直结构和热扩散方式辨识地面热源异常变化，将有助于解决地震遥感热捕捉技术面临的如何识别非地震构造运动“热”的技术难题。

本文研究仍停留在少数的震例分析，该方法对震前热异常识别是否具有普适性，还需开展更广泛研究，在实践中不断检验、修正。