CryoSat-2衛(wèi)星測(cè)高計(jì)劃及其應(yīng)用

魏鑫 李斐 張勝凱 郝衛(wèi)峰

(1武漢大學(xué)測(cè)繪學(xué)院,湖北武漢430079;2武漢大學(xué)中國(guó)南極測(cè)繪研究中心,湖北武漢430079)

0 引言

冰凍圈看似遠(yuǎn)離人們的日常生活,但其變化對(duì)人類的生活環(huán)境有著極其重要的影響。冰蓋和海冰的質(zhì)量變化會(huì)影響海平面變化、地球輻射平衡和海鹽循環(huán)[1]。近10年的衛(wèi)星測(cè)高數(shù)據(jù)加深了科學(xué)家對(duì)極地內(nèi)陸冰蓋物質(zhì)平衡的認(rèn)識(shí)[2-3],但南極大陸和格陵蘭島的沿海區(qū)域(占冰蓋質(zhì)量30%)[3-4]始終難以測(cè)定,海冰的質(zhì)量變化也尚待深入了解。這些問題妨礙了對(duì)氣候變化效應(yīng)以及大氣和海洋潛在反饋的進(jìn)一步推測(cè)。

CryoSat-2是歐空局(ESA)地球探測(cè)計(jì)劃的組成任務(wù)之一。本文介紹歐空局CryoSat-2的科學(xué)目標(biāo)、科學(xué)需求、任務(wù)概況、儀器載荷、數(shù)據(jù)產(chǎn)品以及數(shù)據(jù)的校準(zhǔn),最后探討了CryoSat-2在極地領(lǐng)域的應(yīng)用研究。

1 CryoSat-2概況

1.1 科學(xué)目標(biāo)

冰凍圈在地球的輻射平衡中扮演著很重要的角色,海冰的減少會(huì)導(dǎo)致北極上空的溫室效應(yīng)較地球其他區(qū)域強(qiáng)。同時(shí),冰蓋的變化會(huì)嚴(yán)重影響海平面高度。CryoSat-2任務(wù)的目的在于測(cè)定冰凍圈的變化趨勢(shì),其中最重要的就是驗(yàn)證全球變暖導(dǎo)致的北極海冰逐年減少以及明確南極和格陵蘭冰蓋對(duì)全球海平面上升的貢獻(xiàn)[5]。另外CryoSat-2在季節(jié)性海冰區(qū)域、冰帽和冰川上空進(jìn)行大量測(cè)量,這些測(cè)量將有效地補(bǔ)充近十年來(lái)其他衛(wèi)星對(duì)冰凍圈觀測(cè)數(shù)據(jù)的不足。

1.2 科學(xué)需求

殘余不確定度(residual uncertainty)是衛(wèi)星系統(tǒng)設(shè)計(jì)及衡量測(cè)量結(jié)果質(zhì)量的指標(biāo)。在分析科學(xué)需求時(shí),需要估計(jì)在3年任務(wù)完成時(shí)測(cè)量結(jié)果的殘余不確定度殘余不確定度由厚度自然變化和測(cè)量誤差組成,其中厚度自然變化由已有資料獲取[6],計(jì)算公式如下:

殘余不確定度通過(guò)使測(cè)量結(jié)果接近厚度的自然變化來(lái)確定。測(cè)量誤差一般不超過(guò)由厚度自然變化引起的殘余不確定度的10%,CryoSat-2系統(tǒng)就是根據(jù)這一首要需求來(lái)設(shè)計(jì)。這里近似地給出殘余不確定度的定量表達(dá):

表1 CryoSat-2科學(xué)和測(cè)量需求Table 1.Overview of the CryoSat-2 science and measurement requirements

1.3 CryoSat-2任務(wù)

CryoSat-2于2010年4月在哈薩克斯坦拜科努爾航天發(fā)射場(chǎng)成功發(fā)射,任務(wù)設(shè)計(jì)運(yùn)行期3年[7]。CryoSat-2的科學(xué)軌道在高緯覆蓋和最大數(shù)量軌道交叉點(diǎn)之間做了最好的優(yōu)化。衛(wèi)星的近極地非太陽(yáng)同步軌道平均高度717 km,傾角92°。衛(wèi)星運(yùn)行周期為369天并伴隨著30天的子循環(huán)(369天的周期由連續(xù)變化的30天重復(fù)模式構(gòu)成)。CryoSat-2衛(wèi)星結(jié)構(gòu)簡(jiǎn)單,衛(wèi)星平臺(tái)長(zhǎng)4.5m,寬2.3 m,高2.2 m,重量為720 kg(圖1)。衛(wèi)星底部的S波段的螺旋天線用于接收來(lái)自地面的遙控指令并向地面?zhèn)鬏斝l(wèi)星狀態(tài)和監(jiān)測(cè)信息,X波段天線用于傳送測(cè)量數(shù)據(jù),數(shù)據(jù)量為320 Gbits。

圖1 CryoSat-2衛(wèi)星的儀器載荷[1]Fig.1.The CryoSat-2 satellite viewed from beneath[1]

2 儀器載荷

2.1 SIRAL

SIRAL是CryoSat-2的首要載荷,它在傳統(tǒng)脈沖壓縮雷達(dá)高度計(jì)的基礎(chǔ)上結(jié)合了合成孔徑和干涉技術(shù)。該傳感器使得CryoSat-2能測(cè)量冰蓋邊緣的復(fù)雜地形,并系統(tǒng)地監(jiān)測(cè)北極海冰的干舷高度(浮冰超出海面的高度)變化。

新SIRAL與傳統(tǒng)高度計(jì)主要有三點(diǎn)不同[8]:第一,CryoSat-2的軌道首次覆蓋到南北緯88°;第二,相比先前的ESA任務(wù)約10 km的衛(wèi)星足跡,SIRAL使用改進(jìn)的合成孔徑雷達(dá)技術(shù)將衛(wèi)星地面足跡減少到約沿軌0.3 km和跨軌1.5 km;第三,SIRAL對(duì)地表點(diǎn)進(jìn)行多視處理以減少由雷達(dá)斑點(diǎn)引起的噪聲。

2.1.1 SAR模式

CryoSat-2的沿軌高分辨率得益于在SAR模式中使用了改進(jìn)的合成孔徑雷達(dá)數(shù)據(jù)處理技術(shù)[9]。這使得衛(wèi)星能有效地區(qū)分狹窄冰間水道中的高反向散射海面和低反向散射浮冰表面。處理來(lái)自海面和浮冰表面的不同回波可以得到干舷高度,進(jìn)而得到海冰厚度。

衛(wèi)星移動(dòng)產(chǎn)生的多普勒效應(yīng)會(huì)造成雷達(dá)回波能量中的輕微頻移。同時(shí)處理64個(gè)連續(xù)相關(guān)的雷達(dá)脈沖,利用得到頻移信息可以將脈沖照亮的地面足跡劃分成很多窄的跨軌條帶,每個(gè)條帶沿軌方向上寬約250 m,見圖2a。

圖2 SIRAL高度計(jì)新的操作模式.(a)SAR模式提供250 m的沿軌分辨率;(b)SARIn模式利用隔開1 m基線的天線獲得跨軌干涉數(shù)據(jù)[1]Fig.2.Schematic illustration of key aspects of the new operating modes of the SIRAL altimeter instrument.(a)SAR mode provides 250 m along-track resolution;(b)SAR in mode exploits two antennaswith a baseline separation of 1 m to obtain additional across-track interferometric data[1]

2.1.2 SARIn模式

傳統(tǒng)的脈沖壓縮高度計(jì)測(cè)量的是衛(wèi)星到地表的最近距離,但無(wú)法提供回波來(lái)源的精確信息。在海洋和海冰上空,地表最近點(diǎn)可以認(rèn)為就位于星下點(diǎn)處。然而在面對(duì)復(fù)雜地形(大坡度的冰蓋邊緣)情況時(shí),回波來(lái)源于偏離星下點(diǎn)的未知點(diǎn),因此測(cè)量的距離不是真正的地表高程,這就造成了傳統(tǒng)高度計(jì)無(wú)法測(cè)量格陵蘭和南極冰蓋邊緣的地形變化。

合成孔徑干涉(SARIn)模式旨在獲取多變地形的高程估計(jì)。為實(shí)現(xiàn)SARIn測(cè)量,衛(wèi)星增加一副天線和SAR接收鏈路來(lái)形成跨軌干涉(1 m的基線)。其中一副天線發(fā)射的脈沖回波可以同時(shí)被兩副天線接收,偏離星下點(diǎn)的回波在路徑長(zhǎng)度和相位上存在差別,聯(lián)合1 m的天線間隔距離可以得到基線和回波方向之間的角度,進(jìn)而對(duì)測(cè)距值進(jìn)行修正。圖2b顯示了SARIn模式在冰蓋地形上空的運(yùn)行情況。

2.1.3 低分辨率模式

高度計(jì)在相對(duì)平坦的極地內(nèi)陸冰蓋上空采取傳統(tǒng)脈沖壓縮模式,也稱低分辨率模式。

2.2 其他載荷

衛(wèi)星多普勒定軌定位系統(tǒng)(DORIS),接收器接收全球?qū)⒔?0個(gè)地面臺(tái)站發(fā)射的信號(hào),根據(jù)多普勒頻移計(jì)算出衛(wèi)星運(yùn)行速度,進(jìn)而精確定軌。除了定軌之外,DORIS還可以應(yīng)用到建立和維持地球參考框架,測(cè)定地球定向參數(shù),地殼變形監(jiān)測(cè)以及大氣探測(cè)和研究。

星體跟蹤定位器(Star trackers),用于精確測(cè)定兩個(gè)SIRAL天線的基線方向,直接安裝在SIRAL天線支撐結(jié)構(gòu)上。

激光回射器(LRR),用于衛(wèi)星軌道的精確路基測(cè)量。

3 CryoSat-2數(shù)據(jù)產(chǎn)品

CryoSat-2的首要產(chǎn)品為海冰和陸地冰通量的變化圖。SIRAL的首要輸出為描述表面反射的原始波形。有效載荷數(shù)據(jù)段(PDS)負(fù)責(zé)下列數(shù)據(jù)產(chǎn)品的收集和存儲(chǔ):

(1)0級(jí):原始遙測(cè)數(shù)據(jù)包,誤差已被剔除,包含了時(shí)間和遙測(cè)質(zhì)量信息;

(2)全碼率(FBR):與1b級(jí)相似,但數(shù)據(jù)未經(jīng)SAR和SARIn模式平均處理。預(yù)期數(shù)據(jù)量430 Gbit/d;

(3)1b級(jí):包含儀器回波波形,在SAR和SARIn模式下這些波形被平均處理,數(shù)據(jù)量較FBR級(jí)小很多;

(4)2級(jí):包含高程估計(jì)數(shù)據(jù);

(5)監(jiān)測(cè)數(shù)據(jù):系統(tǒng)地監(jiān)測(cè)SIRAL的運(yùn)行情況。

CryoSat-2提供的數(shù)據(jù)產(chǎn)品需要大量的校準(zhǔn)和改正。不同的數(shù)據(jù)級(jí)描述見表2。

表2 CryoSat-2數(shù)據(jù)產(chǎn)品Table 2.CryoSat-2 data products

3.1 0級(jí)-全碼率數(shù)據(jù)

ESA地面站為了得到CryoSat-2高程估計(jì),需要在不丟失地球物理信息的前提下對(duì)數(shù)據(jù)進(jìn)行大量壓縮。SIRAL試驗(yàn)性傳感器的角色決定著壓縮算法需要鑒于發(fā)射后和Cal/Val(校準(zhǔn)和驗(yàn)證)期間積累的經(jīng)驗(yàn)不斷的改進(jìn)。壓縮借助底級(jí)數(shù)據(jù),FBR是所需數(shù)據(jù)信息完整保留的最高數(shù)據(jù)級(jí)。

3.2 1b級(jí)-單相結(jié)合波形數(shù)據(jù)

1b級(jí)數(shù)據(jù)由衛(wèi)星地面軌跡點(diǎn)的回波組成。數(shù)據(jù)進(jìn)行了壓縮,雖然1b級(jí)是沒有地球物理近似值的最高數(shù)據(jù)級(jí),但影響數(shù)據(jù)質(zhì)量的地球物理改正已經(jīng)得以應(yīng)用,這些改正包括:(1)電離層,干/濕大氣改正;(2)海潮和固體潮影響;(3)電磁(EM)偏差改正;(4)逆氣壓改正。

3.3 2級(jí)-沿軌高程

2級(jí)數(shù)據(jù)由地面點(diǎn)的高程估計(jì)組成,采樣率為沿軌1 s。沿軌高程是被直接用于地球物理研究的最底級(jí)數(shù)據(jù),也是不需要詳盡了解雷達(dá)信號(hào)處理的最底級(jí)數(shù)據(jù)。

3.4 監(jiān)測(cè)數(shù)據(jù)

監(jiān)測(cè)數(shù)據(jù)提供關(guān)于載荷運(yùn)行情況的即時(shí)信息,包括一系列需要快速測(cè)定的日常參數(shù)。當(dāng)發(fā)生異常情況時(shí),它將提供數(shù)據(jù)異常信息。

3.5 數(shù)據(jù)分配

CryoSat-2數(shù)據(jù)無(wú)需實(shí)時(shí)傳輸,使用者可以通過(guò)描述一個(gè)時(shí)空窗口來(lái)獲取1b級(jí)和2級(jí)數(shù)據(jù)。ESA同時(shí)也提供數(shù)據(jù)存取和轉(zhuǎn)換的專業(yè)軟件,數(shù)據(jù)分配通過(guò)CD-ROM光盤和FTP傳輸。表3匯總了分配要求。數(shù)據(jù)等待時(shí)間是數(shù)據(jù)接收和分配之間的時(shí)間延遲。最小和最大持續(xù)時(shí)間表示在軌衛(wèi)星提供有效數(shù)據(jù)產(chǎn)品的時(shí)間間隔。

表3 CryoSat-2數(shù)據(jù)分配要求Table 3.CryoSat-2 data distribution requirements

4 SIRAL數(shù)據(jù)校準(zhǔn)

CryoSat-2校準(zhǔn)和驗(yàn)證的目的在于驗(yàn)證各模式的測(cè)量不確定度SIRAL高度計(jì)數(shù)據(jù)只有通過(guò)校準(zhǔn)和驗(yàn)證才能進(jìn)行科學(xué)應(yīng)用。

4.1 發(fā)射前校準(zhǔn)

發(fā)射前需要大量地面測(cè)量來(lái)驗(yàn)證CryoSat-2的主要誤差源。已經(jīng)開展的機(jī)載、冰上和永久觀測(cè)站測(cè)量活動(dòng)提供了大量有價(jià)值的機(jī)載模擬和實(shí)測(cè)數(shù)據(jù)集。

2002年5月歐空局與美國(guó)航空航天局(NASA)合作開展LaRA(激光雷達(dá)高度計(jì))計(jì)劃。P3飛行器同時(shí)搭載激光高度計(jì)和延遲多普勒處理(D2P)Ku-波段雷達(dá)高度計(jì)飛經(jīng)格陵蘭島和海冰區(qū)域上空,有效地模擬了CryoSat-2的SARIn和 SAR處理,首次獲得了冰面的激光雷達(dá)剖面圖。雷達(dá)波能夠滲透干雪到達(dá)含鹽的冰面[10],激光能識(shí)別干雪表面,LaRA計(jì)劃證明聯(lián)合雷達(dá)和激光能有效認(rèn)識(shí)積雪覆蓋對(duì)CryoSat-2冰厚度反演精度的影響。另一個(gè)飛行計(jì)劃由Danish科學(xué)團(tuán)隊(duì)實(shí)施,獲得了Fram海峽和格陵蘭北部海冰區(qū)域的激光地形剖面圖。類似的驗(yàn)證計(jì)劃在接下來(lái)的幾年將持續(xù)進(jìn)行。

4.2 運(yùn)行中校準(zhǔn)

CryoSat-2衛(wèi)星系統(tǒng)的輔助測(cè)量模式和相關(guān)設(shè)備可以滿足系統(tǒng)誤差的獨(dú)立測(cè)量需要,在某些情況下系統(tǒng)也支持持續(xù)測(cè)量以改正原始數(shù)據(jù)。當(dāng)數(shù)據(jù)下行時(shí),雷達(dá)內(nèi)部的校準(zhǔn)通道會(huì)將傳輸信號(hào)直接導(dǎo)入接收鏈路,所有的儀器時(shí)間和相位通道都可以通過(guò)這種稱為“內(nèi)部”校準(zhǔn)的途徑來(lái)測(cè)量?！皟?nèi)部”校準(zhǔn)無(wú)法測(cè)量雷達(dá)天線系統(tǒng),也不能得到衛(wèi)星的位置和姿態(tài)誤差,這些系統(tǒng)因素的獲取都需要借助由天線生成的外部信號(hào),也就是“外部”校準(zhǔn)。校準(zhǔn)數(shù)據(jù)一般用于改正原始數(shù)據(jù)。

5 CryoSat-2成果在南北極的應(yīng)用

5.1 南極大陸和格陵蘭島的冰面高程變化

南極大陸和格陵蘭冰蓋的表面高程變化揭示了冰蓋對(duì)于冰蓋動(dòng)力學(xué)以及表面質(zhì)量平衡變化的響應(yīng),因此冰蓋表面地形及其變化研究成為冰川學(xué)家關(guān)注的重點(diǎn)。Helm等[11]基于CryoSat-2衛(wèi)星2011年1月到2014年1月,3年的測(cè)高數(shù)據(jù)得到了南極大陸和格陵蘭島冰蓋新的數(shù)字高程模型(DEM),同時(shí)也得到了兩者的高程變化和體積變化估計(jì)。

首先要對(duì)1b級(jí)CryoSat-2數(shù)據(jù)進(jìn)行一些預(yù)處理,Helm團(tuán)隊(duì)使用波形濾波器識(shí)別由于內(nèi)部跟蹤器失鎖或波形沒有清晰的前緣產(chǎn)生的“壞波形”,這些波形來(lái)自復(fù)雜地形區(qū)域。然后使用第一脈沖峰值閾值跟蹤算法(TFMRA)來(lái)測(cè)距。處理過(guò)程因測(cè)量模式不同而異,上述處理改正了大氣折射、固態(tài)地球和極潮引起的延遲誤差。

預(yù)處理得到的數(shù)據(jù),根據(jù)不同的測(cè)量模式采用不同方法以得到DEM。LRM數(shù)據(jù)經(jīng)過(guò)四次迭代來(lái)處理傾斜改正,進(jìn)而得到DEM。SARIn數(shù)據(jù)采用了新的干涉處理算法,具體步驟如圖3。得到的DEM精度與之前的DEM精度相符,與NASA冰、云和陸地高程ICESat數(shù)據(jù)對(duì)比顯示80%的CryoSat-2DEM的不確定度要＜(3±15)m。Helm利用新DEM進(jìn)行傾斜改正,獲得了2011—2014年的南極大陸和格陵蘭島的高程變化,將其與ICESat數(shù)據(jù)2003—2009年的高程變化率的對(duì)比揭示出西南極冰蓋體積減少量增加了3倍,而東南極和毛德皇后地的冰蓋出現(xiàn)了變厚的趨勢(shì)。格陵蘭和南極大陸冰蓋整體體積變化為(-503 ±107)km3·a-1。

南極大陸冰蓋邊緣變化顯著,Galin等[12]利用4年的CryoSat-2干涉數(shù)據(jù)獲取了松島冰川這一典型邊緣區(qū)域的高程變化。與此同時(shí),McMillan等[13]也驗(yàn)證了CryoSat-2在南極大陸復(fù)雜地形區(qū)域的觀測(cè)能力,無(wú)論是數(shù)據(jù)采樣率和測(cè)距精度方面都優(yōu)于傳統(tǒng)測(cè)高衛(wèi)星。

圖3 得到DEM及高程變化的處理過(guò)程Fig.3.Diagram of the processing scheme used to derive the final DEM and elevation change

5.2 北極海冰厚度和體積估計(jì)

近30年的衛(wèi)星記錄顯示北極海冰范圍在不斷縮小,至2012年達(dá)到最低點(diǎn)。泛北極冰與海洋建模和同化系統(tǒng)(PIOMAS)研究指出海冰的范圍縮小同時(shí)伴隨著體積減少[8]。CryoSat-2的高分辨率數(shù)據(jù)可以反演得到北極2010年11月和2011年12月的海冰體積估計(jì),有效地證實(shí)了PIOMAS的研究結(jié)果。

要得到CryoSat-2的海冰體積估計(jì)首要要區(qū)分來(lái)自海面和浮冰表面的回波。浮冰表面的漫反射回波和狹窄冰間水道的鏡面反射回波的反射率存在差異(入射角不同),據(jù)此可以區(qū)分海面和冰面的高程測(cè)量值。表面類型識(shí)別通過(guò)檢查“脈沖峰值(PP)”和“堆棧的標(biāo)準(zhǔn)偏差(SSD)”兩個(gè)參數(shù)完成,SSD是一個(gè)新參數(shù),它使用SIRALs多視來(lái)為測(cè)量表面后向散射變化提供入射角度[14]。來(lái)自冰間水道的回波PP＞18和SSD＜4,來(lái)自浮冰的回波PP＜9和SSD＞4。處理過(guò)程中需要改正距離窗口中心和回波前緣點(diǎn)之間的距離差距。冰間水道回波的處理使用了高斯加指數(shù)模型,而對(duì)于浮冰回波,跟蹤點(diǎn)定位在上升達(dá)到回波中的第一個(gè)峰值振幅50%處。漫反射和鏡面反射回波之間的模型匹配差距會(huì)引起高程偏差,通過(guò)在季節(jié)性海冰區(qū)域比較有冰和無(wú)冰情況下的海面高程可以移除高程偏差。

干舷轉(zhuǎn)換成厚度使用了更精細(xì)的方法,依據(jù)海冰存在的年數(shù)(一年或多年生)來(lái)修改雪荷載和冰密度。對(duì)于一年冰,在將干舷轉(zhuǎn)換成厚度前,雪厚度需要乘以系數(shù)0.5,處理中使用到的一年生冰密度為916.7 kg·m-3,多年生冰密度為882.0 kg·m-3[15-16]。

CryoSat-2海冰厚度估計(jì)受制于不同來(lái)源的偏差,因此需要大量的獨(dú)立數(shù)據(jù)集來(lái)對(duì)其進(jìn)行驗(yàn)證。在此之前,對(duì)干舷平均處理可以減少雷達(dá)斑點(diǎn)誤差。CryoSat-2數(shù)據(jù)分別與“極地-5”[17]數(shù)據(jù)、三波弗特環(huán)流回轉(zhuǎn)計(jì)劃(BGEP)仰視聲納(ULS)得到的冰吃水深度數(shù)據(jù)和NASA冰橋計(jì)劃(OIB)數(shù)據(jù)比較,得到的厚度差距均不超過(guò)0.1 m,ULS顯示出最好的一致性。然而,來(lái)自ULS和OIB的數(shù)據(jù)都不能大規(guī)模驗(yàn)證CryoSat-2厚度反演中的系統(tǒng)偏差,未來(lái)的校準(zhǔn)目標(biāo)在于進(jìn)一步約束CryoSat-2與實(shí)地?cái)?shù)據(jù)之間差異。

CryoSat-2海冰厚度估計(jì)沿襲傳統(tǒng)雷達(dá)高度計(jì)數(shù)據(jù)的處理方法排除了無(wú)冰水域[10],得到的海冰厚度借助ICESat任務(wù)[18-19]獲取的覆蓋面積即得北極海冰體積估計(jì)。對(duì)于超出CryoSat-2緯度限制的格網(wǎng)單元值用相近的包含厚度數(shù)據(jù)的格網(wǎng)單元值替代,海冰密集度取自微波成像儀(SSM/I)。

CryoSat-2海冰體積估計(jì)與 PIOMAS和 ICESat[20-21]任務(wù)的海冰體積估計(jì)對(duì)比進(jìn)一步證實(shí)了北極海冰體積的長(zhǎng)期減少趨勢(shì)。其中,秋季和冬季海冰體積平均減少500 km·a-1,相當(dāng)于0.007 5 m·a-1的厚度減少,與船載記錄得到的峰值變薄率相近。

5.3 南極冰下湖測(cè)定

在覆蓋南極大陸的厚冰蓋下,存在很多和海洋沒有聯(lián)系的淡水湖,它們對(duì)于理解南極冰下水傳輸和冰蓋動(dòng)力學(xué)意義重大。鉆孔可以獲取部分信息,但這在極地極端條件很難實(shí)現(xiàn)。衛(wèi)星測(cè)高或許是唯一能夠大規(guī)模獲取冰下湖信息的方式,聯(lián)合Cryo-Sat-2和ICESat衛(wèi)星測(cè)高數(shù)據(jù)已經(jīng)在南極大陸冰面發(fā)現(xiàn)了一個(gè)由冰下湖排空形成的巨大盆狀凹地,并繪制出了冰下湖的3D圖[22-23]。

CryoSat-2能夠測(cè)定南極冰下湖變化,其攜帶的雷達(dá)高度計(jì)能適應(yīng)惡劣氣象條件和長(zhǎng)時(shí)間黑暗觀測(cè),提供南極大陸的持續(xù)觀測(cè)值。南極大陸冰蓋底部已發(fā)現(xiàn)的冰下湖多達(dá)400個(gè)左右,這些湖的排空會(huì)破壞冰川下的生態(tài)系統(tǒng),引起上層冰蓋加速滑落入海。

5.4 海洋測(cè)量

CryoSat-2的主要目標(biāo)是測(cè)量極地海冰的厚度以及監(jiān)測(cè)覆蓋格陵蘭島和南極洲的冰蓋的變化。但是其收集到的數(shù)據(jù)也被應(yīng)用到其他研究領(lǐng)域,利用SAR模式數(shù)據(jù)獲取高精度的海面高程[24-25]進(jìn)而得到海底地形圖就是其中之一。衛(wèi)星的雷達(dá)測(cè)高儀不僅能夠探測(cè)到冰層厚度的微小變化,也可以測(cè)量海平面高度[26-27]。引力的作用使得海面地形能模擬海底起伏。質(zhì)量更大的區(qū)域,如海底的山脈,會(huì)產(chǎn)生強(qiáng)大的引力,吸引更多的水并在海平面高度上產(chǎn)生一個(gè)微小的增加。因此,儀器測(cè)量海面高度的同時(shí)可繪制出洋底未經(jīng)勘測(cè)區(qū)域的地形。CryoSat-2的雷達(dá)高度計(jì)能在海面檢測(cè)到重力場(chǎng),從而揭示5—10 km尺度的海底特征,這是15年來(lái)首次有高度計(jì)在如此高分辨率下繪制全球海洋重力場(chǎng)。研究表明聯(lián)合CryoSat-2高精度測(cè)距值和三年以上的海洋測(cè)繪數(shù)據(jù)可以得到全球海底地形圖。

6 結(jié)語(yǔ)

傳統(tǒng)測(cè)高衛(wèi)星由于軌道傾角的限制,未覆蓋到兩極相當(dāng)于歐洲面積大小的區(qū)域,CryoSat-2的出現(xiàn)將未覆蓋區(qū)域面積縮小了約10倍。CryoSat-2的沿軌高分辨率能識(shí)別80%左右的海冰,其獨(dú)特的多視處理能得到冰蓋邊緣等復(fù)雜地形的信息。CryoSat-2的高緯覆蓋、高分辨率和多視處理很好地補(bǔ)充了傳統(tǒng)遙測(cè)手段和測(cè)高衛(wèi)星的不足,其他遙感數(shù)據(jù)也能校準(zhǔn)CryoSat-2并彌補(bǔ)CryoSat-2短暫運(yùn)行周期的不足,幫助其發(fā)揮更大的科研價(jià)值。聯(lián)合CryoSat-2的雷達(dá)波數(shù)據(jù)和ICESat的激光測(cè)高數(shù)據(jù)能有效地實(shí)現(xiàn)數(shù)據(jù)校準(zhǔn)和測(cè)量積雪厚度,CryoSat-2的海冰厚度測(cè)量數(shù)據(jù)已將ICESat的北極海冰體積記錄延伸到2012年[28-31]。

CryoSat-2成功發(fā)射至今已三年,其數(shù)據(jù)產(chǎn)品已經(jīng)應(yīng)用到海冰、冰蓋和海洋等測(cè)量領(lǐng)域。隨著全世界CryoSat-2科研團(tuán)隊(duì)的不斷壯大,研究不斷深入,CryoSat-2這顆專為冰凍圈設(shè)計(jì)的衛(wèi)星的應(yīng)用前景會(huì)不斷被拓展,其科研價(jià)值也會(huì)逐漸凸顯。

1 Drinkwater M R,Francis R,Ratier G,etal.The European Space Agency’searth explorermission CryoSat:measuring variability in the cryosphere.Annals of Glaciology,2004,39(1):313—320.

2 Douglas B C,Cheney R E,Miller L,et al.Greenland ice sheet:Is it growing or shrinking.Science,1990,248(4953):288—289.

3 Wingham D J,Ridout A J,Scharroo R,et al.Antarctic elevation change from 1992 to 1996.Science,1998,282(5388):456—458.

4 Thomas R,Csatho B,Davis C,et al.Mass balance of higher-elevation parts of the Greenland ice sheet.Journal of Geophysical Research:Atmospheres(1984—2012),2001,106(D24):33707—33716.

5 Drinkwater M,Rebhan H.CryoSat science report.Noordwijk:European Space Agency,2003:1272.

6 Wingham D.CRYOSAT-2:A Mission to Determine Fluctuations in the Earth’s Ice Fields//EGS-AGU-EUG Joint Assembly.Toronto:UCL Press,2002:1750—1752.

7 Francis C R.CryoSat-2 Mission and data description.European Space Agency:European Space and Technology Centre Document,2001:23—29.

8 Laxon SW,Giles K A,Ridout A L,et al.CryoSat-2 estimates of Arctic sea ice thickness and volume.Geophysical Research Letters,2013,40(4):732—737.

9 Raney R K.CryoSat-2 SAR-mode looks revisited.IEEE Geoscience and Remote Sensing Letters,2012,9(3):393—397.

10 Laxon S,Peacock N,Smith D.High interannual variability of sea ice thickness in the Arctic region.Nature,2003,425(6961):947—950.

11 Helm V,Humbert A,Miller H.Elevation and elevation change of Greenland and Antarctica derived from CryoSat-2.The Cryosphere,2014,8:1539—1559.

12 Galin N,Shepherd A,Wingham D,etal.First results of CryoSat-2 observations of Pine Island Glacier//AGU FallMeeting Abstracts.San Francisco:AGU Publications,2012:C54A-03.

13 McMillan M,Shepherd A,Ridout A,etal.Antarctic Ice Sheet topography from CryoSat-2 interferometricmode data//AGU FallMeeting Abstracts.San Francisco:AGU Publications,2012:C51A—0767.

14 Wingham D J,Francis CR,Baker S,et al.CryoSat-2:A mission to determine the fluctuations in Earth’s land and marine ice fields.Advances in Space Research,2006,37(4):841—871.

15 Kurtz N T,Farrell S L.Large-scale surveys of snow depth on Arctic sea ice from Operation IceBridge.Geophysical Research Letters,2011,38(20):L20505.

16 Alexandrov V,Sandven S,Wahlin J,etal.The relation between sea ice thickness and freeboard in the Arctic.The Cryosphere Discussions,2010,4(2):641—661.

17 Haas C,Hendricks S,Eicken H,et al.Synoptic airborne thickness surveys reveal state of Arctic sea ice cover.Geophysical Research Letters,2010,37(9):L09501.

18 Kwok R,Rothrock D A.Decline in Arctic sea ice thickness from submarine and ICESat records:1958—2008.Geophysical Research Letters,2009,36(15):L15501.

19 Schweiger A,Lindsay R,Zhang J L,et al.Uncertainty in modeled Arctic sea ice volume.Journal of Geophysical Research:Oceans(1978—2012),2011,116(C8):C00D06.

20 Xie H,Ackley SF,YiD,etal.Sea-ice thickness distribution of the Bellingshausen Sea from surfacemeasurements and ICESataltimetry.Deep Sea Research PartⅡ:Topical Studies in Oceanography,2011,58(9-10):1039—1051.

21 李建成,范春波,褚永海,等.ICESAT衛(wèi)星確定南極冰蓋高程模型研究.武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2008,33(3):226—228.

22 European Space Agency(ESA).CryoSatmaps largest-ever flood beneath Antarctica.[2013-11-21].http://www.esa.int/Our-Activities/Observing-the-Earth/CryoSat/CryoSat-maps-largest-ever-flood-beneath-Antarctica.

23 McMillan M,Corr H,Shepherd A,etal.Three-dimensionalmapping by CryoSat-2 of subglacial lake volume changes.Geophysical Research Letters,2013,40(16):4321—4327.

24 Labroue S,Boy F,Picot N,et al.First quality assessment of the CryoSat-2 altimetric system over ocean.Advances in Space Research,2012,50(8):1030—1045.

25 Dibarboure G,Renaudie C,Pujol M I,et al.A demonstration of the potential of CryoSat-2 to contribute to mesoscale observation.Advances in Space Research,2012,50(8):1046—1061.

26 Andersen O B,Stenseng L,Knudsen P,et al.Arctic marine gravity and bathymetry from 3 years of Cryosat-2 SAR altimetry(DTU13 Gravity)//AGU Fall Meeting Abstracts.San Francisco:AGU Publications,2013:G11B-0922.

27 European Space Agency(ESA).CryoSatgoes to sea.[2013-11-21].http://www.esa.int/Our-Activities/Observing-the-Earth/CryoSat/CryoSat-goes-to-sea.

28 楊帆,溫家洪.ICESat與ICESat-2應(yīng)用進(jìn)展與展望.極地研究,2011,23(2):138—148.

29 文漢江,程鵬飛.ICESAT/GLAS激光測(cè)高原理及其應(yīng)用.測(cè)繪科學(xué),2005,30(5):32—35.

30 鄂棟臣,徐瑩,張小紅.ICESat衛(wèi)星及其在南極Dome A地區(qū)的應(yīng)用.武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2007,32(12):1139—1142.

31 鄂棟臣,沈強(qiáng),徐瑩,等.基于ASTER立體數(shù)據(jù)和ICESat/GLAS測(cè)高數(shù)據(jù)融合高精度提取南極地區(qū)地形信息.中國(guó)科學(xué)(D輯):地球科學(xué),2009,39(3):351—359.