 Solution

For fixed $p$, the covariant derivative is linear in its tangent-vector input, so $L_p:T_{p}M\to T_{p}M$, $v\mapsto ({\nabla }_{v}X)(p)$, is a well-defined linear map. Its trace defines the divergence of a Riemannian vector field. In the given Riemannian orthonormal frame, write $X={\sum }_j{f}_j{E}_j$. The Levi-Civita connection product rule givesAt $p$, the second sum vanishes by the normal orthonormal frame hypothesis. Thus the $i$th diagonal coefficient of $L_p$ is $E_i(f_i)(p)$, andEquivalently one can take ${\sum }_{i}g({\nabla }_{E_i}X,E_i)$; the trace does not depend on the orthonormal basis chosen.