和 會 聞洪峰 付利釗

（1.河北省基礎(chǔ)地理信息中心，河北 石家莊 050032；2.河北省第一測繪院，河北 石家莊 050032）

0.引言

GNSS信號在穿過霧霾區(qū)域時，霧霾會對GNSS信號的傳播產(chǎn)生什么樣的影響，對GNSS相對定位基線解算是否會產(chǎn)生影響，對GNSS精密單點定位是否會產(chǎn)生影響，對GNSS反演對流層延遲是否會產(chǎn)生影響，這種影響有多大，這種影響來源哪里？這些方面的研究相對還比較少，也不系統(tǒng)全面。2014年王勇等發(fā)現(xiàn)天頂對流層延遲與可降水量在霧霾天會有升高[1]；2016年田宗彪等研究了霧霾對精密單點定位的影響[2]；2017 年 魏 二 虎 等 以BJFS、LHAZ、CHAN、JFNG、SHAO等5個IGS站秋冬季的觀測數(shù)據(jù)為例，分析得出可吸入顆粒物濃度變化與精密單點定位的精度沒有顯著相關(guān)性[3]。筆者以HBCORS基準(zhǔn)站2014年9月至2015年2月的連續(xù)觀測數(shù)據(jù)為例，結(jié)合CORS站所在地區(qū)AQI、PM2.5、PM10濃度日均值，詳細(xì)分析了大氣顆粒物濃度對GNSS基線解算、精密單點定位和GNSS反演對流層濕延遲的影響，探索霧霾天氣對GNSS對流層濕延遲產(chǎn)生影響的原因。

1.大氣顆粒物濃度對中長基線的影響

GNSS相對定位的基線解算一般是利用兩臺以上接收機的載波觀測值和偽距觀測值，組成雙差觀測值為基礎(chǔ)量，配以組合觀測值為輔助，解算接收機的相對位置關(guān)系。在基線解算過程中，基線長度對于解算結(jié)果會產(chǎn)生影響。對于相距15km以內(nèi)，高差不大的基線，雙差就可以基本消除如對流層延遲誤差在內(nèi)的多項誤差源。而霾天氣就發(fā)生在對流層，霾產(chǎn)生的影響在基線過短時經(jīng)過雙差處理后會被削弱甚至消除，因此本文選擇50km以上的中長基線作為研究對象。

以AQI日均值為霾污染程度判斷標(biāo)準(zhǔn)，假設(shè)基線端點分別為A和B，按照以下4種情況分類統(tǒng)計：

（1）基線一端站點A有重度以上霾發(fā)生（AQI＞200），另一端站點B空氣質(zhì)量良好（AQI＜100）；

（2）基線一端站點B有重度以上霾發(fā)生（AQI＞200），另一端站點A空氣質(zhì)量良好（AQI＜100）；

（3）基線兩端站點A、B均有重度以上霾發(fā)生；

（4）基線兩端站點A、B空氣質(zhì)量均良好。

表1 基線LHAZ-BDZZ各污染程度基線解算結(jié)果 單位：m

LHAZ—BDZZ的基線北方向（N）、東方向（E）、橢球高（U）和基線長度（DIS）的統(tǒng)計情況（如表1所示）?？梢钥闯鰧τ谶@種超過1000km的長基線，基線各個分量以及基線長度均沒有明顯因為霾天氣的影響，在N方向、E方向、U方向都沒有產(chǎn)生系統(tǒng)性的偏差。也就是說無論霾天氣如何變化，基線分量和基線長度的重復(fù)性都是相當(dāng)?shù)摹?/p>

表2 基線BJFS-BDAG各類污染程度基線解算結(jié)果 單位：m

BJFS-BDAG基線各分量及基線長度的統(tǒng)計值（如表2所示），顯然在各種霧霾天氣下，基線精度相當(dāng)，基本沒有變化，完全在誤差范圍內(nèi)波動。無論是2000km以上的長基線還是100km以上的中等基線，霧霾天氣都不會對基線解算精度造成顯著影響，且沒有體現(xiàn)出系統(tǒng)性的偏差。

因此，在相對定位解算基線的過程中，顆粒物濃度對基線解算的精度沒有明顯影響，可以不予考慮。

2.顆粒物濃度對精密單點定位的影響

本文以CORS站的實測數(shù)據(jù)為例，利用GPAS軟件求解各站點的非差精密單點定位日解，并與當(dāng)?shù)氐腁QI日均值進行比較分析，以驗證霾天氣是否會對精密單點定位的精度產(chǎn)生影響。對于沒有準(zhǔn)確的三維坐標(biāo)的CORS站點，將非差精密單點定位日解的均值作為基準(zhǔn)值。然后將精密單點定位日解與基準(zhǔn)值求差，以日解與基準(zhǔn)值的偏差來與AQI指數(shù)進行比較。

圖1 精密單點定位U方向坐標(biāo)值殘差

BDFP和TSZH兩個站2014年9月至2015年3月期間AQI和精密單點定位U方向殘差的比對（如圖1所示），從圖中看不出明顯的相同的趨勢性變化，AQI指數(shù)與精密單點定位的精度都沒有顯著的相關(guān)性。這期間各站坐標(biāo)不同分量偏差與AQI的相關(guān)系數(shù)（如表3所示），也可以明顯看出各個分量與AQI的相關(guān)性都很小，都未通過顯著性檢驗，不具有統(tǒng)計特性。

表3 精密單點定位U方向坐標(biāo)值殘差與大氣顆粒物濃度相關(guān)性

3.大氣顆粒物濃度對GNSS 對流層延遲的影響

從前面分析可以看出，大氣顆粒物濃度的變化并沒有對GNSS基線解算精度和精密單點定位精度產(chǎn)生顯著的影響。大氣顆粒物大量懸浮于對流層底層，隨對流層大氣運動和氣象條件變化而聚集和消散，因此，其濃度變化應(yīng)該與GNSS對流層延遲有一定的關(guān)系。本文分析大氣可吸入顆粒物濃度對GNSS對流層延遲的影響。

GNSS對流層延遲分為干延遲和濕延遲兩部分，干延遲通過對流層延遲模型可以得到高精度改正。對流層濕延遲則需要經(jīng)過數(shù)學(xué)誤差模型在GNSS數(shù)據(jù)處理過程中估計。本文以GNSS對流層濕延遲和大氣顆粒物濃度為研究對象，分析霾天氣持續(xù)過程中，二者變化趨勢。

BDFP和TSQA兩個CORS站點兩次持續(xù)重霧霾的發(fā)生、持續(xù)直至消散的過程期間，天頂對流層濕延遲（ZWD）和PM2.5的變化趨勢（如圖2所示）。以BDFP站點為例，在2014年10月5日至10月13日在保定市阜平縣發(fā)生了一次較為嚴(yán)重的持續(xù)霧霾天氣，從10月6日夜間（25歷元）開始，PM2.5濃度從不到100μg/m3開始逐步攀升，到10月8日達到接近500μg/m3最高點，并在高位持續(xù)徘徊，直到11日夜間（155歷元）開始迅速降低，霧霾消散。在此期間，ZWD也是從10月6日夜間開始隨著PM2.5濃度同步攀升，從不到100mm攀升到160mm左右。之后在整個霧霾持續(xù)過程中，ZWD一直在160mm上下高位震蕩，在10月11日（155歷元）開始下降，一直降至不到40mm。唐山市遷安市2014年10月27日至11月2日期間發(fā)生的一次持續(xù)霧霾天氣過程（如圖2（b）所示），其ZWD和PM2.5濃度的同步變化過程與圖2（a）是高度相似的。

由以上分析可以看出，整個霾天氣持續(xù)過程中，ZWD和PM2.5濃度有著相似的升降過程，但升降的時間點略有偏差，升降速度也略有不同，PM2.5濃度變化曲線比較陡，而ZWD的變化曲線要相對平緩。

圖2 大氣顆粒物濃度與大氣可降水量同步變化趨勢圖

對流層濕延遲ZWD與大氣可降水量，也就是整層大氣的水汽含量是呈一定正比例關(guān)系的，只要知道一個區(qū)域的加權(quán)平均溫度Tm，就可以建立ZWD與大氣可降水量之間的數(shù)學(xué)關(guān)系[5]。在GNSS氣象學(xué)中，GNSS對流層濕延遲通常都被用于反演大氣可降水量。然而在霧霾天氣發(fā)生時，ZWD與顆粒物濃度有著相近的升降趨勢，那顆粒物是否對GNSS反演對流層濕延遲有影響呢？如果有影響，那么在GNSS對流層濕延遲反演大氣可降水量時就需要減去大氣顆粒物的影響。

圖3 對流層濕延遲ZWD與整層大氣可降水量的線性回歸方程

對ZWD與大氣可降水量呈現(xiàn)出很強的線性關(guān)系（如圖3所示），斜率在0.15至0.16之間，與眾多學(xué)者們此前的研究成果相符，也印證了二者計算結(jié)果的可靠性。本地加權(quán)平均氣溫一般需要長時間的探空資料才能計算，而且并非累積時間越長效果越好，最近一年的數(shù)據(jù)分析成果優(yōu)于之前幾年累積數(shù)據(jù)的分析成果。本文沒有研究區(qū)域長時間累積的探空數(shù)據(jù)來統(tǒng)計估算加權(quán)平均溫度。就直接取貝葉斯經(jīng)驗公式計算加權(quán)平均溫度。利用歐洲天氣預(yù)報中心氣象再處理數(shù)據(jù)ERA5的整層大氣水汽含量（IWV）通過轉(zhuǎn)換系數(shù)即可反演水汽造成的對流層濕延遲，用GNSS對流層濕延遲減去通過整層大氣水汽含量估算的對流層濕延遲，就可以得到對流層濕延遲的剩余殘差量，進而分析去掉大氣水汽影響后對流層濕延遲的殘差與大氣顆粒物濃度是否存在相關(guān)性。

9個監(jiān)測站ZWD殘余誤差的平均值、中位值、標(biāo)準(zhǔn)差以及其與AQI、PM2.5、PM10的相關(guān)系數(shù)（如表4所示），可見不同監(jiān)測站的ZWD殘余誤差平均值是有差別的，這表明ERA5提取的大氣水汽含量反演的ZWD與GNSS實測反演的ZWD之間有系統(tǒng)性誤差，且在不同區(qū)域這一系統(tǒng)性誤差存在差異。各個監(jiān)測站的ZWD殘差標(biāo)準(zhǔn)差都在8mm左右，折合成大氣可降水量約在1-2mm。除了TSQA站，其余各站ZWD與PM2.5、PM10、AQI的相關(guān)系數(shù)都低于0.2甚至0.1，且有正有負(fù)，沒有顯著的相關(guān)性。TSQA一個站的微弱正相關(guān)性不具有代表性，綜合來看，顆粒物濃度的變化與ZWD殘余誤差之間沒有顯著的相關(guān)性。大氣顆粒物濃度的上升下降并沒有直接影響到ZWD的升降變化。

表4 ZWD剩余殘差與大氣顆粒物濃度相關(guān)性

GNSS信號在穿過霾區(qū)域時，如果大氣顆粒物會對GNSS信號產(chǎn)生延遲影響，則ZWD去掉水汽影響后的殘差應(yīng)該與大氣顆粒物濃度有一定的相關(guān)性，而且應(yīng)該是正相關(guān)性。而從上述分析中可以看出大氣顆粒物濃度與ZWD剩余殘差并沒有體現(xiàn)出顯著的相關(guān)性。由此可以推斷：在目前的測量精度上，大氣顆粒物對GNSS信號傳播不會產(chǎn)生顯著的直接影響，通過GNSS對流層濕延遲反演大氣可降水量并不必考慮大氣顆粒物濃度的影響；對流層濕延遲與大氣顆粒物濃度在霾天氣有相似變化的主要原因在于大氣可降水量與大氣顆粒物濃度在霾天氣時有著相似變化趨勢。

4.結(jié)論

本文分析了大氣顆粒物濃度對GNSS數(shù)據(jù)處理的影響以及大氣顆粒物濃度與PWV相關(guān)性起因。得出以下結(jié)論：

（1）霧霾天氣，大氣顆粒物的高濃度對于GNSS基線解算不會產(chǎn)生顯著影響，基線解算精度與大氣顆粒物濃度之間沒有顯著相關(guān)性。

（2）霧霾天氣，大氣顆粒物的高濃度對于GNSS精密單點定位日解精度沒有顯著影響，精密單點定位與大氣顆粒物濃度之間也沒有顯著的相關(guān)性。

（3）霧霾天氣，大氣顆粒物濃度與對流層濕延遲會有相似變化趨勢，但二者不是直接相關(guān)關(guān)系，本質(zhì)上是大氣顆粒物濃度與PWV存在相關(guān)性。積累階段，大氣顆粒物濃度比ZWD或者說PWV上升要更快。