I am currently trying to find an expression for the Lagrangian variation of the Christoffel symbols $\Delta \Gamma^\lambda {}_{\mu\nu}$. For the Eulerian variation $\delta \Gamma^\lambda {}_{\mu\nu}$ everything is a lot simpler since the Eulerian variation $\delta$ commutes with the partial derivative $\partial$; but the Lagrangian variation $\Delta = \delta + \mathcal{L}_\xi$ (where $\xi$ is the Lagrangian displacement, i.e. some arbitrary vector field) does not obey such simple commutation relations.
My hope is that the expression for the Lagrangian variation of the Christoffel symbol is $$ \Delta \Gamma^\mu_{\nu\lambda} = \frac{1}{2} g^{\mu\kappa} \left( \nabla_\lambda \Delta g_{\kappa\nu} + \nabla_\nu \Delta g_{\kappa\lambda} - \nabla_\kappa \Delta g_{\nu\lambda} \right), $$ where $\Delta g_{\mu\nu}$ is the Lagrangian variation of the metric.
I can derive this expression if I could show that $$ \partial_\lambda \mathcal{L}_\xi g_{\mu\nu} = \mathcal{L}_\xi \partial_\lambda g_{\mu\nu}. \qquad (1) $$ But $\partial_\lambda g_{\mu\nu}$ is not a tensorial quantity, so I'm not sure how to calculate the Lie derivative of it. If I pretend that these are the components of a tensor and write down the expression for the Lie derivative (right-hand side of (1)), I end up with the additional terms $(\partial_\lambda \partial_\mu \xi^\sigma) g_{\sigma\nu} + (\partial_\lambda \partial_\nu \xi^\sigma) g_{\mu\sigma}$; I believe they appear only because I wrongly use the tensor formula for $\mathcal{L}_\xi$ when applying it to non-tensor objects.
Is (1) correct and if so, how could this be shown?
A colleague of mine gave me the hint, that is as crucial as it is simple, on how to calculate $$ \mathcal{L}_\xi \partial_\lambda g_{\mu\nu}. $$ Express the partial derivative using the covariant derivative (with the handy fact that the metric is compatible) plus some terms including the Christoffel symbols, i.e. $$ \partial_\lambda g_{\mu\nu} = \nabla_\lambda g_{\mu\nu} + \Gamma^\sigma {}_{\mu\lambda} g_{\sigma\nu} + \Gamma^\sigma {}_{\nu\lambda} g_{\mu\sigma}. $$ After using the Leibniz-rule on the Lie derivative, one also needs an expression for the Lie derivative of Christoffel symbols but this can be found also here on SE-Mathematics, and we have $$ \mathcal{L}_\xi \Gamma^\mu {}_{\nu\lambda} = \xi^\sigma R^\mu {}_{\lambda\sigma\nu} + \nabla_\nu \nabla_\lambda \xi^\mu. $$ Using these relations and the properties of the Riemann tensor, it takes only a few lines of algebraic reformulations to arrive at Eq. (1): $$ \mathcal{L}_\xi \partial_\lambda g_{\mu\nu} = \partial_\lambda \mathcal{L}_\xi g_{\mu\nu}. $$