林建民， 方孫珂， 倪四道

1 浙江省海洋觀測(cè)-成像試驗(yàn)區(qū)重點(diǎn)實(shí)驗(yàn)室， 浙江大學(xué)， 浙江舟山 316021 2 浙江大學(xué)海洋學(xué)院， 浙江舟山 316021 3 大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 中科院精密測(cè)量科學(xué)與技術(shù)創(chuàng)新研究院， 武漢 430071

0 引言

每逢夏秋，我國(guó)東南沿海地區(qū)經(jīng)常遭受西北太平洋臺(tái)風(fēng)災(zāi)害的侵?jǐn)_，蒙受巨大的經(jīng)濟(jì)損失和人民生命安全威脅.近年來(lái)，隨著海洋和大氣領(lǐng)域科學(xué)的發(fā)展和技術(shù)的革新，臺(tái)風(fēng)研究與數(shù)值預(yù)報(bào)均取得了長(zhǎng)足的進(jìn)步，但仍無(wú)法滿足防災(zāi)減災(zāi)的切實(shí)需求(陳聯(lián)壽，2006；端義宏等，2020).臺(tái)風(fēng)過(guò)境時(shí)現(xiàn)場(chǎng)觀測(cè)資料的缺乏是目前臺(tái)風(fēng)研究的主要短板，其主要原因在于傳統(tǒng)海洋觀測(cè)儀器布設(shè)困難，很難做到大范圍布陣或提前布設(shè)于臺(tái)風(fēng)行經(jīng)海域；而且，臺(tái)風(fēng)這種極端天氣情況下，海洋儀器的正常工作、觀測(cè)數(shù)據(jù)的獲取與質(zhì)量保障等都會(huì)受到嚴(yán)重干擾；此外，由于操作技術(shù)難度高、風(fēng)險(xiǎn)大、成本高等因素，目前國(guó)內(nèi)如氣象飛機(jī)觀測(cè)等先進(jìn)臺(tái)風(fēng)觀測(cè)技術(shù)應(yīng)用相對(duì)較少，也是導(dǎo)致直接觀測(cè)數(shù)據(jù)相對(duì)缺乏的重要原因.因此，目前臺(tái)風(fēng)監(jiān)測(cè)主要采用衛(wèi)星遙感的方式.以合成孔徑雷達(dá)為代表的衛(wèi)星遙感技術(shù)可以獲得準(zhǔn)確的中低風(fēng)速下海面風(fēng)場(chǎng)，但臺(tái)風(fēng)高風(fēng)速區(qū)域往往伴隨著降雨，會(huì)對(duì)其回波信號(hào)造成一定的影響，導(dǎo)致風(fēng)速反演誤差最高可達(dá)10 m·s-1(Uhlhorn and Black，2003；Klotz and Uhlhorn，2014；Zhang et al.，2014).綜合多學(xué)科研究方法與觀測(cè)數(shù)據(jù)，發(fā)展新的臺(tái)風(fēng)監(jiān)測(cè)方法，對(duì)傳統(tǒng)方法進(jìn)行必要的補(bǔ)充和相互印證，將有助于我們提高臺(tái)風(fēng)監(jiān)測(cè)能力.

近年來(lái)，基于臺(tái)風(fēng)激發(fā)的地震背景噪聲、極低頻海洋噪聲進(jìn)行臺(tái)風(fēng)定位追蹤與強(qiáng)度反演的研究逐漸興起(Gerstoft et al.，2006；Chi et al.，2010；夏英杰等，2011；Sun et al.，2013；Davy et al.，2014；Sufri et al.，2014；Wang and Li，2015；李風(fēng)華等，2016；Lin et al.，2018；Fang et al.，2020).海洋已被證明是地震背景噪聲的最主要源區(qū)，廣袤海洋上的波浪每時(shí)每刻作用于固體地球，從而源源不斷地激發(fā)低頻背景噪聲(Nishida，2013，2017).臺(tái)風(fēng)作為北半球逆時(shí)針旋轉(zhuǎn)的熱帶氣旋，是天然的大尺度“攪拌機(jī)”，能夠顯著加劇海-氣相互作用、海浪-海底相互作用及能量耦合，從而使得期間的地震背景噪聲顯著增強(qiáng)并攜帶著臺(tái)風(fēng)過(guò)程物理信息被遠(yuǎn)處地震臺(tái)所接收, 這為地震學(xué)家利用地震背景噪聲研究臺(tái)風(fēng)提供了可能.

臺(tái)風(fēng)引起的地震背景噪聲根據(jù)激發(fā)機(jī)理及優(yōu)勢(shì)頻段，一般分為“地球嗡鳴”(Hum,約0.003～0.02 Hz；Rhie and Romanowicz, 2004，2006；Maurya et al.，2019)和“地脈動(dòng)”(Microseism,約0.05～0.5 Hz；Longuet-Higgins，1950；Hasselmann，1963；Gerstoft et al.，2006；Bromirski，2009；Ardhuin et al.，2015).地球嗡鳴產(chǎn)生自次重力波與沿岸海底的能量耦合，其強(qiáng)度與海上風(fēng)暴的影響范圍密切相關(guān)(Maurya et al.，2019).近期有研究發(fā)現(xiàn)，地球嗡鳴頻段之上還存在著噪聲級(jí)相對(duì)較低的“風(fēng)暴地震”(Stormquake,約0.02～0.05 Hz)，起源于風(fēng)暴激發(fā)長(zhǎng)周期海浪與陸架特征結(jié)構(gòu)的能量耦合，能量強(qiáng)度可相當(dāng)于3.5級(jí)地震(Fan et al.，2019).地脈動(dòng)則是地震背景噪聲中能量最強(qiáng)的分量，又分為單頻和倍頻地脈動(dòng).單頻地脈動(dòng)由海浪產(chǎn)生的壓力擾動(dòng)直接作用于海底或海岸產(chǎn)生，其頻率與波浪頻率一致.由于海浪產(chǎn)生的垂向壓力擾動(dòng)隨水深呈指數(shù)衰減，所以，一般認(rèn)為單頻地脈動(dòng)源區(qū)主要分布在淺海或近岸海區(qū)(Hasselmann，1963).而倍頻地脈動(dòng)則是地脈動(dòng)信號(hào)中的主導(dǎo)者，由相向傳播的兩列同頻率海浪非線性相干形成駐波后產(chǎn)生的垂向壓力擾動(dòng)作用于海底而激發(fā)，其源區(qū)可分布于沿岸、臺(tái)風(fēng)中心附近或兩個(gè)臺(tái)風(fēng)之間(Cessaro，1994；Bromirski et al.，2005；Ardhuin et al.，2011；Ardhuin and Herbers，2013；Ardhuin et al.，2015；Lin et al.，2017；鄭露露等，2017；陳棟爐等，2018).一般認(rèn)為，倍頻地脈動(dòng)主要以Rayleigh波和P波形式傳播.早在20世紀(jì)中葉，遠(yuǎn)海臺(tái)風(fēng)激發(fā)的P波倍頻地脈動(dòng)就已被學(xué)者觀測(cè)發(fā)現(xiàn)(Ramirez，1940；Gilmore，1946).近年來(lái)，得益于現(xiàn)代化的地震臺(tái)陣部署和技術(shù)發(fā)展，Love波及S波形式的倍頻地脈動(dòng)也已被發(fā)現(xiàn)并確認(rèn)(Gerstoft et al.，2008； Tanimoto et al.，2015；Liu et al.，2016).

近些年來(lái)，P波倍頻地脈動(dòng)被地震學(xué)家初步成功應(yīng)用于近實(shí)時(shí)監(jiān)測(cè)定位臺(tái)風(fēng).自Gerstoft等(2006)利用南加利福尼亞地震臺(tái)網(wǎng)(Southern California Seismic Network，SCSN)成功監(jiān)測(cè)到2000多公里外墨西哥灣2005年超強(qiáng)颶風(fēng)“卡特麗娜”(Katrina)后，利用地震臺(tái)陣定位追蹤臺(tái)風(fēng)的研究逐漸興起.Zhang等(2010)分別基于SCSN以及日本高靈敏度地震臺(tái)網(wǎng)(High Sensitivity Seismograph Network，Hi-net)記錄的倍頻地脈動(dòng)信號(hào)，對(duì)2006年超強(qiáng)臺(tái)風(fēng)“伊歐凱”(Ioke)進(jìn)行定位追蹤，得到與美國(guó)國(guó)家颶風(fēng)中心(National Hurricane Center)公布的臺(tái)風(fēng)最佳路徑吻合良好的定位結(jié)果.Lin等(2018)利用布置于中國(guó)東北的NECESSArray地震臺(tái)陣(NorthEast China Extended SeiSmic Array)監(jiān)測(cè)定位2009年超強(qiáng)臺(tái)風(fēng)“盧碧”(Lupit)，結(jié)果與日本區(qū)域?qū)I(yè)氣象中心(Regional Specialized Meteorological Center，RSMC)公布的臺(tái)風(fēng)軌跡及MODIS衛(wèi)星遙感結(jié)果相印證一致.以上基于地震學(xué)的臺(tái)風(fēng)定位研究未來(lái)有望為傳統(tǒng)臺(tái)風(fēng)監(jiān)測(cè)與研究提供跨學(xué)科的數(shù)據(jù)支持、協(xié)助大氣和海洋相關(guān)數(shù)值模型的校正(Zhang et al.，2010；Ardhuin et al.，2012；Lin et al.，2018；Ward Neale et al.，2018).然而，這些應(yīng)用前景須以對(duì)P波倍頻地脈動(dòng)源區(qū)精確定位為基礎(chǔ)，如何切實(shí)提高臺(tái)風(fēng)激發(fā)P波地脈動(dòng)源區(qū)的定位精度是亟待解決的關(guān)鍵任務(wù).

1 數(shù)據(jù)

本文選擇2014年登陸我國(guó)廣東省南部沿岸的第15號(hào)臺(tái)風(fēng)“海鷗”作為研究對(duì)象，并以RSMC公布的臺(tái)風(fēng)最佳路徑數(shù)據(jù)為參考軌跡(數(shù)據(jù)時(shí)間間隔為6 h，http:∥www.jma.go.jp/jma/jma-eng/jma-center/rsmc-hp-pub-eg/besttrack.html).“海鷗”起源于關(guān)島東部海域的熱帶低壓，于世界協(xié)調(diào)時(shí)(Universal Time Coordinated，UTC)2014年9月12日成長(zhǎng)為熱帶風(fēng)暴，并在次日增強(qiáng)為1級(jí)臺(tái)風(fēng)；14至15日期間，“海鷗”橫穿呂宋島進(jìn)入我國(guó)南海，并在登陸海南島前成長(zhǎng)至巔峰(最大持續(xù)風(fēng)速達(dá)38.6 m·s-1)，最終在廣東湛江沿岸第三次登陸后快速減弱，直至消散于我國(guó)內(nèi)陸(圖1).“海鷗”在其生命周期內(nèi)，徑直穿越于開闊的菲律賓海和相對(duì)較為封閉的南海，并具有三次登陸過(guò)程，為激發(fā)地脈動(dòng)提供了較為豐富的場(chǎng)景，是理想的研究對(duì)象.

圖1 地震臺(tái)陣ChinArray II和Hi-net的地理分布及臺(tái)風(fēng)“海鷗”(Kalmaegi)和“風(fēng)神”(Fengshen)最佳路徑軌跡示意圖紅色實(shí)心小圓點(diǎn)表示地震臺(tái)站；不同顏色的實(shí)心圓表示不同時(shí)刻臺(tái)風(fēng)中心位置，時(shí)間間隔為6 h，近臺(tái)風(fēng)中心最大持續(xù)風(fēng)速(Max sustained wind speed)與圓圈內(nèi)不同顏色相對(duì)應(yīng)，并與圓圈半徑成正比.背景地形圖下載自https:∥visibleearth.nasa.gov/images/147190/explorer-base-map/147192l.Fig.1 Scheme of the deployment of ChinArray II and Hi-net, and the best-tracks of Typhoon Kalmaegi and Fengshen The seismic stations are represented by solid red dots, and the typhoon track is indicated by color-coded circles, which are spaced at a 6 h time interval, with the circle sizes proportional to the max sustained wind speed. The background raster map is from https:∥visibleearth.nasa.gov/images/147190/explorer-base-map/147192l.

本文收集了臺(tái)風(fēng)“海鷗”期間(2014年09月13日至16日)地震臺(tái)陣ChinArray II和Hi-net的垂向連續(xù)波形數(shù)據(jù)，分別下載自中國(guó)地震科學(xué)探測(cè)臺(tái)陣數(shù)據(jù)中心(http:∥www.chinarraydmc.cn)和日本Hi-net臺(tái)網(wǎng)網(wǎng)站(http:∥www.hinet.bosai.go.jp).除去存在數(shù)據(jù)質(zhì)量問(wèn)題(數(shù)據(jù)異常、缺失等)的臺(tái)站，兩個(gè)臺(tái)陣分別具有437和747個(gè)有效臺(tái)站用于本文研究.文中使用的總涌浪平均周期數(shù)據(jù)，屬于歐洲中期氣象預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasts，ECMWF)第五代再分析數(shù)據(jù)集(ECMWF Reanalysis v5，ERA5)，網(wǎng)格分辨率為0.5°×0.5°，下載自哥白尼氣候變化服務(wù)(Copernicus Climate Change Service，C3S，https:∥cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form).文中P波地脈動(dòng)源區(qū)海表等效壓力譜數(shù)據(jù)，網(wǎng)格分辨率為0.5°×0.5°，下載自法國(guó)海洋開發(fā)研究院數(shù)據(jù)庫(kù)(the French Research Institute for the Exploitation of the Sea，IFREMER，ftp:∥ftp.ifremer.fr/ifremer/ww3/HINDCAST/SISMO/).

2 方法

2.1 雙臺(tái)陣聯(lián)合FB定位方法

不同于傳統(tǒng)的頻率-波數(shù)(Frequency-Wave number，f-k)方法，本文采用的FB定位方法基于地震臺(tái)陣中不同臺(tái)站記錄的P波地脈動(dòng)信號(hào)的相干性并結(jié)合相應(yīng)走時(shí)預(yù)測(cè)，反演其源區(qū).首先，地震臺(tái)站記錄數(shù)據(jù)均按以下步驟進(jìn)行預(yù)處理：1)去除儀器響應(yīng)；2)去均值、去趨勢(shì)；3)降采樣至1 Hz；4)0.05～0.45 Hz帶通濾波；5)連續(xù)記錄數(shù)據(jù)分割成時(shí)長(zhǎng)為6 h的小段，計(jì)算每段數(shù)據(jù)的標(biāo)準(zhǔn)差，并截去振幅大于1倍標(biāo)準(zhǔn)差的部分，以去除由地震或異常儀器響應(yīng)等造成的瞬時(shí)振幅異常；6)將各段記錄數(shù)據(jù)進(jìn)行窗口長(zhǎng)度為10 min、50%重疊的短時(shí)傅里葉變換(Short-Time Fourier Transform，STFT)獲得連續(xù)頻譜S(ω)，再除以振幅值以只保留相位信息，避免振幅放大效應(yīng)并消除局地噪聲的影響.

根據(jù)預(yù)處理數(shù)據(jù)，構(gòu)建互譜密度矩陣(Cross-Spectral Density Matrix，CDSM)C，它由N個(gè)臺(tái)站的頻譜S(ω)計(jì)算得到，即為N×N的矩陣，其第i行和第j列的元素包含臺(tái)站i和臺(tái)站j之間的相位延遲，表示為

(1)

其中ω為角頻率，t為數(shù)據(jù)段(t-3h,t+3h)的中值時(shí)刻，?為Hermitian轉(zhuǎn)置運(yùn)算符，〈〉表示對(duì)連續(xù)10-min窗口頻譜的集合平均.

基于網(wǎng)格搜索法，本文將目標(biāo)海域(100°E—140°E, 5°N—25°N)劃分為0.2°×0.2°的網(wǎng)格.當(dāng)頻率為ω時(shí)，網(wǎng)格點(diǎn)(lat,lon)處t時(shí)刻的波束Bt(ω,lat,lon)可表示為

Bt(ω,lon,lat)=P?(ω,lon,lat)Ct(ω)P(ω,lon,lat),

(2)

其中指向各網(wǎng)格點(diǎn)的臺(tái)陣導(dǎo)向矢量P為

(3)

其中Tlon,lat指從目標(biāo)源區(qū)網(wǎng)格點(diǎn)(lon,lat)到地震臺(tái)陣的P波走時(shí)向量，即Tlon,lat=(tp1,…,tpN)?.本文采用LLNL-G3Dv3三維地球速度模型(Myers et al.，2015)預(yù)測(cè)各臺(tái)站對(duì)應(yīng)的P波走時(shí)tpi(i=1,2,3…N)，相對(duì)于PREM等一維地球介質(zhì)模型在一定程度上可以減少淺層介質(zhì)復(fù)雜地質(zhì)結(jié)構(gòu)對(duì)P波走時(shí)預(yù)測(cè)的干擾而提高走時(shí)預(yù)測(cè)精度(Fang et al.，2020).

鑒于單臺(tái)陣FB定位結(jié)果通常受空間分辨能力的限制，本文進(jìn)而提出雙臺(tái)陣聯(lián)合FB定位方法，以提升定位精度.根據(jù)以往工作研究(Fang et al.，2020)，當(dāng)臺(tái)風(fēng)激發(fā)的P波地脈動(dòng)源區(qū)至地震臺(tái)陣距離小于～26°時(shí)，F(xiàn)B定位結(jié)果將受到上地幔三重震相(Triplication)的嚴(yán)重干擾.因此，本文選擇臺(tái)陣至目標(biāo)源區(qū)網(wǎng)格點(diǎn)的距離作為雙臺(tái)陣聯(lián)合FB定位方法的歸一化因子(Ward Neale et al.，2018)，通過(guò)降低近源臺(tái)陣的權(quán)重占比以削弱上地幔三重震相的干擾，定義聯(lián)合臺(tái)陣(Combined array)導(dǎo)向矢量為

P(ω,lon,lat)=w(lon,lat)·e-iω Tlon,lat,

(4)

其中，w(lon,lat)為距離權(quán)重歸一化向量：

w(lon,lat)=(w1(lon,lat)…wN(lon,lat))?,

(5)

(6)

Dj(lon,lat)為當(dāng)前所搜索網(wǎng)格點(diǎn)(lon,lat)與臺(tái)站j之間的距離.

然而，如后文“3.3 臺(tái)陣屬性影響分析”所述，地震臺(tái)陣的組合可能導(dǎo)致聯(lián)合臺(tái)陣的臺(tái)陣響應(yīng)出現(xiàn)顯著的波紋狀割裂，嚴(yán)重干擾FB定位結(jié)果.因此，本文在上述聯(lián)合臺(tái)陣FB方法的基礎(chǔ)上，進(jìn)一步采用反卷積技術(shù)CLEAN-SC(CLEAN based on spatial Source Coherence；Sijtsma, 2007)，以抑制臺(tái)陣響應(yīng)對(duì)聯(lián)合定位結(jié)果的負(fù)面影響.CLEAN-SC技術(shù)本質(zhì)上是一種迭代“去噪”技術(shù)，我們將由式(1)、式(2)得到的Ct(ω)和Bt(ω,lon,lat)作為初始輸入C(0)(ω)和B(0)(ω,lon,lat)，后續(xù)迭代過(guò)程如下：

當(dāng)i≥1時(shí)，首先計(jì)算i-1時(shí)波束能量分布圖中的峰值B(i-1)(ω,lonmax,latmax)，其中(lonmax,latmax)為該峰值對(duì)應(yīng)的經(jīng)緯度，其對(duì)應(yīng)的源我們稱之為峰值源，它在地震臺(tái)陣引起的CDSMG(i)(ω)可表示為

G(i)(ω)=B(i-1)(ω,lonmax,latmax)

·(u(i)(ω)u?(i)(ω)-U(i)(ω)),

(7)

其中u(i)(ω)為峰值源相干分量(即為N×1的矩陣)，U(i)(ω)為u(i)(ω)u?(i)(ω)的對(duì)角線元素(即為N×N的矩陣).u(i)(ω)可設(shè)置初值為e-iω Tlonmax,latmax，并計(jì)算相應(yīng)的U(i)(ω)，然后通過(guò)下式反復(fù)迭代至收斂而得到最終值：

(8)

一般幾十次的迭代就可得到較為收斂的u(i)(ω).其中，我們通過(guò)去除波束能量分布圖中峰值源相干能量來(lái)減少臺(tái)陣響應(yīng)的干擾，即

B(i)(ω,lon,lat)=B(i-1)(ω,lon,lat)-P?(ω,lon,lat)C(i)(ω)P(ω,lon,lat),

(9)

C(i)(ω)=C(i-1)(ω)-φG(i)(ω),

(10)

其中φ為安全因子，取值范圍為(0,1].同時(shí)，為彌補(bǔ)波束能量的缺失，須以“潔凈”的波束取代該峰值源：

Q(i)(ω,lon,lat)=φB(i-1)(ω,lonmax,latmax)·10λ‖(lon,lat)-(lonmax,latmax)‖2,

(11)

其中λ為帶寬系數(shù)，決定“潔凈”波束的帶寬；‖‖表示模長(zhǎng)計(jì)算.

當(dāng)?shù)鶬次，滿足∑|C(i)(ω)|≥∑|C(i-1)(ω)|時(shí)，便可得到足夠“潔凈”的波束能量分布：

(12)

基于上述結(jié)合反卷積技術(shù)的聯(lián)合FB定位方法，本文利用ChinArray II和Hi-net對(duì)臺(tái)風(fēng)“海鷗”期間目標(biāo)海域(100°E—140°E，5°N—25°N)上的P波地脈動(dòng)源區(qū)進(jìn)行雙臺(tái)陣聯(lián)合定位.

2.2 臺(tái)風(fēng)激發(fā)P波地脈動(dòng)源區(qū)數(shù)值模擬

理論上，臺(tái)風(fēng)激發(fā)的P波倍頻地脈動(dòng)源區(qū)可表達(dá)為相向傳播的同頻率波列發(fā)生非線性相互作用后產(chǎn)生的海表等效壓力場(chǎng).根據(jù)Ardhuin和Herbers(2013)和Farra等(2016)的方法，該等效壓力場(chǎng)的功率譜FP可寫為：

FP(x,f)=[2π]2[ρwg]2fE2(x,f/2)I(x,f/2),

(13)

其中f為P波地脈動(dòng)信號(hào)頻率；x表示空間位置；ρw為海水密度(= 1.020 kg·m-3)；g為重力加速度；E(x,f/2)為海浪頻率譜；無(wú)量綱量I(x,f/2)則與海浪能量隨方向θ的分布特征，即方向譜M(x,f/2,θ)，有關(guān)：

(14)

考慮到場(chǎng)地放大效應(yīng)，我們進(jìn)一步引入P波場(chǎng)地放大系數(shù)CP(Gualtieri et al.，2014；Farra et al.，2016)：

(15)

其中θiw(x)為自源區(qū)位置x處傳向地震臺(tái)陣的P波射線出射角，取決于x點(diǎn)與地震臺(tái)陣之間的距離；h(x)為x處水深，取自ETOPO1測(cè)深數(shù)據(jù)(Amante and Eakins，2009)；TP(θiw(x))和RP(θiw(x))分別為海底P波透射和反射系數(shù)：

TP(θiw(x))=

(16)

(17)

(18)

其中ρc(x)表示x點(diǎn)海底地殼密度，αw(x)和αc(x)分別為x點(diǎn)海水及海底地殼中P波傳播速度，βc(x)則為x點(diǎn)海底地殼中S波傳播速度，均取值自全球地殼模型CRUST1.0(http:∥igppweb.uscd.edu/～gabi/rem.html)；p表示基于Snell定律的射線參數(shù).Φw(h(x),f,θiw(x))則是P波在水層中傳播引起的相移，可表示為

(19)

因此，最終P波地脈動(dòng)理論源區(qū)可表示為

(20)

根據(jù)上述方法，本文同時(shí)給出了臺(tái)風(fēng)“海鷗”期間目標(biāo)海域上的P波地脈動(dòng)理論源區(qū).由于目標(biāo)海域至ChinArray II和Hi-net的P波出射角θiw(x)不同，使得二者的CP存在細(xì)微差異，本文取二者CP的均值進(jìn)行計(jì)算.

3 結(jié)果與討論

3.1 雙臺(tái)陣聯(lián)合定位結(jié)果

基于上述雙臺(tái)陣聯(lián)合FB定位方法，本文利用ChinArray II和Hi-net組成的聯(lián)合臺(tái)陣對(duì)臺(tái)風(fēng)“海鷗”生命周期范圍內(nèi)連續(xù)16個(gè)時(shí)刻(UTC 2014/09/13 00∶00至09/16 18∶00，間隔為6 h)的P波地脈動(dòng)源區(qū)進(jìn)行了定位，圖2c1—c8按12 h為時(shí)間間隔對(duì)定位結(jié)果進(jìn)行了顯示，對(duì)應(yīng)頻率為～0.2 Hz(即臺(tái)風(fēng)激發(fā)的P波地脈動(dòng)信號(hào)的峰值頻率；Lin et al., 2018).同時(shí)，為評(píng)估雙臺(tái)陣聯(lián)合FB定位效果，本文對(duì)“海鷗”期間的P波地脈動(dòng)理論源區(qū)進(jìn)行了計(jì)算(圖2d1—d8)，以作對(duì)比驗(yàn)證.雙臺(tái)陣聯(lián)合FB定位成功探測(cè)到“海鷗”激發(fā)的P波地脈動(dòng)源區(qū)，且定位所得源區(qū)分布與理論計(jì)算結(jié)果具有相當(dāng)?shù)囊恢滦裕涸?月15日“海鷗”穿越呂宋島進(jìn)入南海之前，其激發(fā)的P波地脈動(dòng)源區(qū)主要分布于菲律賓海之上(圖2c1—c4，d1—d4)；在“海鷗”深入南海之后，P波地脈動(dòng)源區(qū)主要分布于南海并緊隨臺(tái)風(fēng)中心而移動(dòng)，而菲律賓海上的地脈動(dòng)源區(qū)逐漸衰退(圖2c5—c8，d5—d8).南海海域上定位結(jié)果與理論源區(qū)的一致性相對(duì)菲律賓海域較差：定位源區(qū)較為聚焦地分布于臺(tái)風(fēng)中心后方，而理論源區(qū)則相對(duì)離散地向南漂移.這可能與南海較為封閉的海域環(huán)境有關(guān)，其四周的復(fù)雜海岸環(huán)境使得理論源區(qū)的模擬難度增大.上述結(jié)果表明，本文基于距離權(quán)重的雙臺(tái)陣聯(lián)合FB定位方法在一定程度上實(shí)現(xiàn)了對(duì)“海鷗”激發(fā)的P波地脈動(dòng)源區(qū)的全程有效定位與追蹤.

圖2 臺(tái)風(fēng)“海鷗”期間(9月13日至16日)時(shí)間間隔為12 h的連續(xù)8個(gè)時(shí)刻的基于(a1—a8)ChinArray II、(b1—b8)Hi-net和(c1—c8)聯(lián)合臺(tái)陣(Combined array)的FB定位結(jié)果(對(duì)應(yīng)頻率～0.2 Hz)以及(d1—d8)相應(yīng)時(shí)刻“海鷗”激發(fā)P波地脈動(dòng)理論源區(qū)(單位為log10(Pa2·m2·s)).黑色空心圓圈表示當(dāng)前時(shí)刻臺(tái)風(fēng)中心位置，黑色實(shí)心圓點(diǎn)表示此前不同時(shí)刻臺(tái)風(fēng)中心軌跡，時(shí)間間隔為6 hFig.2 Localized P-wave microseism source regions (at ～0.2 Hz) by FB method with ChinArray II (a1—a8), Hi-net (b1—b8) and Combined array (c1—c8), and corresponding theoretical source regions (d1—d8, in log10(Pa2·m2·s)) at successive 8 time-points during the lifecycle of Typhoon Kalmaegi (09/13—09/16) with 12 h time intervals. The black circles represent the current typhoon center, with the black solid dots indicating the typhoon locations at 6 h time intervals along the typhoon track

同時(shí)，比較有趣的是，雙臺(tái)陣聯(lián)合定位結(jié)果以及理論源區(qū)顯示，臺(tái)風(fēng)“海鷗”行經(jīng)菲律賓海期間(09/13 12∶00至09/15 12∶00)，其激發(fā)的P波地脈動(dòng)源區(qū)較為穩(wěn)定地分布于其軌跡的右上方，并未隨臺(tái)風(fēng)中心移動(dòng)，甚至出現(xiàn)向東北漂移遠(yuǎn)離的趨勢(shì).這表明該P(yáng)波地脈動(dòng)源區(qū)可能并非單純由臺(tái)風(fēng)“海鷗”產(chǎn)生的波浪場(chǎng)激發(fā).鑒于臺(tái)風(fēng)“海鷗”的生成稍晚于臺(tái)風(fēng)“風(fēng)神”的消散，且二者均起源于菲律賓海(圖1)，“風(fēng)神”遺留的波浪場(chǎng)很可能傳播至“海鷗”行經(jīng)海域并與其產(chǎn)生的波浪場(chǎng)相互作用，從而激發(fā)出不隨“海鷗”移動(dòng)的P波地脈動(dòng)源區(qū).為驗(yàn)證這一猜想，本文選用“風(fēng)神”和“海鷗”交替前后連續(xù)12個(gè)時(shí)刻(UTC 2014/09/09 12∶00至09/15 00∶00，間隔為12 h)的ERA5總涌浪平均周期數(shù)據(jù)(圖3a—f，m—r)以及相應(yīng)時(shí)刻的海表等效壓力譜數(shù)據(jù)(圖3g—l，s—x)進(jìn)行對(duì)比分析.當(dāng)“風(fēng)神”沿日本島弧向東北移動(dòng)，分布于其臺(tái)風(fēng)中心附近的波浪場(chǎng)及海表等效壓力場(chǎng)也隨之移動(dòng)(09/09 12∶00后，圖3a—d，g—j)，但有部分風(fēng)生涌浪群卻逐步背向遠(yuǎn)離，向菲律賓群島傳播(圖3a—f)；當(dāng)這些涌浪群入射菲律賓群島沿岸時(shí)，“風(fēng)神”消散而“海鷗”逐漸生成壯大(09/11 12∶00至09/12 00∶00，圖3e—f)；當(dāng)臺(tái)風(fēng)“海鷗”開始向西北移動(dòng)，其產(chǎn)生的波浪場(chǎng)與菲律賓沿岸反射的涌浪群非線性相互作用而激發(fā)海表等效壓力場(chǎng)，且該海表壓力場(chǎng)隨涌浪群傳播而向東北擴(kuò)散漂移(09/13 00∶00后，圖3n—r，t—x).由于場(chǎng)地放大效應(yīng)的存在，源自該海表等效壓力場(chǎng)的P波地脈動(dòng)源區(qū)較穩(wěn)定地分布于“海鷗”軌跡的右上方，并呈現(xiàn)向東北漂移趨勢(shì).上述分析表明，菲律賓海上不隨臺(tái)風(fēng)“海鷗”移動(dòng)的P波地脈動(dòng)源區(qū)，是由“風(fēng)神”產(chǎn)生的涌浪群與“海鷗”產(chǎn)生的波浪場(chǎng)相互作用而產(chǎn)生.

為了評(píng)價(jià)本文所提出的基于距離權(quán)重的雙臺(tái)陣聯(lián)合FB定位方法，以下我們進(jìn)一步通過(guò)與單臺(tái)陣定位結(jié)果的對(duì)比、聯(lián)合臺(tái)陣的臺(tái)陣屬性影響等方面進(jìn)行了分析與討論.

3.2 相對(duì)單臺(tái)陣的定位優(yōu)勢(shì)分析

為評(píng)估上述雙臺(tái)陣聯(lián)合FB技術(shù)對(duì)定位精度的提升，本文分別利用ChinArray II和Hi-net的單臺(tái)陣數(shù)據(jù)對(duì)臺(tái)風(fēng)“海鷗”激發(fā)的P波地脈動(dòng)源區(qū)進(jìn)行了定位.圖2a1—a8和圖2b1—b8分別給出了兩個(gè)臺(tái)陣在“海鷗”生命周期范圍內(nèi)連續(xù)8個(gè)時(shí)刻(UTC 2014/09/13 00∶00至09/16 12∶00，間隔為12 h)的定位結(jié)果，對(duì)應(yīng)頻率為～0.2 Hz.單臺(tái)陣定位亦均成功探測(cè)到“海鷗”激發(fā)的P波地脈動(dòng)源區(qū)，但與理論源區(qū)(圖2d1—d8)相比，單臺(tái)陣的定位結(jié)果均存在一定缺陷：當(dāng)P波地脈動(dòng)源區(qū)位于菲律賓海之上時(shí)(即09/15 12∶00之前)，ChinArray II的定位結(jié)果雖然未能足夠刻畫相對(duì)較弱源區(qū)的分布但卻十分聚焦(圖2a1—a6)，而Hi-net的定位結(jié)果則顯得相對(duì)散亂且向南漂移(圖2b1—b6)；當(dāng)“海鷗”遷移至南海時(shí)(即09/15 00∶00之后)，ChinArray II的定位結(jié)果就變得相當(dāng)零散(圖2a6—a8)，而Hi-net的定位結(jié)果雖然略向北漂移但開始變得相對(duì)聚焦(圖2b6—b8).

對(duì)比上述兩個(gè)單臺(tái)陣的定位結(jié)果，我們發(fā)現(xiàn)一個(gè)有趣的現(xiàn)象：ChinArray II在菲律賓海上的定位結(jié)果比較聚焦，而在南海定位則相對(duì)渙散(圖2a1—a8)；相反，Hi-net在菲律賓海上的定位結(jié)果則比較渙散，而在南海定位相對(duì)聚焦(圖2b1—b8).這可能與地震臺(tái)陣至地脈動(dòng)源區(qū)的距離有關(guān)：當(dāng)臺(tái)陣距離目標(biāo)源區(qū)較近時(shí)，射線路徑分布較淺而受地球淺層介質(zhì)橫向非均勻性的影響較大，致使基于理想地球模型的理論預(yù)測(cè)走時(shí)與實(shí)際走時(shí)之間存在相對(duì)較大的偏差，導(dǎo)致臺(tái)陣的聚束定位結(jié)果不穩(wěn)定、定位能力變差(Fang et al.，2020).當(dāng)然，地震波傳播的“三重震相”(Triplication)也可能是干擾因素之一.P波在固體地球中傳播時(shí)，會(huì)在上地幔的兩個(gè)主要間斷面(約410和660 km)發(fā)生反射，產(chǎn)生不同相位的P波.在15°～26°距離范圍內(nèi)，這些不同相位的P波可能近乎同時(shí)傳播到地震臺(tái)站，造成臺(tái)站記錄的相位信息混雜，這種現(xiàn)象被稱為上地幔三重震相(Chu et al.，2019).當(dāng)“海鷗”激發(fā)的P波地脈動(dòng)源區(qū)與ChinArray II和Hi-net的距離在26°以內(nèi)時(shí)，兩個(gè)臺(tái)陣的FB定位結(jié)果可能因上地幔三重震相干擾而產(chǎn)生額外偏差.而得益于本文引入的距離權(quán)重，降低了聯(lián)合臺(tái)陣定位中近源臺(tái)陣的貢獻(xiàn)比，雙臺(tái)陣聯(lián)合FB定位結(jié)果很大程度上削弱了上述P波走時(shí)預(yù)測(cè)偏差及上地幔三重震相的干擾，為對(duì)“海鷗”激發(fā)的P波地脈動(dòng)源區(qū)進(jìn)行持續(xù)、有效定位與追蹤奠定了基礎(chǔ).

圖3 臺(tái)風(fēng)“風(fēng)神”和“海鷗”交替前后(9月9日至15日)時(shí)間間隔為12 h的連續(xù)12個(gè)時(shí)刻的(a—f，m—r)ERA5總涌浪平均周期(單位s)，和(g—l，s—x)相應(yīng)時(shí)刻海表等效壓力譜(單位為Pa2·m2·s)的空間分布圖.黑色空心圓圈表示當(dāng)前時(shí)刻臺(tái)風(fēng)中心位置，黑色實(shí)心圓點(diǎn)表示此前不同時(shí)刻臺(tái)風(fēng)中心軌跡，時(shí)間間隔為6 hFig.3 Spatial distributions of the mean periods of total swell of ERA5 (a—f, m—r; in s) and calculated equivalent surface pressure spectra (g—l, s—x; in Pa2·m2·s) at successive 12 time-points before and after the alternation of Typhoon Fengshen and Kalmaegi (09/09—09/15) with 12 h time intervals. The black circle represents the current typhoon center, with the black solid dots indicating the typhoon locations at 6 h time intervals along the typhoon track

理論上，地震臺(tái)陣的空間分辨率R取決于其孔徑Φ、地震波波長(zhǎng)λ以及與源區(qū)的距離D，可用R=λD/Φ表示(Rost and Thomas，2002；Ward Neale et al.，2018).例如，對(duì)于3000 km外的P波地脈動(dòng)源區(qū)，假設(shè)其傳播速度為6 km·s-1(0.2 Hz頻率對(duì)應(yīng)的波長(zhǎng)約為30 km)，則ChinArray II的空間分辨率約為62.4 km，而Hi-net則為47.1 km，即Hi-net的理論分辨能力優(yōu)于ChinArray II.該分辨率R一定程度上反映了地震臺(tái)陣對(duì)相鄰源區(qū)的探測(cè)和分辨能力.雖然，實(shí)際應(yīng)用中，基于FB技術(shù)的地震臺(tái)陣定位空間分辨率通常遠(yuǎn)不及理論值(Hillers et al.，2012；Ward Neale et al.，2018)，但雙臺(tái)陣聯(lián)合的空間分辨率仍相對(duì)單臺(tái)陣有了較大提升.例如，聯(lián)合定位結(jié)果清晰地顯示出“海鷗”穿越呂宋島后在其北部沿岸激發(fā)的較弱的近岸P波地脈動(dòng)源區(qū)(圖2c5和2d5)，而ChinArray II和Hi-net則無(wú)法單獨(dú)探測(cè)到(圖2a5和2b5).這進(jìn)一步體現(xiàn)了聯(lián)合臺(tái)陣相對(duì)單臺(tái)陣的空間分辨能力優(yōu)勢(shì).

此外，由于臺(tái)風(fēng)“海鷗”路徑軌跡在空間上具有相當(dāng)大的跨度，在緯向跨越～30°、經(jīng)向跨越～10°(圖1)，上述結(jié)果分析表明，單獨(dú)采用Hi-net或ChinArray II均無(wú)法對(duì)“海鷗”激發(fā)P波地脈動(dòng)源區(qū)進(jìn)行全程有效定位與追蹤.鑒于ChinArray II和Hi-net的單臺(tái)陣FB定位結(jié)果在空間域上具有一定的互補(bǔ)性，本文采用基于距離權(quán)重的雙臺(tái)陣聯(lián)合FB定位方法在一定程度上綜合了兩個(gè)臺(tái)陣的定位優(yōu)勢(shì)，提升了定位源區(qū)的穩(wěn)定性與精度，實(shí)現(xiàn)了對(duì)“海鷗”激發(fā)的P波地脈動(dòng)源區(qū)的全程、穩(wěn)定定位追蹤.

3.3 臺(tái)陣屬性影響分析

地震臺(tái)陣的固有屬性，包括其孔徑、幾何形狀、臺(tái)站數(shù)量及間距等，都會(huì)對(duì)臺(tái)陣定位的空間分辨率產(chǎn)生顯著的影響.由于臺(tái)陣可探測(cè)信號(hào)的最大波長(zhǎng)大致近似于其孔徑(Schweitzer et al.，2012)，能夠有效探測(cè)臺(tái)風(fēng)激發(fā)的P波地脈動(dòng)信號(hào)的臺(tái)陣孔徑至少為30 km(假設(shè)P波波速為6 km/s)；而臺(tái)陣的幾何形狀則決定了定位結(jié)果的優(yōu)勢(shì)空間分辨方向.這種臺(tái)陣固有屬性造成的影響可通過(guò)臺(tái)陣響應(yīng)函數(shù)(Array Response Function, ARF)進(jìn)行表征(Lin et al.，2018)：

(21)

其中k為波數(shù)矢量，k0為參考波數(shù)矢量，rn則表示地震臺(tái)站的位置，N為臺(tái)站數(shù)量.

圖4a和4b分別給出了ChinArray II和Hi-net的ARF.得益于兩個(gè)臺(tái)陣的大孔徑(ChinArray II的孔徑約為1443 km，Hi-net的孔徑約為1910 km)，二者的ARF主瓣都較為聚焦，將對(duì)目標(biāo)源區(qū)具有較強(qiáng)的定位分辨能力.受臺(tái)陣幾何形狀影響，ChinArray II的ARF主瓣較為規(guī)整，略微沿西南-東北向拉伸，表明其對(duì)西北-東南方向的空間分辨能力略優(yōu)；而Hi-net的ARF主瓣呈西北-東南向的梭形，其對(duì)西南-東北方向的空間分辨能力相對(duì)更強(qiáng).因此，臺(tái)風(fēng)“海鷗”在菲律賓海和南海產(chǎn)生的P波地脈動(dòng)源區(qū)正好總體大致處于這兩個(gè)臺(tái)陣的空間分辨優(yōu)勢(shì)方位角方向，即位于ChinArray II的東南方向和Hi-net的西南方向，這為聯(lián)合這兩個(gè)臺(tái)陣進(jìn)行有效定位提供了一定保障.聯(lián)合臺(tái)陣由ChinArray II和Hi-net組合而成，由于ChinArray II的臺(tái)站數(shù)目?jī)H為Hi-net的58.5%(圖1)，聯(lián)合臺(tái)陣的ARF總體呈近似Hi-net的ARF分布，而ChinArray II臺(tái)站的加入一定程度上削弱了Hi-net的ARF旁瓣(圖4d).遺憾的是，聯(lián)合臺(tái)陣的ARF呈現(xiàn)豎向波紋狀割裂，將使FB定位結(jié)果出現(xiàn)類似的波紋干擾.初步判斷，該現(xiàn)象可能與ChinArray II和Hi-net的臺(tái)站數(shù)目差異及空間距離有關(guān).為探究該現(xiàn)象的成因，本文進(jìn)行了聯(lián)合臺(tái)陣組合測(cè)試.

如圖4j所示，我們?cè)贖i-net中部截取了與ChinArray II臺(tái)站數(shù)目相同的臺(tái)陣，記為Array 1，并將其復(fù)制并移動(dòng)[-18°E, -20°N]得到Array 2，由于臺(tái)陣幾何結(jié)構(gòu)相同，二者的臺(tái)陣響應(yīng)如圖4c所示；同時(shí)，我們將ChinArray II分別復(fù)制移動(dòng)[+15°E, 0°N]和[-15°E, 0°N]得到Array 3和Array 4.我們分別計(jì)算這些臺(tái)陣的聯(lián)合臺(tái)陣響應(yīng)并進(jìn)行分析：由ChinArray II和Array 1的聯(lián)合臺(tái)陣響應(yīng)(圖4e)與圖4d比較可知,臺(tái)站數(shù)目并非造成波紋狀割裂的原因；Array 1和Array 3(圖4f)以及Array 1和Array 4(圖4g)的聯(lián)合臺(tái)陣響應(yīng)差異則表明，臺(tái)陣間距是造成波紋狀割裂現(xiàn)象的主要原因之一，且波紋割裂的密集程度與距離成正相關(guān)；而從Array 2和Array 3(圖4h)以及Array 2和Array 4(圖4i)的聯(lián)合臺(tái)陣響應(yīng)差異可見，臺(tái)陣的相對(duì)空間方位決定了波紋狀割裂的分布方向.上述分析表明，聯(lián)合臺(tái)陣的豎向波紋狀割裂由ChinArray II和Hi-net的近緯向間距導(dǎo)致.

圖4 地震臺(tái)陣(a)ChinArray II、(b)Hi-net、(c)Array 1、(d)聯(lián)合臺(tái)陣(Combined array)、(e)ChinArray II+Array 1聯(lián)合、(f)Array 1+Array 3聯(lián)合、(g)Array 1+Array 4聯(lián)合、(h)Array 2+Array 3聯(lián)合及(i)Array 2+Array 4聯(lián)合的響應(yīng)函數(shù)ARFs以及(j)相應(yīng)地理分布Fig.4 Array Response Functions (ARFs) of (a) ChinArray II, (b) Hi-net, (c) Array 1， (d) Combined array, (e) ChinArray II+Array 1, (f) Array 1+Array 3, (g) Array 1+Array 4, (h) Array 2+Array 3 and (i) Array 2+Array 4, and (j) corresponding geographical distribution

由于常規(guī)FB技術(shù)本身無(wú)法消除由地震臺(tái)陣固有屬性(如臺(tái)陣ARF所表征)的影響，本文在雙臺(tái)陣聯(lián)合FB定位技術(shù)的基礎(chǔ)上進(jìn)一步結(jié)合反卷積技術(shù)CLEAN-SC，盡可能地避免聯(lián)合定位結(jié)果出現(xiàn)波紋狀干擾.最終聯(lián)合定位結(jié)果符合預(yù)期，只有少數(shù)時(shí)刻存在淡淡的“波紋”(圖2c1、2c5和2c6)，且并不影響定位結(jié)果的精度，這體現(xiàn)了本文雙臺(tái)陣聯(lián)合FB定位技術(shù)的可靠性與優(yōu)勢(shì).

4 結(jié)論

本文基于ChinArray II和Hi-net的連續(xù)波形數(shù)據(jù)，對(duì)2014年西北太平洋強(qiáng)臺(tái)風(fēng)“海鷗”激發(fā)的P波地脈動(dòng)源區(qū)進(jìn)行定位研究.我們發(fā)現(xiàn)，臺(tái)風(fēng)“海鷗”在菲律賓海激發(fā)的P波地脈動(dòng)主要源自其激發(fā)海浪與臺(tái)風(fēng)“風(fēng)神”遺留涌浪的相互作用，而其本身在不同時(shí)刻激發(fā)海浪之間相互作用而激發(fā)的地脈動(dòng)能量相對(duì)較弱，這導(dǎo)致主要源區(qū)逐漸偏離“海鷗”移動(dòng)軌跡.該現(xiàn)象表明由其他風(fēng)暴引起的涌浪可能“干擾”對(duì)特定臺(tái)風(fēng)的定位追蹤，值得后續(xù)深入研究；同時(shí)，南海作為“相對(duì)封閉”的邊緣海，東部島鏈的存在能有效“攔截”來(lái)自西太平洋風(fēng)暴引起的涌浪的影響，不失為難得“干凈”的臺(tái)風(fēng)激發(fā)地脈動(dòng)觀測(cè)、研究環(huán)境.

由于ChinArray II、Hi-net單臺(tái)陣定位結(jié)果受多重因素的影響，包括：(1)ARF強(qiáng)旁瓣導(dǎo)致定位結(jié)果出現(xiàn)“偽影”、不聚焦；(2)近距離時(shí)(即臺(tái)陣至源區(qū)距離小于26°)基于理論地球模型的P波走時(shí)估計(jì)與實(shí)際走時(shí)之間的偏差以及(3)上地幔三重震相產(chǎn)生的干擾等，本文提出了將臺(tái)陣至源區(qū)的距離作為歸一化因子、結(jié)合反卷積技術(shù)CLEAN-SC的雙臺(tái)陣聯(lián)合FB定位方法.該方法能夠削弱單臺(tái)陣定位應(yīng)用中以上方面的影響，有效提升定位結(jié)果的聚焦性與穩(wěn)定性，實(shí)現(xiàn)了對(duì)臺(tái)風(fēng)“海鷗”在菲律賓海、南海激發(fā)的P波地脈動(dòng)源區(qū)的全程、穩(wěn)定有效定位與追蹤.

然而，本文的雙臺(tái)陣聯(lián)合FB定位研究仍存在一些不足之處.本文的聯(lián)合定位結(jié)果與理論源區(qū)在南海海域的差異較大，暫無(wú)法判別優(yōu)劣與原因，亟待后續(xù)進(jìn)一步研究.另外，雖然本文利用3-D地球介質(zhì)速度模型提高了P波走時(shí)估計(jì)精度，并且能夠通過(guò)雙臺(tái)陣聯(lián)合優(yōu)化近場(chǎng)源區(qū)定位效果，但未能徹底避免上地幔三重震相及近源P波走時(shí)估計(jì)偏差對(duì)定位結(jié)果的影響.上述問(wèn)題的存在仍會(huì)對(duì)雙臺(tái)陣聯(lián)合定位的精確性和可靠性造成影響，有待后續(xù)進(jìn)一步研究.

致謝謹(jǐn)此祝賀陳颙先生從事地球物理教學(xué)科研工作60周年.感謝中國(guó)地震科學(xué)探測(cè)臺(tái)陣數(shù)據(jù)中心和日本Hi-net臺(tái)網(wǎng)為本研究提供連續(xù)地震記錄數(shù)據(jù).感謝日本氣象廳提供臺(tái)風(fēng)最佳路徑數(shù)據(jù)，歐洲中期氣象預(yù)報(bào)中心提供波浪場(chǎng)數(shù)據(jù)，法國(guó)海洋開發(fā)研究院提供的P波地脈動(dòng)源區(qū)海表等效壓力譜數(shù)據(jù).