Recently, I saw the following formula for the non-commutativity of the d'alambert operator acting on the covariant derivative of a scalar field in general relativity, $\Box (\nabla_{\mu}\phi)-\nabla_{\mu}\Box\phi=R_{\mu\nu}\nabla^{\nu}\phi$. How exactly it is derived, considering the metric compatibility and that $\phi$ is a scalar function depending on time?

Recently, I saw the following formula for the non-commutativity of the d'Alembert operator $\Box$ acting on the covariant derivative of a scalar field in general relativity, $\Box (\nabla_{\mu}\phi)-\nabla_{\mu}\Box\phi=R_{\mu\nu}\nabla^{\nu}\phi$. How exactly it is derived, considering the metric compatibility and that $\phi$ is a scalar function depending on time?