We present three recently proposed kernel-based collocation methods in unified notations as an easy reference for practitioners who need to solve PDEs on surfaces S ⊂ R d. These PDEs closely resemble their Euclidean counterparts, except that the problem domains change from bulk regions with a flat geometry of some surfaces, on which curvatures play an important role in the physical processes. First, we present a formulation to solve surface PDEs in a narrow band domain containing the surface. This class of numerical methods is known as the embedding types. Next, we present another formulation that works solely on the surface, which is commonly referred to as the intrinsic approach. Convergent estimates and numerical examples for both formulations will be given. For the latter, we solve both the linear and nonlinear time-dependent parabolic equations on static and moving surfaces. Keywords: kernel-based collocation methods, elliptic partial differential equations on manifolds, convergence estimate. for some surface differential operator L S : H m (S) → H m−2 (S) with W m ∞ (S)-bounded coefficients a S , c S : S → R, and b S : S → R d. We assume the existence [2] of classical solutions u * S to (2) in Hilbert spaces H m (S). 2 EMBEDDING KANSA METHODS Methods in this section are variants of Kansa methods [3], [4] that built upon some constantalong-normal property. They are generalization of the finite difference based closest point method [5] and its meshfree extension [6]. Firstly, we define the closest point mapping cp cp(x) = arg inf ξ∈S ξ − x 2 (R d) ,