地基GNSS大气水汽探测遥感研究进展和展望

张克非，李浩博，王晓明，朱丹彤，何琦敏，李龙江，胡安东，郑南山，李怀展

1. 中国矿业大学环境与测绘学院，江苏 徐州 221116； 2. 清华大学水利水电工程系，北京 100084； 3. 中国科学院空天信息创新研究院，北京 100094； 4. 苏州科技大学地理科学与测绘工程学院，江苏 苏州 215009； 5. 科罗拉多大学博尔德分校环境科学联合研究中心，科罗拉多 博尔德 CO80309； 6. 美国国家海洋和大气管理局空间天气预报中心，科罗拉多 博尔德 CO80309

随着全球气候变暖，极端天气事件增多增强趋势明显，诸如极端强降水、台风、洪涝等突发性强、局地性明显的气象灾害事件对人类社会可持续发展构成了重大威胁。作为表征极端天气事件和气候变化的重要参数，大气水汽在极端天气事件的发展、发生过程中占据着重要的作用。因此，对大气水汽含量进行有效监测和科学分析是实现各类灾害性天气事件精准预报的重要前提与关键所在，也是准确预报极端天气事件、提升防灾减灾应急能力的重要保障。

大气水汽主要分布在对流层内，其复杂的空间分布及快速变化的特性使得大气水汽含量的高精度监测极其困难。目前，大气水汽监测的主要方法有探空气球、大气水汽辐射计、遥感卫星、GNSS等。探空气球利用搭载在氢气球上的探空仪器对大气状态进行监测，能够获得0～36 km不同高度的气温、气压、湿度等多种气象参数，是近百年来最主要大气探测手段，但是探空站数量有限、分布不均匀，且存在时间分辨率低、运行成本高和气球漂移等重大缺陷。水汽辐射计通过测量大气中的水汽辐射量测定水汽及液态水含量，能够提供高精度、高时间分辨率的水汽数据，然而，由于该仪器的造价成本比较昂贵(约150万/台以上)，在实际中难以推广使用。遥感卫星利用搭载在卫星上的传感器接受地表反射的电磁波，基于不同波段的电磁场辐射信息进行大气水汽反演，能够提供高空间连续性的广域(特别是无人区、洋区和极地)水汽数据，但所得结果易受天气条件影响较大，且精度较低。上述水汽监测手段虽然能够提供较为可靠的水汽数据，但是受限于数据处理的时效性，通常难以满足短临极端天气预报的实时性需求。随着新一代GNSS的高速发展，利用GNSS观测值进行大气水汽实时反演能够有效弥补传统水汽探测手段的不足，同时还具有高精度、全球覆盖、全天候、长时序稳定及高时空分辨率等优点，为开展大气水汽监测提供了重要数据支持。根据观测模式的不同，GNSS大气水汽探测技术可以分为地基GNSS和GNSS掩星技术。本文主要针对地基GNSS大气水汽探测遥感技术的研究现状与发展趋势进行概述和总结。

GNSS技术的快速发展给地基GNSS大气水汽探测技术提供了重要的技术支撑，自20世纪70年代美国提出建立全球定位系统(GPS)以来，经历了不断地发展与完善，于1995年达到全面运行能力，能够为全球范围内的用户提供全天候以及高效能的导航定位服务。GPS卫星信号穿过对流层时，受到中性大气的影响会出现路径弯曲和速度衰减，从而使GPS观测值中含有一定的大气延迟，基于对GPS信号大气延迟量的高精度观测便能实现大气水汽含量的监测。自1990年代开始GPS大气反演理论提出以来，开展高精度GPS大气水汽反演及其相关气象应用研究受到了国际学者们的广泛关注。随着我国北斗卫星导航系统(BeiDou navigation satellite system，BDS)的全面建成，以GPS、BDS、GLONASS、Galileo为核心的新一代GNSS观测体系为大气水汽的高精度探测和气象应用奠定了重要的空间设施基础。同时，各级卫星连续运行参考站(continuously operating reference system，CORS)网络的大规模建设，也为GNSS大气水汽探测遥感科学研究和气象创新应用提供了重要的地面观测设施。总体来说，地基GNSS大气水汽探测技术集成了全天候、无漂移、高精度、高时空分辨率及低成本等优点，可较好地与传统大气水汽监测技术形成互补，为大气水汽监测和短临天气预报提供高质量的数据源和基准信息，地基GNSS大气水汽探测技术在大气探测以及气象领域应用中的重要地位已逐步凸显。

1 多类型GNSS大气参数高精度反演

地基GNSS技术不仅可以提供天顶对流层延迟和大气可降水量等二维水汽产品，还可以提供斜路径对流层延迟、斜路径水汽含量、水平梯度和三维水汽密度等三维廓线信息。目前，地基GNSS大气水汽研究主要集中于高精度ZTD/PWV反演和三维水汽层析建模。

1.1 ZTD/PWV反演研究进展

在地基GNSS发展初期，大气延迟作为GNSS定位的主要误差源，通常采用经验模型或气象观测值等方式予以改正或削弱，后来作为独立的未知参数联合求解。文献[1]于1992年首次利用GPS观测数据反演得到接收机上空的PWV，进而形成GNSS气象学这一全新的学科领域。20世纪90年代至今，国内外学者开展了大量试验对GNSS反演ZTD和PWV产品的精度和适应性进行评估，使得GNSS气象学得到了迅猛发展，地基GNSS数据处理主要有双差网解法和精密单点定位法。在双差网解法进行大气水汽反演方面，文献[2—6]先后在不同地区开展了地基GNSS水汽反演试验，并与水汽辐射计等不同水汽产品进行了比较。文献[7—12]分析了GNSS水汽反演的主要误差源，利用探空数据证明了不同天气场景下GNSS水汽产品的精度可达1～2 mm。虽然双差网解法的大气水汽反演精度较高，但是在进行大规模GNSS数据处理时，双差网解法具有较高的复杂性，随着GNSS精密单点定位技术(precise point positioning，PPP)的出现，受其模型简单、处理高效的优势引导，文献[13—15]围绕PPP反演ZTD/PWV开展了大量试验。但是，精密单点定位技术受卫星产品精度的影响较大，在早期难以满足实时水汽产品精度需求。为了克服该问题，2013年国际GNSS服务(international GNSS service，IGS)正式向全球用户播发GNSS实时数据流产品，此后文献[17—22]对基于GNSS实时PPP的ZTD解算和PWV反演精度进行了评估，大量试验证实地基GNSS能够获取高精度的实时PWV产品。此外，事后PWV产品具有高精度、高时空分辨率的特点，鉴于此，文献[23—25]使用GNSS-PWV产品作为定标改善卫星遥感和数值预报模型产品获取的PWV系统偏差。同时，有学者还研究了利用低成本单频GNSS接收机反演水汽的可行性。例如文献[26—27]利用经济型的单频接收机或智能手机代替大地型接收机反演水汽，其结果表明利用低成本的GNSS接收机能够有效获取对流层水汽。由于海洋区域缺乏密集、稳固的GNSS安装平台，文献[28—30]使用船载或浮标GNSS水汽反演技术来获取海洋上空的水汽，为动态GNSS水汽反演理论方法开展了探索性的研究。文献[31—32]利用星载辐射计获取了全球范围内公海区域上空的PWV，并将其结果与多次船载GNSS-PWV进行了对比，结果表明在极地地区均方根最小(1.29 mm)，其他地区为2.78 mm；与NECP再分析产品比较，两者差异为2.96 mm。文献[33—34]利用无线电探空仪、地基GNSS、掩星数据和卫星遥感获取了在海洋区域内水汽时空特征，结果表明上述几种方式的PWV结果具有很好的一致性。此外，中国气象局气象探测中心联合内蒙古自治区与广东省气象局于2017年5月联合起草了国内首部地基导航卫星遥感水汽观测规范[35]，规定了导航卫星遥感观测、值班和数据中心运行监控等要求，为地基GNSS水汽产品的规范化监测与使用奠定了基础。

1.2 GNSS三维水汽建模研究进展

地基GNSS大气水汽的初级产品是天顶方向上的大气水汽总量，无法直接给出信号传播路径上的水汽分布信息，在天气监测和气象预报等领域存在一定局限性。为拓展地基GNSS技术在气象预报中的应用场景，文献[36—40]研究了信号传播路径上的大气水汽反演方法，斜路径水汽(slant water vapor，SWV)精度评估结果表明对流天气条件下的SWV的精度可优于4 mm。然而，SWV仍然无法全面解析目标区域垂直方向上的水汽分布，为了获得三维水汽场信息，文献[41]首次提出了基于GNSS观测值获取对流层水汽结构的设想，文献[42]首次利用层析技术获取了夏威夷基拉韦厄火山上空的湿折射率，并基于欧洲中期天气预报中心(European Center for Medium Range Weather Forecast，ECMWF)的再分析产品对反演结果进行精度评估。随后，文献[43—45]利用三维层析技术成功获取了不同地区的水汽三维分布，并利用探空站和再分析资料验证三为水汽层析的可行性。

目前三维水汽层析模型研究的难点在于无法解决层析方程系数矩阵的病态问题，其主要原因是卫星和地面GNSS站的几何分布不均匀使得卫星信号无法完全覆盖研究区域[46-47]。因此，解决层析方程中的矩阵病态问题是获取高精度水汽三维分布的关键所在，目前的主要处理方法有[48]：①附加约束条件，文献[49]基于水汽垂直变化特征采用高斯指数模型作为约束条件；文献[50]利用探空信息建立了顾及边界信号及垂直约束的方法；文献[51]提出了附加高水平分辨率PWV约束的层析方法。②增加观测值法，文献[52—53]使用探空历史数据建立了优化体素分割方法，从而提高了研究区域测边信号的利用率，研究结果表明射线利用率和有效体素数分别提高了30.32%和12.62%；文献[54]提出了一种基于截断因子模型的对流层层析算法，该方法能够充分利用层析区域内侧边信号的水汽观测值；文献[55]引入中分辨率成像光谱仪(moderate-resolution imaging spectroradiometer，MODIS)的水汽观测值提出了一种融合GNSS、MODIS观测值的水汽层析方法；文献[56]提出了融合GNSS及合成孔径雷达观测值的水汽层析方法。③利用模型改进层析精度，文献[57]利用神经网络技术反演了研究区域内的三维水汽分布；文献[58]提出了一种全新的节点参数化方法(图1)，该方法可以在研究时段内动态调整层析区域的边界和节点位置；文献[59]基于水汽与高度之间的指数变化关系提出了一种自适应非均匀分层层析方法。

图1 基于GNSS信号分布的自适应层析模型[58]

1.3 GNSS水汽反演未来研究及应用展望

当前地基GNSS水汽反演算法已基本成熟，随着GNSS系统的完善和地面监测站的增加，地基GNSS水汽产品已满足气象研究的高时空分辨率、高精度需求。除了GNSS大气水汽产品外，气象卫星、气象站、探空站、雷达、再分析资料等多源数据均可提供高精度大气水汽产品。强化推进测绘、气象等资源交换共享机制，赋能“空-天-地”一体化水汽监测网络体系，提供综合化、多样化和精细化水汽产品是气象业务广泛化的必然要求，当前仍有一些问题尚需研究。

(1) 提高泛在水汽反演能力。低成本GNSS接收机的推广使用可极大节约经济成本，智能时代下基于低成本GNSS接收机的泛在水汽反演已经成为可能，但是，低成本接收机的精度、可靠性在复杂环境、恶劣天气下仍存在一定问题，由此导致的水汽反演精度仍需进一步提高。因此，在未来研究中需要建立基于空间变异的无损水汽反演方法，进一步提高基于低成本设备的大气水汽反演性能，以满足气象、农业、交通等多部门的业务需求。

(2) 提高多源水汽综合利用能力。全球范围内探空气球、水汽辐射仪、数值预报模式等能够提供不同空间分辨率的水汽资料，融合基于多源水汽进行水汽层析联合解算，可有效改善区域GNSS水汽监测网络低密集度的问题，提高水汽层析反演精度。但是，当前研究未充分考虑多源异构水汽产品的时空变异性特征，因此，研究多源水汽观测数据归一化处理方法，解决多源水汽时空不一致等问题，建立基于动态大气背景场的时空水汽融合模型，生成精细化水汽产品是拓展大气水汽气象业务应用场景的必要前提。

2 极端气象事件短临预报应用

极端天气事件是在某个特定地点和具体时段内所发生的一种极其罕见的气象事件。联合国政府间气候变化专门委员会(Intergovernmental Panel on Climate Change，IPCC)的研究报告明确指出，在人类社会快速发展与全球气候变化的大背景下，强降水、热带气旋、洪涝、干旱、寒潮、对流风暴等典型极端气象将向高频率、复合型灾害事件转变，将会严重影响了人们的生命财产和正常的生产生活。随着我国现代化工农业以及航空等诸多行业的快速发展，对天气预报(特别是极端气象事件短临预报)的准确度以及时效性等均提出了更高的要求，天气预报能力亟须进一步发展。地基GNSS大气探测技术能够进行大气水汽的实时监测，可实现面向各类极端天气事件的短时(0～12 h)和临近(0～6 h)预报[60-65]，是改善当前极端天气突发事件预报能力不足的重要手段。此外，GNSS观测站已积累了近30 a的观测数据，为极端天气事件发生的机理及演变规律研究提供了极为重要的数据源，对极端天气的机理研究可进一步提高短临预报预警能力，对减少气象灾害风险、提升防灾减灾应急能力具有重大科学和现实意义。下文将对利用多类型地基GNSS大气观测开展极端气象事件短临预报应用的研究进展及未来展望进行介绍。

2.1 面向气象应用的数据预处理及系统分析

针对面向气象应用的地基GNSS数据预处理，首要是综合利用无线探空仪、水汽辐射计等大气观测数据和甚长基线干涉测量(very long baseline interferometry，VLBI)、IGS等事后产品，分析地基GNSS大气水汽反演过程中的误差特性，严格控制地基GNSS水汽产品的精度，以满足极端天气短临预警预报应用研究的精度要求。表1和表2以PWV为例，分别给出了面向极端天气短临预报以及长期气候分析应用的产品要求[64]。

表1 面向极端天气短临预报应用的PWV产品要求

表2 面向气候分析应用的PWV产品要求

需要指出的是，不同气象应用对输入数据的精度和时空分辨率的要求有所不同。对于天气预报数值(numerical weather predication model,NWP)模式为基础的短临预报，由于NWP自身的非线性较为复杂，对观测值精度的要求较高，具体来说，NWP模型的观测值需满足无偏性，相应的观测误差应是不具有空间互相关性的高斯噪声。非NWP模型下的短临预报更为关注大气参数的相对变化情况，微小误差并不会对预报结果产生较大影响，较之于非NWP模式，非NWP模型对数据绝对精度要求较低。在面向气候分析的应用中，通常需要使用较长时间跨度的数据进行研究，因此长时间序列中存在的强度较小的误差也同样不会影响数据的整体分析结果。

除了对输入数据本身的精度有所要求以外，对多气象参数的系统性分析也是面向极端天气预报中不可或缺的一步。首先，可利用皮尔逊相关系数、斯皮尔曼相关系数、偏相关系数、复相关系数等作为评价因子计算多类型大气参数间的互相关系数从而分析各参数间的相关关系[66-67]；同时，可采用主成分分析、经验正交函数等方法探讨各参数对典型极端气象事件的响应关系及它们在该事件短临预报应用中的有效性和贡献度，从而为各类模型预报因子的选择提供了有效的参考依据。

由于强降水、台风、干旱等事件的发生和演变过程均与GNSS大气水汽信息的时序变化高度相关，故GNSS大气信息广泛地应用于这类事件的预报研究[67-69]。文献[70]分析了阿拉伯海和孟加拉湾逐月及季节性PWV变化同季风的关系，发现了季风的发生同PWV的变化紧密相关，并提出可以使用PWV来预报季风的发生[71]。文献[72]研究了GNSS大气水汽在如暴雨、持续性降水、冰冻雨雪等典型天气中的演变特征，有效揭示了GNSS大气探测技术及产品在天气预报业务中的应用。以暴雨事件为例，文献[73]发现暴雨事件发生前PWV会以9.5 mm/h的速度上升，并将该表现称为PWV的跳跃性变化。文献[74]发现在长时序的PWV中，其快速上升的时序变化特征能有效揭示降水事件的发生。文献[75]则从另一个侧面发现在热带风暴引起的暴雨事件发生前的3 h内，PWV时序数据会出现快速下降。类似地，文献[76]发现在暴雨事件发生前，PWV往往会存在一个以5～10 mm/h速度的下降过程。此外，文献[77]的研究也表明大气水汽在极端降水事件发生前快速增加的现象能准确地在GNSS水汽产品上得以体现。然而由于极端天气的发生具有极为复杂的成因，因此要利用GNSS大气产品进行极端天气短临预报就需要分析各类产品在大量极端天气发生过程中的时序变化特征，再结合其他气象观测信息(如雷达、辐射仪、实时气象观测等)，从而据此利用多种方法建立区域性的基于GNSS大气信息的极端气象事件短临预报模型。

2.2 极端气象事件短临预报的主要方法

上述研究充分表明GNSS大气水汽同极端天气等气象事件具有较强的相关性，通过利用高精度、高时空分辨率的ZTD、PWV等大气资料能够有效地开展气象事件的预报研究。目前，虽然国内外利用GNSS大气反演资料进行极端气象事件的短临预报均进行了一些初步的研究，但尚未形成成熟完备的理论与方法。在当前研究中，所使用的各类极端天气短临预报模型主要包括以下几类：

(1) 气象预报专家系统模型。该类模型主要依赖于所使用的气象预报因子同特定气象事件间的响应关系而构造，最为典型的气象预报专家系统模型是基于大气参数数值异常变化的阈值模型，即当某类极端气象事件发生前，一种或多种大气参数的数值会出现异常增大(或减小)的情况，根据该变化趋势设定预报阈值可开展短临预警预报。在预报阈值的优化构造选取过程中，为了保障所选预报阈值的可行性和普适性，通常需要利用该因子的长时间历史数据以及对应的气象事件观测记录进行统计分析。近年来基于该类模型也形成了面向不同极端天气事件的多种短临预报模型[78-84]。文献[85]通过研究西班牙9年内的GNSS水汽数据同降水间的关系提出了降水预测模型。文献[86]提出了一种结合可降水量季节变化信息的降水预测算法，基于过去30 min的GNSS水汽信息预测未来5 min降水事件的发生概率，结果表明该算法有效提高了探测的准确率。文献[87]基于大气可降水量时变特征与降水之间的响应关系，提出了一种利用GNSS大气可降水量时序变化特征开展短临降雨预报的方法，测试试验表明该方法取得了较高的预测精度和较好的适用性。文献[88]分析了西班牙里斯本地区2010—2012年暴雨事件发生前PWV的时间变化特征，并提出了一种基于6 h时间窗口的降水预报模型，该模型利用PWV最大增长率作为预报因子，经验证其预报正确率达75%，而误报率为60%～70%。基于该研究，文献[89]提出了一个基于PWV、PWV变化量、PWV变化率的三因子降水预报模型，该模型的正确预报率达到了82%，但模型误报率仍为60%左右。文献[90]基于热带地区(新加坡)的数据对于上述方法进行检验，并提出了适用于热带地区的以最大PWV值作为主要评价指标的预报方法，该模型的误报率较上述研究降低了17%。文献[91]提出了一种改进的降水预报模型，该模型利用包括ZTD以及PWV在内的五因子进行预报，结果表明该模型的正确预报率可以达到95%，并且模型误报率降为30%左右。文献[92]首次提出了一种在极端降水发生前顾及包括PWV增长及下降完整时序变化过程的模型，模型的预报性能得到了进一步的提升。除了上述研究外，许多研究还提出利用其他多类型的GNSS观测[66,93-94]、三维/四维层析技术[42,95-98]、更长的PWV时序、更高时空分辨率的时序数据来开展短时降水预报。虽然上述所构造的各类模型方法原理简单且易于操作，具有较强的可行性，但往往也会存在稳定性低、计算量大、可移植性差等缺点，因此该方法通常主要面向某一小范围区域(局域)气象事件的短临预警预报应用。

图2 基于PWV距平的降雨预报示例[79]

(2) 数值天气预报模型。NWP模式是天气预报应用中不可或缺的一种工具，由于其具有完善的大气物理模型，且能够通过数据同化的方式利用多类型的大气参数，因而能够有效加强对于天气和气候系统演变规律的认识，获取高性能的预报表现。NWP模式对于GNSS大气反演信息的利用主要是基于数据同化的方式，常见的是利用变分同化方法开展针对PWV、ZTD、斜路径延迟(slant total delay，STD)及三维层析产品等非常规观测资料的同化应用。因此，为了进一步增加各类模型中使用的气象参数且不断完善模型的物理基础，文献[99—100]对GNSS大气反演资料在NWP模式中的同化开展了一系列卓有成效的研究，目前诸多研究都证明将GNSS大气产品同化到NWP模式中能有效改善模式初始场的精度，进而提升模型对气象事件的预报能力，特别是对极端降雨等突发灾害性天气。文献[101—102]在假定大气湿度廓线不变的条件下，首次应用模拟的PWV资料进行了同化试验，结果表明，同化PWV能有效提高降水预报精度，并进一步验证了利用四维变分资料同化技术同化PWV的可行性和有效性。文献[103]利用梅雨季节长三角地区获取的GNSS大气水汽资料，分析出入梅(雨)过程水汽的空间分布特征，并通过对其的同化有效提高了数值模式预报能力。文献[104—105]进行了ZTD的四维变分同化试验，结果表明其对6 h和12 h的累计降水预报精度分别提升33.15%和25.08%。文献[106]针对2001年地中海西部极端降水事件开展了ZTD同化试验，结果表明同化ZTD比同化PWV得到了更好的预报结果，能够有效降低强对流天气的预报误差。然而由于PWV和ZTD只能反映GNSS测站上方大气水汽的整体变化信息，并不能反映其在水平方位上的梯度变化信息。因此，如何有效同化GNSS反演得到的斜路径产品及三维层析产品成了目前该领域最为前沿的研究热点，并且诸多研究也显示同化STD、三维层析产品比PWV和ZTD在改善大气背景场精度及降水等气象事件的预报结果上有更积极的影响[107-109]。文献[110]研究发现通过同化STD能得到更好的水汽三维分布及降水日变化信息，从而有效改善了实时与短期定量降水预报。文献[111]开展了STD及多普勒雷达观测资料的四维变分同化试验，结果表明相较于同化PWV模式面向降水事件预报取得了更好的效果。文献[112]将三维水汽层析产品同化进数值天气预报模型中，结果显示模式预报的相对湿度及温度等参数的精度均得到了明显提升。文献[113]利用TOMOREF观测算子开展暴雨事件同化试验，结果同样表明通过同化三维层析产品，无论是降水还是湿度场的预报效果都有了显著改善。大量同化试验研究结果均表明，加入GNSS资料对改善数值天气预报模式初始场精度有着较为明显的作用，可有效提高极端气象事件预报精度[84,100,114-124]。然而由于该类模型对于各类参数同化过程的时效性，特别是在处理大规模区域上的多类型观测数据时，因此难以有效地满足0～6 h的临近预报的应用需求。

(3) 基于神经网络技术的模型。目前神经网络技术在地球科学领域中的应用越来越受到人们的关注，由于基于该技术开发的各类模型具有独立学习、并行处理、非线性映射能力、容错能力和自适应能力强等优点，因而为研究多种气象信息间内在关系、开发基于神经网络技术的极端气象事件预报模型提供了前所未有的契机。虽然该类模型同样需要使用大气参数的长时序数据学习拟合出参数同目标气象事件间的耦合关系，然而与阈值模型不同的是，通过利用神经网络技术能够优化使用多种类型的气象参数，通过利用多类型大气参数中蕴含的大气环境变化信息更为全面地揭示极端气象事件的形成与演变过程。近年来，随着各类神经网络模型的逐步成熟以及神经网络集成、遗传神经网络和模糊神经网络等技术的兴起，文献[125—131]通过利用神经网络技术开展了大量的短临天气预报的应用研究。文献[132]采用误差反向传播神经网络建立了汛期降水预测模型，结果表明该模型无论是对历史样本的拟合还是对独立样本的预报检验都取得了较高的精度。文献[133]建立了自适应模糊神经网络天气预报系统模型。文献[134]采用人工神经网络进行短期降水预报，提高了预报的准确率。文献[135]釆用集成预报路线，将低级单模式数值预报输入模型，输出得到统一的集成预报结果。文献[136]在基于卫星遥感图像降雨监测系统中基于遗传优化反向传播神经网络(back propagation neural network，BPNN)对研究区域内云团移动和降水量的预测进行了分析研究。文献[137]提出了一种利用包括PWV在内的七参数的神经网络模型用于开展降水事件预报，结果表明其准确预报率为80.4%，误报率为20.3%。文献[138]提出利用多种气象参数、时间相关参数及可降水量信息，通过多隐层神经网络建立降雨和各种参数的非线性关系对短临降雨进行预测，结果表明，该模型可预测出95%以上的降雨事件。文献[139]基于神经网络对数值预报模型数据进行了处理，并改善了在冰雹天气预报的准确度。文献[140]同样基于BPNN算法，提出了一种改进的极端降水短临预报模型，该模型大大降低了误报率，取得了较好的预报效果。由于GNSS反演的长时序大气资料具有实测独立的特点，通过利用神经网络技术，基于其强大的自主学习及非线性映射能力开展数据分析，能够有效地挖掘多类型大气参数中蕴含的信息，构建针对极端气象事件的短临预警预报方法。由于多种类型的神经网络模型不断被提出，针对特定的气象事件以及所使用的气象参数类型还需要进一步开展分析，从而选择更为合适的模型类型开展短临预报预警应用。此外，模型构造过程中网络结构、预报参数、模型超参数等各类因子的优化选取、对预报模型内部机理的理解以及模型拟合过程的优化也都是当前仍需进一步探究解决的问题。

2.3 GNSS极端天气短临预报未来研究及应用展望

GNSS气象学作为一个跨学科研究方向，通过综合利用GNSS反演的大气资料对提高各类气象模型预报性能，以及开展针对极端气象事件的短临预报应用有着非常积极的作用。当前的研究尽管在多类型GNSS大气信息精度检验、多源大气参数分析、极端气象事件形成与演变过程探究等方面进行了诸多研究，也利用包括GNSS大气产品及常规气象观测在内的多类型大气参数优化构建了各类预警预报模型，充分提升了面向极端气象事件的短临预报性能，但如何充分挖掘海量GNSS观测数据中蕴含的极端气象事件短临预报的前兆信息、如何有效地因地制宜、因时制宜地构造极端气象事件短临预报模型仍存在诸多亟待解决的问题。

(1) 针对气象预报专家系统模型应用。众所周知，各类极端气象事件形成及演变的天气学过程极其复杂。以极端降水事件为例，它不仅需要充足的水汽，还要满足垂直运动、云滴增长等一系列条件，然而当前气象预报专家系统模型中所使用的预报因子仅能有效反映大气中水汽及湿度的相关信息。虽然利用各类衍生变量，也能在一定程度上表征大气动力学变化过程，但从业务预报应用的角度仍远远不够，因此应进一步加强对大气环境动态变化及天气系统形成过程机理的研究，发展实时GNSS水汽资料应用于极端天气短临预报/预警的方法与模型。

(2) 针对数值天气预报模型。虽然当前业务使用的各类数值天气预报模型能用数据同化的手段综合相关信息来提供高时空分辨率的气象产品，也能以较高的精度、较广的覆盖范围提供未来一定时段内的气象预报结果，然而目前在各类业务数值模式中均未优化实现对多类型GNSS大气信息的同化应用，这是该领域目前尚未解决的关键技术。由于不同类型的GNSS大气产品具有不同的误差特性，为了更好地利用多类型GNSS大气产品，应进一步探索构建面向不同探测资料的同化算法，在充分提升模型预报精度、精准确定事件落区的同时，进一步提升面向大规模站网数据的同化效率，并充分延长各类极端天气的预见期，从而有效满足极端天气事件防灾减灾决策响应需求。

(3) 针对基于神经网络技术构造的模型。在模型构造过程中的网络结构确定、权值设定、超参数选取等各类问题目前仍需进一步探究，从而形成一套优化的针对极端气象事件短临预报应用的模型构造方案以增强模型的泛化能力；在完成模型构造过程后，应进一步增强对模型构造机理、内在算法实现等方面的分析解读能力；探索各类成熟的神经网络模型算法在极端天气短临预报应用中的可行性与适用性，从而进一步深化神经网络技术在GNSS气象学领域的应用。

3 气候变化分析应用

大气水汽能够吸收来自地球的红外辐射，贡献了全球60%的温室效应，是大气中最重要的温室气体之一，同时大气水汽与温度之间存在着明显的正相关，这种正相关作用使得大气水汽含量在全球变暖的趋势下不断增加，从而进一步加剧了大气温度变化，严重影响平衡气统敏感性。大气水汽作为大气状况的主要表征之一，分析其长期变化趋势，研究大气水汽与多种气象参数的耦合关系，进而研究大气状态的变化趋势，对深入理解全球气候变化、灾害性天气监测等方面具有十分重要的科学意义。

3.1 PWV时间序列一致性检验方法

随着空间观测技术的完善和数据处理技术的发展，目前已累计了相当长度的地基GNSS-PWV时间序列，可用于全球或区域性的气象分析[141-145]。一致化的基准数据是进行高精度气象研究的前提，但是地基GNSS不可避免地受到如接收机更新、硬件升级、站点环境变化、处理策略更新等多种非气象因素的影响，导致PWV时间序列中含有非气象因素引起的、相对于自然变率无法忽视的系统性偏差，即非一致现象[145-148]。为了减少非气象因素的影响，GNSS-PWV的质量控制和一致化检验是利用GNSS-PWV时间序列进气象研究必不可少的基础。从时空上看，非一致现象包括空间非一致性和时间非一致性两个维度，而PWV时间序列的非一致性主要集中在时间维度。

根据非一致现象产生的原因，非一致性主要分为与接收机有关的、与站点环境有关、与数据处理策略有关的。与接收机有关的非一致现象主要包括仪器更新、天线罩配置、固件升级、地质活动等，这种非一致现象往往表现为某一历元前后PWV时间序列出现跳跃式变化[149-153]。站点日志中记录了大部分接收机硬件相关的变化，因此站点日志文件往往用于非一致性检验中。由于站点日志记录的不完整和误差等原因，站点日志无法实现硬件变化的完整记录[154-156]，且日志文件中不含有与地质活动有关的站点变化[157]。与站点环境有关的非一致现象主要源于站点植被生长、积雪覆盖等因素，相较于硬件更新导致的阶跃点，这种非一致现象在PWV时间序列往往表现为缓慢变化，其大小和位置均不明确，是一致性检验的难点[150,151]。处理策略变化包括国际地球参考框架(international terrestrial reference frame，ITRF)变化、天线相位改正模型、误差改正模型、高度角等，这种非一致现象较为简单，在数据处理时只需采用统一的处理策略即可消除[148,158-160]。PWV时间序列的一致性检验分为非一致性探测和修复，其中前者是一致性检验的重点和难点。针对大气水汽时间序列的非一致性探测，国内外学者展开了研究，提出了一系列基于统计检验的非一致性探测方法，大体可分为基于参考数据的相对探测方法和基于原始数据的绝对探测方法。

相对探测方法是当前使用最多的非一致性探测方法。该方法首先利用参考数据与原始目标数据的空间相关性，利用差分算子消除目标数据的气象信号、周期变化和观测噪声，基于差分序列完成非一致性的探测[147,161-162]。文献[163]基于双相回归模型(two-phase regression，TPR)构造似然比统计值，并利用F检验对温度差分时间序列进行一致性检验。文献[156]在TPR模型的基础上，针对传统T检验和F检验探测结果虚警率较高的缺点，提供了基于经验惩罚函数的惩罚最大T检验(penalized maximal T test，PMT)和惩罚最大F检验(penalized maximal F test，PMF)[164]，并在此基础上提出针对观测噪声的一阶自相关特性的PMTred法和PMFred法[165]。文献[166]利用PMFred检验对并探空站湿度差分数据进行检验，得到了一致化的探空气象数据。文献[167]以ERA-Interim数据、并址GNSS站和13个VLBI站的PWV观测值为参考数据，利用PMTred方法对全球101个IGS站的PWV数据进行了一致性检验，结果表明47个测站有含有62个阶跃点。文献[168]以BIC准则为基础提出了一种顾及阶跃点前后平均值和方差变化特性的异方差阶跃点探测模型，并采用ECMWF再分析数据集ERA-Interim数据作为参考值研究了数个GNSS站PWV时间序列的一致性。

绝对探测方法采用原始数据的时间变化特性拟合数据本身的周期信号，基于去除周期信号的残余时间序列完成非一致性探测，相较于参考数据的相对探测方法而言，该方法不需要外部数据的参考，因此学者们对该类方法进行了研究和探索。文献[169]利用T检验构造了序贯T检验方法(STARS)方法。文献[154]对时间序列进行分段，并将时序的总方差分为外部方差和内部方差，据此以外部方差最大和内部方差最小为参考量提出了一种多阶跃点探测方法。文献[170]基于奇异谱分解方法构造了一种阶跃点探测非差方法，基于仿真数据的结果验证了该方法的探精度为30 d。文献[171]将时间序列分解算法(seasonal-trend decomposition based on LOESS，STL)方法和PMFred(PMF to account for the red noise)方法进行结合提出一种新的绝对自适应探测方法(adaptive absolute homogeneity test，AAHT)。该方法将原始PWV数据分解为周期项、趋势项和噪声项，利用趋势项完成非一致性探测，并基于噪声项构建自适应因子进行探测结果的虚警点检验，保证探测结果的可靠性。

图3 ALBH站大气水汽时间序列基于AAHT法的非一致性探测结果[171]

3.2 PWV时间序列的特征演变分析

长期趋势项和不同尺度的周期信号是基于PWV长时间序列进行气象分析的重要参考之一，很多学者进行了基于全球或区域PWV数据对其长期变化趋势研究，并通过对水汽变化的机制和主要驱动力进行分析。

针对PWV数据的时间演变特征，文献[172]采用全球探空数据和Spearman相关系数，分析了全球1973—1995年和1958—1995年两个时间段PWV的长期趋势，首次发现大部分GPS站的大气水汽与温度呈显著的正相关。文献[173]分析了不同大气水汽资料指出全球水汽分布受全球气温升高影响在整体上呈上升趋势。文献[174]基于全球1997—2007年的探空数据、微波辐射计和GPS观测值，分析了全球的PWV序列的变化特征。结果表明，PWV序列在整体上呈现出上升的趋势，且海洋地区的变化高于陆地地区，夜间变化高于日间变化，并且白天的PWV变化更加符合热力学过程，而夜间的PWV变化主要受动力学过程影响。文献[175—176]顾及PWV观测噪声的一阶自相关特性，分析了德国的PWV时间序列的长期变化趋势。文献[61]指出2000—2014年全球PWV的平均变化趋势为每年1.29‰，但其结果的显著性较差。文献[177]利用GNSS观测值和再分析资料的月平均PWV数据，分析了PWV的长期趋势变化及其与地表温度的相关性。文献[178]利用ERA5和GNSS-PWV分析了印度洋南部区域的日间和季节变化特征。文献[179]分析了地中海东部PWV不同季节的昼夜差异性。文献[180]利用水汽辐射计、傅里叶变换红外光谱仪(Fourier transform infrared spectrometer，FTIR)和GNSS的PWV信息研究了瑞士地区的水汽趋势变化特征。文献[181]利用中国地壳运动观测网络(crustal movement observation network of China，CMONOC)中28个GPS基准站2004—2007年的观测值分析了中国大气水汽的周期变化特征与地理位置之间的关系。文献[182—183]利用经一致化处理的探空数据估算了中国上空大气水汽的长期变化趋势，同时分析了探空数据的不一致性对大气水汽序列长期趋势带来的影响。文献[184]利用成都30 min的PWV时间序列分析了大气水汽的日变化趋势，指出大气水汽日变化趋势与太阳照射强度存在一定的相关性。文献[185]利用CMONOC的观测值数据研究了中国PWV的时空分布特征和热力学过程，结果表明中国的PWV平均值随纬度增减而降低，最大值和最小值分别出现在华南地区和青藏高原地区，在水汽驱动力方面，PWV的热力学变化与中国的气候分布有关。上述研究的结果表明在全球范围内PWV总体上呈现出上升趋势，强度约为1‰～2‰，但其结果的不确定度较大，置信水平较差。

根据克劳修斯-克拉珀龙(Clausius-Clapeyron，C-C)方程，在相对湿度稳定的条件下，温度每升高1°，大气水汽含量会增加7%，但是PWV-温度趋势变化的相关分析结果与C-C方程式的理论值(7%/K)具有一定的差距，因此学者们针对大气水汽变化的其他驱动力，研究了厄尔尼诺现象与大气水汽变化的关联性。文献[186]在2000年首次研究厄尔尼诺与南方涛动(El Nio-Southern Oscillation，ENSO)期间大气水汽的变化情况。文献[173]利用更长跨度的PWV时间序列分析了ENSO对热带、亚热带地区水汽的影响，证明了文献[186]结论的可行性。文献[148]利用1994—2004年195个GPS站的观测值研究了PWV的变化特征与ENSO之间的相关性。文献[187]利用多种气象模型和两种再分析数据(ERA-40和MERRA)计算了ENSO活动期间海表温度对大气水汽反馈强度的影响。文献[188]利用中国南海和东南亚的PWV数据分析了水汽变化与海表温度、厄尔尼诺指数之间的相关性。文献[189]采用2007—2015年GPS站点坐标、重力恢复与气候实验卫星(gravity recovery and climate experiment，GRACE)的等效水深和PWV数据研究了云南干旱事件期间三者的协同关系。文献[190]利用经验正交分解分析了ENSO发生期间中国南部区域PWV和温度的时空演变特征，并基于PWV和温度提出了一种新的ENSO表征指数。文献[65]利用改进的SSA方法对全球56个临海地区的GPS站PWV序列进行信号分解，并对提取得到的非线性趋势与月异常均值进行了相关性分析，基于此研究了PWV与ENSO现象的关系。文献[191]研究了ENSO对俄罗斯PWV年际变化的影响。

图4 地基GNSS大气水汽与海表温度的相关系数[65]

3.3 PWV时序气候应用未来研究及应用展望

关于在非一致性探测方法，相对探测方法是使用最多的探测方法。对PWV而言，常用的参考数据源为并址GNSS站、探空站和再分析数据。但是这几种参考数据具有一定局限性：①受限于地面站点的空间分布和再分析数据的精度，参考数据往往较难获取；②参考数据可能存在非一致现象，故差分序列的探测结果无法确定探测结果的来源；③差分算子可能会掩盖非一致现象，特别是阶跃量较小时，容易受到参考数据和目标数据本身噪声的影响。绝对方法不需要外部参考数据，可避开上述问题，从理论上讲，绝对探测方法是PWV时间序列非一致性探测的最优选择，但受限于当前对PWV时间演变特征认知的局限性，目前的信号分解方法无法给出精确的信号分解结果，故绝对探测方法的准确性和虚警率仍存在一定问题。因此，一致性探测方法需要从两个方面入手：①对于相对探测方法，最佳的参考数据为与目标数据同一系统的观测值，因此增加GNSS观测值的数量，打破不同业务部门GNSS数据的壁垒，同时提高参考数据的可靠性，有利于提高相对探测方法的精度。②对于绝对探测方法，需要提高对PWV时间演变特征和观测噪声特性的认知，进一步提高信号分解的精度，从而改进绝对探测方法的探测效果。

针对PWV的时间演变特征，当前的研究方法基本都是基于长期趋势项和周期振幅、相位为常数的假设。然而，PWV作为大气状态的表征，易受到如短期极端天气，持续数年的厄尔尼诺变化、干旱等不同气象因素的影响，PWV时间变化呈现出明显的非线性变化趋势，故将PWV长时间序列中的趋势项、周期项作为常数项的拟合结果不具有统计意义上的显著性。因此，应充分认识PWV时间序列中的非线性时间演变特征，进一步探测非线性信号分解方法，提高信号分解结果的显著性，并借此分析大气水汽变化与全球变暖、厄尔尼诺等多种气候现象的关联性，从而研究大气水汽变化与不同气象因素的耦合关系，深入认识大气水汽变化的主要驱动力。

4 结束语

GNSS大气水汽探测遥感作为一个跨学科的、新兴的研究方向，通过综合利用GNSS反演的大气资料对于提高各类气象模型的预报性能，开展针对极端气象事件预报及气候变化监测应用有着非常积极的作用。本文系统阐述了地基GNSS大气水汽探测遥感技术及其应用的发展过程，首先介绍了近年来包括对流层延迟、大气可降水量等多类型GNSS大气参数高精度反演的研究进展；然后对GNSS大气反演在极端天气短临预报及气候变化现象解释两个方向的研究工作进行了科学探析；最后，阐明了GNSS大气水汽探测遥感技术面临的主要挑战及未来研究展望。当前，虽然地基GNSS大气反演技术已经有了较大的发展，然而，作为非常规的气象观测资料，各类GNSS大气产品从实时反演到质量控制与误差分析，再到在气象和气候领域的有效利用，都仍存在诸多问题需进一步探究解决。其中，伴随北斗导航卫星系统的全面建成使用，如何进一步完善基于多系统GNSS观测高精度反演大气产品的处理策略，探究不同映射函数、信号频率、对流层约束类型、多系统定权以及大气负载等因素对反演过程及产品精度的影响，寻求一套最优的多系统GNSS反演大气信息处理方法，并利用各类模型方法有效实现其在极端天气短临预警预报与气候变化分析中的应用，将是未来GNSS大气反演的研究重要方向。此外，近年来，全球范围GNSS连续观测站网中大气探测数据的不断累积，为大地测量领域时空大数据分析带来前所未有的新机遇，如何充分挖掘并高效利用GNSS历史大数据中蕴含的信息，也有待进一步开展更为深入系统的研究。