  We prove $L^p\to L^q$ estimates for local maximal operators associated with
dilates of codimension two spheres in Heisenberg groups; these are sharp up to
two endpoints. The results can be applied to improve currently known bounds on
sparse domination for global maximal operators. We also consider lacunary
variants, and extensions to M\'etivier groups.
