姚華建 薛向輝 王寶善 黃顯良 吳華章
1 安徽蒙城地球物理國(guó)家野外科學(xué)觀測(cè)研究站 蒙城 233500
2 中國(guó)科學(xué)技術(shù)大學(xué) 地球和空間科學(xué)學(xué)院 合肥 230026
3 安徽省地震局 合肥 230071
對(duì)地球內(nèi)部及其外部空間的各種物理場(chǎng)的觀測(cè),是探索地球本體及近地環(huán)境空間的物質(zhì)結(jié)構(gòu)、組成、形成和演化,研究各種自然災(zāi)害發(fā)生和演化的規(guī)律,以及開展自然災(zāi)害預(yù)警和預(yù)測(cè)預(yù)報(bào)工作的重要基礎(chǔ)。目前,世界上很多國(guó)家都建有各種地球物理觀測(cè)站點(diǎn),不同區(qū)域的觀測(cè)站點(diǎn)聯(lián)合在一起,形成龐大的地球內(nèi)部和外部空間觀測(cè)網(wǎng)絡(luò),可以實(shí)時(shí)獲取地球與空間的多物理場(chǎng)信息,包括:地球內(nèi)部的振動(dòng)、形變、電磁場(chǎng)、重力場(chǎng)等信息,以及地球外部空間的風(fēng)場(chǎng)、溫度、密度、等離子體、電磁場(chǎng)等多種參數(shù)。安徽蒙城地球物理國(guó)家野外科學(xué)觀測(cè)研究站(以下簡(jiǎn)稱“蒙城站”)是我國(guó)唯一同時(shí)開展地球內(nèi)部和外部空間物理場(chǎng)觀測(cè)和研究的國(guó)家級(jí)野外站,這在國(guó)際上也是非常少見的。
蒙城站于 2007 年在安徽省蒙城地震臺(tái)(中國(guó)地震局國(guó)家基本臺(tái))的基礎(chǔ)上建成,依托中國(guó)科學(xué)技術(shù)大學(xué)與安徽省地震局共同建設(shè),主管單位為中國(guó)科學(xué)院和中國(guó)地震局;2009 年進(jìn)入科學(xué)技術(shù)部國(guó)家野外科學(xué)觀測(cè)研究站。目前,蒙城站是中國(guó)科學(xué)院日地空間環(huán)境觀測(cè)研究網(wǎng)絡(luò)中的重點(diǎn)站,是國(guó)家重大科學(xué)工程“子午工程”(一期)的重要節(jié)點(diǎn)。蒙城站目前開展測(cè)震、地電、地磁、形變等固體地球物理觀測(cè),同時(shí)開展中高層大氣風(fēng)場(chǎng)、溫度、密度、金屬層、大氣氣輝輻射、電離層電子濃度總含量、太陽(yáng)射電等多種近地空間環(huán)境參數(shù)觀測(cè),具備良好的綜合地球物理科學(xué)觀測(cè)、儀器研發(fā)和研究能力。蒙城站還聯(lián)合安徽省內(nèi)其他臺(tái)站開展地球物理場(chǎng)聯(lián)合探測(cè),同時(shí)針對(duì)郯廬—大別構(gòu)造帶開展流動(dòng)地球物理觀測(cè)與研究,為中國(guó)中東部地區(qū)的地震監(jiān)測(cè)、地震災(zāi)害評(píng)估、近地空間環(huán)境監(jiān)測(cè)、空間災(zāi)害(如磁暴)預(yù)警等提供不可或缺的觀測(cè)資料和研究支撐。
蒙城站位于安徽西北部;安徽地處華北地塊、揚(yáng)子地塊和秦嶺—大別造山帶 3 個(gè)大地構(gòu)造單元的接壤地帶,著名的北北東向郯廬斷裂帶穿過安徽境內(nèi)。郯廬斷裂帶是中國(guó)東部最為重要的深大斷裂帶,總長(zhǎng)度超過 2 400 km,郯廬斷裂帶及其兩側(cè)附近先后發(fā)生過 1668 年的山東郯城 81/2級(jí)大地震等多次強(qiáng)震。而秦嶺—大別造山帶附近地震活動(dòng)也比較頻繁,包括造成較大破壞和人員傷亡的 1917 年安徽霍山 61/4級(jí)地震和 2006 年江西九江里氏 5.7 級(jí)地震。此外,蒙城站所在的安徽省也位于長(zhǎng)江中下游成礦帶區(qū)域,礦產(chǎn)資源形成的結(jié)構(gòu)和構(gòu)造背景一直也是熱點(diǎn)研究問題。
為了更好地研究郯廬—大別構(gòu)造帶及其鄰區(qū)的巖石圈結(jié)構(gòu)和變形特征、強(qiáng)震孕育發(fā)生的結(jié)構(gòu)背景和礦產(chǎn)資源形成的深部構(gòu)造環(huán)境,蒙城站科研人員近年來通過固定和流動(dòng)寬頻帶地震臺(tái)陣觀測(cè)資料,采用背景噪聲成像、接收函數(shù)分析等方法,獲得了郯廬斷裂帶中南段及長(zhǎng)江中下游成礦帶區(qū)域的三維巖石圈橫波速度結(jié)構(gòu)和方位各向異性模型[1-4],以及莫霍(Moho)界面和泊松比結(jié)構(gòu)[5],橫向分辨率可達(dá) 50 km。
獲得的結(jié)構(gòu)模型顯示郯廬斷裂帶中南段及鄰近區(qū)域的上地殼速度結(jié)構(gòu)與區(qū)域地質(zhì)構(gòu)造單元具有很好的對(duì)應(yīng)性:盆地區(qū)域表現(xiàn)為明顯的低速異常,在造山帶、隆起區(qū)域、成礦帶區(qū)域表現(xiàn)為明顯的高速異常。郯廬斷裂帶中南段地殼速度結(jié)構(gòu)、地殼厚度和地震活動(dòng)性存在明顯的分段性:其東、西兩側(cè)的地殼厚度存在明顯差異,斷裂帶東側(cè)的揚(yáng)子塊體的地殼厚度不到 30 km。成像結(jié)果發(fā)現(xiàn),在長(zhǎng)江中下游地區(qū)的下地殼和上地幔頂部存在 V 形高速異常帶(圖 1),長(zhǎng)江中下游成礦帶中下地殼普遍存在高速異常體隆起。這可能是來源于古太平洋板塊俯沖洋殼脫水導(dǎo)致上地幔巖石部分熔融,巖漿底侵至地殼下部,冷卻后形成高速異常體。方位各向異性成像模型顯示郯廬斷裂帶中南段地殼的方位各向異性快波方向基本平行于斷裂帶走向,顯示出斷裂帶兩側(cè)華南和華北塊體的相對(duì)運(yùn)動(dòng)對(duì)斷裂帶區(qū)域地殼介質(zhì)變形的控制作用,但巖石圈與下覆軟流圈的各向異性和變形模式顯著不同。
圖1 背景噪聲成像獲得的長(zhǎng)江中下游成礦帶區(qū)域不同深度橫波速度結(jié)構(gòu)模型(a),以及上地幔頂部V型高速異常體地球動(dòng)力學(xué)背景示意圖(b)[2]Figure 1 Shear wave velocity structure at different depths of Middle-Lower Yangtze River Metallogenic Belt from ambient noise tomography (a), and schematic diagram of geodynamic background for the uppermost mantle V-shaped high velocity anomaly (b)[2]
除了巖石圈速度結(jié)構(gòu)外,地下介質(zhì)的密度結(jié)構(gòu)及其變化也是認(rèn)識(shí)地球內(nèi)部構(gòu)造運(yùn)動(dòng)甚至孕震過程的重要參數(shù),而地球重力場(chǎng)是反映地下介質(zhì)密度分布的核心物理場(chǎng),重力場(chǎng)的測(cè)量及其隨時(shí)間的變化反映了地下介質(zhì)密度的分布及其變化。蒙城站科研人員在安徽省內(nèi)每年開展 2 次流動(dòng)重力測(cè)量工作。目前,安徽省重力測(cè)網(wǎng)由相對(duì)、絕對(duì)重力測(cè)量點(diǎn)及連續(xù)重力觀測(cè)站組成,共有 125 個(gè)測(cè)點(diǎn)。2019—2020 年及 2020—2021 年的年度尺重力場(chǎng)變化不大,約在 ±20 μgal(圖 2);整個(gè)觀測(cè)區(qū)域處于相對(duì)穩(wěn)定的狀態(tài)。蒙城站科研人員已經(jīng)發(fā)展了地震數(shù)據(jù)與重力數(shù)據(jù)的聯(lián)合成像算法,正在開展多種地震波數(shù)據(jù)與重力數(shù)據(jù)的聯(lián)合反演,從而更好地約束郯廬—大別構(gòu)造帶巖石圈物性結(jié)構(gòu)及其動(dòng)力學(xué)過程。
圖2 安徽省流動(dòng)重力觀測(cè)獲得的不同時(shí)間段的重力場(chǎng)變化分布Figure 2 Spatial variation of gravity field in different periods obtained from mobile gravity observation in Anhui Province
寬頻帶地震臺(tái)站的空間分布相對(duì)較為稀疏,所得成像結(jié)果的分辨率還難以刻畫斷裂帶的精細(xì)結(jié)構(gòu)。為了更好地獲得郯廬斷裂帶中南段精細(xì)的速度結(jié)構(gòu)模型和變形特征,從 2015 年開始沿著郯廬斷裂帶中南段布設(shè)密集短周期臺(tái)陣 1 600 多臺(tái)次(圖 3),包括密集線性臺(tái)陣和面狀臺(tái)陣;臺(tái)站間距為數(shù)百米到數(shù)千米,觀測(cè)時(shí)間為 1—2 個(gè)月,獲得了一批重要的地震臺(tái)陣觀測(cè)數(shù)據(jù)。
圖3 郯廬—大別構(gòu)造帶布設(shè)的短周期流動(dòng)臺(tái)陣和大地電磁臺(tái)陣Figure 3 Deployed short-period band mobile seismic arrays and MT array in Tanlu-Dabie tectonic belt
基于密集臺(tái)陣數(shù)據(jù)的背景噪聲成像等方法,蒙城站構(gòu)建了橫向分辨率達(dá)數(shù)千米的郯廬斷裂帶中南段不同區(qū)段的精細(xì)地殼淺部速度結(jié)構(gòu)模型[6-8]?;谀P蛯?duì)比分析和模擬得出 4 點(diǎn)科學(xué)結(jié)論:① 郯廬斷裂帶南段(如廬江段、巢湖段、肥東段)的斷裂帶介質(zhì)表現(xiàn)為高速異常(圖 4);由斷裂帶區(qū)域變質(zhì)雜巖和冷卻的巖漿巖所導(dǎo)致,斷裂帶西側(cè)的合肥盆地顯示為明顯的低速異常。② 從張八嶺北部以北的郯廬斷裂帶區(qū)域由 4 條近于平行的北北東向斷層組成,斷裂帶東西寬度達(dá)到 50—70 km;斷裂帶地殼淺部表現(xiàn)出明顯的高低速相間的特征,主要是由于華北克拉通破壞時(shí)期的拉張作用導(dǎo)致一系列斷陷盆地的形成。③ 郯廬斷裂帶巢湖段上地殼三維各向異性模型揭示區(qū)域內(nèi)方位各向異性分布與主要地質(zhì)單元分布相關(guān),主要受形狀優(yōu)勢(shì)方向(如斷裂帶走向)控制;并基于各向同性和各向異性模型提出了郯廬斷裂帶巢湖段的 4 階段構(gòu)造演化模式。④ 郯廬斷裂帶合肥段的多尺度三維地殼模型[9]顯示巢湖西側(cè)的區(qū)域沉積厚度可達(dá) 6—8 km(圖 5),顯著放大地震動(dòng)和造成更大的地震災(zāi)害;該多尺度模型為郯廬斷裂帶設(shè)定地震的強(qiáng)地面運(yùn)動(dòng)模擬和潛在地震災(zāi)害評(píng)估工作提供了重要的基礎(chǔ)。
圖5 郯廬斷裂帶合肥段多尺度速度結(jié)構(gòu)模型[9]Figure 5 Multi-scale velocity model of Hefei segment of Tanlu fault zone[9]
加強(qiáng)大地電磁、流動(dòng)重力等其他地球物理場(chǎng)的觀測(cè),結(jié)合密集地震臺(tái)陣觀測(cè)資料,采用地球物理聯(lián)合成像等方法,能更好地揭示郯廬—大別構(gòu)造帶深淺部精細(xì)結(jié)構(gòu)和變形特征,研究區(qū)域地震孕育發(fā)生與結(jié)構(gòu)的關(guān)系。
中間層頂區(qū)域大氣主要成分為大氣分子;在這一高度范圍,大氣分子間的運(yùn)動(dòng)以碰撞為主,一般認(rèn)為不受地球磁場(chǎng)及地磁活動(dòng)的影響。太陽(yáng)活動(dòng)會(huì)產(chǎn)生強(qiáng)烈的地磁擾動(dòng),同時(shí)伴隨大量高能粒子進(jìn)入地球大氣,在磁層-電離層-熱層產(chǎn)生一系列災(zāi)害性空間天氣現(xiàn)象。近年研究發(fā)現(xiàn)能量粒子注入可以引起中間層痕量成分(如 O3、NOx、HOx等)變化,以及低熱層高度(95 km 以上)大氣溫度變化。但是,對(duì)于中間層的中性大氣是否受到太陽(yáng)活動(dòng)的影響這一問題尚不清楚。
基于流星雷達(dá)探測(cè)的流星尾跡的衰減時(shí)間(雙極擴(kuò)散系數(shù)),通過排除微流星注入大氣速度對(duì)流星燒蝕峰值高度的影響,發(fā)展了中間層頂大氣密度反演算法。該方法應(yīng)用于全球流星雷達(dá)觀測(cè),獲得了不同緯度中間層頂大氣密度季節(jié)變化[10](圖 6)。在此基礎(chǔ)上,發(fā)現(xiàn)極區(qū)中間層頂大氣密度與對(duì)重現(xiàn)性地磁活動(dòng)的響應(yīng)呈現(xiàn)明顯的負(fù)相關(guān);即地磁活動(dòng)Kp指數(shù)增強(qiáng)時(shí),中間層頂大氣密度降低,反之亦然。這與此前報(bào)道的熱層大氣密度與地磁活動(dòng)存在正相關(guān)性的觀測(cè)結(jié)果剛好相反。另外,還發(fā)現(xiàn)在重現(xiàn)性地磁活動(dòng)期間中間層頂風(fēng)場(chǎng)存在明顯西向增強(qiáng)現(xiàn)象,并伴隨顯著的潮汐幅度(特別是半日潮汐幅度)增加。進(jìn)一步分析 2004—2016 年中間層頂中性大氣密度對(duì)強(qiáng)地磁暴(Dst<-100 nT)的緯度依賴特征,蒙城站發(fā)現(xiàn)在強(qiáng)磁暴期間極區(qū)中間層頂大氣密度存在明顯的下降(可超過 10%),兩極地區(qū)的中間層頂大氣密度存在相似的變化并可以延伸至中緯度地區(qū)[11]。這些發(fā)現(xiàn)揭示了強(qiáng)地磁活動(dòng)期間伴隨的高能粒子或者熱層焦耳加熱等過程對(duì)中高層大氣動(dòng)力學(xué)過程的潛在影響。
圖6 北極-南極流星雷達(dá)反演的中間層頂大氣密度季節(jié)變化特征[11]Figure 6 Seasonal variation of relative atmospheric density within mesopause region using meteor radars from Arctic to Antarctica[11]
電離層是距地面約 60—1 000 km 高度的大氣部分電離區(qū)域,是影響無(wú)線電波傳播最重要的媒介。在高空大氣中,流星注入形成的金屬離子是導(dǎo)致高空電離層突發(fā)E層(Es 層)的主要成分。Es 層是 90—130 km 垂直高度區(qū)域中稠密的電離層不均勻體。電離層不均勻體的出現(xiàn)極大地影響著無(wú)線電通信、定位、導(dǎo)航、航空航天等相應(yīng)領(lǐng)域的遠(yuǎn)程空間通訊的準(zhǔn)確性和穩(wěn)定性。
自 20 世紀(jì) 60 年代地球電離層 Es 層被觀測(cè)發(fā)現(xiàn)以來,其形成的物理機(jī)制一直是一個(gè)長(zhǎng)期未解決的問題,而其中最大的一個(gè)難題就是如何解釋 Es 層顯著的季節(jié)變化。通過 COSMIC 衛(wèi)星的 GPS 掩星事件對(duì)全球電離層Es層的強(qiáng)度進(jìn)行了統(tǒng)計(jì)分析,發(fā)現(xiàn)傳統(tǒng) Es 成因的“垂直風(fēng)剪切理論”并不能很好的解釋 Es 層的全球季節(jié)變化特征[12]。在此基礎(chǔ)上,利用國(guó)內(nèi) 120° 子午線和英國(guó)盧瑟福 · 阿普爾頓實(shí)驗(yàn)室 0° 本初子午線附近的電離層數(shù)字測(cè)高儀觀測(cè)鏈,發(fā)現(xiàn) Es 層存在顯著地從冬季半球向夏季半球的子午輸運(yùn)過程,這一全球金屬離子的跨半球輸運(yùn)過程主要受低熱層的大氣經(jīng)向環(huán)流控制(圖 7)。該研究成果對(duì)進(jìn)一步理解 Es 層的形成機(jī)制提供了重要的觀測(cè)和理論依據(jù)[13]。
圖7 我國(guó)東經(jīng)120°子午線和格林威治本初子午線午線的電離層數(shù)字測(cè)高儀觀測(cè)臺(tái)鏈(a);基于2個(gè)臺(tái)鏈電離層數(shù)字測(cè)高儀觀測(cè)擬合的北半球Es水平輸運(yùn)速度(b)Figure 7 Location of the Ionosonde chains along 120oE meridian and Greenwich prime meridian, respectively (a); Fitted meridional transport velocity of Es in northern hemisphere using two meridional chains (b)
在電離層中存在復(fù)雜的不同時(shí)間尺度和空間尺度的變化和結(jié)構(gòu)。一般來說,1 000 km 及以上的大尺度電離層結(jié)構(gòu)和變化主要引起無(wú)線電信號(hào)延遲的差異,而中小尺度的電離層結(jié)構(gòu)和變化還會(huì)造成電離層閃爍,嚴(yán)重時(shí)可導(dǎo)致信號(hào)的中斷。由于觀測(cè)的限制,目前對(duì)中小尺度結(jié)構(gòu)和變化的認(rèn)識(shí)較為缺乏。
基于 GNSS 的電離層觀測(cè)是目前應(yīng)用最廣泛的電離層觀測(cè)技術(shù)。而我國(guó)北斗衛(wèi)星導(dǎo)航系統(tǒng)還包含 5 顆地球靜止軌道同步衛(wèi)星,與美國(guó) GPS 等系統(tǒng)相比擁有獨(dú)特的優(yōu)勢(shì)。北斗同步衛(wèi)星電離層觀測(cè)相當(dāng)于是定點(diǎn)的全天候觀測(cè),能把電離層時(shí)間和空間變化效應(yīng)區(qū)分。借助北斗的觀測(cè)優(yōu)勢(shì),蒙城站創(chuàng)新地建立了中國(guó)中部區(qū)域蒙城北斗電離層短基線觀測(cè)網(wǎng)(圖 8)。該觀測(cè)網(wǎng)絡(luò)臺(tái)站間距基本在 100 km 附近,最小間距在 60 km,可以構(gòu)成兩組三角觀測(cè)網(wǎng)絡(luò)。觀測(cè)網(wǎng)除了可以接收北斗衛(wèi)星信號(hào),還可接收 GPS 和 GLONASS 等導(dǎo)航衛(wèi)星信號(hào)。利用該特色觀測(cè)網(wǎng)絡(luò),結(jié)合蒙城站布設(shè)的光學(xué)氣輝觀測(cè),能有效觀測(cè)和監(jiān)測(cè)區(qū)域中小尺度電離層變化和結(jié)構(gòu)及其引起的閃爍,對(duì)研究和認(rèn)識(shí)電離層區(qū)域結(jié)構(gòu)和變化具有重要作用。
圖8 蒙城北斗電離層短基線觀測(cè)網(wǎng)分布圖Figure 8 Location of Mengcheng ionospheric observation of Beidou Short Baseline Network
利用蒙城北斗電離層短基線觀測(cè)網(wǎng),中國(guó)科學(xué)技術(shù)大學(xué)團(tuán)隊(duì)聚焦電離層中小尺度結(jié)構(gòu)和變化開展了系列工作,取得了創(chuàng)新性成果,包括:首次發(fā)現(xiàn)區(qū)域性 100 km 尺度電離層經(jīng)度梯度(圖 9);黃山站觀測(cè)到不同天北斗同步衛(wèi)星電離層總電子含量(TEC)變化對(duì)比,在 2016 年 3 月 8 日和 11 日不同同步衛(wèi)星 TEC 變化存在明顯不同,3 月 8 日不同同步衛(wèi)星 TEC 較為一致,而 3 月 11 日存在顯著差異,小區(qū)域變化達(dá) 100% 以上。該結(jié)果準(zhǔn)確刻畫既揭示了電離層中尺度行進(jìn)式擾動(dòng)特征及其與不規(guī)則體關(guān)系,發(fā)現(xiàn)低層大氣潮汐活動(dòng)導(dǎo)致電離層變化超過磁暴效應(yīng)的現(xiàn)象等。自 2015 年該網(wǎng)絡(luò)建設(shè)起,利用其觀測(cè)數(shù)據(jù)已發(fā)表多篇重要研究論文[15-17],得到了國(guó)內(nèi)外學(xué)界廣泛關(guān)注和認(rèn)可。
圖9 黃山站觀測(cè)不同天北斗同步衛(wèi)星電離層總點(diǎn)子含量(TEC)[14]Figure 9 Variations of vertical TECs observed by Beidou GEO at Huangshan[14]
通過不斷努力,蒙城站構(gòu)建了從固體地球到高層大氣較為全面的地球物理觀測(cè)網(wǎng)絡(luò),開展了一系列富有成效的地球物理野外觀測(cè)工作,取得了豐碩成果;其在郯廬斷裂帶多尺度地下結(jié)構(gòu)成像與孕震構(gòu)造研究、中緯度區(qū)域中高層大氣和電離層精細(xì)結(jié)構(gòu)和變化機(jī)制等領(lǐng)域的科研成果引起了國(guó)內(nèi)外的關(guān)注,為我國(guó)地球物理野外站從固體地球到行星際空間的綜合觀測(cè)和研究體系的建設(shè)起到了示范作用。
未來蒙城站的觀測(cè)和研究主要從 4 個(gè)方面展開:① 提升地球物理儀器研發(fā)和組網(wǎng)觀測(cè)能力,包括分布式光纖傳感設(shè)備[18]、絕對(duì)重力儀、實(shí)時(shí)傳輸智能地震儀、陣列式大口徑激光雷達(dá)、多站流星雷達(dá)組網(wǎng)觀測(cè)等。② 繼續(xù)開展郯廬—大別構(gòu)造帶多尺度地震、電磁、重力、形變等綜合地球物理觀測(cè),以及多種不同數(shù)據(jù)的聯(lián)合反演及聯(lián)合解釋,深入研究該構(gòu)造帶的深淺部結(jié)構(gòu)和變形、地震活動(dòng)性分布和地震災(zāi)害特征。③ 圍繞近地空間環(huán)境多要素協(xié)調(diào)觀測(cè),加強(qiáng)近地空間環(huán)境主被動(dòng)光學(xué)和無(wú)線電遙感能力建設(shè),為理解近地空間多圈層、多尺度耦合提供觀測(cè)條件。④ 加強(qiáng)固體地球與近地空間的聯(lián)合觀測(cè)與融合研究,更好地約束整個(gè)地球與近地空間系統(tǒng)的物理過程。例如:太陽(yáng)活動(dòng)如何影響地球電離層和地球磁場(chǎng)的變化,大地震活動(dòng)與地球電離層異常之間的關(guān)系,以及電離層變化對(duì)雷達(dá)干涉(InSAR)地表形變測(cè)量的影響等。