A statistical-numerical aerosol parameterization was incorporated into the regional air quality model CMAQ to study the coagulation mixing process focusing on a dust storm event occurred over East Asia. Simulation results show that the coagulation mixing process tends to decrease aerosol mass, surface area, and number concentrations over the dust source areas. Over the downwind oceanic areas, aerosol concentrations generally increased due to enhanced sedimentation as particles became larger upon coagulation. The mixture process can reduce the overall single scattering albedo by up to 10% as a result of enhanced core-with-shell absorption by dust and reduction in the number of scattering particles. The enhanced dry deposition speed also altered the vertical distribution. In addition, the ability of aerosol particles to serve as cloud condensation nuclei (CCN) increased from around 10^7 m^-3 to above 10^9 m^-3 over downwind areas because a large amount of mineral dust particles became effective CCN with solute coating, except over the highly polluted areas where multiple collections of hygroscopic particles by dust in effect reduced CCN number. This CCN effect is much stronger for coagulation mixing than by the uptake of sulfuric acid gas on dust, although the nitric acid gas uptake was not investigated. The ability of dust particles to serve as ice nuclei may decrease or increase at low or high sub-zero temperatures, respectively, due to the switching from deposition nucleation to immersion freezing or haze freezing.