《測繪學報》
構建與學術的橋樑 拉近與權威的距離
金雙根1, 張勤耘1,2, 錢曉東1,3
1. 中國科學院上海天文臺, 上海 200030;
2. 上海大學通信與資訊工程學院, 上海 200444;
3. 中國科學院大學, 北京 100049
收稿日期：2017-05-20；修回日期：2017-09-01
基金專案：國家自然科學基金（11373059）
第一作者簡介：金雙根(1974—)， 男， 博士， 研究員， 研究方向為衛星導航與空間大地測量。 E-mail:sgjin@shao.ac.cn
摘要：全球導航衛星系統（GNSS）具有全天候、近即時、高精度的特點， 可持續發射L波段信號， 廣泛應用於定位、導航和授時（PNT）。 隨著GNSS技術的發展， 最近GNSS反射信號可探測地球表面特徵， 即GNSS反射測量（GNSS+R）。 結合GNSS接收機天線位置和介質資訊， 利用延遲測量值可以確定表面粗糙度和表面特性。 GNSS+R作為當前GNSS和遙感領域的研究熱點， 取得了一些研究進展和成果。 本文詳細介紹了GNSS+R原理和方法及其最新應用進展， 包括各種GNSS+R技術手段和方法， 以及海洋、陸地、水文、植被和冰雪等遙感應用， 特別是最新BeiDou-R和TDS-1研究進展。
New Progress and Application Prospects of Global Navigation Satellite System Reflectometry (GNSS+R)
JIN Shuanggen1, ZHANG Qinyun1,2, QIAN Xiaodong1,3
Abstract: The Global Navigation Satellite System (GNSS) has the characteristics with all-weather, near real-time and high precision observations. It can transmit L band signals continuously, which has been widely used in positioning, navigation and timing (PNT). With the development of GNSS technology, recently GNSS reflected signals can be used to detect Earth's surface characteristics, i.e., GNSS reflectometry (GNSS+R). Combined with GNSS receiver antenna position and medium information, surface roughness and characteristics can be determined by using delay measurements. As a hotspot in GNSS and remote sensing field, GNSS+R has achieved a number of progress and results. In this paper, the principles and methods of GNSS+R are presented as well as their latest advances, including various techniques, methods and remote sensing applications in the ocean, land, hydrological, vegetation and ice and snow, especially the latest progresses in BeiDou-R and TDS-1 GNSS+R. Finally, the future applications and prospective of GNSS+R are given, including multi-GNSS system, new GNSS-R receivers, GNSS+R missions and emerging applications.
Key words: GNSS multipath reflectometry environment remote sensing
全球導航衛星系統(GNSS)， 包括美國全球定位系統(GPS)、中國北斗導航衛星系統(BDS)、俄羅斯格洛納斯系統(GLONASS)和歐洲伽利略導航衛星系統(Galileo)， 以及區域增強系統， 如日本准天頂衛星系統(QZSS)和印度區域導航衛星(IRNSS)等， 具有全天候、近即時、高精度的特點， 可持續發射L波段信號， 廣泛應用於定位、導航和授時(PNT)。 隨著各導航衛星系統的逐漸完善， 星座的增多， 觀測站的增加， 其應用領域越來越廣泛。 GNSS系統不僅僅應用於定位、授時和導航， 還可以利用其表面反射信號進行遙感研究。 GNSS衛星持續向地球播發無線電信號， 其中部分信號會被地球表面反射。 從粗糙表面反射回來的GNSS延遲信號可以提供直射和反射信號路徑的不同資訊。
GPS反射測量是由歐空局(ESA)Martin Neira[1]於1993年首次提出來， 即GPS地表反射信號和直射信號一起被接收機接收， 它們之間延遲可以用於干涉測量， 即PARIS(passive reflectometry and interferometry system)。 1994年， 法國科學家通過進行飛行試驗發現接收機可以接收到海面反射信號， 但由於其對定位精度的影響， 通常將其作為多路徑信號剔除[2]。 1996年， NASA蘭利研究中心的科學家利用雙頻GPS信號進行海面前向散射以獲取並剔除電離層延遲，
本文詳細描述了GNSS+R原理和方法及其最新應用進展， 包括遙感海洋、陸地、水文、植被和冰雪特徵等， 特別是最新BDS-R和TDS-1 GNSS+R研究進展。 最後給出了將來GNSS+R的應用前景和展望， 包括多GNSS系統、GNSS+R接收機、GNSS+R衛星計畫和新興應用， 如中小尺度地表特徵和目標探測等。
1 GNSS+R原理與方法
GNSS+R反射測量技術屬於雙基雷達， 可獲得地表粗糙特徵和地球物理參數， 即利用GNSS測量直射信號與地表鏡面反射的信號之間延遲(時間延遲或相位延遲)， 再根據GNSS衛星、接收機和鏡面反射點之間的幾何位置關係， 可反演地表特徵。按照處理資料的方式，GNSS+R測高可以分為傳統型GNSS+R測高(cGNSS+R測高)和干涉GNSS+R測高(iGNSS+R)。前者是配置左右圓極化天線並利用接收機記錄的直射信號與反射信號的載波相位資料，通過固定模糊度和解算接收機鐘差等方式，確定兩者之間的傳播路徑延遲，進而計算天線至地球表面的高度。後者是利用直射信號與反射信號功率波形相關的原理，測得信號時延，進而計算天線到地球表面的垂直距離，但其涉及複雜的多普勒時延演算法，資料處理方法複雜。
1.1 傳統型GNSS+R
這種測高一般是採用傳統的大地測量型接收機，利用載波相位觀測值作為原始觀測量，得到直射信號與反射信號的傳播路程差，再根據幾何關係得到反射面的高度。如在湖邊、海岸或海島上架設GNSS接收機，進行GNSS+R測高試驗，這就是岸基測高。其所架設的接收機高度都較低，其天線的照射面積決定有效的水面散射面積，可有效代替傳統的驗潮站測量模式。GNSS接收機對所接收到的直射信號與表面反射信號進行相關，並通過測量相關函數的最大值的位置得到時間延遲T，進而可求得接收機到反射面的高度H。在飛機或氣球上假設GNSS接收機進行的GNSS+R測高試驗就是機載測高。與岸基測高相比，機載測高架設的高度較高，觀測的水面面積較大。在LEO衛星上搭載GPS接收機，進行GNSS+R測高試驗，這就是星載測高。在500~800 km的高空可以採用LEO衛星搭載GPS接收機來進行海面高度測量，這種測量方式與傳統衛星測高相比不需要發射機，也可以直接採用多個LEO衛星組成星座，具有較高的時空解析度。
1.2 干涉型GNSS+R
這種測高是利用特製的可以同時接收直射信號和反射信號的接收機，並將接收到的信號在接收機中進行相關處理，利用時延一維相關函數、多普勒一維相關函數或者時延-多普勒二維相關函數得到兩個信號之間時間延遲，再根據幾何關係得到反射面高度。
1.2.1 SNR技術
多路徑效應是GNSS高精度定位的主要誤差，它與反射面的結構和電介質參數密切相關。當衛星高度角低於10°時，GNSS接收到的反射信號是右旋極化。這時具有相同頻率的反射信號與直射信號會發生相干作用。這一相干現象反映在信噪比SNR的變化上，信噪比SNR觀測值是衡量GPS接收機天線接收到的信號的大小的一個量值，反映多路徑與多路徑誤差的大小受衛星信號的發射功率、天線增益、衛星與接收機間的距離及多路徑效應等因素的影響。在高度角較高的情況下，天線增益較大使得SNR得到有效提高；而在高度角較低的情況下，一方面天線增益減小，另一方面多路徑效應影響使得SNR下降較為嚴重。因此，對信噪比SNR進行分析可以評估多路徑效應，進而估計地表環境參數。
1.2.2 DDM技術
GNSS+R海洋遙感主要利用反射信號時延多普勒二維相關功率。其原理是通過計算本地載波信號與散射區域內不同時間延遲和多普勒頻率的接收信號的相關功率值。DDM多用於機載和星載GNSS+R，考慮到其距離地面高度較高，首先介紹GNSS+R的鏡面反射測量幾何關係。GNSS+R幾何關係要用到鏡面反射點，即從反射區域反射的反射信號中路徑延遲最短的理論反射點。根據GNSS衛星、接收機和鏡面反射點的幾何關係建立如圖的本地坐標系。該坐標系的原點為鏡面反射點，z軸為地球切面的法線方向，GNSS衛星T，鏡麵點和接收機R位於yz平面內，x軸按右手定則確定。圖 1中，ht為衛星到地球參考橢球面的高度；hr為接收機到地球參考橢球面的高度；Re為地球半徑；G為衛星到地心的距離；L為接收機到地心的距離；Rt和Rr分別是衛星和接收機到鏡面反射點的距離；β是地心、鏡面反射點和GNSS衛星間的夾角；α是地心、鏡面反射點和接收機間的夾角；θ為鏡麵點處衛星的仰角。
圖 1 GNSS鏡面反射測量幾何關係Fig. 1 Geometric relation of GNSS specular reflectiometry
圖選項
利用鏡面反射點與接收機和GNSS衛星之間的關係，以及GNSS反射信號時延多普勒二維相關功率關係可獲得地表反射面參數。如對於海洋表面，文獻[8]在2000年提出了完整的GNSS信號散射理論模型，其模型是關於時延和多普勒的二維函數，通過此模型可以得到散射區域不同散射點的反射功率，Z-V模型具體表達如下
(1)
式中，Y(τ,f)是接收機處理後的信號相關功率，是關於延遲和多普勒的函數；Pt是衛星發射的信號功率；Gt和Gr是衛星天線和接收機天線增益；Ti是信號處理中相干積分時間；Rt和Rr是衛星和接收機到地面反射點的距離；σ0是與反射面參數有關的項，包括反射係數和粗糙度等；S(f)是多普勒頻移函數；Λ(τ)是GNSS碼相關函數；A是有效的散射區域，大概是閃耀區大小。
2 GNSS+R最新應用進展2.1 海洋測高
目前岸基傳統型大地測量GNSS接收機獲得較好測高結果。如文獻[4]利用反射信號碼相位資訊，測高精度達到數米。文獻[9]在Crater湖面上利用iGNSS+R設備，解算了取樣速率為1 Hz的水面高，精度大約在2 cm。文獻[10—11]利用反射信號載波相位資訊，將精度提高到釐米級。文獻[12]利用GNSS+R SNR資料反演出的測高結果與驗潮站的觀測結果做了對比，結果顯示兩者之差小於3 cm。另外，文獻[13]在飛機上搭載iGNSS+R設備的試驗表明，利用GPS反射信號測高的精度能達到5 cm，其空間解析度能達到5 km。文獻[14]將cGNSS+R設備搭載在飛行高度為1 km的飛機上，利用L1載波相位信號解算得到的平均海面高度(SSH)與Jason-1的測量對比，精度在10 cm左右。
目前國內外利用岸基和機載GNSS+R等獲得了一系列的新結果。例如文獻[15]首次利用岸基BDS觀測獲得了海平面變化，與驗潮站觀測結果有較好的一致(見圖 2)。在星載GNSS+R方面，文獻[16]首次利用TDS-1資料反演得到海表面高度(sea surface height, SSH), 其結果與全球DTU10的海表面高度資料進行了對比，存在一定的一致性但依然也有較大的差異。
圖 2 利用BDS L2、L6、L7 SNR以及三頻組合所監測的海平面變化與潮汐觀測值的對比(10 d期間)[15]Fig. 2 Comparison of the monitored sea level changes with tide observations based on the SNR data of BDS L1, L6, L7 and three-frequency combinations during 10 days[15]
圖選項
2.2 海面測風
海面在風的作用下會引起波浪，從而引起海面粗糙度的變化。根據海面風場反演原理，反射信號的功率峰值隨著風速的增大而減小，波形也相應越扁平。這說明海面散射能量隨風速的變化而分佈到更大的範圍，即閃爍區範圍變大。故可以通過和模擬波形相對比反演海面風速。此外，風速和風向還會影響海平面坡度。風向的變化會對閃爍區的非對稱性產生影響。沿順風方向，海面散射功率曲線函數的包絡較寬，沿測風方向則相對較窄，表現為曲線波形上的非對稱性。
自1997年，NASA蘭利研究中心與科羅拉多大學聯合開展了關於GNSS+R海面風場反演演算法的研究，並進行多次機載飛行試驗[3]。試驗發現，GPS反射信號的相關功率波形與風速有直接關係，且海面風速越大，波形的後延變化越趨於平緩，反之越陡峭。2000年，文獻[8]提出雙尺度表面模型(Z-V模型)在早期的GNSS+R海面風場反演模型中較為系統和成熟，但是該模型並沒考慮閃爍區之外由海水表面小尺度坡度造成的布拉格散射所產生的反射。後來，文獻[17]等對Z-V模型進行了改進。除此之外，NASA/科羅拉多大學、ESA和Starlab研究所等進行了大量機載試驗，並根據反演理論對試驗資料進行了海面風場反演，風速反演精度達±2 m/s，風向精度為±20°。2014年，文獻[18]使用GNSS+R方法獲得的時延多普勒圖像(DDM)的5個觀測量，並使用根據NDBC浮標資料建立的經驗地球物理模型函數(GMF)，建立了基於回歸分析的最小方差風速估算模型，結果表明使用最小方差風速估算模型的均方根誤差比僅僅使用單個觀測值進行反演的均方根誤差要小得多。2016年，文獻[19]又針對熱帶氣旋全球導航衛星系統(CYGNSS)的第二等級產品為資料來源，並從樣本、功率的選擇，使用二維地球物理模型函數等方面，改進之前他所提出的海洋測風演算法，設置理想的CYGNSS在軌參數模擬得到的時延多普勒圖像(DDM)對演算法進行測試。同時，文獻[20]第一次利用TDS-1資料反演得到了海洋風速結果，該結果與實測資料達到了很好的一致性，如圖 3是利用TDS-1衛星資料反演得到的海洋風速結果。
圖 3 TDS-1海洋風速反演結果[20]Fig. 3 TDS-1 inversion results of sea wind speed[20]
圖選項
2.3 海冰探測
近年來，全球氣候發生顯著變化，冰川和凍土加劇融化，從而導致海平面產生明顯變化。早期的海冰監測常通過設立觀測站和雷達站，並利用海上工具進行實地監測，結合使用聲呐技術以及光學測量技術獲取海冰厚度，但解析度低和成本高。隨著遙感技術的發展，利用光學遙感和雷達遙感測冰的技術得以發展，但前者會受雲層和光照等因素的影響，而微波遙感則具有良好的穿透性且可全天候全天時工作等優勢。近些年，全球導航衛星系統的發展與完善，使得GNSS+R反射測量技術在海冰監測方面的研究也取得了顯著的成果。
相較於國外，國內對海冰監測的進展還處於起步階段。2013年，文獻[26]利用歐空局格陵蘭島的資料進行類比分析，得到GNSS反射信號的極化比與海冰密集度的關係，驗證了GNSS+R監測海冰的可行性。文獻[27]在中國渤海進行岸基試驗，對渤海海冰進行監測，通過類比GNSS+R信號在海水和海冰表面的反射，建立3層輻射傳輸模型得到GNSS信號的反射率和海冰厚度及入射角度間的關係，並將結果與2013年天津海冰消融過程的試驗資料進行對比分析，進一步驗證了GNSS+R信號對海冰監測的敏感。
2.4 土壤濕度
土壤濕度是地表土壤的一個重要參數。土壤水分在全球水迴圈中有著重要作用，是地表能量平衡的重要決定因素，影響著地表農作物的生長和洪水泥石流災害的發生，也影響著全球天氣變化。地面GPS反射信號包含有土壤濕度、介質常數、地面粗糙度以及植被覆蓋等方面的資訊。2000年，文獻[28]延伸了Z-V模型，使其適用於土壤濕度測量，發現反射信號的波形峰值雖然與土壤濕度相關，但是地面粗糙度對波形尾部有影響，這給利用波形反演土壤濕度帶來困難，於是產生利用特製雙天線接收機來接收信號從而利用反射信號來反演土壤濕度的構想。2002年以後，美國宇航局(NASA)連續執行了一系列GPS反射信號土壤水分遙感試驗，分析土壤介電常數與反射信號功率的定性關係。西班牙Starlab研究所設計了SMIGOL(soil moisture interference-pattern GNSS observations at L-band reflectometer)探測裝置，用於探測直射信號和反射信號干涉之後的信號，通過分析干涉信號波形特徵與土壤水分的關係，反演地表參數，即干涉模式技術(IPT)。文獻[29]利用IPT技術，從最簡單的裸土模型開始研究土壤水分反演演算法，到有植被覆蓋的複雜模型，證實土壤水分與干涉信號的振幅有關。以及文獻[30]又提出利用傳統的GPS接收機也可以進行土壤濕度反演，並對此進行了長期研究，與其他科研學者一起建立了基於物理的多路徑正演模型[31-32]，很好地分析了多路徑信噪比的相位振幅隨土壤濕度的變化而變化，促進了GNSS+R技術在土壤濕度監測應用。文獻[33]在河南鄭州開展的GNSS+R航空飛行試驗, 給出土壤濕度的估算方法和結果分析, 為中國GNSS+R觀測試驗開展、自主載荷研製、陸面觀測資料處理等提供參考。但由於地面上植物、樹葉及玻璃雜物等因素的存在，GPS多路徑信號包含的資訊是非常複雜的，從地基GPS多路徑信號中提取土壤水分參數的資訊，首先必須剔除其他因素的影響。
2.5 植被探測
植被探測，包括植物含水量，可以確保農作物的正常生長，及時發現存在的問題並解決。植物覆蓋在土壤上，其含水量會影響土壤濕度的監測，而估計植物含水量可以提高土壤濕度的反演精度。隨著GNSS系統的完善，GNSS台站已經遍佈全球各地，同一時刻可觀測的衛星個數、可接收到的衛星信號越來越多，有助於植物監測網的建立。而且GNSS信號屬於微波波段，而微波對地表粗糙度、含水量等敏感。一方面有利於對植物含水量的監測，但另一方面在監測植物含水量的同時，受地表粗糙度和土壤濕度的影響，因此在估計植物含水量時，需要採用一定的方法消除這些因素的影響。
文獻[30]首次提出了利用傳統GPS接收機多路徑觀測可監測土壤濕度，為GPS反射信號監測植物含水量奠定了基礎。文獻[34]首次利用GPS的雜訊統計量MP1 RMS定性地估計了植物生長，指出信噪比會隨著植被的生長而減小。文獻[35]用Larson等提出的用於土壤濕度反演的正演模型對植物含水量與信噪比和實際反射面高度進行了定量分析，研究發現在植被含水量不超過1 kg/m2時，植物含水量與信噪比振幅呈線性關係。文獻[36]提出了一種較為複雜的基於信噪比干涉圖的振幅與頻率分析的方法，以消除植被含水量在土壤濕度中的影響，並獲得較好的結果。
2.6 積雪厚度
積雪是重要的淡水資源，對全球大氣和海洋的熱狀況和區域性氣候有著重要影響。積雪的融化一方面補給了地球水資源，但是另一方面卻會造成海平面上升，淹沒沿海城市，影響人類正常生活。積雪探測十分重要，但是由於積雪的時間、空間變化特性，使得積雪探測十分困難。傳統的地基觀測方法雖然可以觀測到積雪厚度和密度，觀測精度也高，但是依然存在時間解析度較低或缺乏空間動態變化等缺陷，而新型的空基遙感探測，如光學星載感測器雖然可以提供積雪覆蓋資訊，但是積雪的厚度和密度等資訊卻無法獲得。而近年來，逐步發展起來的GNSS技術，利用其優勢在積雪探測方面取得了一定的成果。
文獻[37]利用GNSS+R技術對科羅拉多發生的兩次暴風雪的測量在積雪厚度方面的研究取得了初步結果，將GPS信噪比估計得到的積雪厚度與實際的積雪厚度進行 對比，證明傳統的GPS接收機可以用來測量積雪厚度。文獻[32]建立的基於物理的多路徑正演模型，充分考慮了GPS發射信號的左旋和右旋極化以及天線/地面的回應特性。他們還利用Matlab建立了該模型的模擬器，使用該模擬器可以用來分析反射測量的特徵。並利用多路徑正演模型建立了GPS多路徑測量積雪厚度的反演模型，估計結果與實測值達到了0.97以上的相關性[38]。文獻[39]利用GPS-L4觀測值和非參數的自舉模型來估計積雪厚度的變化，與實測結果符合度良好，進一步推動了積雪探測的發展；文獻[39]利用GPS的L2P觀測值的信噪比資料分析積雪厚度的變化，在高度角為5°—30°時計算得到的相關係數高達0.98，表明當GPS衛星高度角在30°以下時更有利於雪深測量。同時，文獻[41]又利用GLONASS信噪比觀測值與L4相位組合觀測值分析了積雪厚度結果，並與GPS結果進行了對比，同時融合兩者觀測值提高了估計結果精度(見圖 4)。
圖 4 GLONASS與GPS觀測值融合後的結果[41]Fig. 4 Fusion results of GLONASS and GPS observation[41]
圖選項
3 前景與展望
目前，新一代的全球導航衛星系統(GNSS)，包括美國現代化GPS、俄羅斯現代化GLONASS系統、歐盟Galileo系統和中國北斗系統(BDS)，正在不斷發展與完善，GNSS正向多頻、多模的方向發展。結合空基增強系統(QZSS和IRNSS等)，地基GNSS測網可以接收更多頻率的多系統GNSS反射信號，這將大大提高時空解析度，更好地估計地表反射面地球物理參數。另外，隨著全球IGS測站以及區域性GPS測站的不斷增多，可獲得更多的地基GNSS測站周邊環境資訊，為全球水文研究及氣候變化提供了一個新的工具。
3.1 GNSS+R接收機
美國噴氣推進實驗室(JPL)研發的“BlackJack”GPS接收機目前廣泛應用於精密定軌及無線電掩星等。為滿足未來新任務以及准即時資料處理能力的需求，新型GNSS接收機不僅能夠接收和跟蹤GPS信號，還將能夠接收Galileo、GLONASS及中國北斗等系統信號。目前JPL正在研發新一代的GNSS多頻接收機TriG(Tri-GNSS, GPS+Galileo+GLONASS)用於精密定軌及無線電掩星和反射測量。通過跟蹤GPS的L1 C/A、L2、L2C和L5以及Galileo、GLONASS等新型系統的信號，TriG可以進行多頻GNSS信號折射和反射測量。
為利用GPS反射信號推演海面高度、海況及土壤濕度等相關地球物理參數，加泰羅尼亞理工大學(UPC)研發了“高級GNSS無源反射儀(griPAU)”，griPAU設備利用GPS L1波段的C/A碼，可以即時、高精度的得到延遲-多普勒相關值(DDM)，並可根據不同的解析度配置及選擇相關或非相關的積分時間來計算24(32複雜點的DDM。griPAU的高靈敏度提高了地球物理參數的恢復品質。現代化的多頻GPS、未來的Galileo系統等帶來了更多的衛星星座、新的信號和頻段，西班牙的ICE也正在進行GNSS+R設備研發、GNSS信號散射和反射應用等方面的研究。
3.2 GNSS+R低軌衛星
星載GPS反射試驗，如美國國家海洋和大氣管理局(NOAA)和國家太空中心(NSPO)計畫在2017年發射下一代的FORMOSAT-7/COSMIC-2衛星。新一代的衛星系統將擁有12顆裝載GNSS+R接收機的衛星，可接收GPS、GLONAS以及Galileo反射信號。因此它可以收集大量的反射測量資料以用於陸地和海洋等方面的應用。同時，NASA設計的CYGNSS反射測量衛星2016年發射。CYGNSS含有8個繞著運載工具運行的小衛星，在該軌道上，小衛星將接收到來自GPS衛星的直射信號和由反射信號[42]。NASA團隊致力於研究海洋表面特性、潮濕大氣熱力學、判斷熱帶颶風的形成和是否會繼續加強以及加強幅度多少的對流輻射動力學之間的關係，而CYGNSS將會改善預報和記錄方法。由於快速變化的風暴中心在颶風形成初期和颶風的加劇中扮演著十分重要的角色，CYGNSS將首次幫助科學家探測發生在風暴中心的大氣與海洋之間的交互過程。
3.3 新興應用
確定洋面粗糙度和冰凍圈有很多種方法，然而缺少足夠的星載資料來驗證這些方法以及監視更多的細節，而傳統方法即時性不足且空間解析度較低。在未來幾年裡，將有更多搭載下一代GNSS遙感接收機的試驗衛星發射。在環境遙感的各個領域，高時空解析度、准即時的GNSS+R將得到更多的應用，如監測海冰、海況、海洋渦流、海面浮油、地質災害及太空梭測控等。這些准即時的GNSS反射資料將向公共領域開放，而對這些資料的分析將在預報強風、危險海況、洪澇災害、海洋渦流以及風暴潮等方面發揮重要作用。此外，高取樣速率的海面粗糙度數據也將從未來大量的GNSS反射信號中估計出來，這將加深人們對洋面海氣通量、浮冰脊、霜花、破冰以及冰雪交界處的粗糙度的認識，尤其是在那些無法接近且環境惡劣的海冰區域。
合成孔徑雷達(SAR)可生成關於陸地表面的高解析度影像，其發射機和接收機可以在同一平臺即單基SAR，也可分別安裝在不同平臺即雙基或多基SAR。對於GNSS-SAR，由星載發射機和鄰近地表的接收機等組成，其中GNSS衛星作為發射源，這些發射源可以是GPSGLONASSGalileoBDS，接收機可以固定在地表、裝置在航空器上或在陸面車輛上。與其他衛星系統相比，作為發射源的GNSS既有優點也有缺點。首先GNSS擁有眾多的衛星星座群，對於地球表面的同一點，同時有4~8個衛星可見，來自于多顆衛星的信號有助於建立3D射線表面圖，其次GNSS的導航信號具有最佳的遙感同步的特性。但是GNSS的預算功率低而且距離解析度受發射信號的頻寬限制。
4 小結
GNSS+R反射測量遙感技術是利用全球導航衛星系統(GNSS)的反射信號對海洋、陸地或冰川雪地進行遙感探測的新興手段。GNSS衛星作為可提供免費且長期穩定的L波段信號源，可充分發揮GNSS系統自身優勢-全天候、全天時、覆蓋範圍廣、時空解析度高等，在海洋遙感上，可以進行海洋測高、反演海面風場、估計海水鹽度以及海面溢油；在陸地遙感上，利用微波波段對水分的敏感的特性，可以估計土壤濕度和植物生長量；在冰川雪地遙感上，充分利用GNSS在時空解析度上的優勢，可以測量海冰厚度、積雪厚度、密度、粗糙度等。目前，GNSS+R在海洋、土壤濕度和冰雪等遙感監測等取得較好的進展。隨著GNSS+R反射測量技術的發展，將來或可能監測火山、地震形變和滑坡等自然災害。隨著越來越多的多頻GNSS導航衛星星座及空基增強系統，以及空基GNSS反射測量任務的實施(如接下來的FORMOSAT-7/COSMIC-2和CYGNSS任務)，人們將獲取更多的高時空解析度的地表特徵資訊。另外，相關人員正在研發更加先進的GNSS接收機，為了能夠滿足不同的應用需求而改進了演算法，為滿足未來空基高性能任務的需求而具有准即時資料處理能力(例如具備多模GNSS反射和折射技術的下一代Tri-GNSS接收機)。未來幾年裡，公眾使用一些大學等機構研發的低成本衛星將成為可能，GNSS信號在遙感領域的應用也將擴展到全球範圍內。
【引文格式】金雙根，張勤耘，錢曉東。全球導航衛星系統反射測量(GNSS+R)最新進展與應用前景[J]. 測繪學報，2017，46(10)：1389-1398. DOI: 10.11947/j.AGCS.2017.20170282
可反演地表特徵。按照處理資料的方式，GNSS+R測高可以分為傳統型GNSS+R測高(cGNSS+R測高)和干涉GNSS+R測高(iGNSS+R)。前者是配置左右圓極化天線並利用接收機記錄的直射信號與反射信號的載波相位資料，通過固定模糊度和解算接收機鐘差等方式，確定兩者之間的傳播路徑延遲，進而計算天線至地球表面的高度。後者是利用直射信號與反射信號功率波形相關的原理，測得信號時延，進而計算天線到地球表面的垂直距離，但其涉及複雜的多普勒時延演算法，資料處理方法複雜。1.1 傳統型GNSS+R
這種測高一般是採用傳統的大地測量型接收機，利用載波相位觀測值作為原始觀測量，得到直射信號與反射信號的傳播路程差，再根據幾何關係得到反射面的高度。如在湖邊、海岸或海島上架設GNSS接收機，進行GNSS+R測高試驗，這就是岸基測高。其所架設的接收機高度都較低，其天線的照射面積決定有效的水面散射面積，可有效代替傳統的驗潮站測量模式。GNSS接收機對所接收到的直射信號與表面反射信號進行相關，並通過測量相關函數的最大值的位置得到時間延遲T，進而可求得接收機到反射面的高度H。在飛機或氣球上假設GNSS接收機進行的GNSS+R測高試驗就是機載測高。與岸基測高相比，機載測高架設的高度較高，觀測的水面面積較大。在LEO衛星上搭載GPS接收機，進行GNSS+R測高試驗，這就是星載測高。在500~800 km的高空可以採用LEO衛星搭載GPS接收機來進行海面高度測量，這種測量方式與傳統衛星測高相比不需要發射機，也可以直接採用多個LEO衛星組成星座，具有較高的時空解析度。
1.2 干涉型GNSS+R
這種測高是利用特製的可以同時接收直射信號和反射信號的接收機，並將接收到的信號在接收機中進行相關處理，利用時延一維相關函數、多普勒一維相關函數或者時延-多普勒二維相關函數得到兩個信號之間時間延遲，再根據幾何關係得到反射面高度。
1.2.1 SNR技術
多路徑效應是GNSS高精度定位的主要誤差，它與反射面的結構和電介質參數密切相關。當衛星高度角低於10°時，GNSS接收到的反射信號是右旋極化。這時具有相同頻率的反射信號與直射信號會發生相干作用。這一相干現象反映在信噪比SNR的變化上，信噪比SNR觀測值是衡量GPS接收機天線接收到的信號的大小的一個量值，反映多路徑與多路徑誤差的大小受衛星信號的發射功率、天線增益、衛星與接收機間的距離及多路徑效應等因素的影響。在高度角較高的情況下，天線增益較大使得SNR得到有效提高；而在高度角較低的情況下，一方面天線增益減小，另一方面多路徑效應影響使得SNR下降較為嚴重。因此，對信噪比SNR進行分析可以評估多路徑效應，進而估計地表環境參數。
1.2.2 DDM技術
GNSS+R海洋遙感主要利用反射信號時延多普勒二維相關功率。其原理是通過計算本地載波信號與散射區域內不同時間延遲和多普勒頻率的接收信號的相關功率值。DDM多用於機載和星載GNSS+R，考慮到其距離地面高度較高，首先介紹GNSS+R的鏡面反射測量幾何關係。GNSS+R幾何關係要用到鏡面反射點，即從反射區域反射的反射信號中路徑延遲最短的理論反射點。根據GNSS衛星、接收機和鏡面反射點的幾何關係建立如圖的本地坐標系。該坐標系的原點為鏡面反射點，z軸為地球切面的法線方向，GNSS衛星T，鏡麵點和接收機R位於yz平面內，x軸按右手定則確定。圖 1中，ht為衛星到地球參考橢球面的高度；hr為接收機到地球參考橢球面的高度；Re為地球半徑；G為衛星到地心的距離；L為接收機到地心的距離；Rt和Rr分別是衛星和接收機到鏡面反射點的距離；β是地心、鏡面反射點和GNSS衛星間的夾角；α是地心、鏡面反射點和接收機間的夾角；θ為鏡麵點處衛星的仰角。
圖 1 GNSS鏡面反射測量幾何關係Fig. 1 Geometric relation of GNSS specular reflectiometry
圖選項
利用鏡面反射點與接收機和GNSS衛星之間的關係，以及GNSS反射信號時延多普勒二維相關功率關係可獲得地表反射面參數。如對於海洋表面，文獻[8]在2000年提出了完整的GNSS信號散射理論模型，其模型是關於時延和多普勒的二維函數，通過此模型可以得到散射區域不同散射點的反射功率，Z-V模型具體表達如下
(1)
式中，Y(τ,f)是接收機處理後的信號相關功率，是關於延遲和多普勒的函數；Pt是衛星發射的信號功率；Gt和Gr是衛星天線和接收機天線增益；Ti是信號處理中相干積分時間；Rt和Rr是衛星和接收機到地面反射點的距離；σ0是與反射面參數有關的項，包括反射係數和粗糙度等；S(f)是多普勒頻移函數；Λ(τ)是GNSS碼相關函數；A是有效的散射區域，大概是閃耀區大小。
2 GNSS+R最新應用進展2.1 海洋測高
目前岸基傳統型大地測量GNSS接收機獲得較好測高結果。如文獻[4]利用反射信號碼相位資訊，測高精度達到數米。文獻[9]在Crater湖面上利用iGNSS+R設備，解算了取樣速率為1 Hz的水面高，精度大約在2 cm。文獻[10—11]利用反射信號載波相位資訊，將精度提高到釐米級。文獻[12]利用GNSS+R SNR資料反演出的測高結果與驗潮站的觀測結果做了對比，結果顯示兩者之差小於3 cm。另外，文獻[13]在飛機上搭載iGNSS+R設備的試驗表明，利用GPS反射信號測高的精度能達到5 cm，其空間解析度能達到5 km。文獻[14]將cGNSS+R設備搭載在飛行高度為1 km的飛機上，利用L1載波相位信號解算得到的平均海面高度(SSH)與Jason-1的測量對比，精度在10 cm左右。
目前國內外利用岸基和機載GNSS+R等獲得了一系列的新結果。例如文獻[15]首次利用岸基BDS觀測獲得了海平面變化，與驗潮站觀測結果有較好的一致(見圖 2)。在星載GNSS+R方面，文獻[16]首次利用TDS-1資料反演得到海表面高度(sea surface height, SSH), 其結果與全球DTU10的海表面高度資料進行了對比，存在一定的一致性但依然也有較大的差異。
圖 2 利用BDS L2、L6、L7 SNR以及三頻組合所監測的海平面變化與潮汐觀測值的對比(10 d期間)[15]Fig. 2 Comparison of the monitored sea level changes with tide observations based on the SNR data of BDS L1, L6, L7 and three-frequency combinations during 10 days[15]
圖選項
2.2 海面測風
海面在風的作用下會引起波浪，從而引起海面粗糙度的變化。根據海面風場反演原理，反射信號的功率峰值隨著風速的增大而減小，波形也相應越扁平。這說明海面散射能量隨風速的變化而分佈到更大的範圍，即閃爍區範圍變大。故可以通過和模擬波形相對比反演海面風速。此外，風速和風向還會影響海平面坡度。風向的變化會對閃爍區的非對稱性產生影響。沿順風方向，海面散射功率曲線函數的包絡較寬，沿測風方向則相對較窄，表現為曲線波形上的非對稱性。
自1997年，NASA蘭利研究中心與科羅拉多大學聯合開展了關於GNSS+R海面風場反演演算法的研究，並進行多次機載飛行試驗[3]。試驗發現，GPS反射信號的相關功率波形與風速有直接關係，且海面風速越大，波形的後延變化越趨於平緩，反之越陡峭。2000年，文獻[8]提出雙尺度表面模型(Z-V模型)在早期的GNSS+R海面風場反演模型中較為系統和成熟，但是該模型並沒考慮閃爍區之外由海水表面小尺度坡度造成的布拉格散射所產生的反射。後來，文獻[17]等對Z-V模型進行了改進。除此之外，NASA/科羅拉多大學、ESA和Starlab研究所等進行了大量機載試驗，並根據反演理論對試驗資料進行了海面風場反演，風速反演精度達±2 m/s，風向精度為±20°。2014年，文獻[18]使用GNSS+R方法獲得的時延多普勒圖像(DDM)的5個觀測量，並使用根據NDBC浮標資料建立的經驗地球物理模型函數(GMF)，建立了基於回歸分析的最小方差風速估算模型，結果表明使用最小方差風速估算模型的均方根誤差比僅僅使用單個觀測值進行反演的均方根誤差要小得多。2016年，文獻[19]又針對熱帶氣旋全球導航衛星系統(CYGNSS)的第二等級產品為資料來源，並從樣本、功率的選擇，使用二維地球物理模型函數等方面，改進之前他所提出的海洋測風演算法，設置理想的CYGNSS在軌參數模擬得到的時延多普勒圖像(DDM)對演算法進行測試。同時，文獻[20]第一次利用TDS-1資料反演得到了海洋風速結果，該結果與實測資料達到了很好的一致性，如圖 3是利用TDS-1衛星資料反演得到的海洋風速結果。
圖 3 TDS-1海洋風速反演結果[20]Fig. 3 TDS-1 inversion results of sea wind speed[20]
圖選項
2.3 海冰探測
近年來，全球氣候發生顯著變化，冰川和凍土加劇融化，從而導致海平面產生明顯變化。早期的海冰監測常通過設立觀測站和雷達站，並利用海上工具進行實地監測，結合使用聲呐技術以及光學測量技術獲取海冰厚度，但解析度低和成本高。隨著遙感技術的發展，利用光學遙感和雷達遙感測冰的技術得以發展，但前者會受雲層和光照等因素的影響，而微波遙感則具有良好的穿透性且可全天候全天時工作等優勢。近些年，全球導航衛星系統的發展與完善，使得GNSS+R反射測量技術在海冰監測方面的研究也取得了顯著的成果。
相較於國外，國內對海冰監測的進展還處於起步階段。2013年，文獻[26]利用歐空局格陵蘭島的資料進行類比分析，得到GNSS反射信號的極化比與海冰密集度的關係，驗證了GNSS+R監測海冰的可行性。文獻[27]在中國渤海進行岸基試驗，對渤海海冰進行監測，通過類比GNSS+R信號在海水和海冰表面的反射，建立3層輻射傳輸模型得到GNSS信號的反射率和海冰厚度及入射角度間的關係，並將結果與2013年天津海冰消融過程的試驗資料進行對比分析，進一步驗證了GNSS+R信號對海冰監測的敏感。
2.4 土壤濕度
土壤濕度是地表土壤的一個重要參數。土壤水分在全球水迴圈中有著重要作用，是地表能量平衡的重要決定因素，影響著地表農作物的生長和洪水泥石流災害的發生，也影響著全球天氣變化。地面GPS反射信號包含有土壤濕度、介質常數、地面粗糙度以及植被覆蓋等方面的資訊。2000年，文獻[28]延伸了Z-V模型，使其適用於土壤濕度測量，發現反射信號的波形峰值雖然與土壤濕度相關，但是地面粗糙度對波形尾部有影響，這給利用波形反演土壤濕度帶來困難，於是產生利用特製雙天線接收機來接收信號從而利用反射信號來反演土壤濕度的構想。2002年以後，美國宇航局(NASA)連續執行了一系列GPS反射信號土壤水分遙感試驗，分析土壤介電常數與反射信號功率的定性關係。西班牙Starlab研究所設計了SMIGOL(soil moisture interference-pattern GNSS observations at L-band reflectometer)探測裝置，用於探測直射信號和反射信號干涉之後的信號，通過分析干涉信號波形特徵與土壤水分的關係，反演地表參數，即干涉模式技術(IPT)。文獻[29]利用IPT技術，從最簡單的裸土模型開始研究土壤水分反演演算法，到有植被覆蓋的複雜模型，證實土壤水分與干涉信號的振幅有關。以及文獻[30]又提出利用傳統的GPS接收機也可以進行土壤濕度反演，並對此進行了長期研究，與其他科研學者一起建立了基於物理的多路徑正演模型[31-32]，很好地分析了多路徑信噪比的相位振幅隨土壤濕度的變化而變化，促進了GNSS+R技術在土壤濕度監測應用。文獻[33]在河南鄭州開展的GNSS+R航空飛行試驗, 給出土壤濕度的估算方法和結果分析, 為中國GNSS+R觀測試驗開展、自主載荷研製、陸面觀測資料處理等提供參考。但由於地面上植物、樹葉及玻璃雜物等因素的存在，GPS多路徑信號包含的資訊是非常複雜的，從地基GPS多路徑信號中提取土壤水分參數的資訊，首先必須剔除其他因素的影響。
2.5 植被探測
植被探測，包括植物含水量，可以確保農作物的正常生長，及時發現存在的問題並解決。植物覆蓋在土壤上，其含水量會影響土壤濕度的監測，而估計植物含水量可以提高土壤濕度的反演精度。隨著GNSS系統的完善，GNSS台站已經遍佈全球各地，同一時刻可觀測的衛星個數、可接收到的衛星信號越來越多，有助於植物監測網的建立。而且GNSS信號屬於微波波段，而微波對地表粗糙度、含水量等敏感。一方面有利於對植物含水量的監測，但另一方面在監測植物含水量的同時，受地表粗糙度和土壤濕度的影響，因此在估計植物含水量時，需要採用一定的方法消除這些因素的影響。
文獻[30]首次提出了利用傳統GPS接收機多路徑觀測可監測土壤濕度，為GPS反射信號監測植物含水量奠定了基礎。文獻[34]首次利用GPS的雜訊統計量MP1 RMS定性地估計了植物生長，指出信噪比會隨著植被的生長而減小。文獻[35]用Larson等提出的用於土壤濕度反演的正演模型對植物含水量與信噪比和實際反射面高度進行了定量分析，研究發現在植被含水量不超過1 kg/m2時，植物含水量與信噪比振幅呈線性關係。文獻[36]提出了一種較為複雜的基於信噪比干涉圖的振幅與頻率分析的方法，以消除植被含水量在土壤濕度中的影響，並獲得較好的結果。
2.6 積雪厚度
積雪是重要的淡水資源，對全球大氣和海洋的熱狀況和區域性氣候有著重要影響。積雪的融化一方面補給了地球水資源，但是另一方面卻會造成海平面上升，淹沒沿海城市，影響人類正常生活。積雪探測十分重要，但是由於積雪的時間、空間變化特性，使得積雪探測十分困難。傳統的地基觀測方法雖然可以觀測到積雪厚度和密度，觀測精度也高，但是依然存在時間解析度較低或缺乏空間動態變化等缺陷，而新型的空基遙感探測，如光學星載感測器雖然可以提供積雪覆蓋資訊，但是積雪的厚度和密度等資訊卻無法獲得。而近年來，逐步發展起來的GNSS技術，利用其優勢在積雪探測方面取得了一定的成果。
文獻[37]利用GNSS+R技術對科羅拉多發生的兩次暴風雪的測量在積雪厚度方面的研究取得了初步結果，將GPS信噪比估計得到的積雪厚度與實際的積雪厚度進行 對比，證明傳統的GPS接收機可以用來測量積雪厚度。文獻[32]建立的基於物理的多路徑正演模型，充分考慮了GPS發射信號的左旋和右旋極化以及天線/地面的回應特性。他們還利用Matlab建立了該模型的模擬器，使用該模擬器可以用來分析反射測量的特徵。並利用多路徑正演模型建立了GPS多路徑測量積雪厚度的反演模型，估計結果與實測值達到了0.97以上的相關性[38]。文獻[39]利用GPS-L4觀測值和非參數的自舉模型來估計積雪厚度的變化，與實測結果符合度良好，進一步推動了積雪探測的發展；文獻[39]利用GPS的L2P觀測值的信噪比資料分析積雪厚度的變化，在高度角為5°—30°時計算得到的相關係數高達0.98，表明當GPS衛星高度角在30°以下時更有利於雪深測量。同時，文獻[41]又利用GLONASS信噪比觀測值與L4相位組合觀測值分析了積雪厚度結果，並與GPS結果進行了對比，同時融合兩者觀測值提高了估計結果精度(見圖 4)。
圖 4 GLONASS與GPS觀測值融合後的結果[41]Fig. 4 Fusion results of GLONASS and GPS observation[41]
圖選項
3 前景與展望
目前，新一代的全球導航衛星系統(GNSS)，包括美國現代化GPS、俄羅斯現代化GLONASS系統、歐盟Galileo系統和中國北斗系統(BDS)，正在不斷發展與完善，GNSS正向多頻、多模的方向發展。結合空基增強系統(QZSS和IRNSS等)，地基GNSS測網可以接收更多頻率的多系統GNSS反射信號，這將大大提高時空解析度，更好地估計地表反射面地球物理參數。另外，隨著全球IGS測站以及區域性GPS測站的不斷增多，可獲得更多的地基GNSS測站周邊環境資訊，為全球水文研究及氣候變化提供了一個新的工具。
3.1 GNSS+R接收機
美國噴氣推進實驗室(JPL)研發的“BlackJack”GPS接收機目前廣泛應用於精密定軌及無線電掩星等。為滿足未來新任務以及准即時資料處理能力的需求，新型GNSS接收機不僅能夠接收和跟蹤GPS信號，還將能夠接收Galileo、GLONASS及中國北斗等系統信號。目前JPL正在研發新一代的GNSS多頻接收機TriG(Tri-GNSS, GPS+Galileo+GLONASS)用於精密定軌及無線電掩星和反射測量。通過跟蹤GPS的L1 C/A、L2、L2C和L5以及Galileo、GLONASS等新型系統的信號，TriG可以進行多頻GNSS信號折射和反射測量。
為利用GPS反射信號推演海面高度、海況及土壤濕度等相關地球物理參數，加泰羅尼亞理工大學(UPC)研發了“高級GNSS無源反射儀(griPAU)”，griPAU設備利用GPS L1波段的C/A碼，可以即時、高精度的得到延遲-多普勒相關值(DDM)，並可根據不同的解析度配置及選擇相關或非相關的積分時間來計算24(32複雜點的DDM。griPAU的高靈敏度提高了地球物理參數的恢復品質。現代化的多頻GPS、未來的Galileo系統等帶來了更多的衛星星座、新的信號和頻段，西班牙的ICE也正在進行GNSS+R設備研發、GNSS信號散射和反射應用等方面的研究。
3.2 GNSS+R低軌衛星
星載GPS反射試驗，如美國國家海洋和大氣管理局(NOAA)和國家太空中心(NSPO)計畫在2017年發射下一代的FORMOSAT-7/COSMIC-2衛星。新一代的衛星系統將擁有12顆裝載GNSS+R接收機的衛星，可接收GPS、GLONAS以及Galileo反射信號。因此它可以收集大量的反射測量資料以用於陸地和海洋等方面的應用。同時，NASA設計的CYGNSS反射測量衛星2016年發射。CYGNSS含有8個繞著運載工具運行的小衛星，在該軌道上，小衛星將接收到來自GPS衛星的直射信號和由反射信號[42]。NASA團隊致力於研究海洋表面特性、潮濕大氣熱力學、判斷熱帶颶風的形成和是否會繼續加強以及加強幅度多少的對流輻射動力學之間的關係，而CYGNSS將會改善預報和記錄方法。由於快速變化的風暴中心在颶風形成初期和颶風的加劇中扮演著十分重要的角色，CYGNSS將首次幫助科學家探測發生在風暴中心的大氣與海洋之間的交互過程。
3.3 新興應用
確定洋面粗糙度和冰凍圈有很多種方法，然而缺少足夠的星載資料來驗證這些方法以及監視更多的細節，而傳統方法即時性不足且空間解析度較低。在未來幾年裡，將有更多搭載下一代GNSS遙感接收機的試驗衛星發射。在環境遙感的各個領域，高時空解析度、准即時的GNSS+R將得到更多的應用，如監測海冰、海況、海洋渦流、海面浮油、地質災害及太空梭測控等。這些准即時的GNSS反射資料將向公共領域開放，而對這些資料的分析將在預報強風、危險海況、洪澇災害、海洋渦流以及風暴潮等方面發揮重要作用。此外，高取樣速率的海面粗糙度數據也將從未來大量的GNSS反射信號中估計出來，這將加深人們對洋面海氣通量、浮冰脊、霜花、破冰以及冰雪交界處的粗糙度的認識，尤其是在那些無法接近且環境惡劣的海冰區域。
合成孔徑雷達(SAR)可生成關於陸地表面的高解析度影像，其發射機和接收機可以在同一平臺即單基SAR，也可分別安裝在不同平臺即雙基或多基SAR。對於GNSS-SAR，由星載發射機和鄰近地表的接收機等組成，其中GNSS衛星作為發射源，這些發射源可以是GPSGLONASSGalileoBDS，接收機可以固定在地表、裝置在航空器上或在陸面車輛上。與其他衛星系統相比，作為發射源的GNSS既有優點也有缺點。首先GNSS擁有眾多的衛星星座群，對於地球表面的同一點，同時有4~8個衛星可見，來自于多顆衛星的信號有助於建立3D射線表面圖，其次GNSS的導航信號具有最佳的遙感同步的特性。但是GNSS的預算功率低而且距離解析度受發射信號的頻寬限制。
4 小結
GNSS+R反射測量遙感技術是利用全球導航衛星系統(GNSS)的反射信號對海洋、陸地或冰川雪地進行遙感探測的新興手段。GNSS衛星作為可提供免費且長期穩定的L波段信號源，可充分發揮GNSS系統自身優勢-全天候、全天時、覆蓋範圍廣、時空解析度高等，在海洋遙感上，可以進行海洋測高、反演海面風場、估計海水鹽度以及海面溢油；在陸地遙感上，利用微波波段對水分的敏感的特性，可以估計土壤濕度和植物生長量；在冰川雪地遙感上，充分利用GNSS在時空解析度上的優勢，可以測量海冰厚度、積雪厚度、密度、粗糙度等。目前，GNSS+R在海洋、土壤濕度和冰雪等遙感監測等取得較好的進展。隨著GNSS+R反射測量技術的發展，將來或可能監測火山、地震形變和滑坡等自然災害。隨著越來越多的多頻GNSS導航衛星星座及空基增強系統，以及空基GNSS反射測量任務的實施(如接下來的FORMOSAT-7/COSMIC-2和CYGNSS任務)，人們將獲取更多的高時空解析度的地表特徵資訊。另外，相關人員正在研發更加先進的GNSS接收機，為了能夠滿足不同的應用需求而改進了演算法，為滿足未來空基高性能任務的需求而具有准即時資料處理能力(例如具備多模GNSS反射和折射技術的下一代Tri-GNSS接收機)。未來幾年裡，公眾使用一些大學等機構研發的低成本衛星將成為可能，GNSS信號在遙感領域的應用也將擴展到全球範圍內。
【引文格式】金雙根，張勤耘，錢曉東。全球導航衛星系統反射測量(GNSS+R)最新進展與應用前景[J]. 測繪學報，2017，46(10)：1389-1398. DOI: 10.11947/j.AGCS.2017.20170282