邵禮明， 張云， 孟婉婷， 袁國良

（1. 上海海事大學 上海 201306； 2. 上海海洋大學 上海 201306）

目前， 包括美國的GPS、 俄羅斯的GLONASS、 中國的BeiDou 等星座在內(nèi)的全球?qū)Ш叫l(wèi)星系統(tǒng) （Global Navigation Satellite System）， 通過穩(wěn)定的電磁波信號， 可以提供精確的定位、導航和授時信息，GNSS-R（GNSS-Reflection）以其大量的信號資源、低成本、全天時等特點，近年來成為國內(nèi)外的研究熱點。

干雪是處于高緯度、 高海拔地區(qū)擁有的特殊天氣現(xiàn)象，干雪與濕雪的區(qū)別主要在于所含的水分多少，干雪中不包含水分，例如，像南極洲這樣天氣極其干燥、全年降水量甚少的地方，適合干雪的存在。

現(xiàn)有的雪深檢測方法有很多，從傳統(tǒng)上講，可以通過人工測量工具，例如尺子或者標桿，現(xiàn)場對雪深進行測量，這種方法得到的結果非常準確，但是如果需要測量不同區(qū)域范圍的雪深時，該方法將會費時費力。

Proko[1]所采用的激光掃描儀（TLS，terrestrial laser scanner）的方法也可用來檢測雪深。 TLS 測量法獲得的雪深精確度受目標分辨率、光束直徑、掃描持續(xù)時間、雪的反射率和穿入雪層中激光雷達信號限制。 除此之外，假如雪表面是濕的，或者雪顆粒較大， 只有50%左右的TLS 發(fā)射脈沖可以被接收到，這會對最終結果造成很大影響。 但是通過TLS 方法，可以得到比較精確的雪深結果，誤差在10 cm 范圍內(nèi)。

在Lutz[2]等人研究中，對雪深測量采用的是空載雷射掃描法（ALS，airborne laser scanning）。 ALS 方法對大面積（至少幾十平方公里）檢測雪深的精度要求來說，是一個穩(wěn)定、成本劃算的方法。 但是，ALS 方法同樣面臨著一些挑戰(zhàn)：1）每次掃描范圍有限；2）必須對每次掃描的位置單獨進行測量；3）需要對另外區(qū)域進行再掃描，才可以得到雪深的連續(xù)覆蓋范圍。

文中采用基于GNSS-R 技術的方法，對南極洲干雪深度進行定性檢測。相比于其他雪深檢測實驗，GNSS-R 技術有著許多優(yōu)勢：1）由于發(fā)射信號源是衛(wèi)星，采用異源觀測模式，因此有著大量的信號源可以利用，可以在一次實驗中搭載多個接收機進行數(shù)據(jù)采集，接收機成本大大下降；2）反射信號直接攜帶反射面特性信息，反射信號的幅值、相位、極性和頻率等參量的變化，都直接與反射面有關；3）采用的擴頻技術具有干擾小、抗干擾性能強等特點，可以接收到淹沒在噪聲之中的低功率信號；4）對于一些環(huán)境惡劣的條件，例如：南極、北極、海面等地方，GNSS-R 檢測技術將帶來極大的方便性和安全性。因此，GNSS-R 技術對于海洋遙感、陸地遙感、災害預警等方面來說，有著大好前景和深遠意義。

1 實驗原理

在1993 年， 歐空局的科學家Martin-Neira 首次提出，對檢測到的GPS 的反射信號進行利用的概念[3]。 1994 年，法國科學家Auber 在進行機載飛行實驗時，意外接收到GPS 的反射信號，證明了反射信號是可以被接收到的[4]。 自此之后，歐洲、美國等發(fā)達國家陸續(xù)展開了利用衛(wèi)星反射信號反演地面參數(shù)的實驗，主要方向集中在：海面測高[5]、海面風場、海洋鹽度、土壤濕度等領域。

當發(fā)射信號入射到干雪區(qū)域時，由于干雪從表面到深處的幾百米處有著不同的雪層，每層干雪的相對介電常數(shù)不同（ε），信號會在干雪-空氣接觸面以及內(nèi)部不同雪層間交匯面發(fā)生反射， 或者再次入射， 達到更深的雪層。 在MRSR（multiple-ray single-reflection）模型中，如圖1 所示，衛(wèi)星信號平行入射，并遵循斯涅耳定律（n 表示折射率）。MRSR 模型表示從信號發(fā)射到接收的整個過程，每一條信號射線只進行單次反射，向上傳播，并到達接收機。

圖1 MRSR 模型示意圖Fig. 1 Sketch of MRSR model

1.1 第i 層干雪對信號造成的延遲

當信號射線傳播到第i 層，并在此層底部反射，向上傳播到接收機，這時反射信號和直射信號的延遲ρi，可以用迭代方法解決，在文中公式中，第0 層是空氣介質(zhì)，所以，0 層延遲表示雪-空氣表面對信號造成的延遲（圖1 中的S 點）：

下標TS，SR，TR 分別指的是發(fā)射機-鏡面反射點， 鏡面反射點-接收機，發(fā)射機-接收機。 那么，通過遞推可知，第i層的一般表達式是[6]:

1.2 第i 層干雪對反射信號幅度產(chǎn)生的影響

MRSR 模型的下一步是估算信號射線在第i 層內(nèi)部的幅度Ui。 首先，需要考慮以下兩個問題：

1）入射和反射的菲涅爾方程，入射系數(shù)用T 表示，反射系數(shù)用? 表示。

2）信號射線在每層干雪中的傳播衰減常數(shù)α[7]。

衛(wèi)星發(fā)射的GPS 信號是右旋圓極化的（RHCP，right-hand circular polarization）。 為了從圓極化信號獲得反射和入射系數(shù)，要求采用以下運算：

其中下標co 和cross 分別和共極化T‖和?‖交叉極化，和表示菲涅爾方程中的入射系數(shù)T⊥和?⊥反射系數(shù)的平行分量，和表示菲涅爾方程中的入射系數(shù)和反射系數(shù)的垂直分量。另外，要求考慮到信號向上傳播和向下傳播兩種情況：用Tij作為入射系數(shù)，表示衛(wèi)星信號從第i 層介質(zhì)和第j 層介質(zhì)界面穿過，需要指出的是Tij和Tji是兩個不同的概念，數(shù)值上并不相等。 ?ij指的是從第i 層介質(zhì)與第j 層介質(zhì)界面信號的反射系數(shù)。 因此，圖1 中的在接下來的分析中必須理解成?cross，T理解成Tco。

因為發(fā)射的GPS 信號經(jīng)過介質(zhì)反射后，由原來的右旋圓極化信號， 大多轉(zhuǎn)化為左旋圓極化信號 （LHCP，left-hand circular polarization），所以文中采用LHCP 天線接收的信號數(shù)據(jù)進行分析。

類似的，在第0 層空氣介質(zhì)中，經(jīng)過干雪-空氣接觸面反射的信號，相對于入射信號的振幅可以表示成：

那么，用迭代的方式，可以得到第k 層中信號振幅：

2 數(shù)據(jù)分析

2009 年12 月， 歐空局開展了一項關于海冰和干雪的研究項目（GPS－SIDS，Sea Ice and Dry Snow project），實驗目的是為了探究GPS－R 接收機接收的GPS 反射信號用于遙感海冰和積雪特性的可行性。 文中采用歐空局地球觀測中心提出的多路射線單次反射模型（MRSR）反演雪深的方法，對歐空局2009 年12 月20 日8 號衛(wèi)星采集GPS－SIDS 的實驗數(shù)據(jù)進行分析。

實驗所用的GNSS－R 接收機是GPS 開環(huán)差分實時接收機（GOLD－RTR）[8]，能夠通過相關器（lag）接收I、Q 兩路波形。在距離干雪表面高度約45 m 的接收塔上， 直射天線水平放置朝向天頂，接收GPS 發(fā)射的直射信號；兩個反射天線，左旋圓極化（LHCP）天線和右旋圓極化（RHCP）天線水平放置，接收從干雪表面和內(nèi)部反射回來的反射信號。

每個波形分別用64 個同相（I）和正交（Q）延遲器（lag）以1 kHz 速率進行采樣儲存，每個lag 空間分辨率為15 m，將獲得的波形進行1－s 積分，I、Q 支路結果的總波形如圖2 所示（顏色的深淺表示時間的不同），接收到的反射信號不是規(guī)則的三角波形，而是一系列歪曲的三角形，有著拖尾和次峰，這與MRSR 模型中不同雪層對信號造成不同延遲相對應。

圖2 同相和正交信號的總振幅Fig. 2 Total amplitude of in－phase and quadrature signals

從信號總振幅的時間序列可以看出， 波形在不同的lag中是振蕩的。 GPS C/A 碼碼片長度約為300 m，等于20 個lag總分辨率長度，一些來自更深干雪層中的延遲信號會超過碼片長度，也就是說，分析雪深時需要考慮所有相關器接收到的信號，而不僅僅是波峰所在位置，因此，對每個lag 接收到的信號單獨進傅里葉運算，并進行累加，可以大致得到干雪在不同深度的回應值，如圖3 所示，在表面的10 m 附近，內(nèi)部的85 m，135 m，225 m 附近有較明顯峰值出現(xiàn)。

3 實驗結果驗證

為了驗證真實數(shù)據(jù)得到的結果是否正確，可以用真實數(shù)據(jù)與相同幾何條件下的仿真結果進行對比，從而得出結論。

圖3 干雪深度的初步結果Fig. 3 Preliminary result of dry snow depth

通過仿真， 把不同雪層中得到的所有信號強度相加起來，可以得到接收機接收到的I、Q 支路信號波形，每層雪中的信號波形可以用GPS C/A 碼三角函數(shù)近似表示。 在仿真復波形之前，需要知道每層雪相對于直射信號產(chǎn)生的延遲ρi、雪層中信號的振幅Ui，以及接收到的信號相對直射信號產(chǎn)生的相位：

作為參考的直射信號的峰值設置在第22 個lag 處：

考慮到GPS C/A 碼的碼片寬度為20 個lag，也就是說C/A 碼自相關函數(shù)必須限制在±300 m 的范圍內(nèi)，因此，每個lag接收到的復波形可以表示成：

圖4 接收機接收到的仿真波形Fig. 4 The simulated waveform received by receiver

對仿真出的信號波形，按照處理真實數(shù)據(jù)的方法，在每個lag 中單獨進行傅里葉變換，并把所有延遲軸的傅里葉變換結果相加起來，可以近似得到雪深與信號振幅的關系，見圖5。圖中可以看到在表面的10 m 附近， 內(nèi)部的85 m、135 m 和265 m 附近有較多信號反射出去，到達接收機。 仿真結果與真實數(shù)據(jù)有很好的一致性，證明了真實數(shù)據(jù)的正確性。

4 結 論

圖5 仿真得到的雪深與振幅的關系Fig. 5 The relationship between snow depth and amplitude achieved by simulation

通過文中的數(shù)據(jù)分析，運用MRSR 模型初步反演出南極洲干雪的深度，結果與仿真數(shù)據(jù)有很好的一致性，初步驗證了GNSS－R 技術可以用來反演雪深的可行性，這為以后的干雪研究做了基礎性工作。 后續(xù)工作將繼續(xù)分析歐空局的南極洲干雪數(shù)據(jù)，并與其他方法反演雪深結果作比較。

GNSS－R 作為一個新興的、高精度、全方位、實時遙感技術，已經(jīng)顯示出對研究地球表面參數(shù)具有的潛力；將來，隨著國內(nèi)外發(fā)射的衛(wèi)星越來越多，衛(wèi)星的反射信號資源也會越來越豐富，這會給科學研究提供可靠依據(jù)，因此，GNSS－R 技術有著大好前景。

[1] Prokop A. Assessing the applicability of terrestrial laser scanning forspatial snow depth measurements[J]. Cold Reg.Sci. Technol.，2008，45（3）:155－163.

[2] Lutz E R，Geist T，St?tter J. Investigations of Airborne Laser Scanning Signal Intensity on glacial surfaces: utilizing comprehensive Laser Geometry Modeling and Surface Type Modeling；（a case study: Svartisheibreen， Norway）[M]. na，2003.

[3] Martin－Neira M. A passive reflectometry and interferometry system （PARIS）: Application to ocean altimetry [J]. ESA journal，1993（17）:331－355.

[4] Auber J C，Bibaut A，Rigal J M. Characterization of multipath on land and sea at GPS frequencies [C]//Proceedings of the 7th International Technical Meeting of the Satellite Division of The Institute of Navigation （ION GPS 1994）. 1994:1155－1171.

[5] Martín－Neira M，Caparrini M，F(xiàn)ont－Rossello J，et al. The PARIS concept: an experimental demonstration of sea surface altimetry using GPS reflected signals [J]. Geoscience and Remote Sensing， IEEE Transactions on，2001，39（1）:142－150.

[6] Fabra F. GNSS－R as a source of opportunity for sea ice and dry snow remote sensing [D]. Ph. D. Dissertation，Polytechnical University of Catalonia （UPC）， Barcelona，2013.

[7] Ulaby F T，Moore R K，F(xiàn)ung A K. Microwave Remote Sensing: Active and Passive [M]. Dedham， MA: Artech House，1986: 1797－1848.

[8] Nogués－Correig O，Galí E C，Campderrós J S，et al. A GPSreflections receiver that computes Doppler/delay maps in real time [J]. Geoscience and Remote Sensing， IEEE Transactions on，2007，45（1）:156－174.