Radial Riccati equation for distance spheres ($A^{\prime }+A^2+{\mathcal{R}}_{{\partial }_r}=0$)

Before the cut locus, define $A(X)={\nabla }_X{\partial }_r$ on the tangent space of a distance sphere. This is the negative of the shape operator with outward normal under the convention $S=-\nabla \nu$. A nonsingular matrix $Y$ of radial Jacobi fields satisfies $Y^{\prime \prime }+{\mathcal{R}}_{{\partial }_r}Y=0$ and $A=Y^{\prime }{Y}^{-1}$, yielding the displayed matrix Riccati equation. Taking its trace gives $h^{\prime }+\mathrm{tr}(A^2)+\mathrm{Ric}({\partial }_r,{\partial }_r)=0$, where $h=\mathrm{tr}A={\partial }_r\log \det Y$. This is the local differential identity behind the Bishop-Gromov inequality.