• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    聯(lián)合GRACE和ICESat數(shù)據(jù)分離南極冰川均衡調(diào)整(GIA)信號(hào)

    2016-11-24 00:44:21高春春陸洋史紅嶺張子占姜永濤
    地球物理學(xué)報(bào) 2016年11期
    關(guān)鍵詞:冰蓋南極變化

    高春春,陸洋,史紅嶺,張子占,姜永濤

    1 南陽(yáng)師范學(xué)院 環(huán)境科學(xué)與旅游學(xué)院,河南南陽(yáng) 4730612 中國(guó)科學(xué)院測(cè)量與地球物理研究所 大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,武漢 4300773 中國(guó)科學(xué)院大學(xué),北京 100049

    ?

    聯(lián)合GRACE和ICESat數(shù)據(jù)分離南極冰川均衡調(diào)整(GIA)信號(hào)

    高春春1,2,陸洋2,3*,史紅嶺2,張子占2,姜永濤1

    1 南陽(yáng)師范學(xué)院 環(huán)境科學(xué)與旅游學(xué)院,河南南陽(yáng) 4730612 中國(guó)科學(xué)院測(cè)量與地球物理研究所 大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,武漢 4300773 中國(guó)科學(xué)院大學(xué),北京 100049

    2002年發(fā)射的GRACE重力衛(wèi)星為南極冰蓋質(zhì)量平衡提供了一種新的測(cè)量方式,但由于南極GIA模型的不確定較大,進(jìn)而影響GRACE結(jié)果的可靠性.本文聯(lián)合2003—2009年的GRACE和ICESat等數(shù)據(jù)實(shí)現(xiàn)了南極GIA信號(hào)的分離,聯(lián)合方法所分離的GIA不依賴于不確定性很大的冰負(fù)荷等假設(shè)模型,而是直接基于衛(wèi)星觀測(cè)數(shù)據(jù)估算而來(lái)的,具有更大的可靠性.在分離過(guò)程中,本文提出了冰流速度加權(quán)改正法和GPS球諧擬合改正法對(duì)GIA結(jié)果進(jìn)行精化,同時(shí)引入了南極GPS觀測(cè)站的位移數(shù)據(jù)對(duì)分離的GIA進(jìn)行詳細(xì)的評(píng)估和驗(yàn)證,GPS驗(yàn)證表明經(jīng)過(guò)冰流速度加權(quán)和GPS球諧擬合雙改正后的GIA結(jié)果精度明顯得到提高.最后本文利用所分離的GIA對(duì)GRACE和ICESat結(jié)果進(jìn)行了改正,得到2003—2009年南極冰蓋質(zhì)量變化的趨勢(shì)為-66.7±54.5 Gt/a (GRACE)和-77.2±21.5 Gt/a (ICESat),相比采用其他的GIA模型,本文的GIA結(jié)果使GRACE和ICESat這兩種不同觀測(cè)技術(shù)得到的南極冰蓋質(zhì)量變化結(jié)果更加趨于一致.

    南極冰蓋;冰川均衡調(diào)整;GRACE;ICESat;質(zhì)量變化

    The fact that GRACE senses all sources of mass change at the Earth′s surface requires removing non-glacier changes by independent datasets or models.Mass change signals associated with GIA are poorly known in Antarctica due to very sparse geophysical and climatological data to constrain the glacial history and often dominate GRACE error budgets.In this study,we explore an approach to estimate GIA through the combination of GRACE satellite gravity and Ice,Cloud,and Land Elevation Satellite (ICESat) altimetry data sets.The GIA results are completely independent of any previous reconstruction of the Antarctic glacial history and derived from direct observations by satellite techniques.

    Our study developed the approach proposed by Gunter et al.(2014) for empirical estimation of present-day Antarctic GIA and ice mass change using a combination of satellite altimetry and gravimetry.An ICESat surface height change estimated from the crossover analysis approach was combined with three different GRACE solutions for the period between February of 2003 and October of 2009.Based on the previous method of the low-precipitation zone (LPZ) bias correction,weighted ice flow speed correction and GPS spherical fitting methods are created.When the GIA results computed by the combination data were compared to the GPS ground station displacements,the results showed our correction methods helped reduce the impact of the mm-level biases.

    The new GIA results show the presence of strong uplift in the two main ice shelves and the northern Transantarctic Mountains regions in West Antarctica (WA),as well as subsidence in most areas of East Antarctica (EA) except for Kalser WilhelmⅡ Land.The total GIA-related mass change estimates for the entire AIS ranged from 38.4 to 47.5 Gt/a depending on different GRACE solution used.Our estimates are considerably smaller than the values of 75.6 to 146.8 Gt/a obtained from G14,W12a and ICE5G.Over the period between February of 2003 and October of 2009,the corresponding ice mass change for AIS showed the average values of -66.7±54.5 Gt/a (GRACE) and -77.2±21.5 Gt/a (ICESat).Compared to other GIA models,our estimates can minimize the discrepancy between the values of ice mass change measured by the GRACE and ICESat.

    1 引言

    南極冰蓋作為地球上最大的冰體,如果全部融化,海平面將上升約60 m,這無(wú)疑會(huì)給人類帶來(lái)巨大的環(huán)境災(zāi)難,因此精確評(píng)估南極冰蓋質(zhì)量平衡及其相關(guān)的氣候響應(yīng)具有顯著的科學(xué)、社會(huì)和經(jīng)濟(jì)效益(Velicogna and Wahr,2006).但由于南極冰蓋地處偏遠(yuǎn)且地域廣袤,探測(cè)其質(zhì)量平衡實(shí)際上非常困難,近20年來(lái)隨著衛(wèi)星測(cè)高、衛(wèi)星重力等空間大地測(cè)量技術(shù)的迅速發(fā)展,為南極冰蓋質(zhì)量平衡的監(jiān)測(cè)提供了多種強(qiáng)有力的新型觀測(cè)手段.空間大地測(cè)量技術(shù)具有受氣候環(huán)境限制小、覆蓋面積大、監(jiān)測(cè)能力強(qiáng)、可連續(xù)大規(guī)模對(duì)冰蓋進(jìn)行監(jiān)測(cè)并獲取多源數(shù)據(jù)的技術(shù)優(yōu)勢(shì),同時(shí)也是一種高效經(jīng)濟(jì)的觀測(cè)手段和方法.其中2002年發(fā)射的重力衛(wèi)星GRACE(Gravity Recovery and Climate Experiment)可通過(guò)測(cè)量地球重力場(chǎng)隨時(shí)間的變化來(lái)獲取南極冰蓋質(zhì)量的時(shí)空變化,GRACE探測(cè)南極冰蓋質(zhì)量平衡已成為近年來(lái)冰蓋物質(zhì)平衡研究領(lǐng)域的熱點(diǎn)和重點(diǎn),并涌現(xiàn)出一大批重大的研究成果(如:Tapley et al.,2004;Velicogna and Wahr,2006;Jacob et al.,2012;King et al.,2012;Gardner et al.,2013;Shepherd et al.,2012).但這些研究結(jié)果彼此并不一致,甚至差異很大,最新的研究結(jié)果表明差異較大的原因在于南極GIA(Glacial Isostatic Adjustment)模型的不確定性太大(Shepherd et al.,2012;Velicogna and Wahr,2013;高春春等,2015).

    GIA主要是指黏彈地球?qū)δ┐巫钍⒈谝詠?lái)冰退期地表冰負(fù)荷和海水負(fù)荷變化的響應(yīng),GIA能夠引起地球內(nèi)部物質(zhì)的重新分布導(dǎo)致地球重力場(chǎng)的變化(Wang and Wu,2006;賈路路等,2011),因此為了估算冰蓋質(zhì)量變化,需要從GRACE重力變化中扣除GIA信號(hào).但由于目前對(duì)冰蓋歷史變化和地幔黏滯度結(jié)構(gòu)缺乏足夠的認(rèn)識(shí),現(xiàn)有的GIA模型不可避免地存在著很大誤差和不確定性(賈路路,2012).前期的研究結(jié)果表明,不同的GIA改正導(dǎo)致GRACE估算結(jié)果的差異可達(dá)80 Gt/a,量級(jí)甚至超過(guò)了冰蓋質(zhì)量變化本身(高春春等,2015),因此GIA模型本身的極大不確定性很大程度上降低了GRACE結(jié)果的可靠性,南極GIA信號(hào)的精確估算成為GRACE監(jiān)測(cè)南極冰蓋物質(zhì)平衡研究中的重點(diǎn)和難點(diǎn).很多研究已經(jīng)提出采用其他觀測(cè)數(shù)據(jù)對(duì)GIA信號(hào)進(jìn)行約束和分離,如GPS(Global Positioning System)數(shù)據(jù)(Wang et al.,2013;Wu et al.,2010),ICESat(Ice,Cloud,and land Elevation Satellite)數(shù)據(jù)(Riva et al.,2009;Gunter et al.,2014).就南極地區(qū)而言,GPS觀測(cè)站布設(shè)稀少很難形成密集的觀測(cè)網(wǎng),難以單獨(dú)準(zhǔn)確構(gòu)建GIA,而2003年發(fā)射的激光測(cè)高衛(wèi)星ICESat可直接觀測(cè)南極冰蓋表面高程變化,與 GRACE衛(wèi)星有8年的重疊觀測(cè)期(2003—2009年),因此聯(lián)合ICESat和GRACE數(shù)據(jù)分離南極GIA信號(hào)是目前最佳的選擇.

    聯(lián)合GRACE和ICESat數(shù)據(jù)分離南極GIA信號(hào)的想法最早由Wahr等(2000) 提出,他們利用5年的GRACE和ICESat模擬數(shù)據(jù)對(duì)GIA進(jìn)行了分離,指出兩種衛(wèi)星數(shù)據(jù)的聯(lián)合能夠有效減小GIA結(jié)果的不確定性(Wahr et al.,2000).隨后Velicogna和Wahr(2002)提出利用衛(wèi)星重力、衛(wèi)星測(cè)高和GPS數(shù)據(jù)在頻域中通過(guò)迭代計(jì)算實(shí)現(xiàn)南極 GIA和現(xiàn)今冰蓋質(zhì)量變化的分離,并做了數(shù)值模擬(Velicogna and Wahr,2002).GRACE和ICESat衛(wèi)星發(fā)射后,Riva等(2009)利用5年的GRACE和ICESat數(shù)據(jù),采用粒雪和基巖的密度差結(jié)合質(zhì)量守恒和體積守恒原理在空間域?qū)崿F(xiàn)了南極 GIA 和冰蓋質(zhì)量變化的分離,得到南極GIA信號(hào)在南極引起的質(zhì)量變化為100±67 Gt/a(Riva et al.,2009).Gunter等(2014)利用2003—2009年的GRACE和ICESat數(shù)據(jù),對(duì)Riva等的分離方法進(jìn)行了改進(jìn),采用不同數(shù)據(jù)中心公布的GRACE數(shù)據(jù)得到南極GIA信號(hào)引起的質(zhì)量變化為53~103 Gt/a (Gunter et al.,2014).

    本文選取2003—2009年的GRACE和ICESat數(shù)據(jù),并結(jié)合同一時(shí)期的南極粒雪密實(shí)化FDM(Firn Densification Model)模型和表面質(zhì)量平衡SMB (Surface Mass Balance)模型對(duì)南極GIA信號(hào)進(jìn)行了分離,在分離過(guò)程中提出冰流速度加權(quán)法和GPS球諧擬合改正法對(duì)分離的GIA信號(hào)進(jìn)行進(jìn)一步精化,同時(shí)對(duì)分離的GIA信號(hào)進(jìn)行詳細(xì)的精度評(píng)估和GPS驗(yàn)證,最后對(duì)南極GIA信號(hào)和冰蓋質(zhì)量變化及其時(shí)空分布特征進(jìn)行了討論和分析.

    2 數(shù)據(jù)

    2.1 ICESat數(shù)據(jù)

    ICESat衛(wèi)星于2003年1月13日在美國(guó)加利福尼亞范登堡空軍基地成功發(fā)射,是人類歷史上第一顆激光測(cè)高衛(wèi)星,該衛(wèi)星的主要任務(wù)之一就是測(cè)量冰蓋的高程變化,衛(wèi)星軌道傾角94°,覆蓋南極大部分地區(qū),衛(wèi)星攜帶的激光測(cè)高儀GLAS(Geoscience Laser Altimeter System)對(duì)冰蓋表面反射非常敏感,反射信號(hào)不會(huì)深入冰雪面,激光腳點(diǎn)小(約60 m),沿軌采樣密集(水平間距約170 m),測(cè)量精度可達(dá)±10 cm (Schutz et al.,2005),ICESat激光測(cè)高儀的工作任務(wù)期見(jiàn)表1.利用GLAS數(shù)據(jù)采用交叉點(diǎn)分析法可以獲取ICESat任務(wù)期間的南極冰蓋高程變化結(jié)果(見(jiàn)圖1)(史紅嶺,2010;史紅嶺 et al.,2011).ICESat所得南極冰蓋高程變化的空間分辨率約為100 km,時(shí)間跨度為2003年02月至2009年10月,高程變化負(fù)增長(zhǎng)的區(qū)域主要集中在Amundsen海岸的Pine Island/Thwaites冰川周圍、Queen Maud Land北部的Jutulstraumen等冰川、Wilkes Land北部的Totten冰川以及橫貫?zāi)蠘O山脈北部的Mercer冰流等地區(qū),而高程變化正增長(zhǎng)的區(qū)域主要分布在南極半島的Palmer Land以及西南極的Kamb冰流和東南極廣大內(nèi)陸等地區(qū).

    2.2 GRACE數(shù)據(jù)

    為了更好地跟ICESat數(shù)據(jù)進(jìn)行匹配,需要盡可能選擇工作任務(wù)期跟ICESat相一致的GRACE月重力場(chǎng)模型,因此根據(jù)表1所列的ICESat工作任務(wù)期將GRACE GSM RL05產(chǎn)品中2003-02和2003-03平均作為第1任務(wù)期,2003-09、2003-10和2003-11平均作為第2任務(wù)期,2004-02和2004-03平均作為第3任務(wù)期,2004-05和2004-06平均作為第4任務(wù)期,2004-10和2004-11作為第5任務(wù)期,2005-02和2005-03平均作為第6任務(wù)期,2005-05和2005-06平均作為第7任務(wù)期,2005-10和2005-11平均作為第8任務(wù)期,2006-02和2006-03平均作為第9任務(wù)期,2006-05和2006-06平均作為第10任務(wù)期,2006-10和2006-11平均作為第11任務(wù)期,2007-03和2007-04平均作為第12任務(wù)期,2007-10和2007-11平均作為第13任務(wù)期,2008-02和2008-03平均作為第14任務(wù)期,2008-10月作為第15任務(wù)期,2008-11和2008-12平均作為第16任務(wù)期,2009-03和2009-04平均作為第17任務(wù)期,2009-10作為第18任務(wù)期.這樣將GSM數(shù)據(jù)產(chǎn)品歸算到18個(gè)工作任務(wù)期后GRACE數(shù)據(jù)便與ICESat數(shù)據(jù)形成了良好的一致性,有助于后期的聯(lián)合處理.

    圖1 (a) ICESat數(shù)據(jù)估算的南極冰蓋表面高程變化(2003-02—2009-10);(b)對(duì)應(yīng)的誤差Fig.1 (a) The height change trend for AIS estimates by ICESat data over the period between February of 2003 and October of 2009;(b) Corresponding uncertainties

    序號(hào)激光器任務(wù)期起始時(shí)間結(jié)束時(shí)間工作天數(shù)1L12003-02-202003-03-29382L2A2003-09-252003-11-18543L2B2004-02-172004-03-20334L2C2004-05-182004-06-20345L3A2004-10-032004-11-08376L3B2005-02-172005-03-24367L3C2005-05-202005-06-22348L3D2005-10-212005-11-23349L3E2006-02-222006-03-273410L3F2006-05-242006-06-253211L3G2006-10-252006-11-273412L3H2007-03-122007-04-143413L3I2007-10-022007-11-043614L3J2008-02-172008-03-213415L3K2008-10-042008-10-181416L2D2008-11-252008-12-172317L2E2009-03-092009-04-113418L2F2009-09-302009-10-1112

    為了測(cè)試不同GRACE數(shù)據(jù)產(chǎn)品和后處理方法對(duì)GIA信號(hào)的影響情況,本文盡可能選擇差別較大的GRACE結(jié)果,共選取了三種不同的GRACE數(shù)據(jù)產(chǎn)品并對(duì)應(yīng)三種不同的后處理方法,分別為:CSR-60(Slepian法)(高春春,2015)、CSR-96(DDK5濾波)(Kusche et al.,2009)和GFZ(組合濾波法) (高春春等,2015).此外,由于GRACE對(duì)重力場(chǎng)的短波信息相對(duì)不敏感,而且球諧截?cái)嗪秃筇幚硪矔?huì)使重力場(chǎng)短波信號(hào)進(jìn)一步丟失,從而導(dǎo)致GRACE結(jié)果出現(xiàn)泄漏誤差,一方面,南極區(qū)域外的非冰蓋質(zhì)量變化信號(hào)會(huì)泄漏到南極冰蓋內(nèi),形成內(nèi)泄漏(leakage-in)誤差;另一方面,南極區(qū)域內(nèi)的冰蓋質(zhì)量變化信號(hào)會(huì)削弱或泄漏到區(qū)域外,造成外泄漏(leakage-out)誤差.內(nèi)泄漏誤差主要來(lái)源于冰蓋以外的陸地水、非潮汐的海洋及冰蓋附近冰帽等質(zhì)量變化信號(hào);陸地水文信號(hào)引起的泄漏,本文采用GLDAS(Global Land Data Assimilation System)模型(Rodell et al.,2004)進(jìn)行模擬改正;海洋環(huán)流引起的泄漏,本文采用OGCM(Ocean General Circulation Model)數(shù)據(jù)(Dobslaw et al.,2013)進(jìn)行估算扣除;對(duì)于海陸質(zhì)量交換引起的泄漏,采用GLDAS模型計(jì)算的陸地水質(zhì)量平衡將其平鋪到海洋上進(jìn)行估算扣除.而外泄漏誤差本文采用一個(gè)新的信號(hào)恢復(fù)法——自適應(yīng)邊界選取法(高春春,2015)進(jìn)行改正,以避免尺度因子恢復(fù)法(Velicogna and Wahr,2006)引起的信號(hào)變形,泄漏誤差改正的結(jié)果見(jiàn)表2.

    圖2 GRACE數(shù)據(jù)估算的南極冰蓋質(zhì)量變化(2003-02—2009-10)及其不確定性(a) 采用CSR-60數(shù)據(jù)和Slepian法;(b) 采用CSR-96數(shù)據(jù)和DDK5濾波法;(c) 采用GFZ數(shù)據(jù)和組合濾波法;(d—f) 圖a—c對(duì)應(yīng)的誤差.Fig.2 The mass change trends and corresponding uncertainties for AIS estimate from GRACE over the period between February of 2003 and October of 2009(a) CSR-60 solutions and the Slepian method;(b) CSR-96 solutions and the DDK5;(c) GFZ solutions and the combined filter method;(d—f) the errors for Fig.a—c.

    圖3 (a)FDM數(shù)據(jù)估算的南極冰蓋表面高程變化趨勢(shì)(2003-02—2009-10);(b)對(duì)應(yīng)的誤差Fig.3 Surface height change trend estimates from FDM for the period between February of 2003 and October of 2009;(b) corresponding uncertainties

    表2 泄漏誤差改正值(單位:Gt/a)

    GRACE估算的南極冰蓋質(zhì)量變化(見(jiàn)圖2a—2c)顯示冰蓋質(zhì)量增加和減小的區(qū)域跟ICESat所得的冰蓋高程上升和下降的區(qū)域基本一致,但在Queen Maud Land和橫貫?zāi)蠘O山脈北部區(qū)域,ICESat結(jié)果有明顯的負(fù)增長(zhǎng),而GRACE沒(méi)有體現(xiàn).Queen Maud Land北部區(qū)域的不一致可能跟ICESat的不確定性或GRACE的空間分辨率不足有關(guān),而橫貫?zāi)蠘O山脈北部的冰川區(qū)域,可能跟GIA信號(hào)在該區(qū)域?qū)RACE結(jié)果影響較大有關(guān).

    2.3 FDM數(shù)據(jù)

    ICESat所觀測(cè)到的南極冰蓋表面高程變化dH/dt包含了很多因素的影響(Ligtenberg et al.,2011):

    (1)

    vacc為雪積累率,vme為雪融化率,vfc粒雪壓實(shí)率,vice冰流速率,vb基巖運(yùn)動(dòng).其中粒雪壓實(shí)率vfc只影響ICESat所觀測(cè)的冰蓋高程變化,并不影響GRACE所觀測(cè)的冰蓋質(zhì)量變化,如果不加以考慮,會(huì)降低分離GIA信號(hào)的可靠性,本文采用南極粒雪密實(shí)化FDM模型(Ligtenberg et al.,2011)來(lái)消除這種影響,FDM數(shù)據(jù)的時(shí)間分辨率為2天,依然參照ICESat工作任務(wù)期(表1),將FDM的時(shí)間歸算到與ICESat一致的工作任務(wù)期,然后進(jìn)行時(shí)間序列分析,最后獲得FDM冰蓋表面高程變化趨勢(shì)(見(jiàn)圖3a).

    3 GIA信號(hào)分離方法

    (2)

    乘上密度轉(zhuǎn)換為地表質(zhì)量變化為

    (3)

    (4)

    (5)

    其中

    (6)

    圖4 (a) 南極冰蓋表面密度;(b) 南極冰蓋表面質(zhì)量平衡SMB (2003—2009),白線內(nèi)區(qū)域?yàn)長(zhǎng)PZFig.4 (a) Surface density for AIS;(b) Surface mass balance (SMB) for AIS from 2003 to 2009,outline of the LPZ is white line

    圖5 (a) ICESat冰蓋高程變化與FDM之差;(b) 2δhFig.5 (a) Surface height change differences between the ICESat and FDM:firn;(b) 2δh

    圖6 聯(lián)合方法估算的GIA結(jié)果(a—i)以及GIA與GPS速率的比較,圓圈的半徑為,品紅色圓圈代表GPS>0,藍(lán)色圓圈代表Fig.6 Antarctic GIA vertical rates (a—i) computed from the combination approach and comparison of the vertical rates from the GPS stations and the GIA uplift rates,the radius of the circles is magenta circles represent GPS>0,blue circles represent GPS<0

    4 結(jié)果與精度評(píng)估

    4.1 初始結(jié)果

    4.2 LPZ偏差改正

    除了粗差,不同數(shù)據(jù)的聯(lián)合處理不可避免會(huì)造成所得的結(jié)果出現(xiàn)一定系統(tǒng)偏差,本文采用Gunter等(2014)提出的LPZ(Low-Precipitation Zone)偏差改正法來(lái)消除系統(tǒng)偏差的影響(Gunter et al.,2014),LPZ即為圖4b中白線所圈的區(qū)域,由于GIA信號(hào)在LPZ區(qū)域信號(hào)非常小,理論上假設(shè)為0,那么Ga、Gb、Gc中LPZ區(qū)域的信號(hào)即為偏差改正值(見(jiàn)表3),將Ga、Gb、Gc分別減去偏差改正值即可得到經(jīng)過(guò)LPZ偏差改正后的GIA結(jié)果:Ga-LPZ、Gb-LPZ、Gc-LPZ(見(jiàn)圖6d—6f).

    4.3 精度評(píng)價(jià)

    為了評(píng)估所計(jì)算GIA結(jié)果的可靠性,還需對(duì)其進(jìn)行精度評(píng)價(jià),其中ICESat和FDM的不確定性見(jiàn)圖1b和圖3b,SMB數(shù)據(jù)的不確定性為模型本身的10%(Rignot et al.,2008),ρrock的不確定性為100 kg·m-3(Wahr et al.,2000),ρα的不確定性為其本身的10%(Gunter et al.,2014),GRACE估算結(jié)果的誤差(見(jiàn)圖2d—2f)需要綜合考慮泄漏誤差、大氣誤差、混頻誤差、測(cè)量誤差等聯(lián)合估算(高春春等,2015;高春春,2015),各分項(xiàng)誤差獲取后,即可根據(jù)公式(5)按照誤差傳播定律計(jì)算Ga-LPZ、Gb-LPZ和Gc-LPZ的不確定性(見(jiàn)圖7a—7c).

    4.4 冰流速度加權(quán)改正法

    通過(guò)圖7a—7c可知,GIA信號(hào)誤差較大的區(qū)域主要集中在南極沿海地區(qū),因?yàn)檫@些區(qū)域是GRACE結(jié)果泄漏、混頻等誤差較為嚴(yán)重的區(qū)域,也是ICESat結(jié)果和FDM數(shù)據(jù)誤差較大的區(qū)域,而利用這些數(shù)據(jù)聯(lián)合估算的GIA信號(hào)在這些區(qū)域的誤差自然就相對(duì)較大.我們研究發(fā)現(xiàn)誤差較大的區(qū)域一般都集中在冰流速度較大的地區(qū)(見(jiàn)圖8a),因?yàn)楸魉俣仍酱蟮膮^(qū)域冰蓋高程變化越復(fù)雜,ICESat所監(jiān)測(cè)的結(jié)果不確定性就越大,冰流速度大的區(qū)域質(zhì)量變化信號(hào)也相對(duì)劇烈,GRACE結(jié)果中的泄漏誤差也就越嚴(yán)重,而且冰流速度大的區(qū)域一般分布在沿海地區(qū),氣候變化多樣,FDM和SMB數(shù)據(jù)誤差也相對(duì)較大,GRACE的大氣誤差也較為明顯,因此本文提出利用南極地區(qū)冰流速度分布數(shù)據(jù)(Rignot et al.,2011)對(duì)GIA結(jié)果進(jìn)行加權(quán)處理,即冰流速度加權(quán)改正法.

    具體做法:首先將冰流數(shù)據(jù)采用與ICESat等數(shù)據(jù)一樣的高斯濾波進(jìn)行平滑,這樣得到一個(gè)與其他數(shù)據(jù)空間分辨率一致的冰流速度vice flow,然后假定vice flow=100 m·a-1時(shí)權(quán)值為0,vice flow=0 m·a-1時(shí)權(quán)值為1,即

    (7)

    最后計(jì)算的權(quán)值wice flow見(jiàn)圖8b 那么將Ga-LPZ、Gb-LPZ、Gc-LPZ及其對(duì)應(yīng)的誤差都乘上wice flow,即可得到經(jīng)過(guò)加權(quán)改正后的GIA結(jié)果Ga-LPZ-w、Gb-LPZ-w、Gc-LPZ-w(見(jiàn)圖6g—6i)和相應(yīng)的誤差(見(jiàn)圖7d—7f).經(jīng)過(guò)加權(quán)改正之后GIA信號(hào)誤差(圖7d—7f)明顯小了許多,特別是在Amundsen海岸冰蓋質(zhì)量變化最劇烈的地區(qū)誤差得到了有效的控制.

    5 結(jié)果驗(yàn)證與GPS球諧擬合

    (8)

    其中

    (9)

    GIA信號(hào)和ICE5G、W12a、G14的GPS驗(yàn)證結(jié)果列于表3,由驗(yàn)證結(jié)果可知:①在東南極LPZ偏差改正法和冰流速度加權(quán)改正法能夠有效減小WRMS和RMS值,也就是有效提升GIA信號(hào)的可靠性和精度,而在西南極和南極半島只有Ga和Gb模型精度有所提高,而Gc精度則有所下降,就整個(gè)南極而言,兩種方法總體改善了GIA信號(hào)的精度;②本文分離的GIA在西南極和南極半島精度要明顯優(yōu)于ICE5G、W12a和G14這三個(gè)外圍GIA模型,但是在東南極精度要明顯差一些,在整個(gè)南極Ga_LPZ_w、Gb_LPZ_w和Gc_LPZ_w與外圍GIA模型精度相當(dāng).

    圖7 GIA結(jié)果的不確定性Fig.7 The uncertainties of GIA vertical rates

    圖8 (a) 南極冰流速度分布;(b) 利用南極冰流速度計(jì)算的權(quán)值Fig.8 (a) Antarctic ice velocity;(b) Weight computed based on Antarctic ice velocity

    從圖6可以看出:本文分離的GIA信號(hào)與GPS速率的差值在東南極基本都大于0,在南極半島都小于0,在西南極地區(qū)主要是大于0,也有少數(shù)站點(diǎn)小于0,這說(shuō)明GIA結(jié)果存在著區(qū)域性的系統(tǒng)偏差.為了消除這種區(qū)域性系統(tǒng)偏差,本文采用GPS擬合法對(duì)GIA結(jié)果進(jìn)行進(jìn)一步精化,針對(duì)南極地區(qū)的區(qū)域特點(diǎn),采用球諧函數(shù)進(jìn)行曲面擬合是相對(duì)較優(yōu)的選擇,根據(jù)GPS的站點(diǎn)數(shù)量采用4階球諧函數(shù)作為擬合模型進(jìn)行GPS擬合效果最佳,其中球諧擬合的數(shù)學(xué)模型為

    +Bn msinm λ]+ε.

    (10)

    圖9 GPS球諧擬合曲面Fig.9 GPS spherical fitting surfaces

    圖10 GPS球諧擬合改正后的GIA結(jié)果Fig.10 Antarctic GIA vertical rates corrected by GPS spherical fitting

    表3 LPZ偏差改正值和GPS驗(yàn)證結(jié)果

    注:WA+AP:西南極+南極半島,EA:東南極.

    6 GIA結(jié)果討論與應(yīng)用

    經(jīng)過(guò)聯(lián)合方法分離、LPZ偏差改正、冰流速度加權(quán)改正以及GPS球諧擬合改正一系列處理后的GIA結(jié)果Ga_LPZ_w_fit、Gb_LPZ_w_fit和Gc_LPZ_w_fit(見(jiàn)圖10a-10c)顯示:GIA信號(hào)上升速率較大的區(qū)域主要集中在西南極的Ronne冰架、橫貫?zāi)蠘O山脈北部和Ross冰架以及東南極Kalser Wilhelm Ⅱ Land等地區(qū).東南極的廣大內(nèi)陸地區(qū)主要呈現(xiàn)微弱的下降趨勢(shì),這種下降趨勢(shì)在W12a模型中有更為明顯的表現(xiàn)(見(jiàn)圖10e).Ga_LPZ_w_fit、Gb_LPZ_w_fit和Gc_LPZ_w_fit相比同樣采用GRACE和ICESat數(shù)據(jù)分離的G14最大的區(qū)別在于Amundsen海岸和東南極Kalser Wilhelm Ⅱ Land沿海區(qū)域,G14在Amundsen海岸表現(xiàn)出較強(qiáng)的隆升趨勢(shì),Ga_LPZ_w_fit、Gb_LPZ_w_fit和Gc_LPZ_w_fit并沒(méi)有體現(xiàn);Ga_LPZ_w_fit、Gb_LPZ_w_fit和Gc_LPZ_w_fit在Kalser Wilhelm Ⅱ Land沿海區(qū)域表現(xiàn)出較大的上升速率,而G14并沒(méi)有那么強(qiáng)烈.將Ga_LPZ_w_fit、Gb_LPZ_w_fit和Gc_LPZ_w_fit所得的高程變化與巖石密度相乘即可轉(zhuǎn)換為地表質(zhì)量變化并從GRACE估算結(jié)果中扣除,可以得到2003年02月至2009年10月GRACE觀測(cè)的南極冰蓋質(zhì)量變化趨勢(shì)(見(jiàn)圖11a—11c).另外還可利用Ga_LPZ_w_fit、Gb_LPZ_w_fit和Gc_LPZ_w_fit的高程變化作為對(duì)ICESat觀測(cè)的GIA改正,采用以下公式:

    (11)

    可以計(jì)算得到ICESat觀測(cè)的南極冰蓋質(zhì)量變化趨勢(shì)(見(jiàn)圖11d—11f),數(shù)值結(jié)果列于表4.

    從圖11可以看出2003年02月至2009年10月GRACE觀測(cè)到的南極冰蓋質(zhì)量消融最大的區(qū)域集中在Amundsen海岸和南極半島北部的冰川區(qū)域,而ICESat結(jié)果顯示在Amundsen海岸和橫貫?zāi)蠘O山脈北部的冰川區(qū)域,兩者最大的差別在于橫貫?zāi)蠘O山脈北部的Mercer冰流區(qū)域,這可能跟GIA信號(hào)在該區(qū)域較大有關(guān),因?yàn)镚IA對(duì)GRACE結(jié)果的影響相比對(duì)ICESat結(jié)果的影響要顯著的多.從表4可知:①未經(jīng)GPS擬合改正的GIA結(jié)果Ga_LPZ_w、Gb_LPZ_w和Gc_LPZ_w所計(jì)算的地表質(zhì)量變化分別為80.2±51.4 Gt/a、86.9±49.3 Gt/a和87.4±48.0 Gt/a,彼此相差很小,經(jīng)過(guò)GPS擬合后的GIA所計(jì)算的地表質(zhì)量變化分別為46.6±51.4 Gt/a、47.5±49.3 Gt/a和38.4±48.0 Gt/a,擬合改正后量級(jí)減小30~40 Gt/a;②利用Ga_LPZ_w、Gb_LPZ_w和Gc_LPZ_w分別進(jìn)行GIA改正后的GRACE結(jié)果彼此相近,三者的平均質(zhì)量變化趨勢(shì)為-107.4±54.5 Gt/a,與G14的-100±44 Gt/a(Gunter et al.,2014) 非常接近,而利用Ga_LPZ_w_fit、Gb_LPZ_w_fit和Gc_LPZ_w_fit進(jìn)行GIA改正的GRACE結(jié)果平均質(zhì)量變化趨勢(shì)為-66.7±54.5 Gt/a,與G14的-68.9±44 Gt/a更為接近;③不同的GIA模型對(duì)ICESat的結(jié)果影響很小,最大差異不超過(guò)10 Gt/a,而對(duì)GRACE的結(jié)果影響很大,最大差異可達(dá)75 Gt/a;④采用Ga_LPZ_w、Gb_LPZ_w和Gc_LPZ_w進(jìn)行GIA改正后的GRACE和ICESat結(jié)果彼此相差約25Gt/a,但是采用經(jīng)過(guò)GPS擬合改正后的Ga_LPZ_w_fit、Gb_LPZ_w_fit和Gc_LPZ_w_fit,GRACE和ICESat結(jié)果差異縮小為10 Gt/a,可見(jiàn)GPS擬合改正法有效降低了GIA的不確定性,使得不同源的衛(wèi)星數(shù)據(jù)估算的南極冰蓋質(zhì)量平衡趨于一致.

    圖11 GRACE(a—c)和ICESat(d—f)估算的南極冰蓋質(zhì)量變化趨勢(shì)(2003-02—2009-10)Fig.11 Antarctic ice sheet mass change trends estimated from GRACE (a—c) and ICESat (d—f) for the period between February of 2003 and October of 2009

    表4 GRACE和ICESat估算的南極冰蓋質(zhì)量變化(2003-02—2009-10)

    本文的GRACE結(jié)果平均為-66.7±54.5 Gt/a與Gunter等(2014)的結(jié)果-100±40 Gt/a有所差別,主要原因在于本文采用了GPS球諧擬合對(duì)GIA信號(hào)進(jìn)行了改正,而未改正前本文的結(jié)果是-107.4±54.5 Gt/a與Gunter解算的結(jié)果基本一致,與King等(2012)的GRACE結(jié)果-69±18 Gt/a非常接近(King et al.,2012),與賈路路等(2011)的-82±18 Gt/a(賈路路等,2011)、羅志才等(2012)的-80 Gt/a(羅志才等,2012)、Luthcke等(2013)的-81±26 Gt/a(Luthcke et al.,2013)、Velicogna和Wahr(2013)的-83±49 Gt/a(Velicogna and Wahr,2013)都相差不到20 Gt/a,GRACE結(jié)果彼此差異的原因主要在于選擇GIA改正信號(hào)的不同(高春春等,2015).本文的ICESat結(jié)果平均為-77.2±21.5 Gt/a,低于Groh等(2014)的-126±20 Gt/a(Groh et al.,2014),高于李斐等(2016)的-44±21 Gt/a(李斐等,2016),主要原因在于ICESat數(shù)據(jù)解算的方法不同,本文采用的是交叉點(diǎn)分析法,而Groh和李斐采用的是重復(fù)軌跡法,此外ICESat數(shù)據(jù)偏差改正方法的不同也會(huì)導(dǎo)致解算結(jié)果的差異(李斐等,2016).

    7 結(jié)論

    本文聯(lián)合GRACE衛(wèi)星重力和ICESat激光測(cè)高、南極FDM和SMB等數(shù)據(jù)分離了南極GIA信號(hào),所分離的GIA信號(hào)不依賴于不確定性很大的地球模型和冰負(fù)荷模型,是基于觀測(cè)數(shù)據(jù)直接估算的GIA信號(hào),更具有可靠性.分離過(guò)程中,為了有效降低GIA信號(hào)的不確定性,本文在前人LPZ偏差改正法的基礎(chǔ)上提出了冰流速度加權(quán)改正法和GPS球諧擬合改正法,實(shí)現(xiàn)了對(duì)GIA結(jié)果的進(jìn)一步精化,GPS驗(yàn)證結(jié)果顯示雙改正法能夠有效提高GIA結(jié)果的精度.最后本文對(duì)分離的GIA結(jié)果進(jìn)行了詳細(xì)的討論和應(yīng)用分析,利用聯(lián)合數(shù)據(jù)所估算的GIA對(duì)GRACE和ICESat進(jìn)行改正,得到南極冰蓋質(zhì)量變化的趨勢(shì)分別為-66.7±54.5 Gt/a(GRACE)和-77.2±21.5 Gt/a(ICESat),相比采用其他的GIA模型,本文所得的GIA使GRACE和ICESat這兩種不同測(cè)量方式所得的南極冰蓋質(zhì)量變化結(jié)果更加趨于一致.由于目前所能獲得的ICESat數(shù)據(jù)只有18個(gè)工作任務(wù)期,時(shí)間段相對(duì)較短,工作任務(wù)期較少,GRACE數(shù)據(jù)的空間分辨率也相對(duì)較低,南極GPS站點(diǎn)也比較少,聯(lián)合方法所分離的GIA結(jié)果不確定性還相對(duì)較大,但是隨著未來(lái)ICESat-2和GRACE Follow-on任務(wù)的實(shí)施,相信隨著更高質(zhì)量更高精度的時(shí)變重力和激光測(cè)高數(shù)據(jù)用于南極GIA信號(hào)的分離中,南極GIA結(jié)果的精度將會(huì)得到極大的改善.

    致謝 感謝匿名評(píng)審專家和編輯提供的寶貴建議和幫助,感謝CSR和GFZ提供GRACE數(shù)據(jù)產(chǎn)品,感謝荷蘭烏得勒支大學(xué)Lenaerts博士提供南極SMB數(shù)據(jù),感謝荷蘭烏得勒支海洋大氣研究所Ligtenberg博士提供南極FDM數(shù)據(jù),感謝荷蘭代爾夫特理工大學(xué)Gunter博士提供G14數(shù)據(jù),感謝英國(guó)杜倫大學(xué)Whitehouse博士提供W12a模型,感謝美國(guó)科羅拉多大學(xué)Geruo博士提供ICE5G模型,感謝英國(guó)紐卡斯?fàn)柎髮W(xué)Thomas博士提供南極GPS數(shù)據(jù),特別感謝中國(guó)科學(xué)院測(cè)量與地球物理研究所許厚澤院士、王勇研究員和汪漢勝研究員在本文完成中給予的指導(dǎo)和意見(jiàn).

    Dobslaw H,Flechtner F,Bergmann-Wolf I,et al.2013.Simulating high-frequency atmosphere-ocean mass variability for dealiasing of satellite gravity observations:AOD1B RL05.J.Geophys.Res.,118(7):3704-3711.

    Gao C C,Lu Y,Zhang Z Z,et al.2015.Ice sheet mass balance in Antarctica measured by GRACE and its uncertainty.Chinese J.Geophys.(in Chinese),58(3):780-792,doi:10.6038/cjg20150308.

    Gao C C.2015.Theory,method and application for detecting Antarctic ice sheet mass balance from time-variable gravity measurements [Ph.D.thesis] (in Chinese).Wuhan:Institute of Geodesy and Geophysics,Chinese Academy of Sciences.

    Gardner A S,Moholdt G,Cogley J G,et al.2013.A reconciled estimate of glacier contributions to sea level rise:2003 to 2009.Science,340(6134):852-857,doi:10.1126/science.1234532.

    Groh A,Ewert H,Rosenau R,et al.2014.Mass,volume and velocity of the Antarctic Ice sheet:present-day changes and error effects.Surv.Geophys.,35(6):1481-1505,doi:10.1007/s10712-014-9286-y.

    Gunter B C,Didova O,Riva R E M,et al.2014.Empirical estimation of present-day Antarctic glacial isostatic adjustment and ice mass change.The Cryosphere,8(2):743-760,doi:10.5194/tc-8-743-2014.

    Jacob T,Wahr J,Pfeffer W T,et al.2012.Recent contributions of glaciers and ice caps to sea level rise.Nature,482(7386):514-518,doi:10.1038/nature10847.

    Jia L L,Wang H S,Xiang L W,et al.2011.Effects of glacial isostatic adjustment on the estimate of ice mass balance over Antarctica and the uncertainties.Chinese J.Geophys.(in Chinese),54(6):1466-1477,doi:10.3969/j.issn.0001-5733.2011.06.006.

    Jia L L.2012.Effects of glacial isostatic adjustment on the estimate of present-day mass balance on the Earth′s surface [Ph.D.thesis] (in Chinese).Wuhan:Institute of Geodesy and Geophysics,Chinese Academy of Sciences.

    King M A,Bingham R J,Moore P,et al.2012.Lower satellite-gravimetry estimates of Antarctic sea-level contribution.Nature,491(7425):586-589,doi:10.1038/nature11621.

    Kusche J,Schmidt R,Petrovic S,et al.2009.Decorrelated GRACE time-variable gravity solutions by GFZ,and their validation using a hydrological model.J.Geod.,83(10):903-913,doi:10.1007/s00190-009-0308-3.

    Lenaerts J T M,van den Broeke M R,van de Berg W J,et al.2012.A new,high-resolution surface mass balance map of Antarctica (1979—2010) based on regional atmospheric climate modeling.Geophys.Res.Lett.,39(4):L04501,doi:10.1029/2011gl050713.

    Li F,Yuan L X,Zhang S K,et al.2016.Mass change of the Antarctic ice sheet derived from ICESat laser altimetry.Chinese J.Geophys.(in Chinese),59(1):93-100,doi:10.6038/cjg20160108.

    Ligtenberg S R M,Helsen M M,van den Broeke M R.2011.An improved semi-empirical model for the densification of Antarctic firn.The Cryosphere,5:809-819,doi:10.5194/tc-5-809-2011.

    Luo Z C,Li Q,Zhang K,et al.2012.Trend of mass change in the Antarctic ice sheet recovered from the GRACE temporal gravity field.Sci.China:Earth Sci.,55(1):76-82.

    Luthcke S B,Sabaka T J,Loomis B D,et al.2013.Antarctica,Greenland and Gulf of Alaska land-ice evolution from an iterated GRACE global mascon solution.J.Glaciol.,59(216):613-631.

    Rignot E,Bamber J L,Van Den Broeke M R,et al.2008.Recent Antarctic ice mass loss from radar interferometry and regional climate modelling.Nat.Geosci.,1(2):106-110,doi:10.1038/ngeo102.

    Rignot E,Mouginot J,Scheuchl B.2011.Ice flow of the Antarctic ice sheet.Science,333(6048):1427-1430.

    Riva R E M,Gunter B C,Urban T J,et al.2009.Glacial isostatic adjustment over Antarctica from combined ICESat and GRACE satellite data.Earth Planet.Sci.Lett.,288(3-4):516-523,doi:10.1016/j.epsl.2009.10.013.

    Rodell M,Houser P R,Jambor U,et al.2004.The global land data assimilation system.Bull.Amer.Meter.Soc.,85(3):381-394.

    Schutz B E,Zwally H J,Shuman C A,et al.2005.Overview of the ICESat Mission.Geophys.Res.Lett.,32:L21S01,doi:10.1029/2005GL024009.

    Shepherd A,Ivins E R,Geruo A,et al.2012.A reconciled estimate of ice-sheet mass balance.Science,338(6111):1183-1189,doi:10.1126/science.1128102.

    Shi H L.2010.Mass change detection in Antarctica Ice sheet using ICESat altimeter data [Ph.D.thesis].Wuhan:Institute of Geodesy and Geophysics,Chinese Academy of Sciences.

    Shi H L,Lu Y,Du Z L,et al.2011.Mass change detection in Antarctic ice sheet using ICESat block analysis techniques from 2003—2008.Chinese J.Geophys.(in Chinese),54(4):958-965,doi:10.3969/j.issn.0001-5733.2011.04.010.

    Tapley B D,Bettadpur S,Ries J C,et al.2004.GRACE measurements of mass variability in the Earth system.Science,305(5683):503-505,doi:10.1126/science.1099192.

    Thomas I D,King M A,Bentley M J,et al.2011.Widespread low rates of Antarctic glacial isostatic adjustment revealed by GPS observations.Geophys.Res.Lett.,38(22):L22302,doi:10.1029/2011gl049277.

    Velicogna I,Wahr J.2002.A method for separating Antarctic postglacial rebound and ice mass balance using future ICESat Geoscience Laser Altimeter System,Gravity Recovery and Climate Experiment,and GPS satellite data.J.Geophys.Res.,107(B10):ETG 20-1-ETG 20-11,doi:10.1029/2001JB000708.

    Velicogna I,Wahr J.2006.Measurements of time-variable gravity show mass loss in Antarctica.Science,311(5768):1754-1756.

    Velicogna I,Wahr J.2013.Time-variable gravity observations of ice sheet mass balance:Precision and limitations of the GRACE satellite data.Geophys.Res.Lett.,40(12):3055-3063,doi:10.1002/grl.50527.

    Wahr J,Wingham D,Bentley C.2000.A method of combining ICESat and GRACE satellite data to constrain Antarctic mass balance.J.Geophys.Res.,105(B7):16279-16294.

    Wang H S,Wu P.2006.Effects of lateral variations in lithospheric thickness and mantle viscosity on glacially induced relative sea levels and long wavelength gravity field in a spherical,self-gravitating Maxwell Earth.Earth Planet.Sci.Lett.,249(3-4):368-383.

    Wang H S,Jia L L,Steffen H,et al.2013.Increased water storage in North America and Scandinavia from GRACE gravity data.Nat.Geosci.,6:38-42,doi:10.1038/ngeo1652.

    Whitehouse P L,Bentley M J,Milne G A,et al.2012.A new glacial isostatic adjustment model for Antarctica:calibrated and tested using observations of relative sea-level change and present-day uplift rates.Geophys.J.Int.,190(3):1464-1482.

    Wu X P,Heflin M B,Schotman H,et al.2010.Simultaneous estimation of global present-day water transport and glacial isostatic adjustment.Nat.Geosci.,3(9):642-646,doi:10.1038/ngeo938.

    附中文參考文獻(xiàn)

    高春春,陸洋,張子占等.2015.GRACE重力衛(wèi)星探測(cè)南極冰蓋質(zhì)量平衡及其不確定性.地球物理學(xué)報(bào),58(3):780-792,doi:10.6038/cjg20150308.

    高春春.2015.時(shí)變重力探測(cè)南極冰蓋質(zhì)量平衡的理論、方法與應(yīng)用研究[博士論文].武漢:中國(guó)科學(xué)院測(cè)量與地球物理研究所.賈路路,汪漢勝,相龍偉等.2011.冰川均衡調(diào)整對(duì)南極冰質(zhì)量平衡監(jiān)測(cè)的影響及其不確定性.地球物理學(xué)報(bào),54(6):1466-1477,doi:10.3969/j.issn.0001-5733.2011.06.006.

    賈路路.2012.冰川均衡調(diào)整對(duì)現(xiàn)今地表質(zhì)量平衡監(jiān)測(cè)影響的研究[博士論文].武漢:中國(guó)科學(xué)院測(cè)量與地球物理研究所.李斐,袁樂(lè)先,張勝凱等.2016.利用ICESat數(shù)據(jù)解算南極冰蓋冰雪質(zhì)量變化.地球物理學(xué)報(bào),59(1):93-100,doi:10.6038/cjg20160108.羅志才,李瓊,張坤等.2012.利用GRACE時(shí)變重力場(chǎng)反演南極冰蓋的質(zhì)量變化趨勢(shì).中國(guó)科學(xué):地球科學(xué),42(10):1590-1596.

    史紅嶺.2010.ICESat在南極冰蓋質(zhì)量變化中的應(yīng)用研究[博士論文].武漢:中國(guó)科學(xué)院測(cè)量與地球物理研究所.

    史紅嶺,陸洋,杜宗亮等.2011.基于ICESat塊域分析法探測(cè)2003—2008年南極冰蓋質(zhì)量變化.地球物理學(xué)報(bào),54(4):958-965,doi:10.3969/j.issn.0001-5733.2011.04.010.

    (本文編輯 胡素芳)

    Combination of GRACE and ICESat data sets to estimate Antarctica Glacial Isostatic Adjustment (GIA)

    GAO Chun-Chun1,2,LU Yang2,3*,SHI Hong-Ling2,ZHANG Zi-Zhan2,JIANG Yong-Tao1

    1 Department of Environment Science and Tourism,Nanyang Normal University,Nanyang Henan 473061,China2 State Key Laboratory of Geodesy and Earth′s Dynamics,Institute of Geodesy and Geophysics, Chinese Academy of Sciences,Wuhan 430077,China3 University of Chinese Academy of Sciences,Beijing 100049,China

    The Antarctic Ice Sheet (AIS) is the largest single mass of ice on Earth,if fully melted,would contribute more than 60 meters of sea level rise.Accurate quantification of the mass balance of AIS is very important to improve our understanding and prediction of its response and contribution to sea level change.The Gravity Recovery and Climate Experiment (GRACE) mission has provided a new means to detect AIS mass change since its launch in 2002,but there are large uncertainties in applying GRACE to estimate AIS mass balance.Of these,the dominant error comes from the inaccurate knowledge of Glacial Isostatic Adjustment (GIA),which is the ongoing response of the solid Earth due to the changing ice-ocean load during the last ice age.

    Antarctic ice sheet;Glacial Isostatic Adjustment;GRACE;ICESat;Mass change

    高春春,陸洋,史紅嶺等.2016.聯(lián)合GRACE和ICESat數(shù)據(jù)分離南極冰川均衡調(diào)整(GIA)信號(hào).地球物理學(xué)報(bào),59(11):4007-4021,

    10.6038/cjg20161107.

    Gao C C,Lu Y,Shi H L,et al.2016.Combination of GRACE and ICESat data sets to estimate Antarctica Glacial Isostatic Adjustment (GIA).Chinese J.Geophys.(in Chinese),59(11):4007-4021,doi:10.6038/cjg20161107.

    國(guó)家重大科學(xué)研究計(jì)劃(2012CB957703),國(guó)家自然科學(xué)基金(41604009,41674085,41274025,41174064,41204013),大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室自主項(xiàng)目(SKLGED2013-2-2-Z) 和南陽(yáng)師范學(xué)院高層次人才科研啟動(dòng)費(fèi)聯(lián)合資助.

    高春春,男,1987年生,講師,主要從事空間大地測(cè)量技術(shù)及其在極地冰蓋監(jiān)測(cè)中的應(yīng)用等研究.E-mail:gaocc@asch.whigg.ac.cn

    *通訊作者 陸洋,男,1959年生,研究員,主要從事地球重力場(chǎng)精化、衛(wèi)星測(cè)高及南極冰蓋質(zhì)量平衡等研究.E-mail:luyang@whigg.ac.cn

    10.6038/cjg20161107

    P228

    2015-12-17,2016-03-30收修定稿

    猜你喜歡
    冰蓋南極變化
    軍事文摘(2024年6期)2024-04-30 03:13:59
    格陵蘭島的冰蓋悄悄融化
    參花(下)(2022年1期)2022-01-15 00:45:01
    我去南極能住哪兒
    來(lái)到南極要補(bǔ)課
    從9到3的變化
    南極大逃亡
    這五年的變化
    特別文摘(2018年3期)2018-08-08 11:19:42
    難忘的南極之旅
    長(zhǎng)距離輸水工程的冰期冰蓋數(shù)值模擬研究
    鳥的變化系列
    詩(shī)選刊(2015年6期)2015-10-26 09:47:11
    久久精品91无色码中文字幕| 啪啪无遮挡十八禁网站| 免费看十八禁软件| 可以在线观看的亚洲视频| 久久久国产成人免费| 国内精品久久久久久久电影| 色尼玛亚洲综合影院| 日韩欧美国产在线观看| 熟女人妻精品中文字幕| 精品久久久久久成人av| 九色国产91popny在线| 精品国产乱子伦一区二区三区| 黄色日韩在线| 搡老岳熟女国产| 亚洲,欧美精品.| 91字幕亚洲| 国产成人av激情在线播放| 亚洲一区二区三区不卡视频| 美女被艹到高潮喷水动态| 久久热在线av| 国产精品99久久久久久久久| 麻豆成人av在线观看| 一个人看的www免费观看视频| 精品国产亚洲在线| 人人妻人人澡欧美一区二区| 国产三级在线视频| 不卡一级毛片| 午夜免费成人在线视频| 18禁国产床啪视频网站| 国产精品自产拍在线观看55亚洲| 色精品久久人妻99蜜桃| 欧美成人性av电影在线观看| 欧美午夜高清在线| 国产精品久久久av美女十八| 亚洲午夜精品一区,二区,三区| 窝窝影院91人妻| 两个人看的免费小视频| 日本黄色视频三级网站网址| 久久中文看片网| 国产精品久久视频播放| 日本在线视频免费播放| 夜夜看夜夜爽夜夜摸| 免费av毛片视频| 欧美黑人巨大hd| 午夜福利成人在线免费观看| 国产精品98久久久久久宅男小说| 久久99热这里只有精品18| 色播亚洲综合网| 成人亚洲精品av一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品aⅴ一区二区三区四区| 两个人视频免费观看高清| 精品久久蜜臀av无| 一个人看的www免费观看视频| 国内揄拍国产精品人妻在线| 国产成人影院久久av| 日韩欧美国产在线观看| 一进一出抽搐gif免费好疼| 国产精品一区二区三区四区久久| 亚洲成人久久爱视频| 国产1区2区3区精品| 久久久久久人人人人人| 听说在线观看完整版免费高清| 国内精品久久久久久久电影| 成人国产综合亚洲| 亚洲激情在线av| 国模一区二区三区四区视频 | 亚洲精品一卡2卡三卡4卡5卡| 一a级毛片在线观看| 91在线精品国自产拍蜜月 | 一个人看视频在线观看www免费 | 午夜福利在线观看吧| svipshipincom国产片| 久久久久国产精品人妻aⅴ院| 一级a爱片免费观看的视频| а√天堂www在线а√下载| 精品一区二区三区视频在线观看免费| 91麻豆精品激情在线观看国产| 18禁裸乳无遮挡免费网站照片| 亚洲天堂国产精品一区在线| 亚洲人成网站在线播放欧美日韩| 1024香蕉在线观看| 国产一区二区三区在线臀色熟女| 亚洲人成网站在线播放欧美日韩| 国产麻豆成人av免费视频| 性欧美人与动物交配| 国产伦在线观看视频一区| 亚洲国产精品sss在线观看| 一夜夜www| 首页视频小说图片口味搜索| 国产精品1区2区在线观看.| 国产蜜桃级精品一区二区三区| 国产成人影院久久av| 淫妇啪啪啪对白视频| 两个人的视频大全免费| 亚洲九九香蕉| 高清在线国产一区| 精品国产超薄肉色丝袜足j| www.999成人在线观看| 成人性生交大片免费视频hd| 无限看片的www在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 母亲3免费完整高清在线观看| 国产精品久久电影中文字幕| 夜夜看夜夜爽夜夜摸| 免费高清视频大片| 久久这里只有精品19| 亚洲一区二区三区不卡视频| 国产精品免费一区二区三区在线| 亚洲人与动物交配视频| 国内精品一区二区在线观看| 99热只有精品国产| e午夜精品久久久久久久| 一个人看视频在线观看www免费 | 91在线精品国自产拍蜜月 | 一级黄色大片毛片| 亚洲精品色激情综合| 在线观看日韩欧美| 岛国在线观看网站| 国产高清视频在线观看网站| e午夜精品久久久久久久| www日本在线高清视频| 热99在线观看视频| 国产真实乱freesex| 亚洲av片天天在线观看| 少妇丰满av| 十八禁人妻一区二区| 中文字幕久久专区| 美女大奶头视频| 一个人观看的视频www高清免费观看 | 国产黄片美女视频| 哪里可以看免费的av片| 人妻久久中文字幕网| 一个人看的www免费观看视频| 久久天堂一区二区三区四区| 夜夜夜夜夜久久久久| 国产毛片a区久久久久| 免费av毛片视频| 国产精华一区二区三区| 亚洲无线观看免费| 欧美午夜高清在线| 精品熟女少妇八av免费久了| 国产精品亚洲美女久久久| 少妇丰满av| 黄片大片在线免费观看| 黑人操中国人逼视频| 欧美在线黄色| 午夜精品久久久久久毛片777| 一区二区三区国产精品乱码| www日本黄色视频网| 最近最新免费中文字幕在线| 欧美高清成人免费视频www| 身体一侧抽搐| 91老司机精品| or卡值多少钱| 18禁裸乳无遮挡免费网站照片| 精品久久久久久成人av| 黄色日韩在线| 久久亚洲精品不卡| 欧美一区二区国产精品久久精品| 久久九九热精品免费| 俺也久久电影网| 日本一二三区视频观看| 日韩高清综合在线| 亚洲国产精品成人综合色| 免费搜索国产男女视频| 免费人成视频x8x8入口观看| 一夜夜www| 18禁美女被吸乳视频| 两个人视频免费观看高清| 亚洲aⅴ乱码一区二区在线播放| 岛国在线免费视频观看| 日韩欧美免费精品| 午夜亚洲福利在线播放| 一夜夜www| 老鸭窝网址在线观看| 亚洲精品久久国产高清桃花| 又大又爽又粗| 美女cb高潮喷水在线观看 | 看免费av毛片| 成人欧美大片| 波多野结衣巨乳人妻| 国产亚洲精品av在线| 国产成+人综合+亚洲专区| 一个人看视频在线观看www免费 | 人人妻人人看人人澡| 午夜日韩欧美国产| 熟妇人妻久久中文字幕3abv| 亚洲欧美日韩东京热| 成年版毛片免费区| 免费av不卡在线播放| 精品国产亚洲在线| 1000部很黄的大片| 日韩人妻高清精品专区| 欧美日韩瑟瑟在线播放| 亚洲电影在线观看av| 在线观看免费视频日本深夜| 日韩中文字幕欧美一区二区| 不卡一级毛片| 国产乱人视频| av欧美777| 久久精品国产清高在天天线| 久久久久久久久久黄片| 国产黄a三级三级三级人| 久久久久久大精品| 51午夜福利影视在线观看| 精品久久久久久久末码| 免费观看的影片在线观看| 麻豆成人午夜福利视频| 99热6这里只有精品| 亚洲欧美精品综合久久99| 日本在线视频免费播放| 高潮久久久久久久久久久不卡| 亚洲中文日韩欧美视频| 99国产精品一区二区三区| 一个人看的www免费观看视频| 国产一级毛片七仙女欲春2| 男人舔奶头视频| 一a级毛片在线观看| 深夜精品福利| 视频区欧美日本亚洲| 一本一本综合久久| 淫妇啪啪啪对白视频| 精品一区二区三区视频在线观看免费| 色综合亚洲欧美另类图片| 十八禁人妻一区二区| 久久久精品大字幕| 两个人的视频大全免费| 国产精品爽爽va在线观看网站| xxxwww97欧美| 无限看片的www在线观看| 午夜福利18| 欧美成人性av电影在线观看| 日本三级黄在线观看| av天堂中文字幕网| 99精品欧美一区二区三区四区| 精品熟女少妇八av免费久了| 欧美乱码精品一区二区三区| 午夜福利欧美成人| 国产精品久久视频播放| 真实男女啪啪啪动态图| 午夜福利在线观看免费完整高清在 | 啪啪无遮挡十八禁网站| 亚洲精品乱码久久久v下载方式 | 国产一区二区三区视频了| 国产精品一区二区三区四区久久| 在线看三级毛片| 欧美乱色亚洲激情| 18禁黄网站禁片免费观看直播| 天天躁日日操中文字幕| 91麻豆av在线| 欧美日韩黄片免| 日韩欧美一区二区三区在线观看| 午夜精品久久久久久毛片777| 国产亚洲精品综合一区在线观看| 成人三级黄色视频| 性色avwww在线观看| 日韩欧美在线乱码| 亚洲自拍偷在线| 中国美女看黄片| 国产人伦9x9x在线观看| 成人av一区二区三区在线看| 国产精品精品国产色婷婷| 狂野欧美激情性xxxx| 好看av亚洲va欧美ⅴa在| 午夜影院日韩av| 一边摸一边抽搐一进一小说| 观看免费一级毛片| 男人舔奶头视频| 男女视频在线观看网站免费| 天堂网av新在线| 女人高潮潮喷娇喘18禁视频| 欧美激情在线99| 国产 一区 欧美 日韩| 午夜视频精品福利| 日日干狠狠操夜夜爽| 最新美女视频免费是黄的| 精品人妻1区二区| 欧美一区二区国产精品久久精品| 美女高潮的动态| av福利片在线观看| 亚洲av成人av| 嫁个100分男人电影在线观看| 国产69精品久久久久777片 | 欧美日韩乱码在线| 久久午夜综合久久蜜桃| 一级黄色大片毛片| 男女下面进入的视频免费午夜| 亚洲欧美日韩高清专用| 国产伦在线观看视频一区| 91麻豆精品激情在线观看国产| 很黄的视频免费| 欧美日韩瑟瑟在线播放| 国产av不卡久久| 欧美一级毛片孕妇| 中文亚洲av片在线观看爽| 久久久成人免费电影| 午夜精品在线福利| 精品国内亚洲2022精品成人| 身体一侧抽搐| 国内毛片毛片毛片毛片毛片| 真实男女啪啪啪动态图| 日韩欧美精品v在线| 99riav亚洲国产免费| 黑人操中国人逼视频| 国产91精品成人一区二区三区| 国产淫片久久久久久久久 | 国产精品久久久久久久电影 | 黄片小视频在线播放| 校园春色视频在线观看| 欧美在线黄色| 亚洲欧美精品综合久久99| 免费看a级黄色片| 在线观看日韩欧美| 国产爱豆传媒在线观看| 国产淫片久久久久久久久 | 99热这里只有是精品50| 国产精品久久久久久精品电影| 成人亚洲精品av一区二区| 亚洲人成网站高清观看| 免费av毛片视频| 精品久久久久久成人av| 久久这里只有精品中国| 亚洲avbb在线观看| avwww免费| 99久国产av精品| 婷婷精品国产亚洲av在线| 麻豆成人午夜福利视频| 亚洲一区二区三区色噜噜| 性色avwww在线观看| h日本视频在线播放| 岛国在线免费视频观看| 老熟妇乱子伦视频在线观看| 亚洲国产精品成人综合色| 天天一区二区日本电影三级| 成人午夜高清在线视频| 成人高潮视频无遮挡免费网站| 叶爱在线成人免费视频播放| 香蕉久久夜色| 免费大片18禁| 欧美乱码精品一区二区三区| 免费观看人在逋| 在线观看舔阴道视频| 欧美日韩中文字幕国产精品一区二区三区| 真人做人爱边吃奶动态| 一个人免费在线观看的高清视频| 看黄色毛片网站| 国产成人福利小说| 亚洲中文日韩欧美视频| 黄色视频,在线免费观看| 国产成人精品久久二区二区免费| 久久久精品欧美日韩精品| 日日摸夜夜添夜夜添小说| 亚洲av美国av| 看免费av毛片| 亚洲av日韩精品久久久久久密| 丁香六月欧美| 琪琪午夜伦伦电影理论片6080| 成人欧美大片| 国产精品爽爽va在线观看网站| 国产1区2区3区精品| 亚洲精品国产精品久久久不卡| 观看免费一级毛片| 五月玫瑰六月丁香| 一级毛片女人18水好多| 久久精品综合一区二区三区| 亚洲中文日韩欧美视频| 国产精品女同一区二区软件 | 欧美zozozo另类| 欧美激情久久久久久爽电影| 亚洲aⅴ乱码一区二区在线播放| 国产精品亚洲美女久久久| 国内毛片毛片毛片毛片毛片| 成人高潮视频无遮挡免费网站| а√天堂www在线а√下载| 天天一区二区日本电影三级| 夜夜爽天天搞| 在线a可以看的网站| 久久精品91无色码中文字幕| 熟女电影av网| 免费观看的影片在线观看| 免费av不卡在线播放| 免费看十八禁软件| 国产一区二区三区视频了| 1024手机看黄色片| 又粗又爽又猛毛片免费看| 精品一区二区三区视频在线 | 黄片小视频在线播放| 久久热在线av| 国产黄a三级三级三级人| 免费一级毛片在线播放高清视频| 综合色av麻豆| 少妇的逼水好多| or卡值多少钱| 搞女人的毛片| 天天躁日日操中文字幕| 99视频精品全部免费 在线 | 成人国产一区最新在线观看| 国产97色在线日韩免费| 五月伊人婷婷丁香| 日本在线视频免费播放| 国产精品免费一区二区三区在线| 亚洲 欧美一区二区三区| 天堂√8在线中文| 99热精品在线国产| 午夜福利免费观看在线| 久久精品国产综合久久久| 最近视频中文字幕2019在线8| 韩国av一区二区三区四区| 99国产综合亚洲精品| 不卡av一区二区三区| 亚洲欧美一区二区三区黑人| 在线观看66精品国产| 这个男人来自地球电影免费观看| 亚洲中文字幕一区二区三区有码在线看 | 好看av亚洲va欧美ⅴa在| 性色av乱码一区二区三区2| 国产毛片a区久久久久| 色精品久久人妻99蜜桃| 身体一侧抽搐| aaaaa片日本免费| 日韩成人在线观看一区二区三区| 亚洲国产精品久久男人天堂| 在线免费观看的www视频| 啦啦啦韩国在线观看视频| 90打野战视频偷拍视频| 亚洲自偷自拍图片 自拍| 18禁美女被吸乳视频| 欧美日韩精品网址| 高清毛片免费观看视频网站| 999精品在线视频| 午夜激情欧美在线| www日本黄色视频网| 在线播放国产精品三级| 天堂动漫精品| 久久久久久久久中文| 一个人免费在线观看电影 | 久久午夜亚洲精品久久| 日韩精品青青久久久久久| xxxwww97欧美| 国产主播在线观看一区二区| 夜夜躁狠狠躁天天躁| 国产亚洲欧美在线一区二区| 亚洲成av人片在线播放无| 精品国产乱码久久久久久男人| 国产精品一区二区三区四区免费观看 | 久久久色成人| 国产精品影院久久| 成人三级做爰电影| 亚洲av免费在线观看| 可以在线观看毛片的网站| 我要搜黄色片| 岛国在线观看网站| 午夜福利视频1000在线观看| 日本免费一区二区三区高清不卡| 免费看日本二区| 久久亚洲精品不卡| 一级毛片女人18水好多| 日韩欧美免费精品| 日本成人三级电影网站| 日韩有码中文字幕| 少妇人妻一区二区三区视频| 国产精品永久免费网站| 亚洲精品国产精品久久久不卡| 久久久久免费精品人妻一区二区| 免费人成视频x8x8入口观看| 色哟哟哟哟哟哟| 欧美色欧美亚洲另类二区| 男女下面进入的视频免费午夜| 一本精品99久久精品77| 女生性感内裤真人,穿戴方法视频| 国产伦一二天堂av在线观看| 一边摸一边抽搐一进一小说| 91字幕亚洲| 99久国产av精品| 色吧在线观看| 99久久99久久久精品蜜桃| 欧美成人性av电影在线观看| 欧美三级亚洲精品| 国产视频内射| 他把我摸到了高潮在线观看| 麻豆一二三区av精品| 91av网站免费观看| 变态另类成人亚洲欧美熟女| 一级作爱视频免费观看| 成人三级黄色视频| 九色国产91popny在线| 久久久久久大精品| 精品福利观看| 操出白浆在线播放| 国产精品香港三级国产av潘金莲| 中文字幕久久专区| 一夜夜www| 成人特级黄色片久久久久久久| 国产精品美女特级片免费视频播放器 | 麻豆成人午夜福利视频| 亚洲国产精品sss在线观看| 天堂影院成人在线观看| 一个人免费在线观看的高清视频| 国产精品一区二区精品视频观看| 一级作爱视频免费观看| 国产一区在线观看成人免费| 精品久久久久久久久久免费视频| 国产精品精品国产色婷婷| 一二三四在线观看免费中文在| 真人一进一出gif抽搐免费| 村上凉子中文字幕在线| 国产又黄又爽又无遮挡在线| 国产精品自产拍在线观看55亚洲| 麻豆成人av在线观看| 国产麻豆成人av免费视频| 国产不卡一卡二| 日韩大尺度精品在线看网址| 亚洲七黄色美女视频| 国产一级毛片七仙女欲春2| 国产视频内射| 午夜精品久久久久久毛片777| 亚洲美女黄片视频| 国产综合懂色| 亚洲九九香蕉| 动漫黄色视频在线观看| 桃色一区二区三区在线观看| 久久中文看片网| 黑人操中国人逼视频| 麻豆av在线久日| 大型黄色视频在线免费观看| 两人在一起打扑克的视频| 午夜免费观看网址| 亚洲无线观看免费| 免费观看的影片在线观看| 亚洲美女黄片视频| 91老司机精品| 91在线精品国自产拍蜜月 | 麻豆国产97在线/欧美| 亚洲最大成人中文| 国产淫片久久久久久久久 | 成年版毛片免费区| 国产麻豆成人av免费视频| 久久精品aⅴ一区二区三区四区| 国产1区2区3区精品| 母亲3免费完整高清在线观看| 日韩人妻高清精品专区| 久久伊人香网站| 综合色av麻豆| 亚洲人成网站在线播放欧美日韩| 国产97色在线日韩免费| 男女床上黄色一级片免费看| 国产精品乱码一区二三区的特点| 精品一区二区三区视频在线 | 亚洲精华国产精华精| 伊人久久大香线蕉亚洲五| 丰满人妻熟妇乱又伦精品不卡| 露出奶头的视频| 熟女电影av网| 免费看日本二区| 天堂√8在线中文| 激情在线观看视频在线高清| 国产美女午夜福利| 久久婷婷人人爽人人干人人爱| 看片在线看免费视频| 夜夜爽天天搞| 人人妻,人人澡人人爽秒播| 国产人伦9x9x在线观看| 精品国产乱码久久久久久男人| 国产99白浆流出| 精品日产1卡2卡| 亚洲中文av在线| 在线观看66精品国产| e午夜精品久久久久久久| 亚洲精品色激情综合| 欧美精品啪啪一区二区三区| 99精品久久久久人妻精品| 国产主播在线观看一区二区| 亚洲成av人片在线播放无| 美女cb高潮喷水在线观看 | 精品一区二区三区av网在线观看| 51午夜福利影视在线观看| 一进一出抽搐gif免费好疼| 久久伊人香网站| 俄罗斯特黄特色一大片| av在线蜜桃| 法律面前人人平等表现在哪些方面| 欧美性猛交黑人性爽| 亚洲成人免费电影在线观看| 美女高潮的动态| 长腿黑丝高跟| 欧美不卡视频在线免费观看| 免费在线观看亚洲国产| 国产精品,欧美在线| 搞女人的毛片| 日韩有码中文字幕| 香蕉av资源在线| 国产精品99久久久久久久久| 99精品久久久久人妻精品| 日韩欧美一区二区三区在线观看| 亚洲专区字幕在线| 欧美日韩福利视频一区二区| 香蕉av资源在线| 亚洲欧美日韩东京热| 一进一出抽搐动态| 国产精品一及| 搞女人的毛片| 一进一出抽搐动态| 97超视频在线观看视频| 又黄又爽又免费观看的视频| 在线观看免费午夜福利视频| 国产精品一及| 欧美激情久久久久久爽电影| 很黄的视频免费| www.999成人在线观看| 国产精品98久久久久久宅男小说| 美女高潮的动态| 亚洲国产欧美人成| 亚洲 欧美 日韩 在线 免费| 午夜福利成人在线免费观看|