There is no need to be so specific in the question, you can indeed answer this in general for $M$ a (paracompact, finite-dimensional) manifold.
It is well known, that in this setup $\mathrm{Diff}_c(M)$ is an infinite-dimensional Lie group with Lie algebra $\mathfrak{X}_c(M)$ (this was your space of $C^\infty_c(M,M)$.
Before we continue, it needs to be mentioned that there are at least two settings in which the Lie group statement makes sense and is true:
- The convenient setting of global analysis (a la Kriegl and Michor, this is Theorem 43.1. in 1)
- The Bastiani setting of calculus (again Michor, Theorem 11.11 in 2)
In both settings the exponential is of course the flow map zou mention (this is said explicitely in the convenient calculus source 1, but NOT in the Bastiani source 2.
How does this relate to your question? Well $\exp$ is convenient smooth according to 1, but we can not use that! Convenient smooth mappings are in general NOT continuous with respect to the original topology. Note that this phenomenon happens only outside of Frechet spaces and some other nice classes of spaces. However, the model space of compactly supported vector fields is not one of the nice spaces, so we can not deduce it from this result.
In the Bastiani setting I know several sources wich mention Bastiani smoothness of the exponential map (this is a consequence of regularity of the Lie group, here regularity means that the flow map for curves in the Lie algebra exists and is smooth, in your picture this means that you take the flow map of time dependent vector fields) but only feature proofs in more restricted settings (i.e. $M$ compact). To my knowledge a proof of the regularity of $\mathrm{Diff}(M)$ in the Bastiani setting for paracompact $M$ appeared for the first time in print in my thesis 3 (as a special case of the orbifold diffeomorphism group discussed there).
Regularity of a Lie group (in the Bastiani setting) guarantees smoothness (in Bastiani sense) of the exponential map. Now as Bastiani smooth implies continuity, this answers your first question affirmatively.
First of all the second question seems to be off, or at least contain some non-trivial identification since the vector fields map into the tangent bundle and are still comared using the metric induced by the Riemannian metric on $M$. Anyway, a closer analysis of the problem should show that whatever you want is probably a finite-dimensional problem anyway. Therefore it might be misleading to think about your second property as Lipschitzness of the exponential.
To see this, note that the source of the differentiability of the exponential map is the so called exponential law (for function spaces, has nothing to do with Lie group exponentials). Here this translates roughly to the insight that you do not need to solve the defining equation on an infinite-dimensional manifold, but can place yourself in a situation where you are solving a differential equation on the underlying manifold $M$ (see e.g. 4 for more information on this topic). This is usually highly technical but since your question fixes the compact set, it should be doable (I would have to think more closely about how to do this though).