We extend many known results for harmonic maps from the 2-sphere into a Grassmannian to harmonic maps of finite uniton number from an arbitrary Riemann surface. Our method relies on a new theory of nilpotent cycles arising from the diagrams of F.E. Burstall and the second author associated to such harmonic maps; these properties arise from a criterion for finiteness of the uniton number found recently by the authors with A. Aleman. Applications include a new classification result on minimal surfaces of constant curvature and a constancy result for finite type harmonic maps.