孫 楊， 姜效典， 王 燚， 李朝陽(yáng)(中國(guó)海洋大學(xué) .海洋地球科學(xué)學(xué)院；.計(jì)算機(jī)基礎(chǔ)部，山東 青島 6600)

衛(wèi)星測(cè)高反演重力場(chǎng)在西太平洋海域的應(yīng)用?

孫 楊1， 姜效典1， 王 燚2， 李朝陽(yáng)1
(中國(guó)海洋大學(xué) 1.海洋地球科學(xué)學(xué)院；2.計(jì)算機(jī)基礎(chǔ)部，山東 青島 266100)

利用高精度衛(wèi)星測(cè)高數(shù)據(jù)反演重力是獲取海洋重力異常的重要方法。本文介紹了自2008年以來(lái)新一代高精度測(cè)高衛(wèi)星的相關(guān)信息。在此基礎(chǔ)上，聯(lián)合多顆新型衛(wèi)星，進(jìn)行交叉點(diǎn)平差，通過(guò)中心差分代替常規(guī)Sandwell方法計(jì)算垂線偏差，利用改進(jìn)的逆Vening-Meinesz法和Stokes數(shù)值反解法對(duì)西太平洋某區(qū)域進(jìn)行了重力場(chǎng)反演。反演結(jié)果與NGDC提供的船測(cè)重力值比較表明，逆Vening-Meinesz法計(jì)算誤差小于Stokes數(shù)值反解法。通過(guò)添加隨機(jī)噪聲數(shù)值實(shí)驗(yàn)，進(jìn)一步證明逆Vening-Meinesz法對(duì)于測(cè)高數(shù)據(jù)誤差的抗差能力更強(qiáng)，在研究區(qū)的重力異常計(jì)算結(jié)果可信度較高。

測(cè)高衛(wèi)星；海洋重力；逆Vening-Meinesz；垂線偏差

衛(wèi)星測(cè)高技術(shù)作為大地測(cè)量的重要手段，具有全天候，高重復(fù)率，高精度，低成本的特點(diǎn)。近年來(lái)隨著測(cè)高數(shù)據(jù)精度的不斷提高，其在地球重力場(chǎng)推算，尤其是快速推求大區(qū)域重力數(shù)據(jù)過(guò)程中的突出優(yōu)勢(shì)越來(lái)越引起人們的重視。與船測(cè)重力異常相比，衛(wèi)星反演重力數(shù)據(jù)穩(wěn)定性好，所有數(shù)據(jù)都基于同一基準(zhǔn)面觀測(cè)，基本不受參考框架統(tǒng)一性問(wèn)題的限制，避免了不同船測(cè)區(qū)域間需要數(shù)據(jù)融合的困擾。國(guó)內(nèi)外眾多學(xué)者對(duì)聯(lián)合多顆衛(wèi)星測(cè)高數(shù)據(jù)進(jìn)行重力場(chǎng)反演開(kāi)展了大量工作，萬(wàn)劍華利用Cryosat-2衛(wèi)星數(shù)據(jù)在南海進(jìn)行了重力異常反演，分辨率達(dá)到了2’×2’，精度達(dá)到了4.5×10-5m/s2[1]。彭富清，李洋等分析了衛(wèi)星測(cè)高誤差對(duì)海洋重力場(chǎng)反演精度的影響，指出反演誤差與測(cè)高隨機(jī)誤差呈正比變化，與空間分辨率成反比變化[2-3]。Sandwell聯(lián)合Cryosat-2，Envisat和Jason-1衛(wèi)星在墨西哥灣進(jìn)行了重力反演，精度達(dá)到了2×10-5m/s2[4]。在此基礎(chǔ)上，建立了新一代的海洋重力模型，分析了目前衛(wèi)星重力反演存在的主要問(wèn)題，指出了未來(lái)的研究方向[5]。此外，Hwang[6]，Keating[7]和Garcia[8]等也先后開(kāi)展了針對(duì)淺水海域衛(wèi)星重力反演的研究，并取得了卓有成效的結(jié)果。鑒于針對(duì)新型測(cè)高衛(wèi)星數(shù)據(jù)聯(lián)合反演的系統(tǒng)研究較少，本文以近年來(lái)多顆新測(cè)高衛(wèi)星為研究對(duì)象，對(duì)西太平洋某海域進(jìn)行了重力場(chǎng)的計(jì)算，并與NGDC船測(cè)數(shù)據(jù)對(duì)比分析誤差結(jié)果。

1 數(shù)據(jù)來(lái)源與研究區(qū)概況

利用衛(wèi)星測(cè)高數(shù)據(jù)反演得到的重力場(chǎng)精度主要依賴于四個(gè)因素，即軌道的空間密度，測(cè)高計(jì)高程精度，異軌數(shù)據(jù)的方向(方位角)，以及近岸涌浪模型的精度。近幾年，多顆新型測(cè)高衛(wèi)星進(jìn)入太空，同時(shí)有已服役衛(wèi)星進(jìn)行了軌道調(diào)整，大量新觀測(cè)數(shù)據(jù)的獲得,為進(jìn)行高精度，高分辨率的衛(wèi)星重力反演提供了可能[3]。下表列出了這些衛(wèi)星的主要參數(shù)。

其中，Envisat New和Jason-1 Geodetic分別在執(zhí)行了8和11年任務(wù)之后實(shí)施了變軌操作，而Jason-2，Cryosat-2，Altika和HY-2則發(fā)射于2008—2014的6年間。

研究區(qū)位于西太平洋帕里西維拉海盆和馬里亞納弧盆系統(tǒng)的接合部(140.5°E～143.5°E，15°N～21.5°N)，西部為帕里西維拉海盆，水深較深，東部為西馬里亞納島弧，水深較淺。區(qū)中有2006年實(shí)測(cè)的重力觀測(cè)數(shù)據(jù)，由NGDC公布，編號(hào)BD0601，共35條主測(cè)線，3條聯(lián)絡(luò)測(cè)線，214 189個(gè)采樣點(diǎn)。船測(cè)重力異常最小值為-13.29×10-5m/s2，最大為148.55×10-5m/s2，平均值為22.8×10-5m/s2。

多星聯(lián)合進(jìn)行重力計(jì)算是提高分辨率的重要手段。研究區(qū)不同時(shí)期衛(wèi)星地面軌跡如圖2所示。從圖中可以看出，不同時(shí)期衛(wèi)星軌跡密度有較大差別，在一定程度上會(huì)影響最終重力反演的結(jié)果。

2 重力異常計(jì)算

不同衛(wèi)星之間存在參考框架，測(cè)高精度以及各種校正模型的系統(tǒng)偏差，在進(jìn)行重力反演之前需要先進(jìn)行交點(diǎn)平差計(jì)算，減小系統(tǒng)誤差。之后利用測(cè)高數(shù)據(jù)和海面地形數(shù)據(jù)計(jì)算大地水準(zhǔn)面高，從而可以得到垂線偏差數(shù)據(jù)。在此基礎(chǔ)上，利用逆Vening-Meinesz公式和“移去-恢復(fù)”技術(shù)即可反演出研究區(qū)的自由空間重力異常。下面介紹具體的計(jì)算過(guò)程。

2.1 交叉點(diǎn)平差

交點(diǎn)平差可以減弱衛(wèi)星軌道徑向誤差造成的影響[14-16]。若將交叉點(diǎn)海面高的不符值作為觀測(cè)量，采用線性模型模擬軌道誤差，可建立平差觀測(cè)方程

(1)

(2)

2.2 大地水準(zhǔn)面確定

將得到的平差后的海面高數(shù)據(jù)利用Kriging法進(jìn)行網(wǎng)格化，然后減去由AVISO提供的平均海平面MSS_CNES_CLS11模型校正后的動(dòng)態(tài)海面地形網(wǎng)格數(shù)據(jù)，即可得到大地水準(zhǔn)面高。下圖給出了西太平洋地區(qū)10′×10′的大地水準(zhǔn)面分布。

2.3 計(jì)算垂線偏差

Sandwell和Olgiati等先后提出了不同的垂線偏差計(jì)算方法[17-19]。兩種方法得到的數(shù)據(jù)不規(guī)則地分布在觀測(cè)點(diǎn)和交叉點(diǎn)上，須額外將其處理成網(wǎng)格點(diǎn)的平均值。Sandwell方法利用了星下點(diǎn)速度參量，Olgiati方法則需要內(nèi)插垂直于軌跡的垂線偏差，這些步驟都引入了不必要的數(shù)值誤差。此外，褚坤曾嘗試?yán)么蟮厮疁?zhǔn)面的網(wǎng)格數(shù)據(jù)直接計(jì)算垂線偏差[8]。

為了建立快速，準(zhǔn)確的計(jì)算方法，本文從直觀物理意義出發(fā)，利用大地水準(zhǔn)面空間差分形式計(jì)算垂線偏差。網(wǎng)格點(diǎn)(i,j)處的垂線偏差子午和卯酉分量分別為

(3)

(4)

式中：ξ為垂線偏差子午分量；η為垂線偏差卯酉分量；N為大地水準(zhǔn)面高；R為地球平均半徑；φ為大地緯度；λ為大地經(jīng)度。在邊界處，可以相應(yīng)的用前向差分和后向差分代替中心差分。通過(guò)這種方法，可以快速的進(jìn)行垂線偏差的計(jì)算。圖4為利用差分法計(jì)算得到的10′×10′垂線偏差子午和卯酉分量。

2.4 計(jì)算重力異常

根據(jù)衛(wèi)星測(cè)高數(shù)據(jù)計(jì)算重力異常目前常用的有逆Vening-Meinesz方法和Stokes數(shù)值反解法。本文主要利用改進(jìn)的逆Vening-Meinesz方法計(jì)算重力異常。其二維球面譜計(jì)算式為[22]

(5)

其中：H′(φ，λ)為積分核函數(shù)H(φ，λ)對(duì)球面角距φ的導(dǎo)數(shù)；α為積分點(diǎn)到計(jì)算點(diǎn)之間的方位角。當(dāng)Ψ=0時(shí)，核函數(shù)H′(Ψ)將發(fā)生奇異。應(yīng)當(dāng)單獨(dú)考慮計(jì)算點(diǎn)所在方塊數(shù)據(jù)的影響，并對(duì)積分核函數(shù)做出相應(yīng)的改進(jìn)。其內(nèi)環(huán)帶的影響量可以表示為

(6)

其中：S0為內(nèi)環(huán)帶的球冠半徑；ξx和ηy分別表示ξ和η在x和y方向上的變化率。

與逆Vening-Meinesz法不同，Stokes反解法直接將測(cè)高數(shù)據(jù)作為輸入量，其二維譜計(jì)算式為

(7)

其中：N為大地水準(zhǔn)面高；S(φ，λ)為Stokes函數(shù)。

2.5 移去-恢復(fù)技術(shù)

在計(jì)算過(guò)程中，還需要采用“移去-恢復(fù)”技術(shù)，消除長(zhǎng)波長(zhǎng)數(shù)據(jù)造成的影響[4]。通過(guò)EGM2008模型球諧系數(shù)計(jì)算模型垂線偏差[20-24]。從2.3步得到的垂線偏差中減去模型垂線偏差，即為殘余垂線偏差，進(jìn)而計(jì)算得出殘余重力異常。殘余重力異常與EGM2008模型重力異常疊加，即為最終重力異常[17]。

3 實(shí)際資料處理

采用逆Vening-Meinesz法，聯(lián)合圖2右圖所示的Jason-2，Cryosat-2，Altika和HY-2四顆衛(wèi)星對(duì)研究區(qū)進(jìn)行自由空間重力異常的反演。結(jié)果如下圖所示。

從圖中可以看出，研究區(qū)以正異常為主，正異常高值集中在東側(cè)西馬里亞納島弧處，最大值為130.5×10-5m/s2。西馬里亞納島弧西側(cè)異常變化劇烈，沿島弧分布有南北向的異常梯度帶，顯示此處地殼結(jié)構(gòu)變化劇烈。最小值位于研究區(qū)中部， 為-11.6×10-5m/s2。進(jìn)一步選取圖5中的AA’和BB’兩條測(cè)線，并利用Stokes數(shù)值反解法計(jì)算測(cè)線上的重力異常，最后將船測(cè)重力異常與對(duì)應(yīng)點(diǎn)兩種方法的反演結(jié)果進(jìn)行對(duì)比。如下圖所示，圖中縱坐標(biāo)為自由空間重力異常，橫坐標(biāo)分別為測(cè)線上點(diǎn)距A點(diǎn)或B點(diǎn)的水平距離。

從圖中可以看出，兩種方法得到的衛(wèi)星重力異常與船測(cè)重力異常波動(dòng)趨勢(shì)一致，但都較為平滑。測(cè)高衛(wèi)星數(shù)據(jù)的主要信息集中在中長(zhǎng)波段，對(duì)于波長(zhǎng)小于10km的短波特征的表現(xiàn)較弱。整體來(lái)看，逆Vening-Meinesz法結(jié)果與實(shí)測(cè)重力更接近。而且對(duì)于細(xì)節(jié)的反映要好于Stokes數(shù)值反解法。具體到測(cè)線上，AA’測(cè)線上兩種方法的反演結(jié)果要優(yōu)于BB’測(cè)線。在BB’測(cè)線上，Stokes數(shù)值反解的結(jié)果已經(jīng)與實(shí)測(cè)重力有較大的誤差。對(duì)研究區(qū)和兩條測(cè)線的數(shù)據(jù)分別統(tǒng)計(jì)誤差信息如下表所示。

誤差統(tǒng)計(jì)表明，對(duì)于整個(gè)研究區(qū)，逆Vening-Meinesz法精度在5.77×10-5m/s2左右。兩種方法在AA’測(cè)線的結(jié)果要好于BB’測(cè)線的結(jié)果。這主要由兩個(gè)原因造成的，一是兩種方法在重力數(shù)據(jù)變化劇烈的區(qū)域反解精度要低于數(shù)值變化平緩的區(qū)域。另一方面，BB’測(cè)線附近衛(wèi)星測(cè)高數(shù)據(jù)精度較差。由于此處水深較淺，衛(wèi)星測(cè)高數(shù)據(jù)受回波影響精度大大降低，同時(shí)一些校正量如大氣濕對(duì)流層校正、海洋潮汐校正以及大氣高頻因數(shù)校正在淺海不準(zhǔn)確，也會(huì)對(duì)結(jié)果造成影響。AA’測(cè)線處水深在4 500m左右，衛(wèi)星數(shù)據(jù)質(zhì)量高，逆Vening-Meinesz法誤差約為1.45×10-5m/s2?？紤]到船測(cè)重力也存在(1～3)×10-5m/s2左右的偏差，利用衛(wèi)星測(cè)高數(shù)據(jù)計(jì)算得到的重力異常在該區(qū)域已經(jīng)接近于普通科考船的數(shù)據(jù)精度。

通過(guò)兩種方法間誤差對(duì)比可以看出，逆Vening-Meinesz法的精度要高于Stokes數(shù)值反解法。在AA’測(cè)線處逆Vening-Meinesz法精度要高出2.76×10-5m/s2，在BB’測(cè)線處精度更是高出4.26×10-5m/s2，說(shuō)明其對(duì)于測(cè)高數(shù)據(jù)的抗差能力更強(qiáng)。為進(jìn)一步論證該結(jié)論，將BB’測(cè)線的測(cè)高數(shù)據(jù)添加平均值為0，均方差為0.05m的隨機(jī)誤差。分別計(jì)算兩種方法受隨機(jī)誤差的影響。結(jié)果如圖7所示。

從圖7中可以看出，逆Vening-Meinesz法和Stokes數(shù)值反解法都一定程度地受到噪聲的影響。噪聲數(shù)據(jù)通過(guò)逆Vening-Meinesz法產(chǎn)生誤差的標(biāo)準(zhǔn)差為8.263 4×10-5m/s2，Stokes數(shù)值反解法的誤差則為17.330 6×10-5m/s2?？梢?jiàn)，逆Vening-Meinesz法的穩(wěn)定性要優(yōu)于Stokes數(shù)值反解法。

研究區(qū)內(nèi)既存在重力場(chǎng)劇烈變化的區(qū)域，也包含了淺水區(qū)。這些不利影響對(duì)利用衛(wèi)星測(cè)高數(shù)據(jù)反解重力場(chǎng)造成了很大的困難。逆Vening-Meinesz法可以有效地減弱測(cè)高誤差對(duì)于結(jié)果的影響，同時(shí)劇烈變化的重力場(chǎng)對(duì)其造成的干擾也相對(duì)較小。因此，利用逆Vening-Meinesz法得到的研究區(qū)重力異常數(shù)據(jù)具有較高的精度，可以很好地反映區(qū)域內(nèi)重力異常的空間分布。

4 結(jié)語(yǔ)

本文在介紹近年來(lái)新發(fā)射的多顆測(cè)高衛(wèi)星的基礎(chǔ)上，以改進(jìn)的逆Vening-Meinesz方法為基礎(chǔ)，引入一階中間差分算法代替Sandwell垂線偏差法，計(jì)算了西太平洋某海區(qū)重力異常。同時(shí)利用Stokes數(shù)值反解法計(jì)算的測(cè)線結(jié)果和NGDC提供的船測(cè)實(shí)測(cè)數(shù)據(jù)與其進(jìn)行了對(duì)比。結(jié)果顯示逆Vening-Meinesz法受重力場(chǎng)變化劇烈程度影響較小，計(jì)算結(jié)果與實(shí)測(cè)重力異常的符合度更高，很好地反映了研究區(qū)內(nèi)中長(zhǎng)波重力異常分布。

數(shù)值實(shí)驗(yàn)進(jìn)一步證明，相對(duì)于Stokes數(shù)值反解對(duì)測(cè)高數(shù)據(jù)的抗差能力更強(qiáng)。由于目標(biāo)區(qū)域內(nèi)包含淺水海域，衛(wèi)星測(cè)高數(shù)據(jù)誤差增大，故利用逆Vening-Meinesz法得到的重力場(chǎng)分布更加準(zhǔn)確。

衛(wèi)星測(cè)高數(shù)據(jù)反演在獲取大區(qū)域海洋重力場(chǎng)分布中正起到越來(lái)越重要的作用。針對(duì)重力場(chǎng)變化劇烈和淺水海域等特殊研究區(qū)，逆Vening-Meinesz法可以提供更好的重力異常數(shù)據(jù)。

致謝：感謝法國(guó)國(guó)家空間研究中心(CNES)提供衛(wèi)星數(shù)據(jù)。

[1] 萬(wàn)劍華, 李家軍, 劉善偉, 等. 基于Cryosat-2 數(shù)據(jù)的南海 2'× 2'重力異常計(jì)算與分析[J]. 中國(guó)石油大學(xué)學(xué)報(bào), 2015, 39(3): 70-75.WanJH,LiJJ,LiuSW,etal.Calculationandanalysisof2′×2′gravityanomaliesovertheSouthChinaSeabasedonCryosat-2satellitealtimeterdata[J].JournalofChinaUniversityofPetroleum, 2015, 39(3): 70-75.

[2] 彭富清, 陳雙軍, 金群峰. 衛(wèi)星測(cè)高誤差對(duì)海洋重力場(chǎng)反演的影響[J]. 測(cè)繪學(xué)報(bào), 2014, 43(4): 337-340.PengFQ,ChenSJ,JinQF.Influenceofaltimetryerrorsonmarinegeopotentialrecovery[J].ActaGeodaeticaetCartographicaSinica, 2014, 43(4): 337-340.

[3] 李洋, 張潤(rùn)寧. 高度計(jì)測(cè)距精度對(duì)沿軌跡重力異常反演的影響[J]. 測(cè)繪學(xué)報(bào), 2015, 44(4):G363-G369.LiY,ZhangRN.Inferenceofaltimeteraccuracyonalong-trackgravityanomalyrecovery[J].ActaGeodaeticaetCartographicaSinica, 2015, 44(4):G363-G369.

[4]SandwellD,GarciaE,SoofiK,etal.Toward1-mGalaccuracyinglobalmarinegravityfromCryoSat-2,Envisat,andJason-1[J].TheLeadingEdge, 2013, 32(8): 892-899.

[5]SandwellDT,MullerRD,SmithWHF,etal.NewglobalmarinegravitymodelfromCryoSat-2andJason-1revealsburiedtectonicstructure[J].Science, 2014, 346(6205): 65-67.

[6]HwangC,HsuHY.Shallow-watergravityanomaliesfromsatellitealtimetry:CasestudiesintheeastchinaseaandTaiwanstrait[J].JournaloftheChineseInstituteofEngineers, 2008, 31(5): 841-851.

[7]KeatingP,PinetN.Comparisonofsurfaceandshipbornegravitydatawithsatellite-altimetergravitydatainHudsonBay[J].TheLeadingEdge, 2013, 32(4): 450-458.

[8]GarciaES,SandwellDT,SmithWHF.RetrackingCryoSat-2,EnvisatandJason-1radaraltimetrywaveformsforimprovedgravityfieldrecovery[J].GeophysicalJournalInternational, 2014:ggt469.

[9]JiangX,LinM,LiuJ,etal.TheHY-2satelliteanditspreliminaryassessment[J].InternationalJournalofDigitalEarth, 2012, 5(3): 266-281.

[10] 彭海龍, 林明森, 穆博, 等.HY-2A衛(wèi)星雷達(dá)高度計(jì)數(shù)據(jù)的全球統(tǒng)計(jì)評(píng)價(jià)及質(zhì)量分析[J]. 海洋學(xué)報(bào), 2015(7): 54-66.PengHL,LinMS,MuB,etal.HY-2AGlobalstatisticalevaluationandperformanceanalysisofHY-2Asatelliteradaraltimeterdata[J].ActaOceanologicaSinica, 2015(7): 54-66.

[11] 楊磊, 周興華, 彭海龍, 等. 基于Jason-2 的Saral/AltiKa高度計(jì)全球統(tǒng)計(jì)評(píng)估與交叉定標(biāo)[J]. 海洋科學(xué)進(jìn)展, 2014, 32(4): 482-490.YangL,ZhouXH,PengHL,etal.Globalassessmentandcross-calibrationofSaral/AltikabasedonJason-2Altimeter[J].AdvancesinMarineScience, 2014, 32(4): 482-490.

[12]WinghamDJ,FrancisCR,BakerS,etal.CryoSat:AmissiontodeterminethefluctuationsinEarth’slandandmarineicefields[J].AdvancesinSpaceResearch, 2006, 37(4): 841-871.

[13]ChaudharyA,BasuS,KumarR,etal.RetrievingsignificantwaveheightintheIndianOceannearVisakhapatnamusingJason-2altimeterdata[J].RemoteSensingLetters, 2015, 6(4): 286-294.

[14] 褚坤.Geosat高度計(jì)數(shù)據(jù)處理與南海重力異常反演精度評(píng)價(jià)[D]. 青島: 中國(guó)石油大學(xué), 2010.ChuK.GeosatAltimeterDataProcessingandAccuracyAssessmentontheInversionofGravityAnomaliesintheSouthSea[D].Qingdao:ChinaUniversityofPetroleum, 2010.

[15] 王虎彪, 王勇, 陸洋, 等. 用衛(wèi)星測(cè)高和船測(cè)重力資料聯(lián)合反演海洋重力異常[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2005, 25(1): 81-85.WangHB,WangY,LuY,etal.Inversionofmarinegravityanomaliesbycombinatingmulti-altimeterdataandshipbornegravimetricdata[J].JournalofGeodesyandGeodynamics, 2005, 25(1): 81-85.

[16] 黃謨濤, 翟國(guó)君, 歐陽(yáng)永忠, 等. 利用多代衛(wèi)星測(cè)高數(shù)據(jù)反演海洋重力場(chǎng)[J]. 測(cè)繪科學(xué), 2006, 31(6): 37-39.HuangMT,ZhaiGJ,OuyangYZ,etal.Recoveryofmarinegravityfieldusingintegrateddatafrommulti-satellitemissions[J].ScienceofSurveyingandMapping, 2006, 31(6): 37-39.

[17] 王莉娟. 多源衛(wèi)星高度計(jì)中國(guó)近海及其鄰域重力異常反演與精度檢驗(yàn)[D]. 青島: 中國(guó)石油大學(xué), 2012.WangLJ.TheInversionoftheGravityAnomaliesandtheAssessmentofItsAccuracyBasedonMulti-satelliteAltimeterDatainChinaSeaandItsVicinity[D].Qingdao:ChinaUniversityofPetroleum, 2012.

[18]SandwellDT.AdetailedviewoftheSouthPacificgeoidfromsatellitealtimetry[J].JournalofGeophysicalResearch:SolidEarth, 1984, 89(B2): 1089-1104.

[19]OlgiatiA,BalminoG,SarrailhM,etal.Gravityanomaliesfromsatellitealtimetry:comparisonbetweencomputationviageoidheightsandviadeflectionsofthevertical[J].Bulletingéodésique, 1995, 69(4): 252-260.

[20] 黃謨濤. 海洋重力場(chǎng)測(cè)定及其應(yīng)用 [M]. 北京： 測(cè)繪出版社, 2005.HuangMT.DeterminationandApplicationofMarineGravityField[M].Beijing：SurveyingandMappingPress, 2005.

[21] 胡淑梅, 文漢江, 李洪超, 等. 利用多種測(cè)高數(shù)據(jù)反演中國(guó)南海海域重力異常[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2011, 31(4): 56-59.HuSM,WenHJ,LiHC,etal.InversionofgravityanomaliesoverSouthChinaSeabyuseofcombinationofmulti-satellitealtimeterdata[J].JournalofGeodesyandGeodynamics, 2011, 31(4): 56-59.

[22] 孫鳳華, 郭玉良, 陳春旺, 等. 我國(guó)沿海及海域 1'× 1'垂線偏差模型的研究與建立[J]. 海洋測(cè)繪, 2002, 22(5): 12-15.SunFH,GuoYL,ChenCW,etal.Researchandconstructionof1' × 1'deflectionoftheverticalmodelforChinaseaarea[J].HydrographicSurveyingandCharting, 2002, 22(5): 12-15.

[23] 黃謨濤, 翟國(guó)君, 管錚, 等. 利用衛(wèi)星測(cè)高數(shù)據(jù)反演海洋重力異常研究[J]. 測(cè)繪學(xué)報(bào), 2001, 30(2): 179-184.HuangMT,ZhaiGJ,GuanZ,etal.Ontherecoveryofgravityanomaliesfromaltimeterdata[J].ActaGeodaeticaetCartographicaSinica, 2001, 30(2): 179-184.

[24] 楊樂(lè), 林明森, 張有廣,等. 中國(guó)近岸海域高度計(jì)JASON-1測(cè)量數(shù)據(jù)的波形重構(gòu)算法研究[J]. 海洋學(xué)報(bào), 2010, 32(6):91-101.YangL,LinMS,ZhangYG,etal.ImprovingthequalityofJASON-1altimetrydatabywaveformretrackingincoastalwatersoffChina[J].ActaOceanologicaSinica, 2010, 32(6):91-101.

[25]SandwellDT,SmithWHF.RetrackingERS-1altimeterwaveformsforoptimalgravityfieldrecovery[J].GeophysicalJournalInternational, 2005, 163(1): 79-89.

責(zé)任編輯 徐 環(huán)

Application of Gravity Inversion From Satellite Altimeter Data in the West Pacific

SUN Yang1，JIANG Xiao-Dian1，WANG Yi2，LI Chao-Yang1

(Ocean University of China 1. College of Marine Geosciences; 2. Department of Comupter Foudation，Qingdao 266100， China)

The application of gravity anomaly inversion from high quality satellite altimeter data is an important method to acquire marine gravity. The main characteristics of newly available satellite altimeter data since 2008 are reviewed in this paper. On the basis of multi-satellite data, we apply crossover adjustment and use central difference method to calculate deflection of the vertical instead of conventional Sandwell method.After that, improved inverse Vening-Meinesz method and numerical inverse-Stokes approach are adopted to calculate the gravity filed of certain area in the West Pacific.Comparision with the shipboard-based gravity provided by NGDC shows that inverse Vening-Meinesz method has smaller errors than numerical inverse-Stokes approach.With further numerical expericment ofadding stochastic noise, we show that inverse Vening-Meinesz method is more robust so that it can provide reliable gravity estimation for the target area.

Satellinte altimeter; marine gravity; inverse Vening-Meinesz; deflection of the vertical

海洋公益性行業(yè)科研專項(xiàng)(201305029-02)資助 Supported by the Public Science and Technology Research Funds Projects of Ocean (201305029-02)

2015-12-07;

2016-03-07

孫楊(1987-)，男，博士。E-mail：sunyang_ouc@163.com

P312.1

A

1672-5174(2017)06-087-08

10.16441/j.cnki.hdxb.20150411

孫楊，姜效典，王燚,等.衛(wèi)星測(cè)高反演重力場(chǎng)在西太平洋海域的應(yīng)用[J].中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2017, 47(6): 87-94.

SUN Yang，JIANG Xiao-Dian，WANG Yi，et al.Application of gravity inversion from satellite altimeter data in the West Pacific[J].Periodical of Ocean University of China, 2017, 47(6): 87-94.