Nicolas Burq
email: nicolas.burq@universite-paris-saclay.fr

AND Belhassen Dehman
email: belhassen.dehman@fst.utm.tn

J Ér Ôme

L E Rousseau

Rauch-Taylor / Bardos-Lebeau

Jérôme Le

MEASURE PROPAGATION ALONG C 0 -VECTOR FIELD AND WAVE CONTROLLABILITY ON A ROUGH COMPACT MANIFOLD

Keywords: Mathematics Subject Classification. 35L05, 35Q49, 35Q93, 58Jxx, 93B05, 93B07. 1 1

The celebrated

Measure support propagation: proof of Theorem 1.10 5. Exact controllability: proof of Theorem 1.12 6. Lack of continuity of the control operator with respect to coefficients 6.1. Proof of Theorems 1.14 and 1.14 6.2. Proof of Theorem 1. [START_REF] Gérard | Leichtnam Ergodic properties of eigenfunctions for the Dirichlet problem[END_REF]

Introduction

The observability property for the wave equation has been intensively studied during the last decades mainly because of its deep connection with the problem of exact controllability. Until the end of the 80's, most of the positive results of observability were established under a (global) geometric assumptions, the so-called Γ-condition introduced by J.-L. Lions, essentially based on and well adapted to a multiplier method [START_REF] Lions | Contrôlabilité exacte, Stabilisation et Perturbations de Systèmes Distribués. Tome 1. Contrôlabilité exacte[END_REF]. Later, following [START_REF] Rauch | Exponential decay of solutions to hyperbolic equations in bounded domains[END_REF], Bardos, Lebeau and Rauch established in [START_REF] Bardos | Sharp sufficient conditions for the observation, control and stabilization of waves from the boundary[END_REF], boundary observability inequalities under a geometric control condition (GCC in short), linking the set on which the control acts and the generalized geodesic flow. Proofs of this result are based on microlocal tools, such as the propagation in phase space of wavefront sets in [START_REF] Bardos | Sharp sufficient conditions for the observation, control and stabilization of waves from the boundary[END_REF] or the propagation of microlocal defect measures in more modern proofs [START_REF] Burq | Condition nécessaire et suffisante pour la contrôlabilité exacte des ondes[END_REF]. For the latter approach, microlocal defect measures originate from the concentration phenomena for sequences of waves if one assumes that observability does not hold. Away from boundaries one obtains t H p µ = 0, (1.1) yielding the transport of the measure µ along the bicharacteristic flow in phase space. This flow is generated by the hamiltonian vector field H p associated with the symbol of the wave operator p. However note that despite their high efficiency and robustness, these methods present the great disadvantage of requiring too much regularity in the coefficients of the wave operator and the geometry. To define the generalized bicharacteristic flow and prove the propagation properties mentioned above a minimal smoothness of the metric and the boundary domain is needed. To our knowledge, the best result, C 2 metric, was proven in [START_REF] Burq | Contrôle de l'équation des ondes dans des ouverts peu réguliers[END_REF], and barely misses the natural minimal smoothness required to define the geodesic flow (W 2,∞ ) and thus the geometric control condition.

In this context, in the present article, we address the following natural question: how can one derive observability estimate for the wave equation from optimal observation regions in the case of a nonsmooth metric? This problem has already received some attention and answers by E. Zuazua and his collaborators, in [START_REF] Castro | Concentration and lack of observability of waves in highly heterogeneous media[END_REF][START_REF] Castro | Addendum to: "Concentration and lack of observability of waves in highly heterogeneous media[END_REF], and more recently in [START_REF] Fanelli | Weak observability estimates for 1-D wave equations with rough coefficients[END_REF] (see also the result of [START_REF] Dehman | Observability estimates for the wave equation with rough coefficients[END_REF]). More precisely, in [START_REF] Castro | Concentration and lack of observability of waves in highly heterogeneous media[END_REF][START_REF] Castro | Addendum to: "Concentration and lack of observability of waves in highly heterogeneous media[END_REF], the authors prove a lack of observability of waves in highly heterogeneous media, that is, if the density is of low regularity. In [START_REF] Fanelli | Weak observability estimates for 1-D wave equations with rough coefficients[END_REF], the authors establish observability with coefficients in the Zygmund class and also observability with loss when the coefficients are Log-Zygmund or Log-Lipschitz. Furthermore, this result is proven sharp since one observes a infinite loss of derivatives for a regularity lower than Log-Lipschitz. Note that these analyses are carried out in one space dimension. This calls for the following comments. First, in this simplified framework, for smooth coefficients all the geodesics reach the observability region in uniform time: captive geodesics are not an issue. Second, proofs are based on sidewise energy estimate, a technique that is specific to the one-dimensional setting; the underlying idea consists in exchanging the rôles of the time and space variables and, in fine, in proving hyperbolic energy estimates for waves with rough coefficients. Unfortunately, such method does not extend to higher space dimension. Furthermore, for the low regularity considered in these articles, the geodesic flow is not well defined. Proving propagation results for wavefront sets or microlocal defect measure appears quite out of reach in such cases.

The present work is the first in a series of three articles devoted to the question of observability (and equivalently exact controllability) of wave equation with nonsmooth coefficients. Here, we initiate this study on a compact Riemannian manifold with a rough metric, yet without boundary, while the two forthcoming articles will present the counterpart analysis on manifolds with boundary (or bounded domains of R d ) [START_REF] Burq | Measure and continuous vector field at a boundary I: propagation equation and wave observability[END_REF][START_REF] Burq | Measure and continuous vector field at a boundary II: geodesics and support propagation[END_REF]. The presence of a boundary yields a much more involved analysis and in [START_REF] Burq | Measure and continuous vector field at a boundary I: propagation equation and wave observability[END_REF][START_REF] Burq | Measure and continuous vector field at a boundary II: geodesics and support propagation[END_REF] we develop Melrose-Sjöstrand generalized propagation theory in a low regularity framework. In the present article, our main result is the observability of the wave equation with a C 1 -metric, completed with the stability of the observability property for small Lipschitz (W 1,∞ ) perturbations of the metric. More precisely, we first show that if the geometric control condition in time T holds for geodesics associated with a C 1 -metric g, then the observability property holds for the wave equation, and equivalently exact controllability. For this low regularity case one has to carefully consider the meaning of the geometric condition (or more generally the meaning of a geodesic) since the metric does not define a natural geodesic flow: geodesics are not uniquely defined. Only their existence is guaranteed. Second, we consider a reference C 1 -metric g 0 as above and we prove that observability also holds for any Lipschitz metric g chosen sufficiently close to g 0 (in the Lipschitz topology). It has to be noticed that Lipschitz metric are too rough to permit the use of microlocal tools and a direct proof of the observability property. Even worse for such a metric the geometric control condition itself does not seem make sense (as the generating vector field is only L ∞ ), and we have to use a perturbation argument near the (not so) smooth C 1 reference metric.

Following the strategy of [START_REF] Burq | Contrôle de l'équation des ondes dans des ouverts peu réguliers[END_REF], we argue by contradiction and we prove a propagation result for microlocal defect measures in a low regularity setting. We prove that these measures are solutions to the ODE (1.1) with here H p having C 0 -coefficients. Then, we deduce some general properties about their support. Namelly we show that their support is a union of integral curves of the vector field. This latter step also follows from Ambrosio-Crippa's superposition principle [START_REF] Ambrosio | Continuity equations and ODE flows with non-smooth velocity[END_REF]. Yet, we give a completely different proof which is of interest since it can be extended to the case of a domain with a boundary [START_REF] Burq | Measure and continuous vector field at a boundary I: propagation equation and wave observability[END_REF][START_REF] Burq | Measure and continuous vector field at a boundary II: geodesics and support propagation[END_REF]. We have not been able to extend the approach of [START_REF] Ambrosio | Continuity equations and ODE flows with non-smooth velocity[END_REF] to to that case. To derive the ODE fufilled by the microlocal defect measure we heavily relies on some harmonic analysis results due to R. Coifman and Y. Meyer [START_REF] Coifman | Au delà des opérateurs pseudo-différentiels[END_REF]Proposition IV.7], that expresses that the commutator of a pseudo-differential operator of order one and a Lipschitz function is a bounded operator on L 2 .

Finally, going further in the analysis, we investigate another stability property with respect to perturbations of the metric. We prove that the HUM optimal control associated with a fixed initial data is not stable with respect to perturbations of the metric.

In Section 2 we recall some geometric facts and the notions of pseudo-differential calculus and microlocal defect (density) measures on a manifold. In addition, using bicharacteristics we state the geometric control condition of [START_REF] Bardos | Sharp sufficient conditions for the observation, control and stabilization of waves from the boundary[END_REF] in its classical form (C 2 -metric) and generalized form (C 1 -metric).

In Section 3 we recall what microlocal defect measures are and we show how, if associated with sequences of solutions of PDEs, their support can be estimated and how a transport ODE can be derived, in the particular context of low regularity of coefficients.

Section 4 is devoted to our proof of the support propagation for measures solutions of a ODE with C 0 -coefficients, Theorem 1.10.

In Section 5 we use the results of Sections 3 and the propagation result of Theorem 1.10 to prove the observability and controllability results for the wave equation, Theorems 1.11 and 1.12.

Finally, in Section 6 we prove the results related to stability properties of the HUM control process.

1.2. Setting and well-posedness. Throughout the article, we consider M a d-dimensional C ∞ -compact manifold, that is, a manifold without boundary with a topology that makes it compact equipped with a C ∞ -atlas. We assume that the topology is also given by a Riemannian metric g, to be chosen either Lipschitz or of class C k for some value of k to be made precise below 1 .

We denote by µ g the canonical positive Riemannian density on M, that is, the density measure associated with the density function (det g) 1/2 . We also consider a positive Lipschitz or of class C k -function κ and we define the density κµ g .

The L 2 -inner product and norm are considered with respect to this density κµ g , that is,

(u, v) L 2 (M) = M uv κµ g , u 2 L 2 (M) = M |u| 2 κµ g . (1.2)
We denote by L 2 V (M) the space of L 2 -vector fields on M, equipped with the norm

v 2 L 2 V (M) = M g(v, v) κµ g , v ∈ L 2 V (M).
We recall that the Riemannian gradient and divergence are given by

g(∇ g f, v) = v(f ) and M f div g vµ g = - M v(f ) µ g ,
for f a function and v a vector field, yielding in local coordinates

(∇ g f ) i = 1≤j≤d g ij ∂ x j f, div g v = (det g) -1/2 1≤i≤d ∂ x i (det g) 1/2 v i , with (g ij x ) = (g x,ij ) -1 .

We introduce the elliptic operator

A = A κ,g = κ -1 div g (κ∇ g ), that is, in local coordinates Af = κ -1 (det g) -1/2 1≤i,j≤d ∂ x i κ(det g) 1/2 g ij (x)∂ x j f .
1 Note that despite considering C k metrics with k < ∞ we still impose the underlying manifold to be smooth. This is due to our use of pseudo-differential techniques that are simple to introduce on a smooth manifold. See Section 2.3 Its principal symbol is simply a(x, ξ) = -1≤i,j≤d g ij x ξ i ξ j . Note that for κ = 1, one has A = ∆ g , the Laplace-Beltrami operator associated with g on M. Similarly to ∆ g , the operator A is unbounded on L 2 (M). With the domain D(A) = H 2 (M) one finds that A is self-adjoint, with respect to the L 2 -inner product given in (1.2), and negative. Moreover, one has

(Au, v) L 2 (M) = - M g(∇ g u, ∇ g v) κµ g , u ∈ H 2 (M), v ∈ H 1 (M).
Together with A we consider the wave operator P κ,g = ∂ 2 t -A κ,g + m, with m > 0 a constant and the following equation

P κ,g y = f in (0, +∞) × M, y |t=0 = y 0 , ∂ t y |t=0 = y 1 in M. (1.3) It is well-posed in the energy space H 1 (M) ⊕ L 2 (M).
Proposition 1.1. Consider κ and g both of Lipschitz class. Let (y 0 , y 1 ) ∈ H 1 (M) × L 2 (M) and let f ∈ L 2 0, T ; L 2 (M) , for any T > 0. There exists a unique

y ∈ C 0 [0, +∞); H 1 (M) ∩ C 1 [0, +∞); L 2 (M)
that is a weak solution of (1.3), that is, y |t=0 = y 0 and ∂ t y |t=0 = y 1 and

P κ,g y = f in D (0, +∞) × M .
Remark 1.2. At this level of regularity of κ and g, the well-posedness of the wave equation is classical. For less regular coefficients we refer to [START_REF] Colombini | A note on hyperbolic operators with log-Zygmund coefficients[END_REF] and [START_REF] Colombini | Time-dependent loss of derivatives for hyperbolic operators with non regular coefficients[END_REF].

In what follows, for simplicity we shall consider the case m = 1, that is for

P κ,g = ∂ 2 t -A κ,g + 1.
In this case, we denote by

E κ,g (y)(t) = 1 2 y(t) 2 H 1 (M) + ∂ t y(t) 2 L 2 (M) = 1 2 y(t) 2 L 2 (M) + ∇ g y(t) 2 L 2 V (M) + ∂ t y(t) 2 L 2 (M) ,
the energy of this solution at time t. For a weak solution y of (1.3), if f = 0 this energy is independent of time t, that is,

E κ,g (y)(t) = E κ,g (y)(0) = 1 2 y 0 2 H 1 (M) + y 1 2 L 2 (M) .
Remark 1.3. The equation we consider, with the constant m > 0, is often referred to the Klein-Gordon equation. Here, we keep the name wave equation. We choose this equation instead of the classical wave equation that corresponds to the case m = 0. In fact, on a compact manifold without boundary, constants are eigenfunctions of the elliptic operator A κ,g with 0 as an eigenvalue. Hence, constant functions are solutions to the wave equation and are so-called invisible solutions, as far as the observability property we are interested in is concerned. If one considers a manifold with boundary and say, homogeneous Dirichlet conditions, this issue becomes irrelevant. We could have dealt with the case m = 0 (the usual wave equation) to the price of additional technical complications.

1.3. Exact controllability and observability. Let ω be a nonempty open subset of M and T > 0. The notion of exact controllability for the wave equation from ω at time T is stated as follows.

Definition 1.4 (exact controllability in H 1 (M) ⊕ L 2 (M)). One says that the wave equation is exactly controllable from ω at time T > 0 if for any (y 0 , y 1 ) ∈ H 1 (M) × L 2 (M), there exists f ∈ L 2 ((0, T ) × M) such that the weak solution y to (1.4) P κ,g y = 1 (0,T )×ω f, (y |t=0 , ∂ t y |t=0 ) = (y 0 , y 1 ), as given by Proposition 1.1 satisfies (y, ∂ t y) |t=T = (0, 0). The function f is called the control function or simply the control.

Observability of the wave equation from the open set ω in time T is the following notion.

Definition 1.5 (observability). One says that the wave equation is observable from ω at time T if there exists C obs > 0 such that for any (u 0 , u

1 ) ∈ H 1 (M) × L 2 (M) one has (1.5) E κ,g (u)(0) ≤ C obs 1 (0,T )×ω ∂ t u 2 L 2 (L) , for u ∈ C 0 [0, T ]; H 1 (M) ∩ C 1 [0, T ]; L 2 (M)
the weak solution of solution to P κ,g u = 0 with u |t=0 = u 0 and ∂ t u |t=0 = u 1 as given by Proposition 1.1 (see [START_REF] Lions | Contrôlabilité exacte, Stabilisation et Perturbations de Systèmes Distribués. Tome 1. Contrôlabilité exacte[END_REF]).

Proposition 1.6. Let ω be an open subset of M and T > 0. The wave equation is exactly controllable from ω at time T if and only if it is observable from ω at time T . Remark 1.7. In the case m = 0 the energy function is given by

E κ,g (u)(t) = 1 2 ∂ t u(t) 2 L 2 (M) + ∇ g u(t) 2 L 2 V (M) .
It follows that a constant function u, solution to the the wave equation (∂ 2 t -A)u = 0 has zero energy. Since 1 (0,T )×ω ∂ t u 2 L 2 (L) also vanishes, one sees that such solution are invisible for an observability inequality of the form of (1.5). Possibilities to overcome this difficulty are to work in a quotient space or to change the wave operator into the Klein-Gordon operator. Here, we chose for simplicity the latter option.

1.4. Main results. We introduce the following spaces for the coefficients (κ, g) to distinguish various levels of regularity:

X 2 (M) = {(κ, g); κ ∈ C 2 (M) and g is a C 2 -metric on M}, X 1 (M) = {(κ, g); κ ∈ C 1 (M) and g is a C 1 -metric on M}, Y(M) = {(κ, g); κ ∈ W 1,∞ (M) and g is a W 1,∞ -metric on M}.
We start by recalling the controllability result known for regularity higher than or equal to C 2 , under the Rauch-Taylor geometric control condition. Definition 1.8 (Rauch-Taylor, geometric control condition). Let g be a C k metric, k = 1 or 2, and let ω be an open set of M and T > 0. One says that (ω, T ) fulfills the geometric control condition if all maximal geodesics associated with g, travelled at speed one, encounter ω for some time t ∈ (0, T ).

A second formulation of this geometric condition based on the dual notion of bicharacteristics is given in Section 2.2 below. Theorem 1.9 (Exact controllability -C 2 -regularity). Consider (κ, g) ∈ X 2 (M), ω an open subset of M and T > 0 such that (ω, T ) fulfills the geometric control condition of Definition 1.8. Then, the wave equation is exactly controllable from ω at time T .

This result was first proven by Rauch and Taylor [START_REF] Rauch | Exponential decay of solutions to hyperbolic equations in bounded domains[END_REF] for a smooth metric. The case (κ, g) ∈ X 2 (M) was proven by the first author in [START_REF] Burq | Contrôle de l'équation des ondes dans des ouverts peu réguliers[END_REF]. On smooth open sets of R d , or equivalently on manifolds with boundary equipped with smooth (κ, g), for instance in the case of homogeneous Dirichlet boundary conditions, this result is given in the celebrated articles of Bardos, Lebeau and Rauch [START_REF] Bardos | Un exemple d'utilisation des notions de propagation pour le contrôle et la stabilisation de problèmes hyperboliques[END_REF][START_REF] Bardos | Sharp sufficient conditions for the observation, control and stabilization of waves from the boundary[END_REF].

In the present article, we extend the result of Theorem 1.9 to cases of rougher coefficients. Our extension is twofold: (1) we treat the case (κ, g) ∈ X 1 (M) and, (2) we treat small perturbations in Y(M) of some (κ, g) ∈ X 1 (M). Most importantly, these two results rely on the understanding of the structure of the support of a nonegative measure subject to a homogeneous transport equation with continuous coefficients. Consider a continuous vector field X on O and let µ be a nonnegative measure density on O. Assume that µ is such that t Xµ = 0 in the sense of distributions, that is,

(1.6) t Xµ, a 1 D (O),C ∞ c (O) = µ, Xa 1 D ,0 (O),C 0 c (O) = 0, a ∈ C ∞ c (O).
If X is moreover Lipschitz, one concludes that µ is invariant along the flow that X generates. However, if X is not Lipschitz, there is no such flow in general. Yet, integral curves do exist by the Cauchy-Peano theorem. The following theorem provides a structure of the support of µ.

Theorem 1.10. Let X be a continuous vector field on O and µ be a nonnegative density measure on O that is solution to t Xµ = 0 in the sense of distributions. Then, the support of µ is a union of maximally extended integral curves of the vector field X.

In other words, if m 0 ∈ O is in supp(µ), then there exist an interval I in R with 0 ∈ I and a C 1 curve γ : I → O that cannot be extended such that γ(0) = m 0 and d ds γ(s) = X(γ(s)), s ∈ I, and γ(I) ⊂ supp(µ). Theorem 1.10 can actually be obtained as a consequence of the superposition principle of L. Ambrosio and G. Crippa [START_REF] Ambrosio | Continuity equations and ODE flows with non-smooth velocity[END_REF]Theorem 3.4]. Here, we provide an alternative proof that is of interest as it allows one to extend this measure support structure result to the case of an open set or a manifold with boundary [START_REF] Burq | Measure and continuous vector field at a boundary II: geodesics and support propagation[END_REF] as needed for our application to observability and controllability. Ambrosio and Crippa's proof is based on a smoothing-by-convolution argument. Extending this approach does not seem to be straightforward in the context of a boundary.

Theorem 1.10 is proven in Section 4 and its proof is independent of the other sections of the article. A reader only interested in our proof of Theorem 1.10 may thus head to Section 4 directly.

1.4.2. Exact controllability results. If (κ, g) ∈ X 2 (M), x ∈ M and v ∈ T x M there is a unique geodesic originating from x in direction v. In the case (κ, g) ∈ X 1 (M) uniqueness is lost. Existence holds however and maximal (here global, see below) geodesics can still be defined by the Cauchy-Peano theorem. In particular, the geometric control condition of Definition 1.8 still makes sense. As announced above, our first result is the following theorem.

Theorem 1.11 (Exact controllability -C 1 -regularity). Consider (κ, g) ∈ X 1 (M), ω an open subset of M and T > 0 such that (ω, T ) fulfills the geometric control condition of Definition 1.8. Then, the wave equation is exactly controllable from ω at time T .

A second result is the following perturbation result.

Theorem 1.12 (Exact controllability -Lipschitz perturbation). Let (κ 0 , g 0 ) ∈ X 1 (M), ω an open subset of M and T > 0 be such that (ω, T ) fulfills the geometric control condition of Definition 1.8 with respect to the metric g 0 . There exists ε > 0 such that for any (κ, g) ∈ Y(M) satisfying

(κ, g) -(κ 0 , g 0 ) Y(M) ≤ ε,
the wave equation associated with (κ, g) is exactly controllable by ω in time T .

Observe that Theorem 1.11 is a direct consequence of Theorem 1.12. We shall thus concentrate on this second more general result. Its proof relies on the measure support structure result of Theorem 1.10.

The sequence of Theorems 1.9, 1.11, and 1.12 calls for the following important comment. Under the assumption of Theorem 1.9, that is, (κ, g) ∈ X 2 (M), there is a geodesic flow and the geometric condition of Definition 1.8 is actually a condition on the flow. Under the assumption of Theorem 1.11, that is, (κ, g) ∈ X 1 (M), as pointed out above there is no geodesic flow in general. Yet, maximal geodesics are still well defined and, the geometric condition of Definition 1.8 makes sense because it does not refer to a flow. However, under the assumption of Theorem 1.12, that is, (κ, g) ∈ Y(M), geodesics cannot be defined in general. No geometric condition can be formulated. Yet, Theorem 1.12 is a perturbation result and a geometric condition is expressed for a reference pair (κ 0 , g 0 ) ∈ X 1 (M) around which a (small) neighborhood in Y(M) is considered.

The following remark further emphasizes that the perturbation is to be considered around a pair (κ 0 , g 0 ) ∈ X 1 (M) for which the geometric control condition holds and not around a pair (κ 0 , g 0 ) ∈ X 1 (M) for which exact controllability (or equivalently observability) holds.

Remark 1.13 (On the perturbation result). Having both our results, geometric control for C 1 metrics and Lipschitz stability of exact controllability around a reference metric satisfying the geometric control condition, a natural question is whether the exact controllability property is itself stable by perturbation. On the one hand, it is classical that the exact controllability property is stable under lower-order perturbations of the elliptic operator A κ,g but on the other hand, it is possible to show that it is not stable under (smooth) perturbations of the geometry or the metric.

Let us illustrate this instability property with a quite simple example. Consider the wave equation on the sphere

S d = {x ∈ R d+1 ; i x 2 i = 1},
endowed with its standard metric and with control domain the open hemisphere

ω = {x ∈ S d ; x 1 > 0}.
Even though ω does not fulfill the geometric control condition of Definition 1.8 exact controllability holds for this geometry, an unpublished result by G. Lebeau (see [25, Section VI.B] and [START_REF] Zhu | Stabilization of damped waves on spheres and Zoll surfaces of revolution[END_REF] for extensions). Consider now the sphere endowed with the above standard metric, with the smaller control domain

ω ε = {x ∈ S d ; x 1 > ε},
for some ε > 0. This second geometry is ε-close to the Lebeau example in the C ∞ -topology.

Yet, for all ε > 0, exact controllability does not hold, because there exists a geodesic (the equator, {x ∈ S d ; x 1 = 0}) that does not encounter ω ε ). This shows that in Theorem 1.12, the assumption that the reference geometry should satisfy the geometric control condition cannot be replaced by the weaker assumption that it should satisfy the exact controllability property. This also shows that our perturbation argument will have to be performed on the actual proof that geometric control implies exact controllability and not on the final property itself.

1.4.3. Further results on the control operator. We finish this section with results analyzing the influence of some metric perturbations on the control process.

We introduce further levels of regularity for the coefficients by setting for k ∈ N ∪ {+∞},

X k (M) = {(κ, g); κ ∈ C k (M) and g is a C k -metric on M}.
First, we consider k ≥ 2. We recall the notation

P κ,g = ∂ 2 t -A κ,g +1 with A κ,g = κ -1 div g (κ∇ g
), and we assume that (κ, g) ∈ X k (M), and that (ω, T ) satisfies the geometric control condition of Definition 1.8 for geodesics given by the metric g. Then, by Theorem 1.9, given (y 0 , y 1 ) ∈ H 1 (M) × L 2 (M), there exists f ∈ L 2 ((0, T ) × ω) such that the solution to (1.4) satisfies y(T ) = 0 and ∂ t y(T ) = 0. One can prove that among all possible control functions there is one of minimal L 2 -norm. We denote by f y 0 ,y 1 κ,g this control function usualy named HUM control function (cf. for instance [START_REF] Lions | Contrôlabilité exacte, Stabilisation et Perturbations de Systèmes Distribués. Tome 1. Contrôlabilité exacte[END_REF]). Moreover, the map

H κ,g : H 1 (M) ⊕ L 2 (M) → L 2 ((0, T ) × M) (1.7) (y 0 , y 1 ) → f y 0 ,y 1 κ,g , is continuous. Note that f y 0 ,y 1 κ,g
is actually a weak solution of the wave equation with initial data in L 2 (M) × H -1 (M), meaning that one moreover has f y 0 ,y 1 κ,g

∈ C 0 ([0, T ], L 2 (M)).
Theorem 1.14 (Lack of continuity of the HUM-operator -Case k ≥ 2). Let k ≥ 2 and (κ, g) as above. For any neighborhood U of (κ, g) in X k (M), there exist (κ, g) ∈ U and an initial data (y 0 , y

1 ) ∈ H 1 (M) × L 2 (M), with y 0 2 H 1 + y 1 2 L 2 = 1
, such that the respective solutions y and ỹ of (1.8)

P κ,g y = 1 (0,T )×ω f y 0 ,y 1 κ,g in (0, T ) × M, (y, ∂ t y) |t=0 = (y 0 , y 1 ) in M, P κ,g ỹ = 1 (0,T )×ω f y 0 ,y 1 κ,g in (0, T ) × M, (ỹ, ∂ t ỹ) |t=0 = (y 0 , y 1 ) in M, are such that (1.9) E κ,g (ỹ -y)(T ) = E κ,g (ỹ)(T ) ≥ 1/2.
Moreover, there exists C T > 0 such that

(1.10) (H κ,g -H κ,g )(y 0 , y 1 ) L 2 ((0,T )×ω) = f y 0 ,y 1 κ,g -f y 0 ,y 1 κ,g
L 2 ((0,T )×ω) ≥ C T , for (y 0 , y 1 ) as given above.

Remark 1.15. The result of Theorem 1.14 states that starting from the same initial data and solving the two wave equations with the same control vector f κ,g associated with P κ,g , a small perturbation of the metric can induce a large error for the final state (y(T ), ∂ t y(T )). In other words, the two dynamics are no longer close. In particular, the map

X k (M) (κ, g) -→ H κ,g ∈ L H 1 (M) ⊕ L 2 (M), L 2 ((0, T ) × M)
is not continuous.

Remark 1.16. The result of Theorem 1.14 can also be stated on open bounded smooth domains of R n in the case of homogeneous Dirichlet condition. In fact, as can be checked in what follows, its proof only relies on basic properties of microlocal defect measures (support localization and propagation) that are known to be valid in this framework (see [START_REF] Burq | Contrôle de l'équation des ondes dans des ouverts peu réguliers[END_REF]).

Remark 1.17. In the statement of Theorem 1.14 if the neighborhood U of (κ, g) in X k is small enough, the pair (ω, T ) also satisfies the geometric control condition of Definition 1.8 for (κ, g) and therefore f y 0 ,y 1 κ,g is well defined. In particular, this is clear as in the case k ≥ 2 there is a well defined and unique geodesic flow.

The case k = 1 is quite different as there is no geodesic flow, as already mentioned above. However, given (κ, g) ∈ X 1 and (ω, T ) if the Rauch-Taylor geometric control condition of Definition 1.8 holds for (ω, T ) for the geodesics associated with g, given any neighborhood U of (κ, g) in X 1 one can still find (κ, g) ∈ U such that (1) the geometric control condition still holds for the geodesics associated with g, (2) the result of Theorem 1.14 also holds.

Theorem 1.14 (Lack of continuity of the HUM-operator -Case k = 1). Let k = 1 and (κ, g) ∈ X 1 as above. For any neighborhood U of (κ, g) in X 1 (M), there exist (κ, g) ∈ U and an initial data (y 0 , y

1 ) ∈ H 1 (M) × L 2 (M), with y 0 2 H 1 + y 1 2 L 2 = 1
, such that the geometric control condition of Definition 1.8 for geodesics given by the metric g holds and moreover the results listed in Theorem 1.14 hold.

The proofs of Theorems 1.14 and 1.14 are given in Section 6.1.

We finish this section with some remarks and some questions.

Remark 1.18. In all results above we have used 1 (0,T )×ω as a control operator, that is, the characteristic function of an open set. We could have also considered a control operator given by 1 (0,T ) (t)χ(x), with χ a smooth function on M. The controlled wave equation then has the form

(1.11) P κ,g y = 1 (0,T ) χ f, (y |t=0 , ∂ t y |t=0 ) = (y 0 , y 1 ),
In such case, the open set to be used in the geometric control condition is ω = {χ = 0}. This is often done this way, in particular since the smoothness of the function χ allows one to use some microlocal techniques that require regularity in the operator coefficients. The results and proofs of the present article can be written mutatis mutandis for this type of control operator.

1.4.4. Comparison with the smooth case and some open questions. Following on the previous remark, with a smooth in space control operator, one can wonder above the smoothness of the HUM operator. This question is addressed in the work of the second author jointly with G. Lebeau [START_REF] Dehman | Analysis of the HUM control operator and exact controllability for semilinear waves in uniform time[END_REF]. In fact, a gain of regularity in the initial data (y 0 , y 1 ) yields an equivalent gain of regularity in the HUM control function f y 0 ,y 1 κ,g . For instance, for (y 0 , y 1 ) ∈ H 2 (M)×H 1 (M) one finds f y 0 ,y 1 κ,g ∈ C 0 ([0, T ], H 1 (M)). Note that the result of [START_REF] Dehman | Analysis of the HUM control operator and exact controllability for semilinear waves in uniform time[END_REF] is proven in the case of smooth coefficients, that is, (κ, g) ∈ X ∞ . We thus consider this smooth case in the discussion that ends this introductory section. Open questions around the results of Theorems 1.14 and Theorem 1.14 are then raised.

As we shall see in their proof, the result of Theorems 1.14 and Theorem 1.14 relies on the high frequency behavior of the solutions to (1.8). In the case of smooth coefficients and a smooth control operator, if we assume smoother data (y 0 , y 1 ) in the HUM control process, the result of Theorem 1.14 does not hold any more. The HUM control process becomes regular with respect to (κ, g) as expressed in the following proposition.

Proposition 1.19 (HUM control process for smooth data). Consider (κ, g) ∈ X ∞ (M) and let χ ∈ C ∞ (M). Set ω = {χ = 0} and assume that (ω, T ) fulfills the geometric control condition of Definition 1.8 for the geodesics associated with (κ, g). Let α ∈ (0, 1]. There exists C α > 0 such that for any (κ, g) ∈ X ∞ (M) and any (y 0 , y 1 ) ∈ H 1+α (M) × H α (M), the respective solutions y and ỹ to

P κ,g y = 1 (0,T ) χ f y 0 ,y 1 κ,g in (0, T ) × M, (y, ∂ t y) |t=0 = (y 0 , y 1 ) in M, P κ,g ỹ = 1 (0,T ) χ f y 0 ,y 1 κ,g in (0, T ) × M, (ỹ, ∂ t ỹ) |t=0 = (y 0 , y 1 ) in M, satisfy E κ,g (y -ỹ)(T ) 1/2 ≤ C α (κ, g) -(κ, g) α X 1 (M) (y 0 , y 1 ) H 1+α (M)⊕H α (M) .
The proof of Proposition 1.19 is given in Section 6.2.

In the above proposition coefficients are chosen smooth, quite in contrast with the rest of this article. As explained above, and as the reader can check in the proof, this lies in the use of the regularity of the HUM operator with respect to the data (y 0 , y 1 ), a result proven for smooth coefficients in [START_REF] Dehman | Analysis of the HUM control operator and exact controllability for semilinear waves in uniform time[END_REF]. The result of Proposition 1.19 raises the following natural questions:

(1) Does the HUM operator exhibit regularity with respect to the data (y 0 , y 1 ) similar to what is proven in [START_REF] Dehman | Analysis of the HUM control operator and exact controllability for semilinear waves in uniform time[END_REF] in the case of not so smooth coefficients? (2) If so, if one increases the smoothness of the data (y 0 , y 1 ) as in Proposition 1.19, does the HUM control process also become regular with respect of the metric?

Geometric aspects and operators

We define the smooth manifold L = R × M and T * L its cotangent bundle. We denote by π : T * L → L the natural projection. Elements in T * L are denoted by (t, x, τ, ξ). One has π(t, x, τ, ξ) = (t, x).

Setting |ξ| 2 x = g x (ξ, ξ) the Riemannian norm in the cotangent space of M at x, we define

S * L = {(t, x, τ, ξ) ∈ T * L, τ 2 + |ξ| 2 x = 1},
the cosphere bundle of L. We shall also use the associated cosphere bundle in the spatial variables only,

S * M = {(x, ξ) ∈ T * M, |ξ| 2 x = 1/2}. For a C k -metric both S * M and S * L are C k -manifolds. Consider a C ∞ -atlas A M = (C M j ) j∈J of M, #J < ∞, with C M j = (O j , θ j )
where O j is an open set of M and θ j : O j → Õj is a bijection for Õj an open set of R d . For j ∈ J, we define

C j = (O j , ϑ j ) with O j = R × O j and ϑ j : O j → Õj (t, x) → t, θ j (x) , with Õj = R × Õj . Then A = (C j ) j∈J is a C ∞ -atlas for L.
In what follows for simplicity we shall use the same notation for an element of T * L and its local representative if no confusion arises.

2.1. Hamiltonian vector field and bicharacteristics. Let (κ, g) ∈ X k , k = 1 or 2. The principal symbol of the wave operator P κ,g is given by

(2.1) p(t, x, τ, ξ) = p κ,g (t, x, τ, ξ) = -τ 2 + |ξ| 2 x , (t, x, τ, ξ) ∈ T * L.
In local charts, one has

p(t, x, τ, ξ) = -τ 2 + 1≤i,j≤d g ij (x)ξ i ξ j .
Note that (g ij (x)) i,j is the inverse of (g ij (x)) i,j , the latter being the local representative of the metric.

We denote by H p the Hamiltonian vector field associated with p, that is, the unique vector field such that {p, f } = H p f for any smooth function f . Here, {., .} denotes the Poisson bracket, that is, in local chart

{p, f } = ∂ τ p ∂ t f -∂ t p ∂ τ f + 1≤j≤d (∂ ξ j p ∂ x j f -∂ x j p ∂ ξ j f ), (2.2) yielding H p = -2τ ∂ t + ∇ ξ p • ∇ x -∇ x p • ∇ ξ ,
as p is in fact independent of the time variable t. The Hamiltonian vector field H p is of class C k-1 . Observe that, for a function f of the variables (t, x, τ, ξ), one has

t H p f = 2τ ∂ t f -div x (f ∇ ξ p) + div ξ (f ∇ x p),
with which one deduces

t H p = -H p , (2.3) even in the case (κ, g) ∈ X 1 .
First, consider the case k = 2. Thus, H p is a C 1 -vector field. For ∈ T * L one denotes by s → φ s ( ) the unique maximal solution to (2.4) d ds φ s ( ) = H p φ s ( ), s ∈ R, and φ s=0 ( ) = , as given by the Cauchy-Lipschitz theorem. One calls (s, ) → φ s ( ) the Hamiltonian flow map. Let s → γ(s) be an integral curve of H p , that is, γ(s) = φ s ( ) for some ∈ T * L. For any smooth function f on T * L one has

d ds f • γ(s) = H p f γ(s) ,
Note that H p τ = 0, meaning that the variable τ is constant along γ. Note also that the value of p remains constant along γ since H p p = {p, p} = 0. Hence, |ξ| 2 x = g x (ξ, ξ) is also constant. Thus, if γ(0) ∈ S * L then γ(s) remains in S * L, and for ∈ S * L, the vector field H p at is tangent to S * L. consequently, we may consider H p as a tangent vector field on the

C 2 -manifold S * L. In particular H p a makes sense if a ∈ C 1 c (S * L). If moreover a ∈ C 2+ c (S * L), ≥ 0, one has H p a ∈ C 1 c (S * L)
. Since H p p = 0, the flow φ s preserves Char(p) = p -1 ({0}), the characteristic set of p. As is done classically, we call bicharacteristic an integral curve for which p = 0. Observe then that (2.4) defines a flow on the C 2 -manifold

Char(p) ∩ S * L = {(t, x, τ, ξ); τ 2 = 1/2 and |ξ| 2 x = 1/2}.
Second, consider the case k = 1. Then H p is only a continuous vector field. Thus, for any ∈ Char(p) there exists a maximal bicharacteristic s → γ(s) defined on R such that γ(0) = , that is,

d ds γ(s) = H p γ(s) , s ∈ R,
by the Cauchy-Peano theorem. Uniqueness is however not guaranteed and the notion of flow cannot be used in the case k = 1. Since the value of |ξ| x remains constant and the manifold M is compact, maximal bicharacteristics are actually defined globally.

As above, if γ(0) ∈ S * L (resp. Char(p) ∩ S * L) one has γ(s) ∈ S * L (resp. Char(p) ∩ S * L) for all s ∈ R. The hamiltonian vector field H p can be viewed a C 0 -vector field on the C 1manifold S * L (resp. on the

C 1 -manifold Char(p) ∩ S * L). For a ∈ C 1+ c (S * L), ≥ 0, one finds H p a ∈ C 0 c (S * L).
Finally, connection between bicharacteristic and geodesics can be made. For this we recall that if ξ ∈ T *

x M for some x ∈ M one can define v ∈ T x M by v = ξ , that reads in local coordinates v i = j g ij (x)ξ j . In particular |v| 2

x = g x (v, v) = |ξ| x . If now 0 = (t 0 , x 0 , τ 0 , ξ 0 ) ∈ Char(p)∩S * L and let s → (s) = t(s), x(s), τ, ξ(s) be a bicharacteristic such that (0) = 0 . One has τ = τ 0 and t(s) = t 0 -2τ 0 s. The map

X : t → x (t 0 -t)/(2τ 0 ) ,
can be proven to be the geodesic originating from x 0 in the direction given by v 0 = (ξ 0 ) and parameterized by t.

We now compute the speed at which the geodesic is travelled. We have dX dt (t) = -1

2τ 0 dx(s)
ds , which yields

dX dt (t) = -1 2τ 0 ∇ ξ p x(s), ξ(s) = - ξ(s) τ 0 .
It follows that

| dX dt (t)| x = |ξ(s) | x /|τ 0 | = |ξ(s)| x /|τ 0 | = |ξ 0 | x /|τ 0 | = 1,
since 0 ∈ Char(p). Hence, the projection of the bicharacteristic s → γ(s) yields a geodesic travelled at speed one.

2.2. Geometric control condition. As the projections of bicharacteristics onto L yield geodesics, in the case k ≥ 2, we can state the Rauch-Taylor geometric control condition [START_REF] Rauch | Exponential decay of solutions to hyperbolic equations in bounded domains[END_REF] formulated in Definition 1.8 with the notion of Hamiltonian flow introduced above.

Definition 1.8 (geometric control condition, k ≥ 2). Let g be a C 2 metric and let ω be an open set of M and T > 0. One says that (ω, T ) fulfills the geometric control condition if for all ∈ Char(p) one has π φ s ( ) ∈ (0, T ) × ω for some s ∈ R.

In the case k = 1, since g is only C 1 there is no flow in general, one rather writes the geometric control condition by means of maximal bicharacteristics. Definition 1.8 (generalized geometric control condition, k = 1). Let g be a C 1 metric and let ω be an open set of M and T > 0. One says that (ω, T ) fulfills the geometric control condition if for any maximal bicharacteristic s → γ(s) in Char(p) one has π γ(s) ∈ (0, T ) × ω for some s ∈ R.

In other words, for all ∈ Char(p), all bicharacteristics that go through meet the cotangent bundle above (0, T ) × ω.

Naturally, the Definitions 1.8 and 1.8 coincide in the case k = 2 because of the uniqueness of a bicharacteristic going through a point of Char(p).

Symbols and pseudo-differential operators.

Here, we follow [5, Section 1.1] for the notation. We denote by H k (X) or H k loc (X), with X = M or L, the usual Sobolev space for complex valued functions, endowed with its natural inner product and norm. In particular, the L 2 (X)-inner product is denoted by (., .) L 2 (X) .

Classical polyhomogeous symbol classes on T * R n R n × R n are denoted by S m ph (R n × R n ) and the classes of associated operators by Ψ m ph (R n ). We recall that symbols in the class S m ph (R n ×R n ) behave well with respect to changes of variables, up to symbols in S m-1 ph (R n ×R n ) (see [START_REF] Hörmander | The analysis of linear partial differential operators[END_REF]Theorem 18.1.17 and Lemma 18.1.18]).

We define S m c,ph (T * L) as the set of polyhomogeneous symbols of order m on T * L with compact support in the variables (t, x) ∈ L (note that compactness with respect to x ∈ M is obvious). Having the manifold M smooth is important for symbols and following pseudodifferential operators to be simply defined.

For any m, the restriction to the sphere

(2.5) S m c,ph (T * L) → C ∞ c (S * L), a → a |S * L , is
onto. This allows one to identify a homogeneous symbol with a smooth function on S * L with compact support.

We denote by Ψ m c,ph (L) the space of polyhomogeneous pseudo-differential operators of order m on L:

one says that Q ∈ Ψ m c,ph (L) if Q maps C ∞ c (L) into D (L) and (1) its kernel K(x, y) ∈ D (L × L) is such that supp(K) is compact in L × L;
(2) K(x, y) is smooth away from the diagonal ∆ L = {(t, x; t, x); (t, x) ∈ L};

(3) for any local chart C j = (O j , ϑ j ) and all φ [START_REF] Hörmander | The analysis of linear partial differential operators[END_REF]Chapter 18.1]). Note that the principal symbol is uniquely defined in S m c,ph (T * L) because of the polyhomogeneous structure (see the remark following Definition 18.1.20 in [START_REF] Hörmander | The analysis of linear partial differential operators[END_REF]). The application σ m enjoys the following properties.

0 , φ 1 ∈ C ∞ c ( Õj ) one has φ 1 • ϑ -1 j * • Q • ϑ * j • φ 0 ∈ Op S m c,ph (R d+1 × R d+1 ) . For Q ∈ Ψ m c,ph (L), we denote by σ m (Q) ∈ S m c,ph (T * L) the principal symbol of Q (see
(1) The map

σ m : Ψ m c,ph (L) → S m c,ph (T * L) is onto. (2) For all Q ∈ Ψ m c,ph (L), σ m (Q) = 0 if and only if Q ∈ Ψ m-1 c,ph (L). (3) For all Q ∈ Ψ m c,ph (L), σ m (Q * ) = σ m (Q). (4) For all Q 1 ∈ Ψ m 1 c,ph (L) and Q 2 ∈ Ψ m 2 c,ph (L), one has Q 1 Q 2 ∈ Ψ m 1 +m 2 c,ph (L) with σ m 1 +m 2 (Q 1 Q 2 ) = σ m 1 (Q 1 )σ m 2 (Q 2 ).
(

) For all Q 1 ∈ Ψ m 1 c,ph (L) and Q 2 ∈ Ψ m 2 c,ph (L), one has [Q 1 , Q 2 ] = Q 1 Q 2 -Q 2 Q 1 ∈ Ψ m 1 +m 2 -1 c,ph (L) with σ m 1 +m 2 -1 ([Q 1 , Q 2 ]) = 1 i {σ m 1 (Q 1 ), σ m 2 (Q 2 )}. (6) If Q ∈ Ψ m c,ph (L), then Q maps continuously H k loc (L) into H k-m comp (L). In particular, for m < 0, Q is compact on L 2 loc (L). Given an operator Q ∈ Ψ m c,ph (L), one sets Char(Q) = Char σ m (Q) = { ∈ T * L, σ m (Q)( ) = 0}. 5

Microlocal defect measure and propagation properties

A defect measure is used to characterize locally the failure of a sequence to strongly converge, meaning some concentration phenomenum. This characterization can be made finer by further considering microlocal concentration phenomena. for the duality bracket. This notation will also be used for a ∈ S 0 c,ph (T * L) according to the identification map (2.5).

Consider a sequence (u k ) k∈N ⊂ L 2 loc (L) that converges weakly to 0. Here, to define the L 2 -norm and inner product on L we use a fixed (κ 0 , g 0 ) chosen in X 1 (M); see (1.2).

As a consequence of [18, Theorem 1], there exists a subsequence of (u k ) k∈N (still denoted by (u k ) k∈N in what follows) and a density measure µ ∈ M + (S * L), such that

(3.1) lim k→∞ Qu k , u k L 2 comp (L),L 2 loc (L) = µ, σ 0 (Q) S * L , for any Q ∈ Ψ 0 c,ph (L).
Recall that symbols in S 0 c,ph (T * L) are compactly supported in time t here. We also refer to [START_REF] Tartar | H-measures, a new approach for studying homogenisation, oscillations and concentration effects in partial differential equations[END_REF] and [START_REF] Burq | Mesures semi-classiques et mesures de défaut Séminaire Bourbaki[END_REF]. One calls µ a microlocal defect (density) measure associated with (u k ) k∈N .

Similarly, one can use the notion of H 1 -microlocal defect density measure. Consider (u k ) k∈N ⊂ H 1 loc (L) that converges weakly to 0. Then, there exists a subsequence of (u k ) k∈N (still denoted by (u k ) k∈N ) and a density measure µ ∈ M + (S * L) such that for any

Q ∈ Ψ 2 c,ph (L) (3.2) lim k→∞ Qu k , u k H -1 comp (L),H 1 loc (L) = µ, σ 2 (Q) S * L .
Naturally, in either cases, the density measure µ depends on the choice made of (κ 0 , g 0 ) ∈ X 1 (M). In what follows we shall make clear what choice is made.

3.2.

Local representatives. Consider a finite atlas A = (C j ) j∈J on L, as introduced in Section 2, with C j = (O j , ϑ j ). Consider a smooth partition of unity (χ j ) j∈J subordinated to the covering by the open sets (O j ) j . We consider also χj , χj ∈ C ∞ (L) supported in O j such that χj ≡ 1 on a neighborhood of supp(χ j ) and χj ≡ 1 on a neighborhood of supp( χj ). Set also χ C j j = (ϑ -1 j ) * χ j , χC j j = (ϑ -1 j ) * χj , and χC j j = (ϑ -1 j ) * χj . One has χ

C j j , χC j j , χC j j ∈ C ∞ c ( Õj ), with Õj = ϑ j (O j ). Let (u k ) k ⊂ H 1 loc (L) that converges weakly to 0, Q ∈ Ψ 2 c,ph (
L), and j ∈ J. One can write

χ j Q = χ j Q χj + χ j Q(1 -χj ). Since χ j Q(1 -χj ) is a regularizing operator one finds µ, χ j σ 2 (Q) S * L ∼ χ j Qu k , u k H -1 comp (L),H 1 loc (L) ∼ χ j χj Q χj v k j , v k j H -1 comp (L),H 1 loc (L) ,
as k → +∞, for v k j = χj u k . The operator Q j = (ϑ -1 j ) * χj Q χj (ϑ j ) * is a pseudo-differential operator of order 2 on R d+1 with principal symbol q j = χ2 j q C j , where q C j is the local representative of σ 2 (Q). Set also v k,C j j = (ϑ -1 j ) * v k j . It converges weakly to 0 in H 1 (R d+1 ). Associated with this sequence is a microlocal defect measure µ j . If one writes

χ j χj Q χj v k j , v k j H -1 comp (L),H 1 loc (L) = χ C j j Q j v k,C j j , v k,C j j H -1 comp (R d+1 ),H 1 loc (R d+1 ) , one obtains µ, χ j σ 2 (Q) S * L = µ j , χ C j j q j S * Õj = µ j , χ C j j q C j S * Õj .
Note that here, the L 2 and H s -norms on R d+1 are based on the local representative of the density measure κ 0 µ g 0 dt. One thus sees that the local representative of χ j µ is precisely χ C j j µ j , that is, χ j µ = ϑ * j χ C j j µ j = χ j ϑ * j µ j . Summing up, we thus have [START_REF] Coifman | Au delà des opérateurs pseudo-différentiels[END_REF]Proposition IV.7]) and some of its consequences that we list below.

µ = j∈J χ j µ = j∈J χ j ϑ * j µ j . and µ, σ 2 (Q) S * L = j∈J µ, χ j σ 2 (Q) S * L = j∈J µ j , χ C j j q C j S *
Theorem 3.2 (Coifman-Meyer). Let Q ∈ Ψ 1 ph (R n × R n ). If m ∈ W 1,∞ (R n ) the commutator [Q, m] maps L 2 (R n ) into itself continuously. Moreover there exists C > 0 such that [Q, m] L 2 →L 2 ≤ C m W 1,∞ , m ∈ W 1,∞ (R n ).
We deduce the following corollary.

Corollary 3.3. Let Q ∈ Ψ 1 ph (R n ×R n ) be such that its kernel has compact support in R n ×R n . With q ∈ S 1 ph (R n × R n ) its principal symbol. Let m ∈ C 1 (R n ). There exist K 1 and K 2 , compact operators on L 2 (R n ), with compactly supported kernels, such that [Q, m] = 1 i ∇ x m • Op(∇ ξ q) + K 1 = 1 i Op(∇ ξ q) • ∇ x m + K 2 . (3.3) Proof. Consider a sequence (m k ) k∈N ⊂ C ∞ (R n ) such that |α|≤1 ∂ α x (m k -m) L ∞ → 0 as k → +∞. Classical symbolic calculus gives [Q, m k ] = 1 i ∇ x m k • Op(∇ ξ q) + K k 1 , (3.4) with K k 1 = Op(r k 1 ) for some r k 1 ∈ S -1 ph , j = 1, 2. Thus, K k 1 is bounded from L 2 (R n ) into H 1 (R n ). In addition, since K k 1 has a kernel with compact supports in R n × R n , it is compact on L 2 (R n ). Note that the support of the kernel of K k 1 lies in a compact K of R n × R n that is uniform with respect to k.
On the other hand, observe that

∇ x m k • Op(∇ ξ q) → ∇ x m • Op(∇ ξ q) in L (L 2 (R n )).
Moreover, from Theorem 3.2 applied to m k -m, one also has

[Q, m k ] → [Q, m] in L (L 2 (R n )).
Using then (3.4) we deduce that (K k 1 ) n∈N converges to some K 1 in L (L 2 (R n )), and from the closedness of the set of compact operators in L (L 2 (R n )) we find that K 1 is compact. Moreover, K 1 has a kernel supported in K. The limits above give the first equality in (3.3). The second equality follows similarly.

Let Ω be a bounded open set of R n and (κ 0 , g 0 ) ∈ X 1 (Ω), with definition adapted from that of X 1 (M). The L 2 -inner product and norm are given by the density κ 0 µ g 0 . The following result is also a consequence of Theorem 3.2. Proposition 3.4. Let (u k ) k∈N ⊂ H 1 loc (Ω) be a sequence that converges weakly to 0 and let µ be a H 1 -microlocal defect density measure on S * Ω associated with the sequence

(u k ) k . Let b 1 ∈ W 1,∞ (R n ) and b 2 ∈ C 0 (R n ). Consider also Q 1 , Q 2 ∈ Ψ 1 ph (R n ), both with kernels compactly supported in Ω × Ω, with q 1 , q 2 ∈ S 1 ph (R n × R n ) for respective principal symbol. Then, one has b 1 Q 1 b 2 Q 2 u k , u k H -1 comp (Ω),H 1 loc (Ω) -→ k→+∞ µ, b 1 b 2 q 1 q 2 S * Ω . (3.5) More generally, assume that (b k 1 ) k∈N ⊂ W 1,+∞ (R n ) and (b k 2 ) k∈N ⊂ L ∞ (R n ), and (κ k , g k ) k∈N ⊂ Y(Ω) with b k 1 -b 1 W 1,+∞ (R n ) + b k 2 -b 2 L ∞ (R n ) + (κ k , g k ) -(κ 0 , g 0 ) Y(Ω) → 0, as k → +∞. Then b k 1 Q 1 b k 2 Q 2 u k , u k H -1 comp (Ω,κ k µg k ),H 1 loc (Ω,κ k µg k ) -→ k→+∞ µ, b 1 b 2 q 1 q 2 S * Ω . (3.6) Remark 3.5. Note that b 1 is chosen in W 1,∞ (R n ) because
one cannot multiply an element in H -1 by a bounded function. One derivative is needed.

Proof of Proposition 3.4. With Lemma 3.6 below we may replace the density κ k µ g k in the L 2 -inner product by κ 0 µ g 0 and thus in the H -1 comp -H 1 loc duality. We write

b k 1 Q 1 b k 2 Q 2 = b 1 Q 1 b 2 Q 2 + R k , R k = b 1 Q 1 (b k 2 -b 2 ) Q 2 + (b k 1 -b 1 )Q 1 b k 2 Q 2 . Note that R k maps H 1 loc (Ω) into H -1 comp (Ω) continuously.
Moreover because of the convergences of b k 1 and b k 2 , and the boundedness of (u k ) k∈N in H 1 loc (Ω) one finds that R k u k → 0 strongly in H -1 comp (Ω). Thus we can write According to Theorem 3.

b k 1 Q 1 b k 2 Q 2 u k , u k H -1 comp (Ω),H 1 loc (Ω) = b 1 Q 1 b 2 Q 2 u k , u k H -1 comp (Ω)
2 the commutator [b 1 , Q 1 ] is bounded on L 2 (Ω) implying that [b 1 , Q 1 ] b 2 Q 2 u k is bounded in L 2 (Ω) yielding [b 1 , Q 1 ] b 2 Q 2 u k , u k H -1 comp (Ω),H 1 loc (Ω) = ([b 1 , Q 1 ] b 2 Q 2 u k , u k ) L 2 (Ω) -→ k→+∞ 0,
since u k → 0 strongly in L 2 (Ω). We may thus assume that b 1 = 1 without any loss of generality.

Let ε > 0 and let b ε 2 ∈ C ∞ (Ω) be such that b 2 -b ε 2 L ∞ ≤ ε. Write Q 1 b 2 Q 2 = Q 1 b ε 2 Q 2 + R ε , R ε = Q 1 (b 2 -b ε 2 ) Q 2 . One has | R ε u k , u k H -1 comp (Ω),H 1 loc (Ω)
| ≤ Cε, and this leads to

Q 1 b 2 Q 2 u k , u k H -1 comp (Ω),H 1 loc (Ω) = Q 1 b ε 2 Q 2 u k , u k H -1 comp (Ω),H 1 loc (Ω) + o(1) ε→0 + o(1) k→+∞ . (3.7) Since b ε 2 is smooth, by symbolic calculus one has Q 1 b ε 2 Q 2 u k , u k H -1 comp (Ω),H 1 loc (Ω) -→ k→+∞ µ, b ε 2 q 1 q 2 S * Ω . (3.8)
Finally, since µ, b ε 2 q 1 q 2 S * Ω → µ, b 2 q 1 q 2 S * Ω as ε → 0, with (3.7) and (3.8) one concludes that (3.5) holds. Lemma 3.6. Assume that (κ k , g k ) -(κ 0 , g 0 ) Y(Ω) → 0 and consider a sequence

(f k , h k ) k∈N bounded in L 2 comp (Ω) ⊕ L 2 loc (Ω). Then (f k , h k ) L 2 (Ω,κ k µg k ) = (f k , h k ) L 2 (Ω) + o(1) k→+∞ . If (f k , h k ) k∈N is bounded in H -1 comp (Ω) ⊕ H 1 loc (Ω) then f k , h k H -1 comp (Ω,κ k µg k ),H 1 loc (Ω,κ k µg k ) = f k , h k H -1 comp (Ω),H 1 loc (Ω) + o(1) k→+∞ .
Here, Lemma 3.6 is written in the case of a bounded open set of the the Euclidean space but the same result holds in the case of a compact manifold.

Proof. One has µ g 0 = (det g 0 ) 1/2 dx and µ

g k = (det g k ) 1/2 dx. Therefore κ k µ g k = α k κ 0 µ g 0 with α k = κ k κ 0 det g k det g 0 1/2
and α k → 1 in the Lipschitz norm.

Measures and partial differential equations.

Microlocal defect measures associated with sequences of solutions of partial differential equations with smooth coefficients can have properties such as support localization in the characteristic set and invariance along the Hamiltonian flow. With the material developed above, we extend these results to the case of C 1 -coefficients. We focus on the case of wave operators. Proposition 3.7. Let (κ 0 , g 0 ) ∈ X 1 (M) and set p 0 (x, τ, ξ) = -τ 2 + g 0

x (ξ, ξ), that is, the principal symbol of

P 0 = P κ 0 ,g 0 . Let (κ k , g k ) k∈N ⊂ Y(M) be such that (κ k , g k ) -(κ 0 , g 0 ) Y(M) → 0 as k → +∞ and set P k = P κ k ,g k .
Consider a sequence (u k ) k∈N ⊂ H 1 loc (L) that converges to 0 weakly and µ a H 1 -microlocal defect density measure associated with (u k ) k∈N .

Let T 1 < T 2 . The following properties hold.

(

) If P k u k → 0 strongly in H -1 loc (T 1 , T 2 ) × M then (3.9) supp(µ) ∩ S * ((T 1 , T 2 ) × M) ⊂ Char(p 0 ). 1
(

) If moreover P k u k → 0 strongly in L 2 loc (T 1 , T 2 ) × M then one has (3.10) t H p 0 µ = 0 in the sense of distributions on S * (T 1 , T 2 ) × M , that is, µ, H p 0 q S * L = 0 for all q ∈ C ∞ c S * (T 1 , T 2 ) × M . 2
Since H p 0 is a tangent vector field on S * L where µ lives (see Section 2.1) note that t H p 0 µ makes sense in the second item of the proposition. Moreover note that H p 0 is a tangent vector field on S * L ∩ Char(p 0 ) and one has supp(µ) ∩ S * ((T 1 , T 2 ) × M) ⊂ Char(p 0 ) by the first item of the proposition. Finally, notice that for a Hamiltonian vector field, H p 0 = -t H p 0 as recalled in Section 2.1 even in the case (κ 0 , g 0 ) ∈ X 1 (M).

Naturally, Proposition 3.7 and its proof can be adapted to the other energy levels. We shall also need the following result. Let T 1 < T 2 . The following properties hold.

(

) If P k u k → 0 strongly in H -2 loc (T 1 , T 2 ) × M then supp(µ) ∩ S * ((T 1 , T 2 ) × M) ⊂ Char(p 0 ). (2) If moreover P k u k → 0 strongly in H -1 loc (T 1 , T 2 ) × M then one has t H p 0 µ = 0 in the sense of distributions on S * (T 1 , T 2 ) × M . Proof of Proposition 3.7. Consider B ∈ Ψ 0 c,ph (L) with kernel supported in (T 1 , T 2 )×M 2 and 1
b ∈ S 0 c,ph (L) its principal symbol. For the definition of the L 2 -inner product we use (κ 0 , g 0 ). We also use the partition of unity 1 = j∈J χ j with χ j ∈ C ∞ c (O j ) associated with the atlas A and the additional cutoff functions χj , χj ∈ C ∞ c (O j ) that are introduced in Section 3.2 and, as obtained in that section, we write

BP k u k , u k H -1 comp (L),H 1 loc (L) = j∈J χ j BP k u k , u k H -1 comp (L),H 1 loc (L) (3.11) = j∈J χ j χj BP k χj v k j , v k j H -1 comp (L),H 1 loc (L) + o(1) k→+∞ , with v k j = χj u k . Associated with (ϑ -1 j ) * v k j , the local representative of v k j , is a microlocal defect measure µ j in ϑ j (O j ) = Õj = R × Õj and χ C j j µ j is the local representative of χ j µ in this chart. See Section 3.2.
Note that we use local representatives of the operators, functions, measures without introducing any new symbols. Yet to keep clear that the analysis is carried out in a local chart we use the notation L 2 ( Õj ), H s ( Õj ) and not L 2 (L), H s (L). To further lighten notation we set κk = (det g k ) 1/2 κ k . One has

P k = ∂ 2 t -(κ k ) -1 p,q ∂ p κk g pq k ∂ q + 1 = Pk - p,q R p,q k , with Pk = ∂ 2 t -p,q ∂ p g pq k ∂ q + 1 and R p,q k = (κ k ) -1 [∂ p , κk ]g pq k ∂ q .
Note that χj BR p,q k χj defines a sequence of bounded operators from H 1 (L) into L 2 (L), uniformly with respect to k. Consequently, one has

χ j χj BR p,q k χj v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = χ j χj BR p,q k χj v k j , v k j L 2 ( Õj )
→ k→+∞ 0 since v k j converges strongly to 0 in L 2 ( Õj ). This leads to 1) k→+∞ , by Proposition 3.4. Since χ j µ j = χ j µ locally, lifting back the analysis to the manifold level, with (3.11), one finds

χ j χj BP k χj v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = χ j χj B Pk χj v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) + o(1) k→+∞ = µ j , χ j bp 0 S * ( Õj ) + o(
BP k u k , u k H -1 comp (L),H 1 loc (L) = j∈J µ, χ j bp 0 S * (L) = µ, bp 0 S * (L) + o(1) k→+∞ .

Now, one has

BP k u k , u k H -1 comp (L),H 1 loc (L) = P k u k , t Bu k H -1 loc (L),H 1 comp (L) + o(1) k→+∞ ,
with the transpose operator t B bounded from

H 1 loc (L) into H 1 comp (L) since B is itself bounded from H -1 loc (L) into H -1 comp (L). If one assumes that P k u k → 0 strongly in H -1 loc (T 1 , T 2 ) × M one obtains BP k u k , u k H -1 comp (L),H 1 loc (L) → k→+∞ 0,
and thus µ, bp 0 S * (L) = 0, ∀b ∈ S 0 c,ph (L) with supp(b) ⊂ T * (T 1 , T 2 ) × M , and one obtains the support estimation (3.9).

We now prove the second item of the proposition. We assume that P k u k lies in L 2 loc (T 1 , T 2 )× M) and converges strongly to 0 in this space. Consider B ∈ Ψ 1 c,ph (L) with kernel supported in (T 1 , T 2 ) × M 2 and b ∈ S 1 c,ph (L) its principal symbol. We are interested in the limit of

[P k , B]u k , u k H -1 comp (L),H 1 loc (L)
, which makes sense since [P k , B] is of order 2. We have

[P k , B]u k = P k Bu k -BP k u k ∈ H -1 (T 1 , T 2 ) × M). Since P k u k lies in L 2 (T 1 , T 2 ) × M) by assumption then BP k u k lies in H -1 (T 1 , T 2 ) × M
) and the same holds for P k Bu k . We may thus write

[P k , B]u k , u k H -1 comp (L),H 1 loc (L) = P k Bu k , u k H -1 comp (L),H 1 loc (L) -P k u k , B * u k L 2 loc (L),L 2 comp (L)
, where the adjoint is computed with respect to the L 2 -inner product associated with (k 0 , g 0 ) here. As B maps continuously

L 2 loc (T 1 , T 2 ) × M) into H -1 comp (T 1 , T 2 ) × M) then B * maps continuously H 1 loc (L) into L 2 comp (L). Thus, one has P k u k , B * u k L 2 (L) → k→+∞ 0.
By Lemma 3.6 it is asymptotically equivalent to use (κ 0 , g 0 ) or (κ k , g k ) for the definition of the L 2 -inner product and H -1 comp -H 1 loc duality, that is,

P k Bu k , u k H -1 comp (L),H 1 loc (L) = P k Bu k , u k H -1 comp (L,κ k µg k dt),H 1 loc (L,κ k µg k dt) + o(1)
k→+∞ . Since P k is selfadjoint for this latter L 2 -inner product, one obtains

P k Bu k , u k H -1 comp (L),H 1 loc (L) = Bu k , P k u k L 2 comp (L,κ k µg k dt),L 2 loc (L,κ k µg k dt) + o(1) k→+∞ = Bu k , P k u k L 2 comp (L),L 2 loc (L) + o(1) k→+∞ . Using again that P k u k → 0 strongly to 0 in L 2 loc (T 1 , T 2 ) × M) we obtain P k Bu k , u k H -1 comp (L),H 1 loc (L) → k→+∞ 0,
and finally

[P k , B]u k , u k H -1 comp (L),H 1 loc (L) → k→+∞ 0. (3.12)
As above, with the partition of unity 1 = j∈J χ j we write

[P k , B]u k , u k H -1 comp (L),H 1 loc (L) = j∈J χ j [P k , B]u k , u k H -1 comp (L),H 1 loc (L) . (3.13)
For each term in the sum one has

χ j [P k , B]u k , u k H -1 comp (L),H 1 loc (L) = χ j [P k , Bj ]v k j , v k j H -1 comp (L),H 1 loc (L) + o(1)
k→+∞ . with Bj = χj B χj . This allows one to work in a local chart and write

[P k , B]u k , u k H -1 comp (L),H 1 loc (L) = j∈J χ j [P k , Bj ]v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) , (3.14)
with the (manifold-local chart) identifications described above. With

A k = A κ k ,g k , in the local chart C j one writes χ j [P k , Bj ] = χ j [∂ 2 t , Bj ] -χ j [A k , Bj ] = χ j [∂ 2 t , Bj ] - 1≤p,q≤d Q pq 1 + Q pq 2 + Q pq 3 + Q pq 4 ,
with

Q pq 1 = χ j κ-1 k ∂ xp κk g pq k [∂ xq , Bj ], Q pq 2 = χ j κ-1 k ∂ xp [κ k g pq k , Bj ]∂ xq , Q pq 3 = χ j κ-1 k [∂ xp , Bj ]κ k g pq k ∂ xq , Q pq 4 = χ j [κ -1 k , Bj ]∂ xp κk g pq k ∂
xq . We now compute the limit of each term associated with this decomposition of [P k , Bj ] on the right-hand side of (3.14). The principal symbol of χ j [∂ 2 t , Bj ] is iχ j {τ 2 , b} and thus

χ j [∂ 2 t , Bj ]v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = µ j , iχ j {τ 2 , b} S * ( Õj ) + o(1) k→+∞ . Proposition 3.

applies and yields

Q pq 1 v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = µ j , iχ j g 0,pq ξ p ∂ xq b S * ( Õj ) + o(1) k→+∞ , and
Q pq 3 v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = µ j , iχ j g 0,pq (∂ xp b)ξ q S * ( Õj ) + o(1) k→+∞ . With Theorem 3.2 one has [κ k g pq k , Bj ] → [κ 0 g 0,pq , Bj ] in L L 2 ( Õj ) as k → +∞. It follows that Q pq 2 v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = Q pq 2,a v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) + o(1) k→+∞ , with Q pq 2,a = χ j κ-1 k ∂ xp [κ 0 g 0,pq , Bj ]∂ xq . With Corollary 3.3 one writes [κ 0 g 0,pq , Bj ] = - 1 i ∇ x (κ 0 g 0,pq ) • Op ∇ ξ ( χ2 j b) + K 1 ,
with K 1 a compact operator on L 2 (R d+1 ), with compactly supported kernel. One thus obtains

Q pq 2 v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = Q pq 2,b v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) + o(1) k→+∞ , with Q pq 2,b = - 1 i χ j κ-1 k ∂ xp ∇ x (κ 0 g 0,pq ) • Op ∇ ξ ( χ2 j b) ∂ xq . Proposition 3.

applies and yields

Q pq 2 v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = µ j , -iχ j ξ p ξ q (κ 0 ) -1 ∇ x (κ 0 g 0,pq ) • ∇ ξ b S * ( Õj ) + o(1) k→+∞ .
We now treat the term associated with Q pq 4 . Note that one has p,q Q pq 4 = χ j [κ -1 k , Bj ]κ k A k . We write, lifting temporarily the analysis back to the manifold, p,q

Q pq 4 v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = χ j [κ -1 k , B]κ k A k v k j , v k j H -1 comp (L),H 1 loc (L) = χ j [κ -1 k , B]κ k A k u k , u k H -1 comp (L),H 1 loc (L) + o(1) k→+∞ . Setting f k = (∂ 2 t -A k )u k with f k → 0 strongly in L 2 loc (T 1 , T 2 ) × M , we thus find p,q Q pq 4 v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = χ j [κ -1 k , B]κ k ∂ 2 t u k , u k H -1 comp (L),H 1 loc (L) -χ j [κ -1 k , B]κ k f k , u k H -1 comp (L),H 1 loc (L) + o(1) k→+∞ = χ j [κ -1 k , Bj ]κ k ∂ 2 t v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) + o(1)
k→+∞ , bringing again the analysis at the level of the local chart.

Using that κk is independent of t we may write

χ j [κ -1 k , Bj ]κ k ∂ t = χ j ∂ t [κ -1 k , Bj ]κ k + χ j [κ -1 k , E j ]κ k ,
where

E j = [∂ t , Bj ] ∈ Ψ 1 c,ph ( Õj ), with ∂ t b ∈ S 1 c,ph ( Õj ) for principal symbol. With Theo- rem 3.2 we see that [κ -1 k , E j ] maps L 2 ( Õj ) into itself continuously and moreover [κ -1 k , E j ] → [(κ 0 ) -1 , E j ] in L L 2 ( Õj ) . Thus we obtain χ j [κ -1 k , E j ]κ k ∂ t v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = χ j [(κ 0 ) -1 , E j ]κ k ∂ t v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) + o(1) k→+∞ → k→+∞ 0,
arguing as above. Similarly we write

χ j ∂ t [κ -1 k , Bj ]κ k ∂ t v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) ∼ k→+∞ χ j ∂ t [(κ 0 ) -1 , Bj ]κ 0 k ∂ t v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj )
Arguing as we did for the term associated with Q p,q 2 we thus find

χ j ∂ t [κ -1 k , Bj ]κ k ∂ t v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = µ j , -iχ j τ 2 κ0 (∇ x (κ 0 ) -1 ) • ∇ ξ b S * ( Õj ) + o(1) k→+∞ .
Collecting the various estimates we found we obtain

χ j [P k , Bj ]v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = µ j , χ j σ S * ( Õj ) + o(1) k→+∞ . (3.15)
with

σ = i{τ 2 , b} -i p,q g 0,pq ξ p ∂ xq b + g 0,pq (∂ xp b)ξ q -ξ p ξ q (κ 0 ) -1 ∇ x (κ 0 g 0,pq ) • ∇ ξ b + iτ 2 κ0 (∇ x (κ 0 ) -1 ) • ∇ ξ b.
Recalling that p 0 = -τ 2 + p,q g 0,pq ξ p ξ q one finds σ = -i{p 0 , b}

+ ip 0 (κ 0 ) -1 ∇ x (κ 0 ) • ∇ ξ b.
Since µ, and thus µ j , is supported in Char(p 0 ) by the first part of the proposition, one concludes that The strategy we follow is very much inspired by the Melrose and Sjöstrand approach to the propagation of singularities [START_REF] Melrose | Singularities of boundary value problems I[END_REF] and relies on careful choices of test functions allowing one to construct sequences of points in the support of the measure relying on nonnegativity 2 . Then, a limiting procedure leads to the conclusion, in the spirit of the classical proof of the Cauchy-Peano theorem.

χ j [P k , Bj ]v k j , v k j H -1 comp ( Õj ),H 1 loc ( Õj ) = -i µ j , χ j {p 0 , b} S * ( Õj ) + o(1) k→+∞ . Since χ j µ = χ j µ j (see Section 3.2), with (3.13)-(3.14) one obtains [P k , B]u k , u k H -1 comp (L),H 1 loc (L) = -i µ, {p 0 , b} S * (L) + o(1) k→+∞ . With (3.
The proof of Theorem 1.10 is made of two steps that are stated in the following propositions. (1) The set F is a union of maximally extended integral curves of the vector field X.

(2) For any compact K ⊂ Ω where the vector field X does not vanish,

∀ε > 0, ∃δ 0 > 0, ∀x ∈ K ∩ F, ∀δ ∈ [-δ 0 , δ 0 ], B x + δX(x), δε ∩ F = ∅.
Proposition 4.2. Let X be a C 0 -vector field on Ω an open set of R d . Consider a nonnegative measure µ on Ω solution to t Xµ = 0 in the sense of distributions, that is,

(4.1) t Xµ, a D (Ω),C ∞ c (Ω) = µ, Xa D ,0 (Ω),C 0 c (Ω) = 0, a ∈ C ∞ c (Ω).
Then, the closed set F = supp(µ) satisfies the second property in Proposition 4.1.

Proof of Proposition 4.1. First, we prove that Property (1) implies Property (2) and consider a compact set K of R d such that K ⊂ Ω and K ∩ F = ∅.

There exists η > 0 such that K ⊂ K η ⊂ Ω with K η = {x ∈ Ω; dist(x, K) ≤ η}. One has X ≤ C 0 on K η for some C 0 > 0. Let x ∈ K and let γ(s) be a maximal integral curve defined on an interval ]a, b[, a, b ∈ R and such that 0 ∈]a, b[ and γ(0

) = x. If b < ∞ then there exists s 1 ∈ ]0, b[ such that γ(s 1 ) / ∈ K η . Since γ(s) ∈ K η if s < η/C 0 , one finds that b ≥ η/C 0 .
Similarly, one has |a| ≥ η/C 0 . Consequently, there exists S > 0 such that any maximal integral curve γ(s) of the vector field X with γ(0) ∈ K is defined for s ∈ I = (-S, S).

2 of the measure in our case and of some operators for Melrose and Sjöstrand, via the Gårding inequality. By uniform continuity of the vector field X in a compact neighborhood of K we have

γ(s) = γ(0) + s 0 γ(s)ds = γ(0) + s 0 X(γ(s))ds = x + sX(x) + r(s), s ∈ (-S, S)
where lim s→0 r(s) /s = 0, uniformly with respect to x. We deduce that for any ε > 0 there exists 0 < δ 0 < S such that r(s) < sε for any s ∈ (-δ 0 , δ 0 ), which implies

F γ(s) ∈ B x + sX(x), sε .
Second, we prove that Property (2) implies Property [START_REF] Ambrosio | Continuity equations and ODE flows with non-smooth velocity[END_REF]. It suffices to prove that for any x ∈ F there exist an interval I 0 and an integral curve γ : I → F such that γ(s) = X(γ(s)) and γ(0) = x.

Then, the standard continuation argument shows that this local integral curve included in F can be extended to a maximal integral curve also included in F .

If X(x) = 0, then the trivial integral curve γ(s) = x, s ∈ R, is included in F . As a consequence, we assume X(x) = 0 and we pick a compact neighborhood K of x containing B(x, η) with η > 0 and where, for some 0

< c K < C K , c K ≤ X(y) ≤ C K , y ∈ K.
Let n ∈ N * . Set x n,0 = x and ε = 1/n and apply Property (2). One deduces that there exist 0 < δ n ≤ 1/n and a point

x n,1 ∈ F ∩ B x n,0 + δ n X(x n,0 ), δ n /n .
If x n,1 ∈ K one can perform this construction again, yet starting from x n,1 instead of x n,0 . If a sequence of points x n,0 , x n,1 , . . . , x n,L + is obtained in this manner one has

x n, +1 ∈ F ∩ B x n, + δ n X(x n, ), δ n /n , = 0, . . . , L + -1. (4.2)
One can carry on the construction as long as x n,L + ∈ K. We can perform the same construction for ≤ 0, with the property

x n, -1 ∈ F ∩ B x n, -δ n X(x n, ), δ n /n , | | = 0, . . . , L --1. (4.3)
Having X ≤ C K on K and B(x, η) ⊂ K ensures that we can construct the sequence at least for

L + = L -= L n = η δ n (C K + 1) + 1 ≤ η δ n (C K + 1/n) + 1,
where . denotes the floor function. With the points x n, , | | ≤ L n , we have constructed we define the following continuous curve γ n (s) for |s| ≤ L n δ n :

γ n (s) = x n, + (s -δ n ) x n, +1 -x n, δ n for s ∈ [ δ n , ( + 1)δ n ) and | | ≤ L n -1.
This curve and its construction is illustrated in Figure 1(a). Note that γ n (s) remains in a compact set, uniformly with respect to n. In this compact set X is uniformly continuous.

We set S = η/(C K + 1). Since S ≤ L n δ n we shall in fact only consider the function γ n (s) for |s| ≤ S in what follows. Note that since

x n, ∈ F for | | ≤ L n then one has dist γ n (s), F ≤ δ n (C K + 1/n), |s| ≤ S. (4.4)
From (4.2), for ≥ 0 and s ∈ ( δ n , ( + 1)δ n ), we have

γn (s) = x n, +1 -x n, δ n = X(x n, ) + O(1/n).
Similarly, from (4.3), for ≤ 0 and s ∈ (( -1)δ n , δ n ), we have

γn (s) = x n, -x n, -1 δ n = X(x n, ) + O(1/n).
In any case, using the uniform continuity of the vector field X, we find γn (s) = X(γ n (s)) + e n (s),

where the error |e n | goes to zero uniformly with respect to |s| ≤ S as n → +∞.

Since the curve γ n is absolutely continuous (and differentiable except at isolated points), we find We now let n grow to infinity. With (4.5), the family of curves (s → γ n (s), |s| ≤ S) n∈N * is equicontinuous and pointwise bounded; by the Arzelà-Ascoli theorem we can extract a subsequence (s → γ np ) p∈N that converges uniformly to a curve γ(s), |s| ≤ S. Convergence is illustrated in Figure 1(b). Passing to the limit n p → +∞ in (4.5) we find that γ(s) is solution to

γ(s) = x + s 0 X(γ(σ))dσ.
From estimation (4.4), for any |s| ≤ S, there exists (y p ) p ⊂ F such that lim p→+∞ y p = γ(s).

Since F is closed we conclude that γ(s) ∈ F .

Positivity argument and proof of Proposition 4.2. We consider a compact set K where the vector field X does not vanish. By continuity of the vector field there exist 0 < c K ≤ C K such that 0 < c K ≤ X(x) ≤ C K , for all x ∈ K.

Let us consider x 0 ∈ K ∩ supp(µ). By performing a rotation and a dilation of coefficient X(x 0 ) ∈ [c K , C K ], we can assume that X(x 0 ) = (1, 0, . . . , 0) ∈ R d . We shall write x = (x 1 , x ) with x ∈ R d-1 .

Let χ ∈ C ∞ (R) be given by (4.6) χ(s) = 1 s<1 exp(1/(s -1)), and β ∈ C ∞ (R) be such that (4.7)

β ≡ 0 on ] -∞, -1], β > 0 on ] -1, -1/2[, β ≡ 1 on [-1/2, +∞[.
We then set

q ε,δ,x 0 = (χ • v)(β • w), g ε,δ,x 0 = (χ • v)(β • w)Xv, h ε,δ,x 0 = (χ • v)(β • w)Xw, (4.8) with v(x) = 1/2 -δ -1 (x 1 -x 0 1 ) + 8(εδ) -2 x -x 0 2 and w(x) = 2ε -1 1 -δ -1 (x 1 -x 0 1
) , for ε > 0 and δ > 0 both meant to be chosen small in what follows. We have Xq ε,δ,x 0 = g ε,δ,x 0 + h ε,δ,x 0 .

The function q ε,δ,x 0 is compactly supported. Indeed, in the support of β • w one has w ≥ 1 implying

x 1 -x 0 1 ≤ δ(1 + ε/2), while on the support of χ • v one has v ≤ 1 which gives -1/2 + 8(εδ) -2 x -x 0 2 ≤ δ -1 (x 1 -x 0 1
). On the supports of q ε,δ,x 0 and (χ • v)(β • w) one thus finds (4.9)

-δ/2 ≤ x 1 -x 0 1 ≤ δ(1 + ε/2) and 8(εδ) -2 x -x 0 2 ≤ 3/2 + ε/2. Similarly, on the support of β • w one has -1 ≤ w ≤ -1/2

δ(1 + ε/4) ≤ x 1 -x 0 1 ≤ δ(1 + ε/2
), which implies that on the support of h ε,δ,x 0 one has (4.10)

δ(1 + ε/4) ≤ x 1 -x 0 1 ≤ δ(1 + ε/2
) and 8(εδ) -2 x -x 0 2 ≤ 3/2 + ε/2. In particular, in the case ε ≤ 1, one finds (4.11) supp(h ε,δ,x 0 ) ⊂ B x 0 + δX(x 0 ), εδ .

These estimations of the supports of q ε,δ,x 0 and h ε,δ,x 0 are illustrated in Figure 2.

Lemma 4.3. For any 0 < ε ≤ 1 there exists δ 0 > 0 such that for any x 0 ∈ K and 0 < δ ≤ δ 0 , the function g ε,δ,x 0 is nonnegative. Moreover, g ε,δ,x 0 is positive in a neighborhood of x 0 .

Proof. Let 0 < ε ≤ 1. We have g ε,δ,x 0 = (χ • v)(β • w)Xv. Since β ≥ 0 and χ < 0 it suffices to prove that Xv(x) ≤ 0 for x in the support of (χ • v)(β • w) for δ > 0 chosen sufficiently small, uniformly with respect to x 0 ∈ K. We write with α 1 (x, x 0 ) ∈ R and α (x, x 0 ) ∈ R d-1 . By (4.9), for x ∈ supp(χ • v)(β • w) we have x -x 0 δ. From the uniform continuity of X in any compact set we conclude that

X(x) -X(x 0 ) = α 1 (x, x 0 )∂ x 1 + α (x, x 0 ) • ∇ x , x 0 x x 1 X(x 0 ) = ∂ x1 εδ 2 δ -δ/2 εδ 2 (a) Support of q ε,δ,x 0 . x 0 x x 1 X(x 0 ) δ -δ/2 εδ 4 εδ εδ 2 (b) Support of h ε,δ,x 0 .
|α 1 (x, x 0 )| + α (x, x 0 ) = o(1) as δ → 0 + , (4.12
) uniformly3 with respect to x 0 ∈ K and x ∈ supp(χ • v)(β • w). Using that X(x 0 ) = ∂ x 1 and the form of v given above, we write

Xv(x) = X(x)v (x) = ∂ x 1 v + X(x) -X(x 0 ) v (x) = -δ -1 1 + α 1 (x, x 0 ) -16ε -1 (εδ) -1 α (x, x 0 ) • (x -x 0 ) .
Using again (4.9), we thus find for x ∈ supp(χ

• v)(β • w) α 1 (x, x 0 ) -16ε -1 (εδ) -1 α (x, x 0 ) • (x -x 0 ) α 1 (x, x 0 ) + ε -1 α (x, x 0 ) .
With ε fixed above and with (4.12) we find that Xv(x) ∼ -δ -1 as δ → 0 + uniformly with respect to x 0 ∈ K and x ∈ supp(χ • v)(β • w). Finally, we have g ε,δ,x 0 (x 0 ) = -δ -1 χ (1/2)β(2ε -1 ) > 0 and thus g ε,δ,x 0 is positive in a neighborhood of x 0 .

We are now in a position to conclude the proof of Proposition 4.2. Note that it suffices to prove the result for 0 < ε ≤ 1. We choose δ 0 > 0 as given by Lemma 4.3. Let then x 0 ∈ K ∩ supp(µ). We apply (4.1) to the family q ε,δ,x 0 of test functions with 0 < δ ≤ δ 0 :

(4.13) 0 = µ, X(q ε,δ,x 0 ) = µ, g ε,δ,x 0 + µ, h ε,δ,x 0 .
By Lemma 4.3, g ε,δ,x 0 ≥ 0 and g ε,δ,x 0 is positive in a neighborhood of x 0 . As x 0 ∈ supp(µ) we find µ, g ε,δ,x 0 > 0. Consequently, µ, h ε,δ,x 0 = 0. By the support estimate for h ε,δ,x 0 given in (4.11) the conclusion follows: supp(µ) ∩ B x 0 + δX(x 0 ), εδ = ∅.

Exact controllability: proof of Theorem 1.12

Let (κ 0 , g 0 ) ∈ X 1 (M) and assume that (ω, T ) fulfills the geometric control condition of Definition 1.8 .

Let also (κ, g) ∈ Y(M). With Proposition 1.6, the result of Theorem 1.12 follows if we prove that there exists ε > 0 and C obs > 0 such that E κ,g (u)(0) ≤ C obs 1 (0,T )×ω ∂ t u 2 L 2 (L,κµgdt) , for any weak solution u of the wave equation associated with (κ, g) chosen such that

(κ, g) -(κ 0 , g 0 ) Y(M) ≤ ε.
The L 2 -norm on the r.h.s. is associated with (κ, g), that is,

1 (0,T )×ω ∂ t u 2 L 2 (L,κµgdt) = T 0 ω |∂ t u| 2 κµ g dt.
Yet, for ε > 0 chosen sufficiently small one has . L 2 (L,κ 0 µ g 0 ) . L 2 (L,κµg) , where A B means c 1 ≤ A/B ≤ c 2 for some c 1 , c 2 > 0. In other words, we have equivalence with constants uniform with respect to (κ, g). In what follows, L 2 -and more generally H s -norms on M are chosen with respect to κ 0 µ g 0 unless explicitely written. Our goal is thus to prove the following observability inequality

(5.1)

E κ 0 ,g 0 (u)(0) ≤ C obs 1 (0,T )×ω ∂ t u 2 L 2 (L)
. The Bardos-Lebeau-Rauch uniqueness compactness argument reduces the proof of (5.1) to the proof of the weaker estimate (5.2)

E κ 0 ,g 0 (u)(0) ≤ C 1 (0,T )×ω ∂ t u 2 L 2 (L) + C u(0), ∂ t u(0) 2 L 2 (M)⊕H -1 (M)
, that exhibits an additional compact term, and expresses observability for high-frequencies. Low frequencies are dealt with by means of a unique continuation argument.

To prove (5.2) we argue by contradiction and we assume that there exists a sequence

(κ k , g k ) k∈N ⊂ Y(M) such that (5.3) lim k→+∞ (κ k , g k ) -(κ 0 , g 0 ) Y(M) = 0,
and yet for each k ∈ N the associated observability inequality does not hold. Thus, for each k ∈ N, there exists a sequence of initial data v k,p,0 , v k,p,1 p∈N ⊂ H 1 (M) × L 2 (M) with associated solution (v k,p ) p∈N , that is,

P k v k,p = 0 in (0, +∞) × M, v k,p |t=0 = v k,p,0 , ∂ t v k,p |t=0 = v k,p,1 in M, with P k = P κ k ,g k , that moreover has the properties E κ 0 ,g 0 (v k,p )(0) = 1 and 1 (0,T )×ω ∂ t v k,p L 2 (L) + v k,p,0 , v k,p,1 L 2 (M)⊕H -1 (M) ≤ 1 p + 1 .
We take p = k and we set u k,0 , u k,1 = v k,k,0 , v k,k,1 and u k = v k,k , one obtains P k u k = 0 in L and

(5.4) E κ 0 ,g 0 (u k )(0) = 1 and

1 (0,T )×ω ∂ t u k L 2 (L) + u k,0 , u k,1 L 2 (M)⊕H -1 (M) ≤ 1 k + 1 .
From (5.4) one has u k 0 weakly in H 1 loc (L). With (3.1)-(3.2), we can associate with (a subsequence of) (u k ) k a H 1 -microlocal defect measure µ on S * (L). Here, the measure is understood with respect to L 2 (L, κ 0 µ g 0 dt). From the second part of (5.4) one has µ = 0 in S * ((0, T ) × ω).

(5.5)

In fact, for any ψ ∈ C ∞ ((0, T ) × ω) one has ψ∂ t u k L 2 (L) ∼ 0 and thus µ, τ 2 ψ 2 = 0. Hence, supp(µ) ∩ S * ((0, T ) × ω) ⊂ {τ = 0}. Since {τ = 0} ∩ Char(p 0 ) ∩ S * (L) = ∅ with (3.9) one obtains (5.5).

With the first part of (5.4) one has the following lemma.

Lemma 5.1. The measure µ does not vanish on S * L .

A proof is given below.

We now use Proposition 3.7 to obtain a precise description of the measure µ. First, one has supp(µ) ∩ S * ((0, T ) × M) ⊂ Char(p 0 ). Furthermore, one has t H p 0 µ = 0 in the sense of distributions on S * ((0, T ) × M). Since H p 0 is a C 0 -vector field on the manifold S * L, Theorem 1.10 implies that supp(µ) is a union of maximally extended bicharacteristics in S * ((0, T ) × M).

Under the geometric control condition of Definition 1.8 , any maximal bicharacteristic meets S * ((0, T ) × ω) where µ vanishes by (5.5). Thus supp(µ) = ∅ yiedling a contradiction with the result of Lemma 5.1. We thus obtain that (5.1) holds. This concludes the proof of Theorem 1.12.

Proof of Lemma 5.1. Let T 1 < T 2 and φ ∈ C ∞ c (R) nonnegative and equal to 1 on a neighborhood of [T 1 , T 2 ]. On L, consider the elliptic operator Q = -∂ 2 t -A κ 0 ,g 0 + 1 with symbol q = τ 2 + p,q g 0 p,q (x)ξ p ξ q . Taking (3.2) and Lemma 3.6 into account one can write

φ 2 Qu k , u k H -1 comp (L),H 1 loc (L) ∼ k→+∞ µ, φ 2 q S * L . (5.6)
Integrating by parts one obtains

φ 2 Qu k , u k H -1 comp (L),H 1 loc (L) = L φ(t) 2 |∂ t u k | 2 + g 0 (∇ g 0 u k , ∇ g 0 u k ) + |u k | 2 κ 0 µ g 0 dt + 2(φ φ ∂ t u k , u k ) L 2 (L) = R φ(t) 2 E κ 0 ,g 0 (u k )(t)dt + 2(φ φ ∂ t u k , u k ) L 2 (L) .
Since the energy built on κ p , g p is preserved by the evolution given by P p , we have by (5.4)

(5.7)

E κ 0 ,g 0 (u k )(t) = E κ k ,g k (u k )(t) + o(1) = E κ k ,g k (u k )(0) + o(1) = E κ 0 ,g 0 (u k )(0) + o(1) = 1 + o(1)
and since (φ φ ∂ t u k , u k ) L 2 (L) → 0 as u k → 0 strongly in L 2 loc (L), one obtains

φ 2 Qu k , u k H -1 comp (L),H 1 loc (L) ∼ k→+∞ φ 2 L 2 (R) .
With (5.6) this proves that µ = 0.

6. Lack of continuity of the control operator with respect to coefficients 6.1. Proof of Theorems 1.14 and 1.14 . We prove the result of both theorems, that is, in the case k ≥ 1. In the case k = 1 we are simply required to prove additionnally that the geometric control condition of Definition 1.8 is fulfilled for geodesics given by the chosen metric g; see Remark 1.17. Let ε > 0. We set g = (1 + ε)g. Given any neighborhood U of (κ, g) in X k (M), for ε > 0 chosen sufficiently small one has (κ, g) ∈ U.

Moreover, observe that for ε > 0 chosen sufficiently small geodesics associated with g can be made arbitrary close to those associated with g uniformly in t ∈ [0, T ]. Hence, for such ε > 0 the geometric control condition is fulfilled for geodesics associated with g.

Observe that one has Char(p κ,g ) ∩ Char(p κ,g ) ∩ S * L = ∅. (6.1)

We consider a sequence (y k,0 , y k,1 ) (0, 0) weakly in H 1 (M) ⊕ L 2 (M) such that 1 2 ( y k,0 2 H 1 (M) + y k,1 2 L 2 (M) ) = 1.

L 2 -and H 1 -norms are based on the κµ g dt measure on L. Setting f k κ,g = H κ,g (y k,0 , y k,1 ) ∈ L 2 ((0, T ) × M) with H κ,g defined in (1.7), one obtains a sequence of control functions. According to the HUM method [START_REF] Lions | Contrôlabilité exacte, Stabilisation et Perturbations de Systèmes Distribués. Tome 1. Contrôlabilité exacte[END_REF], f k κ,g is itself a (weak) solution to the following free wave equation (6.2) P κ,g f k κ,g = 0, in the energy space L 2 (M) ⊕ H -1 (M), that is, (f k κ,g (0), ∂ t f k κ,g (0)) ∈ L 2 (M) × H -1 (M). Moreover, (f k κ,g (0), ∂ t f k κ,g (0)) depend continuously on (y k,0 , y k,1 ). The function f k κ,g is thus bounded in C 0 ((T 1 , T 2 ), L 2 (M)) uniformly with respect to k for any T 1 < T 2 . Since the map H κ,g is continuous f k κ,g 0 weakly in L 2 loc (L). Up to extraction of a subsequence, it is associated with a L 2 -microlocal defect measure µ f . With Proposition 3.7 one has (6.3) supp(µ f ) ⊂ Char(p κ,g ).

We consider the sequences of solutions (y k ) k and (ỹ k ) k to P κ,g y k = 1 (0,T )×ω f k κ,g in L, (y k , ∂ t y k ) |t=0 = (y k,0 , y k,1 ) in M, P κ,g ỹk = 1 (0,T )×ω f k κ,g in L, (ỹ k , ∂ t ỹk ) |t=0 = (y k,0 , y k,1 ) in M.

Both are bounded and weakly converges to 0 in H 1 loc (L). Up to extraction of subsequences, both are associated with H 1 -microlocal defect density measures µ and μ respectively. Since 1 (0,T )×ω f k κ,g 0 weakly in L 2 loc (L) then 1 (0,T )×ω f k κ,g → 0 strongly in H -1 loc (L) and, with Proposition 3.7, one finds supp(μ) ⊂ Char(p κ,g ). Thus one has supp(μ) ∩ supp(µ f ) = ∅. (6.4) The sequence (∂ t ỹk ) converges to 0 weakly in L 2 loc (L) and can be associated with a L 2microlocal defect density measure whose support is given by supp(μ). Lemma 6.1. One has (1 (0,T )×ω f k κ,g , ∂ t ỹk ) L 2 (L,κµ g dt) → 0 as k → +∞.

A proof is given below. Using the density of strong solutions of the wave equation, with integration by parts, one finds the following classical energy estimate E κ,g (ỹ k )(T ) -E κ,g (ỹ k )(0) = (1 (0,T )×ω f k κ,g , ∂ t ỹk ) L 2 (κµ g dt) .

With Lemma 6.1 one obtains

E κ,g (ỹ k )(T ) ∼ k→+∞ E κ,g (ỹ k )(0).
With the form of g chosen above one has

E κ,g (ỹ k )(t) = (1 + O(ε))E κ,g (ỹ k )(t),
uniformly with respect to t ∈ [0, T ]. Chosing ε > 0 sufficiently small and k sufficiently large, the first part of Theorem 1.14 follows since E κ,g (ỹ k )(0) = 1.

We use the values of ε and k chosen above. To prove (1.10), we write ỹk in the form ỹk = v 1 + v 2 where v 1 and v 2 are solution to (6.5) P κ,g v 1 = 1 (0,T )×ω f k κ,g in L, (v 1 , ∂ t v 1 ) |t=0 = (y k,0 , y k,1 ) in M, P κ,g v 2 = 1 (0,T )×ω (f k κ,g -f k κ,g ) in L, (v 2 , ∂ t v 2 ) |t=0 = (0, 0) in M.

with f k κ,g = H κ,g (y k,0 , y k,1 ). A hyperbolic energy estimation for v 2 solution to the second equation in (6.5) gives

E κ,g (v 2 )(T ) ≤ C T 1 (0,T )×ω (f k κ,g -f k κ,g ) 2 L 2 (L) .
Since one has (v 1 (T ), ∂ t v 1 (T )) = (0, 0) because of the definition of f k κ,g one finds E κ,g (v 2 )(T ) = E κ,g (ỹ k )(T ) ≥ 1/2, which gives the second result of Theorem 1.14.

Proof of Lemma 6.1. The key point in the proof is the following lemma. Lemma 6.2 ([18, Proposition 3.1]). Assume that u k and v k are two sequences bounded in L 2 loc that converge weakly to zero and are associated with defect measures µ and ν respectively. Assume that µ ⊥ ν, that is, µ and ν are supported on disjoint sets. Then, for any ψ ∈ C 0 c , lim k→+∞ (ψu k , v k ) L 2 = 0.

To apply this result, we just need to exchange the rough cutoff 1 (0,T )×ω for a smooth cutoff ψ(t, x). First, note that one has (1 (0,T )×ω f k κ,g , ∂ t ỹk ) L 2 (L,κµ g dt) ∼

k→+∞ (1 (0,T )×ω f k κ,g , ∂ t ỹk ) L 2 (L,κµgdt) .

We may thus simply consider the L 2 -norm and inner product associated with κµ g dt. Second, let δ > 0. Since (f k κ,g ) k and (ỹ k ) k are both bounded in C 0 ((0, T ), L 2 (M)) uniformly with respect to k, there exists 0 < T 1 < T 2 < T and O ω such that With (6.4) and Lemma 6.2, one finds (ψf k κ,g , ∂ t ỹk ) L 2 (L) → k→+∞ 0, (6.6) and the conclusion of the lemma follows.

6.2. Proof of Theorem 1.19. We consider first the case α = 1. As proven in [START_REF] Dehman | Analysis of the HUM control operator and exact controllability for semilinear waves in uniform time[END_REF] one has f y 0 ,y 1 κ,g ∈ C 0 ([0, T ], H 1 (M)) and the estimate f y 0 ,y 1 κ,g L ∞ (0,T ;H 1 (M))

(y 0 , y 1 ) H 2 (M)⊕H 1 (M) .

With this regularity of the source term in the right-hand-side of the wave equations in (1.8), one finds y, ỹ ∈ C 0 ([0, T ], H 2 (M)). Computing the difference in (1.8) one writes (6.7) P κ,g (y -ỹ) = (A κ,g -A κ,g )ỹ.

A hyperbolic energy estimate yields E κ,g (y -ỹ)(T ) 1/2 (A κ,g -A κ,g )ỹ L ∞ (0,T ;L 2 (M)) (κ, g) -(κ, g) X 1 ỹ L ∞ (0,T ;H 2 (M))

(κ, g) -(κ, g) X 1 f y 0 ,y 1 κ,g L ∞ (0,T ;H 1 (M))

(κ, g) -(κ, g) X 1 (y 0 , y 1 ) H 2 (M)⊕H 1 (M) .

In the case α = 0, one writes E κ,g (y -ỹ)(T ) 1/2 E κ,g (y)(T ) 1/2 + E κ,g (ỹ)(T ) 1/2 E κ,g (y)(T ) 1/2 + E κ,g (ỹ)(T ) 1/2 (y 0 , y 1 ) H 1 (M)⊕L 2 (M) .

Finally, the result follows from interpolation between the two cases α = 0 and α = 1.

1. 4 . 1 .

 41 Transport equation and measure support. Let O be an open set of a smooth manifold. We denote by 1 D (O) and 1 D ,0 (O) the spaces of density distributions and density Radon measures on O.

3. 1 .

 1 Microlocal defect density measures. We define M + (S * L) as the set of positive density measures on S * L. For µ ∈ M + (S * L) and a ∈ C 0 c (S * L), we shall write µ, a S * L = S * L a( )µ(d ),

,H 1 loc

 1 (Ω) + o(1) k→+∞ and (3.6) follows if we prove (3.5).

Proposition 3 . 7 .

 37 With the notation of Proposition 3.7, consider a sequence (u k ) k∈N ⊂ L 2 loc (L) that converges to 0 weakly and µ a L 2 -microlocal defect density measure associated with (u k ) k∈N .

4 .

 4 12), this concludes the proof of the second part of the proposition since {p 0 , b} = H p 0 b. Measure support propagation: proof of Theorem 1.10 Theorem 1.10 is stated on an open subset of a smooth manifold. Yet, its result is of local nature. Using a local chart we may assume that we consider an open of set Ω of R d instead without any loss of generality.

Proposition 4 . 1 .

 41 Let X be a C 0 -vector field on Ω an open set of R d . For a closed set F of Ω, the following two properties are equivalent.

  Iterative construction of the curve γn.

  Convergence of γn as n increases.

Figure 1 .

 1 Figure 1. Construction and convergence of the sequence (γ n ) n .

(4. 5 )

 5 γ n (s) = γ n (0) + s 0 γn (σ)dσ = x + s 0 X(γ n (σ))dσ + s 0 e n (σ) dσ.

Figure 2 .

 2 Figure 2. Estimation of the test function supports in the case ε ≤ 1.

K

  |f k κ,g ||∂ t ỹk |κµ g dt ≤ δ. with K = (0, T ) × ω \ (T 1 , T 2 ) × O . Let ψ ∈ C ∞ c ((0, T ) × ω) such that 0 ≤ ψ ≤ 1 and equal to 1 in a neighborhood of [T 1 , T 2 ] × O. One thus has (1 (0,T )×ω f k κ,g , ∂ t ỹk ) L 2 (L) ≤ (ψf k κ,g , ∂ t ỹk ) L 2 (L) + ((1 (0,T )×ω -ψ)f k κ,g , ∂ t ỹk ) L 2 (L)≤ (ψf k κ,g , ∂ t ỹk ) L 2 (L) + δ.

  Õj . Remark 3.1. Local properties of microlocal defect measures like µ can be deduced from the properties of χ C j j µ j . In what follows most results are of local nature. In such cases we shall work in local charts and use Sections 2.3 and 3.2 to bring the analysis to open domains of R d+1 . 3.3. Operators with a low regularity. An important tool we use to handle low regularity terms in what follows is a result due to R. Coifman and Y. Meyer (see

Observe that the change of variables made above for X(x 0 ) = (1, 0, . . . , 0) does not affect uniformity since the dilation is made by a factor in [cK , CK ].

Acknowledgements. This research was partially supported by Agence Nationale de la Recherche through projects ANA É ANR-13-BS01-0010-03, and ISDEEC ANR-16-CE40-0013, Institut universitaire de France (NB), and by the Tunisian Ministry for Higher Education and Scientific Research within the LR-99-ES20 program (BD). Finally, the authors wish to thank the anonymous reviewer for pointing out arguments that needed clarifications.