We report a concentration instability at zero Reynolds number created by hydrodynamically interacting bubbles with surfactant. This instability is driven by Marangoni stresses that force bubbles to migrate in directions perpendicular to gravity. We characterize the lateral motion of a single buoyant bubble when it is subject to a weak, low wavenumber disturbance velocity. We use this result to determine which mean flow wavevectors amplify concentration fluctuations in a dilute suspension. The suspension is linearly unstable at small horizontal wavenumbers by a mechanism similar to the concentration instabilities demonstrated in suspensions of sedimenting nonspherical or deformable particles.It is now well known that hydrodynamic interactions between particles under the action of a body force can greatly alter a suspension's microstructure and even create concentration instabilities, particularly when the particles are characterized by an intrinsic microstructural variable such as orientation or shape deformation. 1-5 These instabilities share a common mechanism in the regime of dilute concentration and low particle Reynolds number, 1-5 and they are very important in describing the rates of sedimentation of fibers in suspension. 2,4,5 We describe another example of this instability that is perhaps not obviously related to the previous work in this area. As described below, buoyant bubbles contaminated with surfactant are subject to a concentration instability due to the Marangoni forces induced by hydrodynamic interactions.The simplest mathematical description of these instabilities is through a mean field theory, where the particles only interact through the mean flow. We assume the microstructural variable ͑deformation, orientation, etc.͒ changes instantaneously compared to the growth rate of concentration fluctuations caused by this variable ͑an assumption we justify at the end of the analysis͒. In this situation, the instability of a dilute, monodisperse suspension at low Reynolds number can be captured by a particle advection-diffusion equation coupled with an averaged momentum equation 2In Eq. ͑1͒, is the suspension concentration, F is the body force on a particle, M is the mobility of a particle, D is a hydrodynamic dispersion tensor, and uЈ is the mean disturbance velocity of the fluid caused by all particles. uЈ satisfies an averaged Stokes flow ͑2͒ which includes a body force proportional to the suspension concentration. 1-3 Hydrodynamic interactions occur between particles via uЈ, and these interactions amplify concentration fluctuations when the advection flux M · F + uЈ is larger than the diffusive flux −D · ٌ, and both point in opposite directions. Because of the microstructural variable describing the mobility of the particle, the particle's mobility is coupled to the mean velocity gradient via M = fٌ͑uЈ͒, and this coupling is what drives instability. For weakly deformable particles, M = M 0 ͑1+E͒, where E is the local rate of strain ͓i.e., 1 2 ٌ͑uЈ + ٌuЈ T ͔͒, and is a relaxation time-scale for the...