The surface gradient of a function defined on a surface $\Gamma \subset \mathbb{R}^n$ is defined $$\nabla_{\Gamma} f = \nabla f - (\nabla f \cdot N)N$$ where $N$ is the unit normal on $\Gamma.$
How do I obtain this from the covariant derivative operator $\nabla_X Y$?