Numerical simulation of tropical cyclone thane: role of boundary layer and surface drag parameterization schemes

In the recent times, several advanced numerical models are utilized for the prediction of the intensity, track and landfall time of a cyclone. Still there are number of issues concerning their prediction and the limitation of numerical models in addressing those issues. The most pertinent question is how intensive a cyclone can become before it makes a landfall and where the cyclone moves under the ambient large-scale flow. In this paper, detailed study has been carried out using Weather Research Forecast model with two boundary schemes to address the above question by considering a recent tropical cyclone in Bay of Bengal region of North Indian Ocean. In addition, the impact of the surface drag effect on the low-level winds and the intensity of the cyclone are also studied. The result reveals that large differences are noted in the ocean surface fluxes between YSU and MYJ with MYJ producing relatively higher fluxes than YSU. It is found that the YSU scheme produced a better simulation for the THANE cyclone in terms of winds, pressure distribution and cloud fractions. Comparison with available observations indicated the characteristics of horizontal divergence, vorticity and vector track positions produced by YSU experiment are more realistic than with MYJ and other experiments. However, when the drag coefficient is changed as 0.5 or 2.0 from the default values, appreciable changes in the surface fluxes are not noticed. A maximum precipitation is reported in YSU as compared to the MYJ PBL scheme for the tropical cyclone THANE.     

Sensitivity of simulated cyclone Gonu intensity and track to variety of parameterizations: Advanced hurricane WRF model application

Domain configuration and several physical parameterization settings such as planetary boundary layer, cumulus convection, and ocean–atmosphere surface flux parameterizations can play significant roles in numerical prediction of tropical cyclones. The present study focuses to improve the prediction of the TC Gonu by investigating the sensitivity of simulations to mentioned configurations with the Advanced Hurricane WRF model. The experiments for domain design sensitivity with 27 km resolution has been shown moving the domains towards the east improve the results, due to better account for the large-scale process. The fixed and movable nests on a 9-km grid were considered separately within the coarse domain and their results showed that despite salient improvement in simulated intensity, an accuracy reduction in simulated track was observed. Increasing horizontal resolution to 3 km incredibly reduced the simulated intensity accuracy when compared to 27 km resolution. Thereafter, different initial conditions were experimented and the results have shown that the cyclone of 1000 hPa sea level pressure is the best simulation initial condition in predicting the track and intensity for cyclone Gonu. The sensitivity of simulations to ocean–atmosphere surface-flux parameterizations on a 9-km grid showed the combination of ‘Donelan scheme’ for momentum exchanges along with ‘Large and Pond scheme’ for heat and moisture exchanges provide the best prediction for cyclone Gonu intensity. The combination of YSU and MYJ PBL scheme with KF convection for prediction of track and the combination of YSU PBL scheme with KF convection for prediction of intensity are found to have better performance than the other combinations. These 22 sensitivity experiments also implicitly lead us to the conclusion that each particular forecast aspect of TC (e.g., track, intensity, etc.) will require its own special design.     

Simulation of atmospheric states for a storm surge on the west coast of Korea: model comparison between MM5, WRF and COAMPS

High-quality informations on sea level pressure and sea surface wind stress are required to accurately predict storm surges over the Korean Peninsula. The storm surge on 31 March 2007 at Yeonggwang, on the western coast, was an abrupt response to mesocyclone development. In the present study, we attempted to obtain reliable surface winds and sea level pressures. Using an optimal physical parameterization for wind conditions, MM5, WRF and COAMPS were used to simulate the atmospheric states that accompanied the storm surge. The use of MM5, WRF and COAMPS simulations indicated the development of high winds in the strong pressure gradient due to an anticyclone and a mesocyclone in the southern part of the western coast. The response to this situation to the storm surge was sensitive. A low-level warm advection was examined as a possible causal mechanism for the development of a mesocyclone in the generating storm surge. The low-level warm temperature advection was simulated using the three models, but MM5 and WRF tended to underestimate the warm tongue and overestimate the wind speed. The WRF simulation was closer to the observed data than the other simulations in terms of wind speed and the intensity of the mesocyclone. It can be concluded that the magnitude of the storm surge at Yeonggwang was dependent, not only on the development of a mesocyclone but on ocean effects as well.     

Effects of various combinations of boundary layer schemes and microphysics schemes on the track forecasts of tropical cyclones over the South China Sea

This study investigates the effects of various combinations of the planetary boundary layer (PBL) schemes and the microphysics schemes on the numerical forecasting of tropical cyclones (TCs). Using different combinations of three PBL schemes (YSU, MYJ and MYNN2) and four microphysics schemes (Ferrier, Goddard, WSM6 and Lin), a number of experiments are carried out for five landed TCs in the South China Sea during 2012. Results show that the combination of the YSU and Ferrier schemes performs the best for the TC track forecasting, although it does not perform the best for the forecast of precipitation. Further analysis reveals that the best performance of the track forecast by the combination of the YSU and Ferrier schemes mainly attributes to a more accurate steering flow as well as TC wind structure produced by this combination. These results provide a valuable reference to the operational numerical forecasting of TC tracks in the future.     

Customization of WRF-ARW model with physical parameterization schemes for the simulation of tropical cyclones over North Indian Ocean

The convection and planetary boundary layer (PBL) processes play significant role in the genesis and intensification of tropical cyclones (TCs). Several convection and PBL parameterization schemes incorporate these processes in the numerical weather prediction models. Therefore, a systematic intercomparison of performance of parameterization schemes is essential to customize a model. In this context, six combinations of physical parameterization schemes (2 PBL Schemes, YSU and MYJ, and 3 convection schemes, KF, BM, and GD) of WRF-ARW model are employed to obtain the optimum combination for the prediction of TCs over North Indian Ocean. Five cyclones are studied for sensitivity experiments and the out-coming combination is tested on real-time prediction of TCs during 2008. The tracks are also compared with those provided by the operational centers like NCEP, ECMWF, UKMO, NCMRWF, and IMD. It is found that the combination of YSU PBL scheme with KF convection scheme (YKF) provides a better prediction of intensity, track, and rainfall consistently. The average RMSE of intensity (13?hPa in CSLP and 11?m?s^?1 in 10-m wind), mean track, and landfall errors is found to be least with YKF combination. The equitable threat score (ETS) of YKF combination is more than 0.2 for the prediction of 24-h accumulated rainfall up to 125?mm. The vertical structural characteristics of cyclone inner core also recommend the YKF combination for Indian seas cyclones. In the real-time prediction of 2008 TCs, the 72-, 48-, and 24-h mean track errors are 172, 129, and 155?km and the mean landfall errors are 125, 73, and 66?km, respectively. Compared with the track of leading operational agencies, the WRF model is competing in 24?h (116?km error) and 72?h (166?km) but superior in 48-h (119?km) track forecast.     

A comprehensive study of surface and upper-air characteristics over two stations on the west coast of India during the occurrence of a cyclonic storm

While qualitative information from meteorological satellites has long been recognized as critical for monitoring weather events such as tropical cyclone activity, quantitative data are required to improve the numerical prediction of these events. In this paper, the sea surface winds from QuikSCAT, cloud motion vectors and water vapor winds from KALPANA-1 are assimilated using three-dimensional variational assimilation technique within Weather Research Forecasting (WRF) modeling system. Further, the sensitivity experiments are also carried out using the available cumulus convective parameterizations in WRF modeling system. The model performance is evaluated using available observations, and both qualitative and quantitative analyses are carried out while analyzing the surface and upper-air characteristics over Mumbai (previously Bombay) and Goa during the occurrence of the tropical cyclone PHYAN at the west coast of Indian subcontinent. The model-predicted surface and upper-air characteristics show improvements in most of the situations with the use of the satellite-derived winds from QuikSCAT and KALPANA-1. Some of the model results are also found to be better in sensitivity experiments using cumulus convection schemes as compared to the CONTROL simulation.     

Skills of different mesoscale models over Indian region during monsoon season: Forecast errors

Performance of four mesoscale models namely, the MM5, ETA, RSM and WRF, run at NCMRWF for short range weather forecasting has been examined during monsoon-2006. Evaluation is carried out based upon comparisons between observations and day-1 and day-3 forecasts of wind, temperature, specific humidity, geopotential height, rainfall, systematic errors, root mean square errors and specific events like the monsoon depressions.It is very difficult to address the question of which model performs best over the Indian region? An honest answer is ‘none’. Perhaps an ensemble approach would be the best. However, if we must make a final verdict, it can be stated that in general, (i) the WRF is able to produce best All India rainfall prediction compared to observations in the day-1 forecast and, the MM5 is able to produce best All India rainfall forecasts in day-3, but ETA and RSM are able to depict the best distribution of rainfall maxima along the west coast of India, (ii) the MM5 is able to produce least RMSE of wind and geopotential fields at most of the time, and (iii) the RSM is able to produce least errors in the day-1 forecasts of the tracks, while the ETA model produces least errors in the day-3 forecasts.     

On the role of the Planetary Boundary Layer in the numerical simulation of a severe cyclonic storm Nargis using a mesoscale model

The very severe cyclonic storm Nargis of 2008 was a strong tropical cyclone that caused the deadliest natural disaster in the history of Myanmar. The time tested NCAR/PSU MM5 model has been used to simulate the Nargis cyclone, which is designed to have two domains covering the Bay of Bengal with horizontal resolutions of 90 and 30?km. The physics options chosen are Kain?CFritsch 2 for convection, Blackadar (BLA), Burk?CThompson, medium range forecast (MRF), Eta Mellor?CYamada (Eta MY) and Gayno?CSeaman (GS) for Planetary Boundary Layer (PBL) and Simple Ice for explicit cloud physics processes. The experiment was conducted with the model integration starting from April 27, 2008, to May 3, 2008. The performance of the five PBL schemes is evaluated in terms of radius height cross-section of the three component winds, surface heat fluxes of sensible heat and latent heat, equivalent potential temperature (?? _e ), precipitation, track and variation of Central Surface Pressure and wind speed with time. The numerical results show a large impact of the PBL schemes on the intensity and movement of the system. The intensity of the storm is examined in terms of pressure drop, strength of the surface wind and rainfall associated with the storm. The results are compared to the India Meteorological Department observations. These experiments indicate that the intensity of the storm is well simulated with the Eta MY and BLA with finer resolution. The simulated track with MRF compared well with the Joint Typhoon Warning Center observation at landfall position both with the 90 and 30?km resolutions.     

The effect of a surface data assimilation technique and the traditional four-dimensional data assimilation on the simulation of a monsoon depression over India using a mesoscale model

The objective of this study is to investigate the impact of a surface data assimilation (SDA) technique, together with the traditional four-dimensional data assimilation (FDDA), on the simulation of a monsoon depression that formed over India during the field phase of the 1999 Bay of Bengal Monsoon Experiment (BOBMEX). The SDA uses the analyzed surface data to continuously assimilate the surface layer temperature as well as the water vapor mixing ratio in the mesoscale model. The depression for the greater part of this study was offshore and since successful application of the SDA would require surface information, a method of estimating surface temperature and surface humidity using NOAA-TOVS satellites was used. Three sets of numerical experiments were performed using a coupled mesoscale model. The first set, called CONTROL, uses the NCEP (National Center for Environmental Prediction) reanalysis for the initial and lateral boundary conditions in the MM5 simulation. The second and the third sets implemented the SDA of temperature and moisture together with the traditional FDDA scheme available in the MM5 model. The second set of MM5 simulation implemented the SDA scheme only over the land areas, and the third set extended the SDA technique over land as well as sea. Both the second and third sets of the MM5 simulation used the NOAA-TOVS and QuikSCAT satellite and conventional upper air and surface meteorological data to provide an improved analysis. The results of the three sets of MM5 simulations are compared with one another and with the analysis and the BOBMEX 1999 buoy, ship, and radiosonde observations. The predicted sea level pressure of both the model runs with assimilation resembles the analysis closely and also captures the large-scale structure of the monsoon depression well. The central sea level pressures of the depression for both the model runs with assimilation were 2–4 hPa lower than the CONTROL. The results of both the model runs with assimilation indicate a larger spatial area as well as increased rainfall amounts over the coastal regions after landfall compared with the CONTROL. The impact of FDDA and SDA, the latter over land, resulted in reduced errors of the following: 1.45 K in temperature, 0.39 m s⁻¹ in wind speed, and 14° in wind direction compared with the BOBMEX buoy observation, and 1.43 m s⁻¹ in wind speed, 43° in wind direction, and 0.75% in relative humidity compared with the CONTROL. The impact of SDA over land and sea compared with SDA over land only showed a further marginal reduction of errors: 0.23 K in air temperature (BOBMEX buoy) and 1.33 m s⁻¹ in wind speed simulations.     

Sensitivity of tropical cyclone intensification to boundary layer and convective processes

This study examines the role of the parameterization of convection, planetary boundary layer (PBL) and explicit moisture processes on tropical cyclone intensification. A high-resolution mesoscale model, National Center for Atmospheric Research (NCAR) model MM5, with two interactive nested domains at resolutions 90 km and 30 km was used to simulate the Orissa Super cyclone, the most intense Indian cyclone of the past century. The initial fields and time-varying boundary variables and sea surface temperatures were taken from the National Centers for Environmental Prediction (NCEP) (FNL) one-degree data set. Three categories of sensitivity experiments were conducted to examine the various schemes of PBL, convection and explicit moisture processes. The results show that the PBL processes play crucial roles in determining the intensity of the cyclone and that the scheme of Mellor-Yamada (MY) produces the strongest cyclone. The combination of the parameterization schemes of MY for planetary boundary layer, Kain-Fritsch2 for convection and Mixed-Phase for explicit moisture produced the best simulation in terms of intensity and track. The simulated cyclone produced a minimum sea level pressure of 930 hPa and a maximum wind of 65 m s⁻¹ as well as all of the characteristics of a mature tropical cyclone with an eye and eye-wall along with a warm core structure. The model-simulated precipitation intensity and distribution were in good agreement with the observations. The ensemble mean of all 12 experiments produced reasonable intensity and the best track.     

A sensitivity study of the WRF model in offshore wind modeling over the Baltic Sea

Accurate wind modeling is important for wind resources assessment and wind power forecasting.To improve the WRF model configuration for the offshore wind modeling over the Baltic Sea,this study per-formed a sensitivity study of the WRF model to multiple model configurations,including domain setup,grid resolution,sea surface temperature,land surface data,and atmosphere-wave coupling.The simu-lated offshore wind was evaluated against LiDAR observations under different wind directions,atmo-spheric stabilities,and sea status.Generally,the simulated wind profiles matched observations,despite systematic underestimations.Strengthening the forcing from the reanalysis data through reducing the number of nested domains played the largest role in improving wind modeling.Atmosphere-wave cou-pling further improved the simulated wind,especially under the growing and mature sea conditions.Increasing the vertical resolution,and updating the sea surface temperature and the land surface infor-mation only had a slight impact,mainly visible during very stable conditions.Increasing the horizontal resolution also only had a slight impact,most visible during unstable conditions.Our study can help to improve the wind resources assessment and wind power forecasting over the Baltic Sea.     

A Study on the Impact of Parameterization of Physical Processes on Prediction of Tropical Cyclones over the Bay of Bengal With NCAR/PSU Mesoscale Model

Prediction of the track and intensity of tropical cyclones is one of the most challenging problems in numerical weather prediction (NWP). The chief objective of this study is to investigate the performance of different cumulus convection and planetary boundary layer (PBL) parameterization schemes in the simulation of tropical cyclones over the Bay of Bengal. For this purpose, two severe cyclonic storms are simulated with two PBL and four convection schemes using non-hydrostatic version of MM5 modeling system. Several important model simulated fields including sea level pressure, horizontal wind and precipitation are compared with the corresponding verification analysis/observation. The track of the cyclones in the simulation and analysis are compared with the best-fit track provided by India Meteorological Department (IMD). The Hong-Pan PBL scheme (as implemented in NCAR Medium Range Forecast (MRF) model) in combination with Grell (or Betts-Miller) cumulus convection scheme is found to perform better than the other combinations of schemes used in this study. Though it is expected that radiative processes may not have pronounced effect in short-range forecasts, an attempt is made to calibrate the model with respect to the two radiation parameterization schemes used in the study. And the results indicate that radiation parameterization has noticeable impact on the simulation of tropical cyclones.     

Surface wind speed probability distribution in the Southeast Pacific of Marine Stratus and Stratocumulus regions

Probability distributions of surface wind speeds (SWS) near coastal regions are needed for applications such as estimating offshore wind power and ocean surface fluxes and for offshore wind risk assessments. Ocean surface wind speed probability distribution (PDF) is characterized using three-year QuikSCAT and AIRS satellite observations in the southeast Pacific of marine stratus and stratocumulus (MSC) regions. Seasonal variation is removed from wind statistics. It was found that the observed SWS standard deviation has a linear positive relationship with its mean SWS; while the SWS skewness decreases with mean SWS in regimes of strong winds and increases with mean SWS in regimes of weak winds. A simple 1D conceptual model is developed near the Peruvian region, which successfully reproduces the observed relationship between higher moments of SWS and its mean value. The model based physical picture among ocean surface winds, SST, and marine boundary clouds are supported by three-year QuikSCAT surface wind observations and fifteen-year ERA40 re-analysis data. Model sensitive tests suggest that large-scale divergence, and strengths of momentum and cloud fluctuations have significant effects on the ocean SWS-PDF in marine stratus and stratocumulus regions.     

Numerical simulation of cyclonic storms FANOOS, NARGIS with assimilation of conventional and satellite observations using 3-DVAR

In this work, the impact of assimilation of conventional and satellite data is studied on the prediction of two cyclonic storms in the Bay of Bengal using the three-dimensional variational data assimilation (3D-VAR) technique. The FANOOS cyclone (December 6?C10, 2005) and the very severe cyclone NARGIS (April 28?CMay 2, 2008) were simulated with a double-nested weather research and forecasting (WRF-ARW) model at a horizontal resolution of 9?km. Three numerical experiments were performed using the WRF model. The back ground error covariance matrix for 3DVAR over the Indian region was generated by running the model for a 30-day period in November 2007. In the control run (CTL), the National Centers for Environmental Prediction (NCEP) global forecast system analysis at 0.5° resolution was used for the initial and boundary conditions. In the second experiment called the VARCON, the conventional surface and upper air observations were used for assimilation. In the third experiment (VARQSCAT), the ocean surface wind vectors from quick scatterometer (QSCAT) were used for assimilation. The CTL and VARCON experiments have produced higher intensity in terms of sea level pressure, winds and vorticity fields but with higher track errors. Assimilation of conventional observations has meager positive impact on the intensity and has led to negative impact on simulated storm tracks. The QSCAT vector winds have given positive impact on the simulations of intensity and track positions of the two storms, the impact is found to be relatively higher for the moderate intense cyclone FANOOS as compared to very severe cyclone NARGIS.     

Atmospheric mesoscale conditions during the Boothbay meteotsunami: a numerical sensitivity study using a high-resolution mesoscale model

The article aims to test the sensitivity of high-resolution mesoscale atmospheric model to fairly reproduce atmospheric processes that were present during the Boothbay Harbor meteotsunami on 28 October 2008. The simulations were performed by the Weather and Research Forecasting (WRF) model at 1-km horizontal grid spacing by varying initial conditions (ICs) and lateral boundary conditions (LBCs), nesting strategy, simulation lead time and microphysics and convective parameterizations. It seems that the simulations that used higher-resolution IC and LBC were more successful in reproduction of precipitation zone and surface pressure oscillations caused by internal gravity waves observed during the event. The results were very sensitive to the simulation lead time and to the choice of convective parameterization, while the choice of microphysics parameterization and the type of nesting strategy (one-way or two-way) was less important for reproducibility of the event. The success of the WRF model appears limited to very short-range forecasting, most advanced parameterizations, and very high-resolution grid spacing; therefore, the applicability of present atmospheric mesoscale models to future operational meteotsunami warning systems still has a lot of room for improvements.     

Investigation of the atmospheric mesoscale circulation patterns and their simulation with WRF-CHEM model of the dust storm occurrence over the southern coast of the Caspian Sea

For the very first time, the mesoscale circulation patterns and synoptic-dynamic structure of the atmospheric systems that led to the dust emission to the south coast region of Caspian Sea (SCRCS) were identified and classified using the region synoptic stations’ observations of 2005–2013. Satellite measurements and images, NCEP/NCAR reanalysis data, and Weather Research and Forecasting (WRF) model coupled with Chemistry (WRF-CHEM) model products were also used in this study. Results showed that in 49 % of cases, cyclonic circulations over the Middle East deserts were the main transporter of dust particles into the atmosphere where then transferred to the SCRCS by southerly winds over the Alborz mountains in the lower troposphere and by westerly waves in the middle and upper troposphere. During the warm seasons, the surface heating lead to the development of mesoscale thermal low pressures over the hot deserts on the eastern regions of the Caspian Sea, like Turkmenistan and Qura Qum. Those heat lows were responsible for the 38 % of the occurred events. Turbulence and instabilities in the lower troposphere were identified as the second important dust emitter to the atmosphere where those dust particles transported to the SCRCS with the strong northeasterly wind. The third pattern by 13 % of cases was belonging to the mesoscale thermal low pressure that was developed over the arid regions of Iran like Dasht-é-Kavir. Because of the nature of the turbulence in the lower troposphere and heat lows, the ascent of dust particles by these two mechanisms was limited to a shallow layer in the troposphere. The results of simulation with the WRF-CHEM model, analysis of moderate resolution imaging spectroradiometer (MODIS) images, and spatial zoning of atmospheric optical depth (AOD) confirmed the results of the synoptic study.     

A multiple scale modeling system for coastal hurricane wind damage mitigation

Hurricane wind damage constitutes the largest percentage of catastrophic insured losses in the US. Yet the complicated wind structures and their changes are not fully understood and, thus, have not been considered in current risk catastrophic models. To obtain realistic landfall hurricane surface winds, a large eddy simulation (LES) framework in a weather forecasting mode has been developed from a multiple nested Weather Research & Forecasting (WRF) model to explicitly simulate a spectrum of scales from large-scale background flow, hurricane vortex, mesoscale organizations, down to fine-scale turbulent eddies in a unified system. The unique WRF-LES enables the high resolution data to be generated in a realistic environment as a hurricane evolves. In this paper, a simulation of the landfalling Hurricane Katrina is presented to demonstrate various features of the WRF-LES. It shows that the localized damaging winds are caused by the large eddy circulations generated in the hurricane boundary layer. With a sufficient computational power, WRF-LES has the potential to be developed into the next generation operational public wind-field model for hurricane wind damage mitigation.     

Wind-induced salt-wedge intrusion in the Tiber river mouth (Rome–Central Italy)

The wind effect on river water quality was illustrated by means of thermohaline measurements carried out in the Tiber River in May 2012. The survey was carried out using a boat, in stations located in the two Tiber branches: Fiumara Grande and Traiano Canal. Thermohaline variables (salinity and temperature) were used to describe the water-type patterns and to define the salt-wedge position. Although the river flow rate was rather high, saltwater intrusion happened. Wind data suggested that the more probable cause of salt-wedge intrusion was the wind action. Especially wind speeds higher than 4 m/s are able to dominate the sea current at surface layers, determining an increase in the sea level. Therefore, westerly winds determined a seawater inflow in the river mouths.     

新一代中尺度预报模式（WRF）国内应用进展

随着中尺度大气模式的不断发展,新一代中尺度天气研究和预报模式WRF因其完全开放、可移植性强、更新快等特点在国内外得到了广泛应用。从物理参数化方案研究,实时个例模拟研究及与中尺度大气模式MM5的对比研究3个方面介绍近10年来WRF模式在国内的发展和应用概况,阐明WRF模式在中尺度模拟中的普适性和优越性,展望WRF模式在国...     

Energetics of Indian winter monsoon

The Indian subcontinent is characterized by complex topography and heterogeneous land use-land cover. The Himalayas and the Tibetan Plateau are spread across the northern part of the continent. Due to its highly variable topography, understanding of the prevailing synoptic weather systems is complex over the region. The present study analyzes the energetics of Indian winter monsoon (IWM) over the Indian subcontinent using outputs of mesoscale model (MM5) forced with National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), US, initial and boundary conditions. MM5 modeling framework, designed to simulate or predict mesoscale atmospheric circulations, is having a limited-area, non-hydrostatic and terrain following 12 sigma levels. The IWM energetics is studied using MM5 model outputs. Prior to this model’s validity and deviation from the corresponding observations (NCEP/NCAR) is assessed. The model’s overestimation/underestimation of wind, temperature and specific humidity at upper troposphere proves that the model has difficulty in picking up corresponding fields at all the model grid points because of terrain complexity over the Himalayas and Tibetan Plateau. Hence, the model fields deviate from the corresponding observations. However, model results match well with the winter global energy budget calculated using reanalysis dataset by Peixoto and Oort (1992). It suggests MM5 model’s fitness in simulating large scale synoptic weather systems. And, thus the model outputs are used for calculation of energetics associated with IWM. It is observed that beyond \(15^{{\circ }}\hbox {N}\) lower as well as upper level convergence of diabatic heating, which represents continental cooling and sinking of heat from atmosphere to land mass (i.e., surface is cooler than surrounding atmosphere) dominates. The diabatic heating divergence (cooling of continents) is found over ocean/sea and whole of the China region, Tibetan and central Himalayas (because of excess condensation than evaporation). The adiabatic generation of kinetic energy depends on the cross isobaric flow (north to south in winter, i.e., the present study shows strong circulation during IWM). It is found that wind divergence of model concludes lower level convergence over study region (i.e., strong winter circulation in the model fields).