<p style='text-indent:20px;'>We consider variational discretization [<xref ref-type="bibr" rid="b18">18</xref>] of a parabolic optimal control problem governed by space-time measure controls. For the state discretization we use a Petrov-Galerkin method employing piecewise constant states and piecewise linear and continuous test functions in time. For the space discretization we use piecewise linear and continuous functions. As a result the controls are composed of Dirac measures in space-time, centered at points on the discrete space-time grid. We prove that the optimal discrete states and controls converge strongly in <inline-formula><tex-math id="M1">\begin{document}$ L^q $\end{document}</tex-math></inline-formula> and weakly-<inline-formula><tex-math id="M2">\begin{document}$ * $\end{document}</tex-math></inline-formula> in <inline-formula><tex-math id="M3">\begin{document}$ \mathcal{M} $\end{document}</tex-math></inline-formula>, respectively, to their smooth counterparts, where <inline-formula><tex-math id="M4">\begin{document}$ q \in (1,\min\{2,1+2/d\}] $\end{document}</tex-math></inline-formula> is the spatial dimension. Furthermore, we compare our approach to [<xref ref-type="bibr" rid="b8">8</xref>], where the corresponding control problem is discretized employing a discontinuous Galerkin method for the state discretization and where the discrete controls are piecewise constant in time and Dirac measures in space. Numerical experiments highlight the features of our discrete approach.