The purpose of this paper is to assess the prediction of particle properties of aggregates and particle size distributions with the Spectral Ice Habit Prediction System (SHIPS) and to investigate the effects of crystal habits on aggregation process. Aggregation processes of ice particles are critical to the understanding of precipitation and the radiative signatures of cloud systems. Conventional approaches taken in cloud-resolving models (CRMs) are not ideal to study the effects of crystal habits on aggregation processes because the properties of aggregates have to be assumed beforehand. As described in Part III, SHIPS solves the stochastic collection equation along with particle property variables that contain information about crystal habits and maximum dimensions of aggregates. This approach makes it possible to simulate properties of aggregates explicitly and continuously in CRMs according to the crystal habits.
The aggregation simulations were implemented in a simple model setup, assuming seven crystal habits and several initial particle size distributions (PSDs). The predicted PSDs showed good agreement with observations after rescaling except for the large-size end. The ice particle properties predicted by the model, such as the mass-dimensional (m-D) relationship and the relationship between diameter of aggregates and number of component crystals in an aggregate, were found to be quantitatively similar to those observed. Furthermore, these predictions were dependent on the initial PSDs and habits. A simple model for the growth of a particle's maximum dimension was able to simulate the typically observed fractal dimension of aggregates when an observed value of the separation ratio of two particles was used. A detailed analysis of the collection kernel indicates that the m-D relationship unique to each crystal habit has a large impact on the growth rate of aggregates through the cross-sectional area or terminal velocity difference, depending on the initial equivalent particle distribution. A significant decrease in terminal velocity differences was found in the inertial flow regime for all the habits but the constant-density sphere. It led to formation of a local maximum in the collection kernel and, in turn, formed an identifiable mode in the PSDs. Remaining issues that must be addressed in order to improve the aggregation simulation with the quasi-stochastic model are discussed.