The Application of Radiosonde Observation Blacklisting Check to Variable Data Assimilation System
-
摘要: 针对探空观测资料使用造成的某些区域GRAPES分析场存在虚假的高、低压系统问题,该文通过对比全球探空资料的位势高度观测与NCEP分析场,统计站点中观测质量较差的时次出现频数,确定探空位势高度观测黑名单。研究表明:500 hPa在印度地区、北大西洋和南极洲附近的探空位势高度观测与NCEP分析场的均方根误差在30 gpm以上的站点较多,且位势高度观测不可靠观测比率为20%以上的站点主要集中这些区域,以上观测站均列入黑名单。文中在GRAPES全球三维变分分析场的质量控制中加入探空位势高度观测黑名单检查,通过6 h分析预报循环试验表明:探空位势高度观测黑名单检查能有效提高分析场质量,GRAPES位势高度分析场在南极洲附近和印度地区有所改善。Abstract: The radiosonde observations are very important for the numerical weather forecasting as they are more reliable data compared to other observations in the assimilation system. The radiosonde observation needs serious quality controlling before used in data assimilation, and the blacklisting check which uses monthly blacklist to delete the abnormal noise and bias data compared with background field is a significant step in observation quality monitoring. At present the radiosonde observations used in GRAPES (Global-Regional Assimilation and PrEdiction System) analysis system are scanned through preliminary checks such as climatological extreme check, consistency checks and background quality control in data assimilation system, but blacklisting check is not included. For the spurious high/low-pressure system structure of GRAPES analysis caused by the radiosonde observation in some region, the radiosonde observations are analyzed through comparing with NCEP analysis from June to August in 2007. Results show that the bias of radiosonde geopotential height observation and NCEP analysis is less than 50 gpm under 100 hPa. The bias is much bigger over 100 hPa and the radiosonde geopotential height is much smaller than NCEP analysis. Stations where the root mean square errors of the geopotential height between radiosonde observations and NCEP analysis for 500 hPa exceed 30 gpm are located in India, North Atlantic Ocean and Antarctic region. In those regions the unreliable observation ratios of some stations are above 20%, so they are listed in the blacklist. The geopotential height blacklist of radiosonde is provided by comparing observations with NCEP analysis and estimating the frequency of bad observation-quality station. The geopotential height blacklisting check of radiosonde is applied to GRAPES global 3D-var analysis system, and the cycle experiments of the analysis and forecast from 1 July to 21 July in 2009 are carried out. Two experiments are performed, Experiment 1 (control test) without blacklisting check and Experiment 2 with blacklisting check. The results indicate that the quality of geopotential height analysis of GRAPES is improved in India and Antarctic regions when considering geopotential height blacklist. The average root mean square errors of geopotential height GRAPES analysis and NCEP analysis for 500 hPa is reduced by 6 gpm in India Region and near 10 gpm in Antarctic Region after using blacklisting check in the 3D-var analysis system. In Antarctic region, the average biases and root mean square errors of geopotential height GRAPES analysis and NCEP analysis both decrease from 1000 hPa to high levels. The blacklisting check of radiosonde observation is a necessary step in the analysis system. The updated blacklist information will be used to assimilate the radiosonde observation more effectively.
-
图 4 站点42369 2007年6月1日—8月31日885~965 hPa位势高度观测与NCEP分析场偏差频数 (黑色) 分布 (灰色为观测偏差绝对值大于阈值部分的频数)
Fig. 4 The frequency (the black) of bias between geopotential height observation and NCEP analysis for 885-965 hPa on Station 42369 from 1 June to 31 August 2007 (the gray denotes the part of absolute observation bias bigger than the threshold value)
图 5 2007年6月1日—8月31日站点42369 450~550 hPa位势高度与NCEP分析场偏差的频数 (黑色) 分布 (灰色为观测偏差绝对值大于阈值部分的频数)
Fig. 5 The frequency (the black) of bias between geopotential height observation and NCEP analysis for 450-550 hPa of Station 42369 from 1 June to 31 August in 2007 (the gray denotes the part of observation bias bigger than the threshold value)
图 9 2009年7月1日12:00—7月21日12:00 500 hPa位势高度GRAPE分析场与NCEP分析场平均均方根误差 (单位:gpm) (试验1与试验2的位势高度分析场均方根误差的差别)
Fig. 9 Root mean square errors of geopotential height between GRAPES analysis and NCEP analysis for 500 hPa at 1200 UTC from 1 July to 21 July in 2009(unit:gpm) (root mean square errors of geopotential height analysis for Experiment 1 minus Experiment 2)
表 1 2009年7月10日12:00探空位势高度观测资料份数
Table 1 Radiosonde observation pieces at 1200 UTC 10 July 2009
项目 试验1 试验2 资料缺省 7607 7607 资料使用 16999 16872 预质量控制剔除 24 24 背景场检查剔除 2571 2578 黑名单检查剔除 0 120 -
[1] Eskridge R E, Luers J K, Redder C R.Unexplained discontinuity in the US radiosonde temperature data.Part Ⅰ:Troposphere.J Climate, 2003, 16:2385-2395. doi: 10.1175/2762.1 [2] Collins W G.The operational complex quality control of radiosonde heights and temperatures at the national centers for environmental prediction.Part Ⅰ:Description of the method.J Appl Meteor, 2001, 40:137-151. doi: 10.1175/1520-0450(2001)040<0137:TOCQCO>2.0.CO;2 [3] Collins W G.The operational complex quality control of radiosonde heights and temperatures at the national centers for environmental prediction.Part Ⅱ:Examples of error diagnosis and correction from operational use.J Appl Meteor, 2001, 40:152-168. doi: 10.1175/1520-0450(2001)040<0152:TOCQCO>2.0.CO;2 [4] 陶士伟, 张跃堂, 陈卫红, 等.全球观测资料质量监视评估.气象, 2006, 32(6):53-58. doi: 10.7519/j.issn.1000-0526.2006.06.009 [5] St-James J S, Laroche S.Assimilation of wind profiler data in the Canadian Meteorological Centre's analysis systems.J Atmos Oceanic Technol, 2005, 22:1181-1194. doi: 10.1175/JTECH1765.1 [6] White P W.IFS Documentation Part Ⅱ:Data Assimilation.ECMWF Research Department, IFS Documentation Cycle CY23r4, 2003. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.153.2466&rep=rep1&type=pdf [7] 张华, 薛纪善, 庄世宇, 等.GRAPES三维变分同化系统的理想试验.气象学报, 2004, 62(1):31-41. doi: 10.11676/qxxb2004.004 [8] 庄世宇, 薛纪善, 朱国富, 等.GRAPES全球三维变分同化系统——基本设计方案与理想试验.大气科学, 2005, 29(6):872-884. http://www.cnki.com.cn/Article/CJFDTOTAL-DQXK200506003.htm [9] 薛纪善, 庄世宇, 朱国富, 等.GRAPES新一代全球/区域变分同化系统研究.科学通报, 2008, 53(20):2408-2417. doi: 10.3321/j.issn:0023-074X.2008.20.003 [10] 薛纪善.新世纪初我国数值天气预报的科技创新研究.应用气象学报, 2006, 17(5):602-610. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=200605103&flag=1 [11] 马旭林, 庄照荣, 薛纪善, 等.GRAPES非静力数值预报模式的三维变分资料同化系统的发展.气象学报, 2009, 67(1):50-60. doi: 10.11676/qxxb2009.006 [12] Chen Dehui, Xue Jishan, Yang Xuesheng, et al.New generation of multi-scale NWP system (GRAPES):General scientific design.Chinese Science Bulletin, 2008, 53(22):3433-3445. http://www.cnki.com.cn/Article/CJFDTOTAL-JXTW200822003.htm [13] Zhang Renhe, Shen Xueshun.On the development of the GRA PES—A new generation of the national operational NWP system in China.Chinese Science Bulletin, 2008, 53(22):3429-3432. [14] Sun Xiaoming, Ana P B.The impact of forcing datasets on the high-resolution simulation of tropical storm Ivan (2004) in the southern Appalachians.Mon Wea Rev, 2012, 140:3300-3326. doi: 10.1175/MWR-D-11-00345.1 [15] 李长青, 闫之辉, 王瀛, 等.NCEP-FNL与T213L31同化资料对比分析中的地形差异影响.气象, 2007, 33(1):62-69. doi: 10.7519/j.issn.1000-0526.2007.01.010 [16] Marshall G J.Trends in Antarctic geopotential height and te-mperature:A comparison between radiosonde and NCEP-NCAR reanalysis data.J Climate, 2002, 15:659-674. doi: 10.1175/1520-0442(2002)015<0659:TIAGHA>2.0.CO;2 [17] 刘红亚, 薛纪善, 沈桐立, 等.探空气球漂移及其对数值预报影响的研究.应用气象学报, 2005, 16(4):518-526. doi: 10.11898/1001-7313.20050413 [18] 马颖, 姚雯, 黄炳勋.59型与L波段探空仪温度和位势高度记录对比.应用气象学报, 2010, 21(2):214-220. doi: 10.11898/1001-7313.20100211 [19] 姚雯, 马颖, 王战, 等.用数值预报场间接对比新疆两种型号探空系统.应用气象学报, 2012, 23(2):159-166. doi: 10.11898/1001-7313.20120204