荊 鳳，申旭輝，張鐵寶，熊 攀，孫 珂

(1.中國地震局地震預(yù)測研究所，北京 100036;2.四川省地震局，成都 610041)

0 引言

20世紀(jì)80年代，前蘇聯(lián)科學(xué)家在分析中亞地區(qū)衛(wèi)星熱紅外遙感圖像時意外發(fā)現(xiàn)了地震前存在熱輻射異?，F(xiàn)象[1]，從此，利用衛(wèi)星熱紅外遙感數(shù)據(jù)開展與地震孕育有關(guān)的輻射異常研究得到了廣泛關(guān)注，并取得了一定成果[2－16]。近年來，陸續(xù)有震例研究表明，一些在震前檢測到的輻射異常與發(fā)震斷層具有一定的關(guān)聯(lián)性[17－18]。眾所周知，地震通常是由斷層活動引發(fā)的，有關(guān)應(yīng)力場與溫度場相關(guān)性的巖石力學(xué)試驗結(jié)果表明，溫度場對變形場具有很好的響應(yīng)，對斷層開展溫度變化監(jiān)測有利于研究斷層活動[19－20];研究熱紅外異常與斷層活動的關(guān)系是解決紅外異常與地震關(guān)系的首要條件[21]。此外，通過對主要活動斷裂帶開展紅外輻射特征分析，研究其紅外輻射變化特征與周邊區(qū)域地震活動的關(guān)聯(lián)性，有助于提高對地震紅外信息識別的準(zhǔn)確率和可靠性，為開展地震監(jiān)測提供技術(shù)支撐。

本文針對2010年4月14日發(fā)生在我國青海省玉樹縣的Ms 7.1級地震，對地震前后發(fā)震斷裂的紅外亮溫圖像進(jìn)行了分析，并與2011年同期該斷裂帶的亮溫變化進(jìn)行比較，希望能為進(jìn)一步利用衛(wèi)星熱紅外遙感數(shù)據(jù)開展地震斷層活動性監(jiān)測提供參考。

1 玉樹地震的構(gòu)造背景

圖1 玉樹地震震中及發(fā)震構(gòu)造分布Fig.1 Locations of epicentre and seismogenic structure of Yushu earthquake

2010年4月14日，青海省玉樹藏族自治州玉樹縣發(fā)生了Ms 7.1級地震，這是繼2008年5月12日汶川地震后我國發(fā)生的又一次破壞性極強(qiáng)的大地震。地震的震中位于 N33.2°，E96.6°(圖 1)，震源深度達(dá)14km。發(fā)震斷裂為左旋走滑的甘孜—玉樹斷裂帶。該斷裂帶與鮮水河斷裂帶在甘孜附近呈左階排列，并與鮮水河斷裂帶共同構(gòu)成巴彥喀拉地塊南邊界。根據(jù)震后野外地質(zhì)調(diào)查結(jié)果[22]，本次地震形成的地表破裂帶總長約51 km，實測最大水平位錯發(fā)生在主破裂北段約1.8 m處，認(rèn)為玉樹地震是巴顏喀拉山地塊和川滇地塊向東擠出而產(chǎn)生的。

2 數(shù)據(jù)源及數(shù)據(jù)預(yù)處理

使用 MODIS 熱紅外波段(10.78 ～11.28 μm)遙感數(shù)據(jù)，圖像空間分辨率為1 km。通過數(shù)據(jù)定標(biāo)及普朗克公式計算，將輻射率轉(zhuǎn)換成亮度溫度數(shù)據(jù);通過閾值限制對圖像進(jìn)行去云處理;對1個月內(nèi)數(shù)據(jù)進(jìn)行均值計算，獲得晴空條件下的月際亮溫圖像。利用該方法獲得了2010年玉樹地震發(fā)震構(gòu)造甘孜—玉樹斷裂帶及其周邊地區(qū)的亮度溫度月值圖像。

3 輻射變化特征分析

分析玉樹地震前后震中區(qū)的亮溫圖像(圖2)發(fā)現(xiàn)，在地震前1個月內(nèi)，甘孜—玉樹斷裂帶北段和南段出現(xiàn)了微弱的亮溫高值區(qū)，斷裂帶的西南方向也出現(xiàn)了零散分布的亮溫高值區(qū)，亮度溫度在285 K左右;地震發(fā)生當(dāng)月，斷裂帶北段和南段的2處亮溫持續(xù)增加，且北段高值的覆蓋面積和溫度高于南段，斷裂帶西南方出現(xiàn)的小面積亮溫高值區(qū)也呈現(xiàn)出北移，向震中靠近，且面積擴(kuò)大、強(qiáng)度增強(qiáng)，亮溫最高值接近290 K;地震后一個月，斷裂帶亮溫條帶南北兩端的高值顯著下降，特別是北段的亮溫高值條帶急劇衰減，衰減程度強(qiáng)于南段，而出現(xiàn)在3月和4月份斷裂帶西南側(cè)的亮溫高值區(qū)域也完全消失。

圖2 2010年3—5月甘孜—玉樹斷裂帶亮溫圖像Fig.2 Image of brightness temperature for Ganzi－Yushu fault from March to May，2010

針對地震發(fā)生當(dāng)月(2010年4月)斷裂帶顯示出的亮溫變化特征，重點分析該月每10 d晴空條件下的亮溫分布情況(圖3)。

圖3 2010年4月甘孜—玉樹斷裂帶亮溫圖像Fig.3 Image of brightness temperature for Ganzi－ Yushu fault in April，2010

從圖3可以看出，在2010年4月1—10日，斷裂帶的亮溫高值出現(xiàn)在玉樹斷裂帶的南北兩端，與月際亮溫圖像所顯示的2處亮溫高值分布相同;4月11—20日，南北兩端的亮溫高值區(qū)范圍擴(kuò)大，且整個區(qū)域的亮溫值都有增高的趨勢;4月21—30日，斷裂帶的條帶亮溫值最高，主要分布在斷裂帶北段以及斷裂帶西南一側(cè)?？梢姡敬蔚卣鸢l(fā)震斷裂帶亮溫的最高值出現(xiàn)在震后7～16 d內(nèi)。據(jù)地震發(fā)生期間水汽和CO變化情況的研究結(jié)果[23]，水汽含量和CO的增多現(xiàn)象都發(fā)生在震后，且震后這一時間段內(nèi)水汽含量和CO含量均表現(xiàn)為上升，認(rèn)為這主要是由于震后大量的地下熱水汽以及CO沿地震產(chǎn)生的地表破裂帶向大氣擴(kuò)散所致。因此認(rèn)為，震后7～16 d出現(xiàn)的亮溫增高可能與震后地下大量氣體釋放有關(guān)。

為了與2010年3—5月份甘孜—玉樹斷裂帶亮溫特征進(jìn)行對比，本研究處理并得到了2011年同時期、同區(qū)域的亮溫圖像(圖4)，發(fā)現(xiàn)在2011年3—4月期間，該地區(qū)亮溫分布格局較為均一，整條斷裂亮溫呈均一的條帶狀分布;到了2011年5月，該斷裂帶的亮溫值又整體升高，西南側(cè)出現(xiàn)了大面積的亮溫高值區(qū)域(圖4(c))。通過檢索地震目錄，發(fā)現(xiàn)2011年6月26日在本區(qū)囊謙縣發(fā)生過Ms 5.3級地震，震中位置為 N32.4°，E96.05°，恰好是甘孜—玉樹斷裂帶西南側(cè)的亮溫高值區(qū)域，因此推測這一亮溫值增高現(xiàn)象可能與一個月后發(fā)生的這次地震活動有關(guān)。

圖4 2011年3—5月甘孜—玉樹斷裂帶亮溫圖像Fig.4 Image of brightness temperature for Ganzi－Yushu fault from March to May，2011

4 討論

通過對2010年4月14日青海省玉樹Ms 7.1級地震前后月際亮度溫度圖像的分析看出，甘孜—玉樹斷裂帶上亮溫變化最顯著的區(qū)域發(fā)生在斷裂帶北段，位于(N33°，E97°)點附近，這一位置與震后野外調(diào)查發(fā)現(xiàn)的玉樹地震產(chǎn)生的地表破裂帶相吻合[22]。而震后4月21—30日的斷裂帶及其周邊的亮度溫度也顯示出這個位置及其西南方向是亮度溫度增高區(qū)。結(jié)合本次地震震后發(fā)現(xiàn)的水汽及CO增強(qiáng)現(xiàn)象，認(rèn)為這一亮度溫度增高現(xiàn)象可能是由于同震破裂引起的摩擦和震后大量地下氣體上涌而引起的。單新建等[21]在用紅外資料分析2001年昆侖山口西8.1級地震時也發(fā)現(xiàn)過類似的現(xiàn)象。

在分析本區(qū)熱紅外亮溫數(shù)據(jù)時還發(fā)現(xiàn)，甘孜—玉樹斷裂帶無論在有震還是無震時都呈現(xiàn)出清晰的線性條帶特征，這是由下墊面的輻射特性所決定的[24－25]，這與紅河斷裂帶的亮溫特征相似;屈春燕等[26]對紅河斷裂帶亮溫條帶的分析結(jié)果顯示，在有震發(fā)生年，紅河斷裂帶的亮溫值要高于無震發(fā)生年，而且斷裂的不同部位也由于活動性的不同顯示出不同的亮溫變化特征。而在本文中所分析的2次地震前1個月，斷裂帶位置均發(fā)生了較為顯著的增溫現(xiàn)象，包括斷裂帶上局部及斷裂帶周邊區(qū)域亮溫的增高，這些變化可能均與斷裂活動有關(guān)。

5 結(jié)論

地震的發(fā)生與活動構(gòu)造關(guān)系密切，利用衛(wèi)星熱紅外遙感數(shù)據(jù)研究主要斷裂帶的紅外輻射特征，能夠達(dá)到監(jiān)測地震構(gòu)造帶是否活動的目的。本文利用MODIS熱紅外遙感數(shù)據(jù)，分析了2010年玉樹地震前后甘孜—玉樹斷裂帶的亮溫變化情況，獲得了以下初步認(rèn)識:

1)由于下墊面的輻射特性，一些活動斷裂帶能夠在亮溫圖像上顯示出清晰的條帶狀影像特征，分析這些條帶的亮溫變化，能夠間接獲得斷裂帶的活動信息。

2)2010年4月14日玉樹Ms 7.1級地震期間和2011年6月26日囊謙縣Ms 5.3級地震前一個月，甘孜—玉樹斷裂帶均出現(xiàn)了明顯的亮溫增高現(xiàn)象，認(rèn)為這與斷裂的活動性增強(qiáng)有關(guān)。

3)2010年玉樹地震前，甘孜—玉樹斷裂帶北段出現(xiàn)的亮溫高值區(qū)與野外考察所發(fā)現(xiàn)的本次地震造成的地表破裂帶位置一致，而亮溫最高值出現(xiàn)的時間也與本次地震水汽含量和CO含量顯著增強(qiáng)的時間吻合。因此認(rèn)為，斷裂帶顯示出的亮溫變化與該區(qū)震前及同震的應(yīng)力集中、摩擦增強(qiáng)以及震后大量地下氣體上涌有關(guān)。

4)利用衛(wèi)星熱紅外遙感數(shù)據(jù)監(jiān)測主要活動斷裂帶亮溫變化，同時有效結(jié)合其他觀測成果的分析，能夠為日常地震監(jiān)測提供輔助決策信息。

志謝:感謝四川省地震局提供研究所用的MODIS數(shù)據(jù)。

[1]Gorny V I，Salman A G，Tronin A A，et al.The earth outgoing IR radiation as an indicator of seismic activity[J].Proc Acad Sci USSR，1988，30(1):67－69.

[2]強(qiáng)祖基，賃常恭.青海共和7級地震衛(wèi)星熱紅外臨震增溫前兆[J].現(xiàn)代地質(zhì)，1992，6(3):297－300.Qiang Z J，Dian C G.Satellite thermal infrared impending temperature increase precursor of Gonghe earthquake with magnitude 7.0 Qinghai province[J].Geoscience，1992，6(3):297－300.

[3]Lu Q Q，Ding J H，Cui C Y.Probable satellite thermal infrared anomaly before the Zhangbei Ms 6.2 earthquake on January 10，1998[J].Acta Seismologica Sinica，2000，13(2):203－209.

[4]康春麗.川滇地區(qū)長波輻射場變化與地震活動關(guān)系研究[J].地震，2008，28(3):43－48.Kang C L.Research on relation between the OLR field change and earthquake activity in Chuan－ Dian region of China[J].Earthquake，2008，28(3):43－48.

[5]Ouzounov D，Bryant N，Logan T，et al.Satellite thermal IR phenomena associated with some of the major earthquakes in 1999－2003[J].Physics and Chemistry of the Earth，2006，31:154－163.

[6]Ouzounov D，F(xiàn)reund F.Midinfrared emission prior to strong earthquakes analyzed by remote sensing data [J].Advances in Space Research，2004，33(3):268－273.

[7]Dey S，Singh R P.Surface latent heat flux as an earthquake precursor[J].Natural Hazards and Earth System Sciences，2003，3(6):749－755.

[8]Cervone G，Kafatos M，Napoletani D，et al.Wavelet maxima curves of surface latent heat flux associated with two recent Greek earthquakes[J].Natural Hazards and Earth System Sciences，2004(4):1－16.

[9]Chen M H，Deng Z H，Yang Z Z，et al.Surface latent heat flux anomalies prior to the Indonesia Mw 9.0 earthquake of 2004[J].Chinese Science Bulletin，2006，51(8):1010－1013.

[10]徐秀登，強(qiáng)祖基，賃常恭.非增溫背景下的熱紅外異常兼機(jī)制討論[J].科學(xué)通報，1991，36(11):841－844.Xu X D，Qiang Z J，Dian C G.Thermal infrared anomalies under the non－temperature－increment background and the mechanism discussion[J].Chinese Science Bulletin，1991，36(11):841－844.

[11]強(qiáng)祖基，孔令昌，郭滿紅，等.衛(wèi)星熱紅外增溫機(jī)制的實驗研究[J].地震學(xué)報，1997，19(2):197－201.Qiang Z J，Kong L C，Guo M H，et al.An Experimental study on temperature increasing mechanism of satellite thermal infrared[J].Acta Seismologica Sinica，1997，19(2):197－201.

[12]Freund F T，Rocks that crackle and sparkle and glow:Strange pre－ earthquake phenomena[J].Journal of Scientific Exploration，2003，17(1):37－71.

[13]強(qiáng)祖基，徐秀登，賃常恭.衛(wèi)星熱紅外異常——臨震前兆[J].科學(xué)通報，1990，35(17)，1324－1327.Qiang Z J，Xu X D，Dian C G.Satellite thermal infrared anomalies:Precursor for impending earthquake[J].Chinese Science Bulletin，1990，35(17)，1324－1327.

[14]劉培洵，劉力強(qiáng)，陳順云，等.地表巖石變形引起熱紅外輻射的實驗研究[J].地震地質(zhì)，2004，26(3):502－511.Liu P X，Liu L Q，Chen S Y，et al.An during deformation experiment on the infrared radiation of surficial rocks[J].Seismology and Geology，2004，26(3):502－511.

[15]吳立新，劉善軍，吳育華，等.遙感－巖石力學(xué)(Ⅰ)——非連續(xù)組合斷層破裂的熱紅外輻射規(guī)律及其構(gòu)造地震前兆意義[J].巖石力學(xué)與工程學(xué)報，2004，23(1):24－30.Wu L X，Liu S J，Wu Y H，et al.Remote sensing－ rock mechanics(Ⅰ):Laws of thermal infrared radiation from fracturing of discontinuous jointed faults and its meanings for tectonic earthquake omens[J].Chinese Journal of Rock Mechanics and Engineering，2004，23(1):24－30.

[16]吳立新，劉善軍，吳育華，等.遙感－巖石力學(xué)(Ⅱ)——斷層雙剪粘滑的熱紅外輻射規(guī)律及其構(gòu)造地震前兆意義[J].巖石力學(xué)與工程學(xué)報，2004，23(2):192－198.Wu L X，Liu S J，Wu Y H，et al.Remote sensing－ rock mechanics(Ⅱ):Laws of thermal infrared radiation from viscosity－sliding of bisheared faults and its meanings for tectonic earthquake omens[J].Chinese Journal of Rock Mechanics and Engineering，2004，23(2):192－198.

[17]陳梅花，鄧志輝，賈慶華.地震前衛(wèi)星紅外異常與發(fā)震斷裂的關(guān)系研究——以2001年昆侖山8.1級地震為例[J].地震地質(zhì)，2003，25(1):100－108.Chen M H，Deng Z H，Jia Q H.The relationship between the satellite infrared anomalies before earthquake and the seismogenic fault[J].Seismology and Geology，2003，25(1):100－108.

[18]郭衛(wèi)英，單新建，馬 瑾.對2001年昆侖山口西8.1級地震前斷層帶紅外增溫異常的討論[J].地震地質(zhì)，2004，26(3):548－556.Guo W Y，Shan X J，Ma J.Discussion on the anomalous increase of ground temperature along the seismogenic fault before the Kunlunshan Ms 8.1 earthquake in 2001[J].Seismology and Geology，2004，26(3):548－556.

[19]馬 瑾，馬少鵬，劉培洵，等.識別斷層活動和失穩(wěn)的熱場標(biāo)志:實驗室的證據(jù)[J].地震地質(zhì)，2008，30(2):363－382.Ma J，Ma S P，Liu P X，et al.Thermal field indicators for identifying active fault and its instability from laboratory experiments[J].Seismology and Geology，2008，30(2):363－382.

[20]陳順云，劉力強(qiáng)，劉培洵，等.應(yīng)力應(yīng)變與溫度響應(yīng)關(guān)系的理論與實驗研究[J].中國科學(xué)(D):地球科學(xué)，2009，39(10):1446－1455.Chen S Y，Liu L Q，Liu P X，et al.Theoretical and experimental study on relationship between strees－strain and temperature variation[J].Science in China(D):Earth Science，2009，39(10):1446－1455.

[21]單新建，屈春燕，馬 瑾.衛(wèi)星熱紅外觀測與發(fā)震斷層不同段落交替活動特征分析[J].大地測量與地球動力學(xué)，2005，25(2):58－62.Shan X J，Qu C Y，Ma J.Satellite thermal infrared observation and analysis on alternate activity of different segments of seismogenic fault[J].Journal of Geodesy and Geodynamics，2005，25(2):58－62.

[22]陳立春，王 虎，冉勇康，等.玉樹Ms 7.1級地震地表破裂與歷史大地震[J].科學(xué)通報，2010，55(13):1200－1205.Chen L C，Wang H，Ran Y K，et al.The Ms 7.1 Yushu earthquake surface ruptures and historical earthquakes[J].Chinese Science Bulletin，2010，55(13):1200－1205.

[23]崔月菊，杜建國，陳 志，等.2010年玉樹Ms 7.1地震前后大氣物理化學(xué)遙感信息[J].地球科學(xué)進(jìn)展，2011，26(7):787－794.Cui Y J，Du J G，Chen Z，et al.Remote sensing signals of atmospheric physics and chemistry related to 2010 Yushu Ms 7.1 earthquake[J].Advances in Earth Science，2011，26(7):787－794.

[24]屈春燕，閆麗莉，溫少妍，等.首都圈地區(qū)衛(wèi)星熱紅外亮溫變化特征研究[J].地球科學(xué)進(jìn)展，2011，26(2):202－211.Qu C Y，Yan L L，Wen S Y，et al.A study of variation characteristics of satellite thermal infrared brightness temperature in capital region[J].Advances in Earth Science，2011，26(2):202－211.

[25]郭衛(wèi)英，馬 瑾，單新建.活動斷裂帶的地形地貌差異與紅外亮溫年變特征的研究[J].地球物理學(xué)進(jìn)展，2008，23(5):1437－1773.Guo W Y，Ma J，Shan X J.Researches on the yearly characteristic of IR brightness temperature and the different hypsography on active faults[J].Progress in Geophysics，2008，23(5):1437－ 1773.

[26]屈春燕，單新建，馬 瑾.紅河斷裂熱紅外高溫條帶的成因及其與地震活動性關(guān)系的探討[J].地震學(xué)報，2006，28(1):91－97.Qu C Y，Shan X J，Ma J.Formation cause of thermal infrared high temperature belt along Honghe fault and its relation to earthquakes[J].Acta Seismologica Sinica，2006，28(1):91－97.