蘇軾鵬，王華

(1.海軍大連艦艇學(xué)院航海系 大連 116018；2.海軍大連艦艇學(xué)院軍事海洋系 大連 116018)

?

風(fēng)廓線雷達聯(lián)合RASS監(jiān)測海岸帶大氣濕度的方法研究

蘇軾鵬1，王華2

(1.海軍大連艦艇學(xué)院航海系 大連 116018；2.海軍大連艦艇學(xué)院軍事海洋系 大連 116018)

海岸帶濕度對人類生產(chǎn)生活有較大影響。文章根據(jù)風(fēng)廓線雷達和RASS (Radio Acoustic Sounding System )測量濕度的相關(guān)理論，結(jié)合實測數(shù)據(jù)計算某海濱觀測站個例的大氣濕度分布，并與探空數(shù)據(jù)進行比對，分析誤差及其原因。結(jié)果表明：風(fēng)廓線雷達聯(lián)合RASS測量濕度的方法得到的濕度結(jié)果與常規(guī)探測結(jié)果較為接近,可進一步研究實現(xiàn)其監(jiān)測海岸帶大氣濕度的業(yè)務(wù)化應(yīng)用。

海岸帶；風(fēng)廓線雷達；大氣濕度；大氣環(huán)境；海洋監(jiān)測

1 引言

海岸帶是海陸的分界線和海陸相互作用的地帶，是人類生存和社會發(fā)展的重要場所。由于海-陸-氣等因素的共同作用，海岸帶大氣海洋環(huán)境極其復(fù)雜，對人類生產(chǎn)生活有較大影響。隨著21世紀(jì)海上絲綢之路戰(zhàn)略的正式實施，深入研究海岸帶大氣海洋環(huán)境特征將為海洋開發(fā)建設(shè)做出貢獻[1-4]。尤其是其降水特征具有極強的局地性和不穩(wěn)定性，難以預(yù)測。降水產(chǎn)生的大氣背景是水汽條件，即大氣的濕度分布狀況，這是降水預(yù)報的關(guān)鍵問題之一。而目前海岸觀測站使用的常規(guī)探空方法只能反映1天2個時次的大氣濕度情況，不能實現(xiàn)實時監(jiān)測；地基GPS(Globe Position System)水汽探測系統(tǒng)可以反演監(jiān)測水汽條件的變化，但要達到一定精度需要組網(wǎng)探測，在海岸帶難以實現(xiàn)，且探測要素單一。近年來，風(fēng)廓線雷達(Wind Profile Radar)開始在國內(nèi)廣泛應(yīng)用，其主要功能是監(jiān)測高空風(fēng)場，此外將風(fēng)廓線雷達嵌入RASS(Radio Acoustic Sounding System) 同時可監(jiān)測氣溫的垂直分布，結(jié)合其他資料可以分析降水過程[5]、強對流天氣過程[6]、邊界層特征[7]等。風(fēng)廓線雷達的觀測要素除風(fēng)外，還有許多其他物理量產(chǎn)品，結(jié)合RASS的溫度產(chǎn)品可擴展其測量內(nèi)容，如濕度[8]、湍流耗散率[9]等；但目前這些伴隨產(chǎn)品的應(yīng)用不足，導(dǎo)致觀測資源的浪費，如能實現(xiàn)其對濕度的監(jiān)測，將有利于海岸帶大氣環(huán)境的研究。鑒于此，本研究利用某海濱觀測站的數(shù)據(jù)個例，研究風(fēng)廓線雷達和RASS測量濕度的理論和方法，將其應(yīng)用于海濱觀測個例的濕度計算中，并與探空資料進行比對，分析其誤差和測量的有效性，從而得到有價值的研究結(jié)論。

2 原理、方法與數(shù)據(jù)

2.1 風(fēng)廓線雷達的探測原理

在對流層區(qū)域中，空氣垂直方向上的對流混合作用很強，會產(chǎn)生大量的湍流，湍流導(dǎo)致空氣不均勻的大氣折射率，這是風(fēng)廓線雷達可以探測到大氣后向散射的理論依據(jù)。湍流的移動可認為與風(fēng)是基本一致的，風(fēng)廓線雷達通過對“湍流”反射的雷達徑向速度的測量，反演“湍流”的平均移動速度，并以此作為風(fēng)的速度。為獲得風(fēng)廓線雷達上空三維風(fēng)場信息，至少需要3個以上正交的波束。由于5波束反演風(fēng)場的精度高于3波束，大部分風(fēng)廓線雷達采用5波束，5波束指向為1個垂直指向波束和4個傾斜指向波束。傾斜波束一般為正東、正西、正南、正北，傾斜波束的天頂夾角一般在15°左右，雷達工作時按東、西、天頂、南、北的順序進行探測，1個周期探測約為5 min。5波束風(fēng)廓線雷達波束指向示意如圖1所示[10]。除三維高空風(fēng)外，風(fēng)廓線雷達測量過程中還有譜寬、Cn2和信噪比等伴隨產(chǎn)品。

圖1 風(fēng)廓線雷達5波束的幾何關(guān)系

2.2 RASS測量虛溫原理

無線電-聲探測技術(shù)(Radio Acoustic Sounding Technology)是一種新興氣象遙感技術(shù)，該技術(shù)基于聲傳播理論可反演大氣虛溫的垂直分布，與風(fēng)廓線雷達結(jié)合可拓展其探測功能。RASS的探測原理是：聲波傳播時會使周圍空氣密度發(fā)生變化，導(dǎo)致其折射率隨之發(fā)生變化，當(dāng)雷達電磁波經(jīng)過這些區(qū)域會產(chǎn)生后向散射，雷達通過接收這部分能量能測量出不同高度層聲傳播的速度，并根據(jù)聲速與虛溫的關(guān)系得到虛溫垂直分布。RASS反演大氣溫度的原理如圖2所示[11]。

圖2 RASS反演大氣溫度的原理

這些聲波的速度與大氣的虛溫存在密切的關(guān)系，其表達式為

(1)

如進行精確計算，可用式(1)將虛溫轉(zhuǎn)換為溫度，得到實際的大氣溫度廓線。

2.3 濕度反演的原理

濕度反演理論涉及的許多物理量都可用風(fēng)廓線雷達直接或間接測量得到，主要包括布維頻率、湍流耗散率、湍流結(jié)構(gòu)常數(shù)和大氣折射梯度等，需首先介紹這些物理量及其計算方法，原理與文獻[8]類似。

布維頻率N(Brunt-Vaisala頻率)是一個與大氣靜力穩(wěn)定度有關(guān)的參數(shù)，N值越大說明大氣穩(wěn)定度越大，計算式為

(2)

式中：θ為位溫；z為高度；g=9.8 m/s2為重力加速度；θ可由虛溫Tv近似計算得到。

湍流耗散率ε是表征湍流強弱的重要參數(shù)，Hocking[12]、Cohn[13]、Worthington[14]指出湍流耗散率ε與譜寬σ2之間的關(guān)系[6]，即

(3)

式中：N為布維頻率；b為1個無量綱常數(shù)，Weinstock指出b≈α1-3/2≈ 0.5，本研究b值采用0.6，與文獻[15]相同；α1為柯爾莫戈洛夫常數(shù)；譜寬σ2可由風(fēng)廓線雷達直接測量得到，描述有效照射體積內(nèi)散射粒子相對運動的劇烈程度。

大氣折射指數(shù)梯度dn/dz與大氣折射率結(jié)構(gòu)常數(shù)之間的關(guān)系，其形式為

(4)

在大氣低層，比濕的垂直分布表達式為

(5)

Γ為干絕熱遞減率，其值為0.009 8 ℃/m，上述方程的解采用位溫θ可表示

(6)

式中：z0為起始高度；q0為起始邊界比濕；θ0為起始邊界位溫。

2.4 具體方法

利用位于大連附近海濱觀測站的CFL-16風(fēng)廓線雷達和RASS數(shù)據(jù)，采用的計算方法為：先由溫度計算出布維頻率，再得到湍流耗散率，加上譜寬和湍流結(jié)構(gòu)常數(shù)計算出大氣折射指數(shù)梯度，最后根據(jù)比濕與大氣折射指數(shù)梯度的關(guān)系計算出大氣濕度廓線(圖3)。以該方法計算海岸帶大氣濕度的垂直分布情況。

圖3 反演大氣濕度廓線的計算方法

2.5 實測數(shù)據(jù)

本研究采用的雷達是航天23所研制的CFL-16型風(fēng)廓線雷達，從探測范圍上屬于對流層風(fēng)廓線雷達，配有4個RASS發(fā)生裝置，位于大連附近海濱觀測站。CFL-16風(fēng)廓線雷達是以大氣湍流為探測目標(biāo)的脈沖多普勒雷達，采用全相參相控陣體制，能夠連續(xù)提供150 m～12 km高度范圍內(nèi)(垂直分辨率可到75 m)的大氣風(fēng)場、譜寬、湍流結(jié)構(gòu)常數(shù)等氣象要素?？砂凑涨笆隼碚摵头椒ǖ玫綕穸葦?shù)據(jù)，并與附近的大連市氣象臺觀測場的標(biāo)準(zhǔn)探空數(shù)據(jù)進行誤差比對分析。大連市氣象臺觀測場與風(fēng)廓線雷達與RASS均位于海岸帶附近，距離約30 km，在允許的比對范圍內(nèi)。所選數(shù)據(jù)個例日期為2015年6月1日，但二者時間不同，無線電探空儀測量數(shù)據(jù)時間為19時15分，風(fēng)廓線雷達和RASS數(shù)據(jù)為17時38分。

3 結(jié)果分析

3.1 研究數(shù)據(jù)

從6月1日20時的地面分析圖(圖4)來看，大連地區(qū)當(dāng)天處于鞍形場中，天氣形勢相對穩(wěn)定，大氣物理特性在一定時間內(nèi)變化不大。因此，可近似地認為1738時刻風(fēng)廓線雷達+RASS探測時，氣象要素的垂直分布與1915時刻基本相同。

圖4 2015年6月1日20時地面分析

3.2 數(shù)據(jù)對比及誤差

可按前述理論和方法計算風(fēng)廓線雷達和RASS的數(shù)據(jù)，計算結(jié)果如圖5所示。布維頻率總的趨勢是隨著高度升高逐漸減小，在370 m附近出現(xiàn)極大值，為6.255×10-4℃/S2(圖5(a))；由此造成湍流耗散率隨高度變化情況與布維頻率基本相同，在295.8 m處出現(xiàn)極大值，為0.21(圖5(b))；大氣湍流結(jié)構(gòu)常數(shù)在370 m處較大(圖5(c))。可見，風(fēng)廓線雷達和RASS的數(shù)據(jù)可以反映大氣中細微結(jié)構(gòu)變化，這是常規(guī)資料無法實現(xiàn)的。

圖5 2015年6月1日1630時風(fēng)廓線雷達+ RASS測量數(shù)據(jù)的計算結(jié)果

濕度方面，2015年6月1日風(fēng)廓線雷達和RASS反演的比濕與無線電探空儀獲得的比濕的比較情況如圖6所示，可見二者吻合較好。

圖6 2015年6月1日1630時風(fēng)廓線雷達+ RASS反演的比濕

比濕垂直分布的總體趨勢是隨高度升高而下降，細微差別較為明顯，數(shù)值差異最大處位于接近1 000 m處、2 000 m附近和2 500 m附近，最大誤差值超過2 g/kg；其對應(yīng)布維頻率、湍流耗散率和湍流結(jié)構(gòu)常數(shù)在該處附近都有1個極大值出現(xiàn)，同時布維頻率、湍流耗散率、湍流結(jié)構(gòu)常數(shù)在300～400 m附近數(shù)值均達到較大值的高度，在2種探空資料中都顯示為濕度變化劇烈。濕度與這些物理量之間的關(guān)系值得深入研究，此外計算濕度采用的溫度值是RASS測量的，同時2個觀測地點和時間的差異也是造成誤差的原因。

4 結(jié)論與展望

本研究利用CEL-16型風(fēng)廓線雷達和RASS的實測數(shù)據(jù)，采用風(fēng)廓線雷達聯(lián)合RASS測量濕度的方法，計算海岸帶實測個例的大氣濕度分布。從本次海岸帶濕度個例計算結(jié)果可得到結(jié)論。

(1)風(fēng)廓線雷達聯(lián)合RASS測量濕度的方法得到的濕度結(jié)果與常規(guī)探測結(jié)果較為接近，誤差較小，可進一步進行誤差訂正研究，在業(yè)務(wù)應(yīng)用中作為常規(guī)探空手段的補充。

(2)風(fēng)廓線雷達聯(lián)合RASS的伴隨產(chǎn)品，包括布維頻率、湍流耗散率和湍流結(jié)構(gòu)常數(shù)等，與海岸帶大氣濕度有一定的相關(guān)性，可用來進一步深入研究海岸帶大氣的微觀結(jié)構(gòu)和機理。

本研究也存在不足和需要進一步深入研究的問題：首先，風(fēng)廓線雷達聯(lián)合RASS測量濕度的垂直分辨率低于探空資料，本個例應(yīng)用的CEL-16風(fēng)廓線雷達垂直分辨率為75 m，低于常規(guī)探空資料；其次，比對數(shù)據(jù)的時間地點不一致以及計算濕度應(yīng)用RASS測量溫度，是帶來誤差的原因；最后，可結(jié)合風(fēng)廓線雷達測量的三維風(fēng)場與其濕度場相結(jié)合，得到海岸帶大氣環(huán)境的綜合狀況。

總之，隨著對風(fēng)廓線雷達聯(lián)合RASS測量方法開展進一步深入研究并實現(xiàn)其對海岸帶大氣濕度的監(jiān)測，將為海岸帶天氣過程演變、海岸帶氣象預(yù)報等相關(guān)研究提供有益參考。

[1] 鄭崇偉，潘靜，孫威，等.經(jīng)略21世紀(jì)海上絲路之海洋環(huán)境特征系列研究[J].海洋開發(fā)與管理，2015，32(7)：4-9.

[2] 鄭崇偉，李崇銀，楊艷，等.巴基斯坦瓜達爾港的風(fēng)能資源評估[J].廈門大學(xué)學(xué)報(自然科學(xué)版)，2016，55(2)：210-215.

[3] 鄭崇偉，黎鑫，陳璇，等.經(jīng)略21世紀(jì)海上絲路：地理概況、氣候特征[J].海洋開發(fā)與管理，2016，33(2)；3-10.

[4] 鄭崇偉，李崇銀.中國南海島礁建設(shè)：重點島礁的風(fēng)候、波候特征分析[J].中國海洋大學(xué)學(xué)報(自然科學(xué)版)，2015，45(9)：1-6.

[5] 耿建軍,李浚河,杜佳,等.2013年北京地區(qū)一次強對流降水天氣成因分析[J].氣象與環(huán)境科學(xué),2016,39(1):52-58.

[6] 曾瑾瑜,夏麗花,陳文鍵.福建省風(fēng)廓線雷達資料在一次強對流天氣過程中的應(yīng)用[J].氣象與環(huán)境科學(xué),2014,37(3):40-48.

[7] 邱粲,王靜,王棟成,等.濟南章丘大氣邊界層風(fēng)溫場特征研究[J].氣象與環(huán)境科學(xué),2015,38(2):36-41.

[8] 孫康遠.風(fēng)廓線雷達反演比濕廓線的試驗研究[D].南京：南京信息工程大學(xué)，2013.

[9] 張彩云，董德保，翁寧泉.對風(fēng)廓線雷達的譜寬分析及湍流耗散率的變化特征初步研究[J].大氣與環(huán)境光學(xué)學(xué)報，2009，6(4):406-413.

[10] 何平.相控陣風(fēng)廓線雷達[M].北京：氣象出版社，2006.

[11] 譚劍波,彭揚.RASS聲發(fā)射系統(tǒng)的設(shè)計[J].現(xiàn)代電子，2002(2)：17-20.

[12] HOCKING W K.Measurement of turbulent eddy dissipation rates in the middle atmosphere by radar techniques：A review[J].Radio Sci，1985(20)：1403-1422.

[13] COHN S A.Radar measurements of turbulent eddy dissipation rate in the troposphere：A comparison of techniques[J].J.Atmoc.Oceanic Technol，1995(12)：85-95.

[14] WORTHINGSTON R M.Tropopausal turbulence caused by the breaking of mountain waves[J].J.Atmos.Sci.Terr.Phys，1998,60：1543-1547.

[15] CAROI,ANNF,CLAYSON.On Turbulence and Mixing in the free atmosphere inferred from high-resolution soundings[J].J.Journal of Atmospheric and Oceanic Technology，2008(6)：833-852.

Theory of Wind Profile Radar and RASS Inverting Humidity of Coast Zone

SU Shipeng1，WANG Hua2

(1.Department of Navigation，the Dalian Naval Academy，Dalian 116018，China； 2.Department of Military Ocean，the Dalian Naval Academy，Dalian 116018，China)

Humidity of coast zone affects people.According to the method of calculating Wind Profile Radar and RASS inverting humidity of coast zone，the single example of station humidity was calculated，and it was compared with WPR and RASS data for analyzing the error and the reasons.It was concluded that the two data are very close，and the method of calculating Wind Profile Radar and RASS inverting humidity of coast zone is a useful way to research atmospheric humidity of coast zone.

Coast zone，Wind Profile Radar，Atmospheric humidity，Atmospheric environment，Marine monitoring

2016-11-11；

2017-03-10

蘇軾鵬，講師，博士研究生，研究方向為海戰(zhàn)場環(huán)境建設(shè)和軍事海洋氣象保障，電子信箱：qxsspgx@163.com

P7

A

1005-9857(2017)04-0061-05