For a compact connected Riemannian $n$-manifold $(\Omega,g)$ with smooth

boundary, we explicitly calculate the first two coefficients $a_0$ and $a_1$ of

the asymptotic expansion of $\sum_{k=1}^\infty e^{-t \tau_k^\mp}= a_0t^{-n/2}

\mp a_1 t^{-(n-1)/2}+a_2^\mp t^{-(n-2)/2} +\cdots+ a_m^\mp t^{-(n-m)/2}

+O(t^{-(n-m-1)/2})$ as $t\to 0^+$, where $\tau^-_k$ (respectively, $\tau^+_k$)

is the $k$-th Navier-Lam\'{e} eigenvalue on $\Omega$ with Dirichlet

(respectively, Neumann) boundary condition. These two coefficients provide

precise information for the volume of the elastic body $\Omega$ and the surface

area of the boundary $\partial \Omega$ in terms of the spectrum of the

Navier-Lam\'{e} operator. This gives an answer to an interesting and open

problem mentioned by Avramidi in \cite{Avr10}. More importantly, our method is

valid to explicitly calculate all the coefficients $a_l^\mp$, $2\le l\le m$, in

the above asymptotic expansion. As an application, we show that an

$n$-dimensional ball is uniquely determined by its Navier-Lam\'{e} spectrum

among all bounded elastic bodies with smooth boundary.