郭相明，康士峰，趙棟梁，張玉生

(1.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東青島 266003；2.中國(guó)電波傳播研究所電波環(huán)境特性及?；夹g(shù)重點(diǎn)實(shí)驗(yàn)室，山東青島 266107)

海上大氣折射率的相似性研究

郭相明1,2，康士峰2，趙棟梁1，張玉生2

(1.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東青島 266003；2.中國(guó)電波傳播研究所電波環(huán)境特性及?；夹g(shù)重點(diǎn)實(shí)驗(yàn)室，山東青島 266107)

基于邊界層相似理論，通過(guò)數(shù)值模擬的方法研究了大氣折射率的相似性，數(shù)值模擬結(jié)果顯示大氣折射率在不穩(wěn)定大氣條件和近中性條件下大氣折射率近似滿足相似理論，在穩(wěn)定大氣條件下大氣折射率不滿足相似理論假設(shè)。其次通過(guò)氣象梯度塔測(cè)量的水文氣象參數(shù)計(jì)算的大氣折射率剖面和開(kāi)路渦度系統(tǒng)測(cè)量的數(shù)據(jù)計(jì)算的相似性特征量，實(shí)驗(yàn)研究了大氣折射率的相似性，研究也表明大氣折射率在穩(wěn)定大氣條件下不滿足相似性。

大氣折射率；相似理論；海洋大氣邊界層

1 引言

大氣中的微波傳播直接依賴(lài)大氣折射率的變化特征，為此通常采用大氣折射率來(lái)描述大氣無(wú)線電波傳播物理機(jī)制[1]。基于大氣折射率的垂直梯度特征可劃分為4類(lèi)，即次折射、正常折射、超折射和波導(dǎo)折射。在微波波段，大氣折射率可以表達(dá)為大氣溫度、濕度和氣壓的非線性函數(shù)[2-3]，因此可以通過(guò)氣象傳感器測(cè)量大氣溫度、濕度和氣壓的方式確定大氣折射率的結(jié)構(gòu)特征，例如氣象上常用的高空探測(cè)系統(tǒng)和低空系留探空系統(tǒng)成為探測(cè)大氣折射率的有效途徑。

從20世紀(jì)50年代邊界層相似理論的提出，對(duì)邊界層的研究具有開(kāi)創(chuàng)性的意義，相似性理論在描述大氣邊界層，特別是近地層研究與應(yīng)用中發(fā)揮了極大的作用[4-5]?；谙嗨评碚?，可以采用近地（海）面一層大氣溫度、濕度、氣壓和風(fēng)速，以及海面溫度來(lái)確定近海面的大氣溫度、濕度和風(fēng)速剖面，為近地海面大氣折射率剖面構(gòu)建提供了理論基礎(chǔ)和方法支撐，同時(shí)認(rèn)為近海面大氣壓近似靜態(tài)平衡，并滿足壓高公式，基于此提出了許多基于相似理論的近海面折射率剖面預(yù)測(cè)方法和蒸發(fā)波導(dǎo)模型[6-10]。由于大氣折射率為大氣溫度、濕度和氣壓的非線性函數(shù)，于是有研究者也認(rèn)為大氣折射率或其位大氣折射率為相似量[11]，不過(guò)從20世紀(jì)90年代以來(lái)，Hill[12]指出，大氣折射率作為相似量是有一定的條件限制，Cook等[13]通過(guò)數(shù)值模擬的方式指出，大氣位折射率僅在不穩(wěn)定和近中性條件下可以看作為相似量，在1993年，Sorbjan[14]指出當(dāng)重新定義折射率特征量的情況下，折射率方差可以看作為相似量。

長(zhǎng)期以來(lái)，雖然有學(xué)者指出大氣折射率不能看作為相似量，但在實(shí)際應(yīng)用中又往往把其看作為相似量，同時(shí)由于條件限制，雖然通過(guò)數(shù)值模擬的方式指出不為相似量，但缺乏有效的實(shí)驗(yàn)手段進(jìn)行驗(yàn)證，為此本文從數(shù)值模擬和實(shí)驗(yàn)的角度全面闡述大氣折射率的相似性。

2 大氣折射率相似性原理

2.1 大氣折射率及其擾動(dòng)

在無(wú)線電氣象學(xué)中大氣折射指的是電磁波在大氣介質(zhì)中傳播的彎曲特性，其折射程度用折射指數(shù)n來(lái)衡量，定義為電波在自由空間中的傳播速度c（光速）與介質(zhì)中的傳播速度v的比值：

地球表面大氣折射指數(shù)的正常值一般在1.000 25和1.000 4之間。對(duì)于電波傳播研究，由于數(shù)值較小，n值不方便實(shí)際應(yīng)用。因此，定義了一個(gè)與折射指數(shù)可換算的量即通常使用的大氣折射率N。在微波及以下頻段，大氣折射指數(shù)n或大氣折射率N可由常規(guī)氣象要素確定，其關(guān)系為：

式中：p為氣壓（hPa），e為水汽壓（hPa），T為氣溫（K）?？紤]水汽壓e和比濕q之間關(guān)系為：

式中：ε為氣體常數(shù)，通常取0.622。把上式代入式（2）得：

因此，大氣溫濕度和氣壓的變化或擾動(dòng)將引起大氣折射率的變化或擾動(dòng)。在固定高度，大氣折射率的擾動(dòng)N'可以表達(dá)為[15]：

忽略氣壓擾動(dòng)，則大氣折射率擾動(dòng)可以表達(dá)為[12]：

2.2 標(biāo)量相似理論和原理

相似性理論是描述大氣邊界層，特別是近地層研究的有力工具。在近地層相似理論中，采用包括相似長(zhǎng)度、摩擦速度、溫濕度的尺度參數(shù)和其它尺度參數(shù)，距離地面海面的高度，來(lái)描述邊界層大氣參數(shù)的變化特征。前蘇聯(lián)的莫寧最早在1946年采用白漢金Π定理[16]，將摩擦速度u*，感熱通量以及浮力因子進(jìn)行組合得到一個(gè)具有長(zhǎng)度量綱的特征量，現(xiàn)稱(chēng)莫寧-奧布霍夫相似長(zhǎng)度來(lái)描述的大氣的狀態(tài)，其定義為[17]：

式中：k=0.4，為卡曼常數(shù)；g為重力加速度；w′和θ′分別為垂直風(fēng)速和位溫脈動(dòng)量；u*和θ*分別為摩擦速度和位溫的尺度參數(shù)，代表近地（海）面湍流垂直輸送的強(qiáng)度；“-”代表整體平均。莫寧-奧布霍夫長(zhǎng)度反映了雷諾應(yīng)力和浮力做功的相對(duì)大小，并定義L＞0時(shí)，大氣穩(wěn)定，L＜0，不穩(wěn)定，

|L|→∞，大氣中性。隨后，1954年莫寧和奧布霍夫最終提出了完善的近地層相似理論。其中尺度參數(shù)的定義：

式中：q′大氣比濕脈動(dòng)量；q*為比濕的尺度參數(shù)。需要說(shuō)明的是，在近地面層由于位溫和溫度相差很小通常假定T*≈θ*來(lái)定義溫度特征尺度。

與溫濕度尺度參數(shù)定義相似，大氣折射指數(shù)尺度參數(shù)或大氣折射率尺度參數(shù)定義為：

在近地（海）面，風(fēng)速、溫度和濕度的無(wú)量綱化普遍廓線方程為：

同樣，通?；诖髿庹凵渎适菧囟群蜐穸群瘮?shù)的原理，假定大氣折射率同樣滿足相似理論，即

3 基于數(shù)值模擬的大氣折射率相似性研究

3.1 數(shù)值模擬方法

首先采用數(shù)值模擬的方式確定大氣折射率是否可以看作為相似量，是否滿足相似理論?；谙嗨评碚?，式（15）—（17）對(duì)應(yīng)的積分形式，即近地層內(nèi)水平風(fēng)速u(mài)、位溫θ、比濕q的垂直廓線表示為[18]：

基于位溫剖面確定的溫度剖面為：

式中：u(z)、θ(z)、q(z)分別是高度z處的水平風(fēng)速、位溫和比濕；z0t和z0q分別為溫度和濕度粗糙度；ψh和ψq分別位溫和濕度普適函數(shù)，通常取二者相等；Γd為干絕熱遞減率，約等于0.009 76 K/m；Tsea、qsea分別是海面溫度和比濕，且考慮海水鹽度對(duì)比濕的影響，qsea=0.98qsat(Tsea)、qsat(Tsea)為基于海面溫度計(jì)算的海面飽和比濕。目前主要采用海洋-大氣耦合響應(yīng)試驗(yàn)（Tropical Oceans Global Atmosphere，TOGA-Coupled Ocean Atmosphere Response Experiment，COARE）發(fā)展起來(lái)的海氣通量整體算法COARE進(jìn)行位溫θ、比濕q的特征尺度θ?、q?的計(jì)算[19]，考慮到目前海氣通量整體算法已更新到3.0版本[20]，因此本文采用最新發(fā)布的COARE3.0算法確定海表層尺度參數(shù)和海面粗糙度，進(jìn)行大氣溫濕剖面的計(jì)算。

為了確定近海面大氣折射率剖面，還要確定大氣壓剖面，通過(guò)聯(lián)立流體靜力學(xué)方程和理想氣體定律并積分得到，即[18]

式中：R為干空氣氣體常數(shù)；g為重力加速度；P(z1)測(cè)量高度z1的氣壓，P(z)為計(jì)算高度z的氣壓；為高度z1和高度z處的平均虛溫。根據(jù)式（20）—（24）確定的溫度、比濕和氣壓剖面，代入折射率公式就可以確定大氣折射率剖面，進(jìn)而可以確定不同的大氣穩(wěn)定度條件下，不同高度的大氣折射率梯度。

3.2 數(shù)值模擬結(jié)果

假定在10 m高度處測(cè)量的空氣溫度從12℃變化到21℃，間隔0.25℃，相對(duì)濕度取80%，風(fēng)速4 m/s，氣壓1 022.07 hPa，海表溫度18℃，按照上述方法分別計(jì)算大氣的穩(wěn)定度和對(duì)應(yīng)的式（18）的左端，結(jié)果見(jiàn)圖1。與圖1條件類(lèi)似，圖2給出了相對(duì)濕度從40%變化到95%時(shí)，其它模擬條件不變時(shí)，大氣折射率相似函數(shù)的變化情況。

由圖1、2可見(jiàn)，所有的計(jì)算曲線并沒(méi)有完全在一條曲線上，但在不穩(wěn)定或近中性大氣條件下，曲線近似在一條直線上，而在穩(wěn)定大氣條件下，特別是強(qiáng)穩(wěn)定條件下曲線更為發(fā)散。因此，可以初步確定大氣折射率在穩(wěn)定大氣條件下不滿足相似理論。

圖1 不同氣-海溫差下的大氣折射率相似函數(shù)ΦN變化

圖2 相對(duì)濕度變化時(shí)大氣折射率相似函數(shù)ΦN變化

4 基于實(shí)驗(yàn)數(shù)據(jù)的大氣折射率相似性研究

4.1 實(shí)驗(yàn)簡(jiǎn)介與數(shù)據(jù)處理

為研究大氣折射率的相似性，在渤海某海域開(kāi)展了實(shí)驗(yàn)[18]。觀測(cè)儀器架設(shè)某小島岸邊的鐵塔上，鐵塔塔基大部時(shí)間被海水所淹沒(méi)，在海風(fēng)作用下近似代表了附近海域大面積海洋下墊面的特性，在維薩拉海用自動(dòng)氣象站MAWS420的基礎(chǔ)上進(jìn)行了傳感器擴(kuò)展[21]，在距離塔基平臺(tái)6、10、15、20、27、34 m處架設(shè)溫濕度傳感器，在6、15、27 m處架設(shè)風(fēng)速傳感器，同時(shí)在第二層觀測(cè)平臺(tái)上（10 m）安裝了氣壓傳感器和紅外海溫傳感器，所有氣象傳感器采集的數(shù)據(jù)存儲(chǔ)在采集卡中，所有數(shù)據(jù)的采樣頻率為每分鐘一次，具有實(shí)時(shí)和連續(xù)的優(yōu)點(diǎn)。在塔基高度10 m處，在鐵塔上架設(shè)渦動(dòng)開(kāi)路協(xié)方差系統(tǒng)，為了盡可能地避免塔體對(duì)測(cè)量的影響，CSAT3三維超聲風(fēng)速儀和LI7500開(kāi)路氣體分析儀盡量遠(yuǎn)離塔體架設(shè)（距塔體約2 m）。渦動(dòng)協(xié)方差系統(tǒng)采樣頻率為10 Hz，采樣時(shí)間每次30 min?；谠加^測(cè)資料，經(jīng)過(guò)去趨勢(shì)、坐標(biāo)旋轉(zhuǎn)等處理后計(jì)算出近海層特征尺度（u*、θ?、q?）、穩(wěn)定度參數(shù)ξ等，進(jìn)而確定折射率特征尺度。

圖3 實(shí)驗(yàn)期間大氣折射率相似函數(shù)ΦN剖面隨大氣穩(wěn)定度ξ的變化

為了確定不同高度的大氣折射率梯度，基于氣象梯度塔測(cè)量的大氣參數(shù)和海面溫度，進(jìn)行同時(shí)間段30 min的平均，在靜力平衡和理想氣體假設(shè)下，利用壓高公式將10 m處氣壓換算到其他高度層上，然后利用式（2）和（3），可以計(jì)算不同高度層上的M值，假定海面相對(duì)濕度為100%，確定海面大氣折射率，然后對(duì)7層修正折射率進(jìn)行非線性最小二乘擬合，可以確定大氣修正折射率剖面所用的擬合公式為：

式中：M代表修正折射率；z代表高度；f0、f1、f2為擬合系數(shù)。在此需要指出的是在近海面大氣折射率的計(jì)算中，采用了海面為飽和水汽的假設(shè)，同時(shí)在實(shí)際的計(jì)算過(guò)程中，必須考慮潮位的影響，本文在數(shù)據(jù)處理過(guò)程中，利用潮汐數(shù)據(jù)對(duì)氣象傳感器的實(shí)際距離海面的高度進(jìn)行了訂正。根據(jù)大氣折射率和修正折射率的關(guān)系則可確定對(duì)應(yīng)的大氣折射率剖面，可確定不同高度的大氣折射率梯度值。

4.2 實(shí)驗(yàn)結(jié)果

取對(duì)應(yīng)時(shí)間段渦動(dòng)開(kāi)路協(xié)方差采集的數(shù)據(jù)和氣象梯度儀采集的數(shù)據(jù)計(jì)算大氣折射率特征量、大氣穩(wěn)定性和對(duì)應(yīng)高度的大氣折射率梯度值，計(jì)算不同穩(wěn)定度條件下大氣折射率普適函數(shù)。圖3給出了實(shí)驗(yàn)期間，海風(fēng)影響下大氣折射率普適函數(shù)隨大氣穩(wěn)定度變化情況。

由圖3可見(jiàn)，所有的曲線并沒(méi)有完全在一條曲線上，在不穩(wěn)定大氣條件下與數(shù)值模擬結(jié)果不同，表現(xiàn)出更加的分散性，這可能受實(shí)驗(yàn)條件的限制，在海岸區(qū)域受來(lái)自陸地氣流的影響，不能完全代表廣闊海域的實(shí)際海洋環(huán)境，雖然在數(shù)據(jù)處理過(guò)程中考慮了測(cè)量誤差，但是不可避免的也存在測(cè)量誤差和環(huán)境代表性問(wèn)題，使得在不穩(wěn)定大氣條件也顯示出很強(qiáng)的分散性。在穩(wěn)定大氣條件下所有的曲線相對(duì)模擬結(jié)果表現(xiàn)出更強(qiáng)的分散性。不過(guò)在近中性條件下，曲線還是表現(xiàn)出一定的匯聚趨勢(shì)。因此，也可以初步確定大氣折射率在海上實(shí)驗(yàn)環(huán)境下表現(xiàn)出不滿足相似量的假設(shè)，不滿足標(biāo)量相似理論，特別是在穩(wěn)定大氣條件下。

5 結(jié)論與討論

長(zhǎng)期以來(lái)，考慮大氣折射率為溫度、濕度和氣壓的函數(shù)，以及大氣溫濕度在近地（海）面層滿足相似理論，將其看作為相似量，并在許多方面得到了廣泛的應(yīng)用。例如折射率結(jié)構(gòu)常數(shù)的確定和近海面蒸發(fā)波導(dǎo)模型的構(gòu)建等。本文通過(guò)數(shù)值模擬和觀測(cè)實(shí)驗(yàn)兩種途徑論述了大氣折射梯度是否滿足相似理論，大氣折射率是否可以看作為相似量，基于數(shù)值模擬和實(shí)驗(yàn)結(jié)果顯示大氣折射率梯度在穩(wěn)定大氣條件下不滿足相似理論，大氣折射率不為相似量。

在此需要指出，（1）由于受實(shí)驗(yàn)條件的限制，試驗(yàn)結(jié)果顯示在不穩(wěn)定大氣條件下大氣折射率梯度也不滿足相似理論，與數(shù)值模擬結(jié)果不一致，需要開(kāi)展海上實(shí)驗(yàn)進(jìn)一步驗(yàn)證；（2）大氣折射率不僅是溫、濕度的函數(shù)，同時(shí)大氣壓強(qiáng)也是一個(gè)重要的量，在上述定義大氣折射率特征尺度量時(shí)忽略了大氣壓強(qiáng)的影響；（3）在折射率結(jié)構(gòu)常數(shù)的定義中，考慮了溫度和濕度的相關(guān)性對(duì)折射率結(jié)構(gòu)常數(shù)的影響，因此在定義大氣折射率特征尺度時(shí)考慮大氣溫度和濕度在不同的穩(wěn)定度條件下相關(guān)性對(duì)大氣折射率的影響，大氣折射率梯度有可能滿足相似理論；（4）在相似函數(shù)的定義中，僅考慮了大氣穩(wěn)定度對(duì)大氣參數(shù)梯度的影響，F(xiàn)airall等[19]指出，除了大氣穩(wěn)定度外，邊界層高度、地海面粗糙度等也是一個(gè)重要的因素，因此發(fā)展更為完善的相似理論也有可能使大氣折射率滿足相似理論。

致謝：感謝中科院大氣物理所王介民研究員編著的《渦動(dòng)相關(guān)通量觀測(cè)指導(dǎo)手冊(cè)》2012版對(duì)實(shí)驗(yàn)和數(shù)據(jù)處理的指導(dǎo)。

[1]Atkinson B W,Li J G,Plant R S.Numerical modeling of the propagation environment in the atmospheric boundary layer over the Persian Gulf[J].Journalof Applied Meteorology and Climatology,2000,40(3):586-603.

[2]Babin S M,Dockery G D.LKB-based evaporation duct model comparison with buoy data[J].Journal of Applied Meteorology and Climatology,2001,41(4):434-446.

[3]Bean B R,Dutton E J.Radio meteorology[M].Washington D.C.: United States,Government Printing Office,1966:435.

[4]Businger J A,Wyngaard J C,Izumi Y,et al.Flux-profile relationshipsin theatmospheric surface layer[J].Journalofthe Atmospheric Sciences,1971,28(2):181-189.

[5]Gryning S E,Batchvarova E,Brümmer B,et al.On the extension of the wind profile over homogeneous terrain beyond the surface boundary layer[J].Boundary-Layer Meteorology,2007,124(2): 252-268.

[6]Paulus R A.Practical application of an evaporation duct model[J]. Radio Science,1985,20(4):887-896.

[7]Babin S M,Young G S,Carton J A.A new model of the oceanic evaporation duct[J].Journal of Applied Meteorology,1997,36(3): 193-204.

[8]Ivanov V K,Shalyapin V N,Levadnyi Y V.Determination of the evaporation duct height from standard meteorological data[J]. Izvestiya,Atmospheric and Oceanic Physics,2007,43(1):36-44.

[9]Newton D A.COAMPS modeled surface layer refractivity in the roughness and evaporation duct experiment 2001[D].California: Naval Postgraduate School,2003:7-9.

[10]Ding J L,Fei J F,Huang X G,et al.Development and validation of an evaporation duct model.Part II:evaluation and improvement of stability functions[J].Journal of Meteorological Research, 2014,29(3):482-495.

[11]Gossard E E.The height distribution of refractive index structure parameter in an atmosphere being modified by spatial transition at its lower boundary[J].Radio Science,1978,13(3):489-500.

[12]Hill R J.Implications of Monin-Obukhov similarity theory for scalar quantities[J]Journal of the Atmospheric Sciences,1989,46 (14):2236-2244.

[13]Cook J,Burk S.Potential refractivity as a similarity variable[J]. Bound-Layer Meteorology,1992,58(1-2):151-159.

[14]Sorbjan Z.Monin-Obukhov similarity for refractive index revisited [J].Journal of the Atmospheric Sciences,1993,50(21):3677-3684.

[15]Gossard E E.Refractive index variance and its height distributionin different air masses[J].Radio Science,1977,12(1):89-105.

[16]郭相明,康士峰,張玉生,等.蒸發(fā)波導(dǎo)模型特征及其適用性研究[J].海洋預(yù)報(bào),2013,30(5):75-83.

[17]盛裴軒,毛節(jié)泰,李建國(guó),等.大氣物理學(xué)[M].北京:北京大學(xué)出版社,2003:247-250.

[18]郭相明,康士峰,韓杰,等.蒸發(fā)波導(dǎo)數(shù)據(jù)庫(kù)及中國(guó)海域的統(tǒng)計(jì)分析[J].電波科學(xué)學(xué)報(bào),2013,28(6):1152-1157,1182.

[19]Fairall C W,Bradley E F,Rogers D P,et al.Bulk parameterization of air-sea fluxes for tropical ocean-global atmosphere coupledocean atmosphere response experiment[J].Journal of Geophysical Research:Oceans,1996,101(C2):3747-3764.

[20]Fairall C W,Bradley E F,Hare J E,et al.Bulk parameterization of air-sea fluxes:updates and verification for the COARE algorithm [J].Journal of Climate,2003,16:571-591.

[21]范建業(yè),高山紅,康士峰,等.基于氣象梯度儀的蒸發(fā)波導(dǎo)探測(cè)[J].無(wú)線電工程,2012,42(11):32-33.

Study on the similarity of ocean atmospheric refractivity

GUO Xiang-ming1,2，KANG Shi-feng2，ZHAO Dong-liang1，ZHANG Yu-sheng2
（1.College of Oceanic and Atmospheric Sciences,Ocean University of China,Qingdao 266003 China; 2.National Key Laboratory of Electromagnetic Environment，China Research Institute of Radiowave Propagation,Qingdao 266107 China）

Based on Monin-Obukhov similarity theory,using the numerical method to validate whether the radio refractivity behave as a similarity variable,numerical results show that the vertical gradient of radio refractivity does not collapse to a single atmospheric stability universal function,but under unstable and near neutral conditions the radio refractivity approximately behave as a similarity variable.Using the hydrometeorology parameters measured by the meteorology grads tower to determine the radio refractivity profile,and with the turbulence data observed using the open path eddy covariance system,experiment results show that the vertical gradient of radio refractivity also does not collapse to a single atmospheric stability universal function when the atmosphere condition is stable.

radio refractivity;similarity theory;marine atmospheric boundary layer

P732.6

A

1003-0239（2017）03-0071-06

10.11737/j.issn.1003-0239.2017.03.009

2016-09-04；

2016-10-18。

國(guó)家自然科學(xué)基金（41175012）。

郭相明（1981-），男，工程師，博士在讀，主要從事海洋環(huán)境與電波傳播方面的研究。E-mail:xmguo08@163.com