We know that the standard result derived in this post
$$\delta \Gamma_{bc}^{a} = \frac{1}{2}g^{da}(\nabla_{b}(\delta g_{dc}) + \nabla_{c}(\delta g_{db}) - \nabla_{d}(\delta g_{bc})) \, .$$
is a tensor. Now, I want to derive how it would look by lowering the index $a$. I am expecting the result to look similar to this
$$g_{ad}\delta \Gamma^{d}_{bc}=\delta \Gamma_{abc} = \frac{1}{2} \left( \nabla_b \delta g_{ac} + \nabla_c \delta g_{ab} - \nabla_a \delta g_{bc} \right)$$
based on the answer presented here. I was expecting it to be a tensor based on the following arguments.
\begin{align} \Gamma_{abc}' &= \frac{\partial x^d}{\partial x'^a}\frac{\partial x^e}{\partial x'^b}\frac{\partial x^f}{\partial x'^c}\Gamma_{def}+\frac{\partial^2 x^d}{\partial x'^b\partial x'^c}\frac{\partial x_d}{\partial x'^a} \end{align} Under variation \begin{align} \Gamma_{abc}'+\delta \Gamma_{abc}'&= \frac{\partial x^d}{\partial x'^a}\frac{\partial x^e}{\partial x'^b}\frac{\partial x^f}{\partial x'^c}(\Gamma_{def}+\delta \Gamma_{def})+\frac{\partial^2 x^d}{\partial x'^b\partial x'^c}\frac{\partial x_d}{\partial x'^a}\\ \delta \Gamma_{abc}' &=\frac{\partial x^d}{\partial x'^a}\frac{\partial x^e}{\partial x'^b}\frac{\partial x^f}{\partial x'^c}\delta \Gamma_{def} \end{align}
However, when I tried to derive $ \delta \Gamma_{abc} $ using the Leibniz rule from $ \delta (g_{ad} \Gamma^d_{bc})$, I got:
$$ \delta (g_{ad} \Gamma^d_{bc}) = (\delta g_{ad}) \Gamma^d_{bc} + g_{ad} (\delta \Gamma^d_{bc}).$$
Substituting $ \delta \Gamma^d_{bc} $,
$$ \delta \Gamma_{abc} = (\delta g_{ad}) \Gamma^d_{bc} + \frac{1}{2} \left( \nabla_b \delta g_{ac} + \nabla_c \delta g_{ab} - \nabla_a \delta g_{bc} \right). $$
There is an extra term $ (\delta g_{ad}) \Gamma^d_{bc} $ from what I was expecting. One possibility that cross my mind is that $\delta \Gamma_{abc}$ is not a tensor at all and the transformation property I used when considering variation is wrong. Could anyone provide clarification as to where I went wrong?