Abstract. Below-cloud scavenging (BCS) is the process of aerosol removal from the atmosphere between cloud base and the ground by precipitation (e.g., rain or snow), and affects aerosol number or mass concentrations, size distribution, and lifetime. An accurate representation of precipitation phases is important in treating BCS as the efficiency of aerosol scavenging differs significantly between liquid and solid precipitation. The impact of different representations of BCS on existing model biases was examined through implementing a new aerosol BCS scheme in the Environment and Climate Change Canada (ECCC) air quality prediction model GEM-MACH and comparing it with the existing scavenging scheme in the model. Further, the current GEM-MACH employs a single-phase precipitation for BCS: total precipitation is treated as either liquid or solid depending on a fixed environment temperature threshold. Here, we consider co-existing liquid and solid precipitation phases as they are predicted by the GEM microphysics. GEM-MACH simulations, in a local-area domain over the Athabasca oil sands areas, Canada, are compared with observed precipitation samples, with a focus on the particulate base cation NH4+; acidic anions NO3-, SO4=, HSO3- in precipitation; and observed ambient particulate sulfate, ammonium, and nitrate concentrations. Overall, the introduction of the multi-phase approach and the new scavenging scheme enhances GEM-MACH performance compared to previous methods. Including a multi-phase approach leads to altered SO4= scavenging and impacts the BCS of SO2 into the aqueous phase over the domain. Sulfate biases improved from +46 % to −5 % relative to Alberta Precipitation Quality Monitoring Program wet sulfate observations. At Canadian Air and Precipitation Monitoring Network stations the biases became more negative, from −10 % to −30 % for the tests carried out here. These improvements contrast with prior annual average biases of +200 % for SO4=, indicating enhanced model performance. Improvements in model performance (via scores for correlation coefficient, normalized mean bias, and/or fractional number of model values within a factor of 2 of observations) could also be seen between the base case and the two simulations based on multi-phase partitioning for NO3-, NH4+, and SO4=. Whether or not these improvements corresponded to increases or decreases in NO3- and NH4+ wet deposition varied over the simulation region. The changes were episodic in nature – the most significant changes in wet deposition were likely at specific geographic locations and represent specific cloud precipitation events. The changes in wet scavenging resulted in a higher formation rate and larger concentrations of atmospheric particle sulfate.