Numerical simulations to investigate the e ects of shear ow on wakes behind a sphere were carried out using a three-dimensional ÿnite element method. The time-dependent Navier-Stokes equations were solved by a modiÿed explicit time integration scheme. The Reynolds number, deÿned by the sphere diameter and oncoming velocity at its center, ranged from 20 to 500. Linear shear is assumed, and the dimensionless shear rate, deÿned as the vertical gradient of oncoming shear ow, was varied from 0.0 to 0.2. Extensive comparisons were made between the present computed results and available experimental and numerical investigations, and showed that they are in close agreement. Computed lift force was directed toward the low velocity side for Re ¿ 50 whereas it was toward the high velocity side for Re ¡ 50. In addition, drag forces showed slight increases compared to those in a uniform ow. The Reynolds number, at which vortex shedding begins, was found to be lower than that in a uniform ow and was a function of the shear rate. Vortex shedding frequency also increases with the Reynolds number and the shear rate. The irregularity of the vortex shedding and the oscillation of the vortex detachment location appearing in uniform and low shear ows at Re = 500 vanishes when high shear was present.