Microphysical Process Comparison of Three Microphysics Parameterization Schemes in the WRF Model for an Idealized Squall-Line Case Study

Abstract

Abstract Three bulk microphysics schemes with different complexities in the Weather Research and Forecasting Model are compared in terms of the individual microphysical process terms of the hydrometeor mass and number mixing ratio tendency equations in an idealized 2D squall-line case. Through evaluation of these process terms and of hydrometeor size distributions, it is shown that the differences in the simulated population characteristics of snow, graupel, and rainwater are the prominent factors contributing to the differences in the development of the simulated squall lines using these schemes. In this particular case, the gust front propagation speed produced by the Thompson scheme is faster than in the other two schemes during the first 2 h of the simulation because it has a larger dominant graupel size. After 2 h into the simulation, the initially less intense squall lines in the runs using the WSM6 and Morrison schemes start to catch up in intensity and development to the run using the Thompson scheme. Because the dominant size of graupel particles in the runs using the WSM6 and Morrison schemes is smaller, these particles take more time to fall below the freezing level and enhance the rainwater production and its evaporative cooling. In the run using the Thompson scheme, the graupel production slows down at later times while the snow particle growth increases, leading to more snow falling below the freezing level to melt and surpass graupel particle melting in the production of rainwater.