成印河 徐劍輝 張玉生 郭相明 游志偉

（1. 江蘇海洋大學(xué)海洋技術(shù)與測(cè)繪學(xué)院，連云港 222005；2. 江蘇海洋大學(xué)海洋遙感大數(shù)據(jù)應(yīng)用中心，連云港 222005；3. 中國(guó)電波傳播研究所 電波環(huán)境特性及?；夹g(shù)重點(diǎn)實(shí)驗(yàn)室，青島 266107）

引 言

大氣中異常大氣折射率結(jié)構(gòu)會(huì)導(dǎo)致電磁波被陷獲在一定厚度的大氣層內(nèi)，這種現(xiàn)象被稱為大氣波導(dǎo)傳播[1]. 大氣波導(dǎo)環(huán)境的存在對(duì)無(wú)線電通信、雷達(dá)探測(cè)等系統(tǒng)產(chǎn)生重要影響，導(dǎo)致超視距傳播[2]和探測(cè)盲區(qū)[3]. 大氣波導(dǎo)環(huán)境參數(shù)可以用大氣折射率垂直梯度小于零時(shí)來(lái)描述. 而大氣折射率是氣溫、大氣壓和水汽壓等物理量表示的經(jīng)驗(yàn)函數(shù)[4]，一般采用修正大氣折射率M表示，對(duì)于微波波段，通常采用以下形式：

式中：M是無(wú)量綱單位，數(shù)值上采用M單位；T為氣溫，K；P為大氣壓，hPa；e為水汽壓，hPa；h為海拔高度，m.

由式(1)可知濕度隨高度銳減或者溫度隨高度遞增，或兩者同時(shí)發(fā)生是大氣波導(dǎo)形成的必要條件，這種大氣結(jié)構(gòu)往往與大氣邊界層結(jié)構(gòu)密切相關(guān). 已有研究表明，海洋下墊面的海表面溫度(sea surface temperature, SST)分布對(duì)大氣邊界層數(shù)值模擬具有重要的影響. Doyle等[5]利用中尺度模式MM5(mesoscale model5)分析了不同分辨率的SST對(duì)大西洋海域的海洋大氣邊界層分布的影響；孟憲貴等[6]研究了夏季黃海SST對(duì)大氣邊界層的影響；霍志麗等[7]根據(jù)部分海上觀測(cè)SST分析了東海黑潮海洋鋒過(guò)程中大氣邊界層高度的變化. 在海洋大氣邊界層頂部往往伴隨著逆溫及濕度銳減的大氣結(jié)構(gòu)，容易發(fā)生大氣波導(dǎo)事件. 由此可知SST可改變海洋大氣邊界層高度及結(jié)構(gòu)，進(jìn)而影響大氣波導(dǎo)環(huán)境特征.

除了海上大氣邊界層結(jié)構(gòu)變化會(huì)直接影響到海上大氣波導(dǎo)環(huán)境，實(shí)際上還有多種因素的影響如地形、風(fēng)場(chǎng)、海岸線、SST等. 劉桂艷[8]利用中尺度天氣預(yù)報(bào)模式天氣研究與預(yù)報(bào)(weather research and forecaoting, WRF)分析了海陸風(fēng)過(guò)程中發(fā)生的大氣波導(dǎo)；丁軒茹等[9]研究表明夏季風(fēng)爆發(fā)對(duì)南海蒸發(fā)波導(dǎo)產(chǎn)生影響；成印河等[10]分析了夏季風(fēng)爆發(fā)對(duì)南海南北部低空大氣波導(dǎo)的影響；Burk和Thompson[11]使用NORAPS (Navy Operational Regional Atmospheric Prediction System)模擬研究了加州南部海陸風(fēng)對(duì)海洋大氣折射率環(huán)境的影響；Brooks等[12]指出波斯灣表面波導(dǎo)厚度是沿著風(fēng)向變化的，跟離岸距離成正比；Atkinson和Zhu[13-14]使用MM5模式研究了海風(fēng)、沿岸結(jié)構(gòu)、地形等因素對(duì)大氣波導(dǎo)的影響；Reddy等[15]研究表明在熱帶近岸海域大氣波導(dǎo)環(huán)境同樣受海陸風(fēng)環(huán)流的影響. 由此可知海洋大氣折射率環(huán)境受到多種因素的影響，風(fēng)場(chǎng)是其中一個(gè)重要的因素，然而SST是否是影響大氣波導(dǎo)形成的一個(gè)重要因素呢? Haack等[16]2010年利用4種中尺度數(shù)值模式研究分析了美國(guó)東海岸的SST對(duì)大氣波導(dǎo)數(shù)值模擬的影響，但其模擬結(jié)果是否適用于南海，有待于進(jìn)一步研究.

南海是西太平洋最大的邊緣海，其大氣環(huán)境受大陸海岸線、季風(fēng)及下墊面SST分布等因素影響，已有的研究表明海洋邊界層高度存在日變化[17-18]，也是低空大氣波導(dǎo)的高發(fā)區(qū)[19-23]，究竟每種因子對(duì)南海大氣波導(dǎo)環(huán)境影響怎樣，需要有針對(duì)性的研究. 本文基于中尺度天氣數(shù)值模式WRF，擬開(kāi)展不同SST分布對(duì)南海低空大氣波導(dǎo)數(shù)值模擬影響的研究.

1 研究方法

1.1 WRF模式配置

為了開(kāi)展下墊面不同SST對(duì)南海低空大氣波導(dǎo)數(shù)值模擬的影響，本文利用中尺度天氣數(shù)值模式ARW(Advanced Research WRF)3.6版本開(kāi)展數(shù)值模擬研究. WRF模式是新一代中尺度天氣預(yù)報(bào)模式，已經(jīng)被廣泛地應(yīng)用于降水、臺(tái)風(fēng)等方面的研究. 在WRF模式中，以美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心提供的氣候預(yù)報(bào)再分析系統(tǒng)(climate forecast system reanalysis,CFSR)再分析資料作為模式的初始場(chǎng)，積分步長(zhǎng)為60 s，以(15°N, 115°E)為中心設(shè)置了四重嵌套區(qū)域，分為南、北兩部分共6個(gè)區(qū)域，具體見(jiàn)圖1. 模式垂向采用不等間距分層，底層加密，共分為42層，模擬時(shí)間為1998年5月到6月.

圖1 研究海區(qū)(表示“科學(xué)1號(hào)”和“實(shí)驗(yàn)3號(hào)”觀測(cè)站)Fig. 1 Map of the study area and the cross-shaped stars show the “Kexue 1”and “Shiyan 3” observation stations

為了區(qū)別不同下墊面SST分布對(duì)低空大氣波導(dǎo)的影響，本文數(shù)值模擬過(guò)程中采用相同的比較成熟的積云、輻射、微物理等參數(shù)化方案，具體見(jiàn)表1.

表1 WRF模式通用的參數(shù)化方案Tab. 1 Universal parameterization schemes applied in the WRF model

1.2 試驗(yàn)方案設(shè)計(jì)

在下墊面SST驅(qū)動(dòng)模式過(guò)程中，主要設(shè)計(jì)了3種數(shù)據(jù)設(shè)計(jì)試驗(yàn)方案，分析其對(duì)南海低空大氣波導(dǎo)的影響. 3種數(shù)據(jù)分別為國(guó)家環(huán)境預(yù)報(bào)中心(National Centers for Environmental Prediction, NCEP)提供的CFSR再分析資料中的SST數(shù)據(jù)、歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasts, ECMWF)提供的ERA-I (ERA-Interim)再分析資料中的SST數(shù)據(jù)和美國(guó)國(guó)家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)提供的OISST(optimum interpolation sea surface temperature)數(shù)據(jù)[24]，具體時(shí)間、空間分辨率見(jiàn)表2.同時(shí)表2給出了4組下墊面SST驅(qū)動(dòng)試驗(yàn)方案，其中：Exp-N方案僅保留通用設(shè)置，無(wú)下墊面驅(qū)動(dòng)場(chǎng)；Exp-C方案和Exp-E方案每6 h更新一次下墊面驅(qū)動(dòng)場(chǎng)SST；而Exp-O方案每24 h更新一次.

表2 SST數(shù)值試驗(yàn)方案設(shè)計(jì)Tab. 2 SST experiment schemes setting

為了研究不同SST對(duì)南海低空大氣波導(dǎo)的影響，首先我們分析不同SST的差異. 以1998年5月5日為例，我們給出了0時(shí)(UTC，下同)的3種SST數(shù)據(jù)以及ERA-I和 CFSR提供的SST在6時(shí)、12時(shí)和18時(shí)相比0時(shí)刻的變化情況，如圖2～4所示，并以此估計(jì)SST對(duì)低空大氣波導(dǎo)的影響.

圖2 1998年5月5日0時(shí)3種SST數(shù)據(jù)((a)表示ERA-I, (b)表示CFSR, (c)表示OISST)Fig. 2 Distributions of three SST datasets at 0 UTC on May 5, 1998((a) ERA-I SST, (b) CFSR SST, (c) OISST)

圖3 1998年5月5日ERA-I SST數(shù)據(jù)日變化((a)、(b)和(c)分別表示6時(shí)、12時(shí)、18時(shí)數(shù)據(jù)減掉0時(shí)的差值)Fig. 3 Diurnal variation of the ERA-I SST dataset on May 5, 1998((a)，(b) and (c) show the difference between each time(6，12,18 UTC) and 0 UTC, respectively)

圖4 1998年5月5日CFSR SST數(shù)據(jù)日變化((a)、(b)和(c)分別表示6時(shí)、12時(shí)、18時(shí)數(shù)據(jù)減掉0時(shí)的差值)Fig. 4 Diurnal variation of the CFSR SST dataset on May 5, 1998((a)，(b) and (c) show the difference between each time(6，12,18 UTC) and 0 UTC, respectively)

從圖2～4可知，SST區(qū)別非常明顯. 由于空間分辨率的不同，3種數(shù)據(jù)都能夠表現(xiàn)出SST從南到北的區(qū)域空間變化，但ERA-I給出的SST比較平滑，溫度跨度比較大，而其他場(chǎng)分辨率較高，能夠表征一些海洋現(xiàn)象細(xì)節(jié)特征. 如以(15°N, 118°E)為中心的區(qū)域均有大小不等暖渦，但也有區(qū)別，OISST 給出的暖渦范圍大，靠近呂宋島，同時(shí)南部有多個(gè)的冷渦. 在日變化方面，從低分辨率的ERA-I給出的SST幾乎無(wú)日變化情況，而CFSR給出的SST存在較為明顯的日變化情況，如南海中西部的冷暖渦變化從小變大，然后變?。煌瑫r(shí)研究區(qū)域中陸地西海岸海域存在明顯的溫差變化，如中南半島南部海岸附近6時(shí)溫差由2.8 ℃減弱為凌晨18時(shí)的?2.4 ℃以上. 由此推測(cè)這種SST分布的不同可能會(huì)對(duì)南海低空大氣波導(dǎo)模擬具有重要的影響.

1.3 探空數(shù)據(jù)與分析方法

為了對(duì)比分析不同試驗(yàn)方案的數(shù)值模擬結(jié)果，我們與1998年5月5—25日和6月5—25日南海季風(fēng)試驗(yàn)期間“實(shí)驗(yàn)3號(hào)”和“科學(xué)1號(hào)”科考船的GPS探空觀測(cè)數(shù)據(jù)作對(duì)比. 科考船分別位于南海北部和南部，見(jiàn)圖1，北部為“實(shí)驗(yàn)3號(hào)”科考船位置(20.36°N, 116.84°E)，南部為“科學(xué)1號(hào)”科考船位置(6.23°N,109.83°E). 在南海季風(fēng)試驗(yàn)期間分別采集在0、6、12和18時(shí)4個(gè)時(shí)次的氣溫、相對(duì)濕度、露點(diǎn)溫度、海拔高度等廓線數(shù)據(jù)，采樣頻率為1 s或2 s.該數(shù)據(jù)經(jīng)過(guò)質(zhì)量控制[25]后有296個(gè)探空廓線用于對(duì)比WRF模式的輸出結(jié)果，該探空站數(shù)據(jù)幾乎平均分布在南海南北兩個(gè)區(qū)域. WRF模式輸出的氣溫、比濕等以兩科考船觀測(cè)時(shí)間、空間為中心，時(shí)間上相同，空間上采用周?chē)徑逯捣椒ǐ@得匹配數(shù)據(jù)，用于開(kāi)展對(duì)比研究.

分析方法主要采用列聯(lián)表方法用于統(tǒng)計(jì)波導(dǎo)發(fā)生情況，平均誤差和標(biāo)準(zhǔn)差用于波導(dǎo)特征參數(shù)定量分析. 列聯(lián)表是一種矩陣格式的表格，用于多元分類變量的頻率分布. Haack等[26]曾利用列聯(lián)表對(duì)美國(guó)海軍COAMPS(Coupled Ocean-Atmosphere Mesoscale Prediction System)模擬結(jié)果與無(wú)線電探空儀采集的數(shù)據(jù)對(duì)比分析，研究其預(yù)報(bào)準(zhǔn)確率、虛報(bào)率等. 針對(duì)低空大氣波導(dǎo)特征參數(shù)采用平均誤差和標(biāo)準(zhǔn)差等公式進(jìn)行量化，具體見(jiàn)公式(2)和(3).

2 結(jié)果分析與討論

2.1 低空大氣波導(dǎo)統(tǒng)計(jì)分析

為分析不同SST對(duì)低空大氣波導(dǎo)發(fā)生事件的影響，利用各模擬方案的第四重嵌套輸出數(shù)據(jù)分析大氣波導(dǎo)發(fā)生情況，并與南海兩個(gè)觀測(cè)站探空數(shù)據(jù)進(jìn)行對(duì)比分析，結(jié)果見(jiàn)表3.

表3 各試驗(yàn)方案中數(shù)值模擬與觀測(cè)的低空波導(dǎo)事件統(tǒng)計(jì)Tab. 3 Ducting events statistics for numerical simulation based on each experiment scheme and GPS radiosonde data

由表3可知，各方案的模擬正確率為對(duì)角線上兩組數(shù)據(jù)之和(都無(wú)波導(dǎo)和都有波導(dǎo))與觀測(cè)總數(shù)的比值，虛報(bào)率為各方案僅數(shù)值模擬有波導(dǎo)數(shù)與觀測(cè)中無(wú)波導(dǎo)總數(shù)的比值，例如Exp-N方案模擬正確率為61.5%，即 (120+62)/296，虛報(bào)率為26.2%，即22/84.

從各方案模擬正確率數(shù)據(jù)可知：Exp-O試驗(yàn)方案最好，正確率最高為68.2%；其次是Exp-C方案(正確率為65.5%)和Exp-N方案(正確率為61.5%)，正確率最低的為Exp-E方案(正確率為58.1%)，即NOAA提供的每24 h更新一次的SST驅(qū)動(dòng)場(chǎng)最好，每6 h更新一次的ERA-I數(shù)據(jù)中SST的方案模擬最差. 從虛報(bào)率數(shù)據(jù)可知，Exp-O和Exp-C兩種試驗(yàn)方案虛報(bào)率相當(dāng)，為21%左右，其他兩種方案虛報(bào)率都超過(guò)26%. 綜合正確率和虛報(bào)率的結(jié)果，在設(shè)計(jì)的4種試驗(yàn)方案中，Exp-O方案最好，即精細(xì)化的SST對(duì)數(shù)值模擬的低空大氣波導(dǎo)事件具有重要的影響，其次是下墊面的更新周期.

為了定量地研究不同SST對(duì)低空大氣波導(dǎo)模擬參數(shù)的影響，基于式(2)和(3)，我們開(kāi)展模擬的波導(dǎo)參數(shù)與探空觀測(cè)的大氣波導(dǎo)參數(shù)的對(duì)比分析，結(jié)果見(jiàn)表4.

標(biāo)準(zhǔn)差用來(lái)表征一組數(shù)據(jù)中離散程度，值越小說(shuō)明數(shù)據(jù)越集中；平均誤差和其標(biāo)準(zhǔn)差越小，表明該方案越好. 從表4可知，在波導(dǎo)底高中Exp-O方案平均誤差與標(biāo)準(zhǔn)差最小，其次為Exp-C方案，Exp-N方案的平均誤差最大而Exp-E方案的標(biāo)準(zhǔn)差最大. 在波導(dǎo)強(qiáng)度中Exp-N方案平均誤差最大而標(biāo)準(zhǔn)差最小，Exp-E方案平均誤差最小而標(biāo)準(zhǔn)差最大，Exp-C與Exp-O方案平均誤差與標(biāo)準(zhǔn)差整體較小. 在波導(dǎo)厚度中Exp-N方案平均誤差與標(biāo)準(zhǔn)差最小，Exp-C方案平均誤差與標(biāo)準(zhǔn)差最大，而Exp-E與Exp-O方案整體較小. 綜合波導(dǎo)特征參數(shù)平均誤差與標(biāo)準(zhǔn)差的模擬結(jié)果，在設(shè)計(jì)的4種試驗(yàn)方案中，Exp-O方案最好，其模擬的波導(dǎo)底高最好，波導(dǎo)強(qiáng)度與波導(dǎo)厚度相對(duì)較好，即精細(xì)化的SST對(duì)數(shù)值模擬的低空大氣波導(dǎo)事件具有重要的影響，其次是下墊面的更新周期.

表4 數(shù)值模擬波導(dǎo)特征參數(shù)與探空數(shù)據(jù)的對(duì)比分析Tab. 4 Contrastive analysis of the duct parameter characteristics between simulated results and radiosonde data

為了直觀地了解大氣波導(dǎo)誤差存在位置及原因，以1998年5月13日6時(shí)南海南部和1998年6月8日12時(shí)南海北部為例，做出了各方案數(shù)值模擬和GPS探空觀測(cè)的修正大氣折射率、相對(duì)濕度和氣溫隨高度分布，見(jiàn)圖5～6.

圖5 1998年5月13日6時(shí)南海南部模擬和實(shí)測(cè)參數(shù)隨高度分布((a)大氣折射率, (b)相對(duì)濕度, (c)氣溫)Fig. 5 Profiles of the simulated results and observed parameters in south SCS at 06 UTC, May 13, 1998((a)atmospheric modified refractivity index, (b) relative humidity, (c) air temperature)

從圖5(a)與圖6(a)不同日期的兩種數(shù)據(jù)的修正大氣折射率隨高度分布可知，Exp-O方案模擬修正大氣折射率隨高度分布與GPS探空觀測(cè)最接近. 由于濕度銳減與逆溫是形成大氣波導(dǎo)的必要條件，從圖5(b)與圖6(b)可知，Exp-O方案模擬的相對(duì)濕度隨高度分布與GPS探空觀測(cè)整體上最接近，有明顯的濕度銳減現(xiàn)象. 從圖5(c)與圖6(c)可知，各個(gè)方案模擬的氣溫隨高度分布相對(duì)接近但是在圖5(c)中僅Exp-O方案模擬出逆溫現(xiàn)象. 綜上所述，NOAA提供的每24 h更新一次的SST驅(qū)動(dòng)場(chǎng)最好，其模擬相對(duì)濕度與氣溫的變化更準(zhǔn)確，該結(jié)果表明精細(xì)化的SST對(duì)相對(duì)濕度和氣溫?cái)?shù)值模擬的重要性，同時(shí)也影響了低空大氣波導(dǎo)發(fā)生概率以及特征參數(shù).

圖6 1998年6月8日12時(shí)南海北部模擬和實(shí)測(cè)參數(shù)隨高度分布((a)大氣折射率, (b)相對(duì)濕度, (c)氣溫)Fig. 6 Profiles of the simulated results and observed parameters in the north SCS at 12 UTC, June 8, 1998((a)atmospheric modified refractivity index, (b) relative humidity, (c) air temperature)

2.2 子、母網(wǎng)格模擬結(jié)果分析

由于嵌套網(wǎng)格層數(shù)的增加，模式的運(yùn)算時(shí)間也隨之增加. 在使用WRF模式進(jìn)行數(shù)值模擬研究時(shí)，時(shí)間成本也是需要考慮的一個(gè)重要因素. 以最優(yōu)的SST數(shù)據(jù)為驅(qū)動(dòng)場(chǎng)，利用Exp-O方案對(duì)比分析了子、母網(wǎng)格數(shù)值模擬結(jié)果，見(jiàn)表5.

表5 第三、四重嵌套網(wǎng)格模擬結(jié)果的低空波導(dǎo)統(tǒng)計(jì)Tab. 5 Ducting events statistics for the third and fourth nest grid simulation results.

由表5可知：在都有波導(dǎo)分類中，四重嵌套樣本數(shù)明顯提高，而在都無(wú)波導(dǎo)分類中變化不大，這導(dǎo)致四重嵌套網(wǎng)格模擬正確率提高了11.8%；同樣地，四重嵌套網(wǎng)格模擬虛報(bào)率降低了2.4%. 此外，從三、四重嵌套的有波導(dǎo)總計(jì)中可知，模式在四重嵌套網(wǎng)格下模擬發(fā)生低空大氣的概率為52%，比三重嵌套41.6%提高了10.4%. 表明子網(wǎng)格雖然增加了運(yùn)算時(shí)間但其提高了模擬正確率、降低了模擬虛報(bào)率并且提高了模擬發(fā)生低空大氣波導(dǎo)的概率.

為進(jìn)一步討論該方案模擬結(jié)果，對(duì)比分析三、四重嵌套網(wǎng)格下波導(dǎo)特征參數(shù)，見(jiàn)表6.

表6 第三、四重嵌套網(wǎng)格模擬波導(dǎo)特征參數(shù)統(tǒng)計(jì)Tab. 6 Statistical characteristics of duct parameters for the third and fourth nest grid simulation results

由表6可知，四重嵌套網(wǎng)格模擬的波導(dǎo)底高平均誤差與標(biāo)準(zhǔn)差大于三重嵌套而波導(dǎo)厚度平均誤差與標(biāo)準(zhǔn)差則小于三重嵌套，四重嵌套網(wǎng)格模擬的波導(dǎo)強(qiáng)度平均誤差小于三重嵌套而標(biāo)準(zhǔn)差大于三重嵌套. 表明子網(wǎng)格由于提高了模擬發(fā)生低空大氣波導(dǎo)的概率，在一定程度上提高了模擬波導(dǎo)底高的誤差，但同時(shí)也降低了模擬波導(dǎo)厚度的誤差.

3 結(jié) 論

本文基于先進(jìn)的中尺度大氣模式WRF，利用不同SST開(kāi)展了南海低空大氣波導(dǎo)數(shù)值模擬試驗(yàn)，并與GPS探空觀測(cè)數(shù)據(jù)進(jìn)行了對(duì)比分析，初步得到了以下結(jié)論：

1)在設(shè)計(jì)的4種方案中，SST精度最高的Exp-O方案模擬結(jié)果最好，其正確率為68.2%，虛報(bào)率為21.4%，波導(dǎo)底高平均誤差以及標(biāo)準(zhǔn)差最?。?/p>

2)SST分辨率對(duì)南海低空大氣波導(dǎo)的模擬影響最大，較高的分辨率可以提高模擬正確率、更準(zhǔn)確模擬波導(dǎo)特征參數(shù)，其次是SST的更新周期，在一定程度上能改進(jìn)模擬結(jié)果；

3)在最優(yōu)的SST驅(qū)動(dòng)下，子網(wǎng)格提高了模擬正確率和發(fā)生概率，降低了虛報(bào)率和波導(dǎo)強(qiáng)度、厚度誤差，但在波導(dǎo)底高上優(yōu)勢(shì)不明顯.

上述結(jié)論對(duì)南海低空大氣波導(dǎo)數(shù)值預(yù)報(bào)具有重要的科學(xué)意義和應(yīng)用價(jià)值.