## STABILITY IN CONDUCTIVITY IMAGING FROM PARTIAL MEASUREMENTS OF ONE INTERIOR CURRENT

## CARLOS MONTALTO AND ALEXANDRU TAMASAN

Abstract. We prove a stability result in the hybrid inverse problem of recovering the electrical conductivity from partial knowledge of one current density field generated inside a body by an imposed boundary voltage. The region where interior data stably reconstructs the conductivity is well defined by a combination of the exact and perturbed data.

## 1. Introduction

We consider the problem of recovering the electrical conductivity of a body from partial interior measurements of one current density field generated by an imposed boundary voltage. Originally, one sought to recover the conductivity solely from boundary measurements. Most generally formulated by Calder´n [9], the problem has been extensively o studied: the unique determination results [5,8,30,31,35] are known to have weak (of logarithmic type) stability estimates [2], which cannot be, in general, improved [21]. The need for determining the conductivity by more robust methods has lead to considering new problems in the conductivity imaging, which employ some interior data. Such interior knowledge is obtained by coupling the electromagnetic model with the physics of some high resolution imaging method. For example, magnetic resonance measurements are employed in [18,19] to determine (one component) of the magnetic field inside, or in [13, 20, 40] to obtain the current density field inside, whereas ultrasound measurements are employed in [6, 10, 11] to obtain the interior knowledge of power densities; see [29,34,36,39] for some recent reviews.

This work concerns the problem of stability in determining the conductivity from partial knowledge of one current density inside. The interior data can be obtained from Magnetic Resonance measurements

2000 Mathematics Subject Classification. 35R30, 35J65, 65N21.

Key words and phrases. hybrid inverse problems, current density imaging, magnetic resonance, electrical impedance tomography, 1-Laplacian.

A. Tamasan was supported by the NSF Grant DMS 1312883.

as discovered in [33]. A first attempt to image the conductivity from the interior data of current density fields goes back to [40]. Since then several mathematical methods have been proposed, see [15, 19, 26-28] for the isotropic case, and [7,22,23] for the anisotropic case.

Of interest here, the work in [26, 27] shows that interior knowledge of the magnitude of one current density field generated by an imposed boundary voltage uniquely determine the conductivity inside. A general approach in [17] shows that using multiple measurements (two in the planar case) the linearized problem is stable, up to a finite dimensional kernel. Stability for the linear and non-linear problem from knowledge of the magnitude of one current density field was proved in [24]. The second author would like to point out a small omission in [24, Theorem 1.2], were the assumption that the voltage potentials must be equal at the boundary needs to be added. This is a natural assumption and it was clearly used in Equation (21) in there, but unfortunately omitted in the final version of the paper.

When the entire current density field is known, a stability result has been obtained in [16]. For the anisotropic conductivities we mention the work in [7] which a proved stability in two dimensions from full knowledge of four current density fields. All these stability results assume perturbations within the range of the data, locally near a given conductivity. If one assumes general perturbations in the interior data, the recent uniqueness result in [25] combined with the structural stability in [32] shows that the voltage potential can still be stably recovered inside. However, the class of smoothness is insufficient to further yield stability for the conductivity.

In the case of partial interior data, unique determination of planar conductivity have been showed to be possible in certain subdomains (see the injectivity region defined below) in [29], but the stability question has been left open. While the voltage potential can be stably recovered in certain subregions from partial interior data the stability estimates obtained in [37] do not guarantee the necessary regularity to extend to estimates for the conductivity.

In this paper we assume that knowledge of the current density field is only partially available inside. In fact we will show that knowledge of a certain component of the field suffices. We identify the specific subregions where conductivity can be stably recovered from partial interior and boundary data. The boundary voltage potential imposed to generate the current density field is assumed unperturbed, but need not be known everywhere. In general the stability region is a subset of the injectivity region as illustrated in Figure 2. However, if the accessible part of the boundary (where the voltage potential is known)

is simply connected, then the injectivity region and the stability region coincide under some geometric (strict convexity) assumptions on the shape of the boundary.

Different from the work in [24] we use a non-linear version of the decomposition of the data operator and identify a natural component of the current density field sufficient to yield stability. To authors' knowledge, this is the first result where either injectivity or stability is established from knowledge of just one component of the current density. The required component depends both on the measured and on the exact data.

Similar to all the works in the hybrids methods, we also work with voltages free of singular points inside the domain. In dimension two, the condition of no critical points can be satisfied under the assumption that the boundary illumination is almost two-to-one [1,26]. In dimensions n ≥ 3, it is still unclear if similar conditions exists; the recent results in [4, 14] describing the set of singular points assume nonzero frequencies.

## 2. Statement of the results

Let Ω ⊂ R n be a bounded domain C 2 ,α -diffeomorphic with the unit ball, 0 &lt; α &lt; 1. This assumption can be relaxed as explained in the proofs. Let σ be a positive function in C 2 ,α (Ω) and let f ∈ C 2 ,α ( ∂ Ω). From Schauder's regularity theory we know that the unique solution u of the Dirichlet problem

$$\left ( 1 \right ) \quad \nabla \cdot \sigma \nabla u = 0 \ \text{ in } \Omega, \quad \ \ u | _ { \partial \Omega } = f.$$

is in C 2 ,α (Ω), see e.g. [12, Theorem 6.18]. We refer to such a u as being σ -harmonic. The corresponding current density vector field J : C 1 ,α (Ω) → C 1 ,α ( ¯ ) Ω is defined by

$$\mathbf J ( \sigma ) = - \sigma \nabla u.$$

In order to state our partial data results we introduce the following notation. For some u ∈ C 1 (Ω) arbitrary (not necessarily σ -harmonic) with |∇ u | &gt; 0 in Ω, and p an arbitrary point in Ω, let γ p denote the integral curve starting at p in the direction of ∇ u . Since |∇ u | &gt; 0, there exist a first time, denoted by t + p (resp. t -p ), such that the integral curve starting at p and moving in the direction ±∇ u hits the boundary. We denote by

$$l _ { p } ^ { \pm } = \{ \gamma _ { p } ( t ) \, \colon t \in [ 0, t _ { p } ^ { \pm } ] \}$$

the segment of the integral curve from t = 0 to t = t ± p . Also, let

$$\Sigma _ { p } = \{ x \in \overline { \Omega } \, \colon u ( x ) = u ( p ) \}$$

denote the level curve of u passing through p , see Figure 1.

Figure 1. Illustration of the segment l p and the surface Σ p .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0423v2-figure-4f064b92-0b37-4d72-a7a4-b2b74faf56f9.jpg)

Let Γ be an open subset of the boundary to denote the accessible part. The following definitions specify the interior subdomains, where the unique and stable determination of the conductivity can be guaranteed.

Definition 2.1. Let Ω ⊂ R n be a bounded domain C 2 ,α -diffeomorphic to the unit ball and let u ∈ C 2 ,α (Ω) , where 0 &lt; α &lt; 1 .

Injectivity region: We say that a point p ∈ Ω is visible from Γ along the equipotential set if

$$\Sigma _ { p } \cap \partial \Omega \subset \Gamma.$$

The set I (Γ , u ) of points in Ω that are visible from Γ along equipontential sets is called the injectivity region ,

$$\mathcal { I } ( \Gamma, u ) \, \coloneqq \, \{ p \in \overline { \Omega } \, \colon \Sigma _ { p } \cap \partial \Omega \subset \Gamma \}.$$

Stability region: We say that the trajectory through p along ∇ u is visible from Γ if either l + p ⊂ I (Γ , u ) or l -p ⊂ I (Γ , u ) . The set S (Γ , u ) of points in Ω for which the corresponding trayectories along ∇ u are visible from Γ is called the stability region , i.e.,

$$( 4 ) \quad \mathcal { S } ( \Gamma, u ) \coloneqq \{ p \in \mathcal { I } ( \Gamma, u ) \colon l _ { p } ^ { + } \subset \mathcal { I } ( \Gamma, u ) \text{ or } l _ { p } ^ { - } \subset \mathcal { I } ( \Gamma, u ) \}.$$

Clearly S (Γ , u ) ⊆ I (Γ , u ) and, if Γ is connected, the equality can hold as illustrated in Figure 3. In Figure 2 we illustrate the injectivity and stability regions, when Γ is not connected.

Figure 2. The injectivity region I (Γ , u ) is the light grey region that contains the stability region S (Γ , u ) in dark grey.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0423v2-figure-c87c2bed-3e67-4a08-9ad2-baf79d3b44df.jpg)

Figure 3. When Γ is connected the visible region and the trajectory region can be the same.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0423v2-figure-7dc53b88-c117-48f7-816a-4dd7329be7f2.jpg)

Let ˜ be the voltage potential corresponding to some perturbed conu ductivity ˜ σ ∈ C 1 ,α (Ω),

$$\nabla \cdot \tilde { \sigma } \nabla \tilde { u } = 0 \quad \text{in } \Omega, \quad \ \tilde { u } | _ { \partial \Omega } = f,$$

and let J (˜) := σ - ∇ ˜ σ ˜ be the corresponding current density field. u Let

$$\delta \sigma \coloneqq \sigma - \tilde { \sigma } \quad \text{and} \quad \delta \mathbf J \coloneqq \mathbf J - \tilde { \mathbf J }$$

denote the corresponding perturbations.

We prove that δσ is controlled by the component of δ J in the direction ∇ ( u + ˜) throughout u S (Γ , u + ˜). u To the author's knowledge the result below is the first stability result from partial interior data. It is

also the first time the recovery is done from just one component of the current density field.

For a vector field w 0 in Ω we denote the orthogonal projection onto w 0 by Π w 0 , and onto the orthogonal complement (both in the Euclidean metric) by Π ⊥ w 0 := Id -Π w 0 the orthogonal projection; i.e., for an arbitrary vector field w ,

$$\Pi _ { w _ { 0 } } \mathbf w = \frac { \mathbf w _ { 0 } \cdot \mathbf w } { | \mathbf w _ { 0 } | ^ { 2 } } \mathbf w _ { 0 } \quad \text{and} \quad \Pi _ { w _ { 0 } } ^ { \perp } \mathbf w = \mathbf w - \Pi _ { w _ { 0 } } \mathbf w.$$

Figure 4. From knowledge to the scalar component of w = Π ∇ ( u +˜) u ( δ J ) we can stably recover the difference in the conductivities σ -˜. σ Here δ J represents the difference J ( σ ) -J (˜). σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0423v2-figure-1518e6aa-d753-430b-8abb-a9fdc12420c5.jpg)

Theorem 2.1 (Stability with Partial Data) . Let Ω ⊂ R n be C 2 ,α -diffeomorphic to the unit ball, and σ, σ ˜ ∈ C 1 ,α (Ω) be positive on Ω , for some 0 &lt; α &lt; 1 . Let u (respectively ˜ u ) in C 2 ,α (Ω) be σ (respectively ˜ σ )-harmonic. Let Γ ⊆ ∂ Ω be the union of finitely many open connected components, and Γ ′ ⊂⊂ Γ be a open subset compactly contained in Γ . Assume that

$$\sigma | _ { \Gamma } = \tilde { \sigma } | _ { \Gamma }, \, u | _ { \Gamma } = \tilde { u } | _ { \Gamma },$$

and

$$| \nabla ( u + \tilde { u } ) | > 0, \quad i n \ \overline { \mathcal { I } ( \Gamma ^ { \prime }, u + \tilde { u } ) }.$$

Then there exist C &gt; 0 dependent on a lower bound ot |∇ ( u + ˜) u | in I (Γ ′ , u + ˜) u , the domain Ω , and the C 1 (Ω) - norm of σ and ˜ σ , such that (10)

$$\stackrel { \dots } { \left \| \sigma - \tilde { \sigma } \right \| } _ { L ^ { 2 } ( \mathcal { S } ( \Gamma ^ { \prime }, u + \tilde { u } ) ) } \leq C \left | \left | \nabla \cdot \left ( \Pi _ { \nabla ( u + \tilde { u } ) } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \right ) \right | \right | _ { L ^ { 2 } ( \mathcal { I } ( \Gamma ^ { \prime }, u + \tilde { u } ) ) } ^ { \frac { \alpha } { 2 + \alpha } }.$$

In particular,

$$( 1 1 ) \ \| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \mathcal { S } ( \Gamma ^ { \prime }, u + \tilde { u } ) ) } \leq C \left | \left | \Pi _ { \nabla ( u + \tilde { u } ) } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \right | \right | _ { H ^ { 1 } ( \mathcal { I } ( \Gamma ^ { \prime }, u + \tilde { u } ) ) } ^ { \frac { \alpha } { 2 + \alpha } }.$$

The proof of the Theorem 2.1 is presented in Section 4.

A few remarks are in order. Note first, that the estimates (10) and (11) concern norms over (possibly) different domains. However, in many interesting situation S (Γ ′ , u + ˜) = u I (Γ ′ , u + ˜) as illustrated in u Figure 3. In particular, in the full data case with Γ = ∂ Ω, we have S (Γ ′ , u + ˜) = u I (Γ ′ , u + ˜) = Ω, and the estimate (11) yields u

$$\| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) } \leq C \| \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) \| | _ { H ^ { 1 } ( \Omega ) } ^ { \frac { \alpha } { 2 + \alpha } }.$$

The stability estimate requires knowledge of a component of J in a direction ∇ ( u +˜) that is well defined by the exact and perturbed data u due to the uniqueness results in [27]. However, if ˜ is σ a priori close to σ then the level curves of ˜ will be close to the level curves u u , and it will suffice to project δ J onto the ∇ u with a penalty term that will depend on the apriori closeness assumption, as illustrated in Figure 5. More precisely, we have the following local stability result.

Figure 5. We illustrate how we can controlled the visible region and the projection of the current density. The underlying idea is to chose voltages potential f at the boundary to induced specific level curves u 0 =const. By doing so, we get a better understanding of the required direction of the current density to obtain an stable reconstruction. Here, δ J 0 denote J ( σ ) -J ( σ 0 ) and w 0 = Π ∇ ( u 0 ) ( δ J 0 ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0423v2-figure-f5bb2897-4b96-4361-bc75-ccaac57a051f.jpg)

Theorem 2.2. Let σ ∈ C 1 ,α (Ω) , 0 &lt; α &lt; 1 , be positive in Ω and u be σ -harmonic with |∇ u | &gt; 0 in Ω . There exists an glyph[epsilon1] &gt; 0 depending on Ω and some C &gt; 0 depending on glyph[epsilon1] , such that the following holds: If ˜ σ ∈ C 1 ,α (Ω) with

$$\| \sigma - \tilde { \sigma } \| _ { C ^ { 1, \alpha } ( \overline { \Omega } ) } < \epsilon \ \ a n d \ \sigma | _ { \partial \Omega } = \tilde { \sigma } | _ { \partial \Omega },$$

and ˜ u is the ˜ σ -harmonic map with

$$\tilde { u } | _ { \partial \Omega } = u | _ { \partial \Omega },$$

then

$$\| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) } \leq C \, \| \nabla \cdot \left ( \Pi _ { \nabla u } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \| _ { L ^ { 2 } ( \Omega ) } ^ { \frac { \alpha } { 2 + \alpha } } \,.$$

In particular,

$$\| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) } \leq C \, | | \Pi _ { \nabla u } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) | | _ { H ^ { 1 } ( \Omega ) } ^ { \frac { \alpha } { 2 + \alpha } } \,.$$

As a consequence of Theorem 2.2 let us consider the problem of determining a perturbation from a constant conductivity: If the boundary voltage is uniform in the z -direction, the equipotential surfaces are planes perpendicular to the z -direction and we only need two components of the magnetic field B = ( B , B , B x y z ) to stably recover the conductivity, as exemplified below.

Example 1 . For simplicity assume that σ ≡ 1 and f ( x, y, z ) = z . Then u x, y, z ( ) = z and, by Amp` ere's law,

$$\delta J ^ { z } = \nabla u \cdot \delta \mathbf J = \frac { 1 } { \mu _ { 0 } } \nabla u \cdot ( \nabla \times \delta \mathbf B ) = \frac { 1 } { \mu _ { 0 } } \left ( \frac { \partial } { \partial x } \delta B ^ { y } - \frac { \partial } { \partial y } \delta B ^ { x } \right ).$$

By Theorem 2.2, we can stably recover a sufficiently small pertubation of the conductivity from δB x and δB y .

## 3. Preliminaries

The stability result relies on the following decomposition of the perturbation δ J in the data.

glyph[negationslash]

Proposition 3.1. Let u and ˜ u be σ -harmonic and ˜ σ -harmonic, respectively. If ∇ ( u + ˜) = 0 u in V , for some open subset V of Ω , then

$$2 \nabla \cdot \Pi _ { \nabla ( u + \tilde { u } ) } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) = L ( u - \tilde { u } ) \quad i n \ V.$$

where L is the differential operator defined by

$$( 1 7 ) \ L v \coloneqq - \nabla \cdot ( \sigma + \tilde { \sigma } ) \nabla v + \nabla \cdot \left ( ( \sigma + \tilde { \sigma } ) \frac { \nabla ( u + \tilde { u } ) \cdot \nabla v } { | \nabla ( u + \tilde { u } ) | ^ { 2 } } \nabla ( u + \tilde { u } ) \right ).$$

Proof. Consider the difference of two internal measurements

$$( 1 8 ) \quad \ 2 ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) = ( \sigma - \tilde { \sigma } ) \nabla ( u + \tilde { u } ) + ( \sigma + \tilde { \sigma } ) \nabla ( u - \tilde { u } ).$$

If we dot product with ∇ ( u + ˜) we obtain u

(19)

$$\dot { 2 } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \cdot \nabla ( u + \tilde { u } ) = ( \sigma - \tilde { \sigma } ) | \nabla ( u + \tilde { u } ) | ^ { 2 } + ( \sigma + \tilde { \sigma } ) \nabla ( u - \tilde { u } ) \cdot \nabla ( u + \tilde { u } ).$$

Solving for σ -˜ we have σ

$$( 2 0 ) \ \sigma - \tilde { \sigma } = \frac { 2 ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \cdot \nabla ( u + \tilde { u } ) } { | \nabla ( u + \tilde { u } ) | ^ { 2 } } - ( \sigma + \tilde { \sigma } ) \frac { \nabla ( u + \tilde { u } ) \cdot \nabla ( u - \tilde { u } ) } { | \nabla ( u + \tilde { u } ) | ^ { 2 } }$$

On the other hand it follows easily from (1)

$$( 2 1 ) \quad \nabla \cdot ( \sigma + \tilde { \sigma } ) \nabla ( u - \tilde { u } ) = - \nabla \cdot ( \sigma - \tilde { \sigma } ) \nabla ( u + \tilde { u } ) \ \text{ in } \Omega.$$

Using (20) we substitute σ -˜ in (21) and we get σ

$$\text{Using} \, ( 2 0 ) \, \text{we substitute} \, \sigma - \tilde { \sigma } \, \text{ in} \, ( 2 1 ) \, \text{and} \, \text{we get} \\ L ( u - \tilde { u } ) & = - \nabla \cdot ( \sigma + \tilde { \sigma } ) \nabla ( u - \tilde { u } ) \\ & \quad + \nabla \cdot ( \sigma + \tilde { \sigma } ) \frac { \nabla ( u + \tilde { u } ) \cdot \nabla ( u - \tilde { u } ) } { | \nabla ( u + \tilde { u } ) | ^ { 2 } } \nabla ( u + \tilde { u } ) \\ & = \nabla \cdot \left ( \frac { 2 ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \cdot \nabla ( u + \tilde { u } ) } { | \nabla ( u + \tilde { u } ) | ^ { 2 } } \nabla ( u + \tilde { u } ) \right ) \\ & = 2 \nabla \cdot \Pi _ { \nabla ( u + \tilde { u } ) } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \\ \quad \text{We remark that} \, ( 1 6 ) \, \text{is the non-linear version of the decomposition}$$

We remark that (16) is the non-linear version of the decomposition obtained in the linearized case in [24]. The representation of the operator L in local coordinates can then be obtained similarly as in [24]:

glyph[negationslash]

Proposition 3.2. Let u, u ˜ ∈ C 2 (Ω) , with ∇ ( u + ˜)( u x 0 ) = 0 for x 0 ∈ Ω . Then locally near Σ x 0 , the operator L in (17) is the restriction of ∇ ( σ + ˜) σ ∇ onto the level surfaces ( u + ˜) = u const. Moreover, if y n = u + ˜ u , and y ′ = ( y , . . . , y 1 n -1 ) are any local coordinates of the level set Σ x 0 , then in a neighborhood of Σ x 0 the Euclidean line element d s is given by

$$( 2 2 ) \quad \text{d} s ^ { 2 } = c ^ { 2 } ( \text{d} y ^ { n } ) ^ { 2 } + g _ { \alpha \beta } \text{d} y ^ { \alpha } \text{d} y ^ { \beta }, \ \ g _ { \alpha, \beta } \coloneqq \sum _ { i } ^ { n - 1 } \frac { \partial x ^ { i } } { \partial y ^ { \alpha } } \frac { \partial x ^ { i } } { \partial y ^ { \beta } },$$

where c = |∇ ( u + ˜) u | -1 . In this coordinates the operator L becomes

$$L = - \sum _ { \alpha, \beta } \frac { 1 } { \sqrt { \det g } } \frac { \partial } { \partial y ^ { \beta } } ( \sigma + \tilde { \sigma } ) g ^ { \alpha \beta } \sqrt { \det g } \frac { \partial } { \partial y ^ { \alpha } },$$

where the Greek super/subscripts run from 1 to n -1 .

Proof. Let Σ x 0 be the level set of u + ˜ passing through u x 0 . On Σ x 0 choose local coordinates y ′ = ( y , . . . , y 1 n -1 ) and set y n = ( u + ˜)( u x ). Then y = ( y , y ′ n ) = ϕ x ( ) define local coordinates in a neighborgood of Σ x 0 . In the new coordinates, the Euclidean metric becomes

$$( 2 4 ) \quad \text{d} s ^ { 2 } = c ^ { 2 } ( \text{d} y ^ { n } ) ^ { 2 } + g _ { \alpha \beta } \text{d} y ^ { \alpha } \text{d} y ^ { \beta }, \text{ where } g _ { \alpha \beta } \coloneqq \sum _ { i = 1 } ^ { n - 1 } \frac { \partial x ^ { i } } { \partial y ^ { \alpha } } \frac { \partial x ^ { i } } { \partial y ^ { \beta } },$$

where c = |∇ ( u + ˜) u | -1 . This follows easily by noticing that u + ˜ u trivially solves the eikonal equation c 2 |∇ φ | 2 = 1 for the speed c = |∇ ( u + ˜) u | -1 . Then, near x 0 , u + ˜ is the signed distance from the u x to the level surface Σ x 0 and the Euclidean metric has block structure given by (24). Denote by g be metric after the change of variables, i.e.,

$$g = \left [ \begin{array} { c c } [ g _ { \alpha, \beta } ] & 0 \\ 0 & c ^ { 2 } \end{array} \right ]$$

Let φ ∈ C ∞ 0 (Ω) be a test function supported near some ˆ x 0 ∈ Σ x 0 . We have

$$\langle L v, \phi \rangle & = \langle ( \sigma + \tilde { \sigma } ) \nabla v, \nabla \phi \rangle - \langle ( \sigma + \tilde { \sigma } ) \Pi _ { \nabla ( u + \tilde { u } ) } \nabla v, \nabla \phi \rangle, \\ & = \langle ( \sigma + \tilde { \sigma } ) \Pi _ { \nabla ( u + \tilde { u } ) } ^ { \perp } \nabla v, \Pi _ { \nabla ( u + \tilde { u } ) } ^ { \perp } \nabla \phi \rangle \,,$$

where 〈· , ·〉 denotes the scalar product in L 2

(Ω).

From (24), we get that for any function v ∈ C 1 (Ω)

$$\Pi _ { \nabla ( u + \tilde { u } ) } \nabla _ { x } v = \left ( 0, \dots, 0, c ^ { 2 } \frac { \partial v } { \partial y ^ { n } } \right ).$$

and

$$\Pi ^ { \perp } _ { \nabla ( u + \tilde { u } ) } \nabla _ { x } v = \left ( \frac { \partial v } { \partial y ^ { 1 } }, \dots, \frac { \partial v } { \partial y ^ { n - 1 } }, 0 \right ).$$

thus in ( y , y ′ n ) coordinates (25) becomes

(26)

$$\langle L v, \phi \rangle & = \sum _ { \alpha, \beta } \int _ { \varphi ^ { - 1 } ( \Omega ) } ( \sigma + \tilde { \sigma } ) \left ( g ^ { \alpha \beta } \frac { \partial v } { \partial y ^ { \alpha } } \frac { \partial \phi } { \partial y ^ { \beta } } \right ) \sqrt { \det g } \, d y \\ & = - \sum _ { \alpha, \beta } \int _ { \varphi ^ { - 1 } ( \Omega ) } \frac { 1 } { \sqrt { \det g } } \left ( \frac { \partial } { \partial y ^ { \beta } } ( \sigma + \tilde { \sigma } ) g ^ { \alpha \beta } \sqrt { \det g } \frac { \partial } { \partial y ^ { \alpha } } \right ) \phi ( y ) d y \\ \circ \dots \dots \lambda \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \sim \lambda \not \, \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

Since φ and ˆ x 0 ∈ Σ x 0 were arbitrarily we conclude that

$$L = - \sum _ { \alpha, \beta } \frac { 1 } { \sqrt { \det g } } \left ( \frac { \partial } { \partial y ^ { \beta } } ( \sigma + \tilde { \sigma } ) g ^ { \alpha \beta } \sqrt { \det g } \frac { \partial } { \partial y ^ { \alpha } } \right ).$$

near Σ . 0 glyph[square]

We note that Proposition 3.2 shows L to be an elliptic differential operator acting onto the level sets of u +˜, for which the normal comu ponent becomes a parameter.

## 4. Stability estimates for partial data

In this section we prove Theorem 2.1 and its corollary. The proof is based on establishing separate estimates for u -˜ and for the operator u L in (16). For points in S (Γ ′ , u + ˜) the visibility from Γ u ′ parallel to ∇ ( u + ˜) will allow us to control u u -˜, while the the visibility along u the equipotential sets will be used for estimates on L .

Recall that S (Γ ′ , u + ˜) = u S + (Γ ′ , u + ˜) u ∪ S -(Γ ′ , u + ˜), where u

$$\mathcal { S } ^ { \pm } ( \Gamma ^ { \prime }, u ) = \{ p \in \overline { \Omega } \, \colon \Sigma _ { p } \cap \partial \Omega \subset \Gamma ^ { \prime } \text{ and } l _ { p } ^ { \pm } \subset \mathcal { I } ( \Gamma, u ) \}.$$

We will show that

(28)

$$\stackrel { \llcorner } { \sigma } - \tilde { \sigma } \| _ { L ^ { 2 } ( \mathcal { S } ^ { \pm } ( \Gamma ^ { \prime }, u + \tilde { u } ) ) } \leq C \left | \left | \nabla \cdot \left ( \Pi _ { \nabla ( u + \tilde { u } ) } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \right ) \right | ^ { \frac { \alpha } { 2 + \alpha } } _ { L ^ { 2 } ( \mathcal { I } ( \Gamma ^ { \prime }, u + \tilde { u } ) ) }.$$

We prove (28) for S + (Γ ′ , u + ˜). u The corresponding inequality for S -(Γ ′ , u + ˜) follows similarly. u

To simplify notation, let

$$( 2 9 ) \ \mathcal { S } \coloneqq \mathcal { S } ^ { + } ( \Gamma, u + \tilde { u } ), \quad \mathcal { S } ^ { \prime } \coloneqq \mathcal { S } ^ { + } ( \Gamma ^ { \prime }, u + \tilde { u } ), \text{and} \ \mathcal { I } ^ { \prime } \coloneqq \mathcal { I } ( \Gamma ^ { \prime }, u + \tilde { u } ).$$

We assume w.l.o.g. that S ′ and S are connected. Notice that, since there are only finitely many components, we can add the estimates that we obtain for the single-component case, to the more general case.

We will make use of the following two technical estimates in the Lemmas below.

Lemma 4.1. There exist C &gt; 0 dependent on the domain Ω such that

$$\| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \mathcal { S } ^ { \prime } ) } \leq C \| u - \tilde { u } \| _ { H ^ { 2 } ( \mathcal { S } ^ { \prime } ) },$$

where for simplicity we denote S ′ := S + (Γ ′ , u + ˜) u .

glyph[negationslash]

Proof. Since ∇ ( u + ˜) = 0 in Ω, without loss of generality (otherwise u work with a diffeormorphic image of Ω, and, possibly, with -( u + ˜)), u we may assume that

$$\min _ { \overline { \Omega } } \frac { \partial ( u + \tilde { u } ) } { \partial x ^ { n } } > 0.$$

Since the n -th coordinate now plays a special role, to simplify notation, let x ′ := ( x , . . . , x 1 n -1 ), so that x = ( x , x ′ n ).

Let Z be the set (possibly empty) of boundary points , where ∇ ( u +˜) u is tangential to the boundary,

$$Z \coloneqq \{ x \in \partial \Omega \, \colon \, n ( x ) \cdot \nabla ( u + \tilde { u } ) ( x ) = 0 \}.$$

The regularity assumptions on the boundary of the domain, and on u, u ˜ yield that Z is closed and confined to a co-dimension 1 variety of the boundary, in particular, Z is negligible with respect to the induced area measure on the boundary.

Let Γ ′ ⊂⊂ Γ \ Z , x 0 ∈ S ′ = S + (Γ ′ , u + ˜), u and Σ 0 := { x ∈ Ω : ( u + ˜)( u x ) = ( u + ˜)( u x 0 ) } be the level set passing through x 0 .

glyph[negationslash]

↦

Since ∂ u ( +˜) u ∂x n ( x , x ′ 0 n 0 ) = 0, by the implicit function theorem, in a neighborhood O ⊂ R n of x ′ 0 there is a local representation of Σ : 0 O glyph[owner] x ′ → ( x , g ′ ( x ′ )), for some C 1 ,α -smooth coordinate map g : O ⊂ R n -1 → R with g x ( ′ 0 ) = x n 0 .

We define next a tubular neighborhood T x 0 by flowing along ∇ ( u +˜) u the points on O ⊂ Σ 0 until the boundary is met. More precisely, we consider the family (indexed in x ′ ∈ O ) of the initial value problems

$$( 3 3 ) \quad \frac { d } { d t } \gamma ( t ; x ^ { \prime } ) = \frac { \nabla ( u + \tilde { u } ) } { | \nabla ( u + \tilde { u } ) | } ( \gamma ( t ; x ^ { \prime } ) ), \ \ \gamma ( 0 ) = ( x ^ { \prime }, g ( x ^ { \prime } ) ), \ x ^ { \prime } \in O.$$

Since x 0 ∈ S ′ , for the solution γ ( ; · x ′ 0 ) of (33) with parameter x ′ = x ′ 0 meets the boundary transversally, hence there exists a smallest time t x ′ 0 such that γ t ( x ′ 0 ; x ′ 0 ) ∈ Γ , and ′ γ t ( ; x ′ 0 ) ∈ Ω for 0 ≤ t &lt; t x ′ 0 .

↦

The continuous dependence with parameters of solutions of initial value problems for ODE's shows that, by possibly shrinking the neighborhood O , we have that t → γ t, x ( ′ ) are defined on a maximal interval with γ ( · , x ′ ) : [0 , t x ′ ) → Ω and γ t ( x ′ ; x ′ ) ∈ Γ . ′

We consider the following change of coordinates ( x , x ′ n ) → ( y , y ′ n ), where y ′ = x ′ and x n = γ y ( n ; y ′ ). In the new coordinates, the tubular neighborhood T x 0 rewrites

$$T _ { x _ { 0 } } = \{ ( y ^ { \prime }, y ^ { n } ) \in \mathbb { R } ^ { n } \, \colon y ^ { \prime } \in O, \, 0 \leq y ^ { n } \leq t _ { y ^ { \prime } } \}.$$

Let J y , y ( ′ n ) denotes the absolute value of the Jacobian determinant of the change of coordinates above at some ( y , y ′ n ) ∈ T x 0 , and let β &gt; 0 denote an upper bound on the quotient

$$\beta \coloneqq \frac { \max _ { T _ { x _ { 0 } } } J } { \min _ { T _ { x _ { 0 } } } J } < \infty$$

Recall the equation (1) interpreted as a transport equation for δσ = ( σ -˜): σ

$$( 3 5 ) \quad \ \nabla ( u + \tilde { u } ) \cdot \nabla ( \delta \sigma ) + ( \delta \sigma ) \Delta ( u + \tilde { u } ) = G \ \text{ in } \Omega,$$

where, for brevity, we let

$$G \coloneqq - \nabla \cdot ( \sigma + \tilde { \sigma } ) \nabla ( u - \tilde { u } ).$$

In the local coordinates in the interior of T x 0 , the equation (35) becomes

$$\frac { 1 } { c ^ { 2 } } \frac { \partial ( \delta \sigma ) } { \partial y ^ { n } } + ( \delta \sigma ) \Delta ( u + \tilde { u } ) = G,$$

where c = |∇ ( u + ˜) u | -1 . By (9),

$$\max _ { \overline { \Omega } } c < \infty.$$

By hypothesis (13) we also have that

$$( \delta \sigma ) ( y ^ { \prime }, t _ { y ^ { \prime } } ) = 0, \quad \forall y ^ { \prime } \in O.$$

Using an integral factor we can solve (37) and (39) to get

$$( 4 0 ) \quad \ \ ( \delta \sigma ) ( y ^ { \prime }, y ^ { n } ) = \left ( \frac { 1 } { \mu ( y ^ { \prime }, y ^ { n } ) } \int _ { t _ { y ^ { \prime } } } ^ { y ^ { n } } c ^ { 2 } G ( y ^ { \prime }, t ) \mu ( y ^ { \prime }, t ) d t \right ),$$

where

$$\mu ( y ^ { \prime }, y ^ { n } ) = \exp \left ( - \int _ { y ^ { n } } ^ { t _ { y ^ { \prime } } } c ^ { 2 } \Delta ( u + \tilde { u } ) ( y ^ { \prime }, t ) d t \right ).$$

Let m,M (not necessarily positive) be such that

$$m \leq c ^ { 2 } \Delta ( u + \tilde { u } ) ( x ) \leq M, \quad \forall x \in \overline { \Omega }.$$

Then, it is easy to see that, for y ′ ∈ O and 0 ≤ t ≤ y n ≤ t y ′ , we have

$$e ^ { - | M | d i a m ( \Omega ) } \leq \frac { \mu ( y ^ { \prime }, t ) } { \mu ( y ^ { \prime }, y ^ { n } ) } \leq e ^ { | m | d i a m ( \Omega ) }.$$

For a constant ˜ C &gt; 0 depending on the bounds in (42) and (38), and β in (34), we estimate

(43)

$$\beta \text{ in (34), we estimate } \cdot \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \
 \| \delta \sigma \| _ { L ^ { 2 } ( T _ { x _ { 0 } } ) } & = \int _ { y ^ { \prime } \in O } \int _ { 0 } ^ { t _ { y ^ { \prime } } } \left | \int _ { 0 } ^ { y ^ { n } } c ^ { 2 } G ( y ^ { \prime }, t ) \frac { \mu ( y ^ { \prime }, t ) } { \mu ( y ^ { \prime }, y ^ { n } ) } d t \right | ^ { 2 } J ( y ^ { \prime }, y ^ { n } ) d y ^ { n } d y ^ { \prime } \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \	 \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \cdots \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \end{array} \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ = \, \, \,$$

Since S ′ can be covered by finitely many tubular neighborhoods, and the norm in the right hand side of (43), we conclude that

$$( 4 4 ) \quad \| \delta \sigma \| _ { L ^ { 2 } ( \mathcal { S } ^ { \prime } ) } ^ { 2 } \leq C \| \nabla \cdot ( \sigma + \tilde { \sigma } ) \nabla ( u - \tilde { u } ) \| _ { L ^ { 2 } ( \mathcal { S } ^ { \prime } ) } \leq C \| u - \tilde { u } \| _ { H ^ { 2 } ( \mathcal { S } ^ { \prime } ) },$$

for a constant C depending on u, u ˜ and the domain Ω, and a priori bounds on C 1 -norm of σ and ˜. This finishes the proof of (30). σ

glyph[square]

Lemma 4.2. There exists C &gt; 0 dependent on Ω and a lower bound on σ + ˜ σ , such that,

$$\| u - \tilde { u } \| _ { L ^ { 2 } ( \mathcal { I } ^ { \prime } ) } \leq C \| L ( u - \tilde { u } ) \| _ { L ^ { 2 } ( \mathcal { I } ^ { \prime } ) },$$

where, for simplicity we denote by I ′ := I (Γ ′ , u + ˜) u .

Proof. Let x 0 be arbitrarily fixed in the injectivity region I ′ . Assume first that the level set Σ 0 of u + ˜ passing through u x 0 intersects the boundary in Γ , transversally at any point of intersection. ′ Consider the local coordinates ( y , y ′ n ) for the equipotential set Σ , introduced in 0 Proposition 3.2. By continuity of y n = u + ˜ w.r.t. u y ′ , for sufficiently small glyph[epsilon1] , the neighborhood U glyph[epsilon1] of Σ 0 defined by

$$U _ { \epsilon } \coloneqq \{ ( y ^ { \prime }, y ^ { n } ) \in \mathcal { I } ^ { \prime } \colon y ^ { \prime } \in \Sigma _ { 0 }, \text{and} \ y ^ { n } \in ( - \epsilon, \epsilon ) \}$$

is visible from Γ ′ along the equipotential sets of u + ˜. u

Using the local representation of L obtained in (***) and local representation the metric g on U glyph[epsilon1] we estimate L δu ( ). We denote by g ′ the restriction of the metric to y ′ , i.e., g ′ = [ g α,β ] for 1 ≤ α, β ≤ n -1. Let δu = u -˜. u Note that in U glyph[epsilon1] , ( u -˜) u | ∂ Ω = 0, hence Poincar´'s inequality e is valid in Σ . 0 We estimate

$$\begin{array} { l } \text{is valid in $\Sigma_{0}$. We estimate} \\ \\ \langle L \delta u \cdot \delta u \rangle | _ { U _ { \epsilon } } = & \int _ { - \epsilon } ^ { \epsilon } \left ( \frac { 1 } { \sqrt { \det g } } \frac { \partial \delta u } { \partial y ^ { \beta } } ( \sigma + \vec { \sigma } ) g ^ { \alpha \beta } \sqrt { \det g } \frac { \partial \delta u } { \partial y ^ { \alpha } } \right ) \cdot \delta u \sqrt { \det g } \, d y ^ { \prime } d y ^ { n } \\ \\ = & \int _ { - \epsilon } ^ { \epsilon } \int ( \sigma + \vec { \sigma } ) g ^ { \alpha \beta } \frac { \partial \delta u } { \partial y ^ { \alpha } } \frac { \partial \delta u } { \partial y ^ { \beta } } \sqrt { \det g } \, d y ^ { n } \\ \\ \geq \delta \int _ { - \epsilon } ^ { \epsilon } \int _ { \Sigma _ { 0 } } | \nabla _ { y } \delta u | ^ { 2 } \, \ d y ^ { \prime } d y ^ { n }, \\ \\ \geq \delta \int _ { - \epsilon } ^ { \epsilon } \int _ { \Sigma _ { 0 } } | \delta u | ^ { 2 } \, \ d y ^ { \prime } d y ^ { n } \\ \\ ( 4 6 ) & = \delta \| \delta u \| _ { L ^ { 2 } ( U _ { \epsilon } ) }, \\ \text{where} \, \delta > 0 \, \text{depends on an a priori lower bound of $\sigma + \vec { \sigma }$ and of $g^{\alpha,\beta}$ on} \\ \Omega \text{ for $1\leq\alpha$, $\beta\leq n$. The second inequality is the Poincaré inequality.} \end{array}$$

where δ &gt; 0 depends on an a priori lower bound of σ +˜ and of σ g α,β on Ω for 1 ≤ α, β ≤ n . The second inequality is the Poincar´ inequality. e Since U glyph[epsilon1] ⊂ I ′ we obtain from (46)

$$\| \delta u \| _ { L ^ { 2 } ( U _ { \epsilon } ) } ^ { 2 } \leq C \langle L ( \delta u ), \delta u \rangle _ { L ^ { 2 } ( \mathcal { I } ^ { \prime } ) },$$

By compactness, we can get a finite covering of I ′ with neighborhoods U glyph[epsilon1] . Summing (47) over the covering of U glyph[epsilon1] we obtain

$$\| u - \tilde { u } \| _ { L ^ { 2 } ( \mathcal { I } ^ { \prime } ) } ^ { 2 } \leq C ( L ( u - \tilde { u } ), u - \tilde { u } ) _ { L ^ { 2 } ( \mathcal { I } ^ { \prime } ) },$$

for some constant C which depends on Ω, and the lower bound on σ + ˜. Finally from Cauchy-Schwartz inequality (45) follows. σ

Finally the level set Σ 0 intersects tangentially the boundary of Ω as some point then the analysis presented above could fail depending on the degeneracy of such interseccion. In particular it is not clear how to justify integration by parts on U glyph[epsilon1] in the intersection is degenerate. To avoid this difficulty we will extend σ, σ, u ˜ and ˜ in two different u ways to an open domain Ω 1 containing I ′ . The first extension takes advantage of the boundary information while the second extension uses the regularity of the functions to extend the geometry and the operators in decomposition (16).

Since Γ ′ is compactly supported in Γ it follows that I ′ ⊂ I . Thus there exists an open bounded Ω 1 with smooth boundary, such that I ′ ⊂ Ω 1 and Ω 1 ∩ Ω ⊂ I . First, using the fact that ( σ -˜) σ | Γ = 0 we

can find H 1 (Ω ) extensions 1 σ 1 and ˜ σ 1 of σ and ˜, respectively, σ such that ( σ 1 -˜ ) = 0 in Ω σ 1 1 \ Ω; (e.g., extend first σ to some σ 1 ∈ C 2 (Ω ), 1 and define ˜ σ 1 := σ 1 in Ω 1 \ Ω. Since ˜ σ 1 coincide with σ 1 on ∂ Ω, the difference ˜ σ 1 -σ 1 , and thus ˜ σ 1 lie in H 1 (Ω ). 1 Moreover, since σ and ˜ are positive in Ω, we may assume that σ σ 1 and ˜ σ 1 are positive in Ω . 1 Second, denote by σ , σ 2 ˜ 2 ∈ C 2 (Ω ) extensions of 1 σ and ˜, respectively. σ We can take Ω 1 close enough to the boundary of Ω so that σ 2 and ˜ σ 2 are positive in Ω . 1

glyph[negationslash]

In a similar way we extend u and ˜ as follows. u Since ( u -˜) u | Γ = 0, we can find H 1 (Ω ) extensions 1 u 1 and ˜ u 1 of u and ˜, respectively, such u that ( u 1 -˜ ) u 1 | Ω 1 \ Ω = 0. Also denote by u , u 2 ˜ 2 ∈ C 2 (Ω ) extensions of 1 u and ˜ respectively. u We can take Ω 1 close enough to the boundary of Ω so that ∇ ( u 2 + ˜ ) = 0 on Ω . u 2 1 Finally we denote by g 1 a C 1 (Ω ) 1 extension to Ω 1 of the metric g obtained in (22).

We then define U glyph[epsilon1] in Ω 1 instead of I ′ by

$$U _ { \epsilon } \coloneqq \{ ( y ^ { \prime }, y ^ { n } ) \in \Omega _ { 1 } \, \colon y ^ { \prime } \in \Sigma _ { 0 }, \text{and} \, y ^ { n } \in ( - \epsilon, \epsilon ) \}$$

Now, even if the intersection of Σ 0 with ∂ Ω 1 is tangential the function that we are integrating is compactly supported, so we can justify integration by parts in (47). glyph[square]

Using the estimates in Lemmas 4.1 and 4.2 we proceed to proving (10). Recall from Schauder theory that u, u ˜ ∈ C 2 ,α (Ω). By using in order the estimate (30), an interpolation estimate in [38, Sec. 4.3.1 ], (45), and Proposition 3.1 we obtain

(49)

$$\| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( V ^ { \prime } ) } & \leq C \| u - \tilde { u } \| _ { H ^ { 2 } ( V ^ { \prime } ) } \\ & \leq C \| u - \tilde { u } \| _ { L ^ { 2 } ( V ^ { \prime } ) } ^ { \frac { \alpha } { 2 + \alpha } } \cdot \| u - \tilde { u } \| _ { H ^ { 2 } + \alpha } ^ { \frac { 2 } { 2 + \alpha } } ( V ^ { \prime } ) \leq C \| u - \tilde { u } \| _ { L ^ { 2 } ( \mathcal { T } ^ { \prime } ) } ^ { \frac { \alpha } { 2 + \alpha } } \\ & \leq C \| L ( u - \tilde { u } ) \| _ { L ^ { 2 } ( \mathcal { T } ^ { \prime } ) } ^ { \frac { \alpha } { 2 + \alpha } } \\ & = C \left | \left | 2 \nabla \cdot \Pi _ { \nabla ( u + \tilde { u } ) } ( \mathbf J ( \sigma ) - \mathbf J ( \tilde { \sigma } ) ) \right | _ { L ^ { 2 } ( \mathcal { T } ^ { \prime } ) } ^ { \frac { \alpha } { 2 + \alpha } } \\ \tau \infty \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \ e \cdot \end{array}$$

This finishes the proof of Theorem 2.1.

Proof of Theorem 2.2. Recall glyph[epsilon1] &gt; 0 in (13). The Schauder stability estimates for solutions of elliptic equations with C α (Ω)-coefficients and C 2 ,α traces on the boundary, yield

$$\| u - \tilde { u } \| _ { C ^ { 2, \alpha } ( \overline { \Omega } ) } \leq M \epsilon,$$

for a constant M dependent only on Ω and a lower bound of σ . In particular, for glyph[epsilon1] &lt; 1 /M we have that

$$\| \tilde { u } \| _ { C ^ { 2 } ( \overline { \Omega } ) } & \leq \| u \| _ { C ^ { 2 } ( \overline { \Omega } ) } + 1.$$

Moreover, since m := min Ω |∇ u &gt; | 0, for glyph[epsilon1] &lt; m/M , we also have

$$( 5 2 ) \quad \ | \nabla ( u + \tilde { u } ) | \geq 2 | \nabla u | - | \nabla ( u - \tilde { u } ) | \geq 2 m - M \epsilon \geq m > 0.$$

Denote by δ J = J ( σ ) -J (˜). σ From Theorem 2.1 we have that there exist a C &gt; 0 independent of glyph[epsilon1] , such that

$$\| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) } \leq C \left | \left | \frac { \delta \mathbf J \cdot \nabla ( u + \tilde { u } ) } { | \nabla ( u + \tilde { u } ) | } \right | \right | _ { H ^ { 1 } ( \Omega ) }.$$

Write

$$\frac { \delta \mathbf J \cdot \nabla ( u + \tilde { u } ) } { | \nabla ( u + \tilde { u } ) | } = 2 \frac { | \nabla u | } { | \nabla ( u + \tilde { u } ) | } \left ( \frac { \delta \mathbf J \cdot \nabla u } { | \nabla u | } \right ) + \frac { \delta \mathbf J \cdot \nabla ( \tilde { u } - u ) } { | \nabla ( u + \tilde { u } ) | }.$$

Note that all the terms in (54) are C 1 (Ω)-regular.

By using the identity (19), the second term in the right hand side of (54) rewrites as

$$( 5 5 ) \quad \frac { \delta \mathbf J \cdot \nabla ( \tilde { u } - u ) } { | \nabla ( u + \tilde { u } ) | } = \frac { 1 } { 2 } ( \sigma - \tilde { \sigma } ) | \nabla ( u - \tilde { u } ) | + \frac { 1 } { 2 } ( \sigma + \tilde { \sigma } ) \frac { | \nabla ( u - \tilde { u } ) | ^ { 2 } } { | \nabla ( u + \tilde { u } ) | }.$$

Using (21), (51), and the Cauchy-Schwartz inequality we obtain

$$\dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \\ \langle ( \sigma + \tilde { \sigma } ) \nabla ( \tilde { u } - u ), \nabla ( \tilde { u } - u ) \rangle _ { L ^ { 2 } ( \Omega ) } & = - \langle \nabla \cdot ( \sigma + \tilde { \sigma } ) \nabla ( \tilde { u } - u ), \tilde { u } - u \rangle _ { L ^ { 2 } ( \Omega ) } \\ & = \langle \nabla \cdot ( \tilde { \sigma } - \sigma ) \nabla ( u + \tilde { u } ), \tilde { u } - u \rangle _ { L ^ { 2 } ( \Omega ) } \\ & \leq \| ( \tilde { \sigma } - \sigma ) \nabla ( u + \tilde { u } ) \| _ { L ^ { 2 } ( \Omega ) } \| \nabla ( \tilde { u } - u ) \| _ { L ^ { 2 } ( \Omega ) } \\ ( 5 6 ) & \leq C \| \tilde { \sigma } - \sigma \| _ { L ^ { 2 } ( \Omega ) } \| \nabla ( \tilde { u } - u ) \| _ { L ^ { 2 } ( \Omega ) }, \\ \epsilon \quad \sim \cdots \quad \cdots \quad \cdots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \circ \dots \end{array}$$

for C dependent on the C 1 -norm of u but independent of glyph[epsilon1] . From (56), for

$$\delta \coloneqq \min _ { \overline { \Omega } } ( \sigma + \tilde { \sigma } ) > 0,$$

we obtain

$$\left ( 5 7 \right ) \right ) \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \right ) \right ) \right ) \right ) \right ) } \right ) \right ) } \right ) \right ) } \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \ \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \bottom \right ) \right ) \right ) \right ) \right ) } \right ) \right ) } \right ) } \right ) \right ) } \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \right ) \right ) \left ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) } \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right ) \right )$$

We claim that

$$\left \| \frac { \delta \mathbf J \cdot \nabla ( u - \tilde { u } ) } { | \nabla ( u + \tilde { u } ) | } \right \| _ { H ^ { 1 } ( \Omega ) } \leq C \epsilon \| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) }.$$

for some C independent of glyph[epsilon1] . Indeed: one differentiation in (55) in any variable, which for simplicity we denote by subscript x , yields four terms. We treat each term individually, the appearing constants are

independent of glyph[epsilon1] . To bound the first term we use (13), the Cauchy inequality and (57) to obtain

$$\int _ { \Omega } | ( \sigma - \tilde { \sigma } ) _ { x } | \, | \nabla ( u - \tilde { u } ) | d x \leq C _ { 1 } \epsilon \| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) }.$$

To bound the second term we use the Cauchy inequality and uniform bounds (dependent on u only via (51)) on the second derivatives of u and ˜ and obtain u

$$\int _ { \Omega } | \sigma - \tilde { \sigma } | \left ( | \nabla ( u - \tilde { u } ) | \right ) _ { x } d x \leq C _ { 2 } \epsilon \| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) }.$$

To bound the third term we use uniform bounds on the second derivatives of u (and implicitly on ˜ via (51)) and on the first derivatives of u σ (and implicitly of ˜ via (13)), (52) and (57) and obtain σ

$$\int _ { \Omega } \left ( \frac { \sigma + \tilde { \sigma } } { | \nabla ( u + \tilde { u } ) | } \right ) _ { x } | \nabla ( u - \tilde { u } ) | ^ { 2 } d x \leq C _ { 3 } \epsilon \| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) }.$$

To bound the fourth term we use the Cauchy inequality and uniform bounds on σ, σ ˜, and C 2 -bounds on u (and implicitly via on ˜ (51)), u (52), and (57) and obtain

$$\int _ { \Omega } \frac { \sigma + \tilde { \sigma } } { | \nabla ( u + \tilde { u } ) | } ( | \nabla ( u - \tilde { u } ) | ) _ { x } | \nabla ( u - \tilde { u } ) | d x \leq C _ { 4 } \epsilon \| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) }.$$

This finishes the proof of (58).

From Equations (53), (54) and (58) we get

$$\| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) } \leq C \left \| \frac { \delta \mathbf J \cdot \nabla u } { | \nabla u | } \right \| _ { H ^ { 1 } ( \Omega ) } + C \epsilon \| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) },$$

for a constant C independent of glyph[epsilon1] .

Finally, using (59) for glyph[epsilon1] &lt; 1 /C we obtain

$$\| \sigma - \tilde { \sigma } \| _ { L ^ { 2 } ( \Omega ) } \leq \frac { C } { 1 - C \epsilon } \left \| \frac { \delta \mathbf J \cdot \nabla u } { | \nabla u | } \right \| _ { H ^ { 1 } ( \Omega ) }$$

which implies (15) and proves Theorem 2.2.

glyph[square]

## Acknowledgments

We would express our gratitute to Plamen Stefanov for his valuable suggestions during the writing of this paper.

## References

- 1. Giovanni Alessandrini, Critical points of solutions of elliptic equations in two variables , Annali della Scuola Normale Superiore di Pisa-Classe di Scienze 14 (1987), no. 2, 229-256. 3
- 2. , Stable determination of conductivity by boundary measurements , in Appl. Anal. [3], 153-172. MR 922775 (89f:35195) 1
- 3. , Stable determination of conductivity by boundary measurements , Appl. Anal. 27 (1988), no. 1-3, 153-172. MR 922775 (89f:35195) 19
- 4. , Global stability for a coupled physics inverse problem , arXiv preprint arXiv:1404.1275 (2014). 3
- 5. Kari Astala and Lassi P¨ aiv¨rinta, a Calder´n's o inverse conductivity problem in the plane , Ann. of Math. (2) 163 (2006), no. 1, 265-299. MR 2195135 (2007b:30019) 1
- 6. Guillaume Bal, Eric Bonnetier, Fran¸ cois Monard, and Faouzi Triki, Inverse diffusion from knowledge of power densities , Inverse Problems and Imaging 7 (2013), no. 2, 353-375. 1
- 7. Guillaume Bal, Chenxi Guo, and Fran¸ cois Monard, Inverse anisotropic conductivity from internal current densities , Inverse Problems 30 (2014), no. 2, 025001. 2
- 8. AL Bukhgeim, Recovering a potential from cauchy data in the two-dimensional case , Journal of Inverse and Ill-posed Problems jiip 16 (2008), no. 1, 19-33. 1
- 9. A. P. Calder´n, o On an inverse boundary value problem , in Seminar on Numerical Analysis and its Applications to Continuum Physics, Soc. Brasiliera de Matemarica, Rio de Janeiro, 1980, article republished in Comp. Appl. Math 25 (2006), no. 2-3. 1
- 10. Y. Capdeboscq, J.e Fehrenbach, F. De Gournay, and O. Kavian, Imaging by modification: numerical reconstruction of local conductivities from corresponding power density measurements , SIAM Journal on Imaging Sciences 2 (2009), 1003. 1
- 11. B. Gebauer and O. Scherzer, Impedance-acoustic tomography , SIAM J. Appl. Math. 69 (2008), no. 2, 565-576. MR 2465856 (2009j:35381) 1
- 12. David Gilbarg and Neil S. Trudinger, Elliptic Partial Differential Equations of Second Order , Springer, 2001. 3
- 13. K.F. Hasanov, A.W. Ma, R.S. Yoon, A.I. Nachman, and M.L. Joy, A new approach to current density impedance imaging , Engineering in Medicine and Biology Society, 2004. IEMBS '04. 26th Annual International Conference of the IEEE, vol. 1, Sept 2004, pp. 1321-1324. 1
- 14. N. Honda, J. McLaughlin, and G. Nakamura, Conditional stability for single interior measurement , Inverse Problems 30 (2013), no. 5. 3
- 15. Sungwhan Kim, Ohin Kwon, Jin Keun Seo, and Jeong-Rock Yoon, On a nonlinear partial differential equation arising in magnetic resonance electrical impedance tomography , SIAM journal on mathematical analysis 34 (2002), no. 3, 511-526. 2
- 16. Yong-Jung Kim and Min Gi Lee, Well-posedness of the conductivity reconstruction from an interior current density in terms of Schauder theory , Quart. Appl. Math. 2
- 17. Peter Kuchment and Dustin Steinhauer, Stabilizing inverse problems by internal data , Inverse Problems 28 (2012), no. 8, 084007, 20. MR 2956563 2

- 18. O. Kwon, E.J. Woo, J. R. Yoon, and J. K. Seo, Magnetic resonance electrical impedance tomography (MREIT): Simulation study of J-substitution algorithm , IEEE Transactions on Biomedical Engineering 49 (2002), no. 2, 160-167. 1
- 19. Ohin Kwon, June-Yub Lee, and Jeong-Rock Yoon, Equipotential line method for magnetic resonance electrical impedance tomography , Inverse Problems 18 (2002), no. 4, 1089. 1, 2
- 20. June-Yub Lee, A reconstruction formula and uniqueness of conductivity in mreit using two internal current distributions , Inverse Problems 20 (2004), no. 3, 847. 1
- 21. Niculae Mandache, Exponential instability in an inverse problem for the Schr¨dinger equation o , Inverse Problems 17 (2001), no. 5, 1435. 1
- 22. Fran¸ cois Monard and Guillaume Bal, Inverse anisotropic diffusion from power density measurements in two dimensions , Inverse Problems 28 (2012), no. 8, 084001. 2
- 23. , Inverse anisotropic conductivity from power densities in dimension n 3 , Communications in Partial Differential Equations 38 (2013), no. 7. 2
- 24. Carlos Montalto and Plamen Stefanov, Stability of coupled-physics inverse problems with one internal measurement , Inverse Problems 29 (2013), no. 12, 125004. 2, 3, 9
- 25. A. Moradifam, A. Nachman, and A. Tamasan, Uniqueness of minimizers of weighted least gradient problems arising in conductivity imaging , preprint (2014). 2
- 26. Adrian Nachman, Alexandru Tamasan, and Alexandre Timonov, Conductivity imaging with a single measurement of boundary and interior data , Inverse Problems 23 (2007), no. 6, 2551-2563. MR 2441019 (2009k:35325) 2, 3
- 27. , Recovering the conductivity from a single measurement of interior data , Inverse Problems 25 (2009), no. 3, 035014, 16. MR 2480184 (2010g:35340) 2, 7
- 28. , Reconstruction of planar conductivities in subdomains from incomplete data , SIAM Journal on Applied Mathematics 70 (2010), no. 8, 3342-3362. 2
- 29. , Current density impedance imaging , Tomography and inverse transport theory, Contemp. Math., vol. 559, Amer. Math. Soc., Providence, RI, 2011, pp. 135-149. MR 2885199 1, 2
- 30. Adrian I Nachman, Reconstructions from boundary measurements , Annals of Mathematics (1988), 531-576. 1
- 31. , Global uniqueness for a two-dimensional inverse boundary value problem , Annals of Mathematics (1996), 71-96. 1
- 32. M. Z. Nashed and A. Tamasan, Structural stability in a minimization problem and applications to conductivity imaging , Inverse Probl. Imaging 4 (2010). 2
- 33. G. C. Scott, M. L. G. Joy, R. L. Armstrong, and R. M. Henkelman, Measurement of nonuniform current density by magnetic resonance , IEEE Transactions on Medical Imaging 10 (1991), no. 3, 362-374, cited By (since 1996) 141. 2
- 34. Jin Keun Seo and Eung Je Woo, Magnetic resonance electrical impedance tomography (MREIT) , SIAM Rev. 53 (2011), no. 1, 40-68. MR 2785879 (2012d:35412) 1
- 35. John Sylvester and Gunther Uhlmann, A global uniqueness theorem for an inverse boundary value problem , Ann. of Math. (2) 125 (1987), no. 1, 153-169. MR 873380 (88b:35205) 1

- 36. A. Tamasan and Timonov A., Coupled physics electrical conductivity imaging , Eurasian J. Math. Comp. Appl. 2 (2014), no. 3-4, 5-29. 1
- 37. A. Tamasan, A. Timonov, and J. Veras, Stable reconstruction of regular 1harmonic maps with a given trace at the boundary , Appl. Anal 94 (2015), no. 6, 1098-1115. 2
- 38. Hans Triebel, Interpolation theory, function spaces, differential operators , 2nd ed., Johann Ambrosius Barth, Heidelberg, 1995. MR 1328645 (96f:46001) 16
- 39. Thomas Widlak and Otmar Scherzer, Hybrid tomography for conductivity imaging , Inverse Problems 28 (2012), no. 8, 084008. 1
- 40. Nanping Zhang, Electrical impedance tomography based on current density imaging , 1992. 1, 2

Department of Mathematics, University of Washington, Seattle, WA 98195-4350, USA

E-mail address : montcruz@math.washington.edu

Department of Mathematics, University of Central Florida, Orlando, Florida 32816

E-mail address : tamasan@math.ucf.edu
Carlos Montalto Alexandru Tamasan
Stability in Conductivity Imaging from Partial Measurements of One Interior Current
We prove a stability result in the hybrid inverse problem of recovering the
electrical conductivity from partial knowledge of one current density field
generated inside a body by an imposed boundary voltage. The region where
interior data stably reconstructs the conductivity is well defined by a
combination of the exact and perturbed data.