2. 上海市空间导航与定位技术重点实验室, 上海 200030;
3. 中国科学院行星科学重点实验室, 上海 200030;
4. 中国科学院国家空间科学中心空间环境探测研究室, 北京 100190;
5. 北京市天基空间环境探测重点实验室, 北京 100190;
6. 中国科学院大学天文学院, 北京 100190;
7. 中国科学院国家授时中心, 陕西 西安 710600
2. Key Laboratory of Space Navigation and Position Technology, Shanghai 200030, China;
3. Key Laboratory of Planetary Sciences, Chinese Academy of Sciences, Shanghai 200030, China;
4. National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China;
5. Beijing Key Laboratory of Space Environment Exploration, Beijing 100190, China;
6. School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100190, China;
7. National Time Service Center, Chinese Academy of Sciences, Xi'an 710600, China
季节性冻土与永久性冻土约占地球陆地总面积的35%,主要分布在高纬度和高海拔地区。陆地表层土壤冻融状态转换随季节每年重复发生,与人类生活环境密切相关。由于土壤中水的相态变化,这一过程强烈地影响着地表辐射能量的转换过程和地表径流的强度,是地表能量平衡和水分平衡的重要影响因素, 也是气候变化的灵敏指示器。因此有效监测地表冻融状态的时空分布及其相关物理参数是冰冻圈、地学、水文学等研究领域的重要课题。
可见光和热红外遥感受天气条件限制,微波遥感却可以全天时全天候观测。主/被动微波遥感(雷达/辐射计)是地表冻融状态监测的有效手段之一。星载观测将空前提高地表冻融监测的空间分辨率,但其时间分辨率(每3天全球重复覆盖)与实际监测的科学需求之间存在一定的差距。
近年来,利用导航卫星的反射信号(GNSS-R)或直射信号与反射信号的相干信号(GNSS-IR)进行陆面参数遥感研究成为一种新兴的对地观测方式。GNSS-R/IR显著的低成本、小功耗、轻体积以及高时空分辨率等诸多优点可以成为现有地面站点观测和机载/星载遥感监测的有效有益补充[1-2]。
GNSS-R的工作模式分为散射计模式和高度计模式。针对GNSS-R散射计模式,需要研制专门GNSS反射信号接收机对陆面参数进行遥感监测,在该种接收机的基础上利用其相关功率值遥感地物参数。目前陆地表面研究领域涉及土壤水分、植被含水量和生物量以及积雪厚度研究。
GNSS-IR遥感不需要研制专门的接收机,可以直接使用地球测绘或者地球物理的接收机对地表土壤水分、植被参数、积雪厚度和潮位变化进行遥感监测[3-6]。GNSS-IR遥感的空间分辨率约为1 km,介于传统站点式传感器(< 1 m2)和星载观测(>100 km2)之间。利用从多路径观测数据中的有效反射计高度、幅度及相位信息可以提取相关地物参数。美国PBO H2O研究团队将上述环境遥感产品每隔12 h公布一次[7],其土壤水分反演结果已被用于SMAP(soil moisture active passive)卫星的校验[8]。
目前,利用GNSS-R/IR进行地表冻融状态监测的研究相对较少,利用GNSS-R/IR对地表冻融状态监测具有可行性[9],是一种新兴的监测手段。研究者已经利用相关理论和IGS(international GNSS service)站点数据分析中进行了一些探索性研究[10-11]。本文详细阐述了利用GNSS-R/IR遥感方式进行地表冻融监测的理论机理。
针对星载观测,GNSS卫星群和专门的GNSS-R反射信号接收机之间形成一种双站雷达的工作模式,接收机得到的DDM波形是用于反演地球物理参数的基本观测量。开发相应的仿真软件是开展卫星观测、数据仿真、试验设计,内插外推数据和反演等工作的关键机理工具[13-14]。目前,国内外的研究团队重点针对海况参数展开仿真分析,已经发展建立CYGNSS E2ES(cyclone global navigation satellite systems mission end to end simulator, E2ES)软件[15]、P2EPS(PAU/PARIS end-to-end performance simulator)仿真软件[16]和GREEPS(GNSS-R end-to-end performance simulator)仿真软件[17]、SAVERS(simulator of GNSS reflections from bare and vegetated soils)仿真软件则针对陆面土壤和植被特性研究而建立[18]。近期,研究人员在此Z-V模型基础上[19],发展建立了地表冻融转换DDM仿真分析模型。
大多数的前向GPS多路径模型关注码调制,对于反射能量、相位和延迟采用任意值,或者基于给定几何计算反射延迟,亦或根据经验值定义反射能量。基于物理模型的前向GPS多路径模型,最早的是文献[20]提出的模型。该模型完全考虑了测绘接收机和天线的直射信号和反射信号的极化特性,采用分层的办法计算反射信号,并利用小波扰模型计算反射率,将介电常数转换为反射系数,最后将反射系数和指定的天线辐射样图耦合得到SNRmpi(multipath signal noise ratio)干涉图。基于文献[20]的物理模型,文献[21]给出了一个全极化的前向GPS多路径模型。该模型考虑了GPS广播信号的右旋和左旋圆极化信息、天线朝向、地表响应、伪码测距和噪声能量。本文在此模型基础上,发展建立了地表冻融状态转换模型。
本文将已有的DDM仿真分析模型和前向GPS多路径模型拓展到适宜进行地表冻融状态监测的理论模型,并以此模型为工具,模拟分析了地表冻融状态转换对DDM波形和多路径数据的影响,从散射机理上揭示了利用GNSS-R/IR遥感进行地表冻融特性监测的理论依据,为后续利用GNSS-R/IR方法进行地表冻融监测提供理论依据。
1 地表冻融转换时的电磁散射特性变化物质的复介电常数是其固有属性,实部反映了不同介质表面发生的波的折射和反射现象;虚部则表示入射电磁波在介质中的衰减(吸收和转化)程度。本部分利用混合介电常数模型模拟分析了地表冻融转换时,介电常数以及其对应的各种极化条件下的反射信号的相干和非相干部分的散射特性变化情况。
1.1 冻融转换导致介电常数差异基于电磁波传输理论,冻结土壤可看作是由空气、固体颗粒、自由水、束缚水和冰5种物质组成的介电混合体,其复介电常数理论上受入射电磁波的频率、温度、土壤含水量等因素影响。最终介电常数可看做是各组分介电常数之和,如式(1)所示
式中,α是形状常数因子;Vi是第i个组分的体积百分含量。本文土壤介电常数的实部和虚部可通过普遍采用的半经验模型——Dobson介电常数模型计算得到[22-24]。图 1给出利用该模型模拟的不同土壤水分含量下,GPSL1载波频率,土壤介电常数实部和虚部随土壤温度的变化情况。从中可以看出当土壤从冻结状态转换到融化状态过程中,在各种土壤水分含量下,介电常数的实部和虚部均变化明显。在土壤温度大于0℃时,土壤水分含量对介电常数影响(实部,虚部)明显,在土壤温度小于0℃时,由于冰的存在,土壤水分对实部有影响但不明显,对于虚部的影响完全可以忽略。
对于不同的GPS载波频率,土壤的介电常数如表 1所示。由表 1可知在冻融转换时,GPS载波频率(均为低频L波段)对介电常数的影响可以忽略,因此在后续模拟分析中只采用L1载波频率时的介电常数开展分析。
土壤温度/℃ | L1载波频率 | L2载波频率 | L5载波频率 | |||||
实部 | 虚部 | 实部 | 虚部 | 实部 | 虚部 | |||
-1 | 8.667 69 | 1.034 03 | 8.707 93 | 1.168 76 | 8.712 93 | 1.197 62 | ||
+1 | 21.109 70 | 3.544 15 | 21.277 01 | 3.498 27 | 21.297 76 | 3.510 16 |
1.2 冻融转换时圆极化反射率变化情况
在考虑相位的情况下,利用极化合成的方法,计算得到RR和LR极化下的反射率[25-26]
式中,σ为镜像双站雷达散射截面,下角标R、L分别代表RHCP(right hand circular polarization)和LHCP(left hand circular polarization)极化状态;r为镜像反射系数,下标v和h分别代表垂直极化和水平极化。
图 2给出了当土壤温度分别为-1℃(冻结)和1℃(融化)时,圆极化和圆线极化在同极化和交叉极化下的BRCS(bistatic radar cross section)变化情况。由图 2可知,在布鲁斯诺角附近,VR极化存在BRCS最低点,VR极化在小仰角,即1℃布鲁斯诺角之前,-1℃的BRCS高于1℃,而在-1℃布鲁斯诺角之后,1℃时的BRCS高于-1℃的BRCS;在图 2中给出的仰角范围内,-1℃时的BRCS高于1℃的BRCS(两条BRCS交叉线之间),在BRCS交叉线之后和第二个布鲁斯诺角之前,1℃的BRCS高于-1℃;RR极化随着仰角的增加,BRCS增加,1℃的BRCS高于-1℃的BRCS;LR极化的BRCS随着仰角增加而降低,-1℃的BRCS高于1℃的BRCS。
1.3 冻融转换时BRCS变化
实际地表为随机粗糙表面,入射的导航卫星信号经地表散射后,散射能量中的非相干部分不可忽视。而常用的随机粗糙地表的面散射模型有KA(Kirchhoff approach)模型、SPM(small perturbation method)模型、IEM(integrated equation model)模型以及AIEM(advanced integrated equation model)模型。AIEM模型是在IEM模型的基础上的改进模型,形式更为复杂,但是模型的精度得到了提高。AIEM模型的基本形式可以概括为3项之和:基尔霍夫项σrtk、基尔霍夫补偿项σrtc和二者的交叉项σrtkc,如式(4)所示。式中σ0为双站雷达散射截面,q、p分别为接收和发射时的极化状态[27]
图 3是利用极化合成方法计算得到的地表冻融转换时各种极化BRCS的变化情况。在地表从冻结到融化转换的过程中,BRCS在各种极化和角度时,变化差异明显,即观测几何中的天顶角和方位角以及观测的各种极化都会导致冻融转换时BRCS的变化。如何从角度和极化信息中有效提取冻融转换敏感参数是后续研究重点内容之一。
2 冻融转化时GNSS-R/IR接收机信号变化 2.1 多路径信息变化
当地物参数,如积雪厚度、植被含水量、土壤水分变化时,会引起地球测绘或者地球物理的GPS接收机的多路径信息的变化(SNR、相位和伪距)。因此,可以利用直射信号和反射信号的相干多路径数据对地物参数进行反演。
文献[21]发展建立了可以同时考虑GPS信号极化、天线和地表响应的全极化前向GPS多路径模型,如式(5)和式(6)所示
式中,P为电场能量;G为天线增益;W为Woodward模糊函数;X是地表和天线耦合系数;下标d和r分别代表直射和反射分量,计算公式如下
式中,上标R和L分别为RHCP和LHCP极化;Φ为干涉相位;Rs和RX分别为同极化和交叉极化发射系数,是H极化和V极化的线性组合。
图 4给出了当地表冻融转换时,天线高度为2.0 m,载波频率为GPS L1时,采用TRM29659.00天线模型,模拟分析在上半部分为空气,下半部分为冻融土壤的多径反射情况下,多路径信息(SNR,相位和伪距)的变化情况,以及其引起的相应的有效反射计高度的变化。从图 4可以看出,当土壤温度发生冻结-融化转换时,会引起多路径幅度增加,以及有效反射计高度的增加。而这种冻融转换引起的多路径信息的变化,可以用于地表冻融转换监测。文献[12]中给出了实测GPS(AB33)和SNOTEL (ID 958)站点分析结果,并探讨了相应影响因素及局限性。
2.2 冻融转化时DDM波形变化
当地表发生冻融转换时,介电常数的变化会导致BRCS的变化,进而引起DDM波形的变化。本文在Z-V模型的基础上[19],针对BRCS采用2.3节的随机粗糙面模型,建立了适用于冻融土壤的DDM模型,并模拟分析了冻融转换前后,DDM波形的差异。Z-V模型本质上为双站雷达的积分形式,计算公式为
式中,PT为GNSS卫星发射能量,其波长为λ;天线增益及散射点到发射机/接收机之间的距离分别用G和R表示;下标T和R分别代表发射机和接收机;σo则是积分区域A内的冻融地表的双站雷达散射截面;三角函数和sin c函数分别用Λ和S表示,二者平方的积是双站雷达的模糊度函数。
当地表发生冻融转换时,DDM波形的变化如图 5所示(注:未进行归一化处理)。从中可以看出当冻融转换时,不同极化状态的DDM马蹄状波形存在着明显的差异。极化信息是进行冻融监测的重要指标信息,不同极化时,DDM波形的最大值和最小值如表 2所示。
DDM differences | 极化方式 | ||||
VV | HH | VH | RR | RL | |
max | 8.85e-18 | 7.28e-18 | 0 | 0 | 0 |
min | -1.54e-17 | -1.40e-17 | -4.57e-18 | -8.40e-17 | -9.93e-18 |
3 讨论
本文从仿真模型角度分析了利用GNSS-R/IR技术进行地表冻融监测的理论依据。
本文采用的冻融土介电常数模型为经过实际验证的较为经典的模型,理论模拟的结果与文章中的模拟完全一致[22-24]。各种极化的相干和非相干部分的双站散射特性,在模型验证时,采用的方法是将模型中椭倾角和椭率角修改为线极化角度,与原有模型进行比较,结果一致,因此本文用圆极化/圆线极化模拟分析冻融转换时BRCS差别。另外冻融GPS多路径模型和DDM波形的验证方法可参考文献[25]。
极化是电磁波的重要特性,地表反射信号的极化信息携带了地表的重要信息。极化比是土壤水分反演和植被状态研究的重要信息参数[28-30]。地表冻融过程中,不同极化的反射信号的双站散射特性和DDM波形信息的差异性表明,极化/极化比信息在地表冻融状态的监测领域具有重要前景。
GNSS-R/IR本质上是双站雷达。由于电磁波散射的空间异质性,在各个散射天顶角和方位角下,包含地物的不同信息。有效利用不同角度(天顶角和方位角)信息是提高地物监测时后向反演的重点和难点问题。
GNSS+R遥感是指除了可以利用GNSS导航卫星群的反射信号进行遥感监测外,也可以利用通信卫星或数字广播卫星的反射信号遥感地物参数[31-32]。本文模拟分析时重点考虑GPSL1载波频率下的散射特性、多路径以及DDM差异,但同时本文所提及模型,同样为利用其他数字通讯卫星开展地表冻融监测提供了基本理论模型和依据。
土壤水分、积雪和植被都是复杂寒区冻融地表的重要地表参数,直接影响着相应多路径和DDM波形的变化。本文中给出了裸土情况下,地表冻融状态对GNSS-R/IR接收机信号的影响,复杂寒区地表的积雪、植被覆盖等是后续模型中发展的重点。
4 结论GNSS-R和GNSS-IR是介于微波遥感和卫星导航技术的新兴多元交叉对地观测技术,近年来在国内外引起了广泛关注。本文将该种遥感方式的应用领域从现有的土壤水分、植被和积雪厚度研究拓展到地表冻融状态监测上。理论研究表明当地表发生冻融转换时,土壤介电常数存在较大差别,进而导致了各种极化下反射信号的相干分量和非相干分量在冻融前后变化差异明显。
[1] | ZAVOROTNY V U, GLEASON S, CARDELLACH E, et al. Tutorial on remote sensing using GNSS bistatic radar of opportunity[J]. IEEE Geoscience and Remote Sensing Magazine, 2014, 2(4): 8–45. DOI:10.1109/MGRS.2014.2374220 |
[2] | CARDELLACH E, FABRA F, NOGUÉS-CORREIG O, et al. GNSS-R ground-based and airborne campaigns for ocean, land, ice, and snow techniques:application to the GOLD-RTR data sets[J]. Radio Science, 2011, 46(6): RS0C04. |
[3] | LARSON K M. GPS interferometric reflectometry:applications to surface soil moisture, snow depth, and vegetation water content in the western United States[J]. Wiley Interdisciplinary Reviews:Water, 2016, 3(6): 775–787. DOI:10.1002/wat2.1167 |
[4] | CHEW C. Soil moisture remote sensing using GPS-interferometric reflectometry[D]. Boulder: University of Colorado, 2015. |
[5] |
汉牟田, 张波, 杨东凯, 等.
利用GNSS干涉信号振荡幅度反演土壤湿度[J]. 测绘学报, 2016, 45(11): 1293–1300.
HAN Mutian, ZHANG Bo, YANG Dongkai, et al. Soil moisture retrieval utilizing GNSS interference signal amplitude[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(11): 1293–1300. DOI:10.11947/j.AGCS.2016.20160145 |
[6] |
张双成, 南阳, 李振宇, 等.
GNSS-MR技术用于潮位变化监测分析[J]. 测绘学报, 2016, 45(9): 1042–1049.
ZHANG Shuangcheng, NAN Yang, LI Zhenyu, et al. Analysis of tide variation monitored by GNSS-MR[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(9): 1042–1049. DOI:10.11947/j.AGCS.2016.20150498 |
[7] | CHEW C C, SMALL E E, LARSON K M, et al. Effects of near-surface soil moisture on GPS SNR data:development of a retrieval algorithm for soil moisture[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(1): 537–543. DOI:10.1109/TGRS.2013.2242332 |
[8] | JACKSON T, COLLIANDER A, KIMBALL J, et al. Science data calibration and validation plan-SMAP[J]. 2016, 1, 6 |
[9] |
吴学睿, 金双根, 陶学林.一种地表冻融状态的监测方法及装置: 中国, 201410326132.3[P]. 2014-11-19. WU Xuerui, JIN Shuanggen, TAO Xuelin. Method and device for monitoring freeze-thaw state of earth surface: China, 201410326132.3[P]. 2014-11-19. |
[10] | WU Xuerui, JIN Shuanggen. Can we monitor the bare soil freeze-thaw process using GNSS-R?: a simulation study[C]//Proceedings Volume 9264, Earth Observing Missions and Sensors: Development, Implementation, and Characterization Ⅲ. Beijing: SPIE, 2014. |
[11] | WU Xuerui, CHANG Liang, JIN Shuanggen, et al. Initial results for near surface soil freeze-thaw process detection using GPS-interferometric reflectometry[C]//Proceedings of 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). Beijing: IEEE, 2016: 1989-1992. |
[12] | WU Xuerui, JIN Shuanggen, CHANG Liang. Monitoring bare soil freeze-thaw process using GPS-interferometric reflectometry:simulation and validation[J]. Remote Sensing, 2018, 10(1): 14. |
[13] | FUNG A K, CHEN K S. Microwave scattering and emission models and their applications[R]. Boston: Artech House, 1994. |
[14] | SHI Jiancheng, DU Yang, DU Jinyang, et al. Progresses on microwave remote sensing of land surface parameters[J]. Science China Earth Sciences, 2012, 55(7): 1052–1078. DOI:10.1007/s11430-012-4444-x |
[15] | O'BRIEN A, JOHNSON T J. Comparing the CYGNSS simulator forward scattering model with TDS-1 and CYGNSS on-orbit DDMS[C]//Proceedings of 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). Fort Worth, TX: IEEE, 2017: 2657-2658. |
[16] | PARK H, CAMPS A, PASCUAL D, et al. A generic level 1 simulator for spaceborne GNSS-R missions and application to GEROS-ISS ocean reflectometry[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(10): 4645–4659. DOI:10.1109/JSTARS.2017.2720625 |
[17] | BAI Weihua, XIA Junming, ZHAO Danyang, et al. GREEPS: an GNSS-R end-to-end performance simulator[C]//2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). Beijing: IEEE, 2016: 4831-4834. https://www.researchgate.net/publication/309773667_GREEPS_An_GNSS-R_End-to-End_Performance_Simulator |
[18] | PIERDICCA N, GUERRIERO L, GIUSTO R, et al. SAVERS:a simulator of GNSS reflections from bare and vegetated soils[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(10): 6542–6554. DOI:10.1109/TGRS.2013.2297572 |
[19] | ZAVOROTNY V U, VORONOVICH A G. Scattering of GPS signals from the ocean with wind remote sensing application[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(2): 951–964. DOI:10.1109/36.841977 |
[20] | ZAVOROTNY V U, LARSON K M, BRAUN J J, et al. A physical model for GPS multipath caused by land reflections:toward bare soil moisture retrievals[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2010, 3(1): 100–110. DOI:10.1109/JSTARS.2009.2033608 |
[21] | NIEVINSKI F G, LARSON K M. Forward modeling of GPS multipath for near-surface reflectometry and positioning applications[J]. GPS Solutions, 2014, 18(2): 309–322. DOI:10.1007/s10291-013-0331-y |
[22] | DOBSON M C, ULABY F T, HALLIKAINEN M T, et al. Microwave dielectric behavior of wet soil-part Ⅱ:dielectric mixing models[J]. IEEE Transactions on Geoscience and Remote Sensing, 1985, GE-23(1): 35–46. DOI:10.1109/TGRS.1985.289498 |
[23] | HALLIKAINEN M T, ULABY F T, DOBSON M C, et al. Microwave dielectric behavior of wet soil-part 1:empirical models and experimental observations[J]. IEEE Transactions on Geoscience and Remote Sensing, 1985, GE-23(1): 2534. |
[24] | ZHANG Lixin, SHI Jiancheng, ZHANG Zhongjun, et al. The estimation of dielectric constant of frozen soil-water mixture at microwave bands[C]//Proceedings of 2003 IEEE International Geoscience and Remote Sensing Symposium. Toulouse: IEEE, 2003: 2903-2905. https://www.researchgate.net/publication/224748467_The_estimation_of_dielectric_constant_of_frozen_soil-water_mixture_at_microwave_bands |
[25] | FUNG A K. Microwave scattering and emission models and their applications[M]. Boston, London: Artech House, 2009. |
[26] | ULABY F T, ELACHI C. Radar polarimetry for geoscience applications[M]. Norwood, MA: Artech House Inc., 1990. |
[27] | CHEN K S, WU T D, TSANG L, et al. Emission of rough surfaces calculated by the integral equation method with comparison to three-dimensional moment method simulations[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(1): 90–101. DOI:10.1109/TGRS.2002.807587 |
[28] | ZAVOROTNY V U, VORONOVICH A G. Bistatic GPS signal reflections at various polarizations from rough land surface with moisture content[C]//Proceedings of IEEE 2000 International Geoscience and Remote Sensing Symposium. Taking the Pulse of the Planet: the Role of Remote Sensing in Managing the Environment. Honolulu, HI: IEEE, 2000: 2852-2854. |
[29] | ZAVOROTNY V, MASTERS D, GASIEWSKI A, et al. Seasonal polarimetric measurements of soil moisture using tower-based GPS bistatic radar[C]//Proceedings of 2003 IEEE International Geoscience and Remote Sensing Symposium. Toulouse: IEEE, 2003: 781-783. https://www.researchgate.net/publication/4064154_Seasonal_polarimetric_measurements_of_soil_moisture_using_tower-based_GPS_bistatic_radar |
[30] | WU Xuerui, JIN Shuanggen. GNSS-Reflectometry:Forest canopies polarization scattering properties and modeling[J]. Advances in Space Research, 2014, 54(5): 863–870. DOI:10.1016/j.asr.2014.02.007 |
[31] | TORRES O, LAWRENCE R. Retrieval of reflected direct broadcast satellite (DBS) signals for Earth science applications[C]//IGARSS 2008-2008 IEEE International Geoscience and Remote Sensing Symposium. Boston, MA: IEEE, 2008, 5: 240-243. https://www.researchgate.net/publication/224383874_Retrieval_of_Reflected_Direct_Broadcast_Satellite_DBS_Signals_for_Earth_Science_Applications |
[32] | SHAH R, GARRISON J L, GRANT M S. Demonstration of bistatic radar for ocean remote sensing using communication satellite signals[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(4): 619–623. DOI:10.1109/LGRS.2011.2177061 |