We establish global-in-time frequency localized local smoothing estimates for Schrödinger equations on hyperbolic space H d . In the presence of symmetric first and zeroth order potentials, which are possibly timedependent, possibly large, and have sufficiently fast polynomial decay, these estimates are proved up to a localized lower order error. Then in the timeindependent case, we show that a spectral condition (namely, absence of threshold resonances) implies the full local smoothing estimates (without any error), after projecting to the continuous spectrum. In the process, as a means to localize in frequency, we develop a general Littlewood-Paley machinery on H d based on the heat flow. Our results and techniques are motivated by applications to the problem of stability of solitary waves to nonlinear Schrödingertype equations on H d .As a testament of the robustness of approach, which is based on the positive commutator method and a heat flow based Littlewood-Paley theory, we also show that the main results are stable under small time-dependent perturbations, including polynomially decaying second order ones, and small lower order nonsymmetric perturbations.