I want to show
$$
(\delta S)Z = \delta (S(\cdot, Z))+\frac{1}{2}\langle S, \mathcal L_Zg\rangle
\tag{1}
$$
where $\delta= -tr_{12}\nabla$ is divergence operator, $\mathcal L$ is Lie derivative, $S\in\Gamma(Sym^2 T^*M)$
What I try: I use abstract index notation. First $$ (\delta S)Z =-(\nabla^a S_{aj})Z = -(\nabla^a S_{aj})Z^j $$ Besides, $$ \delta (S(\cdot,Z)) =-\nabla^a[S(\cdot, Z)]_a $$ since $S(\cdot, Z)= S_{ij}Z^j$, I think $[S(\cdot, Z)]_a=S_{aj}Z^j$. Therefore $$ \delta(S(\cdot, Z)) =-\nabla^a(S_{aj}Z^j)=-(\nabla ^a S_{aj})Z^j - S_{aj}(\nabla^a Z^j) =(\delta S)Z- S_{aj}(\nabla^a Z^j) \tag{2} $$ On the other hand, when $\delta$ is restricted on $\Gamma (Sym^2 T^*M)$, its formal adjoint is $\omega\rightarrow \frac{1}{2}\mathcal L_{\omega^\sharp}g$. Therefore, I have $$ \frac{1}{2}\langle S, \mathcal L_Z g \rangle = \langle \delta(S), Z^\flat \rangle = \langle -\nabla^a S_{aj} , Z^i {g_{ik}} \rangle = -(\nabla^a S_{aj}) Z^j \tag{3} $$ However, I can't get (1) from (2) and (3).
If the answer is complex, a picture of handwritten draft is enough for me. Thanks very much.