gnss-r在海冰數(shù)據(jù)中的應(yīng)用

gnss-r在海冰數(shù)據(jù)中的應(yīng)用

0基于模型數(shù)據(jù)的渤海海冰遙感探測系統(tǒng)作為世界海洋和氣候系統(tǒng)的成員，海冰在海洋和氣候的相互作用上對世界的氣候變化有重大影響。海冰的生成、發(fā)展和消融過程,反映并影響了海洋與大氣的交換過程,進而影響海洋與大氣的環(huán)流模式。中國渤海和黃海北部歷年來出現(xiàn)了不同程度的凍結(jié),其冰情變化與全球氣候變化有密切關(guān)系。自1969年出現(xiàn)特大冰封以來,中國已全面開展了海冰監(jiān)測、預(yù)報和研究工作。在海冰常規(guī)觀測基礎(chǔ)上,進一步開展了衛(wèi)星和航空遙感技術(shù)的應(yīng)用研究。2010年,國家海洋局北海預(yù)報中心和國家海洋局第一海洋研究所針對渤海海冰和MODIS(中分辨率成像光譜輻射計)數(shù)據(jù)的特點,構(gòu)建了基于MODIS數(shù)據(jù)的渤海海冰遙感探測系統(tǒng)。2011年,中國海事局煙臺溢油應(yīng)急中心利用SAR(衛(wèi)星合成孔徑雷達)數(shù)據(jù)分析2010年渤海海域海冰冰緣線、海冰類型以及探測海冰運動過程,為有效應(yīng)對海冰災(zāi)害提供了科學(xué)參考。傳統(tǒng)的遙感手段如MODIS和SAR雖然具有較高的時空分辨率但是可用的信號有限,而且是單基工作模式,反射信號接收機復(fù)雜度高,遙感成本高。GNSS-R(GlobalNavigationSatelliteSystems-Reflectometry)遙感技術(shù)是20世紀(jì)90年代以來逐漸發(fā)展起來的GNSS的一個新型分支,是國內(nèi)外遙感探測和導(dǎo)航技術(shù)領(lǐng)域研究的熱點之一。國外采用GNSS-R技術(shù)進行海冰探測相關(guān)研究有如下幾個方面。2000年,Komjathy等利用GNSS反射信號的峰值功率獲得海冰表面有效介電常數(shù)的經(jīng)驗值(Permittivitybypeak-power)。2007年Belmonte擬合反射信號波形得到了海冰表面的粗糙度。2012年,Fabra等用歐洲空間局于2008—2009年在格林蘭島迪斯科海岸采集的數(shù)據(jù),采用雙極化相位測高的方法反演了海冰表面的絕對橢球高,取得了很好的測量精度。目前,國內(nèi)使用GNSS-R技術(shù)研究的方向有:海面?zhèn)蕊L(fēng)、海面測高、海洋鹽度探測、土壤濕度探測、海冰探測和移動目標(biāo)探測,而基于GNSS-R的海冰探測技術(shù)研究還處于起步階段,目前還沒有GNSS-R的海冰遙感實驗,因此我們采用歐空局網(wǎng)站上公布的海冰實驗數(shù)據(jù)進行理論方法的驗證和結(jié)果分析。1gnss—GNSS-R技術(shù)特點歐空局科學(xué)家Martin-Neira于1993年首次提出對GNSS反射信號進行利用的概念———PARIS(PassiveReflectometryandInterferometrySystem);1994年法國科學(xué)家Auber在進行機載飛行實驗時意外的發(fā)現(xiàn)了反射信號,證明反射信號是可以被接收并檢測的。GNSS接收機在接收導(dǎo)航衛(wèi)星直射信號的同時,也將接收反射面的反射信號。從電磁波傳播基本理論看來,該反射信號中攜帶反射面的特性信息,反射信號波形的變化、極化特征的變化,幅值、相位和頻率等參量的變化都直接反應(yīng)了反射面的物理特性,或者說直接與反射面相關(guān)。對反射信號的精確估計和接收處理,可以實現(xiàn)對反射面物理特性的估計與反演。目前GNSS-R在理論、技術(shù)和數(shù)據(jù)反演等方面逐步完善,GNSS衛(wèi)星可以在未來幾十年為人類提供精確、無償?shù)奶綔y信號。該信號具有長期的穩(wěn)定性,而且測量技術(shù)簡單,設(shè)備費用低,易于推廣。因此,這項技術(shù)為監(jiān)測海平面的變化提供了一種非常理想的手段。如圖1所示,如果只有一顆大型衛(wèi)星進行觀測,能夠檢測到的海面面積相對來說非常狹窄,所以存在時間分解能力不高的問題。如圖2所示,GNSS反射信號接收機同時可以接收多個GNSS衛(wèi)星的信號,所以相對于傳統(tǒng)的單個觀測衛(wèi)星,使用GNSS反射信號可以觀測到更大范圍的海面,并且能夠?qū)崿F(xiàn)高時間分辨率的海洋觀測。2電極關(guān)系的演化GNSS信號自導(dǎo)航衛(wèi)星發(fā)射后成為在空間傳輸?shù)碾姶挪?。在時變電磁場中,場矢量和場源既是空間位置的函數(shù),又是時間的函數(shù)。在正弦穩(wěn)態(tài)條件下,由場源所激勵的場矢量的各個分量仍是同頻率的正弦時間函數(shù)。時變電磁場中的任一坐標(biāo)分量隨時間做正弦變化時,其振幅和初始相位也都是空間坐標(biāo)的函數(shù)。以電場強度E為例,以一定的頻率ω隨時間t和空間r按正弦規(guī)律變化,可表示為由于電場強度E、磁場強度H和傳播方向K三者之間的關(guān)系是確定的,所以一般用電場強度E的矢量端點在空間任意固定點上隨時間變化所描述的軌跡來表示電磁波的極化。假設(shè)均勻平面波沿著z軸方向傳播,電場強度和磁場強度均在垂直于z軸的平面內(nèi),令電場強度E分解為兩個相互正交的分量Ex和Ey,其頻率和傳播方向均相同,E矢量端點的軌跡方程可以經(jīng)由三角運算獲得,當(dāng)滿足條件Ex0=Ey0=E0,φy-φx=±π/2時,矢量E端點的軌跡方程為這是半徑為E0的圓的方程,故而稱為圓極化。當(dāng)Ey滯后于Ex,π/2時,電場矢量的旋向和波的傳播方向滿足右手螺旋關(guān)系,稱為右旋圓極化(RHCP);反之稱為左旋圓極化(LHCP)。GNSS的導(dǎo)航信號(直射信號)是右旋圓極化的,當(dāng)直射信號照射到物體表面時極化方式會發(fā)生改變,部分轉(zhuǎn)變?yōu)樽笮龍A極化信號。研究反射信號的極化類型的變化可以提取出反射物體的物質(zhì)特性信息,這是GNSS-R技術(shù)的理論基礎(chǔ)。3數(shù)據(jù)收集3.1ceifr接收本文分析的海冰數(shù)據(jù)是歐洲空間局(ESA)于2009年1—5月在格林蘭島迪斯科海岸,采用ICE/IEEC-CSIC制造的GPS開環(huán)實時差分接收機(GPSOpenLoopDifferentialReal-TimeReceiver,GOLDRTR)采集。接收機原理如圖3所示。GOLD-RTR接收機的Link-1接口與射頻前端的天頂右旋圓極化(RHCP)天線連接收集直射信號,用來提供時間、衛(wèi)星仰角和接收機姿態(tài)等參數(shù)。另外的兩個射頻前端可以連接其他任意方向的左旋圓極化(LHCP)和右旋圓極化(RHCP)天線,接收反射信號。接收機有10個可編程的相關(guān)通道,可以根據(jù)實驗需要用多個通道以不同的多普勒頻移處理同一顆衛(wèi)星信號,也可以用單個相關(guān)通道以不同的幾何參數(shù)處理多達10個不同的衛(wèi)星。GOLD-RTR接收機的原始數(shù)據(jù)采樣頻率為40MHz,相關(guān)通道采樣頻率為20MHz,內(nèi)置的FPGA芯片可以實時地在每個通道每毫秒產(chǎn)生64個復(fù)數(shù)相關(guān)值,波形序列的延遲間隔為50ns(約15m),相關(guān)波形采用1bit量化。3.2海冰和堆積物理特性本次實驗是歐空局在格林蘭島開展的GNSSSIDS(SeaIceandDrySnow)項目的一部分,實驗的目的是驗證星載接收機接收GNSS反射信號研究海冰和積雪物理特性的可行性。如圖4所示,接收機天線架設(shè)在迪斯科海岸懸崖邊的一個電信信號塔上,距離海面高度約650m,直射天線朝向天頂,兩個反射天線水平放置朝向海面,可以通過網(wǎng)絡(luò)遠程控制。實驗從2008年10月下旬開始,到2009年5月初結(jié)束,連續(xù)觀測7個月,采集到了觀測海域的海冰從開始形成到完全結(jié)冰再到消融過程的全部數(shù)據(jù)。但是2008年10月下旬至當(dāng)年底的數(shù)據(jù)由于設(shè)備等原因?qū)嶒灁?shù)據(jù)不能夠使用,所以本文分析的是自2009年1月1日起的實驗數(shù)據(jù)。4實驗數(shù)據(jù)分析4.1極化比測定GNSS反射信號的不同極化類型提供了反射面的物理特性信息。為了序列除以直射右旋圓極化信號的峰值序列,并對處理結(jié)果取平均值,得到了如圖5所示的極化比均值。如圖5所示,反射LHCP/直射RHCP的極化比值有一個逐漸減小,然后逐漸增加的過程。4.2海水淡化后各參數(shù)f的放電特性海水在結(jié)冰之初,鹽分來不及完全排出,有一部分被凍結(jié)在海冰中;在形成海冰后,鹽分會逐漸排出去。如初生海冰的鹽度可以達到20左右,經(jīng)過一個冬天,一年冰的鹽度降到4—6。由于海冰在形成過程中表面包含海水,因此可以用海水介電常數(shù)估計海冰表面的介電常數(shù)。海水介電常數(shù)通常用德拜(Debye)公式表示為:其中:ε0=8.854×10-12F/m,為自由空間介電常數(shù);頻率無限大時ε∞=4.9;ω=2πf為電磁波角頻率;S,T分別為海水的鹽度和溫度(℃);電導(dǎo)率σ;靜態(tài)介電常數(shù)ε1和弛豫時間τ隨溫度T和鹽度S變化。根據(jù)德拜公式和海冰形成過程中海冰表面鹽度和溫度的變化,GPSL1波段的頻率為1575.42MHz,溫度變化范圍-20—0℃,鹽度變化范圍4—20,近似計算出2009年格林蘭島迪斯科海域1—5月的海冰表面介電常數(shù)變化曲線(如圖6所示)。對比GPS14、17、18、22、28號衛(wèi)星極化比均值(圖5)和海冰形成過程中復(fù)介電常數(shù)(實部)的變化趨勢圖(圖6)可知,海冰形成過程中海冰表面的復(fù)介電常數(shù)逐漸減小,而對應(yīng)的反射LHCP/直射RHCP極化比均值逐漸減小,反之,海冰融化過程中,表面的復(fù)介電常數(shù)逐漸增加,而對應(yīng)的反射LHCP/直射RHCP極化比均值逐漸增加,我們可以總結(jié)出反射LHCP/直射RHCP的極化比均值與海冰表面的介電常數(shù)具有相同的變化趨勢,隨著海冰的形成,反射LHCP/直射RHCP極化比均值會不斷減小,隨著海冰的消融,反射LHCP/直射RHCP極化比均值會不斷增加,所以我們可以通過GNSS-R技術(shù)觀測海面GPS反射信號的極化比值檢測海冰的形成和消融過程。由于圖6是我們通過根據(jù)德拜公式模擬出的2009年1—5月海冰形成過程中復(fù)介電常數(shù)(實部)的變化趨勢,在以后的研究中,我們將利用該海域的遙感圖像分析等手段,更進一步證實我們的研究理論。5實驗數(shù)據(jù)的使用利用歐洲空間局(ESA)于2008年10月下旬開始,到2009年5月初結(jié)束,在格林蘭島迪斯科海岸的GNSS反射信號實驗數(shù)據(jù),討論了基于GNSS反射信號用來檢測海冰的形成和消融過程的基本方法。由于2008年10月下旬至2008年底的數(shù)據(jù)不能夠使用,所以著重分析了從2009年1月1日起的實驗數(shù)據(jù)。跟蹤了20°左右的GPS衛(wèi)星,并且計算出反射LHCP/直射RHCP的極化比值,并與德拜公式模擬的海冰形成過程中復(fù)介電常數(shù)(實部)的變化趨勢進行了對比,初步得出了利