We introduce an adaptive multilevel method for solving elliptic eigen-problems on surfaces. Using techniques from adaptive finite element methods, Rayleigh quotient minimization (RQMG) and bootstrap algebraic multigrid we design an efficient algorithm for computing any eigenpair. Numerical results of the proposed approach applied to the Laplace-Beltrami eigen-problem on the sphere are presented.