Higher derivatives than Jacobi fields

The higher derivatives of the exponential map satisfy the corresponding higher derivative of the Jacobi equation (because the first derivative satisfies the Jacobi equation itself), which is just an inhomogeneous Jacobi equation, where the homogeneous part is the original Jacobi equation, and the inhomogeneous term involves lower order covariant derivatives of the Jacobi field and covariant derivatives of the curvature tensor. So you would proceed recursively, bootstrapping pointwise bounds on lower derivatives, as well as pointwise bounds on the curvature tensor and its covariant derivatives, into a pointwise bound of the derivative of the Jacobi field. You'll need to figure out how get pointwise bounds for a solution to an inhomogeneous self-adjoint linear second order ODE. I'm sure this has been done before, probably for exactly the same purpose as here, but I don't know or remember where.

You might be interested in the Jacobi flow on $TTM$ whose flow lines project to geodesics, velocity fields of geodesics, and Jacobi fields. You can continue to higher order.

• Peter W. Michor: The Jacobi Flow. Rend. Sem. Mat. Univ. Pol. Torino 54, 4 (1996), 365-372 (pdf)