胡亞軒 趙凌強 宋尚武 莊文泉

1 中國地震局第二監(jiān)測中心，西安市西影路316號，710054

中國東北地塊現(xiàn)今地殼運動主要受西太平洋板塊俯沖及日本東北地震的影響。太平洋板塊自北向南以不同的速度俯沖至歐亞板塊/北美板塊和菲律賓板塊之下，其中向歐亞板塊俯沖的速度約為80～90 mm/a，從東至西依次經(jīng)過日本島弧東部海域-日本列島-日本海-長白山-松遼盆地，以一個復(fù)雜的海溝-島弧-弧后系統(tǒng)作用于中國大陸東部[1]。自晚中生代以來，由于太平洋板塊西向俯沖作用，中國東北地區(qū)巖石圈經(jīng)歷多期擠壓和拉張演化過程，期間發(fā)生不同時期的地震及火山活動，引起構(gòu)造變形，形成一系列裂谷帶、盆地、火山帶等[2]。日本列島受太平洋板塊、北美板塊、歐亞大陸板塊以及菲律賓板塊相互作用影響，構(gòu)造運動活躍，地震、火山等自然災(zāi)害頻發(fā)。太平洋板塊向北西方向俯沖至日本島之下，引發(fā)日本海溝多次7級及以上板內(nèi)地震。2011-03-11日本本州東海岸發(fā)生MW9.0地震(3·11大震)，日本近海溝處最大同震滑移量超過50 m[3]。中國大陸東北和華北地區(qū)觀測站觀測到的同震水平位移為cm級，張性應(yīng)變相對較明顯，最大處約為40 nanostrain[4]。

中國東北地區(qū)發(fā)育眾多第四紀陸內(nèi)火山，廣義的長白山火山在中國境內(nèi)包括天池火山、望天鵝火山、圖們江火山和龍崗火山。根據(jù)長白山天池火山噴發(fā)的旋回性和周期性規(guī)律以及長期的火山地震、溫泉和形變監(jiān)測結(jié)果可知[5]，長白山天池火山未來具有再次發(fā)生大規(guī)模爆炸式噴發(fā)的危險性。與長白山天池火山毗鄰的龍崗火山群，具有高密度分布、多中心爆炸式噴發(fā)特征[6]，噴發(fā)活動從早更新世持續(xù)到全新世，具有多旋回、多期次和多階段噴發(fā)特點。其中，火山群西北部的金龍頂子火山為最年輕的近代活動火山，在距今約1 600 a發(fā)生過噴發(fā)，因此也具有潛在的噴發(fā)危險[7]。

本文基于日本GNSS地球觀測網(wǎng)絡(luò)(GEONET)資料和中國東北地區(qū)GNSS及水準觀測資料，綜合分析區(qū)域三維地殼運動特征以及3·11大震前后水平速度場變化，探討太平洋板塊俯沖及殼幔對流對長白山火山區(qū)現(xiàn)今地殼運動及火山活動的影響。

1 日本東北地區(qū)水平運動特征

日本GNSS連續(xù)觀測臺網(wǎng)始建于1993年，至1999年已建成由1 000余個連續(xù)觀測站組成的GEONET網(wǎng)，平均點間距約為20 km(圖1)，密集的觀測臺網(wǎng)可為認識大地震變形的動力學機理和區(qū)域介質(zhì)性質(zhì)提供良好的觀測依據(jù)。選取GEONET網(wǎng)中分布在日本東北地區(qū)、距震中不同距離的部分觀測站進行位移分量時間序列分析，研究數(shù)據(jù)來自內(nèi)華達大學大地測量實驗室。

圖1 主要研究區(qū)域及GNSS觀測站分布

日本東北地區(qū)的奧羽山脈地區(qū)(東北地區(qū)中央山脈)是運動較為強烈的區(qū)域，在太平洋板塊近東西向推擠下，山脈地區(qū)向南北兩個方向擠出，奧羽山脈地殼隆升和火山活動會釋放部分能量。對日本東北地區(qū)GNSS觀測站位移時空變化特征進行分析，總體來看，奧羽山脈為地殼變形的主要分界帶[8]。從2014～2019年水平位移得到的面膨脹率可以看出，受太平洋板塊NWW向俯沖擠壓作用影響，現(xiàn)今主壓應(yīng)變?nèi)砸越鼥|西向壓縮為主(圖2)，范圍為距震中500 km區(qū)域，也是大震94%能量的釋放區(qū)域[9]。2009年至2011年3·11大震前，西部海岸至分界帶區(qū)域GEONET連續(xù)觀測站主要為東向運動。以西部觀測站J252、J232、J194、J192為例(圖3)，大震前主要為東向運動，推測在上新世已停止擴張作用的日本海以沿其東緣西向俯沖和擴展運動為主。日本海具有較高的熱流和海底擴張?zhí)卣?，日本海盆東側(cè)多發(fā)地震，表明日本海在繼續(xù)擴展。中部觀測站水平變形較小，圖3中J931、J917觀測站運動速度變化較小。分界帶至東部海岸區(qū)域觀測站則以西向運動為主，從圖3中J172、MIZU、J550觀測站可以看出，西向位移量均大于20 mm，西向運動也是日本東部地區(qū)的優(yōu)勢運動方向，基本上與太平洋板塊的俯沖方向一致。3·11大震發(fā)生后，J550測站東向同震位移達5.4 m(圖4)。大震后日本北部地區(qū)發(fā)生明顯的東向位移，從距離震中最近的J550測站到最西部的J252測站均可以看出震后的指數(shù)衰減。3·11大震影響范圍廣，震后影響隨時間推移而逐漸衰減，地震造成的粘彈性松弛在未來50～100 a尺度上對中國東北地區(qū)具有拉張效應(yīng)，但相對日本海擴展作用的時間影響范圍仍為很小尺度。

圖2 面膨脹率場

圖3 3·11地震前GNSS連續(xù)觀測站位移E分量時間序列

圖4 3·11地震后GNSS連續(xù)觀測站位移E分量時間序列

2 長白山火山區(qū)三維運動特征

2.1 水平運動特征

日本東北部海溝地震頻發(fā)是中國東北區(qū)域應(yīng)力變化調(diào)整的原因之一。日本海溝發(fā)生地震后，日本島、日本海以及中國東北地區(qū)地殼及地幔會出現(xiàn)相應(yīng)活動。大地測量觀測資料主要是反映地殼運動特征，在一定條件下也可反映地幔頂部的水平流動。日本3·11大震在中國105°E以東區(qū)域產(chǎn)生同震變形，震后變形是震后余滑和粘滯性松弛等作用的綜合反映，地震造成的粘彈性松弛對中國東北地區(qū)會產(chǎn)生拉張效應(yīng)。對東北地區(qū)2011～2019年GNSS觀測資料進行解算，結(jié)果表明，相對于歐亞板塊，水平運動以東南向為主，說明受海溝-島弧-弧后體系拉張作用影響，大陸具有向大洋方向運動的趨勢，速率在10 mm/a以內(nèi)(圖5)，點位運動速率較大。3·11大震引起GNSS連續(xù)觀測站JLCB的同震位移近30 mm(圖6)。從CHUN、SUIY測站斜率變化可知，大震后各測站東向運動速度相對增大。從圖5和圖6可以看出，敦化-密山斷裂以東測站的運動速率較大，這也是拉張量較大的原因。JLCB、JLYJ測站的速度略大于CHUN、SUIY測站，且仍大于各點震前的長期運動速度，表明日本大震對中國東北地區(qū)地殼運動趨勢的影響仍在持續(xù)。

藍色表示2011～2019年同震位移；紅色表示3·11大震同震位移[4]

圖6 GNSS連續(xù)觀測站位移E分量時間序列

2.2 垂直運動特征

中國東北地區(qū)在新生代構(gòu)造運動強烈，伊舒地塹和敦密地塹在新生代沉降達數(shù)千米，同時在部分地區(qū)還發(fā)育褶皺構(gòu)造和大量斷裂。1951～1982年一等水準資料表明，吉林-通化地區(qū)上升速率較大[10]。1970年代至1990年代多期水準資料表明，東北地區(qū)火山區(qū)垂直運動明顯，龍崗火山區(qū)表現(xiàn)為相對上升運動。收集1980年代至2010年代一、二等水準路線觀測資料，選用GNSS連續(xù)站垂直速率作為先驗值，通過與水準數(shù)據(jù)進行聯(lián)合解算，獲取統(tǒng)一框架下的垂直形變場(圖7)。結(jié)果表明，龍崗火山區(qū)垂直速率為-0.39～2.16 mm/a，火山區(qū)東部渾江-撫松一帶隆升速率較大。撫松地區(qū)是渾江斷裂在北部的延伸，也是地震多發(fā)地區(qū)，仙人橋溫泉群是長白山火山區(qū)地熱主要分布區(qū)域。

圖7 中國東北地區(qū)垂直運動速度場及龍崗火山區(qū)MT測點分布

3 長白山火山區(qū)動力學環(huán)境分析

3.1 太平洋板塊俯沖與長白山火山區(qū)現(xiàn)今地殼運動

西太平洋俯沖帶以在日本海溝俯沖為主，其以高速度、小傾角向歐亞大陸強烈俯沖推擠，熱物質(zhì)上涌、地幔對流以及巖石層底部拖曳等多種驅(qū)動力共同決定著整個中國東北地區(qū)的構(gòu)造運動狀態(tài)以及地表運動格局，對斷裂帶運動及鄰區(qū)構(gòu)造運動具有重要影響[11]，板塊的西向俯沖在日本海及中國東北地區(qū)形成NWW向主壓應(yīng)力。地震震源深度從日本島向西逐步加深，并延伸至最西端約1 100 km外的中國東北大陸之下，在吉林琿春一帶發(fā)生中國唯一的中深源地震。中國東北深震帶上的深源地震還時有發(fā)生，說明西太平洋板塊仍在向中國東北大陸俯沖，俯沖方向、速度和角度隨時間而發(fā)生變化?，F(xiàn)今板塊的運動方向決定著區(qū)域應(yīng)力變化及地震時空分布特征。受3·11大震震后應(yīng)力調(diào)整影響，現(xiàn)今中國東北地區(qū)地殼運動以粘滯性松弛效應(yīng)影響為主，YYF和DMF兩條斷裂以拉張為主，走滑分量較小，DMF以東觀測站的水平運動速度明顯大于西側(cè)站點。

日本海-中國東北深震區(qū)前緣地震帶是歐亞板塊與西太平洋板塊之間部分區(qū)段的深部分界，日本海的擴張作用可能在上新世后已停止并轉(zhuǎn)為沿其東緣發(fā)生俯沖作用或日本海與日本本島的聚合作用，日本海淺震的發(fā)生以及現(xiàn)今GNSS觀測結(jié)果表明兩者可能仍在發(fā)生作用。歐亞板塊與太平洋板塊碰撞主要出現(xiàn)在日本東北地區(qū)西側(cè)，太平洋板塊俯沖作用影響范圍主要為從海岸到陸地200 km范圍內(nèi)。板塊俯沖對中國東北地塊的影響可能主要表現(xiàn)為與深俯沖有關(guān)的深部熱作用和深源地震，對地殼運動的直接影響較小。中國東北長白山火山是與太平洋俯沖板塊在地幔轉(zhuǎn)換帶內(nèi)的滯留和深部脫水等過程密切相關(guān)的弧后板內(nèi)火山。2002年吉林汪清發(fā)生7.2級深震，受同一區(qū)域應(yīng)力場影響，距離深震區(qū)約250 km以外的長白山天池火山在2002～2006年發(fā)生巖漿擾動，地震事件增多，在火山口附近出現(xiàn)最大14.3 mm/a的水平運動和46.0 mm/a的垂直運動[12]。3·11大震后，中國東北地震主要發(fā)生在松遼盆地和大興安嶺地區(qū)，大震引起的天池火山區(qū)面膨脹為46 nanostrain，體膨脹為14 nanostrain[13]，但地震未明顯增多，表明淺層巖漿未受到擾動[14]。

3.2 巖石圈物質(zhì)流動與長白山火山區(qū)現(xiàn)今地殼運動

中國東北地區(qū)新構(gòu)造變形及板內(nèi)運動明顯受到地幔上隆/火山活動與裂谷盆地擴張等弧后活動的強烈影響，中國東北地區(qū)地殼與上地幔具有較強的橫向不均勻性，較薄的巖石圈與熱物質(zhì)上涌及殼內(nèi)存在巖漿囊有關(guān)，中國東部新生代玄武巖一般被認為是上地幔部分熔融的產(chǎn)物。由于板塊俯沖導(dǎo)致東亞大陸擴張，引起上地幔上拱及局部物質(zhì)熔融并沿巖石圈斷裂帶上涌和噴溢至地表。長白山火山區(qū)的低速異常反映出區(qū)域性熱構(gòu)造信息，上涌的軟流圈可能已經(jīng)到達(或接近)莫霍面位置[15]，同時火山區(qū)的隆升可能說明局部殼幔物質(zhì)活動對地殼活動的影響。郭良遷[10]研究認為，東北斷塊區(qū)現(xiàn)代地形變形態(tài)與地殼厚度變化形態(tài)及其等值線走向截然不同，這可能表明現(xiàn)代地形變未完全反映地殼整體變化，可能主要反映地殼中部和上部地層的變化。自1972年吉林省地震臺網(wǎng)建立后，長白山龍崗火山區(qū)記錄到200余次地震，其中4.0級以上地震4次，且在1997和2009年多次發(fā)生震群活動，震源深度在1985～2002年約為20 km，2003年以后平均深度小于10 km。大地電磁方法獲取的地下電阻率分布特征表明，地下介質(zhì)物性參數(shù)存在差異，火山區(qū)分布的低阻體可能與深部巖漿系統(tǒng)相對應(yīng)。2020年在火山區(qū)布設(shè)由7條剖面、近100個測點組成的大地電磁觀測網(wǎng)，點位分布見圖7，三維電磁反演得到的電性結(jié)構(gòu)結(jié)果揭示殼內(nèi)存在低阻體(圖8)。25 km以上的低阻結(jié)構(gòu)主要分布在水準資料揭示的隆升區(qū)和現(xiàn)今的小震多發(fā)區(qū)，推測巖漿從深部向淺部運移引起火山區(qū)持續(xù)的相對隆升及小震活動。

圖8 龍崗火山區(qū)電性結(jié)構(gòu)及地震分布

4 結(jié) 語

太平洋板塊向西強烈俯沖，使得日本海溝大震多發(fā)，NWW向主壓應(yīng)力從東至西依次作用于日本列島、日本海及中國東北地區(qū)。奧羽山脈是日本東北地殼運動的分界帶，大致為太平洋板塊俯沖作用的影響范圍。3·11大震前，奧羽山脈以西點位以東向位移為主，中部點位位移量較小，西部點位以西向位移為主。3·11大震后，地殼運動以東向位移為主，但奧羽山脈以西點位東向運動速率明顯大于東側(cè)點位。中國東北地區(qū)現(xiàn)今地殼運動特征與太平洋板塊的俯沖、太平洋板塊和菲律賓板塊與歐亞板塊的碰撞和擠壓、日本海溝大震、地幔對流以及斷層調(diào)整等作用有關(guān)；現(xiàn)今地殼運動活躍，存在深、淺源地震；受日本3·11大震震后粘滯性松弛效應(yīng)影響，研究區(qū)東部仍以拉張運動為主。太平洋板塊的深部俯沖引起區(qū)域應(yīng)力變化，2002年在吉林汪清發(fā)生深震。長白山火山區(qū)由于存在殼幔巖漿，應(yīng)力變化引起天池火山巖漿擾動，導(dǎo)致火山區(qū)產(chǎn)生明顯地殼變形。2011年日本3·11地震引起東北地區(qū)應(yīng)力發(fā)生變化和調(diào)整，但并未引起長白山火山區(qū)淺層巖漿擾動。龍崗火山區(qū)地殼的長期隆升和小震發(fā)生，推測是由火山區(qū)巖漿運移所引起。

致謝：本研究中GEONET數(shù)據(jù)的應(yīng)用及文獻資料的積累得益于與東京大學地震研究所Akito Araya教授的合作及探討，工作期間也得到研究所有關(guān)師生的協(xié)助和管理人員的支持，在此一并表示感謝。