The transport and removal of particles suspended in spherical bubbles rising in a liquid pool are modeled based on a hybrid Eulerian-Monte-Cario approach. The bubble internal circulation is assumed to be similar to Hill's vortex flow, and equations of motion are solved for a large number of particles accounting for Brownian dispersion, sedimentation, convective, and inertial mechanisms. Particles that impact the bubble surface are assumed to be removed from the bubble. Parametric calculations are presented, and are compared with relevant Eulerian theoretical models. The strength of the bubble internal circulation is shown to be an important parameter affecting the bubble's particle removal rate. Slower internal circulation, which can result from surface active contaminants, will significantly reduce the removal rate of particles by weakening the inertial particle transport mechanism. (C) 2003 Elsevier Ltd. All rights reserved.