 Solution

Choose an orthonormal basis $(e_1,e_2)$ of $T_{p}S$ and write $e(\theta )=\cos \theta e_1+\sin \theta e_2$. On a punctured disc on which the exponential map is a diffeomorphism, the mapdefines geodesic polar coordinates centred at $p$. The curves $r\mapsto X(r,\theta )$ are unit-speed geodesics, so the first coefficient of the first fundamental form isThe Gauss lemma says that radial and angular coordinate vectors are orthogonal, hence

The angular variation $J=X_{\theta }$ is a Jacobi field along each radial geodesic, with initial dataSince $G=⟨X_{\theta },X_{\theta }⟩=|J{|}^2$, smooth dependence of the geodesic flow givesConsequentlyas required.