As pointed out already by the comments, Sogge has indeed made a lot of contributions in this area.
Consider a 2-dimensional compact Riemannian manifold without boundary, then the $L^2$ normalized eigenfunction of Laplacian  $e_{\lambda}$ which sastisfy
$$
-\Delta e_{\lambda}=\lambda^2e_{\lambda}
$$
has the following estimate
$$
\|e_{\lambda}\|_{L^p}\leq C\lambda^{\sigma(p)},~~\lambda\ge 1,
$$
here $\sigma(p)=\frac{1}{2}(\frac{1}{2}-\frac{1}{p})$, if $2\leq p\leq 6$, and $\sigma(p)=2(\frac{1}{2}-\frac{1}{p})-1$, if $6\leq p\leq \infty$. For dimension $n>2$, there are similar results. 

However, for manifold with boundary, the problem becomes more difficult, for dimension 2, see also the Acta paper by Smith and Sogge http://arxiv.org/pdf/math/0605682.pdf.

It's also interesting to improve the bound above. The above estimates is sharp in general. The $L^{\infty}$ norm is saturated by the Zonal function (which is concentrated near the pole) on the round sphere. For lower $p$(below the critical point, which is 6 here), it's saturated by the highest weight hamonics (which is concentrated near the equator ). However, for manifold with nonpositive curvature, on can get some improvement($\log \lambda$), and for flat torus $\mathbb{T}^n$, one can get a further improvement ($\lambda^{\epsilon(n)})$.