The parametric finite element methods of the Barrett–Garcke–Nürnberg (BGN) type have been successful in preventing mesh distortion/ degeneration in approximating the evolution of surfaces under various geometric flows, including mean curvature flow, Willmore flow, Helfrich flow, surface diffusion, and so on. However, the rigorous justification of convergence of the BGN-type methods and the characeterization of the particle trajectories produced by these methods still remain open since this class of methods was proposed in 2007. The main difficulty lies in the stability of the artificial tangential velocity implicitly determined by the BGN methods. In this paper, we give the first proof of convergence of a stabilized BGN method for curve shortening flow, with optimal-order convergence in
L
2
L^2
norm for finite elements of degree
k
≥
2
k \geq 2
under the stepsize condition
τ
≤
c
h
k
+
1
\tau \leq c h^{k+1}
(for any fixed constant
c
c
). Moreover, we give the first rigorous characterization of the particle trajectories produced by the BGN-type methods for one-dimensional curves, i.e., we prove that the particle trajectories produced by the stabilized BGN methods converge to the particle trajectories determined by a system of geometric partial differential equations which differs from the standard curve shortening flow by a tangential motion. The characterization of the particle trajectories also rigorously explains, for one-dimensional curves, why the BGN-type methods could maintain the quality of the underlying evolving mesh.