Pregunta sobre ecuaciones de movimiento "diferentes" en dependencia de índices

Tengamos la acción

$$S = \int (\partial_{\mu}h^{\mu \sigma}\partial^{\nu}h_{\nu \sigma} - \Lambda h^{\mu \nu}T_{\mu \nu}) d^{4}x.$$ Por definición, $$h_{\mu \nu} = h_{\nu \mu} , \quad T_{\mu \nu} = T_{\nu \mu} \neq f(h_{\mu \nu}), \quad h^{\mu \nu} = \eta^{\mu \alpha}\eta^{\nu \beta}h_{\alpha \beta}, \quad \eta^{\mu \nu} = diag (1, -1, -1, -1).$$ Tomemos la variación de $S$por $h_{\mu \nu}$y establecer los términos de la superficie a cero: $$\delta S = \int \partial_{\mu} \delta h^{\mu \sigma} \partial^{\nu}h_{\nu \sigma} d^{4}x + \int \partial_{\mu}h^{\mu \sigma}\partial^{\nu}\delta h_{\nu \sigma} d^{4}x - \Lambda \int \delta h^{\mu \nu} T_{\mu \nu}d^{4}x =$$ $$= -\int (\partial_{\mu}\partial^{\nu}h_{\nu \sigma}\delta h^{\mu \sigma} + \partial_{\mu}\partial^{\nu}h^{\mu \sigma} \delta h_{\nu \sigma} - \Lambda \delta h^{\mu \sigma} T_{\mu \sigma})d^{4}x.$$ Entonces es hora de mi pregunta estúpida. Las expresiones debajo de la integral son de forma escalar y podemos cambiar el nombre de los índices. Pero cuando sacamos la variación $\delta h^{\mu \sigma}$, todo ha cambiado. Puedo hacerlo de las dos "maneras": $$\delta S = \int \delta h^{\mu \sigma}(\partial_{\mu}\partial^{\nu}h_{\nu \sigma} + \partial^{\alpha}\partial_{\mu}h_{\alpha \sigma} - \Lambda T_{\mu \sigma})d^{4}x = \int \delta h^{\mu \sigma}(2\partial_{\mu}\partial^{\nu}h_{\nu \sigma} - \Lambda T_{\mu \sigma})d^{4}x \qquad (1)$$ y $$\delta S = \int \delta h^{\mu \sigma}(\partial_{\mu}\partial^{\nu}h_{\nu \sigma} + \partial^{\alpha}\partial_{\sigma}h_{\alpha \mu} - \Lambda T_{\mu \sigma})d^{4}x . \qquad (2)$$ Poniendo la variación a cero obtengo diferentes ecuaciones de movimiento $(1), (2)$. ¿Dónde cometí el error?