The authors prove that the singular set of a harmonic map from a smooth Riemammian domain to a Riemannian DM-complex is of Hausdorff codimension at least two. They also explore monotonicity formulas and an order gap theorem for approximately harmonic maps. These regularity results have applications to rigidity problems examined in subsequent articles.