In this paper, we study inhomogeneous gravitational collapse of viscous matter with spherically symmetric configuration in an universe with positive cosmological constant. The matter fields admit a wide variety of initial density and velocity profiles. We study the evolution of spherical marginally trapped surfaces in this setting, since this allows us to track the local black hole horizons during the collapse process.