On the regularity of weak solutions to parabolic systems in two spatial dimensions