Using Auton's force law for the unsteady motion of a spherical bubble in inhomogeneous unsteady flow, two key dimensionless groups are deduced which determine whether isolated vortices or shear-layer vortices can trap bubbles. These groups represent the ratio of inertial to buoyancy forces as a relaxation parameter [tcy ] = ΔU2/2gx and a trapping parameter [Gcy ] = ΔU/VT where ΔU is the velocity difference across the vortex or the shear layer, x is streamwise distance measured from the effective origin of the mixing layer and VT is the terminal slip speed of the bubble or particle. It is shown here that whilst buoyancy and drag forces can lead to bubbles moving in closed orbits in the vortex flows (either free or forced), only inertial forces result in convergent trajectories. Bubbles converge on the downflow side of the vortex at a location that depends on the inertial and lift forces. It is important to note that the latter have been omitted from many earlier studies.A discrete-vortex model is used to simulate the large-scale unsteady flows within horizontal and vertical mixing layers between streams with velocity difference ΔU. Trajectories of non-interacting small bubbles are computed using the general force law. In the horizontal mixing layer it is found that Γ needs to have a value of about 3 to trap about 50% of the bubbles if Π is about 0.5 and greater if Π is less. The pairing of vortices actually enhances their trapping of bubbles. In the vertical mixing layer bubbles are trapped mainly within the growing vortices but bubbles are concentrated on the downflow side of the vortices as Γ and Π increase. In a companion paper we show that lateral dispersion of bubbles can be approximately described by an advective diffusion equation with the diffusivity about equal to the eddy viscosity, i.e. rather less than the diffusivity of heat or other passive scalars.