基于地表亮溫的山東省地震活動性分析

馬俊飛李 煒史 雯

1）中國山東266300 膠州市地震局

2）中國山東266034 青島市地震監測中心

基于地表亮溫的山東省地震活動性分析

馬俊飛1）李 煒2）史 雯2）

1）中國山東266300 膠州市地震局

2）中國山東266034 青島市地震監測中心

利用MODIS數據研究山東省地表亮溫時空變化及與地震震例和斷裂帶活動的相關關系，結果表明，山東省2015年存在4個明顯的地表亮溫增溫帶，可能預示著未來地震活動性的增強區域；長島震群與乳山震群在陸地上均有亮溫增溫帶與之對應；風儀店、乳山—海陽斷裂帶在2015年地表亮溫明顯升高，沂沭、聊城—蘭考斷裂帶則無明顯變化?？傊?，120°E以東的膠東半島地區地表亮溫升高明顯，預示著未來膠東半島地區地震活動性可能強于魯中、魯西地區。

地震活動性；地表亮溫；MODIS；遙感

0 引言

臨震前往住出現大面積增溫異常，利用衛星熱紅外技術可以確定溫度異常的范圍和幅度，從而進行地震預測。這一科學思想于20世紀80年代末由前蘇聯科學家提出，ГOPHЫЙ BИ（1990）在分析宇航熱紅外衛片時發現，1984年在蘇聯塔姆得—托克勞與塔拉斯—費爾干斷裂帶交匯部位，多數4級以上地震伴隨熱紅外異常現象。隨后，Tronin等（1996，2004）用10 000景NOAA衛星熱紅外影像對中亞地震活動帶進行研究，得出熱紅外異常與中亞地震活動帶之間具有統計相關性的結論。中國熱紅外遙感地震前兆監測始于20世紀90年代初。強祖基等（1990）用NOAA衛星紅外觀測數據研究1990年4月24日青海共和MS7.0等地震前出現的大規模地面增溫現象，發現地震前7—9天孕震區及周圍可出現幾萬至100多萬km2、累積2℃—10℃的持續增溫異常，往往沿斷裂構造帶向四周成帶狀、片狀展布，在多組斷裂構造帶的交匯部位，異常顯著增強?？琢畈龋?003）進一步對熱紅外增溫的理論基礎及與其他異常的關系進行探討，論述了震前衛星熱紅外增溫的普遍性。根據此溫度異常，目前成功實現了某些地區的地震活動性監測（呂琪琦等，2000；康春麗等，2003）。監測溫度為根據普朗克公式計算出的亮溫，亮溫非地物真實溫度，但能反映地物熱輻射能的相對強度，由于計算簡單，在震前溫度異常研究中廣為使用。本研究以MODIS數據為主數據源，結合地震目錄，研究山東省行政區地表亮溫時空變化，以期預測山東省地震活動性。

1 數據處理

采用紅外遙感及地震目錄數據，其中紅外遙感數據采用美國地球觀測衛星系統（EarthObservation System，EOS） 的中分辨率成像光譜儀（Moderate-resolution Imaging Spectroradiometer，MODIS）數據。美國地球觀測衛星系列發射于1999年，現包括Terra和Aqua太陽同軌道衛星，分別用于陸地和海洋觀測，其中Terra在地方時上午過境，而Aqua在地方時下午過境，二者結合可實現晝夜觀測。MODIS傳感器具有36個波段，覆蓋可見光到遠紅外的光譜范圍（0.4 μm—14 μm），分辨率為250 m（2個波段）、500 m（5個波段）和1 km（29個波段），掃描寬度為2 330 km，1—2天可覆蓋全球。本研究使用下載自美國地質調查局地表數據存儲分發中心網站（https://lpdaac.usgs.gov/）的MOD11A2數據集，該數據集利用MODIS 31、32波段以分裂窗算法獲取，空間分辨率1 km，時間分辨率8天；包含白天和夜間2套地表亮溫數據，并配有相應的數據質量控制、觀測時間、觀測角等配套數據。在晝夜數據的選擇上，考慮到膠東半島三面被海包圍而魯西、魯中地區處于內陸，三者受海陸風影響明顯不同，夜間由海至陸的海風使膠東半島降溫幅度超過魯西、魯中地區，從而導致與地震無關的地表亮溫變化，因此本研究僅使用晝間數據。本研究采用MTR工具，將原始數據投影至UTM，并將晝夜數據分開，而后用最大值算法（Max Value Calculation, MVC），對每4景8日數據進行合成取得月數據。通過MVC算法可以較大程度地減小云、雪等的影響，但仍有某些像元（5%左右）質量達不到要求并被設為空值（NoData）。為避免在后續計算中相互影響導致空值區過大，本研究統一將空值設為0。而后各像元乘以0.02，將數據由象元值轉化為開爾文溫度，減去273最終轉化為攝氏度。經此方法獲取2008年、2009年、2012年、2015年山東省逐月日間地表亮溫數據，并將各年度逐月數據進行平均計算，最終獲取年度地表亮溫數據。

2 地表亮溫變化計算

2.1 山東半島地區地震活動背景年確定

為獲取山東省地表亮溫的時間變化，需要確定一個地震活動相對平靜的年份作為背景值。衡量某一年份地震活動程度有多種指標，本研究采用地震活動度（S）作為指標，該指標采用模糊數學方法，綜合考量地震頻度、平均震級或平均釋放能量、最大震級、地震空間分布集中度及記憶效應等因素，可以比較全面的對震活動性進行定量化（谷繼成等，1987）?？紤]到EOS衛星首發于1999年，2000年才形成全年觀測數據，為實現數據時間匹配，本研究選取2000—2015年山東省陸域地震目錄進行時間掃描，計算逐年地震活動度（S），結果見圖1。

圖1 山東省2000—2015年地震活動度Fig.1 Seismic activity of Shandong Province between 2000 and 2015

由圖1可見，在2000—2015年期，山東省陸域地震活動度在2000年、2006年、2008年、2009年呈現低值，代表較低的地震活動性?？紤]到時間間隔較久會因地表土地利用變化等因素對地表亮溫產生影響，并且2008年、2009年連續兩年地震活動度均處低值，代表地震活動性在一段較長時間內處于穩定低值，而2009年值更低于2008年，因此選定2009年作為基準年，同時為了進行比較，選擇2008年及2012年做為參考年。

2.2 山東省地表亮溫變化

在ArcGIS 10軟件環境中使用區間統計（zonal statistic）方法，對山東省2008年、2009年、2012年、2015年每個季節的代表月份（1月、4月、7月、10月）及年平均地表亮溫進行計算，獲得此4年度代表月份及年平均地表亮溫空間統計均值，見表1。

表1 山東省2008年、2009年、2012年、2015年代表月份及年平均地表亮溫Table 1 Monthly land surface temperature in the symbol month and average annual land surface temperature of Shandong Province in the year of 2008, 2009, 2012 and 2015

地震孕育、發生過程中伴隨的地殼活動，如：斷層面摩擦、地層蠕動、巖石微隙破裂乃至巖層應力聚集，產生放熱效應，熱能對地表加熱從而加強地表長波輻射，并表現為地表亮溫的增加。這種亮溫增加受氣候因子影響較小，而更多反映地殼活動。由表1可見，2009年7月、10月及年平均地表亮溫在選擇的對比年中均較低，與使用地震目錄計算獲得的該年度地震活動性低可以互相印證。而2015年地震活動性的增強也被地表亮溫增加所反映。因此，以2015年度地表亮溫數據與2009年數據進行相減運算，獲取2015年相對于2009年的地表亮溫變化，以此揭示山東地震活動性，見圖2。

圖2 山東省2015年地表亮溫變化及地震分布Fig.2 Land surface bright temperature variation and earthquakes of Shandong Province, 2015

由圖2可見，2015年山東省相對于2009年地表亮溫有4個顯著升高帶（見圖中標注數字的橢圓區域），分別為：①膠州市北部往東經城陽、即墨延伸至海陽市；②蓬萊與龍口交界處；③青州、淄博交界處；④濟寧、兗州區域。其中，以膠州北部至海陽市的增溫條帶規模最大，呈北東東走向，延伸約130 km，最寬處約30 km，升溫幅度普遍在3℃以上，部分地區達4℃—6℃；位于120°E以東的膠東半島地區地表亮溫明顯升高，而120°E以西的魯中、魯西大部分地區地表亮溫變化很小甚至下降。

3 地表亮溫變化與地震活動性

3.1 地表亮溫變化與地震關系

為確定地表亮溫變化與山東省地震活動性關系，將山東省2015年ML1.0以上地震與地表亮溫變化圖疊加（圖2），可以發現，增溫區域3、4在空間上無重合或相近的地震震例，說明該增溫區可能由非地震活動引起，但也可能是本區域地震活動增強的前兆。位于膠東半島北部的蓬萊—龍口增溫帶在陸域上有多處地震與之對應，如：2015年5月21日位于煙臺的2次ML2.0以上地震、6月6日位于煙臺的2.1級地震，均與該增溫帶距離不超過30 km；多處2.0級以下地震發生在該區域周邊；特別是該區附近海域，2015年發生7月22日ML4.5長島地震及多處ML≥2.0地震，均位于該區域周邊50 km范圍內。從時間與空間上判斷，該增溫帶應該與長島震群有關，是長島震群活動在陸地上的反映。膠州北部—海陽增溫帶規模最大、增溫幅度最高，但在時間與空間上只有其北端的乳山震群與之鄰近，其中2015年5月22日乳山ML5.0地震與該增溫帶相距不過26 km。

3.2 地表亮溫變化與斷裂帶關系

為進一步判斷增溫帶與地震活動性關系，將山東省主要斷裂帶與地表亮溫變化圖疊加，見圖3。

圖3 山東省斷裂帶分布及2015年地表亮溫變化Fig.3 Fracture zone distribution and land surface bright temperature variation of Shandong Province, 2015

由圖3 可以發現，膠州北部—海陽增溫帶與乳山—海陽斷裂在走向、分布空間上基本保持一致，與唐家莊—即墨斷裂雖有重合但分布并不一致。通過空間分布關系判斷，該增溫帶應更多與乳山—海陽斷裂活動有關，可以作為乳山震群地震活動的一個反映。同樣，蓬萊—龍口增溫帶與風儀店斷裂的空間分布大致吻合。青州—淄博增溫帶與濟寧—兗州增溫帶均處于多組斷裂帶交匯處，如：青州—淄博增溫帶位于益都斷裂、上五井斷裂、淄河斷裂的交匯處；濟寧—兗州增溫帶位于嘉祥斷裂與鄆城斷裂交匯處。這與強祖基（1990）關于地震前多組斷裂構造帶的交匯部位經常有亮溫顯著增強的論述吻合，此兩處增溫區與地震活動性有關的可能性較大，可能是孕震區。由圖3還可以發現，相對于2009年，沂沭斷裂帶和聊城—蘭考斷裂帶2015年的地表亮溫無明顯升高，可能表明山東省兩條主要斷裂帶在2015年地震活動性無明顯增強趨勢。

4 結論

地表亮溫作為地殼內熱釋放的一個重要指標，與地震活動性具有密切關系，與地震震例、斷裂帶分布在時空上關系復雜。通過對山東省2015年地表亮溫變化與地震活動性的時空分布關系進行分析，發現部分震例，如：膠東半島北部及長島震群震例，與地表亮溫增溫帶在空間上接近，而乳山震群與相應地表亮溫增溫帶在空間分布上并不完全一致，而是通過斷裂帶分布聯系在一起。就震例震級而言，ML2.0以上地震與地表亮溫增溫帶的空間對應關系好于ML2.0以下地震，可能與較小震級的能量釋放不足以大幅增強地表亮溫有關。2015年120°E以東膠東半島地區地表亮溫升高明顯，而魯中、魯西地區則只有兩個較小規模的增溫帶，特別是沂沭斷裂帶與聊城—蘭考斷裂帶均無明顯增溫，可能預示著未來膠東半島地區地震活動性強于魯中、魯西地區。應該承認，由于影響地表亮溫的因子較多，氣候、土地利用等因子變化也可能導致地表亮溫發生變化，如何將非映震因子影響剔除出去，應該是熱紅外遙感地震監測的重要研究方向。

谷繼成，魏富勝.論地震活動性的定量化：地震活動度 [J].中國地震，1987，（3）：12-22.

康春麗，陳正位，陳立澤，等.昆侖山口西8.1級地震的衛星熱紅外前兆特征分析[J].西北地震學報，2003，25（1）：15-18.

孔令昌，王政哲，王志敏，等.衛星熱紅外觀測進行地震預報的可能性分析[J].地震地磁觀測與研究，2003，24（5）：7-12.

呂琪琦，丁鑒海，崔承禹.1998年1月10日張北6.2級地震前可能的衛星熱紅外異?，F象[J].地震學報，2000，22（2）：183-188.

強祖基，徐秀登，賃常恭.衛星熱紅外異?！R震前兆[J].科學通報，1990，35（17）：1 324-1 327.

ГOPHЫЙ BИ.地球熱紅外輻射：地震活動性的標志[J].徐秀登，等譯.地震地質譯叢，1990，12（2）：38-39.

Tronin A A.Satellite thermal survey a new tool for the study of seismo-active regions[J].INT J Remote Sensing, 1996, 17（8）:1 439-1 455.

Tronin A A, Biagi P F, Molchanov O A.Temperature variations related to earthquakes from simultaneous observation at the ground stations and by satellites in Kamchatka area [J].Physics and Chemistry of the Earth, 2004,（29）:501-506.

Analysis of seismic activity in Shandong Province based on land surface bright temperature

Ma Junfei1)，Li Wei2)and Shi Wen2)
1)Earthquake Administration of Jiaozhou County,Shandong Province266300,China
2)Qingdao Earthquake Monitor Center,Shandong Province266034,China

Based on MODIS data, the temporal and spatial variation of Shandong province’s land surface bright temperature are studied to find out the relationship between it and the seismic activity.Results show that 4 obvious bright temperature rising zones where seismic activity may enhance exist in Shandong Province in the year of 2015.Both the Rushan and Changdao earthquake swarm can be reflected by land bright temperature rising zone.The Fengyidian, Rushan and Haiyang fault got obvious increasement in bright temperature, whereas the Yishu, Liaocheng—Lankao fault did not.On the whole, the surface temperature of the east part of Shandong province which located east to the 120° E longitude increased more than that of the central and the west part, which probably foreshow more seismic activities in the Jiaodong Peninsula in the future.

seismic activity，land surface bright temperature，MODIS，remote sensing

10.3969/j.issn.1003-3246.2016.06.013

馬俊飛（1976—），男，理學博士，主要從事遙感與地震監測研究工作

青島市地震局地震監測、防御、應急專項課題“衛星遙感多源數據與多元方法在地震監測預測中的應用”；山東省地震局合同制科研項目“地震前后遙感紅外亮溫變化研究”

本文收到日期：2015-10-22