北極地區(qū)海霧數(shù)值模擬及其電波衰減研究

張玉生 郭相明 郝曉靜 趙強 張蕊

(電波環(huán)境特性及?；夹g(shù)重點實驗室, 中國電波傳播研究所, 山東 青島 266107)

0 引言

海霧是海上一種常見的天氣現(xiàn)象, 由于海霧降低了海上能見度, 會使航行的船只迷失方向或造成擱淺、碰撞等重大事故, 因此即使是在現(xiàn)代先進(jìn)航海技術(shù)條件下, 在海霧中航行仍要百倍警惕[1]。海霧內(nèi)濕度大, 在海霧頂以上區(qū)域濕度隨高度迅速減小, 并且在海霧頂部易形成逆溫結(jié)構(gòu),易形成大氣波導(dǎo)結(jié)構(gòu), 易引發(fā)超短波、微波頻段大氣波導(dǎo)超視距傳播或者盲區(qū)現(xiàn)象[2-4]; 海霧內(nèi)液滴對于30 GHz 以上的毫米波頻段以及紅外等波段的傳播會產(chǎn)生嚴(yán)重的衰減效應(yīng)。

國內(nèi)外對海霧進(jìn)行了不少觀測、機理、特征和預(yù)報研究, 例如彭邦安[5]對南極麥克韋爾灣的海霧分析, 楊清華等[6]對南極長城站海霧的特征分析等。對海霧的預(yù)報方法主要有天氣學(xué)方法、統(tǒng)計方法和數(shù)值預(yù)報方法, 例如許淙等[7]利用南極長城站1985—2006 年的氣象水文觀測資料和美國國家環(huán)境預(yù)報中心(NCEP)再分析數(shù)據(jù), 對南極長城站附近海霧的氣象特征進(jìn)行了統(tǒng)計分析,基于數(shù)據(jù)生成了500 hPa 氣象要素場和海平面氣象要素場, 初步建立了長城站夏季海霧預(yù)報系統(tǒng),預(yù)報有霧或無霧結(jié)果; 林曉能[8]基于馬爾柯夫鏈原理建立了南海海霧的短期預(yù)報方法, 并得到了很好的應(yīng)用; 楊棋等[9]將洋山海域不同范圍的氣溫、相對濕度、氣溫與海溫之差、風(fēng)向、風(fēng)速等因子進(jìn)行組合, 建立了統(tǒng)計預(yù)報方程, 實現(xiàn)了海霧客觀預(yù)報方法。數(shù)值預(yù)報方法也先后經(jīng)歷了一維數(shù)值、二維數(shù)值和三維數(shù)值等階段, 隨著各種尺度大氣數(shù)值模式的發(fā)展, 海霧的數(shù)值模擬精度也隨之提高[10-11]。

北美洲、歐洲北部和亞洲北部國家之間交通往來直線距離最短的路徑需穿越北極地區(qū), 因此北極地區(qū)戰(zhàn)略位置重要。北極地區(qū)是全球氣候變化的敏感地區(qū), 氣候惡劣, 基本上處于多冰雪、多霧、極寒低溫狀態(tài), 暴風(fēng)雪多發(fā), 在北極圈內(nèi)有極晝極夜現(xiàn)象。北極上空臭氧含量最多, 對于擴散到對流層的臭氧不僅能吸收地-氣系統(tǒng)的長波輻射而加熱大氣, 還可以參與大氣光化學(xué)反應(yīng), 進(jìn)而改變其他溫室氣體的含量和分布, 影響地-氣系統(tǒng)的輻射平衡。北極處于高緯度地區(qū), 在研究大氣中尺度運動時, 考慮的柯氏力等受力也不一樣。海冰的存在也阻礙了海洋與大氣之間的直接動量和熱量交換, 使其成為氣候系統(tǒng)中除大氣、海洋、陸地之外的又一重要組成部分, 海冰、融池、冰間水道、淺薄積雪等特征增加了北極地表復(fù)雜性, 也增加了冰雪下墊面物理參數(shù)表述的復(fù)雜性, 并且由于極區(qū)的下墊面覆蓋大量冰雪, 太陽能量的反照率和吸收率明顯與中低緯度不同。近年來人們持續(xù)研究北極地區(qū)的海、氣、冰等的能量交換關(guān)系, 陳立奇等[12]利用聯(lián)合冰站觀測數(shù)據(jù), 模擬了北冰洋夏季大氣邊界層結(jié)構(gòu)和下墊面能量平衡的變化特征, 定量得出了北冰洋夏季海/氣和冰/氣之間湍流通量和邊界層參數(shù)的差異。

由于北極特殊地理環(huán)境和惡劣氣候多變, 目前北極地區(qū)氣象觀測和研究投入不足, 只有有限的氣象觀測站, 無法有效獲取更大區(qū)域的北極低層大氣精細(xì)結(jié)構(gòu)和海霧液滴數(shù)據(jù)。由于缺乏足夠的數(shù)據(jù)和機理分析, 很難建立準(zhǔn)確描述北極地區(qū)實際狀況的氣象方案, 因此北極海霧的生消機理、數(shù)值模擬和預(yù)報還有待進(jìn)一步研究。在極區(qū)的海霧研究主要集中在特征分析以及基于統(tǒng)計意義上的客觀預(yù)報, 預(yù)報主要是定性預(yù)報, 三維數(shù)值模擬主要集中在南海、黃海海域的海霧模擬和預(yù)報[13], 北極地區(qū)的海霧數(shù)值模擬(預(yù)報)研究很少。

本文根據(jù)全球再分析數(shù)據(jù), 基于適合于極區(qū)的中尺度大氣數(shù)值模式, 模擬了2010 年7 月27日北冰洋的濃霧發(fā)展過程, 并基于模擬的海霧液態(tài)含水量和大氣參數(shù)進(jìn)行了相應(yīng)的電波傳播評估,進(jìn)一步評估了海霧對毫米波和遠(yuǎn)紅外波段上的單位距離衰減值。通過北極海霧的數(shù)值模擬及其電波傳播特性研究, 研究其發(fā)生區(qū)域、強度和時空變化, 對航海氣象安全以及遠(yuǎn)距離導(dǎo)航、通訊和探測、編隊間的通訊和探測保障都有重大意義。

1 極區(qū)中尺度大氣數(shù)值模式介紹

數(shù)值預(yù)報模式是目前研究天氣過程和模擬電波環(huán)境的重要方法, 海霧屬于中小尺度對象, 因此預(yù)報北極海霧采用中尺度數(shù)值模式比較合適。由于中尺度數(shù)值模式既考慮了大尺度天氣系統(tǒng)對所預(yù)報區(qū)域大氣的整體影響, 也考慮了大氣邊界層、積云等小微系統(tǒng)對所預(yù)報區(qū)域大氣特別是低層大氣的局部影響, 因此, 對于海霧的伴隨逆溫、液水密度等氣象環(huán)境參數(shù)會有很好的數(shù)值模擬和預(yù)報[14-16]。

數(shù)值預(yù)報的歷史很短, 真正意義上的數(shù)值預(yù)報開始于20 世紀(jì)初期, 隨著計算機技術(shù)的不斷發(fā)展, 后期數(shù)值預(yù)報的時效性和準(zhǔn)確率也不斷提高,先后成功研制了 Mesoscale Model(MM)系列、Advanced Regional Prediction System(ARPS)、Regional Atmospheric Modeling System(RAMS)、Weather Research and Forecasting Model(WRF)等中尺度數(shù)值預(yù)報模式。隨著資料同化效果的改善以及計算性能的提高, 中尺度數(shù)值模擬結(jié)果越來越準(zhǔn)確, 中尺度數(shù)值模擬逐漸成為氣象環(huán)境、電波環(huán)境模擬和預(yù)報技術(shù)發(fā)展的重要方向[17-19], 它的優(yōu)勢在于可以獲取相對較高精度的隨時間和空間變化的網(wǎng)格數(shù)據(jù)。

WRF 模式是基于業(yè)務(wù)預(yù)報和天氣研究需要所發(fā)展的新一代中尺度數(shù)值天氣預(yù)報模式, 該模式是在MM5 基礎(chǔ)上發(fā)展起來的。WRF 適用于幾米到數(shù)千千米尺度大氣現(xiàn)象的研究, 可應(yīng)用在參數(shù)化模型研究、數(shù)據(jù)同化研究、業(yè)務(wù)預(yù)報研究、耦合模式應(yīng)用和教學(xué)等領(lǐng)域, 因此也適合于海霧的數(shù)值模擬、機理分析和預(yù)報[17]。但由于普通型WRF 模式中的各種氣象方案極少考慮到極區(qū)這種氣候、地理環(huán)境極端的地區(qū), 一般只應(yīng)用于中低緯度的數(shù)值模擬。目前針對極地區(qū)域大氣環(huán)流的數(shù)值模擬已取得很大進(jìn)展, 如先后研制了Polar MM5、Polar WRF(PWRF)等適合極區(qū)大氣數(shù)值模擬的模式, 這些模式具有多個可供選擇的物理過程參數(shù)化方案, 已經(jīng)成功應(yīng)用于南極和格陵蘭島大氣研究, 并為美國南極科學(xué)考察提供氣象服務(wù)。

在極區(qū)大氣數(shù)值模式中對參數(shù)化方案的改進(jìn)主要表現(xiàn)為以下幾個方面: 改善冰相微物理過程的顯式方案, 修改云輻射相互作用, 優(yōu)化湍流參數(shù)化, 增加海冰下墊面類型并改變通過雪蓋/冰面熱量輸送過程的處理方法。在利用普通型中尺度大氣數(shù)值模式過程中, 研究發(fā)現(xiàn)由于模擬過多向下的長波輻射通量而導(dǎo)致冬季南極冰蓋區(qū)地面溫度存在偏高的問題, 對于南極[20-21]和美國大陸[22]也存在模擬云量過多的問題, 這是因為原模式中采取的云量、云水路徑分別是相對濕度、溫度的簡單函數(shù)處理方法。為了改善模式在極區(qū)的適用性, 在Reisner 顯式微物理參數(shù)化中, 采用了考慮冰核濃度的Meyers 等[23]方案代替Fletcher[24]的相應(yīng)處理方法, 并采用云水和冰混合比方法計算云量以改進(jìn)大氣傳輸?shù)拈L波和短波輻射, 通過這個改進(jìn)減小了過大云水路徑, 避免過大的云水路徑在冬季南極冰蓋區(qū)導(dǎo)致過多的向下的長波輻射通量而導(dǎo)致的地面溫度偏高。

極地海冰帶是海-氣交界面具有熱力絕緣效應(yīng)的區(qū)域, 在海-氣相互作用方面, 因為海冰的反照率反饋、熱力屏障等一系列作用和影響而區(qū)別于其他的無冰區(qū)域。因此, 與普通型大氣數(shù)值模式中使用的13 類地表分類相比, 極區(qū)大氣數(shù)值模式中增加了第14 類地表特征, 以刻畫極地的海冰地表特征。在極區(qū)大氣數(shù)值模式初始化中,根據(jù)海水格點處海水表面溫度計算該處海冰分布, 且其分布特征在模式積分過程中不再發(fā)生改變。此外, 分別對海冰、開闊水域和陸地點處的通量和地表溫度進(jìn)行處理, 然后再使其與大氣相互作用[25]。

2 北極海霧數(shù)值模擬研究

2010 年7 月27 日我國第四次北極科學(xué)考察隊搭乘“雪龍”號極地科學(xué)考察船進(jìn)行海洋綜合調(diào)查時, 在72°N、153°W 附近的北冰洋海域遇到了霧鎖北冰洋的景觀。雖然限于條件無法搜集到行駛途中“雪龍”號船載觀測數(shù)據(jù)的準(zhǔn)確能見度數(shù)據(jù), 但鑒于從當(dāng)時新華網(wǎng)隨船記者張建松7 月28 日電報道[26]: “乘坐‘雪龍’號極地科學(xué)考察船來到了《幽靈船》的故事發(fā)生地北冰洋, 遭遇到影片中那樣的彌天大霧, 周圍景象看上去甚至更為詭異?！?7 日, 考察隊乘坐‘雪龍’號來到北緯72 度18 分、西經(jīng)152 度34 分的加拿大海盆深水站位, 進(jìn)行8 個小時的海洋綜合調(diào)查……午后時分, 突然飄來一陣濃霧, 霎時間天色變得陰暗昏沉。冒著濃濃海霧, ……無邊的霧氣緊緊籠罩著‘雪龍’號, 似乎已不知身在何處, 路在何方, 那感覺甚是奇妙?！灰? ‘雪龍’號以每小時三四海里的航速, 在冰海濃霧中緩緩前行, 終于順利抵達(dá)了北緯72 度42 分、西經(jīng)153 度29 分的海洋調(diào)查站位……”的描述可看出當(dāng)時“雪龍”號極地科學(xué)考察船遇到的是一場嚴(yán)重影響航行的濃霧, 當(dāng)時能見度很小, 已嚴(yán)重影響航速。

在北極進(jìn)行海霧的數(shù)值模擬和預(yù)報必須考慮和采用符合北極地區(qū)的優(yōu)化湍流參數(shù)的邊界層、輻射等氣象方案和極地海冰等特征地表。為了研究北極海霧的特性, 數(shù)值模擬和再現(xiàn)這次北極典型海霧, 本文根據(jù)該時段全球再分析數(shù)據(jù), 采用了考慮諸多極區(qū)因素的極區(qū)數(shù)值模式PWRF, 通過改進(jìn) Noah 陸面過程方案、長波輻射方案RRTMG 方案等設(shè)置, 以及特征高度、積分步長的優(yōu)化, 數(shù)值模擬(預(yù)報)了2010 年7 月27 日“雪龍”號在北冰洋遇見的濃霧液態(tài)含水量和其他氣象參數(shù), 并模擬(預(yù)報)了其發(fā)展過程。此外, 選擇時間段為2014 年12 月29 日世界00 時至30 日世界00 時, 采用基于全球再分析數(shù)據(jù)的背景數(shù)據(jù)及相同的模式水平分辨率、各類氣象方案、特征高度、積分步長等參數(shù), 對163°W—100°W, 60°N—88°N 的包含美國阿拉斯加、波弗特海及北冰洋的廣大區(qū)域的液態(tài)含水量和其他氣象參數(shù)進(jìn)行了模擬, 以側(cè)面驗證對“雪龍”號所經(jīng)霧區(qū)的模擬情況。

在兩例大氣數(shù)值模式中采用的各類方案如下。

(1)模擬區(qū)域: 第 1 例(海霧)為 120°W—180°W, 65°N—90°N; 第 2 例(云)為 100°W—163°W, 60°N—88°N。

(2)模式水平分辨率: 30 km。

(3)垂直分層: 38 η, 分別為1.000、0.999、0.9983、0.998、0.997、0.9965、0.996、0.9955、0.995、0.9945、0.994、0.993、0.9922、0.991、0.989、0.987、0.978、0.960、0.950、0.935、0.924、0.914、0.904、0.880、0.848、0.800、0.747、0.700、0.646、0.545、0.444、0.343、0.242、0.141、0.091、0.061、0.020、0.000。

(4)微物理過程方案: WSM 3 類簡單冰方案。

(5)長波輻射方案: RRTMG 方案。

(6)短波輻射方案: Dudhia 方案。

(7)調(diào)用輻散物理方案的時間間隔: 30 min。

(8)近地面層方案: Monin-Obukhov 方案。

(9)陸面過程方案: 改進(jìn)的 Noah 陸面過程方案。

(10)邊界層方案: YSU 方案。

(11)積云參數(shù)化方案采用: 淺對流 Kain-Fritsch 方案。

(12)積云參數(shù)化方案調(diào)用時間間隔: 5 min。

RRTMG 方案包含覆蓋230 nm～56 μm 的絕大部分16 個波段, 覆蓋輻射波, 重點考慮了水汽、臭氧、溫室氣體、氣溶膠以及云水和云冰等大氣成分的作用, 基于隨機云重疊的Monte Carlo Independent Column Approximation (McICA)方法,可以反映不同高度層上和次網(wǎng)格尺度上的云物理作用。針對北極的海冰和永久性冰川/陸地冰、雪表面等特色地貌地形, 修改了Noah LSM 中海冰和永久性冰川/陸地冰表面的地表能量平衡和熱量傳輸。例如改進(jìn)海冰和冰川上的地表徑流處理,修改海冰放射率, 將雪的放射率改為0.98, 修改雪的密度, 熱容量和熱傳導(dǎo)率, 增加蒸發(fā)融化潛熱, 調(diào)整熱容量和擴散率分別適合海冰和陸地冰川, 修改土壤濕度通量等。使改進(jìn)后的Noah 能處理可變的海冰厚度和海冰積雪厚度, 并且擁有可隨季節(jié)變化的海冰反照率, 計算表面交換系數(shù)和熱通量時, 在每個格點上對開闊水域和海冰分別進(jìn)行計算, 然后再根據(jù)海冰密集度求權(quán)重平均。本文第1 個數(shù)值模擬案例從2010 年7 月27 日世界時00 時開始作為模擬(預(yù)報)起始時間, 模擬(預(yù)報)出間隔1 h 的霧區(qū)變化情況, 包含24 個時次數(shù)值模擬的153°W 處的海霧液態(tài)含水量(由液態(tài)含水量混合比表征)分布情況。液態(tài)含水量混合比是指一定體積空氣中含有的液態(tài)水與干空氣的質(zhì)量比, 在云中指云液態(tài)水混合比, 在霧中則指霧液態(tài)水混合比。圖1 為從初始時間7 月27 日00 世界時后分別經(jīng)過3、12、18、24 h 時刻153°W 處的海霧經(jīng)向剖面區(qū)域分布, 從經(jīng)向剖面區(qū)域分布可以說明海霧的高度和剖面結(jié)構(gòu)分布。已有研究常用液態(tài)含水量大于0.1 g·kg-1且高度小于600 m作為判斷海霧的判據(jù)[27-28], 并用于中緯度地區(qū)例如黃海、渤海的海霧模擬中, 本研究采用了相同的條件對極區(qū)海霧進(jìn)行判定。從每隔幾小時數(shù)值模擬圖(圖2)中可看出, 液態(tài)含水量較大區(qū)域液態(tài)水的高度都在200 m 以下并且接觸地海面, 液態(tài)含水量大于0.1 g·kg-1, 符合海霧的判據(jù)判定為霧區(qū)。事實上根據(jù)中國第四次北極科學(xué)考察資料可以證實為能見度較低的海霧, 如圖1 所示, 無論從強度還是時間變化, 從高度還是水平區(qū)域變化, 北極海霧的生消變化都很清楚, 隨著時間從弱到強演變發(fā)展, 與“雪龍”船遇到海霧從弱到強的實地景象吻合, 因此基于數(shù)值模擬技術(shù)獲取北極霧區(qū)的特性是可行的。

圖1 數(shù)值模擬初始時間2010 年7 月27 日00 世界時后不同時間的153°W 處的海霧經(jīng)向剖面區(qū)域分布Fig.1. Regional distribution of sea fog meridional profile at 153°W at different times after 00 UTC on July 27, 2010, the initial time of numerical simulation

圖2 數(shù)值模擬2010 年7 月27 日00 世界時+24 小時海拔3 m 處的海霧區(qū)域分布Fig.2. Numerical simulation of regional distribution of sea fog at an altitude of 3 m at 00:00 UTC + 24 hour on July 27, 2010

為了說明本次海霧的水平分布區(qū)域, 本文數(shù)值模擬給出27 日00 世界時后各時刻的經(jīng)緯度液態(tài)水的水平分布, 這些分布皆為海拔3 m 處的液態(tài)含水量分布, 即海拔3 m 處的海霧區(qū)域分布(如圖2 為27 日00 世界時+24 小時的液態(tài)含水量分布即海霧區(qū)域分布圖)。從圖中可以看出這次海霧水平出現(xiàn)區(qū)域很廣, 圖中A 站點為數(shù)值模擬時段航跡第1 站位(72°18'N、152°34'W), B 站點為第2站位(72°42'N、153°29'W), 圖中顯示這兩個站位均在霧區(qū)。從不同時刻的數(shù)值模擬看, 此次海霧區(qū)域、強度隨時間變化有所變化, 之所以選取海拔高度3 m 處, 是因為霧頂高度一般在幾十米到幾百米高度, 海拔高度3 m 處的液態(tài)水水平區(qū)域足以能夠代表海霧的覆蓋區(qū)域。

本文第2 個數(shù)值模擬案例如圖3 所示, 從圖中可看出, 在高空3～4 km 高度處有液態(tài)含水量較大的氣團(tuán), 離海面高。云和霧的主要區(qū)別之一是云離地/海面有一定高度, 和地/海面不相接觸,而霧則接觸地/海面。在該例數(shù)值模擬過程中, 在接觸地/海面處未見有液態(tài)含水量較大的氣團(tuán), 說明該時段在本地的天氣現(xiàn)象是云而非霧現(xiàn)象, 而在第1 個霧區(qū)模擬案例中, 液態(tài)含水量較大的氣團(tuán)接觸地/海面, 是明顯霧團(tuán), 這些都從側(cè)面說明第1 個霧區(qū)模擬案例合理可信?；谒侥芤姸扰c平流霧液態(tài)含水量經(jīng)驗公式, 利用模擬的平流霧液態(tài)含水量計算出水平能見度很小, 為0.076 km, 此外結(jié)合中國第四次北極科學(xué)考察資料, 證實了第1個模擬案例所在區(qū)域當(dāng)時發(fā)生大霧現(xiàn)象, 能見度較差, 讓“雪龍”號考察船艱難航行, 現(xiàn)場真實場景進(jìn)一步說明了第1 個霧區(qū)模擬案例合理可信。

圖3 2014 年12 月29 日00 世界時+9 小時72°N 處的液態(tài)含水量經(jīng)向剖面區(qū)域分布Fig.3. Regional distribution of water content of longitude direction at 72°N at 00:00 UTC + 9 hour on December 29, 2014

3 北極海霧電波衰減研究

海霧對微波、毫米波、紅外等波段電波傳播影響很大, 在微波頻段有可能引發(fā)大氣波導(dǎo)超視距現(xiàn)象, 在毫米波、紅外等波段引起嚴(yán)重的電波衰減現(xiàn)象。本文根據(jù)電磁輻射與介質(zhì)相互作用原理, 基于獲取的海霧液態(tài)含水量和大氣參數(shù)進(jìn)行北極海霧的毫米波、紅外波段電波傳播衰減評估,海霧液滴對傳播的電磁波能量衰減主要是通過瑞利(Rayleigh)散射、米氏(Mie)散射等作用造成的。由于霧滴粒徑一般在1～100 μm 之間, 當(dāng)霧滴的尺度與入射波的波長相當(dāng)時, 需利用米氏散射理論計算霧滴的消光截面, 這時在衰減過程中霧滴散射起主要作用, 因此主要計算霧滴的散射截面。當(dāng)霧滴的尺度遠(yuǎn)小于入射波的波長時, 可使用瑞利散射近似計算霧滴的散射特性, 要比利用米氏散射理論和其他數(shù)值算法簡單方便, 其已在計算較低頻段的降雨和濃霧傳播特性和雷達(dá)氣象上得到廣泛應(yīng)用[29-32]。

毫米波段兩代表頻率(頻段高低端100 GHz 和30 GHz), 對應(yīng)波長為3～10 mm, 由于霧滴的尺度遠(yuǎn)小于其入射波波長, 因此使用瑞利散射近似計算霧滴的散射特性。遠(yuǎn)紅外波段兩代表波長的波段高低端分別為0.1 cm 和50 μm, 該波段低端入射波波長與霧滴的尺度相當(dāng), 因此使用米氏散射理論計算霧滴的散射特性; 而該波段高端入射波波長遠(yuǎn)大于霧滴的尺度, 因此使用瑞利散射近似計算霧滴的散射特性。圖4 為毫米波段兩代表頻率(頻段高低端分別為100 GHz 和30 GHz)在不同海霧液態(tài)含水量(此處表征量為液態(tài)水密度)下的單位距離衰減值?？梢钥闯? 在30 GHz, 單位距離衰減值隨海霧含水量增大逐漸增大, 但變化平緩, 均小于2 dB·km-1;在100 GHz, 單位距離衰減值隨海霧含水量增大迅速增大, 從小于2 dB·km-1增大到大于10 dB·km-1。

圖4 毫米波段兩代表頻率在不同海霧液態(tài)含水量下的單位距離衰減值Fig.4. Unit distance attenuation values of two representative frequencies in millimeter wave band under different sea fog water content

圖5 為遠(yuǎn)紅外波段兩代表波長(波段高低端分別為0.1 cm 和50 μm)在不同海霧含水量下的單位距離衰減值?？梢钥闯? 在0.1 cm 波長上, 單位距離衰減值隨海霧含水量增大逐漸增大, 但變化平緩; 在50 μm 波長上, 單位距離衰減值隨海霧含水量增大而迅速增大, 可達(dá)700 dB·km-1; 兩波長的衰減值相差很大。從圖4、圖5 可看出, 在相同液態(tài)含水量的情況下, 遠(yuǎn)紅外波的單位距離衰減值遠(yuǎn)遠(yuǎn)大于毫米波, 從2010 年7 月27 日數(shù)值模擬的北極海霧案例可知, 海霧區(qū)域中液態(tài)水含量混合比最大為0.5 g·kg-1左右, 換算成液態(tài)水密度為0.6151 g·m-3, 毫米波單位距離衰減值根據(jù)毫米波頻率上下限(30～100 GHz)計算大概在0.8～3.5 dB·km-1左右, 遠(yuǎn)紅外波段最低端(波長50 μm)單位距離衰減值約為270 dB·km-1, 而其最高端(波長 0.1 cm)單位距離衰減值僅僅在幾dB·km-1, 兩者相差甚大。在傳播距離較長的情況下(如10 km、100 km 等), 單位距離衰減值較大的其總傳播衰減值就特別大, 所以在北極地區(qū)的海霧狀態(tài)下, 對所使用的電波無線傳輸系統(tǒng)的頻段(波段)方面有很高的要求。

圖5 遠(yuǎn)紅外波段兩代表波長在不同海霧液態(tài)含水量下的單位距離衰減值Fig.5. Attenuation value per unit distance of two representative wavelengths in far-infrared band under different water content of sea fog

綜上可知, 基于數(shù)值模擬技術(shù)可以獲取北極海霧的發(fā)生區(qū)域和變化特性。從數(shù)值模擬中可看出通過對液態(tài)水含量混合比的強度、高度變化分析及其時間、水平區(qū)域變化分析, 均能清楚地展示北極海霧的生消變化, 因此基于數(shù)值模擬技術(shù)獲取北極霧區(qū)的特性是可行的。在此基礎(chǔ)之上, 根據(jù)海霧特性可以進(jìn)一步評估北極海霧對電波的衰減效應(yīng)。

4 結(jié)論

本文針對北極地區(qū)的海霧預(yù)報可行性和特性研究, 基于考慮諸多極區(qū)因素如海冰地表類型的數(shù)值模式PWRF, 通過改進(jìn)Noah 陸面過程方案、長波輻射方案(RRTMG 方案)等設(shè)置, 以及特征高度、積分步長的優(yōu)化, 成功地數(shù)值模擬出了北極海霧及其演變發(fā)展過程, 和現(xiàn)場情況對比證明采用數(shù)值模擬獲取北極地區(qū)海霧基本信息并研究其變化機理是可行的。此外根據(jù)霧滴尺度和毫米波、遠(yuǎn)紅外波段上波段高低端入射波波長的比較,通過米氏散射和瑞利散射理論計算了不同海霧液態(tài)含水量下的霧區(qū)電波傳播衰減特性, 評估了海霧對毫米波和遠(yuǎn)紅外波段上的單位距離衰減值, 對于毫米波段兩代表頻率(頻段高低端分別為100 GHz和30 GHz), 在頻率30 GHz, 單位距離衰減值隨海霧含水量增大逐漸增大, 但變化平緩, 衰減值均小于2 dB·km-1; 在頻率100 GHz, 單位距離衰減值隨海霧含水量增大迅速增大, 衰減值從小于2 dB·km-1增大到大于10 dB·km-1。對于紅外波段兩代表波長(波段高低端分別為0.1 cm 和50 μm), 在波長0.1 cm, 衰減值隨海霧含水量增大逐漸增大, 但變化平緩; 在波長50 μm, 衰減值隨海霧含水量增大迅速增大, 衰減值可迅速增大到700 dB·km-1, 兩波長的衰減值相差很大。對于相同液態(tài)含水量霧區(qū), 遠(yuǎn)紅外波的單位距離衰減值遠(yuǎn)遠(yuǎn)大于毫米波的單位距離衰減值, 毫米波頻段高低端和紅外波段波長最高端0.1 cm 單位距離衰減值都很小, 一般在幾十dB·km-1以下, 而紅外波段波長最低端50 μm 衰減值在幾百dB·km-1,與其他波段相比其單位距離衰減值甚大, 對于較長距離的傳播, 在整個路徑上不同頻段(波段)電波傳播衰減值差別很大, 有的衰減值特別大甚至造成該頻段(波段)電波無線傳輸系統(tǒng)無法使用或中斷, 所以在北極地區(qū)的海霧狀態(tài)下, 對所使用的電波無線傳輸系統(tǒng)的頻段(波段)方面有很大的要求, 通過數(shù)值模擬(預(yù)報)北極海霧可以評估其對電波傳播衰減以及無線電信息系統(tǒng)隨時間和空間變化的影響效應(yīng)。