4 [3 _3 D! l p) k5 K2 d 点击蓝字,关注我们
% Y# R, K& x; H) V2 ] 全文链接:http://www.sdkx.net/article/2024/1002-4026/1002-4026-37-2-1.shtml
, w W: x6 o( u; C# R" Z 基于全球卫星导航系统的海洋环境监测方法研究综述
. Q: h' q1 A) T3 c" w! Z/ } 摘要:实时、准确、可靠地监测海洋环境信息,对于海洋灾害预警预报、防灾减灾、海洋资源开发、保障海洋安全等具有至关重要的作用。近年来,随着全球卫星导航系统(GNSS)不断发展与升级,基于GNSS导航信号进行大气和海洋环境信息探测发展成为新技术新方法,并成为海洋环境监测技术研究热点,已在海洋气象监测、数值预报、科学研究等领域得到了广泛应用。系统综述了GNSS技术在海面有效波高、风速、降雨强度、水汽及潮位监测等海洋环境监测领域的应用研究现状,对新技术新方法进行了系统总结,并对未来发展趋势进行了展望,以期对相关领域的研究提供参考。 , M6 |# c$ H3 h/ |; ?
关键词: 全球卫星导航系统; 海洋环境; 有效波高; 降雨强度; 海面风; 水汽监测 , p- @: ]: R' e8 I
我国是一个海洋大国,海洋在经济发展和国家安全中的作用愈加突显 [1-2]。海洋环境监测是关心海洋、认识海洋和经略海洋的重要保障,为政府和业务部门提供决策依据,保障海洋活动安全利益和人员生命财产安全 [3-4]。海洋资料浮标等海洋环境原位监测手段,具有可靠性高、实时准确等特点,可基本满足常规海洋监测需求[5],微波辐射计、雷达高度计和无线电探空仪等新技术也被广泛用于海洋环境监测中,但是这些技术却存在空间监测范围有限等不足,卫星遥感海洋环境监测虽然具有空间范围广等优点,但其空间分辨率较低[6-8]。 2 n2 @# Y) Y6 Z$ X, R/ ~
近30年,随着全球卫星导航系统(global navigation satellite system,GNSS)飞速发展,在用于全球定位、导航和授时服务的同时,开始用于监测海洋环境信息[9-13]。二十世纪九十年代,欧洲空间局学者Martín-Neira [14]首次提出利用GNSS反射信号(GNSS reflectometry,GNSS-R)反演海况信息的概念。美国及我国科学家也都相继开展了GNSS-R的实验验证和理论研究,并已将其成功应用于海面高度、有效波高和海面风场等估算 [15-20]。此外,美国科学家Bevis [21-22]提出利用GNSS直射导航信号反演大气水汽信息,形成了GNSS气象学概念,其后基于GNSS的大气水汽监测研究得到国内外众多学者的关注。
. O: g, {1 z/ _( I$ y( u4 z6 Z 鉴于GNSS信息在海洋区域环境信息监测中的重要作用,本文综述了基于GNSS的海面有效波高、风速、降雨强度、水汽及潮位监测等方面的研究进展,并对未来发展趋势进行了展望,以期对本领域的研究人员有所借鉴。 5 K; b; D1 z* Q* K% V
01
0 e+ K2 K8 Z; ?+ y3 L1 F. U/ ? GNSS有效波高反演 " X& p/ @1 F. w, v/ Z2 q2 Y
有效波高(significant wave height,HSWH)是多个海浪波高值从小到大排序后,前三分之一海浪波高值的平均值[23],如图1所示。有效波高是重要的海洋环境监测信息,它在船舶调度、航行安全及实时信息感知中具有极其重要的作用。
0 k, r# u& L; e' { 图1 有效波高的定义
! ]! a& l+ R \6 R$ d: b3 u
* |/ C) k" y7 o. S; x Fig.1 Definition of significant wave height   
5 q* P) l5 ~/ j. v( M" @7 I, ? 图2 基于GNSS-R的有效波高估算结果 3 C5 }' s2 S8 j. R
9 z ^( R) g0 u1 T
Fig.2 Estimation of significant wave height based on GNSS-R technology
. u: ?8 d! d- \9 S. B+ S 02 / v% y4 v7 W' }# Z" X+ e4 i Q
GNSS海面风场反演 
3 C0 _9 Q* H4 U8 g. K 围绕Z-V模型,国内外科研机构学者陆续开展了机载平台、星载平台和船载平台的GNSS-R反演海面风速研究[29-35]。此外,基于单天线GNSS-IR技术,王笑蕾等[36]详细介绍了地基GNSS-IR估算海面风速原理,证实了可以利用SNR估算海面风速。近年来随着人工智能技术的飞速发展,人工神经网络、卷积神经网络和多模态深度学习网络等机器学习算法被逐渐应用于GNSS-R海面风速反演当中,并取得了精度较高的估算结果[37-41]。 
8 R" x+ l8 }6 [/ B) y 图3 捕风一号卫星海面风速反演结果
& q& T" A7 J$ e* u5 j
J5 F; F0 I/ G( w7 ]* m Fig. 3 Sea-surface wind speed retrieved from BuFeng-1
8 a+ b7 \) h' E& \( O 03 & r: x4 R7 V" V4 ^( ?6 z
GNSS海面降雨强度反演
! D1 w# r6 U! d 降雨强度信息是重要的海洋环境参数之一,也是天气预报、气象研究和航行安全的重要支撑。近年来国内外相关学者逐步开展了GNSS-R反演降雨信息研究。 @' x! @" q( o! q
Asgarimehr等[43]开创了基于GNSS-R技术获取海面降雨信息研究的先河,该研究分析了GNSS-R时延多普勒图(delayed Doppler map,DDM)数据中的降雨信息特征,证明了由于雨滴落到海表会改变海洋表面波浪,GNSS-R能够反演弱风条件下海面降雨信息。随后,Balasubramaniam等[44]利用CYGNSS DDM数据也开展了类似研究。图4给出了基于GNSS-R估算降雨强度的示意图[43]。 ' @4 q5 F/ [0 x! W
图4 海面降雨强度反演示意图
( @' ^: _* s6 [% o/ K0 g+ T b6 s3 r* t
Fig.4 Diagram of sea-surface rainfall intensity retrieval
, t9 `' f3 ~# W8 f2 I# P7 }% ^ 此外,布金伟[45]对星载GNSS-R技术反演降雨强度方法进行了系统研究,提出了三维快匹配滤波算法的DDM降噪模型,利用支持向量机、随机森林和卷积神经网络等机器学习算法实现了海面降雨强度检测,并取得了良好的估算精度。
( T* @* h6 A2 a0 k) [ 04 % c; D# I4 c. z: h
GNSS海面水汽反演
$ |3 v8 @& Q7 U& W7 @ 基于GNSS技术解算大气可降水量(precipitable water vapor, PWV)参数的技术主要有双差网解法和精密单点定位(precise point positioning,PPP)两种[8]。由于海面多为移动平台,无法满足双差网解的观测条件,因此其暂时无法满足海面PWV信息探测要求。而PPP技术仅需单台设备即可实现全球范围灵活作业[46],不需要额外布设基准站,模型简单,解算效率高,因而自其概念提出以来,就被广泛应用到了海面PWV信息探测 [47-48]。基于PPP技术的PWV估算及研究进展简述如下。
& r4 \8 n( I+ }, g% c* S0 ? GNSS伪距和载波相位观测值可以写为:  - Z5 E6 s! j, N4 t( u- z3 c
初期的PPP技术估算PWV只能采用国际GNSS服务(international GNSS service,IGS)发布的精密卫星产品进行事后解算。为满足实时应用需求,IGS发布了PPP实时服务,许多学者针对基于实时PPP的PWV反演进行了试验验证与评估 [51-52]。随着美国GPS、俄罗斯GLONASS、欧盟Galileo和我国BDS-3星座的发展与升级,相关学者围绕多系统GNSS水汽监测开展了研究,结果表明采用多系统能够有效克服单系统解算稳定性的问题,提升水汽精度[53-55]。 3 d( O0 m, f" I2 x. F ]& V: w! ^
05 - @" C6 P/ a; U
潮位变化监测
. O. x$ h, i8 e* C1 { 潮位变化信息是海洋观测系统和气候变化研究的关键要素。近年来,利用GNSS-IR新型技术遥感潮位变化已成为GNSS应用领域新的研究热点之一。图5给出了GNSS-IR反演潮位原理[56],该技术通过对岸基GNSS台站SNR观测数据进行谱分析,可以估算出潮位变化信息,其算法详细介绍如下。 ! U q$ x l2 Q1 q7 i! I2 H- c1 a4 W
图5 GNSS-IR反演潮位原理
8 J* D( h, @. j- ], m$ @4 q& S3 [9 P5 j3 t; }& r
Fig.5   rinciple behind tidal retrieval based on GNSS-IR technique . k, F9 g. F" r7 B
在实际GNSS-IR数据处理时,仅能采用GNSS卫星低高度角SNR观测值,导致潮位监测的时间分辨率和估算精度受限。鉴于此,有学者通过多星GNSS融合解算或采用机器学习算法克服此不足,使得分辨率和精度得到提高 [57-59]。
* L6 {$ `2 t% a" m+ U( f 06
# @% P" l- Q# u! |+ U$ f 海冰探测
/ [6 P& J z2 T% e2 m) \ 海冰探测对冰情评估、气候变化研究、气候预测和海上交通决策等具有重要意义。近些年, 使用GNSS测量技术在海冰探测方面的研究取得了显著的成果。海冰探测的基本原理是利用GNSS 传输信号的反射入射特性获取海冰特性参数[60],GNSS卫星发射的直射信号极化方式通常为右旋圆极化,当GNSS卫星直射信号入射到海面时,经过海面发生反射,在较高的高度角时反射回来的信号极化方式主要为左旋圆极化。在空气、海冰和海水的界面处,电磁波的反射与入射的能量关系由菲涅尔反射系数确定。根据海水和海冰的介电常数不同,GPS 卫星信号经过不同厚度的海冰时反射信号的右旋分量和高度角的变化存在一定的关系,由此可得到海冰厚度。 7 `- q6 s9 m( k3 U# m8 ~
早在1998年,研究者通过机载GPS反射信号得到了向前散射的GPS反射与RADARSAT的向后散射测量相关性是一致的结果, 根据发射信号特征的差异得到了海冰特征参数,在北冰洋西北部巴罗地区进行机载监测海冰实验,表明GPS信号可以提供海冰信息[61]。2003年, Wiehl等[62]建立了海冰GPS反射信号模型, 利用模型及冰表面反射的GPS信号,提取海冰内部结构及结冰过程等信息。张云等[63]利用GNSS信号的极化模型得到幅度极化比值的变化与海冰的形成和消融过程中海冰的密集程度之间的关系。文献[64-65]分别用GNSS反射信号实现了海冰厚度和海冰边界的探测,表明海冰探测已经到了技术稳定的阶段。 7 v3 J. h+ x- b2 F& \0 R
07 8 {& n$ a+ `1 e5 K4 |3 r! B+ ^
总结与展望
; J0 v$ R6 ^, J e0 c0 |1 b 围绕GNSS在海面环境要素遥感中的应用,本文综述了基于GNSS直射/反射信号的有效波高估算、海面风速反演、降雨强度监测、大气水汽反演以及潮位监测研究现状与进展。总结与展望如下: 9 h" v! S$ \ k) Z' A! z5 m- V
基于GNSS-R/IR的有效波高估算、海面风速反演、降雨强度监测和潮位变化监测技术较为成熟,取得了较好的反演精度。但是现有研究大多基于岸基和星基GNSS平台,后续针对海洋资料浮标、船舶等海洋观测平台上利用GNSS进行海洋环境信息探测有望成为研究热点,将突破性地拓展传统的业务化现场观测平台的观测能力,将在海洋环境业务化监测中具有广泛的应用前景。此外,联合多GNSS星座和多频率观测值提升反演分辨率和性能也应重点关注。由于利用机器学习的GNSS海面环境信息反演结果会受到训练数据样本分布的影响,因此如何进一步提升反演模型泛化能力也将是未来研究的重点。
8 U" |. F9 A4 c/ I 当前基于实时精密单点定位的海面水汽反演已较为成熟,且在短临降雨预警及台风等领域进行了大量研究探索。但在海洋中实现测站空间加密困难,导致海面水汽信息的时空分辨率受限。因此,未来应联合多种水汽探测技术,解决高失控分辨率水汽获取问题。此外,GNSS水汽信息在海洋领域的应用场景依然偏少,如何进一步拓展水汽应用场景至关重要。 # o- t% f4 K' f8 p& V
此外,我国的北斗三号全球系统正式开通,除提供全球导航定位和授时服务外,还为用户提供短报文通信和精密单点定位增强服务等多种特色服务。因此如何充分利用北斗三号的优势,挖掘北斗系统在海面环境信息遥感中的应用将是未来科学的重要工作。 5 y, F2 C5 a/ L; P5 T1 }
08
4 f7 Q" Z. I U9 x$ F: I6 v 参考文献
0 `; g0 `1 K; M) \% e 略,详见全文:http://www.sdkx.net/article/2024/1002-4026/1002-4026-37-2-1.shtml I/ Y& M7 J. K+ F
山东科学
$ [4 U0 R; b' J) \. T' p H) F+ ` 《山东科学》是齐鲁工业大学(山东省科学院)主管/主办的自然科学综合性学术期刊,现为中国科技核心期刊、RCCSE核心期刊、中国核心期刊(遴选)、中国高校优秀科技期刊。 . E. q7 t& m# N* F4 v' `
主管/主办单位:齐鲁工业大学(山东省科学院)
! w; F; C" U" f 地址:山东省济南市历下区科院路19号 6 ~$ P# ?4 R/ y9 W" u# h, \
电话:0531-82605310
w! u- Q: ~1 g' W' f 官方主页:http://www.sdkx.net ! ~+ z u. N, b
E-mail:xuxy@qlu.edu.cn # g# C& R+ r' D( \: n# W5 e7 ~$ u
山东科学作者交流QQ群:151903851 % }0 f: W0 ~. x3 @1 K* ?1 Z: f# v
2 p; y, }1 W& V# O' D7 D
4 w8 l* o- }3 G9 e0 b7 P$ y) F
) g8 c6 h4 k R% N
' B/ c# ~- k" g |