往返式探空观测资料的质量控制及不确定性分析

为了推动新型探测资料在数值预报模式中的应用,本文进行了往返式探空资料同化应用前重要的基础性研究工作。基于国内首次往返式探空观测资料,首先建立了面向业务化应用的往返式探空资料质量控制方案,通过对比和分析质量控制前后观测样本的统计特征,论证了质量控制方案的合理性,质量控制后探测要素抽样分布更为合理,要素间一致性得到提高。进而以数值天气预报高时间分辨率的模式预报场和同站址业务常规探空观测资料为参考,分析质量控制后资料的不确定性,结果表明往返探空探测精度达到了世界气象组织WMO（World Meteorological Center）规定的突破目标,部分探测要素甚至实现了理想目标,探测资料具有可用性。最后结合数值模式背景场探讨往返探空资料的可同化性,研究表明往返探空的风场观测和夜间温度观测满足变分同化系统的高斯、无偏假定,可直接同化;气压、湿度和日间温度观测在资料同化前需要开展偏差订正工作,从而更有效的发挥资料价值。本文的研究工作为今后往返探空资料在模式中的同化应用奠定了基础。     

Influence of Surface Temperature and Emissivity on AMSU-A Assimilation over Land

AMSU-A (Advanced Microwave Sounding Unit-A) measurements for channels that are sensitive to the surface over land have not been widely assimilated into numerical weather prediction (NWP) models due to complicated land surface features. In this paper, the impact of AMSU-A assimilation over land in Southwest Asia is investigated with the Weather Research and Forecasting (WRF) model. Four radiance assimilation experiments with different land-surface schemes are designed, then compared and verified against radiosonde observations and global analyses. Besides the surface emissivity calculated from the emissivity model and surface temperature from the background field in current WRF variational data assimilation (WRF-VAR) system, the surface parameters from the operational Microwave Surface and Precipitation Products System (MSPPS) are introduced to understand the influence of surface parameters on AMSU-A assimilation over land. The sensitivity of simulated brightness temperatures to different surface configurations shows that using MSPPS surface alternatives significantly improves the simulation with reduced root mean square error (RMSE) and allows more observations to be assimilated. Verifications of 24-h temperature forecasts from experiments against radiosonde observations and National Centers for Environmental Prediction (NCEP) global analyses show that the experiments using MSPPS surface alternatives generate positive impact on forecast temperatures at lower atmospheric layers, especially at 850 hPa. The spatial distribution of RMSE for forecast temperature validation indicates that the experiments using MSPPS surface temperature obviously improve forecast temperatures in the mountain areas. The preliminary study indicates that using proper surface temperature is important when assimilating lower sounding channels of AMSU-A over land.     

Wind Speed and Altitude Dependent AMDAR Observational Error and Its Impacts on Data Assimilation and Forecasting

Aircraft Meteorological Data Relay (AMDAR) observations have been widely used in numerical weather prediction (NWP) because of its high spatiotemporal resolution. The observational error of AMDAR is influenced by aircraft flight altitude and atmospheric condition. In this study, the wind speed and altitude dependent observational error of AMDAR is estimated. The statistical results show that the temperature and the observational error in wind speeds slightly decrease as altitude increases, and the observational error in wind speed increases as wind speed increases. Pseudo single AMDAR observation assimilation tests demonstrate that the wind speed and altitude dependent observational error can provide more reasonable analysis increment. Furthermore, to assess the performance of wind speed and altitude dependent observational error on data assimilation and forecasting, two-month 3-hourly cycling data assimilation and forecast experiments based on the Weather Research and Forecasting Model (WRF) and its Data Assimilation system (WRFDA) are performed for the period during 1 September-31 October, 2017. The results of the two-month 3-hourly cycling experiments indicate that new observational error improves analysis and forecast of wind field and geo-potential height, and has slight improvements on temperature. The Fractions Skill Score (FSS) of the 6-h accumulated precipitation shows that new wind speed and altitude dependent observational error leads to better precipitation forecast skill than the default observational error in the WRFDA does.     

风廓线雷达资料对华南区域模式预报的影响

设计基于GRAPES_Meso的不同试验模拟2014年3月28日-4月8日的广东前汛期降水过程,评估风廓线资料对同化和预报的影响。对资料同化后分析增量的分析表明:相比同化时仅使用自动气象站资料,风廓线雷达资料对1000 hPa到850 hPa纬向风增量均有贡献,在850 hPa,700 hPa高度以上贡献迅速减小。应用3个试验的预报结果计算探空站、风廓线雷达站预报值与观测值的11 d均方根误差发现,同化加入风廓线雷达资料对各预报要素的改善在850 hPa高度最明显,其中风速预报误差显著降低,为0.7 m·s-1。此外,风廓线雷达资料对700 hPa风速预报有一定改善,而在925 hPa高度模拟效果反而降低。通过对2014年3月30日12:00（世界时）的个例分析发现,同化加入风廓线雷达资料的风速预报均方根误差在大雨级别以上的降水落区更大,其原因还有待于进一步研究。     

多普勒雷达风廓线的反演及变分同化试验

为了将雷达风场资料更好地应用到数值预报模式中, 使用VAD方法反演多普勒雷达风廓线并处理成标准的探空资料进行变分同化试验。结果表明: VAD方法反演的风廓线与探空实况对应较好, 验证了用VAD技术反演风廓线的可行性。用GRAPES-Meso模式的三维变分同化系统对雷达风廓线资料进行同化后, 风场的初始场明显改善, 降水强度和落区预报也有不同程度的改善。其中, 对6 h降水预报的改善明显优于对24 h的预报改善。另外, 在短时强降水预报中, 雷达风场资料的同化频率和同化窗口的不同, 对降水预报的改善情况也有所差异。在个例研究中, 同化间隔为1 h的方案6 h降水预报要优于同化间隔为3 h和6 h的方案, 同化窗口为3 h的试验方案6 h降水预报要好于同化窗口为6 h的试验方案。     

BJ-RUC系统对北京夏季边界层的预报性能评估

以北京市观象台2010年8月、2011年8月每日3次 (08:00, 14:00, 20:00,北京时,下同) L波段探空秒间隔数据为实况,对BJ-RUC系统 (rapid updated cycle system for the Beijing area) 分析和预报边界层性能进行了初步评估。结果表明:BJ-RUC系统对北京地区夏季白天边界层的细致特征具有较好的预报能力,但也存在明显的系统性误差。08:00边界层偏冷; 14:00和20:00 1 km以下的边界层则显著偏暖, 边界层内明显偏湿。整体上模式对边界层内温度、湿度的预报误差均高于自由大气。该系统对北京地区边界层内早晨 (08:00) 从夜间山风向白天谷风环流过渡、午后 (14:00) 到日落后 (20:00)1500 m以下盛行西南偏南气流的日变化特征具有较强的预报能力。系统预报的14:00边界层顶高度与评估时段内实际对流边界层高度的变化趋势一致。但预报的对流边界层顶偏高,这与BJ-RUC系统采用YSU边界层参数化方案的垂直混合更强有关。     

Impact of single-point GPS integrated water vapor estimates on short-range WRF model forecasts over southern India

Specifying physically consistent and accurate initial conditions is one of the major challenges of numerical weather prediction (NWP) models. In this study, ground-based global positioning system (GPS) integrated water vapor (IWV) measurements available from the International Global Navigation Satellite Systems (GNSS) Service (IGS) station in Bangalore, India, are used to assess the impact of GPS data on NWP model forecasts over southern India. Two experiments are performed with and without assimilation of GPS-retrieved IWV observations during the Indian winter monsoon period (November–December, 2012) using a four-dimensional variational (4D-Var) data assimilation method. Assimilation of GPS data improved the model IWV analysis as well as the subsequent forecasts. There is a positive impact of ～10 % over Bangalore and nearby regions. The Weather Research and Forecasting (WRF) model-predicted 24-h surface temperature forecasts have also improved when compared with observations. Small but significant improvements were found in the rainfall forecasts compared to control experiments.     

复合翼无人机不同传感器探测大气温湿度对比北大核心CSCD

无人机为大气探测的重要平台,为克服固定翼起飞降落条件要求高和旋翼机飞行航时短的问题,中国科学院大气物理研究所中层大气和全球环境探测实验室自主研制了一款新型复合翼无人机。为检验其在近地面探测大气温湿度的能力,于2020年7月28日—8月6日及2021年8月1—6日,在内蒙古自治区正镶白旗无人机综合验证基地开展了两期无人机搭载不同传感器探测温湿度的比对试验。结果显示:机载自动站与GPS探空仪所测温度绝对偏差为2.00℃~2.35℃,系统偏差可订正;两者所测相对湿度绝对偏差为4.28%;2021年搭载维萨拉温湿探头,测量对比表明维萨拉温湿探头与GPS探空仪测量结果一致性较好,机载自动站与两者差异较大。飞行探测试验表明:长航时复合翼无人机在近地面大气层探空方面机动性强,与常规旋翼无人机相比,可获取更大垂直与水平范围的气象要素信息。     

南极探空与两套再分析资料的对比分析

基于南极18个站点探空气象观测数据对欧洲中期天气预报中心的再分析数据(ERA-Interim)和美国国家环境预报中心的再分析数据(NECP)在南极地区高层大气的适用性进行验证。结果表明:在南极上空,随着高度抬升,探空气象观测数据与两套再分析数据中四个气象要素的差值均逐渐变大,再分析数据数值愈加偏离实际观测数值。两套再分析数据的位势高度和温度与探空观测数据偏差较小;风向则和探空观测数据相差甚远;两套再分析数据的风速与探空观测数据在300 hPa偏差较大。在季节变化中,南极的春季,再分析数据中的位势高度和温度与探空观测数据相差较大,在其他季节相差相对较小。再分析数据中的风速与探空观测数据在南极的夏季相差较小。再分析数据中的风向与探空观测数据存在较大偏差,且差值没有明显的季节变化。尽管两套再分析数据都存在很大偏差,但ERA-Interim数据整体上优于NCEP数据。对比分析也表明,采用这些再分析资料作为初始条件和边界条件驱动南极区域大气模式将带来较大的误差。未来需要加强南极探空观测,改进再分析资料同化和数值模拟系统。     

Intercomparison of the performance of MM5/WRF with and without satellite data assimilation in short-range forecast applications over the Indian region

The present study is conducted to verify the short-range forecasts from mesoscale model version5 (MM5)/weather research and forecasting (WRF) model over the Indian region and to examine the impact of assimilation of quick scatterometer (QSCAT) near surface winds, spectral sensor microwave imager (SSM/I) wind speed and total precipitable water (TPW) on the forecasts by these models using their three-dimensional variational (3D-Var) data assimilation scheme for a 1-month period during July 2006. The control (without satellite data assimilation) as well as 3D-Var sensitivity experiments (with assimilating satellite data) using MM5/WRF were made for 48 h starting daily at 0000 UTC July 2006. The control run is analyzed for the intercomparison of MM5/WRF short-range forecasts and is also used as a baseline for assessing the MM5/WRF 3D-Var satellite data sensitivity experiments. As compared to the observation, the MM5 (WRF) control simulations strengthened (weakened) the cross equatorial flow over southern Arabian sea near peninsular India. The forecasts from MM5 and WRF showed a warm and moist bias at lower and upper levels with a cold bias at the middle level, which shows that the convective schemes of these models may be too active during the simulation. The forecast errors in predicted wind, temperature and humidity at different levels are lesser in WRF as compared to MM5, except the temperature prediction at lower level. The rainfall pattern and prediction skill from day 1 and day 2 forecasts by WRF is superior to MM5. The spatial distribution of forecast impact for wind, temperature, and humidity from 1-month assimilation experiments during July 2006 demonstrated that on average, for 24 and 48-h forecasts, the satellite data improved the MM5/WRF initial condition, so that model errors in predicted meteorological fields got reduced. Among the experiments, MM5/WRF wind speed prediction is most benefited from QSCAT surface wind and SSM/I TPW assimilation while temperature and humidity prediction is mostly improved due to latter. The largest improvement in MM5/WRF rainfall prediction is due to the assimilation of SSM/I TPW. The assimilation of SSM/I wind speed alone in MM5/WRF degraded the humidity and rainfall prediction. In summary the assimilation of satellite data showed similar impact on MM5/WRF prediction; largest improvement due to SSM/I TPW and degradation due to SSM/I wind speed.     

基于WRF-3DVAR同化多源融合数据对近海风模拟的改进试验

本文利用WRF模式及其3DVAR同化系统,以2008年4月20日00时—23日00时的江苏近海10 m风场为研究个例,对Quik SCAT、Wind SAT、多源测风融合数据进行同化试验,通过比较WRF-3DVAR同化系统对模拟风场初始场和预报场的改进,检验了同化不同类型资料后WRF模式对研究区域内单点及区域近地层风速的预报效果。结果表明:同化试验对初始场有改进,且对预报场的改进较FNL资料明显;不同资料对风场模拟的影响不同,同化星星、星地多源融合资料效果最佳,Quik SCAT次之,Wind SAT最差。此外,在模式分辨率一定的情况下,提高观测资料的分辨率并不一定能够改善模拟效果,资料的稀疏分辨率存在最佳选择。     

TAMDAR Observation Assimilation in WRF 3D-Var and Its Impact on Hurricane Ike (2008) Forecast

This study evaluates the impact of atmospheric observations from the Tropospheric Airborne Meteorological Data Reporting (TAMDAR) observing system on numerical weather prediction of hurricane Ike (2008) using three-dimensional data assimilation system for the Weather Research and Forecast (WRF) model (WRF 3D-Var). The TAMDAR data assimilation capability is added to WRF 3D-Var by incorporating the TAMDAR observation operator and corresponding observation processing procedure. Two 6-h cycling data assimilation and forecast experiments are conducted. Track and intensity forecasts are verified against the best track data from the National Hurricane Center. The results show that, on average, assimilating TAMDAR observations has a positive impact on the forecasts of hurricane Ike. The TAMDAR data assimilation reduces the track errors by about 30 km for 72-h forecasts. Improvements in intensity forecasts are also seen after four 6-h data assimilation cycles. Diagnostics show that assimilation of TAMDAR data improves subtropical ridge and steering flow in regions along Ike’s track, resulting in better forecasts.     

时间加密探空观测对分析和预报质量的敏感性研究

探空观测是气象资料同化中最基本的常规观测资料,对同化分析和预报的有效改善具有重要作用。由于现有探空观测站的空间分辨率较低,分布不均匀,且每日仅有两次观测,数量偏少,限制了其分析场对中小尺度大气状态的准确再现能力。自我国L波段雷达-数字探空仪更新换代以来,探空观测具备了获取每日4次、垂直分辨为秒级和分钟级的大气廓线资料。本文利用WRF中尺度数值模式,通过06时（世界时,下同）加密探空资料和12时常规探空资料的有效同化,研究分析了时间加密探空观测资料对同化分析和预报质量的敏感性影响。结果表明：同化06时的时间加密探空资料的午后暴雨预报质量优于12时常规探空观测。具体而言,同化06时的时间加密探空资料预报的大雨和暴雨的预报技巧高于12时常规探空资料;位势高度、温度和风场等预报场的均方根误差在高层的改进效果更加明显;06时的时间加密探空资料的同化对高层的高空急流和低层的水汽通量散度的预报质量贡献更大。批量试验进一步证实了有效同化时间加密探空资料对分析和数值预报效果改进的积极意义。     

Simulation of heavy rainfall events over Indian monsoon region using WRF-3DVAR data assimilation system

We present the results of the impact of the 3D variational data assimilation (3DVAR) system within the Weather Research and Forecasting (WRF) model to simulate three heavy rainfall events (25–28 June 2005, 29–31 July 2004, and 7–9 August 2002) over the Indian monsoon region. For each event, two numerical experiments were performed. In the first experiment, namely the control simulation (CNTL), the low-resolution global analyses are used as the initial and boundary conditions of the model. In the second experiment (3DV-ANA), the model integration was carried out by inserting additional observations in the model’s initial conditions using the 3DVAR scheme. The 3DVAR used surface weather stations, buoy, ship, radiosonde/rawinsonde, and satellite (oceanic surface wind, cloud motion wind, and cloud top temperature) observations obtained from the India Meteorological Department (IMD). After the successful inclusion of additional observational data using the 3DVAR data assimilation technique, the resulting reanalysis was able to successfully reproduce the structure of convective organization as well as prominent synoptic features associated with the mid-tropospheric cyclones (MTC). The location and intensity of the MTC were better simulated in the 3DV-ANA as compared to the CNTL. The results demonstrate that the improved initial conditions of the mesoscale model using 3DVAR enhanced the location and amount of rainfall over the Indian monsoon region. Model verification and statistical skill were assessed with the help of available upper-air sounding data. The objective verification further highlighted the efficiency of the data assimilation system. The improvements in the 3DVAR run are uniformly better as compared to the CNTL run for all the three cases. The mesoscale 3DVAR data assimilation system is not operational in the weather forecasting centers in India and a significant finding in this study is that the assimilation of Indian conventional and non-conventional observation datasets into numerical weather forecast models can help improve the simulation accuracy of meso-convective activities over the Indian monsoon region. Results from the control experiments also highlight that weather and regional climate model simulations with coarse analysis have high uncertainty in simulating heavy rain events over the Indian monsoon region and assimilation approaches, such as the 3DVAR can help reduce this uncertainty.     

自动站资料在WRF 3DVAR中的同化敏感性试验

为了了解国家级地面自动观测站不同观测要素对数值预报作用,利用WRF三维变分同化系统对自动站资料的不同观测要素开展同化＆预报试验。两个月的数值预报检验结果表明：1）只同化自动站温度观测对位势高度、温度和相对湿度场预报的影响作用最大,且主要为负作用,而同化气压或风向风速观测对上述要素预报影响作用比同化温度小,但会改善位势高度预报和部分温度预报,对相对湿度预报改善作用不明显;2）对于降水预报,同化风向风速方案的评分结果最差,同化温度的方案次之,同化气压的方案可以改善24h降水预报结果。由此看出用WRF 3DVAR同化地面观测资料时,气压观测十分重要,同时也说明要想用好地面观测资料,让其为数值预报提高发挥最大作用仍需开展很多研究工作。     

基于西南涡加密探空资料同化的一次奇异路径耦合低涡大暴雨数值模拟研究

利用WRF模式及WRFDA同化系统,引入业务探空资料和西南涡加密探空资料,对一次四川盆地奇异路径低涡耦合大暴雨过程进行了数值试验,对比检验不同同化试验对本次过程降水和低涡移动路径的模拟能力,分析了加密探空资料同化对西南涡结构及其降水演变的影响。结果表明:在同化业务探空资料的基础上,引入西南涡加密探空资料能改善模式对本次降水和低涡移动路径的模拟,而仅同化业务探空资料对模拟结果的改善作用有限;引入西南涡加密探空资料,一方面能在初始风场上产生气旋式扰动,增加初始高原涡和西南涡的强度,另一方面通过调整初始四川盆地上空大气温、湿度结构,使模式在积分初期就能产生出实况量级的降水;西南涡加密探空资料的同化试验揭示了仅靠高层的高位涡不足以激发和维持700 hPa的西南涡,需要通过低层水平辐合引起正涡度增加并向上输送来增强700 hPa的气旋式环流,进而促进西南涡的移动和发展,而模拟初期降水的潜热释放也起重要作用,加深了对西南涡及其降水成因的认识。      

MWHS/FY-3资料同化在四川盆地暴雨预报中的应用研究

为了研究同化风云三号B星(FY-3B)和C星(FY-3C)的微波湿度计(MWHS及MWHS-2)观测资料在四川暴雨数值预报中的影响,本文基于Weather Research and Forecasting Model(WRF)及其三维变分同化系统Weather Research Forecast Variatinal Data Assimilation System(WRFDA),实现了对MWHS/FY-3B和MWHS-2/FY-3C观测资料的直接同化。针对2018年7月的一次四川盆地区域性暴雨过程的同化试验结果表明:同化风云三号系列卫星的微波湿度计观测资料对试验开始时刻均有改善,对相对湿度和矢量风场等物理量场有一定的正向调整作用,尤其是同化MWHS-2/FY-3C资料对风场的调整较为明显。同化试验对龙门山北部降水有较明显的改善作用,改善了降水的分布与落区,其中同化MWHS/FY-3B对盆地中部到东北部的降水量级的预报更接近实况,雨区更为连续。同化试验证明了同化风云三号系列卫星的微波湿度计观测资料对于四川盆地暴雨数值预报有一定的业务应用价值。      

GPS/PWV资料同化在强降水过程中的定量作用评估

基于WRF模式及其三维变分同化系统,分别采用江苏GPS观测网的大气可降水量资料及区域内探空和地面气象站资料,对2011年8月25日江苏苏南地区的一次区域性暴雨和2008年苏皖地区的一次局地特大暴雨过程进行了同化试验,用以比较分析GPS/PWV、探空和地面观测资料同化对强降水预报的定量作用。结果表明:探空和地面资料的同化通过对模式动力和热力场的影响,在强降水中心附近形成了强烈的辐合上升运动和热力不稳定条件,直接改进了强降水中心分布结构和强度特征,对数值模拟的成功与否起着决定性作用。而GPS/PWV在探空和地面资料同化的基础上,将使水汽条件得到增强且更有组织性,无论对降水中心强度还是位置都具有较为显著的改进作用。     

Analysis of Wind Power Assessment Based on the WRF Model

Assessing wind energy is a key step in selecting a site for a wind farm. The accuracy of the assessment is essential for the future operation of the wind farm. There are two main methods for assessing wind power： one is based on observational data and the other relies on mesoscale numerical weather prediction（NWP）. In this study, the wind power of the Liaoning coastal wind farm was evaluated using observations from an anemometer tower and simulations by the Weather Research and Forecasting（WRF） model, to see whether the WRF model can produce a valid assessment of the wind power and whether the downscaling process can provide a better evaluation. The paper presents long-term wind data analysis in terms of annual, seasonal, and diurnal variations at the wind farm, which is located on the east coast of Liaoning Province. The results showed that, in spring and summer, the wind speed, wind direction, wind power density, and other main indicators were consistent between the two methods. However, the values of these parameters from the WRF model were significantly higher than the observations from the anemometer tower. Therefore, the causes of the differences between the two methods were further analyzed. There was much more deviation in the original material, National Centers for Environmental Prediction（NCEP） final（FNL） Operational Global Analysis data, in autumn and winter than in spring and summer. As the region is vulnerable to cold-air outbreaks and windy weather in autumn and winter, and the model usually forecasted stronger high or low systems with a longer duration, the predicted wind speed from the WRF model was too large.     

The impact of the assimilation of scatterometer winds on surface wind and wave forecasts

Recent work has demonstrated that surface marine winds from the Bureau of Meteorology's operational Numerical Weather Prediction (NWP) systems are typically underestimated by 5 to 10%. This is likely to cause significant bias in modelled wave fields that are forced by these winds. A simple statistical adjustment of the wind components is shown to reduce the observed bias in Significant Wave Height considerably. The impact of increasing the vertical resolution of the NWP model and assimilating scatterometer data into the model is assessed by comparing the resulting forecast wind and waves to observations. It is found that, in general, the inclusion of scatterometer observations improves the accuracy of the surface wind forecasts. However, most of the improvement is shown to arise from the increased number of vertical levels in the atmospheric model, rather than directly from the use of the observations. When the wave model is forced with surface winds from the NWP model that includes scatterometer data, it is found that the scatterometer assimilation does not reduce the systematic bias in surface wave forecasts, but that the random errors are reduced.