中国大陆地区GNSS ZTD对超强东部型厄尔尼诺事件的响应

于腾丽 王 勇,2 黄 靖 杨 军 占 伟

1 天津城建大学地质与测绘学院，天津市津静路26号，300384 2 中国科学院精密测量科学与技术创新研究院大地测量与地球动力学国家重点实验室，武汉市徐东大街340号，430077 3 中国地震局第一监测中心，天津市耐火路7号，300180

水汽(precipitable water vapor, PWV)作为地球大气的重要组成成分，是引发暴雨洪涝灾害的要素之一[1]。地基GNSS反演水汽具有全覆盖、高精度、全天候等特点[2],被广泛应用于降雨、中尺度恶劣天气和全球气候变化预测[3-4]。来自西太平洋地区的水汽输送是厄尔尼诺事件影响中国气象因素的重要纽带[5]。研究表明，中国大陆GNSS水汽变化对厄尔尼诺事件的响应存在区域性规律，与气候类型、地理区划有关，且GNSS水汽可作为监测指标，指示厄尔尼诺事件的发生及演变[6-7]。2014-10～2016-04发生的超强东部型厄尔尼诺事件导致我国长江以南大部分地区降水偏多，部分地区甚至出现“冬汛”，反映出超强厄尔尼诺事件对中国大陆气候的重要影响[8-9]。本文将基于CMONOC数据分析GNSS ZTD代替水汽的可行性，并利用FFT和WT分析中国大陆GNSS ZTD空间、时域、频域变化对超强东部型厄尔尼诺事件的响应规律。

1 研究数据与方法

1.1 研究数据

1.1.1 GNSS 水汽

CMONOC连续观测站分布如图1所示。GNSS 水汽数据由GNSS ZTD结合气压、温度数据反演而来[1]。由于获取的CMONOC气象资料有限，为保证GNSS 水汽时间序列的连续性，以及可分析站点的覆盖率，需避免由于同期CMONOC气象数据缺失导致GNSS水汽数据不完整的问题。因此，本文将首先分析利用GNSS ZTD代替GNSS 水汽的可行性。由于设站时间不同导致各站点积累的观测资料的长度不一，本文分析站点的最长观测时间序列为2008-01-01～2021-07-24，最短观测时间序列为2012-01-01～2021-07-24。

图1 CMONOC站点分布Fig.1 Distribution of CMONOC stations

1.1.2 厄尔尼诺事件

1.2 研究方法

FFT可从频域上提取信号在时域上无法体现的周期特征。本文将利用加汉宁窗的FFT分析GNSS ZTD和IEP的频率/周期的变化规律，从时域和频域上探究二者的关系。

FFT可提取信号的显著变化周期，而WT可提取显著周期的对应层数，获取周期的具体变化情况。本文选择紧支撑标准正交小波db6小波基函数对GNSS ZTD进行分析。小波变换分层原理是采用2，4，8，16，…，2J作为尺度，将信号f(t)(频段0～F)分解为J+1个信号频段，其中0～F/2J频段为低频项(AJ)，F/2n～F/2n-1(n=1，2，3，…，J)为高频项(D1,D2, …,DJ)。分层合格的标准为低频信号AJ曲线呈单一递增或递减趋势。经试算，将GNSS ZTD时间序列划分为17层可满足分层合格标准，其中D1～D17层为高频项，A17层为低频项。

2 GNSS ZTD代替GNSS 水汽可行性分析

本文获取的CMONOC气象资料时间跨度为2013-01～2020-12，结合该时段的GNSS ZTD数据反演出各站点的GNSS水汽序列，并开展GNSS ZTD与GNSS水汽序列的相关性分析。利用FFT对GNSS ZTD与GNSS 水汽时间序列进行频域分析，并探究二者的频域相关性。图2为BJFS、WUHN站分析结果，图中**表示相关性在0.01层上显著。

图2 GNSS ZTD与GNSS水汽时频相关性Fig.2 Time-frequency correlation of GNSS ZTD and GNSS PWV

由图2(a)、2(c)可知，GNSS ZTD与GNSS水汽原始序列吻合度较高，相关系数在0.9以上。从图2(b)、2(d)可以看出，GNSS ZTD与GNSS 水汽频域分析结果较为吻合，均存在极显著的12个月变化周期(对应频率约为0.085)以及较显著的9～10个月、5～6个月周期，且在频率区间(0.3, 1.25)内也存在较显著的变化周期。综上可知，GNSS ZTD与GNSS水汽序列在时域和频域上均存在极强相关性，吻合度较高，利用GNSS ZTD代替GNSS 水汽开展研究具有可行性。

3 GNSS ZTD时空变化对超强东部型厄尔尼诺事件的响应分析

本节对此次超强厄尔尼诺事件期间GNSS ZTD与IEP的相关性进行分析。GNSS ZTD原始时间序列存在显著的年周期、季节性变化周期以及短期的随机波动与噪声。而利用FFT分析IEP的显著变化周期发现，其存在显著的80个月(约7 a)、40个月(约3.5 a)、18个月(1.5 a)、14个月(约1 a)周期，这与GNSS ZTD的显著变化周期存在较大差异。故对二者进行相关性分析前，需要先排除GNSS ZTD随机波动与噪声的干扰。利用WT方法对GNSS ZTD时间序列进行分层处理。IEP显著变化周期分别对应小波分层的D15(4～8 a)、D14(2～4 a)、D13(1～2 a)层，故将GNSS ZTD进行WT后的D13、D14、D15层进行重构。分析GNSS ZTD重构高频项与IEP的相关性，并将其与二者原始序列之间的相关性进行对比，结果如图3所示，图中sig值代表Pearson相关性双尾检验的显著性系数，sig<0.01表示相关性在0.01层上显著，sig<0.05表示相关性在0.05层上显著。

图3 GNSS ZTD与IEP相关性分析Fig.3 Correlation analysis between GNSS ZTD and IEP

由图3(a)可知，中国大陆各区域GNSS ZTD与IEP存在正相关性，且均可通过0.01层显著性检验，基本呈弱或中度正相关。图3(b)显示，GNSS ZTD重构高频项与IEP的相关性相比于原始序列有显著提高，其中1个站点可通过0.05层显著性检验，其余站点均通过0.01层显著性检验。大部分区域相关系数在0.4～0.6范围内，呈中度正相关；热带、亚热带、川滇地区以及黑龙江部分站点的相关系数在0.6～0.8范围内，呈强正相关；个别站点相关系数大于0.8，呈现出极强正相关。因此，超强东部型厄尔尼诺事件会导致中国大陆GNSS ZTD(水汽)增多，且对热带季风和亚热带季风区域以及川滇地区影响较为显著。可能是因为超强东部型厄尔尼诺事件导致西太平洋副热带高压偏强偏西，促使大量低纬度的水汽向我国输送，使得中国大陆区域GNSS ZTD(水汽)整体增多；又因受西北太平洋对流层低层菲律宾反气旋性异常环流的影响，中国东部出现偏南风异常，从而导致南方地区水汽含量剧烈增加[9]。同时，中国大陆热带和亚热带季风气候区域以及川滇地区所处纬度较低，GNSS ZTD(水汽)比其他地区含量高[10]，变化较为活跃。

4 GNSS ZTD显著变化的周期对超强东部型厄尔尼诺事件的响应分析

GNSS ZTD存在显著的年、9～10个月、半年周期以及(0.3, 1.25)频率区间内的较显著周期。本节将开展超强东部型厄尔尼诺事件对GNSS ZTD显著变化周期的影响分析。为确保结果的准确性和可靠性，选择多个正常气候时段的数据与之进行对比。此次厄尔尼诺事件持续时间长达19个月，而本文涉及时间序列中只包含一段与之对应的正常气候时段。因此，以该事件峰值时间为基准向前后共截取12个月的时间长度(2015-05～2016-04)开展分析，并选取2012-05～2013-04、2013-05～2014-04、2016-05～2017-04三个正常气候时段与之进行对比。利用FFT分析4个时间段内GNSS ZTD的显著变化周期，探究其在超强东部型厄尔尼诺事件期间与正常气候时段的差异及规律。由于截选时长为12个月，故只能对除年周期外的其他显著周期进行分析。本文截取各站点GNSS ZTD FFT结果中存在显著周期的时段进行对比，由于篇幅限制，仅展示温带季风气候地区2个站点的分析结果(图4)。

图4 温带季风气候地区GNSS ZTD显著变化周期对比Fig.4 Comparison of significant change periods of GNSS ZTD in temperate monsoon climate region

由图4可知，超强东部型厄尔尼诺事件发生期间，温带季风气候区域GNSS ZTD约9个月的显著变化周期振幅明显低于3个正常气候时段；在(0.2, 1.4)频率区间内，超强东部型厄尔尼诺事件期间的GNSS ZTD显著变化周期振幅较正常时段有不同程度的增长。

综合各区域站点的分析结果发现，在超强东部型厄尔尼诺事件的影响下，中国大陆区域GNSS ZTD约9个月的显著变化周期振幅均存在不同程度的降低，而0.8个月、1个月、1.5个月、2～3个月的较小显著变化周期振幅均有不同程度的增长。表明超强东部型厄尔尼诺期间，GNSS ZTD较长的显著变化周期减弱，较短的显著变化周期增强，GNSS ZTD变化更为活跃。

以3个正常气候时段GNSS ZTD 9个月显著变化周期对应振幅(amplitude,A)的平均值为基准，计算超强东部型厄尔尼诺期间GNSS ZTD 9个月显著变化周期振幅的变化量(ΔA)；并以相同方法分别计算各站点0.8～3个月内的显著变化周期对应振幅的增量，然后求其平均值，结果如图5所示。

图5 GNSS ZTD显著变化周期对超强东部型厄尔尼诺事件的响应分析Fig.5 Response analysis of GNSS ZTD significant periods to super eastern Pacific El Nio event

由图5(a)可知，在超强东部型厄尔尼诺事件的影响下，中国大陆绝大部分站点GNSS ZTD 9个月显著变化周期振幅均有不同程度的降低，热带地区和温带大陆西部地区振幅变化量最为显著，个别显示无变化规律的站点主要集中在中国大陆西北部地区。从图5(b)可以看出，中国大陆绝大部分站点GNSS ZTD 0.8～3个月显著变化周期振幅均有不同程度的增加，热带地区和亚热带地区振幅增量最为显著，个别显示无变化规律的站点主要集中在温带大陆地区北部。由此可知，在中国大陆西北部高原山地气候区域的个别站点9个月显著周期和中国大陆北部温带大陆气候区域的个别站点0.8～3个月显著周期对超强东部型厄尔尼诺事件响应较差。可能是因为这些站点位于中国大陆纬度较高的地区，GNSS ZTD(水汽)含量较低，变化不活跃[10]，而且在此次超强东部型厄尔尼诺事件期间，北方地区偏北风气流受到抑制[9]，影响部分区域水汽的正常流通。

5 结 语

本文首先基于CMONOC提供的GNSS连续观测数据分析GNSS ZTD代替GNSS水汽的可行性，然后开展中国大陆GNSS ZTD对2014-10～2016-04超强东部型厄尔尼诺事件的响应分析，得到以下结论：

1)中国大陆各区域GNSS ZTD与IEP呈正相关，即超强东部型厄尔尼诺事件会导致中国大陆GNSS ZTD增大，且对热带季风和亚热带季风区域以及川滇地区影响较为显著。

2)在超强东部型厄尔尼诺事件的影响下，中国大陆区域GNSS ZTD约9个月的显著变化周期振幅均存在不同程度的降低；而0.8～3个月内的较小显著变化周期振幅均有不同程度的增长。其中，热带地区与温带大陆西部地区GNSS ZTD 9个月显著变化周期、热带与亚热带地区GNSS ZTD 0.8～3个月显著变化周期对此次厄尔尼诺事件的响应最为显著，表明超强东部型厄尔尼诺事件期间GNSS ZTD变化更为活跃。