Let $\Phi : S \times [0,\infty) \to M$ be Struwe's weak global solution to the heat equation with smooth initial data $\phi : S \to M$, where $S$ is a compact surface and $M$ is a compact three-dimensional Riemannian manifold. Except at finitely many points in space-time, $\Phi$ is smooth. At these singular points, non-constant harmonic maps $S^2 \to M$ "bubble off" in the precise sense described by Struwe in 1985.
My question is about the stability of these bubbles under small perturbations of the initial data. Specifically, let $\phi_k : S \to M$ be a sequence of maps such that $\phi_k \to \phi$ in the $C^\infty$ norm. Suppose that $\Phi_k : S \times [0,\infty) \to M$ is a sequence of weak solutions with $\Phi_k(\cdot,0) = \phi_k(\cdot)$ and that $(x_k,t_k) \in S \times (0,T]$ are singularities of the respective $\Phi_k$ at which bubbles $u_k : S \to M$ form. Suppose also that the $u_k$ converge in $C^\infty$ to a harmonic map $u$. If $(x_k,t_k) \to (x,T)$, is it true that $\Phi$ has a bubble of $u$ at $(x,T)$?