土壤湿度对卫星辐射率资料直接同化的影响

本研究利用WRF模式及其三维变分同化系统实现了对NOAA-16 AMSU-A微波资料的直接同化,针对2010年6月19日江西地区的一次强降水过程开展模拟与同化试验,并利用中国区域土壤湿度同化系统(CLSMDAS—China Land Soil Moisture Data Assimilation System)输出的土壤湿度值替换NCEP(National Centers for Environmental Prediction)资料中的土壤湿度,研究土壤湿度初值对辐射率资料直接同化中观测场与背景场偏差调整的影响。结果表明:采用CLSMDAS输出土壤湿度初值条件下模拟的亮温值与实际观测值更为接近,经过质量控制和偏差订正后更多的观测资料能够进入到同化系统中,说明改进的土壤湿度初值条件下观测算子的计算值得到正的调整,对低层地表通道的改进效果明显,尤其以50.3 GHz的窗区通道3的结果最为理想;针对此次强降水过程中24 h累积降水分布的模拟结果,CLSMDAS输出土壤湿度初值条件下同化AMSU-A资料,能够较为准确的把握整个雨带的走向、大雨以上级别降水的落区范围、降水中心落区及强度等。说明准确的土壤湿度初值能够改进卫星辐射率资料的同化结果,进而提高数值模式的模拟预报能力。     

EnSRF雷达资料同化在一次飑线过程中的应用研究

本文利用包含复杂冰相微物理过程的WRF（Weather Research and Forecasting）模式,针对2007年4月23日发生在我国华南地区的一次典型飑线天气过程,分别进行了确定性预报和集合预报试验,发现确定性预报能大致捕捉到飑线系统的发生发展过程,但对飑线后部的层云区模拟效果较差。集合预报能够有效地减少模式的不确定性,大部分集合成员对飑线的模拟效果优于确定性预报。进一步将集合预报得到的40个成员作为背景场,采用EnSRF（Ensemble Square Root Filter）同化多普勒天气雷达资料,并将分析得到的集合作为初始场进行集合预报,通过与未同化雷达资料的集合对比,考察了EnSRF同化多部雷达资料对飑线系统的影响。结果表明:EnSRF雷达资料同化增加了模式初始场的中小尺度信息,大部分集合成员的分析场能够较准确地再现飑线的热力场、动力场和微物理场的细致特征,并且模拟出飑线后部的层云结构。通过对EnSRF分析的集合进行模拟发现,大部分集合成员较未同化雷达资料时模拟效果有明显改善。同化后的集合预报ETS（Equitable Threat Score）评分最高,其次是未同化的集合预报,确定性预报的最低。     

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.     

A comprehensive model inter-comparison study investigating the water budget during the BALTEX-PIDCAP period

Summary A comparison of 8 regional atmospheric model systems was carried out for a three-month late summer/early autumn period in 1995 over the Baltic Sea and its catchment area. All models were configured on a common grid using similar surface and lateral boundary conditions, and ran in either data assimilation mode (short term forecasts plus data assimilation), forecast mode (short term forecasts initialised daily with analyses) or climate mode (no re-initialisation of model interior during entire simulation period). Model results presented in this paper were generally post processed as daily averaged quantities, separate for land and sea areas when relevant. Post processed output was compared against available analyses or observations of cloud cover, precipitation, vertically integrated atmospheric specific humidity, runoff, surface radiation and near surface synoptic observations. The definition of a common grid and lateral forcing resulted in a high degree of agreement among the participating model results for most cases. Models operated in climate mode generally displayed slightly larger deviations from the observations than the data assimilation or forecast mode integration, but in all cases synoptic events were well captured. Correspondence to near surface synoptic quantities was good. Significant disagreement between model results was shown in particular for cloud cover and the radiative properties, average precipitation and runoff. Problems with choosing appropriate initial soil moisture conditions from a common initial soil moisture field resulted in a wide range of evaporation and sensible heat flux values during the first few weeks of the simulations, but better agreement was shown at later times. Received September 8, 2000 Revised April 3, 2001     

Assessment of Data Assimilation Approaches for the Simulation of a Monsoon Depression Over the Indian Monsoon Region

A low pressure system that formed on 21 September 2006 over eastern India/Bay of Bengal intensified into a monsoon depression resulting in copious rainfall over north-eastern and central parts of India. Four numerical experiments are performed to examine the performance of assimilation schemes in simulating this monsoon depression using the Fifth Generation Mesoscale Model (MM5). Forecasts from a base simulation (with no data assimilation), a four-dimensional data assimilation (FDDA) system, a simple surface data assimilation (SDA) system coupled with FDDA, and a flux-adjusting SDA system (FASDAS) coupled with FDDA are compared with each other and with observations. The model is initialized with Global Forecast System (GFS) forecast fields starting from 19 September 2006, with assimilation being done for the first 24 hours using conventional observations, sounding and surface data of temperature and moisture from Advanced TIROS Operational Vertical Sounder satellite and surface wind data over the ocean from QuikSCAT. Forecasts are then made from these assimilated states. In general, results indicate that the FASDAS forecast provides more realistic prognostic fields as compared to the other three forecasts. When compared with other forecasts, results indicate that the FASDAS forecast yielded lower root-mean-square (r.m.s.) errors for the pressure field and improved simulations of surface/near-surface temperature, moisture, sensible and latent heat fluxes, and potential vorticity. Heat and moisture budget analyses to assess the simulation of convection revealed that the two forecasts with the surface data assimilation (SDA and FASDAS) are superior to the base and FDDA forecasts. An important conclusion is that, even though monsoon depressions are large synoptic systems, mesoscale features including rainfall are affected by surface processes. Enhanced representation of land-surface processes provides a significant improvement in the model performance even under active monsoon conditions where the synoptic forcings are expected to be dominant.     

Impact of 4DVAR Assimilation of AIRS Total Column Ozone Observations on the Simulation of Hurricane Earl

The Atmospheric Infrared Sounder(AIRS) provides twice-daily global observations of brightness temperature, which can be used to retrieve the total column ozone with high spatial and temporal resolution.In order to apply the AIRS ozone data to numerical prediction of tropical cyclones, a four-dimensional variational(4DVAR) assimilation scheme on selected model levels is adopted and implemented in the mesoscale non-hydrostatic model MM5. Based on the correlation between total column ozone and potential vorticity(PV), the observation operator of each level is established and five levels with highest correlation coefficients are selected for the 4DVAR assimilation of the AIRS total column ozone observations. The results from the numerical experiments using the proposed assimilation scheme for Hurricane Earl show that the ozone data assimilation affects the PV distributions with more mesoscale information at high levels first and then influences those at middle and low levels through the so-called asymmetric penetration of PV anomalies.With the AIRS ozone data being assimilated, the warm core of Hurricane Earl is intensified, resulting in the improvement of other fields near the hurricane center. The track prediction is improved mainly due to adjustment of the steering flows in the assimilation experiment.     

基于WRFDA-Chem三维变分同化系统的中国PM2.5和PM10资料同化试验北大核心CSCD

将大气化学三维变分同化系统WRFDA_Chem引入睿图—化学环境气象数值预报系统(RMAPS-Chem),利用2016年11月地面观测细颗粒物(PM2.5)和颗粒物(PM10)逐小时质量浓度资料进行同化预报试验:6 h循环同化结果表明,WRFDA-Chem对初始场PM2.5和PM10的模拟偏差和相关性有显著改善,均方根误差(RMSE)减小40%左右,相关性提高0.27~0.37;同化对预报改进能持续24 h以上,PM2.5(PM10)浓度预报RMSE降低25%(10%),相关性提升14%(25%);加密同化频次(逐小时循环同化)进一步改进预报效果。未来需要进一步开展同化数据质量控制方案研究以优化业务预报效果,并在深入理解模式不确定性和偏差来源的情况下,进一步开展模式和同化系统的协同发展。     

中尺度区域模式层顶高度对高空风场数值模拟的影响分析

本文基于中尺度区域模式WRF,开展模式层顶高度变化对高空气象要素,特别是高空风场数值模拟影响的研究。通过设计模式顶高45、5 hPa两个试验,同化来源于NOAA-15、NOAA-18、NOAA-19和METOP-2的AMSU-A辐射计高通道数据,表明提高模式层顶能够使卫星更多的高通道样本数量进入同化系统,达到减小背景场误差,同时减小高于层顶通道辐射能量对低层通道影响的目的,一定程度上改进了同化效果,从而改善高空气象要素,特别是风场的模拟效果,与观测值的均方根误差减小了约0.4~0.5 m·s^-1。     

雷达反射率资料的三维变分同化研究

应用天气研究和预报模式(WRF)三维变分系统中一种新的雷达反射率资料间接同化方法来进行反射率资料的三维变分同化研究,评估雷达反射率资料对夏季短时定量降水预报的作用.该方法不直接同化雷达反射率资料,而是同化由反射率资料反演出的雨水和估计的水汽.以2009年夏季北京地区发生的4次强降水过程为例,考察了北京市气象局业务运行的快速更新循环同化预报系统对京津冀地区雷达网的雷达反射率资料的同化性能以及雷达反射率资料和径向风资料同时同化的效果.数值试验结果表明:(1)同化反演雨水或水汽都能改善降水预报,但同化反演水汽对降水预报效果的改善起了更重要的作用;(2)同化反射率资料能极大地提高短时降水预报的效果,其稳定的正面效果可以延伸到6h的预报时效,而同化径向风资料不能得到稳定的正效果;(3)同化雷达资料时,应用快速更新循环同化预报系统是提高短时定量降水预报的一个有效途径.     

ATOVS资料同化方案对暴雨模拟效果的影响

本文利用WRF 模式及其3DVar 同化系统,以2009 年6 月29 日00 时到30 日00 时的湖北鹤峰暴雨为研究个例,对ATOVS 探测器的AMSU-A、AMSU-B 和HIRS 三类资料进行了不同的同化方案试验。结果表明:同化ATOVS 三类资料对暴雨模拟的影响不同,HIRS 资料对暴雨模拟效果改进最大,AMSU-B 次之,AMSU-A 最小。同时同化ATOVS 三类资料时,AMSU-A 资料起主要作用,其同化效果与同时同化ATOVS 三类资料相近,优化组合同化AMSU-B 和HIRS 资料的效果最好。同化ATOVS 不同资料对初始要素场的影响不同,AMSU-A 资料主要影响大气温度场,AMSU-B 资料对中高空要素场的影响较大,HIRS 资料对低空湿度场及风场的协同改变最有利于降水模拟的改善。同时ATOVS 资料的稀疏分辨率也是影响同化效果的一个因子,在模式分辨率不变时,同化资料稀疏分辨率可能存在最佳选择。     

不同卫星微波遥感资料同化对台风路径模拟的影响

本文以2015年13号超强台风“苏迪罗”为个例,利用WRF模式及其3DVar同化系统对NOAA15、NOAA18和NOAA19的AMSU-A微波遥感资料分别同化及组合同化,探究同化不同卫星的同一种微波遥感资料对于台风路径模拟效果的影响。结果表明:同化不同卫星的同一种微波遥感资料对于台风路径模拟具有不同的调整,本文中NOAA15的同化效果最好,其次是NOAA18,最后是NOAA19;同时同化NOAA15、NOAA18和NOAA19的AMSU-A资料并没有取得最好的同化效果,而组合NOAA15和NOAA18则取得了最好的同化效果,即不是同化的卫星数量越多,同化效果越好;同化试验3个时刻的增量场表明同化不同卫星的同一种微波遥感资料对物理场具有不同的调整,这与其对台风模拟路径的调整有着较好的对应,相对于温度场、海平面气压场、位势高度场和风场的增量结构与模拟台风路径的调整更为密切。     

雷达热动力反演资料的高分辨率同化模拟试验

目前,北京地区的天气预报系统对局地对流性定量降水预报能力较弱,远不能满足人们生产、生活和防灾、减灾工作的需要。针对北京地区对提高0-12 h短时临近天气,尤其是夏季局地对流性降水预报能力的需求,基于中国气象局北京城市气象研究所变分多普勒雷达分析系统（VDRAS）的雷达热动力反演资料,建立了WRF模式初始化模块,采用四维资料同化（FDDA）方法,将VDRAS系统高时空分辨率三维热动力结构分析场资料同化到WRF模式中,实现了北京地区VDRAS分析场资料在WRF中尺度模式系统中的应用。通过降水个例的高分辨率同化模拟试验分析了雷达热动力反演资料同化对模式预报结果的影响。研究结果表明:雷达热动力反演资料的同化能够提高模式系统对近地面温、湿、风大气要素和降水过程的模拟能力,改善2 m比湿、降水落区、降水量级、降水时间的预报效果,减少降水漏报的现象。温度和比湿的同化比风的同化对模拟降水结果的改善更重要。虽然研究表明雷达热动力反演资料在WRF模式中的同化能够明显改善模式对选取降水个例的模拟效果,但其对模式尤其是数值业务模式系统预报效果的影响需要进一步更全面、更系统的检验,为业务化应用奠定更坚实的基础。      

A Multivariate Empirical Orthogonal Function-Based Scheme for the Balanced Initial Ensemble Generation of an Ensemble Kalman Filter

The initial ensemble perturbations for an ensemble data assimilation system are expected to reasonably sample model uncertainty at the time of analysis to further reduce analysis uncertainty. Therefore, the careful choice of an initial ensemble perturbation method that dynamically cycles ensemble perturbations is required for the optimal performance of the system. Based on the multivariate empirical orthogonal function (MEOF) method, a new ensemble initialization scheme is developed to generate balanced initial perturbations for the ensemble Kalman filter (EnKF) data assimilation, with a reasonable consideration of the physical relationships between different model variables. The scheme is applied in assimilation experiments with a global spectral atmospheric model and with real observations. The proposed perturbation method is compared to the commonly used method of spatially-correlated random perturbations. The comparisons show that the model uncertainties prior to the first analysis time, which are forecasted from the balanced ensemble initial fields, maintain a much more reasonable spread and a more accurate forecast error covariance than those from the randomly perturbed initial fields. The analysis results are further improved by the balanced ensemble initialization scheme due to more accurate background information. Also, a 20-day continuous assimilation experiment shows that the ensemble spreads for each model variable are still retained in reasonable ranges without considering additional perturbations or inflations during the assimilation cycles, while the ensemble spreads from the randomly perturbed initialization scheme decrease and collapse rapidly.     

单多普勒天气雷达反演风场高原低纬区域同化应用

为了改善低纬高原地区天气预报水平,利用WRF(Weather Research and Forecasting)模式及其变分同化系统进行雷达VAD(Velocity Azimuth Display)反演风场资料同化试验。通过设计不同的试验方案,对2009年6月30日00:00至7月1日00:00发生在云南的一次强降水过程进行数值模拟和对比分析,结果表明:同化VAD反演风场资料后对区域模式的风矢量初始场有明显影响。同化系统能把雷达反演风场信息有效地引入模式初始场,改善强降水区域的水汽输送和风场辐合强度;同化VAD反演风场资料后对区域模式累计降水预报有一定改进作用。从长时间累计降水量定量检验结果看,具体表现为25mm以上量级的降水准确率明显提高、漏报率下降,预报偏差更趋合理。不同的同化试验方案之间的模拟结果差异较大。同化频率越高、同化持续时间越长,对区域模式初始场和预报场的影响越明显。但同化持续时间不宜过长,否则可能导致系统移速过快、降水强度偏大、空报率增加等异常。     

bbGPS/PWV资料三维变分同化改进MM5降水预报连续试验的评估

利用区域地基GPS网反演的高时空密度的大气垂直方向水汽总量，也称为可降水量（PWV），可大大弥补常规探空探测水汽资料的不足。为了全面评估区域GPS网PWV资料同化对业务数值天气预报改进程度的目的，在个例研究分析的基础上，进行了连续38天的GPS/PWV资料三维同化（3D-Var）改进数值业务预报的试验。研究方法是根据长江三角洲地区GPS气象网在2002年梅雨和盛夏季节观测的刖资料，通过三维变分同化建立中尺度数值预报模式MM5的初始场，逐日作出长江三角洲地区24小时的降水量预报。以6小时累积雨量为对象，与未同化GPS／刖资料的MM5的相应预报比较，通过多种评分方法，评估了GPS／PWV资料改进MM5降水预报的效果。结果表明GPS／PWV资料同化后的MM5降水预报能力在大部分时间和大部分地区都有所提高，主要是伪击率有较明显的下降，对小范围降水预报的改进更为明显。预报明显改进的区域恰好位于GPS站填补常规探空站间距较大的地区。     

Assimilation of the multisatellite data into the WRF model for track and intensity simulation of the Indian Ocean tropical cyclones

Weather Research and Forecasting (WRF-ARW) model and its three-dimensional variational data assimilation (3D-Var) system are used to investigate the impact of the Quick Scatterometer (QuikSCAT) near surface winds, Special Sensor Microwave/Imager (SSM/I)-derived Total Precipitable Water (TPW), and Meteosat-7-derived Atmospheric Motion Vectors (AMVs) on the track and intensity prediction of tropical cyclones over the North Indian Ocean. The case of tropical cyclone, Gonu (June 2007; Arabian Sea), is first tested and the results show significant improvements particularly due to the assimilation of QuikSCAT winds. Three other cases, cyclone Mala (April 2006; Bay of Bengal), Orissa super cyclone (October 1999; Bay of Bengal), and Very Severe Cyclonic storm (October 1999; Bay of Bengal), are then examined. The prediction of cyclone tracks improved significantly with the assimilation of QuikSCAT winds. The track improvement resulted from the relocation of the initial cyclonic vortices after the assimilation of QuikSCAT wind vectors. After the assimilation of QuikSCAT winds, the mean (for four cyclone cases) track errors for first, second, and third day forecasts are reduced to 72, 101, and 166?km, respectively, from 190, 250, and 381?km of control (without QuikSCAT winds) runs. The assimilation of QuikSCAT winds also shows positive impact on the intensity (in terms of maximum surface level winds) prediction particularly for those cyclones, which are at their initial stages of the developments at the time of data assimilation. The assimilation of SSM/I TPW has significant influence (negative and positive) on the cyclone track. In three of the four cases, the assimilation of the SSM/I TPW resulted in drying of lower troposphere over cyclonic region. This decrease of moisture in TPW assimilation experiment resulted in reduction of cyclonic intensity. In three of the four cyclones, the assimilation of Meteosat-7 AMVs shows negative impact on the track prediction.     

ANALYSIS OF THE EFFECT OF 3DVAR AND ENSRF DIRECT ASSIMILATION OF RADAR DATA ON THE FORECAST OF A HEAVY RAINFALL EVENT

The present study designs experiments on the direct assimilation of radial velocity and reflectivity data collected by an S-band Doppler weather radar (CINRAD WSR-98D) at the Hefei Station and the reanalysis data produced by the United States National Centers for Environmental Prediction using the Weather Research and Forecasting (WRF) model, the WRF model with a three-dimensional variational (3DVAR) data assimilation system and the WRF model with an ensemble square root filter (EnSRF) data assimilation system. In addition, the present study analyzes a Meiyu front heavy rainfall process that occurred in the Yangtze -Huaihe River Basin from July 4 to July 5, 2003, through numerical simulation. The results show the following. (1) The assimilation of the radar radial velocity data can increase the perturbations in the low-altitude atmosphere over the heavy rainfall region, enhance the convective activities and reduce excessive simulated precipitation. (2) The 3DVAR assimilation method significantly adjusts the horizontal wind field. The assimilation of the reflectivity data improves the microphysical quantities and dynamic fields in the model. In addition, the assimilation of the radial velocity and reflectivity data can better adjust the wind fields and improve the intensity and location of the simulated radar echo bands. (3) The EnSRF assimilation method can assimilate more small-scale wind field information into the model. The assimilation of the reflectivity data alone can relatively accurately forecast the rainfall centers. In addition, the assimilation of the radial velocity and reflectivity data can improve the location of the simulated radar echo bands. (4) The use of the 3DVAR and EnSRF assimilation methods to assimilate the radar radial velocity and reflectivity data can improve the forecast of precipitation, rain-band areal coverage and the center location and intensity of precipitation.     

The effect of three-dimensional variational data assimilation of QuikSCAT data on the numerical simulation of typhoon track and intensity

In this paper, the three-dimensional variational data assimilation scheme (3DVAR) in the mesoscalemod el version 5 (MM5) of the US Pennsylvania State University/National Center for Atmospheric Research is used to study the effect of assimilating the sea-wind data from QuikSCAT on the prediction of typhoon track and intensity. The case of Typhoon Dujuan (2003) is first tested and the results show appreciable improvements. Twelve other cases in 2003 are then evaluated. The assimilation of the QuikSCAT data produces significant impacts on the structure of Dujuan in terms of the horizontal and vertical winds, sealevel pressure and temperature at the initial time. With the assimilation, the 24-h (48-h) track prediction of 11 (10) out of the 12 typhoons is improved. The 24-h (48-h) prediction of typhoon intensity is also improved in 10 (9) of the 12 cases. These experiments therefore demonstrate that assimilation of the QuikSCAT sea-wind data can increase the accuracy of typhoon track and intensity predictions through modification of the initial fields associated with the typhoon.     

复杂地形下地面观测资料同化I. 模式地形与观测站地形高度差异对地面资料同化的影响研究

中国地形复杂，模式地形与实际观测地形存在一定高度差异，因此设计合理的复杂地形下地面观测资料的同化方案有利于使我国目前仅用作探测手段的地面观测资料（常规地面观测站和地面自动站）在中尺度数值模式中得到充分利用。作者在MM5_3DVAR同化系统中利用近地层相似理论将地面观测资料进行直接三维变分同化分析，并对地面资料同化方案设计中是否需要考虑模式与实际观测站地形高度差异进行探讨研究。研究结果表明：通过近地层相似理论将地面观测资料同化到数值模式能起到一定的作用，并且地面观测资料（温度、 湿度、 风场、 地面气压）中各物理量同化到数值模式都能影响24小时降水数值结果，但各物理量起的作用大小不一样，其中影响最大的是温度，其次为湿度；地面观测资料同化方案设计有必要考虑模式地形与实际观测站地形高度差异，适当考虑这种高度差异能取得较好的结果。