基于MM5V3的大气波导预报可靠性和准确度分析

郝晓静 李清亮 郭立新 韩杰 郭相明

(1.西安电子科技大学物理与光电工程学院,西安 710071; 2. 中国电波传播研究所 电波环境特性及模化技术重点实验室,青岛 266107)

引 言

大气波导是影响电波传播的一种反常大气环境. 大气波导的存在会导致超视距现象和雷达盲区的出现,对雷达、通信等电子信息系统性能产生重要的影响. 因此,大气波导的研究,尤其是大气波导的预报得到了广泛的关注．

利用数值天气预报模式进行低空大气波导模拟和预报研究是近年来大气波导研究领域的一个热点.国内外多所院校和科研机构都开展了此方面的研究.在国外,美国和欧洲从1996—2011利用中尺度数值模式对大气折射环境和大气波导进行了模拟和预报研究[1-5]. 研究结果表明,数值模式已经能够描述大气波导的存在和位置特征,但预报的波导参数(高度、强度等)不准确. 在国内,也有很多院校和科研机构开展了这方面的研究工作[6-12],主要集中在初步实现了中尺度模式预报大气波导及在一些典型天气系统影响下大气波导生消变化机理的个例研究. 分析国内外研究进展发现,利用中尺度数值模式预报大气波导虽然能够描述波导的存在和位置特征,但对波导高度和强度的预报准确度较差. 虽然基于中尺度数值模式预报大气环境还存在许多问题尚待解决,但它仍是预报三维对流层大气折射环境和大气波导的有效手段.

本文利用全球电信系统 (Global Telecommunications System,GTS)探空数据计算的波导参数作为真实值,验证了第五代大气中尺度模式(Mesoscale Model 5 Version 3, MM5V3)预报表面波导和悬空波导的有效性;利用美国环境预报中心最终分析 (Final Operational Analysis,FNL)数据和美国海军研究生学院(Naval Postgraduate School, NPS)蒸发波导模型(下文简称“NPS模型”)计算的蒸发波导参数作为真实值,对MM5V3中尺度模式预报的蒸发波导进行验证. 验证所在区域为东海关注区、南海关注区和亚丁湾关注区．验证中尺度数值模式在这些关注区预报大气波导的可行性及其存在的问题,为日后提高数值模式预报大气波导可靠性和准确度的研究工作提供了参考.

1 MM5V3预报平台建立

MM5V3中尺度模式是美国宾夕法尼亚州立大学和美国国家大气研究中心从20世纪80年代以来共同开发的第5代区域中尺度数值模式的第三个版本. 本文模式初始场和侧边界条件选用美国环境预报中心得到的全球预报系统(Global Forecasting System, GFS)数据,实现对大气波导未来24小时的预报．利用MM5V3模式进行波导预报的参数化方案如表1所示.

表1 MM5V3模式设置

2 验证数据与方法

2.1 验证数据

用于验证数值模式预报结果的数据有两种.第一,2014年每天08:00时、20:00时两个时次的GTS探空数据. GTS是世界气象组为了迅速、准确地在全世界范围内传递气象观测数据而设立的一个数据通信系统. 本文使用的是由GTS提供的全球探空数据(下文简称“GTS数据”). 用GTS数据对预报的表面波导、悬空波导进行验证. 第二,2014年每天08:00时、20:00时两个时次的FNL最终分析数据. FNL最终分析数据是由美国环境预报中心提供的全球气象环境格点场数据,水平分辨率为1°×1°．因为FNL最终分析数据加入了高空探测等各种水文气象观测数据,包括卫星数据等,且已有大量研究证实,在缺少实测数据的情况下,FNL最终分析数据可以作为实测数据来使用[13]．因此,在缺少海洋实测数据的情况下,本文用FNL最终数据(下文简称“FNL数据”)对蒸发波导的预报效果进行验证. 预报区域中选择离验证点最近的格点进行验证,验证点和格点的距离<0.05°,大约为5 km左右,波导的水平范围一般在几km到几百km,所以,验证点基本能用来对预报结果进行验证．

2.2 诊断和验证方法

2.2.1 GTS数据诊断表面和悬空波导方法

利用GTS数据的温、压、湿参数计算修正折射指数(以下简称M)剖面,然后通过分析M的梯度判断波导及其参数. 具体方法和步骤如下.

某高度大气修正折射指数M可表示为

(1)

式中:P为大气压强;T为大气温度;e为水汽压;Z为该层大气的高度. 这些参数都可以直接或间接地从GTS数据中获取. 当修正折射率梯度满足条件

(2)

ΔM=Mmax-Mmin.

(3)

2.2.2 FNL数据诊断蒸发波导方法

目前,通常基于近海面测量的水文气象参数,利用蒸发波导模型获取蒸发波导高度. 国内外存在多个蒸发波导模型,PJ(Paulus-Jeske,)模型[14]、MGB(Musson-Genon-Gauthier-Bruth)模型[15]、BYC(Babin-Young-Carton,)模型[16]、NPS模型[17]、俄罗斯国家水文气象大学(Russian State Hydrometeorological University,RSHMU)模型[18]等. 虽然所有模型均基于莫宁-奥布霍夫相似理论,但美国的Babin[17]和乌克兰的V. K. Ivanov[18]等对不同模型进行了理论分析和试验验证,均推荐了采用NPS模型进行蒸发波导的预测. NPS模型使用海面上一定高度或不同高度上大气温、压、湿、风及海表温度作为输入,基于莫宁-奥布霍夫相似理论,获得大气温、压、湿剖面,然后计算出M剖面,把M最小值所在的高度确定为波导高度.

近地层内温度T、比湿q的垂直剖面如下式:

(4)

(5)

式中:T(z)、q(z)分别为高度z处大气温度和比湿;Tsea、qsea分别为海表温度和比湿;θ*和q*分别为位温θ和比湿q的特征尺度;k为卡曼常数;zot为温度粗糙度高度;ψθ、Γd分别为温度普适函数和干绝热递减率;L为相似长度．

在NPS模型中,气压剖面可以通过联立流体静力学方程和理想气体定律得到:

(6)

式中:p(z1)、p(z2)分别为测量高度z1、z2处的气压;TM为高度z1、z2处的虚温平均值．联立式(4)～(6)即可求得M剖面,剖面最小值对应的高度为波导高度.

2.2.3 验证方法

为了对预报结果的可靠性进行验证,分别计算一般强度波导(下文简称“一般波导”)预报可靠度、漏报率、虚警率和强波导预报可靠度. 这里需要强调,蒸发波导不区分波导强度,只分析波导预报可靠度. 为了对预报结果的准确度进行验证,分别计算了波导顶高和波导强度的均方根误差. 具体方法步骤为:

1)对于表面、悬空波导,把GTS数据同预报值都有效的数据选入样本;对于蒸发波导,把FNL数据同预报值都有效的数据选入样本,并都用样本总数TS记录.

2)把实测值和预报值结果一致(实测值和预报值都存在波导或都不存在波导)的次数记为准确次数TA.

3)统计一般波导预报可靠度RA, 用于验证一般强度波导预报有无的可靠性.可靠度RA为准确次数TA与样本总数TS之比,即

(7)

4)统计漏报率RM. 当实测值存在波导,而预报没有波导时,即为漏报．漏报率RM为漏报次数TM与实测存在波导次数TSY之比,即

(8)

5)统计虚警率RF. 当预报存在波导,而实测值没有波导时,即为虚警,也称错报．虚警率RF为虚警次数TF与实测没有波导次数TSN之比,即

(9)

6)统计强波导预报可靠度RAS. 当波导强度大于20 M时,即为强波导,强波导预报可靠度RAS为强波导次数TAS与实测存在强波导次数TSS之比,即

(10)

7)统计预报结果的准确度

当实测值与预报值波导都存在时,通过计算预报值与实测值之间的偏差,验证预报的准确度.准确度的统计量为均方根误差,分别统计波导顶高(用h表示)的均方根误差ERMS(h)和波导强度(用s表示)的均方根误差ERMS(s),表达式如下:

(11)

(12)

式中:hfi和hsi分别为波导顶高预报值和实测值;sfi和ssi分别为波导强度预报值和实测值;n为实测结果与预报结果都存在波导的次数.

3 验证结果分析

3.1 表面波导和悬空波导验证

3.1.1 东海关注区验证结果

模式预报区域为图1(a)所在区域,分别选取47909和47945两个站点为代表进行统计验证,结果见表2．由表2可知:两个站点一般波导预报可靠度均大于70%;强波导预报可靠度均大于80%;漏报率分别为40%和39%;虚警率分别为19%和26%;波导顶高的均方根误差约为500 m;波导强度均方根误差为12 M.

图1(b)和(c)为波导顶高预报值和实测值的散点图.图中显示:47909预报值高于实测值;47945预报值与实测值较接近,波导顶高所在的高度层较一致.

表2 东海关注区波导预报可靠性和准确度

(a) 站点分布

(b) 站点47909

(c) 站点47945

3.1.2 南海关注区验证结果

模式预报区域为图2(a)所在区域,分别选取48855和96471两个站点为代表进行统计验证,结果见表3．由表3可知:两个站点一般波导预报可靠度均大于80%;强波导预报可靠度较一般波导高;漏报率分别为39%和24%;虚警率分别为14%和7%;波导顶高均方根误差分别为719 m和767m;波导强度均方根误差分别为13 M和9 M.

图2(b)和(c)为波导顶高预报值和实测值的散点图.图中显示:48855实测值高于预报值;96471样本较少. 两个站点有一个共同趋势,预报值低于实测值.

表3 南海关注区波导预报可靠性和准确度

(a) 站点分布

(b) 站点48855

(c) 站点96471图2 南海关注区站点分布以及波导顶高预报值与实测值散点图

3.1.3 亚丁湾关注区验证结果

模式预报区域为图3(a)所在区域,分别选取43333和43369两个站点为代表进行统计验证. 从表4可以看出:两个站点一般波导预报可靠度均大于60%;漏报率分别为23%和15%;虚警率分别为34%和40%;43333站没有强波导样本,43369站的强波导预报可靠度为100%;波导顶高均方根误差分别为887 m和511 m;波导强度均方根误差分别为10 M和19 M.

图3(b)和(c)为波导顶高预报值和实测值散点图,图中显示:43333站样本较少,预报值低于实测值;43369站大部分预报值与实测值波导顶高在相同高度层内.

表4 亚丁湾关注区波导预报可靠性和准确度

(a) 站点分布

(b) 站点43333

(c) 站点43369图3 亚丁湾关注区站点分布以及波导顶高预报值与实测值散点图

3.2 蒸发波导验证

3.2.1 东海关注区验证结果

模式预报区域为图4(a)所在区域,分别选取两个区域进行统计验证:区域1为125°E～127°E,27°N～30°N;区域2为128°E～133°E,21°N～25°N,结果见表5．由表5可知:两个区域波导预报可靠度均大于95%;区域1的漏报率为5%,区域2无漏报;区域1无虚警,区域2没有实测数据无波导的情况,没有虚警样本,所以没有得到虚警率;两个区域波导顶高的均方根误差分别为2 m和3 m.

图4(b)和(c)为波导顶高预报值和实测值的散点图,图中显示,大部分样本预报值和实测值波导顶高值很接近,说明验证期内该海域蒸发波导高度预报效果较好.

表5 东海关注区蒸发波导预报可靠性和准确度

(a) 站点分布

(b) 区域1

(c) 区域2图4 东海关注区站点分布以及波导顶高预报值与实测值散点图

3.2.2 南海关注区验证结果

模式预报区域为图5(a)所在区域,选取经纬度范围为112°E～117°E,10°N～20°N的区域进行统计验证,结果见表6．由表6可知:预报可靠度为100%;无漏报;没有虚警样本;波导顶高的均方根误差为2 m.

图5(b)为波导顶高预报值和实测值散点图.图中显示:大部分样本预报值和实测值波导顶高值很接近,说明验证期内南海关注区蒸发波导预报效果较好.

表6 南海关注区蒸发波导预报可靠性和准确度

(a) 站点分布

(b) 南海区域图5 南海关注区站点分布以及预报值与实测值散点图

3.2.3 亚丁湾关注区验证结果

模式预报区域为图6(a)所在区域,分别选取两个区域进行统计验证:区域1为60°E～70°E,6°N～16°N;区域2为84°E～92°E,6°N～14°N,结果见表7．由表7可知:波导预报可靠度均为100%;且均无漏报;没有虚警样本,所以没有得到虚警率;波导顶高的均方根误差分别为3 m和2 m.

波导预报值和实测值散点图(图6(b)和(c))显示,验证期内,亚丁湾关注区波导预报效果较好.

表7 亚丁湾关注区蒸发波导预报可靠性和准确度

(a) 站点分布

(b) 区域1

(c) 区域2图6 亚丁湾关注区站点分布以及预报值与实测值45°散点图

3.3 验证结果总结

通过利用GTS数据对MM5V3模式预报的表面和悬空波导结果进行了验证. 验证结果表明:

1)一般波导预报可靠度在东海关注区和南海关注区约为70%～95%;而在亚丁湾关注区稍低,为60%～70%.说明MM5V3模式在东海、南海关注区的预报可靠度较亚丁湾好.

2)强波导预报可靠度在上述三个关注区都为80%以上,说明MM5V3模式在强波导的预报上可靠度较高.

3)东海、南海关注区的预报效果比较结果为:南海关注区漏报率和虚警率较东海关注区低;南海关注区和东海关注区的虚警率都较漏报率低．说明从虚警率和漏报率的角度来看,MM5V3模式在南海的预报可靠性较东海高,且两个关注区对无波导预报可靠性较有波导存在时预报可靠性高.

4)预报表面、悬空波导高度均方根误差为500～800 m,强度均方根误差为10～20 M,预报波导高度和强度的均方根误差较大,说明MM5V3模式对波导参数预报准确度较差.

更具体地说，基于本研究的数据和分析，针对在英留学生提出了四个相当重要和有用的建议，忽略了复杂的交易过程和复杂的互联网知识。

通过利用4天的FNL数据对MM5V3模式预报的蒸发波导结果进行了验证. 验证结果表明:因为选用的区域为蒸发波导高发区,且根据历史统计数据分析,验证时间段验证区域蒸发波导出现概率都为80%以上,所以模式预报和FNL数据基本都出现了蒸发波导,因此,初步得出MM5V3对蒸发波导高发区内的波导能够实现预报．

4 结 论

通过以上结果分析可得出以下结论:MM5V3 数值模式能够对强度强的表面和悬空波导做出预报,但对强度较弱的波导预报效果较差;预报虚警率低说明当预报出现波导时,实际存在波导的概率较高;但在预报准确度方面效果不好,有待提高;亚丁湾关注区因为选取的站点在陆地上,而东海关注区和南海关注区站点在岛屿上,相对岛屿上陆地的下垫面不均匀性增强,进而通过背景场进一步影响预报结果,是可能导致亚丁湾关注区预报可靠度相对较低的原因. 而且,对于不同的区域,适用的模式也不一样,以上是可能导致亚丁湾关注区预报可靠度和准确度低的原因.对于蒸发波导,MM5V3对蒸发波导高发区波导能够实现预报,但因使用的数据有限,因此MM5V3对蒸发波导的预报性能还需要进一步详细验证．

目前,中尺度数值预报模式是国内外实现大气波导预报行之有效的方法途径. 但是通过上述分析,模式对大气波导预报可靠性和准确度上还需要很大的改进. 主要考虑以下几方面:

1) 优化初始场. 中尺度气象模式的初始场是使用全球预报模式场和观测数据同化形成的,而全球模式网格粗,海上观测数据少,不能很好地描述海上大气边界层和中尺度模式的初始温度场与湿度场,可通过同化更多的海上观测数据和卫星数据改进海上初始场. 尤其在目前海上观测数据严重缺乏的情况下,同化卫星遥感数据到海上初始场,是改进海上初始场最可行的途径.

2) 中尺度模式的物理机制和参数化方案有待改进. 中尺度模式海洋大气边界层参数化方案需要准确地描述真实海洋大气边界层结构并预报大气折射条件.

[1] MARTIN O J, SEAMAN N L, STAUFFER D R, et al. A mesoscale model for EM ducting in the marine boundary layer[R]. Pennsylvania: The Pennsylvania State University, 1996.

[2] BURK S D, THOMPSON W. Mesoscal modeling of summertime refractive conditions in the southern California bight[J]. Journal of applied meteorology, 1997, 36: 522-31.

[3] 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]. Journal of applied meteorology, 2001, 40(3): 586-603.

[4] KUCAS M E. Forecasting the nighttime evolution of radiowave ducting in complex terrain using the MM5 numerical weather model[D]. Pennsylvania: Pennsylvania State University, 2003.

[5] HAACK T, WANG C G, GARRETT S, et al. Mesoscale modeling of boundary layer refractivity and atmospheric ducting[J]. Journal of applied meteorology and climatology, 2010, 49: 2437-2457.

[6] 胡晓华, 费建芳, 张翔, 等. 一次大气波导过程的数值模拟[J]. 气象科学, 2008, 28(3): 294-300.

HU X H, FEI J F, ZHANG X, et al. Numerical simulation of an atmospheric duct[J]. Scientia meteorologica sinica, 2008, 28(3): 294-300.(in Chinese)

[7] 焦林, 张永刚. 基于中尺度模式MM5下的海洋蒸发波导预报研究[J]. 气象学报, 2009, 67(3): 382-386.

JIAO L, ZHANG Y G. An evaporation duct prediction model coupled with the MM5[J]. Scientia meteorologica sinica, 2008, 67(3): 382-386. (in Chinese)

[8] 成印河, 赵振维, 何宜军, 等. 大气波导过程数值模拟研究[J]. 电波科学学报, 2009, 24(2): 260-263.

CHENG Y H, ZHAO Z W, HE Y J, et al. Numerical simulation of atmospheric duct process[J]. Chinese journal of radio science, 2009, 24(2): 260-263.(in Chinese)

[9] 王振会, 王喆, 康士峰, 等. 利用WRF模式对大气波导的数值模拟研究[J]. 电波科学学报, 2010, 25(5): 913-919.

WANG Z H, WANG Z, KANG S F, et al. Numerical simulation of atmospheric duct process with WRF model[J]. Chinese journal of radio science, 2010, 25(5): 913-919.(in Chinese)

[10] 陈莉, 高山红, 康士峰, 等. 中国近海蒸发波导的数值模拟与预报研究[J]. 中国海洋大学学报, 2011, 41(1/2): 1-8.

CHEN L, GAO S H, KANG S F, et al. Numerical modeling and forecasting of evaporation ducts over China Seas[J]. Periodical of Ocean University of China, 2011, 41(1/2): 1-8.(in Chinese)

[11] 张玉生, 韩杰, 郭相明, 等. 基于湿位涡场的悬空波导预报研究[J]. 电波科学学报, 2011, 26(3): 610-614.

ZHANG Y S, HAN J, GUO X M, et al. Elevated duct forecast based on moist potential vorticity fields[J]. Chinese journal of radio science, 2011, 26(3): 610-614.(in Chinese)

[12] 袁夏玉, 高山红, 王永明, 等. 一次海雾过程大气波导形成机理的数值研究[J]. 中国海洋大学学报, 2013, 43(1): 17-26.

YUAN X Y, GAO S H, WANG Y M, et al. Numerical modeling study on the formation mechanism of an atmospheric duct during a sea fog event[J]. Periodical of Ocean University of China, 2013, 43(1): 17-26. (in Chinese)

[13] 张蕾, 邱崇践. 利用卫星红外遥感资料反演大气温湿廓线的三维物理统计方法[J]. 红外与毫米波学报, 2008, 27(2): 109-113.

ZHANG L, QIU C J. Three-dimensional statistical-physical method of atmospheric parameters retrievals from satellite infrared measurements[J]. Journal of infrared and millimeter waves, 2008, 27(2): 109-113.(in Chinese)

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

[15] MUSSON-GENON L. A simple method to determine evaporation duct height in the sea surface boundary layer[J]. Radio science, 1992, 27(5): 635-644.

[16] BABIN S M. A new model of the oceanic evaporation duct[J]. Journal of applied meteorology, 1997, 36: 193-204.

[17] BABIN S M. LKB-Based evaporation duct model comparison with buoy data[J]. Applied meteorology, 2002, 41: 434-446.

[18] IVANOV V K. Determination of the evaporation duct height from standard meteorological data[J]. Atmospheric and oceanic physics, 2007, 43(1): 36-44.

[19] 邹立维, 周天军. 区域海气耦合模式研究进展[J]. 地球科学进展, 2012, 27(8): 857-864.

ZOU L W, ZHOU T J. A review of development and application of regional ocean-atmosphere coupled model[J]. Advances in Earth sciences, 2012, 27(8): 857-864.(in Chinese)