Abstract Titan’s stratosphere has been observed in a superrotation state, where the atmosphere rotates many times faster than the surface does. Another characteristic of Titan’s atmosphere is the presence of a thick haze layer. In this paper, we performed numerical experiments using a general circulation model to explore the effects of the haze layer on the stratospheric superrotation. We employed a semigray radiation model of Titan’s atmosphere following McKay et al., which takes account of sunlight absorption by the haze particles. The phase change of methane or seasonal changes were not taken into account. Our model with radiation parameters tuned for Titan yielded a global eastward wind around the equator with larger velocities at higher altitudes, except at around 70 km, after 10⁵ Earth days. Although the atmosphere is not in an equilibrium state, the zonal wind profiles are approximately consistent with the observed one. By changing the parameters of the radiation model, we found that the intensity and the location of the maximum zonal wind velocity highly depended on the optical thickness and the altitude of the haze layer, respectively. Analysis of our experiments suggests that the quasi-stationary stratospheric superrotation is maintained by the balance between the meridional circulation decoupled from the surface and the eddies that transport angular momentum equatorward. This is different from, but similar to, the so-called Gierasch mechanism, in which momentum is supplied from the surface. This structure may explain the no-wind region at about 80 km in altitude.