For a broad class of complete hybrid systems evolving on Riemannian manifolds and satisfying mild regularity conditions on the data we introduce a notion of multistability based on the existence of a finite number of compact, globally attractive, and weakly invariant sets. Such notion not only generalizes the standard global uniform asymptotic stability requirement, but can also be characterized in terms of equivalent asymptotic properties, existence of smooth Lyapunov functions, and intrinsic robustness to small perturbations.