## Abstract

The present study focuses on determining the best combination of microphysics (MP) and cumulus parameterization (CP) schemes for the simulation of Tropical Cyclones (TCs) in the Indian subcontinent region, using the Weather Research and Forecasting (WRF) model. From the available schemes, four CP schemes, namely Kain–Fritsch, Betts–Miller–Janjic, New Simplified Arakawa–Schubert, and Grell–Devenyi, and four MP schemes, namely WSM6, Purdue Lin, Thompson, and Morrison, are selected for the sensitivity study. Seven TCs are simulated using all combinations of the chosen physics schemes. The simulated tracks and intensities are compared against the India Meteorological Department (IMD) observations. The results show that the Kain–Fritsch scheme, in combination with all MP schemes, predicts the tracks best among all the available CP schemes, but the performance of microphysics schemes is indistinguishable. A further study is conducted using the Global Precipitation Measurement (GPM) radar data to identify the best MP scheme by comparing the reflectivities. An existing radar simulator is modified to calculate the simulated reflectivities from the WRF model output corresponding to the MP schemes. The simulated reflectivities are compared against the GPM radar reflectivities, and the results show that the Thompson scheme reproduces the reflectivities closest to observations. The performance of the best set of schemes obtained from this study is compared with a random set of schemes for cyclone Bulbul, and the best set of schemes outperformed the random schemes in every aspect.

This is a preview of subscription content, access via your institution.

## References

Alam M M, Hossain M A and Shafee S 2003 Frequency of Bay of Bengal cyclonic storms and depressions crossing different coastal zones;

*Int. J. Climatol.***23(9)**1119–1125.Balaguru K, Taraphdar S, Leung L R and Foltz G R 2014 Increase in the intensity of post-monsoon Bay of Bengal tropical cyclones;

*Geophys. Res. Lett.***41(10)**3594–3601.Beljaars A C M 2004 The parameterization of surface fluxes in large-scale models under free convection;

*Quart. J. Roy. Meteor. Soc.***121(522)**255–270.Chandrasekar R and Balaji C 2012 Sensitivity of tropical cyclone Jal simulations to physics parameterizations;

*J. Earth Syst. Sci.***121(4)**923–946.Chen S H and Sun W Y 2002 A one-dimensional time dependent cloud model;

*J. Meteorol. Soc. Japan***80(1)**99–118.Choudhury D and Das S 2017 The sensitivity to the microphysical schemes on the skill of forecasting the track and intensity of tropical cyclones using WRF-ARW model;

*J. Earth Syst. Sci.***126****(4)**57.Elsberry R L 2014 Advances in research and forecasting of tropical cyclones from 1963–2013;

*Asia-Pacific J. Atmos. Sci.***50(1)**3–16.Franklin C N, Holland G J and May P T 2005 Sensitivity of tropical cyclone rainbands to ice-phase microphysics;

*Mon. Weather Rev.***133(8)**2473–2493.Grell G A and Dévényi D 2002 A generalized approach to parameterizing convection combining ensemble and data assimilation techniques;

*Geophys. Res. Lett.***29(14)**.Grell G A and Freitas S R

*et al*. 2014 A scale and aerosol aware stochastic convective parameterization for weather and air quality modeling;*Atmos. Chem. Phys.***14(10)**5233–5250.Han J and Pan H L 2011 Revision of convection and vertical diffusion schemes in the NCEP global forecast System;

*Weather Forecast.***26(4)**520–533.Hong S Y, Dudhia J and Chen S H 2004 A revised approach to ice microphysical processes for the bulk parameterization of clouds and precipitation;

*Mon. Weather Rev.***132(1)**103–120.Hong S Y, Noh Y and Dudhia J 2006 A new vertical diffusion package with an explicit treatment of entrainment processes;

*Mon. Weather Rev.***134(9)**2318–2341.Hou A Y, Kakar R K, Neeck S, Azarbarzin A A, Kummerow C D, Kojima M, Oki R, Nakamura K and Iguchi T 2014 The global precipitation measurement mission;

*Bull. Am. Meteorol. Soc.***95(5)**701–722.Huang Y, Wang Y, Xue L, Wei X, Zhang L and Li H 2020 Comparison of three microphysics parameterization schemes in the WRF model for an extreme rainfall event in the coastal metropolitan city of Guangzhou, China;

*Atmos. Res.***240**104939.Iacono M J, Delamere J S, Mlawer E J, Shephard M W, Clough S A and Collins W D 2008 Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models;

*J. Geophys. Res. Atmos.***113(D13)**.IMD 2015 Extremely Severe Cyclonic Storm, ‘CHAPALA’ over the Arabian Sea (28 October–4 November, 2015): A Report–India Meteorological Department, New Delhi.

IMD 2016 Very Severe Cyclonic Storm, ‘Vardah’ over the Bay of Bengal (06–13 December 2016): A Report – India Meteorological Department, New Delhi.

IMD 2017a Cyclonic Storm, ‘Maarutha’ over the Bay of Bengal (15–17 April 2017): A Report – India Meteorological Department, New Delhi.

IMD 2017b Severe Cyclonic Storm, ‘Mora’ over the Bay of Bengal (28–31 May 2017): A Report – India Meteorological Department, New Delhi.

IMD 2018a Severe Cyclonic Storm, ‘\(\text{P}_{{\rm Hethai}}\)’ over southeast Bay of Bengal (13–18 December 2018): A Report – India Meteorological Department, New Delhi.

IMD 2018b Very Severe Cyclonic Storm ‘Titli’ over east central Bay of Bengal (08–13 October 2018): Summary – India Meteorological Department, New Delhi.

IMD 2019a Extremely Severe Cyclonic Storm ‘F

_{ani}’ over east central equatorial Indian Ocean and adjoining southeast Bay of Bengal (26 April–04 May, 2019): Summary – India Meteorological Department, New Delhi.IMD 2019b Very Severe Cyclonic Storm ‘Bulbul’ over the Bay of Bengal (05th–11th November 2019): Summary – India Meteorological Department, New Delhi.

Janjić Z I 1994 The step-mountain eta coordinate model: Further developments of the convection, viscous sublayer, and turbulence closure schemes;

*Mon. Weather Rev.***122(5)**927–945.Kain J S 2004 The Kain–Fritsch convective parameterization: An update;

*J. Appl. Meteorol.***43(1)**170–181.Kanase R D and Salvekar P 2015 Effect of physical parameterization schemes on track and intensity of cyclone LAILA using WRF model;

*Asia-Pacific J. Atmos. Sci.***51(3)**205–227.Koch S E, Ferrier B, Stoelinga M T, Szoke E, Weiss S J and Kain J S 2005 The use of simulated radar reflectivity fields in the diagnosis of mesoscale phenomena from high-resolution WRF model forecasts; In: Preprints, 11th Conf. on Mesoscale Processes, Albuquerque, NM,

*Am. Meteorol. Soc.*J4J,**7**1–9.Kumar S and Bhat G 2016 Vertical profiles of radar reflectivity factor in intense convective clouds in the tropics;

*J. Appl. Meteorol. Climatol.***55(5)**1277–1286.Lim K S S and Hong S Y 2010 Development of an effective double-moment cloud microphysics scheme with prognostic cloud condensation nuclei (CCN) for weather and climate models;

*Mon. Weather Rev.***138(5)**1587–1612.Ma L M and Tan Z M 2009 Improving the behavior of the cumulus parameterization for tropical cyclone prediction: Convection trigger;

*Atmos. Res.***92(2)**190–211.Min K H, Choo S, Lee D and Lee G 2015 Evaluation of WRF cloud microphysics schemes using radar observations;

*Weather Forecast.***30(6)**1571–1589.Mishchenko M I and Travis L D 1998 Capabilities and limitations of a current FORTRAN implementation of the T-matrix method for randomly oriented, rotationally symmetric scatterers;

*J. Quant. Spectrosc. Radiat. Transfer***60(3)**309–324.Mlawer E J, Taubman S J, Brown P D, Iacono M J and Clough S A 1997 Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave;

*J. Geophys. Res. Atmos.***102(D14)**16,663–16,682.Morrison H, Thompson G and Tatarskii V 2009 Impact of cloud microphysics on the development of trailing stratiform precipitation in a simulated squall line: Comparison of one-and two-moment schemes;

*Mon. Weather Rev.***137(3)**991–1007.Mukhopadhyay P, Taraphdar S and Goswami B 2011 Influence of moist processes on track and intensity forecast of cyclones over the North Indian ocean;

*J. Geophys. Res. Atmos.***116(D5)**.Osuri K K, Mohanty U, Routray A, Kulkarni M A and Mohapatra M 2012 Customization of WRF-ARW model with physical parameterization schemes for the simulation of tropical cyclones over north Indian ocean;

*Nat. Hazards***63(3)**1337–1359.Oue M, Tatarevic A, Kollias P, Wang D, Yu K and Vogelmann A M 2019 The Cloud Resolving Model Radar Simulator (CR-SIM) Version 3.2: Description and Applications of a Virtual Observatory;

*Geoscientific Model Development Discussions***2019**1–31, https://doi.org/10.5194/gmd-2019-207; https://www.geosci-model-dev-discuss.net/gmd-2019-207.Pattanaik D and Rao Y R 2009 Track prediction of very severe cyclone ‘Nargis’ using high resolution weather research forecasting (wrf) model;

*J. Earth Syst. Sci.***118(4)**309–329.Pattanayak S, Mohanty U and Osuri K K 2012 Impact of parameterization of physical processes on simulation of track and intensity of tropical cyclone Nargis (2008) with WRF-NMM model;

*Sci. World J.***2012**671437.Pennelly C, Reuter G and Flesch T 2014 Verification of the wrf model for simulating heavy precipitation in Alberta;

*Atmos. Res.***135**172–192.Rao D V B, Srinivas D and Satyanarayana G C 2019 Trends in the genesis and landfall locations of tropical cyclones over the Bay of Bengal in the current global warming era;

*J. Earth Syst. Sci.***128(7)**1–10.Reddy P J, Sriram D, Gunthe S and Balaji C 2021 Impact of climate change on intense Bay of Bengal tropical cyclones of the post-monsoon season: A pseudo global warming approach;

*Clim. Dyn.***56(9)**2855–2879.Roberts N M and Lean H W 2008 Scale-selective verification of rainfall accumulations from high-resolution forecasts of convective events;

*Mon. Weather Rev.***136(1)**78–97.Rogers R F, Black M L, Chen S S and Black R A 2007 An evaluation of microphysics fields from mesoscale model simulations of tropical cyclones. Part I: Comparisons with observations;

*J. Atmos. Sci.***64(6)**1811–1834.Ryzhkov A, Pinsky M, Pokrovsky A and Khain A 2011 Polarimetric radar observation operator for a cloud model with spectral microphysics;

*J. Appl. Meteorol. Climatol.***50(4)**873–894.Sandeep C, Krishnamoorthy C and Balaji C 2018 Impact of cloud parameterization schemes on the simulation of cyclone Vardah using the WRF model;

*Curr. Sci.***115(6)**.Skamarock W C, Klemp J B, Dudhia J, Gill D O, Barker D M, Wang W and Powers J G 2008 A description of the Advanced Research WRF version 3 NCAR Technical note-475+ STR.

Srinivas C, Bhaskar Rao D, Yesubabu V, Baskaran R and Venkatraman B 2013 Tropical cyclone predictions over the Bay of Bengal using the high-resolution advanced research weather research and forecasting (arw) model;

*Quart. J. Roy. Meteorol. Soc.***139(676)**1810–1825.Starzec M, Mullendore G L and Kucera P A 2018 Using radar reflectivity to evaluate the vertical structure of forecast convection;

*J. Appl. Meteorol. Climatol.***57(12)**2835–2849.Sun Y, Dong X, Cui W, Zhou Z, Fu Z, Zhou L, Deng Y and Cui C 2020 Vertical structures of typical meiyu precipitation events retrieved from GPM-DPR;

*J. Geophys. Res. Atmos.***125(1)**e2019JD031466.Taraphdar S, Mukhopadhyay P, Leung L R, Zhang F, Abhilash S and Goswami B 2014 The role of moist processes in the intrinsic predictability of Indian ocean cyclones;

*J. Geophys. Res. Atmos.***119(13)**8032–8048.Taylor K E 2001 Summarizing multiple aspects of model performance in a single diagram;

*J. Geophys. Res. Atmos.***106(D7**) 7183–7192.Tewari M, Chen F, Wang W, Dudhia J, LeMone M A, Mitchell K, Ek M, Gayno G, Wegiel J and Cuenca R H 2004 Implementation and verification of the unified NOAH land surface model in the WRF model; 20th conference on weather analysis and forecasting/16th conference on numerical weather prediction

**1115(6)**.Thompson G, Field P R, Rasmussen R M and Hall W D 2008 Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part II: Implementation of a new snow parameterization;

*Mon. Weather Rev.***136(12)**5095–5115.Yamaguchi M, Ishida J, Sato H and Nakagawa M 2017 Wgne intercomparison of tropical cyclone forecasts by operational NWP models: A quarter century and beyond;

*Bull. Am. Meteorol. Soc.***98(11)**2337–2349.Yano J I, Geleyn J F, Köhler M, Mironov D, Quaas J, Soares P M, Phillips V T, Plant R S, Deluca A and Marquet P

*et al*. 2015 Basic concepts for convection parameterization in weather forecast and climate models: COST Action ES0905 final report;*Atmosphere***6(1)**88–147.Zhu T and Zhang D L 2006 Numerical simulation of hurricane bonnie (1998). Part II: Sensitivity to varying cloud microphysical processes;

*J. Atmos. Sci.***63(1)**109–126.

## Acknowledgements

The authors are thankful to the international WRF community for their tremendous efforts in developing the WRF model and supporting us through providing the model. The model simulations are performed on Aqua High-Performance Computing (HPC) system at the Indian Institute of Technology Madras, Chennai, India. The figures are plotted using the python language and NCAR Command Language (NCL). The authors would like to thank the ‘Centre of Excellence (CoE) in climate change impact on coastal infrastructure and the adaptation strategies’ for their suggestions. Department of Science and Technology is funding the research grant DST/CCP/CoE/141/2018(C) for this project (CIE1819265DSTXSACI) under SPLICE – Climate Change Programme.

## Author information

### Affiliations

### Contributions

Harish Baki and Sandeep Chinta conceptualized the study, designed the methodology, performed the formal analysis and edited the draft of the paper. Harish Baki performed the investigation and prepared the original draft. C Balaji and Balaji Srinivasan provided critical revisions for the manuscript and supervised the work.

### Corresponding author

## Additional information

Communicated by Parthsarathi Mukhopadhyay

## Appendix

### Appendix

### A1. Drop size distributions of hydrometeors in the Lin scheme

The Purdue Lin scheme adopted classical exponential drop size distribution for hydrometeors rain, snow, and graupel as:

Here, \(n(D)\) is the drop size distribution as a function of the diameter \(D\), \(dD\) is the differential diameter, \(N_0\) is the intercept parameter, and \(\lambda \) is the slope parameter. The intercept and slope parameter values will vary for different hydrometeors.

#### A1.1 Rain

#### A1.2 Snow

#### A1.3 Graupel

In the above equations, \(Q_R,Q_S,Q_G\) are the mixing ratios (in kg/kg) of the hydrometeors rain, snow, and graupel, \(\rho _R, \rho _S, \rho _G\) are the densities, and \(\rho _{\text {dair}}\) is the dry air density (in kg/m\(^3\)).

### A2. Drop size distributions of hydrometeors in the WSM6 scheme

The drop size distributions of the hydrometeors adopted in the WSM6 scheme are based on the Lin scheme, and thus follow the same exponential distribution, with a little modification adopted for snow.

#### A2.1 Rain

#### A2.2 Snow

Here, \(T\) is the temperature (K) and \(T_0\) is the reference temperature (273.16 K).

#### A2.3 Graupel

### A3. Drop size distributions of hydrometeors in the Thompson scheme

The Thompson scheme assumes generalized gamma distribution for each hydrometeor (except for snow) as:

Here, \(N_t\) is the total number concentration, \(N_x\) is the prognostic number concentration (in 1/kg) of hydrometeor \(x\), \(Q_x\) is the mixing ratio of hydrometeor \(x\), \(a_m \;{\text and}\;b_m\) are the constant and exponents of mass-size relationship, \(\Gamma \) is the gamma function, and \(\mu \) is the shape parameter. When \(\mu =0\), the distribution becomes the classical exponential distribution. In the Thompson scheme, along with the mixing ratios, the number concentrations of rain and ice particles are also simulated. Thus, the scheme is known to be a double-moment scheme. Further information can be found at Oue *et al.* (2019).

#### A3.1 Rain

For rain particles, the shape parameter is 0. Thus, the drop size distributions becomes classical exponential distribution. The scheme uses the model simulated prognostic rain number concentration in the distribution.

to

Here, \(N_R\) is the number concentration (in 1/kg).

#### A3.2 Snow

The Thompson scheme specially build for snow. Unlike the exponential or generalized gamma distribution, it follows the combination of both.

where \(K_0 = 490.6, K_1 = 17.56, \wedge _0 = 20.78, \wedge _1 = 3.29, \mu _s = 0.6357, \text{ and } \, T_c \) is the air temperature in °C, \(\rho _S = 100\,\text{kg/m}^3. \)

#### A3.3 Graupel

For graupel also, the shape parameter is 0.

### A4. Drop size distributions of hydrometeors in the Morrison double moment scheme

The Morrison scheme also adopted generalized gamma distribution for every hydrometeor. The scheme simulated number concentrations for all the hydrometeors, except cloud water. The shape parameter \(\mu \) is considered zero for rain, snow, and graupel. Thus, the generalized gamma distribution simplified to an exponential distribution, with \(N_t\lambda \) being considered as intercept. In addition, the constant and exponents of mass-diameter parameters are \(a_m = \frac{\pi \rho _x}{6}\) and \(b_m = 3\).

#### A4.1 Rain

#### A4.2 Snow

#### A4.3 Graupel

In the above equations, \(N_R, N_S, N_G\) are the number concentrations (in 1/kg) for rain, snow, and graupel, respectively.

## Rights and permissions

## About this article

### Cite this article

Baki, H., Chinta, S., Balaji, C. *et al.* A sensitivity study of WRF model microphysics and cumulus parameterization schemes for the simulation of tropical cyclones using GPM radar data.
*J Earth Syst Sci* **130, **190 (2021). https://doi.org/10.1007/s12040-021-01682-3

Received:

Revised:

Accepted:

Published:

### Keywords

- Weather Research and Forecasting (WRF) model
- microphysics schemes
- cumulus physics schemes
- GPM radar data
- forward radar operator
- cloud resolving model
- radar simulator (CR-SIM)