We study the local regularity of p-caloric functions or more generally of ϕ-caloric functions. In particular, we study local solutions of non-linear parabolic systems with homogeneous right hand side, where the leading terms has Uhlenbeck structure of Orlicz type. This paper closes the gap of [22] where Liebermann proved that if the gradient of a solution is bounded, it is Hölder continuous.The crucial step is a novel local estimates for the gradient of the solutions, which generalize and improve the pioneering estimates of DiBenedetto and Friedman [12,10] for the p-Laplace heat equation.