陳秋杰，沈云中，張興福，陳　武，許厚澤

1. 同濟(jì)大學(xué)測(cè)繪與地理信息學(xué)院，上海 200092； 2. 香港理工大學(xué)土地測(cè)量及地理資訊學(xué)系，香港； 3. 同濟(jì)大學(xué)空間信息與可持續(xù)發(fā)展應(yīng)用中心，上海 200092； 4. 廣東工業(yè)大學(xué)測(cè)繪工程系，廣東 廣州 510006； 5. 大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430077

?

基于GRACE衛(wèi)星數(shù)據(jù)的高精度全球靜態(tài)重力場(chǎng)模型

陳秋杰1,2,3，沈云中1，張興福4，陳武2，許厚澤5

1. 同濟(jì)大學(xué)測(cè)繪與地理信息學(xué)院，上海 200092； 2. 香港理工大學(xué)土地測(cè)量及地理資訊學(xué)系，香港； 3. 同濟(jì)大學(xué)空間信息與可持續(xù)發(fā)展應(yīng)用中心，上海 200092； 4. 廣東工業(yè)大學(xué)測(cè)繪工程系，廣東 廣州 510006； 5. 大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430077

摘要：應(yīng)用GRACE衛(wèi)星數(shù)據(jù)反演高精度靜態(tài)地球重力場(chǎng)是大地測(cè)量學(xué)界的熱點(diǎn)之一?？紤]到經(jīng)典動(dòng)力學(xué)法線性化誤差隨弧長(zhǎng)拉長(zhǎng)而迅速增長(zhǎng)，本文以GRACE衛(wèi)星軌道觀測(cè)值為初值的線性化方法，建立了應(yīng)用GRACE衛(wèi)星軌道和星間距離變率反演地球重力場(chǎng)的改進(jìn)動(dòng)力學(xué)法理論模型。利用2003年1月至2010年12月的GRACE衛(wèi)星姿態(tài)、軌道、星間距離變率和非保守力加速度等觀測(cè)數(shù)據(jù)，解算了一個(gè)180階次的無約束全球靜態(tài)重力場(chǎng)模型Tongji-Dyn01s和一個(gè)采用Kaula規(guī)則約束的全球重力場(chǎng)模型Tongji-Dyn01k。與國(guó)際不同機(jī)構(gòu)最新發(fā)布的純GRACE數(shù)據(jù)解算的重力場(chǎng)模型(包括AIUB-GRACE03S、GGM05S、ITSG-Grace2014k和Tongji-GRACE01)進(jìn)行比較，并利用DTU13海洋重力異常和GPS/水準(zhǔn)高程異常進(jìn)行外部檢核，結(jié)果表明，Tongji-Dyn01s與國(guó)際最新模型精度處于同一水平，然而Tongji-Dyn01k模型總體上更加靠近EIGEN6C2重力場(chǎng)模型。

關(guān)鍵詞：GRACE靜態(tài)重力場(chǎng)；動(dòng)力學(xué)法；線性化；正則化

地球重力場(chǎng)的分布源于地球質(zhì)量、密度分布以及質(zhì)量的重新分布(包括地球內(nèi)部運(yùn)動(dòng)、大氣運(yùn)動(dòng)、水循環(huán)、潮汐)。地球內(nèi)部物質(zhì)的質(zhì)量分布可劃分為穩(wěn)態(tài)分布與非穩(wěn)態(tài)分布，穩(wěn)態(tài)分布決定了靜態(tài)地球重力場(chǎng)，非穩(wěn)態(tài)的質(zhì)量重分布決定了時(shí)變地球重力場(chǎng)的特性。CHAMP(Challenging Mini-satellite Payload)[1]、GRACE(Gravity Recovery and Climate Experiment)[2]和GOCE(Gravity Field and Steady State Ocean Circulation Explorer)[3]是當(dāng)今衛(wèi)星重力探測(cè)的3類典型重力衛(wèi)星。其中，CHAMP主要探測(cè)低階地球重力場(chǎng)[1]；GRACE的最大貢獻(xiàn)在于探測(cè)非穩(wěn)態(tài)地球重力場(chǎng)并用于全球氣候時(shí)變監(jiān)測(cè)[4]；GOCE的主要目標(biāo)是探測(cè)高階靜態(tài)地球重力場(chǎng)信號(hào)[5]，尤其是100階后的重力場(chǎng)信息。對(duì)于探測(cè)前100階的重力場(chǎng)信號(hào)，GRACE衛(wèi)星的優(yōu)勢(shì)仍然明顯大于GOCE衛(wèi)星。盡管GARCE衛(wèi)星的預(yù)期壽命為5 a，但該衛(wèi)星自2002年3月發(fā)射后，至今仍在運(yùn)行，已經(jīng)提供了約12 a的觀測(cè)數(shù)據(jù)，因此可利用GRACE衛(wèi)星的長(zhǎng)時(shí)間數(shù)據(jù)序列反演高精度的靜態(tài)地球重力場(chǎng)。目前國(guó)際上有多個(gè)研究機(jī)構(gòu)利用GRACE衛(wèi)星長(zhǎng)時(shí)間觀測(cè)序列反演靜態(tài)地球重力場(chǎng)模型。其中，用純GRACE數(shù)據(jù)解算的最新無約束重力場(chǎng)模型包括GGM05S[6]、AIUB-GRACE03S[7]、ITSG-Grace2014s[8]以及Tongji-GRACE01[9]，約束重力場(chǎng)模型僅僅只有ITSG-Grace2014k[8]。我國(guó)很多學(xué)者對(duì)應(yīng)用衛(wèi)星重力數(shù)據(jù)反演靜態(tài)地球重力場(chǎng)作出了卓有成效的工作與貢獻(xiàn)，取得了許多研究成果[9-18]。如文獻(xiàn)[18]利用經(jīng)典動(dòng)力學(xué)法恢復(fù)了100階次的GRACE靜態(tài)重力場(chǎng)模型。

動(dòng)力學(xué)法是衛(wèi)星重力反演的常用方法[19]，其核心思想是根據(jù)衛(wèi)星初始狀態(tài)參數(shù)、先驗(yàn)重力場(chǎng)模型以及其他力模型確定的參考軌道為初值，對(duì)衛(wèi)星重力反演的非線性觀測(cè)方程進(jìn)行線性化。動(dòng)力學(xué)法的線性化觀測(cè)方程的設(shè)計(jì)矩陣需通過變分方程數(shù)值積分，且需多次迭代解算。文獻(xiàn)[20]提出以衛(wèi)星軌道觀測(cè)值為初值的線性化方法，然而涉及四重積分計(jì)算，模型比較復(fù)雜。文獻(xiàn)[21]采用類似的思想改進(jìn)加速度法，并利用GRACE衛(wèi)星軌道數(shù)據(jù)反演重力場(chǎng)。本文將該線性化方法用于改進(jìn)動(dòng)力學(xué)法，并利用2003年1月至2010年12月共8 a的GRACE衛(wèi)星軌道、星間距離變率、姿態(tài)和非保守力加速度數(shù)據(jù)，解算一個(gè)180階次的Tongji-Dyn01s全球靜態(tài)重力場(chǎng)模型。為改善高階位系數(shù)精度，進(jìn)一步使用Kaula規(guī)則約束位系數(shù)，解算了一個(gè)180階次的靜態(tài)重力場(chǎng)模型Tongji-Dyn01k。利用現(xiàn)有的高精度重力場(chǎng)模型，GPS水準(zhǔn)和DTU13海洋重力異常數(shù)據(jù)對(duì)Tongji-Dyn01s和Tongji-Dyn01k靜態(tài)重力場(chǎng)模型的精度進(jìn)行了檢核，結(jié)果表明Tongji-Dyn01s模型與國(guó)際最新的純GRACE重力場(chǎng)模型精度相當(dāng)，而且Tongji-Dyn01k模型總體上優(yōu)于國(guó)際最新的純GRACE重力場(chǎng)模型。

1理論方法

(1)

(2)

式中，a(t″)為t″時(shí)刻衛(wèi)星單位質(zhì)量所受合力，可由相鄰歷元加速度的K階多項(xiàng)式內(nèi)插而得

(3)

式中，λj為多項(xiàng)式系數(shù)。將式(3)代入式(1)和式(2)，可將式(1)和式(2)離散化為

(4)

(5)

(6)

式中，保守力ag是衛(wèi)星軌道位置r(ti)和重力位系數(shù)u的函數(shù)；非保守力與衛(wèi)星姿態(tài)qi和星固系下非保守力加速度觀測(cè)值facc有關(guān)；C為星固系到慣性系間的轉(zhuǎn)換矩陣。由于加速度計(jì)觀測(cè)值為有偏觀測(cè)量，一般需要估計(jì)加速度計(jì)尺度S和偏差p。由于尺度參數(shù)與重力場(chǎng)位系數(shù)存在較強(qiáng)相關(guān)性，因此本文與GFZ一樣，只估計(jì)偏差參數(shù)p。將式(4)—式(6)中的衛(wèi)星軌道以幾何軌道rk與為初值展開

(7)

i=1,2,…,N

(8)

i=1,2,…,N

(9)

(10)

i=1,2,…,N

(11)

假設(shè)有M個(gè)弧段，對(duì)于第k個(gè)弧段，為聯(lián)合衛(wèi)星軌道與星間距離變率數(shù)據(jù)求解重力場(chǎng)，該弧段所有歷元的衛(wèi)星軌道和星間距離變率的觀測(cè)方程可用矩陣形式表示為

Akδxk+Bkvk=lkk=1,2,…,M

(12)

(13)

式中，Qk為軌道和星間距離變率觀測(cè)值方差-協(xié)方差陣，可根據(jù)軌道和星間距離變率觀測(cè)值先驗(yàn)精度確定(軌道精度約為2 cm;星間距離變率精度約為0.2 μm/s)[23]。消去局部參數(shù)γk后得到只與位系數(shù)有關(guān)的法方程，進(jìn)一步累加各個(gè)弧段與位系數(shù)有關(guān)的法方程并求得重力場(chǎng)參數(shù)，具體解算方法可參考文獻(xiàn)[21—23]。當(dāng)解算的重力位模型達(dá)到180階次時(shí)，其法方程嚴(yán)重病態(tài)。為改善法方程的病態(tài)性，本文采用Kaula約束的正則化方法進(jìn)行位系數(shù)解算，相關(guān)算法參見文獻(xiàn)[24]。

2靜態(tài)地球重力場(chǎng)模型反演

本文采用2003年1月至2010年12月共8a的GRACE衛(wèi)星A和B的觀測(cè)數(shù)據(jù)反演了180階次的Tongji-Dyn01s靜態(tài)地球重力場(chǎng)。所用的數(shù)據(jù)包括由JPL提供的K波段星間距離變率、衛(wèi)星姿態(tài)、非保守力加速度數(shù)據(jù)以及武漢大學(xué)GNSS研究中心提供的約化動(dòng)力學(xué)軌道[25]。幾何軌道直接由GNSS觀測(cè)值確定，不受先驗(yàn)重力場(chǎng)制約，更適合于重力場(chǎng)反演；而約化動(dòng)力學(xué)軌道無疑會(huì)受先驗(yàn)重力場(chǎng)模型影響，但若與距離變率數(shù)據(jù)一起組合求解，先驗(yàn)重力場(chǎng)的影響就大為減少[22]，因此本文反演重力場(chǎng)時(shí)以約化動(dòng)力學(xué)軌道作為觀測(cè)值，且在線性化時(shí)也以約化動(dòng)力學(xué)軌道為初值。原始觀測(cè)數(shù)據(jù)需進(jìn)行預(yù)處理，比如衛(wèi)星姿態(tài)和非保守力加速度數(shù)據(jù)存在數(shù)據(jù)間斷，尤其是衛(wèi)星姿態(tài)的間斷十分頻繁，需進(jìn)行內(nèi)插處理；距離變率數(shù)據(jù)為GRACE衛(wèi)星最重要的觀測(cè)量，對(duì)中長(zhǎng)波重力信號(hào)十分敏感，然而其存在數(shù)據(jù)間斷和粗差，本文直接忽略掉間斷的距離變率數(shù)據(jù)；約化動(dòng)力學(xué)軌道十分平滑，不存在數(shù)據(jù)間斷和粗差，可用于剔除星間距離變率數(shù)據(jù)粗差。此外，還需顧及日月等行星三體擾動(dòng)、海潮、海洋極潮、固體潮、固體極潮、大氣與海洋的非潮汐變化和相對(duì)論效應(yīng)等保守力以及非保守力攝動(dòng),具體參見文獻(xiàn)[23]。

圖1　不同弧長(zhǎng)所得重力場(chǎng)解相對(duì)于EIGEN6C2大地水準(zhǔn)面階誤差Fig.1　Geoid degree errors of solutions computed by using different arc lengths with respect to EIGEN6C2

圖2　不同正則化因子所對(duì)應(yīng)的MSE的跡(取對(duì)數(shù))Fig.2　Trace of MSE (log10 scale) from different regularized factors

3Tongji-Dyn01s與Tongji-Dyn01k靜態(tài)地球重力場(chǎng)的精度分析

為了檢核Tongji-Dyn01s和Tongji-Dyn01k模型的精度水平，本文將Tongji-Dyn01s、Tongji-Dyn01k與最新發(fā)布的GGM05S、AIUB-GRACE03S、ITSG-Grace2014k、Tongji-GRACE01等純GRACE重力場(chǎng)模型進(jìn)行比較分析。其中，GGM05S為美國(guó)德克薩斯大學(xué)空間中心(CSR)利用經(jīng)典動(dòng)力學(xué)法基于10a的GRACE數(shù)據(jù)所解算的180階次無約束模型；AIUB-GRACE03S為瑞士伯尼爾大學(xué)基于天體力學(xué)方法(其本質(zhì)是動(dòng)力學(xué)法，區(qū)別在于采用了不同的參數(shù)化方法，如采用開普勒軌道根數(shù)并引入經(jīng)驗(yàn)加速度參數(shù))利用6a的GRACE數(shù)據(jù)所解算的160階次無約束模型；ITSG-Grace2014k為德國(guó)波恩大學(xué)基于經(jīng)典短弧長(zhǎng)積分法利用約10.5a的GRACE數(shù)據(jù)所解算的200階次正則化約束模型；Tongji-GRACE01為同濟(jì)大學(xué)利用改進(jìn)的短弧長(zhǎng)積分法基于4aGRACE數(shù)據(jù)所解算的160階次無約束模型。GFZ推出的EIGEN6C2模型采用多種觀測(cè)數(shù)據(jù)求得，包括GRACE、GOCE、Lageos(LaserGeodynamicsSatellite)[29]、海洋測(cè)高、地面重力以及航空重力數(shù)據(jù)，具有非常高的精度，本文將其作為基準(zhǔn)模型評(píng)價(jià)不同GRACE重力場(chǎng)模型內(nèi)符合精度。

各模型相對(duì)于EIGEN6C2的大地水準(zhǔn)面階誤差如圖3所示。由圖3可見，Tongji-Dyn01s和Tongji-Dyn01k在前100階次比當(dāng)前其他所有的GRACE重力場(chǎng)模型更加靠近EIGEN6C2。尤其是在前60階，Tongji-Dyn01s和Tongji-Dyn01k的大地水準(zhǔn)面誤差比GGM05S、AIUB-GRACE03S、ITSG-Grace2014k和Tongji-GRACE01等模型約小一個(gè)量級(jí)。Tongji-Dyn01s在100階次后與GGM05S、AIUB-GRACE03S和Tongji-GRACE01等模型非常接近。Tongji-Dyn01k在101到160階次間大地水準(zhǔn)面誤差稍大于ITSG-Grace2014k。然而在161到180階次之間，Tongji-Dyn01k大地水準(zhǔn)面誤差小于ITSG-Grace2014k。

圖3　不同模型相對(duì)于EIGEN6C2大地水準(zhǔn)面階誤差Fig.3　The geoid degree errors of different models with respect to EIGEN6C2

為了進(jìn)一步比較Tongji-Dyn01s和Tongji-Dyn01k模型與其他GRACE重力場(chǎng)模型間的內(nèi)符合精度，本文利用這些模型截?cái)嗟讲煌A次計(jì)算全球重力異常信號(hào)。圖4給出了不同模型扣除EIGEN6C2系數(shù)后截?cái)嘀?80階次的重力異常信號(hào)。由于EIGEN6C2使用了GOCE梯度數(shù)據(jù)，高階次系數(shù)精度優(yōu)于純GRACE模型，因此圖4反映不同純GRACE模型高階位系數(shù)的誤差。圖4表明，無約束解GGM05S和Tongji-Dyn01s模型高階位系數(shù)誤差均大于正則化解ITSG-GRACE2014k和Tongji-Dyn01k，且高緯度地區(qū)的誤差明顯大于低緯度地區(qū)；其主要原因是衛(wèi)星重力反演模型是病態(tài)的，且GRACE軌道在兩極存在1°左右空白，導(dǎo)致高階次位系數(shù)和高緯度區(qū)域的誤差偏大；正則化解能夠抑制高階次位系數(shù)的誤差，因此正則化解的結(jié)果要明顯優(yōu)于無約束解。Tongji-Dyn01k的誤差小于所有模型，其180階次位系數(shù)與EIGEN6C2重力異常差值標(biāo)準(zhǔn)差僅為3.6mGal(1Gal=10-2m/s2)。為了分析Tongji-Dyn01s和Tongji-Dyn01k的外部精度，筆者利用丹麥技術(shù)大學(xué)發(fā)布的DTU13產(chǎn)品(分辨率為2′×2′，精度約為2mGal)[30-31]中太平洋中部64 980個(gè)海洋重力異常數(shù)據(jù)，對(duì)上述不同模型在太平洋中部海洋重力異常進(jìn)行檢核。具體區(qū)域?yàn)?[17°S,23°S],[314°E,326°E])[32]，面積約為83.8萬km2。各模型分別截?cái)嘀?20、150和180階次，并由EIGEN6C2填補(bǔ)高階次系數(shù)(完全至1949階)，依次計(jì)算太平洋中部重力異常。表1給出了不同模型截?cái)嘀敛煌A次的重力異常差值統(tǒng)計(jì)結(jié)果?？梢钥闯觯瑹o論是截?cái)嗟?20還是150階次，不同模型所計(jì)算的重力異常與DTU13海洋重力異常差值的標(biāo)準(zhǔn)差均十分接近，可以說這些模型在前150階次的精度水平相當(dāng)。對(duì)于截?cái)嗟?80階次的結(jié)果，GGM05S精度最差，其次是Tongji-Dyn01s，Tongji-Dyn01k表現(xiàn)最好，其截?cái)嗟?80階所計(jì)算的太平洋中部海洋重力異常與DTU13重力異常之差標(biāo)準(zhǔn)差僅為4.1mGal，與截?cái)嗟?20和150階次所計(jì)算的標(biāo)準(zhǔn)差差別不大，然而其他模型截?cái)嗟?80階次與截?cái)嗟?20或150相比，其標(biāo)準(zhǔn)差均發(fā)生突變，有力地說明了Tongji-Dyn01k模型180階次位系數(shù)的可靠程度。可以說，截?cái)嘀?80階次計(jì)算重力異常差值能夠有效區(qū)分GRACE高階重力位系數(shù)的可靠程度。

圖4　不同純GRACE模型相對(duì)于EIGEN6C2的全球重力異常(截?cái)嘀?80階次)Fig.4　The difference of global gravity anomalies between GRACE-only models and EIGEN6C2 (truncated to degree and order 180)

表1不同模型截?cái)嘀敛煌A次的重力異常差值統(tǒng)計(jì)結(jié)果

Tab.1Statistics of the differences of gravity anomaly from DTU13 and that computed by 6 models truncated to different degrees

mGal

為了進(jìn)一步客觀評(píng)價(jià)Tongji-Dyn01s與Tongji- Dyn01k模型的精度水平，本文采用美國(guó)、墨西哥和江西3個(gè)地區(qū)的GPS/水準(zhǔn)網(wǎng)觀測(cè)數(shù)據(jù)進(jìn)行外部檢核。其中美國(guó)和墨西哥數(shù)據(jù)可從網(wǎng)站(http:∥www.ngs.noaa.gov/GEOID/GEOID12)下載，其水準(zhǔn)點(diǎn)數(shù)據(jù)本文分別選取8006個(gè)和536個(gè)，江西選取117個(gè)。各模型分別截?cái)嘀?0、120、150和180階，并采用EIGEN6C2填補(bǔ)高階位系數(shù)(完全至1949階)，依次計(jì)算美國(guó)、墨西哥和江西等地區(qū)的大地水準(zhǔn)面高程異常。統(tǒng)計(jì)計(jì)算值與實(shí)測(cè)值的最大值、最小值與標(biāo)準(zhǔn)差，結(jié)果如表2所示。表2中，對(duì)于截?cái)嘀料嗤A次(90和120)，從所計(jì)算的模型高程異常與實(shí)測(cè)結(jié)果之差的標(biāo)準(zhǔn)差來看，所有模型在美國(guó)、墨西哥和江西的表現(xiàn)相當(dāng)。對(duì)于截?cái)嗟?50階次，Tongji-Dyn01s的表現(xiàn)僅僅略遜于ITSG-Grace2014k，然而Tongji-Dyn01k優(yōu)于其他模型。尤其對(duì)于截?cái)嗟?80階次，無論是利用美國(guó)、墨西哥還是江西的數(shù)據(jù)進(jìn)行外部檢核，Tongji-Dyn01k模型精度均遠(yuǎn)優(yōu)于其他GRACE重力場(chǎng)模型。

4結(jié)論

本文以衛(wèi)星軌道觀測(cè)值為初值對(duì)經(jīng)典動(dòng)力學(xué)法衛(wèi)星重力反演的非線性觀測(cè)方程進(jìn)行線性化，建立了衛(wèi)星軌道和星間距離變率反演地球重力場(chǎng)的線性化觀測(cè)模型。與傳統(tǒng)動(dòng)力學(xué)法相比，本文所建立的動(dòng)力學(xué)法不需要解算變分方程與迭代計(jì)算，且不依賴于先驗(yàn)重力場(chǎng)模型。利用2003年1月至2010年12月GRACE衛(wèi)星的姿態(tài)、軌道、星間距離變率和非保守力加速度等觀測(cè)數(shù)據(jù)，解算了一個(gè)180階次的Tongji-Dyn01s全球靜態(tài)重力場(chǎng)模型；進(jìn)一步對(duì)重力位系數(shù)進(jìn)行了Kaula規(guī)則約束，通過正則化解法求得了Tongji-Dyn01k模型。利用多種方法與數(shù)據(jù)對(duì)Tongji-Dyn01s和Tongji-Dyn01k模型和國(guó)際不同機(jī)構(gòu)最新的純GRACE重力場(chǎng)模型AIUB-GRACE03S、ITSG-Grace2014k、GGM05S和Tongji-GARCE01的精度進(jìn)行了分析。

從模型的大地水準(zhǔn)面誤差看，Tongji-Dyn01s和Tongji-Dyn01k前100階位系數(shù)比國(guó)際不同機(jī)構(gòu)的最新GRACE重力場(chǎng)模型更加靠近EIGEN6C2。從全球重力異?？?，Tongji-Dyn01s僅次于ITSG-Grace2014k，而Tongji-Dyn01k相比于其他GRACE無約束模型更加靠近EIGEN6C2。DTU13太平洋中部海洋重力異常以及美國(guó)、墨西哥和江西區(qū)域GPS/水準(zhǔn)點(diǎn)的檢核結(jié)果，證實(shí)了Tongji-Dyn01s和Tongji-Dyn01k模型在前150階次與國(guó)際不同機(jī)構(gòu)最新的GRACE重力場(chǎng)模型精度處于同一水平，然而Tongji-Dyn01k在180階次的位系數(shù)精度優(yōu)于當(dāng)前所有GRACE數(shù)據(jù)反演的重力場(chǎng)模型，其所計(jì)算的太平洋海洋中部重力異常與DTU13數(shù)據(jù)之差標(biāo)準(zhǔn)差僅為4.1 mGal?？傮w來說，本文所建立的軌道與星間距離變率理論模型能夠有效地提取重力場(chǎng)信號(hào)，Kaula規(guī)則約束的正則化解法能夠顯著改善高階位系數(shù)解算精度。

致謝：感謝美國(guó)JPL和德國(guó)GFZ所提供的GRACE Level-1b觀測(cè)數(shù)據(jù)以及武漢大學(xué)GNSS研究中心趙齊樂教授所提供的約化動(dòng)力學(xué)軌道數(shù)據(jù)。

表2　各模型高程異常與實(shí)測(cè)GPS/水準(zhǔn)結(jié)果之差統(tǒng)計(jì)結(jié)果

參考文獻(xiàn)：

[1]REIGBER C, BALMINO G, SCHWINTZER P, et al. A High-quality Global Gravity Field Model from CHAMP GPS Tracking Data and Accelerometry (EIGEN-1S)[J]. Geophysical Research Letters, 2002, 29(14): 1-4.

[2]TAPLEY B D, BETTADPUR S, WATKINS M, et al. The Gravity Recovery and Climate Experiment: Mission Overview and Early Results[J]. Geophysical Research Letters, 2004, 31(9): 278-282.

[3]DRINKWATER M R, FLOBERGHAGEN R, HAAGMANS R, et al. GOCE: ESA’s First Earth Explorer Core Mission [C]∥BEUTLER G, DRINKWATER M R, RUMMEL R, et al. Proceedings of the Earth Gravity Field from Space — From Sensors to Earth Sciences. The Netherlands: Springer, 2006, 6-8.

[4]TAPLEY B D,BETTADPUR S,RIES J C, et al. GRACE Measurements of Mass Variability in the Earth System[J]. Science, 2004, 305(5683): 503-505.

[5]BOUMAN J, FIOROT S, Fuchs M, et al. GOCE Gravitational Gradients along the Orbit[J]. Journal of Geodesy, 2011, 85(11): 791-805.

[6]TAPLEY B D, FLECHTNER F, BETTADPUR S V, et al. The Status and Future Prospect for GRACE after the First Decade[C]∥AGU Fall Meeting. [S.l]: AGU, 2013: 9-12.

[8]MAYER-GüRR T,ZEHENTNER N,KLINGER B, et al. ITSG-Grace 2014: A New GRACE Gravity Field Release Computed in Graz[C]∥Oral Presentation at the GRACE Science Team Meeting. Potsdam: [s.n.], 2014, 29.

[9]SHEN Yunzhong,CHEN Qiujie,HSU H,et al. A Modified Short Arc Approach for Recovering Gravity Field Model[C]∥Oral Presentation at the GRACE Science Team Meeting. Austin: University of Texas, 2013: 22-25.

[10]沈云中. 應(yīng)用CHAMP衛(wèi)星星歷精化地球重力場(chǎng)模型的研究[D]. 武漢: 中國(guó)科學(xué)院測(cè)量與地球物理研究所, 2000.

SHEN Yunzhong. Study of Recovering Gravitational Potential Model from the Ephemerides of CHAMP[D]. Wuhan:The Institute of Geodesy and Geophysics, Chinese Academy of Sciences, 2000.

[11]王正濤. 衛(wèi)星跟蹤衛(wèi)星測(cè)量確定地球重力場(chǎng)的理論與方法[D]. 武漢: 武漢大學(xué), 2005.

WANG Zhengtao. Theory and Methodology of Earth Gravity Field Recovery by Satellite-to-Satellite Tracking Data[D]. Wuhan: Wuhan University, 2005.

[12]徐天河. 利用CHAMP衛(wèi)星軌道和加速度計(jì)數(shù)據(jù)推求地球重力場(chǎng)模型[D]. 鄭州: 信息工程大學(xué), 2004.

XU Tianhe. Gravity Field Recovery from CHAMP Orbit and Accelerometer Data[D]. Zhengzhou: The PLA Information Engineering University, 2004.

[13]肖云. 基于衛(wèi)星跟蹤衛(wèi)星數(shù)據(jù)恢復(fù)地球重力場(chǎng)的研究[D]. 鄭州: 信息工程大學(xué), 2006.

XIAO Yun. Analysis of Earth Gravity Field Recovery by Satellite-to-Satellite Tracking Data[D]. Zhengzhou: The PLA Information Engineering University, 2006.

[14]游為. 應(yīng)用低軌衛(wèi)星數(shù)據(jù)反演地球重力場(chǎng)模型的理論和方法[D]. 成都: 西南交通大學(xué), 2011.

YOU Wei. Theory and Methodology of Earth’s Gravitational Field Model Recovery by LEO Data[D]. Chengdu: Southwest Jiaotong University, 2011.

[15]張興福. 應(yīng)用低軌衛(wèi)星跟蹤數(shù)據(jù)反演地球重力場(chǎng)模型[D]. 上海: 同濟(jì)大學(xué), 2007.

ZHANG Xingfu. The Earth’s Field Model Recovery on the Basis of Satellite-to-Satellite Tracking Missions[D]. Shanghai: Tongji University, 2007.

[16]周旭華. 衛(wèi)星重力及其應(yīng)用研究[D]. 武漢: 中國(guó)科學(xué)院測(cè)量與地球物理研究所, 2005.

Zhou Xuhua. Study of Satellite Gravity and Its Application[D]. Wuhan: Institute of Geodesy and Geophysics, Chinese Academy of Sciences, 2005.

[17]蘇勇, 范東明, 游為. 利用GOCE衛(wèi)星數(shù)據(jù)確定全球重力場(chǎng)模型[J]. 物理學(xué)報(bào), 2014, 63(9): 099101.

SU Yong, FAN Dongming, YOU Wei. Gravity Field Model Calculated by Using the GOCE Data[J]. Acta Physica Sinica, 2014, 63(9): 099101.

[18]肖云, 夏哲仁, 王興濤. 用GRACE星間速度恢復(fù)地球重力場(chǎng)[J]. 測(cè)繪學(xué)報(bào), 2007, 36(1): 19-25.

XIAO Yun, XIA Zheren, WANG Xingtao. Recovering the Earth Gravity Field from Inter-satellite Range-rate of GRACE[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(1): 19-25.

[19]ZHU S, REIGBER C, K?NIG R. Integrated Adjustment of CHAMP, GRACE, and GPS Data[J]. Journal of Geodesy, 2004, 78(1-2): 103-108.

[20]XU Peiliang. Position and Velocity Perturbations for the Determination of Geopotential from Space Geodetic Measurements[J]. Celestial Mechanics and Dynamical Astronomy, 2008, 100(3): 231-249.

[21]陳秋杰, 沈云中, 張興福. 基于重力衛(wèi)星幾何軌道線性化的地球重力場(chǎng)反演方法[J]. 地球物理學(xué)報(bào), 2013, 56(7): 2238-2244.

Chen Qiujie, Shen Yunzhong, Zhang Xinfu. Linearization Method of Recovering Earth’s Gravity Field with Respect to Gravity Satellite’s Kinematic Orbits[J]. Chinese Journal of Geophysics, 2013, 56(7): 2238-2244.

[22]CHEN Qiujie, SHEN Yunzhong, ZHANG Xingfu, et al. Global Earth’s Gravity Field Solution with GRACE Orbit and Range Measurements Using Modified Short Arc Approach[J]. Acta Geodaetica et Geophysica, 2014, 50(2): 173-185.

[23]CHEN Qiujie, SHEN Yunzhong, ZHANG Xingfu, et al. Monthly Gravity Field Models Derived from GRACE Level 1B Data Using a Modified Short-arc Approach[J]. Journal of Geophysical Research, 2015, 120(3): 1804-1819.

[24]KUSCHE J, KLEES R. Regularization of Gravity Field Estimation from Satellite Gravity Gradients[J]. Journal of Geodesy, 2002, 76(6-7): 359-368.

[25]ZHAO Qile, GUO Jing, HU Zhigang, et al. GRACE Gravity Field Modeling with an Investigation on Correlation Between Nuisance Parameters and Gravity Field Coefficients[J]. Advances in Space Research, 2011, 47(10): 1833-1850.

[26]BETTADPUR S.Recommendation for A-Priori Bias & Scale Parameters for Level-1B ACC Data (Version 2)[EB/OL]. [2015-05-01]. http:∥podaac.jpl.nasa.gov/gravity/grace-documentation.

[27]DAHLE C, FLECHTNER F, GRUBER C, et al. GFZ GRACE Level-2 Processing Standards Document for Level-2 Product Release 0005[EB/OL]. [2015-05-01]. http:∥gfzpublic.gfz-potsdam.de/pubman/item/escidoc:61197:9/component/escidoc:65055/1202.

[28]F?RSTE C,BRUINSMA S,FLECHTNER F,et al.EIGEN-6C2-A New Combined Global Gravity Field Model Including GOCE Data up to Degree and Order 1949 of GFZ Potsdam and GRGS Toulouse[C]∥EGU General Assembly Conference. Vienna: [s.n.], 2013.

[29]DIEDRICH R, GENDT G. A Gravity Field Model from LAGEOS Based on Point Masses (POEM-L1)[C]∥Proceedings of the 6th International Symposium “Geodesy and Physics of the Earth”. Potsdam: GDR, 1989: 22-27.

[30]ANDERSEN O B, KNUDSEN P, KENYON S, et al. Global and Arctic Marine Gravity Field from Recent Satellite Altimetry (DTU13)[C]∥76th EAGE Conference and Exhibition 2014, Amsterdam, Netherlands: [s.n.], 2014: 16-19.

[31]ANDERSEN O B, JAIN M, KNUDSEN P. The Impact of Using Jason-1 and Cryosat-2 Geodetic Mission Altimetry for Gravity Field Modeling[J]. International Association of Geodesy Symposia, 2015: 1-6.

(責(zé)任編輯：陳品馨)

修回日期： 2015-12-31

First author： CHEN Qiujie(1987—), male, PhD candidate, majors in satellite gravimetry.

E-mail： chenqiujie2009@163.com

GRACE Data-based High Accuracy Global Static Earth’s Gravity Field Model

CHEN Qiujie1,2,3,SHEN Yunzhong1,ZHANG Xingfu4,CHEN Wu3,XU Houze5

1. College of Surveying and Geo-Informatics, Tongji University, Shanghai 200092, China; 2. Department of Land Surveying and Geo-Informatics, Hong Kong Polytechnic University, Hong Kong,China; 3. Center for Spatial Information Science and Sustainable Development, Shanghai 200092, China; 4. Departments of Surveying and Mapping, Guangdong University of Technology, Guangzhou 510006, China; 5. State Key Laboratory of Geodesy and Earth’s Dynamics, Institute of Geodesy and Geophysics, Wuhan 430077, China

Abstract：To recover the highly accurate static earth’s gravity field by using GRACE satellite data is one of the hot topics in geodesy. Since linearization errors of dynamic approach quickly increase when extending satellite arc length, we established a modified dynamic approach for processing GRACE orbit and range-rate measurements in this paper, which treated orbit observations of the twin GRACE satellites as approximate values for linearization. Using the GRACE data spanning the period Jan. 2003 to Dec. 2010, containing satellite attitudes, orbits, range-rate, and non-conservative forces, we developed two global static gravity field models. One is the unconstrained solution called Tongji-Dyn01s complete to degree and order 180; the other one is the Tongji-Dyn01k model computed by using Kaula constraint. The comparisons between our models and those latest GRACE-only models (including the AIUB-GRACE03, the GGM05S, the ITSG-Grace2014k and the Tongji-GRACE01) published by different international groups, and the external validations with marine gravity anomalies from DTU13 product and height anomalies from GPS/levelling data, were performed in this study. The results demonstrate that the Tongji-Dyn01s has the same accuracy level with those of the latest GRACE-only models, while the Tongji-Dyn01k model is closer to the EIGEN6C2 than the other GRACE-only models as a whole.

Key words：GRACE static gravity field; dynamic approach; linearization; regularization

第一作者簡(jiǎn)介：陳秋杰(1987—)，男，博士生，主要從事衛(wèi)星重力方面的研究。

收稿日期：2015-08-13

基金項(xiàng)目：國(guó)家973計(jì)劃(2012CB957703)；國(guó)家自然科學(xué)基金(41474017；41274035)；大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重力實(shí)驗(yàn)室開放基金(SKLGED2014-1-3-E)；地理信息工程國(guó)家重點(diǎn)實(shí)驗(yàn)室開放基金(SKLGIE2014-M-1-2)

中圖分類號(hào)：P223

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1001-1595(2016)04-0396-08

Foundation support： The National Basic Research Program of China (973 Program) (No. 2012CB957703)；The National Natural Science Foundation of China (Nos. 41474017; 41274035)； State Key Laboratory of Geodesy and Earth's Dynamics (No. SKLGED2014-1-3-E)； State Key Laboratory of Geo-information Engineering (No. SKLGIE2014-M-1-2)

引文格式：陳秋杰，沈云中，張興福,等.基于GRACE衛(wèi)星數(shù)據(jù)的高精度全球靜態(tài)重力場(chǎng)模型[J].測(cè)繪學(xué)報(bào)，2016,45(4)：396-403. DOI:10.11947/j.AGCS.2016.20150422.

CHEN Qiujie,SHEN Yunzhong,ZHANG Xingfu,et al.GRACE Data-based High Accuracy Global Static Earth’s Gravity Field Model[J]. Acta Geodaetica et Cartographica Sinica,2016,45(4):396-403. DOI:10.11947/j.AGCS.2016.20150422.