We calculate the damping of low-lying collective modes of a trapped Bose gas in the hydrodynamic regime, and show that this comes solely from the shear viscosity, since the contributions from bulk viscosity and thermal conduction vanish. The hydrodynamic expression for the damping diverges due to the failure of hydrodynamics in the outer parts of the cloud, and we take this into account by a physically motivated cutoff procedure. Our analysis of available experimental data indicates that higher densities than have yet been achieved are necessary for investigating hydrodynamic modes above the Bose-Einstein transition temperature.PACS numbers: 03.75. Fi, 05.30.Jp, 67.40.Db In recent experiments on magnetically-trapped atomic vapors, alkali atoms [1][2][3] have been cooled to temperatures at which they are degenerate and indeed BoseEinstein condensation has been observed in them. Frequencies and damping rates of collective modes in these systems have been investigated, both above and below the Bose-Einstein transition temperature, T c [4][5][6]. In this Letter we shall focus on properties above T c . One can distinguish two regimes, the hydrodynamic one, for which the characteristic mode frequency is small compared with the collision frequency and the wavelength of the mode is large compared with the atomic mean free path, and the opposite limit, the collisionless one, for which collisions are relatively unimportant. The frequencies of modes in the hydrodynamic regime have been calculated in Ref. [7], and here we calculate their damping.We begin by giving a simple derivation of the basic hydrodynamic equations. Our treatment is essentially that of Ref. [8] generalized to take into account the potential of the trap, and, since we are interested in small oscillations, we shall consider the linearized equations. The Euler equation for the fluid velocity v(r, t) iswhere m is the mass of the atoms, n(r, t) is the particle density, n 0 (r) is the equilibrium particle density, p(r, t) is the pressure and f is the force per unit mass due to the external potential U 0 (r), f = −∇U 0 (r)/m. In equilibrium, where the pressure is p 0 (r), Eq. (1) implies that ∇p 0 (r) = mn 0 (r)f . Taking the time-derivative of Eq.(1) and using the continuity equation, one findsWe calculate the first term on the right hand side of Eq.(2) by using the energy conservation condition [8],where ρ is the mass density and ǫ and w are, respectively, the internal energy and the enthalpy of the fluid per unit mass. Since we assume local thermodynamic equilibrium and neglect contributions to the energy due to interparticle interactions, we may use the results p = ρ(w − ǫ) and ρǫ = 3p/2, and find thatCombining Eqs. (2) and (4), and using the fact that ∇n 0 (r) is proportional to ∇U 0 (r), we obtain for the equation of motion for v(r, t),Equation (5) has previously been derived by Griffin et al.[7] using kinetic theory. In our present discussion we assume that the potential is axially symmetric,where ω 0 is the frequency of the trap in the x − y plane; ...