  In this paper we consider sequences $u_j:B_2\subseteq M\to N$ of stationary
harmonic maps between smooth Riemannian manifolds with uniformly bounded energy
$E[u_j]\equiv \int |\nabla u_j|^2\leq \Lambda$ . After passing to a subsequence
it is known one can limit $u_j\to u:B_1\to N$ with the associated defect
measure $|\nabla u_j|^2 dv_g \to |\nabla u|^2dv_g+\nu$, where $\nu = e(x)\,
H^{m-2}_S$ is an $m-2$ rectifiable measure \cite{lin_stat}. For a.e. $x\in
S=\operatorname{supp}(\nu)$ one can produce a finite number of bubble maps
$b_j:S^2\to N$ by blowing up the sequence $u_j$ near $x$.
  We prove the energy identity in this paper. Namely, we have at a.e. $x\in S$
that $e(x)=\sum_j E[b_j]$ for a complete set of such bubbles. That is, the
energy density of the defect measure $\nu$ is precisely the sum of the energies
of the bubbling maps.
