The Application to Temperature Advection Retrieval Based on Wind Profile Radar Data
-
摘要: 利用北京延庆站风廓线雷达水平风廓线资料进行大气温度平流的反演,详细分析2014年11月15日夜间冷空气入侵过程,并统计分析2015年9—11月6次冷空气入侵过程,同T639L60模式的预报风场及温度平流预报产品进行对比。结果表明:在一定预报时效内 (约6~12 h),风廓线雷达获取的水平风廓线与模式给出的预报风场有较好的一致性;由风廓线雷达反演的温度平流与模式给出的温度平流量级相同,温度平流属性一致;风廓线雷达6 min完成1次垂直高度分辨率为120 m的探测,高时空分辨率使风廓线雷达获取的温度平流较T639L60模式更能反映大气温度平流的细节。Abstract: Temperature advection is a basic physical variable in the weather prediction. It directly leads to the change of atmospheric thermal structure and further brings about the alteration of atmospheric physics field, and can reflect the development of weather systems better. In recent years, the technology of wind profile radar (WPR) improves greatly, and WPR data has both high precision and high spatial-temporal resolution and could continuously provide the distribution of horizontal winds over time. Based on the concept of temperature advection and the principle of thermal wind, retrieving the temperature advection with high precision and spatial-temporal resolution is feasible by algorithms established.The temperature advection information is retrieved from WPR data at Yanqing Station of Beijing, whose results are compared to temperature advection prediction products of Global Medium Range Forecast System of T639L60 (hereinafter, T639L60 model). A cold air invasion process on 15 Nov 2014 in the night is analyzed in detail and 6 more samples occurring from Sep to Nov in 2015 are statistically analyzed for evaluation too. Results show that the horizontal wind of the WPR is highly consistent with the wind prediction data of T639L60 model within a certain amount of prediction time (about 6 to 12 h). The order of magnitude is consistent and the value is close between the retrieval of temperature advection by the WPR data and the initial data of T639L60 model. The quality of retrieving temperature advection depends on the quality of horizontal wind profile data. Temperature advection prediction products of T639L60 model and the real-time retrieval products are accordant within a period of 6-12 h depending on the temporal scale of different weather system. The derivation gradually increases as the prediction time lengthens, sometimes even reverse prediction advection attribute may appear. Temperature advection with real-time, continuous attribute and high precision can be provided through WPR observations whose time and vertical resolution is 6 min and 120 m, respectively. And some prediction results can be calculated using the measured results by linear extrapolation in a short period. Thus, the precision of numerical forecast can be improved with WPR data assimilated in the future. Preliminary results show that the real-time 3D products and forecast products of temperature advection with high precision, spatial and temporal resolution could be generated through high quality WPR network data assimilated.
-
Key words:
- wind profile radar;
- temperature advection;
- T639L60 model
-
图 4 2014年11月15日20:00 850 hPa和500 hPa T639L60模式京津冀地区温度平流
(a)850 hPa温度平流分析场,(b)500 hPa温度平流分析场,(c)850 hPa 3 h温度平流预报场,(d)500 hPa 3 h温度平流预报场 (e)850 hPa 6 h温度平流预报场,(f)500 hPa 6 h温度平流预报场,(g)850 hPa 9 h温度平流预报场, (h)500 hPa 9 h温度平流预报场
Fig. 4 Temperature advection by T639L60 model at 850 hPa and 500 hPa at 2000 BT 15 Nov 2014 in Beijing-Tianjin-Hebei region
(a) at 850 hPa in the analysis field, (b) at 500 hPa in the analysis field, (c) at 850 hPa in the forecast field after 3 h, (d) at 500 hPa in the forecast field after 3 h, (e) at 850 hPa in the forecast field after 6 h, (f) at 500 hPa in the forecast field after 6 h, (g) at 850 hPa in the forecast field after 9 h, (h) at 500 hPa in the forecast field after 9 h
图 6 北京延庆站风廓线雷达与T639L60模式的温度平流对比图
(a)2015年9月4日08:00 850 hPa,(b)2015年9月11日08:00 850 hPa,(c)2015年10月20日20:00 850 hPa (d)2015年10月6日08:00 500 hPa,(e)2015年10月15日08:00 500 hPa,(f)2015年11月17日20:00 500 hPa
Fig. 6 Comparison of temperature advection by WPR and T639L60 model at Yanqing Station of Beijing
(a)850 hPa at 0800 BT 4 Sep 2015, (b)850 hPa at 0800 BT 11 Sep 2015, (c)850 hPa at 2000 BT 20 Oct 2015, (d)500 hPa at 0800 BT 6 Oct 2015, (e)500 hPa at 0800 BT 15 Oct 2015, (f)500 hPa at 2000 BT 17 Nov 2015
表 1 700 hPa北京延庆站实时风廓线和对应预报风场
Table 1 WPR of 700 hPa real-time wind and corresponding forecasted wind at Yanqing Station of Beijing
时间 WPR T639L60 预报时效/h 风速/(m·s-1) 风向/(°) 风速/(m·s-1) 风向/(°) 2014-11-15T 20:00 8.8 247 9.6 269 0 2014-11-15T 23:00 9.3 244 8.3 278 3 2014-11-16T 02:00 9.6 261 4.7 313 6 2014-11-16T 05:00 8.3 287 3.0 333 9 表 2 850 hPa北京延庆站实时风廓线和对应预报风场
Table 2 WPR of 850 hPa real time wind and corresponding forecasted wind at Yanqing Station of Beijing
时间 WPR T639L60 预报时效/h 风速/(m·s-1) 风向/(°) 风速/(m·s-1) 风向/(°) 2014-11-15T 20:00 7.1 270 5.6 235 0 2014-11-15T 23:00 6.8 270 6.9 282 3 2014-11-16T 02:00 5.7 299 9.3 311 6 2014-11-16T 05:00 4.7 323 10.1 331 9 表 3 500 hPa北京延庆站实时温度平流和对应预报温度平流
Table 3 WPR of 500 hPa real-time temperature advection and corresponding forecasted temperature advection at Yanqing Station of Beijing
时间 WRP T639L60 预报时效/h 强度/(10-4 K·s-1) 属性 强度/(10 -4 K·s -1) 属性 2014-11-15T20:00 -2.50 冷 -2.10 冷 0 2014-11-15T23:00 0.55 暖 0.30 暖 3 2014-11-16T02:00 0.90 暖 0.90 暖 6 2014-11-16T05:00 1.50 暖 0.20 暖 9 表 4 850 hPa北京延庆站实时温度平流和对应预报温度平流
Table 4 WPR of 850 hPa real-time temperature advection and corresponding forecasted temperature advection at Yanqing Station of Beijing
时间 WRP T639L60 预报时效/h 强度/(10-4 K·s-1) 属性 强度/(10-4 K·s-1) 属性 2014-11-15T20:00 0.85 暖 0.85 暖 0 2014-11-15T23:00 0.55 暖 0.25 暖 3 2014-11-16T02:00 -0.20 冷 -0.20 冷 6 2014-11-15T05:00 0.60 暖 -0.50 冷 9 -
[1] 李厚楹, 孙承旬.相对广义温度平流在强对流天气分析和预报中的应用.高原气象, 1996, 8(3):363-368. http://www.cnki.com.cn/Article/CJFDTOTAL-GYQX603.013.htm [2] 胡明宝.风廓线雷达数据处理与应用研究.南京:南京信息工程大学, 2012:2-6. [3] 何平, 朱小燕, 阮征, 等.风廓线雷达探测降水过程的初步研究.应用气象学报, 2009, 20(4):465-470. doi: 10.11898/1001-7313.200904011 [4] 何平.相控阵风廓线雷达.北京:气象出版社, 2006. [5] 林晓萌, 何平, 黄兴友.一种抑制降水对风廓线雷达水平风干扰的方法.应用气象学报, 2015, 26(1):66-75. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20150107&flag=1 [6] 杨新蕊, 马舒庆, 吴蕾.UHF风廓线雷达降水数据判别方法的比较与评价.大气科学学报, 2010, 33(5):576-581. http://www.cnki.com.cn/Article/CJFDTOTAL-NJQX201005009.htm [7] 何平, 李柏, 吴蕾, 等.确定风廓线雷达功率谱噪声功率方法.应用气象学报, 2013, 24(3):297-303. doi: 10.11898/1001-7313.20130305 [8] Iwan Holleman.Doppler Radar Wind Profiles.KNMI, 2003:19-21. http://www.knmi.nl/kennis-en-datacentrum/publicatie/doppler-radar-wind-profiles-131b9505-42e3-4cd5-90b6-daa342ee29a5 [9] Iwan Holleman.Quality control and verification of weather radar wind profiles.J Atmos Ocean Technol, 2005, 22:1544-1550. doi: 10.1175/JTECH1781.1 [10] Paul J Neiman, Shapiro M A.Retrieving horizontal temperature gradients and advections from single-station wind profiler observations.Wea Forecasting, 1989, 4(2):222-233. doi: 10.1175/1520-0434(1989)004<0222:RHTGAA>2.0.CO;2 [11] 董丽萍.风廓线雷达探测资料评估与产品研究及应用.成都:成都信息工程大学, 2014. [12] 翟亮, 王令, 刘文军.两次降雪天气过程预报中边界层风廓线雷达资料的应用.气象科技, 2012, 40(5):783-788. http://www.cnki.com.cn/Article/CJFDTOTAL-QXKJ201205018.htm [13] 杨波, 魏东.二次开发风廓线探测资料及其在天气预报中的应用.气象科技, 2010, 38(4):413-417. http://www.cnki.com.cn/Article/CJFDTOTAL-QXKJ201004004.htm [14] 朱乾根, 林锦瑞, 寿绍文, 等.天气学原理和方法.北京:气象出版社, 2007:11-49. [15] 范东光.用测风资料预报强降水的一种方法.气象, 1980, 6(4):11-13. doi: 10.7519/j.issn.1000-0526.1980.04.004 [16] 张明席.计算温度平流的另一表达式.气象, 1981, 7(3):8. doi: 10.7519/j.issn.1000-0526.1981.03.004 [17] 李洪勣.对温度平流计算公式的探讨.气象, 1984, 10(11):17-20. doi: 10.7519/j.issn.1000-0526.1984.11.005 [18] 陶祖钰, 周小刚, 郑永光, 等.天气预报的理论基础——准地转理论概要及其业务应用.气象科技进展, 2012, 2(3):6-16. http://www.cnki.com.cn/Article/CJFDTOTAL-QXKZ201203003.htm [19] 盛裴轩, 毛节泰, 李建国, 等.大气物理学.北京:北京大学出版社, 2009:36-41. [20] 周秀骥.高等大气物理学.北京:气象出版社, 1991:226-239. [21] 杨玉华, 徐祥德, 翁永辉.北京城市边界层热岛的日变化周期模拟.应用气象学报, 2003, 14(1):61-68. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20030107&flag=1 [22] 陈炯, 王建捷.北京地区夏季边界层结构日变化的高分辨模拟对比.应用气象学报, 2006, 17(4):403-410. doi: 10.11898/1001-7313.20060403 [23] 卞林根, 陆龙骅, 程彦杰.青藏高原东南部昌都地区近地层湍流输送的观测研究.应用气象学报, 2001, 12(1):1-13. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20010101&flag=1