  In this paper, the elastic Dirichlet-to-Neumann map $\Xi_g$ is studied for
the stationary elasticity system in a compact Riemannian manifold $(\Omega,g)$
with smooth boundary $\partial \Omega$. By overcoming methodological
difficulties, we explicitly get matrix-valued full symbol for the elastic
Dirichlet-to-Neumann map $\Xi_g$. We prove that for a strong convex or
extendable real-analytic manifold with boundary, the elastic
Dirichlet-to-Neumann map $\Xi_g$ uniquely determines the metric $g$ of $\Omega$
in the sense of isometry, thereby solving an open problem for the uniqueness of
the metric under real-analytic setting. Furthermore, by calculating the symbol
representation of the resolvent operator $(\Xi-\tau I)^{-1}$ we can explicitly
obtain all coefficients $a_0, a_1 \cdots, a_{n-1}$ of the asymptotic expansion
$\sum_{k=1}^\infty e^{-t \tau_k}\sim \sum_{m=0}^{n-1} a_m t^{m+1-n} +o(1)$ as
$t\to 0^+$, where $\tau_k$ is the $k$-th eigenvalue of the elastic
Dirichlet-to-Neumann map $\Xi_g$ (i.e., $k$-th elastic Steklov eigenvalue).
These coefficients (spectral invariants) provide important geometric
information for the manifold, which give an answer to another open problem for
the elastic Steklov spectral asymptotics.
