This paper adapts the multivariate optimal control theory to a Riemannian setting. In this sense, a coherent correspondence between the key elements of a standard optimal control problem and several basic geometric ingredients is created, with the purpose of generating a geometric version of Pontryagin’s maximum principle. More precisely, the local coordinates on a Riemannian manifold play the role of evolution variables (“multitime”), the Riemannian structure, and the corresponding Levi–Civita linear connection become state variables, while the control variables are represented by some objects with the properties of the Riemann curvature tensor field. Moreover, the constraints are provided by the second order partial differential equations describing the dynamics of the Riemannian structure. The shift from formal analysis to optimal Riemannian control takes deeply into account the symmetries (or anti-symmetries) these geometric elements or equations rely on. In addition, various submanifold integral cost functionals are considered as controlled payoffs.