Conditional stability for an inverse source problem and an application to the estimation of air dose rate radioactive substances by drone data

We consider the density field $f(x)$ generated by a volume source $\mu(y)$ in $D$ which is a domain in $\R^3$. For two disjoint segments $\gamma, \Gamma_1$ on a straight line in $\R^3 \setminus \ooo{D}$, we establish a conditional stability estimate of H\"older type in determining $f$ on $\Gamma_1$ by data $f$ on $\gamma$. This is a theoretical background for real-use solutions for the determination of air dose rates of radioactive substance at the human height level by high-altitude data. The proof of the stability estimate is based on the harmonic extension and the stability for line unique continuation of a harmonic function.


Motivation
The Fukushima Daiichi Nuclear Disaster in March 2011 has released radioactive substances such as cesium-137 into environments. In particular, radioactive substances fall on the ground and some of them diffuses in the soil. These radioactive substances are a source and generate the air dose rate of radiation, which can be considerd as influences to inhabitants. Possible sources for such air dose rates and possible sources are in the air and on the ground, in the underground. The sources in the air could be created by floating radioactive substances, but we can assume that such sources in the air can be neglected after 9 years passed after the disaster in 2011 and the susbstances have sufficiently diffused in the air to a very low level. Thus we assume that the air dose rates can be generated by sources on the ground and the underground. As for related works on the several kinds of sources, we refer to Malins, Okumura, Machida, Takeyama and Saito [8], Saito [9] for example.
Moreover these substances on and in the ground can move for example by refloating on dry days and may run-off into rivers, so that the sources may not be stationary. However for a feasible model, we can assume that we can neglect also the time change of sources.
In some areas of Fukushima prefecture, inhabitants have already returned to daily lives, and the estimation of the air dose rate of radioactive substances in towns and villages is crucial for the sake of the health of them. The air dose rate at the human height level (e.g., 1m) is considered as very influential factor to the health through exposures by breathing the air. The observation of air dose rates can be done by drones containing measurement equipments. However direct observations of air dose rates at the human height level is not realistic because of many obstacles against the flights of drones such as houses and other many artificial structures such as walls. Moreover we do not a priori know the source, that is, the density disribution of radioactive substances on the ground and in the shallow underground. Thus, not knowing sources themselves, we can observe only data at higher altitude (e.g., 30m ∼ 50m).
Therefore our task is an inverse problem where we are required to determine the dose rate at the human height leval by means of high-altitude data. This motivates our theoretical research for the correponding inverse problem and our theoretical result shows not bad stability under adequete conditions and suggests how to control flight orbits for a reasonable accuracy in solving an inverse problem.

Mathematical model and the main result.
We describe a mathematical model, by which we can discuss our inverse problem from a general point of view.
Let D ⊂ R 3 be an open domain. Henceforth we set x = (x 1 , x 2 , x 3 ), y = (y 1 , y 2 , y 3 ) and We set (2.1) Here 4πµ(y), y ∈ D describes the density of e.g., radioactive substances. We consider f (x) as the air dose rate at x ∈ R 3 . For an arbirarily fixed constant M > 0, we set Our inverse problem is: Here we note that we assume µ ∈ M but we do not know µ.
Example: Taking into consideration the observation by drones, we choose curves Γ 1 , γ on the same curve Γ which is an orbit of a drone. For example, we choose two segments γ and Γ 1 . Let x 3 = 0 and x 3 > 0 correspond to the flat ground and the underground respectively. Let In particular, the case of D ⊂ {0 < x 3 < δ} with small δ > 0 and µ = 0 in D, means that the substances concentrate on the ground and the shallow underground. Let (a, b) ∈ R 2 be a fixed planar location and δ < h 1 < h 2 < H 1 < H 2 and This setting describes that we are requested to determine the dose rate at the heights h 1 ∼ h 2 by high-altitude data at H 1 ∼ H 2 .
We can expect that θ is smaller when dist (γ, Γ 1 ) is larger. The theorem implies the uniqueness, that is, if f = 0 on γ, then f = 0 on Γ 1 . As can be seen by the proof in Section 3, data f on γ cannot give any information of f outside of the straight line L where γ and Γ 1 are included. In other words, we can determine f only in the extended direction of γ by data on γ, provided that the extended segment is in E. In particular, we assume that 0 < δ < h 1 < h 2 < H 1 < H 2 and recall (2.2). Then γ, Γ 1 ⊂ Γ := (δ, ∞). The theorem asserts a stability estimate of Hölder type in determining x ∈ γ, and the Hölder exponent θ depends only on a geometric configuration of the two segments γ and Γ 1 . Our main theorem guarantees a rather good rate of the stability when we use drone data at high-altitude for determining the dose rate at lower-level.
Here we do not discuss the determination of the source µ(x) itself, and we do not know whether we can extract some information of µ from f | γ . On the other hand, as data we choose f in a domain D 1 ⊂ R 3 such that D ∩ D 1 = ∅, and Cheng, Prössdorf and Yamamoto [2] proves some stability for an inverse problem of determining µ| D by f | D 1 which are not data on a segment. The stability is conditional under assumption that unknown µ satisfies (2.3) and is of logarithmic rate which is much weaker than (2.4). We emphasize that data f on γ cannot determine the source µ but can determine f on Γ with Hölder stability rate which is much better than the logarithmic rate.

Proof of Theorem
The proof is based on a harmonic extension of f and the line unique continuation for a harmonic function.
We note that the harmonic extension G defined by (3.1) is applied for the determination of µ (e.g., Cheng, Prössdorf and Yamamoto [2], Cheng and Yamamoto [4]).

Numerical examples
In the numerical computation, we here propose a method via the determination of µ(x). More precisely, we choose {p 1 , p 2 , ..., p M } ⊂ γ ⊂ R 3 and {q 1 , q 2 , ..., q N } ⊂ D ⊂ R 3 as collocation points, and construct the matrix Then find the least-norm solutionμ to the equation Aμ = b, where b = (f (p 1 ), ..., f (p M )). Taking into consideration observation data f meas with errors, the Tikhonov regularization is adopted andμ is then the minimizer to the cost functional Here α > 0 is a regularizing parameter which we should choose suitably according to the noise level. A discretized form is where l i are the line integral coefficients. The regularization parameter α is chosen as α ∼ δ 2 , where δ is the noise level of observation data measured by the L 2 (γ)-norm and we refer to Cheng and Yamamoto [5] as for a choice strategy of α. The approximation solution can be constructed byμ asf In the following numerical examples, we assume that the source distribution has a compact support in B := {x = (x 1 , x 2 , x 3 ) ∈ R 3 ; |x| < 1, x 3 < 0}. The exact source distribution is given as µ( with finer grid resolution.
In order to further illustrate the performance of the method, a random noise is added on the exact observation. Figure 2 gives the results with observation error of level 1% and 3% respectively, which coincide well with the exact solution with the relative error less than 5% at x 3 = 0.1.

1.
In this paper, we dicuss the estimation of air dose rates of radioactive substances at the human height level by high-altitude data by means by drone. We establish a Hölder stability estimate for this inverse problem on the basis of a mathematical model. For determination of air dose rates in more directions, the main theorem suggests that a drone orbit should include the corresponding many directions. More precisely, let γ be a curve composed of N-segments γ 1 , ..., γ N ⊂ E. Let ℓ k , 1 ≤ k ≤ N be the extended straight lines of γ k such that ℓ k and γ k are in the same connected component of E. The line ℓ k may be a finite segment, an infinite straight line or a half straight line. The theorem means that f | γ determines f on ∪ N k=1 ℓ k uniquely. In other words, when a drone flies zigzag, drone data can determine f on all the extended directions of the components of the zigzag. Moreover the stability is rather good in determining the dose rate along the orbit direction.

2.
The key of the proof is the harmonic extension and the line unique continuation, which asserts stable continuation (3.4) of a harmonic function along the segment and is different from usual unique continuation of solutions to elliptic equations (e.g., Isakov [6], Lavrent'ev, Romanov and Shishatskii [7]). The uniqueness corresponding to (3.4) can be easily proved.
Indeed let x(τ ) with τ ∈ J: some open interval, be a parametrization of the straight line Γ. Since ∆G = 0 in some domain E ⊂ R 4 , the interior real analyticity of a harmonic function yields that G is real analytic in (x, ξ) ∈ E. This analyticity can be proved also by the representation (3.4) of G. Therefore G(τ ) := G(x(τ ), 0) is analytic in τ ∈ J. Consequently with some non-empty open interval J 0 ⊂ J, the unicity theorem of an analytic function implies that if G(τ ) = 0 for τ ∈ J 0 , then G(τ ) = 0 for τ ∈ J.
For more details on the line unique continuation, we refer to [1] and [3] 3. We give numerical examples by a method by first determining a source. Our numerical results indicated acceptably good accuracy also with noisy data. Our rather good stability can support our numerical results.