  This work concerns $L^p$ norms of high energy Laplace eigenfunctions,
$(-\Delta_g-\lambda^2)\phi_\lambda=0$, $\|\phi_\lambda\|_{L^2}=1$. In 1988,
Sogge gave optimal estimates on the growth of $\|\phi_\lambda\|_{L^p}$ for a
general compact Riemannian manifold. The goal of this article is to give
general dynamical conditions guaranteeing quantitative improvements in $L^p$
estimates for $p>p_c$, where $p_c$ is the critical exponent. We also apply
previous results of the authors to obtain quantitative improvements in concrete
geometric settings including all product manifolds. These are the first results
improving estimates for the $L^p$ growth of eigenfunctions that only require
dynamical assumptions. In contrast with previous improvements, our assumptions
are local in the sense that they depend only on the geodesics passing through a
shrinking neighborhood of a given set in $M$. Moreover, the article gives a
structure theorem for eigenfunctions which saturate the quantitatively improved
$L^p$ bound. Modulo an error, the theorem describes these eigenfunctions as
finite sums of quasimodes which, roughly, approximate zonal harmonics on the
sphere scaled by $1/\sqrt{\log \lambda}$.
