On the realism of the rain microphysics representation of a squall line in the WRF Model. Part II: Sensitivity studies on the rain drop size distributions

A comparison between retrieved properties of the rain drop size distributions (DSDs) from multifrequency cloud radar observations and WRF Model results using either the Morrison or the Thompson bulk microphysics scheme is performed in order to evaluate the model's ability to predict the rain microphysics. This comparison reveals discrepancies in the vertical profile of the rain DSDs for the stratiform region of the squall-line system observed on 12 June 2011 over Oklahoma. Based on numerical sensitivity analyses, this study addresses the bias at the top of the rain layer and the vertical evolution of the DSD properties (i.e., of D-m and N0). In this way, the Thompson scheme is used to explore the sensitivity to the melting process. Moreover, using the Thompson and Morrison schemes, the sensitivity of the DSD vertical evolution to different breakup and self-collection parameterizations is studied. Results show that the DSDs are strongly dependent on the representation of the melting process in the Thompson scheme. In the Morrison scheme, the simulations with more efficient breakup reproduce the DSD properties with better fidelity. This study highlights how the inaccuracies in simulated D-m and N0 for both microphysics schemes can impact the evaporation rate, which is systematically underestimated in the model.