Let $\Omega$ be a bounded open subset of ${ℝ}^n$, let $X=(x,t)$ be a point of ${ℝ}^n\times {ℝ}^N$. In the cylinder $Q=\Omega \times (-T,0)$, $T>0$, we deduce the local differentiability result $$u\in L^2(-a,0,H^2(B\left(\sigma \right),{ℝ}^N))\cap H^1(-a,0,L^2(B\left(\sigma \right),{ℝ}^N))$$ for the solutions $u$ of the class $L^q(-T,0,H^{1,q}(\Omega ,{ℝ}^N))\cap C^{0,\lambda }(\overline{Q},{ℝ}^N)$ ($0<\lambda <1$, $N$ integer $\ge 1$) of the nonlinear parabolic system $$-\sum _{i=1}^n{D}_i{a}^i(X,u,Du)+{\displaystyle \frac{\partial u}{\partial t}}=B^0(X,u,Du)$$ with quadratic growth and nonlinearity $q\ge 2$. This result had been obtained making use of the interpolation theory and an imbedding theorem of Gagliardo-Nirenberg type for functions $u$ belonging to $W^{1,q}\cap C^{0,\lambda }$.