On a compact Riemannian manifold $(M,g)$ without boundary, there are two ways to define a Sobolev norm on $M$. Assume that $f\in C^\infty(M)$ in the following.
Use pseudo-differential operators on $M$: $$\|f\|_{H^{s,p}(M)}=\|(I-\Delta_g)^{s/2}f\|_{L^p(M)}.$$
Use local coordinate and Fourier transform in $\mathbb{R}^n$
$$\|f\|_{W^{s,p}(M)}=\sum_{\nu}\|D^sf_\nu\|_{L^p(\mathbb{R}^n)},$$ Here $f_\nu=(\phi_\nu\cdot f)\circ \kappa_\nu^{-1}$, where $\{\phi_\nu\}$ is a partition of unity subordinate to a finite covering $M=\cup \Omega_\nu$, $\kappa_\nu:\Omega_\nu\to \tilde \Omega_\nu\subset \mathbb{R}^n$ is the coordinate map, and $D^sf_\nu$ is defined by the Fourier transform in $\mathbb{R}^n$ $$(D^s g)^\wedge(\xi)=(1+|\xi|^2)^{s/2}\hat g(\xi).$$
It is well-known that these two norms are equivalent when $p=2$. It can be proved by using the $L^2$-boundedness of zero-order pseudo-differential operators. However, I cannot find any reference on the equivalence for $p\ne2$. Any help will be appreciated.