In this talk, we derive global gradient estimates for parabolic $p$-Laplace type equations with measure data, whose model is $$u_t - \textrm{div} \left(|Du|^{p-2} Du\right) = \mu \quad \textrm{in} \ \Omega \times (0,T) \subset \mathbb{R}^n \times \mathbb{R},$$ where $\mu$ is a signed Radon measure with finite total mass. We consider the singular case $$\frac{2n}{n+1}