This paper investigates the free vibration analysis of a cylindrical shell based on first-order shear deformation theory (FSDT) and nonlocal strain gradient theory (NSGT). For the first time, the surface effects are taken into account for the proposed model. For this purpose, governing equations of motion are derived using Hamilton's principle and Gurtin-Murdoch surface elasticity model. Then, governing equations and corresponding boundary conditions are discretized by the generalized differential quadrature method and frequency responses are obtained. In the results section, differences between the obtained natural frequencies with and without considering the surface effects are reported to show their influence on results. Furthermore, the effects of each surface properties including surface elastic properties, residual surface stress, and surface mass density are investigated in different values of size-dependent parameters.