Sensitivity of Typhoon Vamei (2001) Simulation to Planetary Boundary Layer Parameterization Using PSU/NCAR MM5

This paper investigates the sensitivity of the numerical simulations of a near equatorial Typhoon Vamei (2001) to various planetary boundary layer (PBL) parameterization schemes in the Pennsylvania State University (PSU)/National Centre for Atmospheric Research (NCAR) non-hydrostatic mesoscale model (MM5). The numerical simulations are conducted on two domains at 45 and 15 km grids nested in a one-way fashion. Four different PBL parameterization schemes including the Blackadar (BLK) scheme, the Burk–Thompson (BURKT) scheme, the NCEP Eta model scheme and the NCEP medium range forecast (MRF) model scheme are investigated. Results indicate that the intensity and propagation track of the simulated near equatorial typhoon system is not very sensitive to the different PBL treatments. The simulated minimum central pressures and the maximum surface wind speeds differ by only 5–6 hPa and 6–8 ms⁻¹, respectively. Larger variations between the simulations occur during the weakening phase of the typhoon system. While all schemes simulated the typhoon with reasonable accuracy, the ETA scheme produces the strongest storm intensity with the largest heat exchanges over the marine environment and the highest warm moisture air content in the PBL around the core of the storm.     

两种边界层参数化方案模拟热带气旋Megi(2010)路径差异的机理分析

本文利用中尺度数值模式WRF,分别采用YSU和MYJ两种边界层参数化方案对2010年超强台风Megi的移动路径进行了模拟,研究了热带气旋(TC)路径模拟对边界层方案的敏感性,并从模拟TC尺度差异所造成的影响角度揭示了模式边界层方案影响TC路径的机理.结果表明:由于两种方案对边界层垂直混合作用过程的描述不同,两个试验模拟的低层水汽垂直输送存在差异.相对于能很好模拟出Megi路径的MYJ方案,YSU方案模拟的TC外围螺旋雨带更活跃,造成TC尺度增大,引起TC中心北侧外围气压梯度和径向风速增加,使得由副高向TC中心输送更多的质量,造成副高异常减弱,从而导致由副高主导的引导气流发生改变,最终使得采用YSU方案模拟的Megi路径出现提前转向.     

Evaluation of Short-Range Forecasts from a Mesoscale Model Over the Indian Region During Monsoon 2006

This study examines the short-range forecast accuracy of the Pennsylvania State University-National Center for Atmospheric Research Mesoscale Model (MM5) as applied to the July 2006 episode of the Indian summer monsoon (ISM) and the model's sensitivity to the choice of different cumulus parameterization schemes (CPSs), namely Betts-Miller, Grell (GR) and Kain-Fritsch (KF). The results showed that MM5 day 1 (0–24 h prediction) and day 2 (24–48 h prediction) forecasts using all three CPSs overpredicted monsoon rainfall over the Indian landmass, with the larger overprediction seen in the day 2 forecasts. Among the CPSs, the rainfall distribution over the Indian landmass was better simulated in forecasts using the KF scheme. The KF scheme showed better skill in predicting the area of rainfall for most of the rainfall thresholds. The root mean square error (RMSE) in day 1 and day 2 rainfall forecasts using different CPSs showed that rainfall simulated using the KF scheme agreed better with the observed rainfall. As compared to other CPSs, simulation using the GR scheme showed larger RMSE in wind speed prediction at 850 and 200 hPa over the Indian landmass. MM5 24-h temperature forecasts at 850 hPa with all the CPSs showed a warm bias of the order of 1 K over the Indian landmass and the bias doubled in 48-h model forecasts. The mean error in temperature prediction at 850 hPa over the Indian region using the KF scheme was comparatively smaller for all the forecast intervals. The model with all the CPSs overpredicted humidity at 850 hPa. The improved prediction by MM5 with the KF scheme is well complemented by the smaller error shown by the KF scheme in vertical distribution of heat and mean moist static energy in the lower troposphere. In this study, the KF scheme which explicitly resolve the downdrafts in the cloud column tended to produce more realistic precipitation forecasts as compared to other schemes which did not explicitly incorporate downdraft effects. This is an important result especially given that the area covered by monsoon-precipitating systems is largely from stratiform-type clouds which are associated with strong downdrafts in the lower levels. This result is useful for improving the treatment of cumulus convection in numerical models over the ISM region.     

人为热源对城市边界层结构影响的数值模拟研究

用南京大学多尺度模式系统在不同区域进行了多种人为热源引入方案的研究,结果表明：考虑时空变化的人为热源按比例分别引入到地表能量平衡方程和大气热量守恒方程是将人为热源引入模式的最优方案.人为热通量密度变化的敏感性试验结果发现：人为热源的存在对城市热岛的生成有重要作用.计算结果表明,南京现有的人为热源排放量对该地区的城市热岛贡献率约为296％,若人为热通量密度在现有量值的基础上增大1倍,则其热岛贡献率可达429％;此外,人为热的排放对清晨城市边界层逆温结构有一定程度的破坏作用,能明显升高夜间近地层气温达05～10℃,并能使白天湍流活动的影响范围增大,混合层高度抬高,使夜间城市热岛环流的影响范围扩大.     

Influence of mantle convection to the crustal movement pattern in the northeastern margin of the Tibetan Plateau based on numerical simulation

Based on the recent observations about the movement and rheological structure of the lithosphere and deformation pattern of the crust, we developed a three-dimensional finite element model for the northeastern margin of the Tibetan Plateau. The model considered the impacts of both external and internal conditions, including mantle convection, gravitational potential energy and block interactions. We compared the simulated surface movement rates to the observed GPS velocities, and the results revealed that crustal movement gradually decreased toward the edge of the plateau. The factors controlling this pattern are the interactions of adjacent blocks, gravitational potential energy of the plateau, and also mantle convection as well. Additionally, according to the observation that there was an apparent difference between the horizontal movement rate of the lithosphere and convective velocity of the underlying mantle, and also based on the results of seismic anisotropy studies that suggest different strengths and deformation regimes of the lithosphere in different tectonic blocks, we proposed that the impact of mantle convection on the lithosphere may have varied in space, and introduced a parameter named mantle convection intensity factor in numerical simulations. Our simulation results show consistent surface movement rates with GPS observations, which further supports the viewpoint of seismic anisotropy studies, i.e., the degree of coupling between the crust and mantle varies significantly among different blocks.     

A Study on the Impact of Observation Assimilation on the Numerical Simulation of Tropical Cyclones JAL and THANE Using 3DVAR

In this work, the impact of assimilation of conventional and satellite remote sensing observations (Oceansat-2 winds, MODIS temperature/humidity profiles) is studied on the simulation of two tropical cyclones in the Bay of Bengal region of the Indian Ocean using a three-dimensional variational data assimilation (3DVAR) technique. The Weather Research and Forecasting (WRF)-Advanced Research WRF (ARW) mesoscale model is used to simulate the severe cyclone JAL: 5–8 November 2010 and the very severe cyclone THANE: 27–30 December 2011 with a double nested domain configuration and with a horizontal resolution of 27 × 9 km. Five numerical experiments are conducted for each cyclone. In the control run (CTL) the National Centers for Environmental Prediction global forecast system analysis and forecasts available at 50 km resolution were used for the initial and boundary conditions. In the second (VARAWS), third (VARSCAT), fourth (VARMODIS) and fifth (VARALL) experiments, the conventional surface observations, Oceansat-2 ocean surface wind vectors, temperature and humidity profiles of MODIS, and all observations were respectively used for assimilation. Results indicate meager impact with surface observations, and relatively higher impact with scatterometer wind data in the case of the JAL cyclone, and with MODIS temperature and humidity profiles in the case of THANE for the simulation of intensity and track parameters. These relative impacts are related to the area coverage of scatterometer winds and MODIS profiles in the respective storms, and are confirmed by the overall better results obtained with assimilation of all observations in both the cases. The improvements in track prediction are mainly contributed by the assimilation of scatterometer wind vector data, which reduced errors in the initial position and size of the cyclone vortices. The errors are reduced by 25, 21, 38 % in vector track position, and by 57, 36, 39 % in intensity, at 24, 48, 72 h predictions, respectively, for the two cases using assimilation of all observations. Simulated rainfall estimates indicate that while the assimilation of scatterometer wind data improves the location of the rainfall, the assimilation of MODIS profiles produces a realistic pattern and amount of rainfall, close to the observational estimates.     

Simulations of Severe Tropical Cyclone Nargis over the Bay of Bengal Using RIMES Operational System

The Regional Integrated Multi-Hazard Early Warning System (RIMES), an international, intergovernmental organization based in Thailand is engaged in disaster risk reduction over the Asia–Pacific region through early warning information. In this paper, RIMES’ customized Weather Research Forecast (WRF) model has been used to evaluate the simulations of cyclone Nargis which hit Myanmar on 2 May 2008, the most deadly severe weather event in the history of Myanmar. The model covers a domain of 35oE to 145oE in the east—west direction and 12oS to 40oN in the north—south direction in order to cover Asia and east Africa with a resolution of 9?km in the horizontal and 28 vertical levels. The initial and boundary conditions for the simulations were provided by the National Center for Environmental Prediction-Global Forecast System (NCEP-GFS) available at 1o lon/lat resolution. An attempt is being made to critically evaluate the simulation of cyclone Nargis by seven set of simulations in terms of track, intensity and landfall time of the cyclone. The seven sets of model simulations were initialized every 12?h starting from 0000 UTC 28 April to 01 May 2008. Tropical Rainfall Measurement Mission (TRMM) precipitation (mm) is used to evaluate the performance of the simulations of heavy rainfall associated with the tropical cyclone. The track and intensity of the simulated cyclone are compared by making use of Joint Typhoon Warning Center (JTWC) data sets. The results indicate that the landfall time, the distribution and intensity of the rainfall, pressure and wind field are well simulated as compared with the JTWC estimates. The average landfall track error for all seven simulations was 64?km with an average time error of about 5?h. The average intensity error of central pressure in all the simulations were found out to be approximately 6?hPa more than the JTWC estimates and in the case of wind, the simulations under predicted it by an average of 12?m?s^?1.     

Estimation of tropical cyclone parameters and wind fields from SAR images

The traditional method of Synthetic Aperture Radar(SAR)wind field retrieval is based on an empirical relation between the near surface winds and the normalized radar backscatter cross section to estimate wind speeds,where this relation is called the geophysical model function(GMF).However,the accuracy rapidly decreases due to the impact of rainfall on the measurement of SAR and the saturation of backscattered intensity under the condition of tropical cyclone.Because of no available instrument synchronously monitoring rain rate on the satellite platform of SAR,we have to derive the precipitation of the SAR observation time from non-simultaneous passive microwave observations of rain in combination with geostationary IR images,and then use the model of rain correction to remove the impact of rain on SAR wind field measurements.For the saturation of radar backscatter cross section in high wind speed conditions,we develop an approach to estimate tropical cyclone parameters and wind fields based on the improved Holland model and the SAR image features of tropical cyclone.To retrieve the low-to-moderate wind speed,the wind direction of tropical cyclone is estimated from the SAR image using wavelet analysis.And then the maximum wind speed and the central pressure of tropical cyclone are calculated by a least square minimization of the difference between the improved Holland model and the low-to-moderate wind speed retrieved from SAR.In addition,wind fields are estimated from the improved Holland model using the above-mentioned parameters of tropical cyclone as input.To evaluate the accuracy of our approach,the SAR images of typhoon Aere,typhoon Khanun,and hurricane Ophelia are used to estimate tropical cyclone parameters and wind fields,which are compared with the best track data and reanalyzed wind fields of the Joint Typhoon Warning Center(JTWC)and the Hurricane Research Division(HRD).The results indicate that the tropical cyclone center,maximum wind speed,and central pressure are generally consistent with the best track data,and wind fields agree well with reanalyzed data from HRD.     

双模态风垂直切变对台风结构和强度的作用分析

本文应用WRF(v 3.4)模式输出资料,揭示了风垂直切变(Vertical Wind Shear:VWS)在垂直方向上的波状变化特征,这种波状变化在台风不同发展时期又有不同形态,其中在持续强盛期呈双模态分布.应用VWS引起的次级垂直环流影响台风对流分布和强度变化的基本原理,用模式资料分析发现:对流层中层具有的VWS是整层VWS的主要部分,台风强度变化滞后VWS的形态突变6h左右;双模态波状变化的VWS产生的次级环流和台风垂直环流的配置不同使台风强对流带结构变得不对称及眼墙区对流强度在垂直分布上变得不均匀,随着持续强盛期涡旋运动的增强,强对流带分布又趋于对称.又根据VWS形成的垂直方向上涡度力分布不均匀引起台风内中尺度滚轴状对流带不稳定发展原理,分析表明:对流层中、低层的涡度力有利于对流不稳定增强,垂直速度的最大值与风垂直廓线的拐点在同一高度上,这与理论模型的结论一致.因此,VWS的波状变化分布特征不仅影响台风强对流带中尺度结构的改变,也对台风持续强盛具有重要作用;同时也是台风内滚轴状对流带不稳定的可能启动机制.     

Air–Sea Interaction Processes Influencing the Development of A Shapiro–Keyser Type Cyclone over the Subtropical South Atlantic Ocean

This study analyzes the impacts of latent and sensible heat exchanges between the atmosphere and the ocean in a non-explosive Shapiro–Keyser type cyclogenesis event that occurred over the southwestern South Atlantic Ocean. The synoptic evolution shows a relatively strong warm front and a cold frontal fracture during the system’s development and a warm seclusion in its mature stage, characterizing a Shapiro–Keyser type cyclone. Numerical experiments with the ARW-WRF Model version 3.3 were used to investigate the influences of sensible and latent fluxes on the track of the surface low, intensity of the fronts and coupling of the lower and upper troposphere. The simulations indicate that in the presence of these fluxes the cyclone underwent greater intensification, had a longer life time and longer trajectory, and presented a typical southeastward movement. In the absence of these fluxes, the cyclone developed a weaker warm front with consequent reduction of diabatic heating due to grid scale precipitation along it. This reduced the negative pressure tendency southeast of the cyclone center and the surface cyclone moved northeastward, showing a decoupling of the lower- and upper-level waves. A consequence of this anomalous tracking is the location of the surface cyclone beneath the upper-level trough axis, where there is no upper-level divergence associated with cyclonic vorticity advection contributing to the further system intensification. Numerical experiments suggest that for this Shapiro–Keyser type cyclone the air–sea interaction processes are crucial to obtain a cyclone with features similar to the observations.     

Assimilation of Doppler Weather Radar Data in WRF Model for Simulation of Tropical Cyclone Aila

For the accurate and effective forecasting of a cyclone, it is critical to have accurate initial structure of the cyclone in numerical models. In this study, Kolkata Doppler weather radar (DWR) data were assimilated for the numerical simulation of a land-falling Tropical Cyclone Aila (2009) in the Bay of Bengal. To study the impact of radar data on very short-range forecasting of a cyclone's path, intensity and precipitation, both reflectivity and radial velocity were assimilated into the weather research and forecasting (WRF) model through the ARPS data assimilation system (ADAS) and cloud analysis procedure. Numerical experiment results indicated that radar data assimilation significantly improved the simulated structure of Cyclone Aila. Strong influences on hydrometeor structures of the initial vortex and precipitation pattern were observed when radar reflectivity data was assimilated, but a relatively small impact was observed on the wind fields at all height levels. The assimilation of radar wind data significantly improved the prediction of divergence/convergence conditions over the cyclone's inner-core area, as well as its wind field in the low-to-middle troposphere (600–900 hPa), but relatively less impact was observed on analyzed moisture field. Maximum surface wind speed produced from DWR–Vr and DWR–ZVr data assimilation experiments were very close to real-time values. The impact of radar data, after final analysis, on minimum sea level pressure was relatively less because the ADAS system does not adjust for pressure due to the lack of pressure observations, and from not using a 3DVAR balance condition that includes pressure. The greatest impact of radar data on forecasting was realized when both reflectivity and wind data (DWR–ZVr and DWR–ZVr00 experiment) were assimilated. It is concluded that after final analysis, the center of the cyclone was relocated very close to the observed position, and simulated cyclone maintained its intensity for a longer duration. Using this analysis, different stages of the cyclone are better captured, and cyclone structure, intensification, direction of movement, speed and location are significantly improved when both radar reflectivity and wind data are assimilated. As compared to other experiments, the maximum reduction in track error was noticed in the DWR–ZVr and DWR–ZVr00 experiments, and the predicted track in these experiments was very close to the observed track. In the DWR–ZVr and DWR–ZVr00 experiments, rainfall pattern and amount of rainfall forecasts were remarkably improved and were similar to the observation over West Bengal, Orissa and Jharkhand; however, the rainfall over Meghalaya and Bangladesh was missed in all the experiments. The influence of radar data reduces beyond a 12-h forecast, due to the dominance of the flow from large-scale, global forecast system models. This study also demonstrates successful coupling of the data assimilation package ADAS with the WRF model for Indian DWR data.     

Numerical simulation on the evolution of cloud particles in 3-D convective cloud

A 3-D convective cloud model with compressible non-hydrostatic dynamics and the spectral bin microphysics of a 2-D slab-symmetric model has been used to simulate an observed supercell storm occurring on 29 June, 2000 near Bird City, Kansas, USA. The main objective of this paper is to study the evolution of particles in this convective storm with bin spectral microphysics scheme. Graupels form and grow through two mechanisms, deposition and riming, with the riming process dominant on top of the inflow and in the upper portion of main updraft. Over the outflow and during the developing and mature stages of the storm, graupel particles mainly grow through deposition with dominant unimodal spectra. Most fall out after growing up. Reducing initial relative humidity disturbance (increasing initial potential temperature disturbance) has negative impact on the formation and growth of graupels over the inflow (outflow). This study shows that large graupel and hail could be suppressed by altering the deposition and coalescence process over the inflow and main updraft. At different locations of the convective cells and with different initial humidity and potential temperature disturbance, the graupel formation and growth mechanisms are different, so as to the feasible hail suppression locations and methods. Supported by National Natural Science Foundation of China (Grant Nos. 40537034, 40805057), and Foundation of Key Laboratory of Meteorological Disaster of Ministry of Education in Nanjing University of Information Science & Technology (Grant No. KLME060202)     

中尺度模式积云参数化方案的改进及其在暴雨模拟中的应用

本文针对中国暴雨发生发展天气特征,改进和发展了一种适合于描述东亚暴雨的中尺度积云参数化方案.首先基于近年来(1990—2010)江淮流域汛期降水合成分析的基础上,诊断出组织化对流降水环境的动力参数;其次利用该动力参数作为动力控制条件,改进了Kain-Fritsch Eta中尺度积云对流参数化方案;最后利用改进的中尺度积云参数化方案对梅雨期暴雨、华南前汛期暴雨过程进行了数值模拟,结果表明:改进后的中尺度积云参数化方案对上述两次暴雨过程的落区及强度的模拟,均有明显改进.     

Improving NWP through radar rainfall-driven land surface parameters: A case study on convective precipitation forecasting

In this study we investigate the effect of forcing the land surface scheme of an atmospheric mesoscale model with radar rainfall data instead of the model-generated rainfall fields. The goal is to provide improved surface conditions for the atmospheric model in order to achieve accurate simulations of the mesoscale circulations that can significantly affect the timing, distribution and intensity of convective precipitation. The performance of the approach is evaluated in a set of numerical experiments on the basis of a 2-day-long mesoscale convective system that occurred over the US Great Plains in July 2004. The experimental design includes multiple runs covering a variety of forcing periods. Continuous data integration was initially used to investigate the sensitivity of the model’s performance in varying soil state conditions, while shorter time windows prior to the storm event were utilized to assess the effectiveness of the procedure for improving convective precipitation forecasting. Results indicate that continuous integration of radar rainfall data brings the simulated precipitation fields closer to the observed ones, as compared to the control simulation. The precipitation forecasts (up to 48 h) appear improved also in the cases of shorter integration periods (24 and 36 h), making this technique potentially useful for operational settings of weather forecasting systems. A physical interpretation of the results is provided on the basis of surface moisture and energy exchange.     

Internal airflow of a convective storm from dual-Doppler radar measurements

The dual-Doppler radar coplane method of scanning and data reduction has been used to determine the internal airflow and radar reflectivity structure of a convective storm. Cumulus convection growing in a moderately sheared wind environment resulted in a nonsteady, moderate intensity thunderstorm. Precipitation fallout and downward moving air are found downshear of an updraft inclined in the downshear direction. Rapid storm translation, vertical shear of the ambient wind, and slow subcloud ascent of inflow air act to establish this observed draft configuration. The absence of significant cold air outflow and its attendant gust front at the surface is attributed to (i) appreciable inflow of slow moving air into the downdraft at the middle layers and (ii) the fact that the potentially coldest air was located too low to contribute significantly to a deep downdraft circulation.     

Simulation of the time-space evolution of an extratropical cyclonic precipitation field over USA

The time-space evolution of an extratropical cyclonic precipitation field over U S A is simulated in a stochastic setting as outlined in Kavvas et al. (1988). The birth of a cyclonic storm is characterized by the simultaneous birth of a cyclone center and births of subsynoptic precipitation areas (SPA) at preferred locations around the cyclone center. The precipitation cores and cells which are used as the fundamental building blocks of the SPAs are approximated by circular precipitation areas (CPA) of different sizes. The time space evolution of the precipitation field after the birth is governed by (1) the movement of the synoptic cyclone described by the cyclone center trajectory, (2) independent nonidentically distributed random velocities of the individual CPAs relative to the cyclone center, (3) the births of new CPAs in time and space relative to the cyclone center, (4) the independent evolution in time of the individual spatially uniform intensities of the existing CPAs, (5) the expansion and shrinkage of the existing CPAs in the course of movement and (6) the dissipation (death) of a random number of existing CPAs within the cyclonic system. The computer simulation, the results of which are presented in this paper, successfully reproduced the general mesoscale and synoptic scale features of the radar detected cyclonic rain fields as observed by Austin and Houze (1972), Houze et al. (1976), Hobbs (1978), Hobbs and Locatelli (1978), Houze (1981), Houze and Hobbs (1982) and others.     

基于CloudSat资料的北上江淮气旋暴雪云系结构特征

2007年3月3-5日和2013年11月24-25日,受江淮气旋北上影响,我国北方大部地区遭遇罕见暴风雪天气,2次暴雪过程有很多相似之处.利用常规观测、CloudSat卫星云廓线雷达的探测资料和NECP/NCAR再分析资料,分析了这2次暴雪过程江淮气旋云系结构和微物理特征.结果表明：（1）北上江淮气旋的冷锋云系较窄,以深厚对流云为主,回波核心在2~7 km,其结构在气旋发展的不同阶段变化不大;（2）逗点头云系范围宽广,在气旋的不同发展阶段,结构和强度有显著差异.气旋初始锋面波动和锋面断裂阶段,逗点头云系有两个降水区：北部为由多个单体组成的大范围层状云区,强回波从地表向上伸展,上空有高空对流泡,建立了播撒云-供水云机制,有利于下部冰晶粒子长大;南部有对流云柱发展.逗点头西部的冷输送带云系主要集中在6 km以下,强度弱,冰粒子含量少;（3）气旋暖锋后弯阶段,干侵入加强,冷锋后部的无云区或少云区范围扩大,逗点头云系南北范围收缩、变窄,云系的高度、强度和含水量减弱,冷锋云系也减弱;（4）气旋冷锋云系和逗点头南部的对流云柱以降雨为主,位于高纬度陆地上的逗点头云系以降雪为主,当逗点头云系处于海上有对流不稳定发展,以降雨为主.冷锋云系北部和逗点头云系南部均有由层积云或高积云组成的低云,以毛毛雨为主.冷锋云系和逗点头云系北部100-200 km的范围为随高度和距离逐渐变薄的高层云,无降水对应.     

Development and evaluation of a storm surge warning system in Taiwan

An operational storm surge forecasting system aimed at providing warning information for storm surges has been developed and evaluated using four typhoon events. The warning system triggered by typhoon forecasts from Taiwan Cooperative Precipitation Ensemble Forecast Experiment (TAPEX) has been executed with two storm surge forecasting scenarios with and without tides. Three numerical experiments applying different meteorological inputs have been designed to assess the impact of typhoon forcing on storm surges. One uses synthetic wind fields, and the others use realistic wind fields with and without adjustments to the initial wind fields for the background circulation. Local observations from Central Weather Bureau (CWB) weather stations and tide gauge stations are used to evaluate the wind fields and storm surges from our numerical experiments. The comparison results show that the accuracy of the storm surge forecast is dominated by the track, the intensity, and the driving flow of a typhoon. When the structure of a typhoon is disturbed by Taiwan’s topography, using meteorological inputs from real wind fields can result in a better typhoon simulation than using inputs from synthetic wind fields. The driving flow also determines the impact of topography on typhoon movement. For quickly moving typhoons, storm forcing from TAPEX is reliable when a typhoon is strong enough to be relatively unaffected by environmental flows; otherwise, storm forcing from a sophisticated typhoon initialization scheme that better simulates the typhoon and environmental flows results in a more accurate prediction of storm surges. Therefore, when a typhoon moves slowly and interacts more with the topography and environmental flows, incorporating realistic wind fields with adjustments to the initial wind fields for the background circulation in the warning system will obtain better predictions for a typhoon and its resultant storm surges.     

On the Prediction of Tropical Cyclones over the Indian Region Using a Synthetic Vortex Scheme in a Mesoscale Model

The tropical cyclones form over the oceanic regions where conventional meteorological observations are not available. This contributes to a poor initial analysis of the cyclonic vortex and hence inadequate forecast. One way of overcoming the above problem is to modify the initial analysis by replacing the weak and ill-defined vortex in the initial analysis with a synthetic vortex having the correct size and intensity at the correct location. In this study we are investigating the effect of inclusion of a synthetic vortex based on Rankine as well as on Holland wind profiles, using NCAR-AFWA bogussing scheme for the prediction of four tropical cyclones, which formed over the Bay of Bengal during November 2002 and 2005, December 2005 and over the Arabian Sea during May 2004, using the MM5 model. Two numerical experiments are designed in this study for each of the above four cyclones. In the first experiment the model is integrated with a synthetic vortex based on Rankine wind profile while in the second experiment we utilize the Holland wind profile. For the November 2002 cyclone, in both the experiments the model is integrated from 10 November 2002 18 UTC to 12 November, 2002 12 UTC with the synthetic vortex inserted at the initial time. The results of the study for the November 2002 cyclone show that the model simulation with the Holland vortex has produced a stronger cyclone in terms of minimum sea-level pressure and maximum wind speed. Also, the results for the November 2002 cyclone with the Holland vortex showed a better longitudinal height section of the horizontal wind speed across the center of the cyclone. The track error of the cyclone for the November 2002 cyclone is less in the model simulation with the Holland vortex at the initial time and at 24 hours of forecast. The results for the November 2002 cyclone with the Rankine vortex showed greater vertical wind speed as compared to the Holland vortex. However, for the November 2002 cyclone there were no significant differences in the spatial distribution of precipitation for both the experiments. In order to provide an adequate number of case studies for a good statistical sample, the present study is extended for three additional cyclones over the Indian region. All four cyclones studied here show that the Holland vortex has produced a stronger cyclone in terms of the minimum sea-level pressure and maximum wind speed. The Holland vortex showed a better vertical structure of wind speed in the longitudinal height section at 24 hours of forecast for the November 2005 cyclone while the structure was better for the Rankine vortex for the remaining two cyclones. There were no significant differences in the spatial distribution of precipitation for the two experiments corresponding to all four cyclones. Some statistical results pertaining to all four cyclones are provided such as the average track error as well as the average difference between the observed and the model minimum sea-level pressure and the maximum wind speed. The statistical results corresponding to the average of all the four cyclones are at only a slight variance with the results corresponding to the November 2002 cyclone.     

Impacts of boundary layer parameterization schemes and air-sea coupling on WRF simulation of the East Asian summer monsoon

The planetary boundary layer(PBL)scheme in the regional climate model(RCM)has a significant impact on the interactions and exchanges of moisture,momentum,and energy between land,ocean,and atmosphere;however,its uncertainty will cause large systematic biases of RCM.Based on the four different PBL schemes(YSU,ACM2,Boulac,and MYJ)in Weather Research and Forecasting(WRF)model,the impacts of these schemes on the simulation of circulation and precipitation during the East Asian summer monsoon(EASM)are investigated.The simulated results of the two local turbulent kinetic energy(TKE)schemes,Boulac and MYJ,are more consistent with the observations than those in the two nonlocal closure schemes,YSU and ACM2.The former simulate more reasonable low-level southwesterly flow over East China and west pacific subtropical high(WPSH)than the latter.As to the modeling of summer monsoon precipitation,both the spatial distributions and temporal evolutions from Boulac and MYJ are also better than those in YSU and ACM2 schemes.In addition,through the comparison between YSU and Boulac experiments,the differences from the results of EASM simulation are more obvious over the oceanic area.In the experiments with the nonlocal schemes YSU and ACM2,the boundary layer mixing processes are much stronger,which lead to produce more sea surface latent heat flux and enhanced convection,and finally induce the overestimated precipitation and corresponding deviation of monsoon circulation.With the further study,it is found that the absence of air-sea interaction in WRF may amplify the biases caused by PBL scheme over the ocean.Consequently,there is a reduced latent heat flux over the sea surface and even more reasonable EASM simulation,if an ocean model coupled into WRF.