玉樹地震前NCEP異?，F(xiàn)象*

周 波，馬未宇，陳 聰

(1.南京師范大學(xué)虛擬地理環(huán)境教育部重點實驗室，江蘇南京210046;2.中國地震臺網(wǎng)中心，北京100045;3.浙江大學(xué)計算機學(xué)院，浙江杭州310058)

玉樹地震前NCEP異?，F(xiàn)象*

周 波1，馬未宇2，陳 聰3

(1.南京師范大學(xué)虛擬地理環(huán)境教育部重點實驗室，江蘇南京210046;2.中國地震臺網(wǎng)中心，北京100045;3.浙江大學(xué)計算機學(xué)院，浙江杭州310058)

利用天體引潮力周期變化計算模型和NCEP多源綜合溫度數(shù)據(jù)資料，分析了2010年4月14日青海玉樹MS7.1地震的誘因。結(jié)果表明:在玉樹地震前，天體引潮力經(jīng)歷低谷—高峰—低谷—高峰的連續(xù)變化，NCEP遙感數(shù)據(jù)呈現(xiàn)出起始增溫—加強增溫—高峰增溫—增溫衰減的規(guī)律，由此證明以NCEP溫度為主導(dǎo)，以天體引潮力構(gòu)造附加應(yīng)力變化為參考，是一種有希望的預(yù)測思路。

玉樹地震;天體引潮力;NCEP遙感數(shù)據(jù);異常增溫

0 引言

地震是地球內(nèi)部快速而劇烈的構(gòu)造運動，是一種力學(xué)過程。不同學(xué)者從天體引潮力的誘震作用 (Tramutoli et al，2005)和震前溫度變化角度對地震進行了大量的研究 (Gorny et al，1988;Friedemann，2002;Qiang et al，1999;Tronin et al，2002;Milne，John，1913)，結(jié)果表明:當(dāng)震源系統(tǒng)巖石中的構(gòu)造應(yīng)力達到臨界狀態(tài)時，天體引潮力變化是觸、誘發(fā)地震的重要外部因素之一 (馬未宇 等，2006，2008;Ma et al，2007)。Heaton(1975)從理論上證實:由于地球介質(zhì)的不均勻性，潮汐應(yīng)力在地球內(nèi)部的分布也是不均勻的，地殼中的實際潮汐應(yīng)力值往往大于平均值，而潮汐應(yīng)力的積累速度通常比構(gòu)造應(yīng)力大兩個數(shù)量級，當(dāng)構(gòu)造應(yīng)力處于臨界狀態(tài)時，迭加上這樣快速的脈沖應(yīng)力，從力學(xué)觀點來看完全有可能加速破裂過程而觸發(fā)地震。但天體引潮力具有明顯的周期性，不是每個周期都會發(fā)生地震，如何判斷地應(yīng)力強度是否達到發(fā)震的臨界狀態(tài)，是這一問題的關(guān)鍵 (馬未宇等，2008)。另一方面，早期地面站點記錄向衛(wèi)星遙感技術(shù)觀測發(fā)展，大震前地溫的異常變化研究而更趨廣域、快捷、高效，但存在許多問題 (Kalnay et al，1996)，如:紅外線不能有效穿透云層，無法獲得云下地面溫度;地震引起的增溫過程和氣象增溫疊加，直接、單一利用衛(wèi)星紅外技術(shù)將無法有效獲取異常增溫信息;尤其是采用統(tǒng)計方法獲取異常溫度，長期平均溫度可能掩蓋地震引起的短臨增溫波動，而且由于不同研究者采用的統(tǒng)計年時間域長短不一樣等原因也會造成增溫背景選擇的不確定性，引起異常增溫判識的歧義性。探究震前增溫與天體引潮力變化，其所揭示的物理特性與地震活動表現(xiàn)出的增溫異常現(xiàn)象，在本質(zhì)上是一致的，反映的都是構(gòu)造運動達到一定程度發(fā)生突變—短臨地震活動發(fā)生這一臨界點的判定問題，二者具有明顯的互補性。天體引潮力作為目前唯一能夠預(yù)先計算出的地球形變現(xiàn)象，在時間域上具有一定指示作用，而熱異常表現(xiàn)的構(gòu)造運動為地應(yīng)力強度判識提供支持。因此本文將天體引潮力計算和多源遙感數(shù)據(jù)結(jié)合提取異常增溫信息，對青海玉樹地震進行初步研究。

1 天體引潮力附加構(gòu)造應(yīng)力的作用分析

2010年4月14日青海省玉樹縣 (33.1°N，96.6°E)發(fā)生MS7.1地震，該地震發(fā)生在印度板塊向北東方向運動的前緣和受歐亞板塊、揚子板塊阻擋形成的斷裂破碎帶前沿，為構(gòu)造應(yīng)力集中區(qū)域(錢曉東，秦嘉政，2010)。圖1所示為中國西南地區(qū)的地殼斷裂圖。從圖中可以看到，玉樹MS7.1地震發(fā)生于甘孜—玉樹—風(fēng)火山斷裂上。在一系列近SEE走向的平行斷裂中，甘孜—玉樹—風(fēng)火山斷裂是一條走向南偏東約70°的走滑斷裂，斷層傾角較陡，為晚更新世—全新世、距今約10～20萬年的活動斷裂，該斷裂以左旋水平運動為主 (聞學(xué)澤等，2003)，伴有各處不等但量值較小的垂直運動，且斷裂的南西側(cè)地塊相對抬升 (張希等，2010)。

圖1 玉樹地區(qū)附近地質(zhì)構(gòu)造Fig.1 Geological structure in Yushu and its adjacent area

本文根據(jù)美國USGS提供震源機制解，采用Ma等 (2007)的方法，計算了天體引潮力沿發(fā)震構(gòu)造主壓應(yīng)力軸?P與主張應(yīng)力軸?T的引潮力分量并繪制成圖2，圖中T1表示熱異常起始觀測時間，T2表示玉樹地震發(fā)震時間。由圖2可以看出，玉樹地震過程中，引潮力沿主壓應(yīng)力軸方向分量變化不明顯，而其沿主張應(yīng)力軸分量變化劇烈，且有很強的周期性:每個周期經(jīng)歷低谷—高峰—低谷的過程。地震發(fā)生的周期B中，天體引潮力由高向低連續(xù)變化接近最低值，表明此次地震過程中引潮力對發(fā)震構(gòu)造的作用主要表現(xiàn)為加速斷層面的滑動，與該地區(qū)構(gòu)造滑動為主的左旋走滑運動特征吻合，使斷層失穩(wěn)，釋放已有的構(gòu)造應(yīng)變能而誘發(fā)地震。

圖2 玉樹地震前后天體引潮力變化圖Fig.2 Variation of astro-tidal-triggering before and after Yushu earthquake

2 NCEP增溫異常狀況分析

為了解決紅線外不能穿透云層、以及地震引起的增溫過程和氣象增溫疊加干擾而無法有效獲取異常增溫信息的問題，消減干擾因子的作用，本文采用NCEP全球再分析溫度數(shù)據(jù) (由Global Rawinsonde Data、COADS Surface Marine Data、Aircraft Data、 Surface Land Synoptic Data、 Satellite Sounder Data、SSM/I surface wind speeds 6類數(shù)據(jù)組成)，并對幾乎所有的觀測資料進行同化處理，能夠較綜合、準確反映多因素作用下的實際溫度狀況。

同時為了避免由統(tǒng)計算法可能帶來背景選擇時間域的不確定性問題，本文根據(jù)天體引潮力周期提取時間信號，以2010年4月8日 (天體引潮力主張應(yīng)力分量開始減小)為背景溫度，對2010年4月8～16日共9天的NCEP逐日數(shù)據(jù)減背景溫度，獲得了玉樹地震前后的地區(qū)逐日增溫圖像，如圖3所示。

玉樹附近地區(qū)于2010年4月8日出現(xiàn)異常增溫現(xiàn)象，較背景值高5℃ ～6℃，4月9日震中持續(xù)增溫，甘孜—玉樹斷裂帶方向增溫明顯，增溫6℃～7℃，表明溫度場在空間上正在朝東特別是東南方向，且孕震體中變形的范圍正在變大;4月10日異常增溫面積有所縮小，但異常增溫幅度明顯，且主要向甘孜—玉樹斷裂遷移匯集，在震中以東方向增溫幅度高達9℃，異常在空間上呈現(xiàn)間斷分布;4月11日震中附近保持7℃ ～8℃增溫，而震中往東方向異常增溫帶沿甘孜—玉樹斷裂走向繼續(xù)擴大并沿斷裂連續(xù)貫通;4月12日，震中增溫達到峰值接近10℃，此時由于羌塘地塊向東運動的速率加大而巴顏喀拉地塊速率不變，故2個地塊交匯處的甘孜—玉樹—風(fēng)火山斷裂表現(xiàn)為帶有左旋性質(zhì)的走滑斷裂，且增溫帶面積達到最大，與東南方向的增溫帶合并;4月13日，震中附近出現(xiàn)異常減溫現(xiàn)象，與前一天相比減溫達到5℃～6℃。增溫面積迅速減小;4月14日異常增溫面積減小到震中附近，溫度保持在前一天的異常減溫溫度，表明最終的斷層破裂起始點集中在一個較小的空間范圍內(nèi)，此時天體引潮力張力方向分量達到低谷，對板塊運動起促滑作用，地震爆發(fā);4月15日和16日，異常增溫趨于正?；⒅饾u消失。玉樹地震前異常增溫過程與遙感巖石力學(xué)試驗中巖石受力破裂的紅外熱像變化過程 (Wu，Cui，2000)相似 (加載初期紅外輻射上升—加載另一端遷移擴展—主破裂發(fā)生前增溫衰減—破裂—趨于平靜)，體現(xiàn)了構(gòu)造運動擠壓—巖石微破裂—能量積累—巖石破裂加強—能量釋放—發(fā)震的物理過程。更值得一提的是異常前端，西藏自治區(qū)那曲地區(qū)聶榮縣 (32.5°N，92.8°E)，在北京時間2010年4月17日發(fā)生MS5.2地震，該區(qū)域也是熱異常延伸的最西端，再次表明震前熱異常區(qū)域可能是強構(gòu)造運動的一種熱表現(xiàn)。

天體引潮力變化A、C周期相同但未出現(xiàn)與B周期類似的溫度變化過程，也說明天體引潮力是使壓應(yīng)力處于高臨界狀態(tài)的構(gòu)造帶誘發(fā)地震的外部因素，溫度變化的過程正是天體引潮力使巖石連續(xù)積累地應(yīng)力的過程。而臨失穩(wěn)前出現(xiàn)了短暫的降溫，是巖石破裂失穩(wěn)的重要晚期前兆，其出現(xiàn)意味巖石即將發(fā)生宏觀破裂和失穩(wěn)災(zāi)變，可以做出臨期災(zāi)變預(yù)報。

圖3 玉樹地震前震中附近地區(qū)異常增溫變化特征Fig.3 Variation characteristic of abnormal warming near the epicenter region before Yushu earthquake

3 結(jié)論和探討

(1)通過對反映地應(yīng)力強度的NCEP增溫異常與天體附加構(gòu)造應(yīng)力時序變化關(guān)系的研究，可以得出地應(yīng)力是內(nèi)因根據(jù)，天體引潮力附加構(gòu)造應(yīng)力是觸發(fā)地震外因條件的發(fā)震機理。

(2)震前異常增溫一般呈現(xiàn)起始—加強—高峰—衰減—發(fā)震—平靜的過程，與巖石受力破裂的過程具有相似性。尤其當(dāng)?shù)卣鸩皇前l(fā)生在異常增溫達到高峰時，而是發(fā)生在異常增溫衰退的時段，將對地震的短臨預(yù)測具有積極作用。

(3)從玉樹地震的發(fā)生時段與天體引潮力附加構(gòu)造應(yīng)力時序變化圖的所在位置看，該次地震發(fā)生在低谷時段。表明天體引潮力觸發(fā)、誘發(fā)高應(yīng)力活動斷裂帶的地震機理有減壓促滑現(xiàn)象存在，其本質(zhì)是使處于高應(yīng)力狀態(tài)的活動斷層失穩(wěn)破裂，說明研究引潮力與地震關(guān)系不能脫離構(gòu)造環(huán)境。

(4)地震短臨期間，紅外增溫異常實質(zhì)是地應(yīng)力急劇增加的表現(xiàn)，通過紅外增溫異常可較直觀地顯示地震構(gòu)造活動狀況，圈定增溫異常區(qū);根據(jù)天體引潮力附加構(gòu)造應(yīng)力在不同方向、深度的斷裂帶上的大小和時間演變與發(fā)震時刻的關(guān)系，可以推斷可能發(fā)震的斷裂和發(fā)震時間。因此以NCEP溫度為主導(dǎo)，以構(gòu)造附加應(yīng)力變化為誘導(dǎo)，以地震地質(zhì)為基礎(chǔ)這種面中求線、線中求點的地震捕捉戰(zhàn)術(shù)，是一種很有價值的預(yù)測思路。

馬未宇，徐秀登，徐保華.2006.印度尼西亞9.0地震序列與增溫異常和天體引潮力的相關(guān)關(guān)系研究［J］.西北地震學(xué)報，28(2):129－133.

馬未宇，趙虎，周波.2008.汶川8.0地震天體引潮力和溫度變化過程［J］.國際地震動態(tài)，(11):104.

錢曉東，秦嘉政.2010.2010年青海玉樹地震及震后青藏高原強震趨勢分析［J］.地震研究，33(4):255－264.

聞學(xué)澤，徐錫偉，鄭榮章，等.2003.甘孜—玉樹斷裂的平均滑動速率與近代大地震破裂［J］.中國科學(xué)(D輯:地球科學(xué))，33(S1)，199－208.

張希，張曉亮，張四新.2010.青藏塊體東北緣近期GPS水平運動特征與汶川大震影響［J］.地震研究，33(4):265－268.

Friedemann F.2002.Charge generation and propagation in igneous rocks［J］.Journal of Geodynamics，33(4)，543 －570.

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

Heaton T H.1975.Tidal triggering of earthquakes［J］.Geophys J R astr Soc，43(307):307 ～326.

Kalnay E，Kanamitsu M，Kistler R，et al.1996.The NCEP/NCAR 40 years reanalysis project［J］.Bulletin of the American Meteorological Society.77(3):437－471.

Ma W Y，Zhang X C，Dai X F，et al.2007.A preliminary study on the use of NCEP temperature images and Astro-Tidal-Triggering to forecast short-impending earthquake［J］.Earthquake Research in China，21(1):85－93.

Milne，John.1913.Earthquakes and Other Earth Movements［M］.New York:DAppleton.

Qiang ZJ，Dian CG，Li LZ.1999.Satellite thermal infrared brightness temperature anomaly image-short-term and impending earthquake precursors［J］.Science in China(Science D)，42(3):313 － 324.

Tramutoli V，Cuomob V，F(xiàn)ilizzolab C.2005.Assessing the potential of thermal infrared satellite surveys for monitoring seismically active areas:The case of Kocaeli()earthquake，August 17，1999［J］.Remote Sensing of Environment，96(3)，409 －426.

Tronin，Hayakawa M，Molchanov O A.2002.Thermal IR satellite data application for earthquake research in Japan and China［J］.Geodyn，33(4－5)，519－534.

Wu L，Cui C.2000.Remote sensing rock mechanics(RSRM)and associated experimental studies［J］.International Journal of Rock Mechanics and Mining Sciences，37(6):879 －888.

Abnormal Phenomenon of NCEP Before Yushu Earthquake

ZHOU Bo1，MA Wei-Yu2，CHEN Cong3
(1.Key Laboratory of Virtual Geographic Environment，Ministry of Education，Nanjing Normal University，Nanjing 210046，Jiangsu，China)
(2.China Earthquake Networks Center，Beijing 100045，China)
(3.College of Computer Science and Technology，Zhejiang University，Hangzhou 310058，Zhejiang，China)

Using the cycle process model of the astro-tidal-triggering and according to the NCEP integrated multisource temperature data，we analyzed the triggering factors of Yushu MS7.1 earthquake in Qinghai occurred on Apr.14 in 2010.The result indicated that astro-tidal-triggering went through a continuous change between peak and trough，such as“trough-peak-trough-peak”，and the NCEP remote sensing data showed the law of“original temperature rise-enhancement-reaching peak-attenuation”before Yushu MS7.1 earthquake.Thus，under the dominant of temperature variation in the NCEP remote sensing data and taking the tectonic additional stress variation of astro-tidal-triggering as the reference is a promising idea to predict the earthquake.

Yushu earthquake; astro-tidal-triggering; NCEP remote sensing data; abnormal temperature rise

P315.728

A

1000－0666(2012)03－0330－05

2011－05－13.

中國博士后基金項目 (20090460403)、江蘇省普通高校研究生科研創(chuàng)新計劃資助項目 (CXZZ11_0878)和南京師范大學(xué)研究生科研創(chuàng)新計劃資助項目聯(lián)合資助.