鄭增記范麗紅陜西省地震局，西安 710068長安大學(xué)地質(zhì)工程與測繪學(xué)院，西安 710054

GRACE衛(wèi)星觀測到的尼泊爾MS8.1地震震前重力場變化特征

鄭增記1）范麗紅2）
1）陜西省地震局，西安 710068
2）長安大學(xué)地質(zhì)工程與測繪學(xué)院，西安 710054

強(qiáng)震的孕育和發(fā)生必然引起震源區(qū)和外圍地區(qū)一定范圍內(nèi)地球重力場的變化，而GRACE衛(wèi)星則為探測重力場變化提供了一種新的方式。GRACE衛(wèi)星自2002年3月發(fā)射升空以來，已成功記錄到了2004年蘇門答臘MW9.1地震、2010年智利MW8.8地震以及2011年日本MW9.0地震的同震和震后重力場變化。2015 年4月25日，尼泊爾（28.147°N，84.708°E）發(fā)生了MS8.1大地震。這是繼1934年1月15日尼泊爾8級大地震后，尼泊爾遭受的最強(qiáng)烈地震，目前已造成至少8 219人死亡，17 866人受傷。

孫文科指出，理論上大于M9.0的剪切型或者大于M7.5的張裂型地震所產(chǎn)生的同震重力變化均可被GRACE衛(wèi)星探測到。2015年尼泊爾MS8.1地震震前重力場變化在GRACE衛(wèi)星中是否有所反應(yīng)呢？這里，采用GRACE衛(wèi)星RL05月重力場模型數(shù)據(jù)提取了2015年尼泊爾MS8.1地震的震前重力場變化。

采用美國德克薩斯大學(xué)空間研究中心（UTCSR）最新發(fā)布的RL05月重力場模型數(shù)據(jù)，時間范圍取為2012年1月至2015年3月，該數(shù)據(jù)為處理后的規(guī)格化重力場位系數(shù)，最大階數(shù)為60階。該模型已扣除了極潮、固體潮、海潮的影響以及大氣、海洋的非潮汐部分影響，并且相較于RL04模型，RL05模型在信噪比方面有顯著提高。

為了獲取震前重力場隨時間變化的特征，以平均重力場作為背景場，對GRACE月重力場模型進(jìn)行300 km高斯濾波以及P3M6方法進(jìn)行去相關(guān)處理。由于重力場變化時間序列中包含漂移信號，對震前數(shù)據(jù)利用最小二乘擬合去除其中包含的年、半年周期項信號以及S2潮汐波的影響，獲得了2012—2015年的重力場變化速率圖（圖1）。同時，為了壓制降水等季節(jié)性變化信號對時變重力場的影響，還以2012年1～3月月重力場模型數(shù)據(jù)的平均值作為背景場，計算了2013—2015年1～3月相對于2012年的重力場變化圖像（圖略）。最后得出結(jié)論。

（1）尼泊爾及周邊研究區(qū)域重力場具有明顯的季節(jié)性變化特征，震中以南大部分區(qū)域上半年呈現(xiàn)負(fù)變化，在4～5月份達(dá)到最?。s－4～－6×10-8m/s2），下半年呈現(xiàn)正變化，在9～10月份達(dá)到最大（約+6～+8×10-8m/s2）。

（2）2012—2015年期間，震中以東緬甸和泰國及西北方向吉爾吉斯斯坦和塔吉克斯坦地球重力場持續(xù)下降，變化率約為每年－0.5～－1.0×10-8m/s2，震中西南方向印度中部重力場呈現(xiàn)持續(xù)增加，變化率約為每年+0.3～+0.6×10-8m/s2，反映了印度板塊與歐亞板塊長期碰撞運(yùn)動所造成的重力場變化。

（3）2013年，震中的西北方向吉爾吉斯斯坦和塔吉克斯坦開始呈現(xiàn)重力場負(fù)變化，震中以北青藏高原及西南方向印度中部呈現(xiàn)大范圍的重力場正變化；2014年，吉爾吉斯斯坦和塔吉克斯坦地區(qū)的重力場負(fù)變化持續(xù)加強(qiáng)并且范圍擴(kuò)大，震中以北青藏高原的正變化減弱并且范圍減小，但西南方向印度中部的正變化持續(xù)加強(qiáng)；2015年震中及以西附近區(qū)域重力場出現(xiàn)大范圍的反轉(zhuǎn)變化，并且在反轉(zhuǎn)過程中，發(fā)生了MS8.1地震。

圖1　震前2012—2015年重力場變化速率圖