Abel transforms with low regularity with applications to X-ray tomography on spherically symmetric manifolds

We study ray transforms on spherically symmetric manifolds with a piecewise $C^{1,1}$ metric. Assuming the Herglotz condition, the X-ray transform is injective on the space of $L^2$ functions on such manifolds. We also prove injectivity results for broken ray transforms (with and without periodicity) on such manifolds with a $C^{1,1}$ metric. To make these problems tractable in low regularity, we introduce and study a class of generalized Abel transforms and study their properties. This low regularity setting is relevant for geophysical applications.


Introduction
Our aim is to study geodesics and ray transforms on a spherically symmetric manifold. More specifically, our manifold is the Euclidean annulus M =B(0, 1) \B(0, R) ⊂ R n , R ∈ (0, 1) and n 2, with the metric g(x) = c −2 (|x|)e(x), where e is the standard Euclidean metric and c : (R, 1] → (0, ∞) is a piecewise C 1,1 function satisfying the Herglotz condition (see definition 1). A piecewise C 1,1 function may have a finite number of jump discontinuities.
By a geodesic we mean a maximal unit speed geodesic on the Riemannian manifold (M, g) with endpoints at the outer boundary ∂M := ∂B(0, 1). One of the problems we study is to recover a function f ∈ L 2 (M) from its integrals over all geodesics. The same problem with f ∈ C ∞ and c ∈ C ∞ was considered by Sharafutdinov [35]. For f ∈ L 2 and c ∈ C ∞ the result follows from the local support theorem by Uhlmann and Vasy [40] since the Herglotz condition is equivalent with the foliation condition in spherical symmetry. Reducing regularity introduces technical difficulties but makes the problem more relevant for seismic imaging. 3 Author to whom any correspondence should be addressed.
In addition to the x-ray transform, we consider the broken ray transform and the periodic broken ray transform. In order to prove injectivity results for these transforms, we need various integral transforms: Abel-type transforms, the Fourier series, the Funk transform, and the planar average ray transform.
See section 1.2 for an overview of our methods and results.

The Herglotz condition
We now define precisely what we mean by the Herglotz condition (named after Herglotz [13]) and geodesics in the low regularity setting. We will also introduce so-called countable and finite conjugacy conditions. Definition 1. By a piecewise C 1,1 function satisfying the Herglotz condition we mean a function c : (R, 1] → (0, ∞) satisfying the following: • The interval (R, 1] is a finite disjoint union of intervals (a, b] so that c is C 1,1 regular on each such interval. • On each subinterval (a, b] the function satisfies the Herglotz condition (The function c is upper semicontinuous from the right.) The jump condition (2) simply requires the inequality (1) when the left-hand side is interpreted as a signed measure or a distribution rather than a function. The Herglotz condition for regular c means that all geodesics reach the boundary, and with less regularity we have to exclude trapping due to total internal reflection. In other words, the Herglotz condition is equivalent with the manifold being non-trapping.
Consider a subinterval (a, b] ⊂ (R, 1] with a = R. By (1) the function r/c(r) increases on this interval, so it has a limit at a in [0, ∞). Combining this with the jump condition (2) at r = a shows that lim r→a+ c(r) exists in (0, ∞). Using the conserved quantities introduced in section 5.1 below, we conclude that any geodesic approaching the surface r = a from either direction has a well-defined limit-for both position and direction.
The above argument actually shows that c is bounded away from zero and infinity on (R + ε, 1] for any ε > 0. The definition does permit lim r→R c(r) = ∞.
Total internal reflection from the outside is possible, but we exclude all geodesics that reflect where the wave speed c jumps. In practice such geodesics would give valuable data for the inverse problem, but that data are unnecessary for uniqueness and we therefore omit it. We also exclude geodesics that are tangent to a surface where c jumps. Geodesics that meet such a surface are therefore assumed to traverse it according to Snell's law. Our geodesics do not branch; they are completely transmitted. This is not merely a matter of technical convenience. If one wants to take branching into account, the whole concept of x-ray transform needs to be redefined. This generalization is particularly inobvious for the periodic broken ray transform.
We do not know whether the Herglotz condition and piecewise C 1,1 regularity are necessary conditions for our injectivity results. Without these assumptions the geometric framework starts breaking apart and the problem may need a reformulation. There are geodesic x-ray tomography results on various manifolds without symmetries (see section 1.3), but our method of proof relies heavily on symmetry. In general, one may extend injectivity results to nearby geometries by a perturbation argument if one has stability in addition to injectivity. However, we have little or not restrictions for conjugate points and cannot therefore expect stability.
In addition to the Herglotz condition, our results for the broken ray transforms (periodic and not) require an additional assumption. This assumption may in fact follow from the Herglotz condition, but we state it separately since we have no proof of the implication. Definition 2. We say that a C 1,1 wave speed c satisfies the countable conjugacy condition if there are only countably many maximal geodesics whose endpoints are conjugate, up to rotational symmetry.
For some negative results (see theorem 36) we need a stronger condition. Definition 3. We say that a C 1,1 wave speed c satisfies the finite conjugacy condition if there are only finitely many maximal geodesics whose endpoints are conjugate, up to rotational symmetry.
The countable conjugacy condition holds for most wave speeds, and assuming it is unlikely to be an issue in applied problems.

Overview of results and methods
Our main results concern ray transforms on spherically symmetric manifolds. See the theorems mentioned here for more details and related results.
• Assuming a piecewise C 1,1 wave speed satisfying the Herglotz condition, the attenuated geodesic x-ray transform is injective when the attenuation is radially symmetric and Lipschitz continuous. (See theorem 29.) • Assuming a C 1,1 wave speed satisfying the Herglotz condition and the countable conjugacy condition, the broken ray transform is injective on functions that are quasianalytic in the angular variable(s). The broken rays have endpoints on the boundary (they are not periodic), and the endpoints can be limited to any open set. In two dimensions, for example, a function of the form satisfies the regularity assumptions if each a k is Hölder continuous. (See theorem 33.) • Assuming a C 1,1 wave speed satisfying the Herglotz condition and the countable conjugacy condition, the periodic broken ray determines the even part of the function in three or more dimensions. In two dimensions averages over circles centered at the origin are determined. Very little other information can be recovered in any dimension. (See theorem 36.) The tools used in these proofs are similar. They all rely on analysis disc by disc and Abel transforms. The piecewise C 1,1 metric is dealt with using layer stripping. Fourier series in the angular variable(s) plays a crucial role.
The periodic broken ray transform is closely related to a new integral transform we introduce, the planar average ray transform. The kernel of this new transform is precisely the space of odd functions (see theorem 41), and the mentioned result for the periodic broken ray transform follows. The proof of this theorem again relies on Abel transforms, but also on the Funk transform.
In order to prove the mentioned results for ray transforms, we give several results for various other integral transforms in a low regularity setting. For injectivity of Funk-type transforms on distributions, see theorem 23. In addition, we prove that a large class of Abel-type integral transforms is injective (see theorem 12). The proof is based on finding a local approximate inverse and using Neumann series together with layer stripping. Section 2 contains several results regarding continuity, injectivity and differentiability properties of these integral transforms.

Applications and related results
The most important spherically symmetric manifold we have in mind is the Earth. In the preliminary reference earth model (PREM) both pressure and shear wave speeds are piecewise nice and satisfy the Herglotz condition in the mantle [7]. The Herglotz condition is violated at the core-mantle boundary (CMB) but is satisfied at all jump discontinuities above it. Furthermore, the Herglotz condition for the pressure wave speed only fails at the CMB and for shear wave speed only in the liquid outer core. Geometrically, failure of the Herglotz condition at a jump discontinuity means that rays can be trapped by total internal reflection.
We consider a piecewise C 1,1 satisfying the Herglotz condition in an annulus. In PREM, both wave speeds satisfy these assumptions, both regularity and geometry, in the whole mantle.
The attenuated x-ray transform has a direct application in an imaging method known as SPECT. It also appears in the anisotropic Calderón's inverse boundary value problem [6,23]. For a review of attenuated x-ray tomography we refer to [9]. For other results for attenuated ray transforms on manifolds, see e.g. [1,24,25,32,34].
The standard x-ray transform without attenuation has applications in imaging methods such as CT and PET. The Euclidean version of this problem is the starting point of the study of inverse problems in integral geometry (see e.g. the classical works [4,12,27,30]), but our main interests are in non-Euclidean geometry. It appears in the linearizations of some inverse problems, such as the travel time tomography problem where one attempts to recover a manifold from the pairwise distances of boundary points. Solving the linearized problem leads to results known as boundary rigidity; see [37] for a review. See the book [33] for the geodesic x-ray transform and its applications.
A less studied variant of x-ray tomography is broken ray tomography, where instead of geodesics one considers broken rays which reflect at the boundary of the manifold. One may either consider broken rays with endpoints in a given subset of the boundary, or periodic broken rays which do not terminate at all. For recent results in broken ray tomography, we refer to [8,[14][15][16][17][18][19][20][21]. 'Broken ray tomography' may also refer to problems where the reflections take place in the interior rather than the boundary of the domain; see [18] for a discussion of terminology.
The broken ray transform appears in similar ways as the x-ray transform when reflections from the exterior boundary or inaccessible interior boundaries are present. For inverse boundary value problems for PDEs with partial data, we refer to [8,22]. Differentiating the length of a broken ray with respect to the metric leads to the broken ray transform [21, theorem 17] in analogy to the more familiar situation without reflections.
Seismic travel time tomography for USArray [2] (there are numerous subsequent improvements as more data have become available) is based on broken ray tomography. We will study broken ray tomography with measurements in a small boundary set in section 7. The prevalent inversion method is to make a least squares fit, but a better understanding of broken ray tomography is likely to lead to advances in seismic tomographic inversion. For a comparison of linearized tomographic travel time inversion and the non-linear Herglotz-Wiechert inversion method in radial symmetry, and a discussion of Abel transforms in this context, see [28]. Earlier work on the Abel transform will be discussed in section 2.
The periodic x-ray transform on closed manifolds is related to spectral rigidity problems. This relation is encoded by the length spectrum, the set of lengths of all periodic geodesics. For a prominent example, see [29]. In a similar fashion, the periodic broken ray transform is related to the rigidity of the spectrum and the length spectrum of a manifold with boundary. The first results in this direction can be found in [5] where the problem is studied on spherically symmetric manifolds. 1.4. The structure of the paper Sections 2-4 are dedicated to tools we need for ray transforms. This includes Abel transforms (section 2), Funk transforms (section 4), and other tools (section 3). A reader whose sole interest lies in ray transforms may skip these sections. Some of our auxiliary results are well known, but we present them with proofs for completenees and ease of reading.
We will discuss geodesics and broken rays in radial symmetry in section 5. With the key tools ready for our disposal, we will focus on ray transforms in sections 6-8.

Abel transforms
We start with studying Abel transforms. We will need to understand these transforms to prove injectivity results for ray transforms in subsequent sections, but we look at more general Abel transforms than is needed for ray tomography. For a review of the Abel transform, see [3], and for more about general Abel-like transforms, see [4,10,16,26,[41][42][43][44]. We are not sure of originality of the results in this section but we were unable to find the exact results we need in the literature. The case of K ≡ 1 is well known.
The notations of this section are somewhat different from those elsewhere in this paper.

Definitions
Let ∆ = {(x, y); 0 x y 1} and let K : ∆ → R be any bounded measurable function. Fix a number α ∈ [0, 1). For a function f : [0, 1] → R we define its integral transform whenever this integral is defined, and we let I α K f (1) = 0. These transforms generalize the classical Abel transform I The kernel K is often assumed to be in the space Lip(∆) of Lipschitz functions from Δ to R, but some results hold for K ∈ L ∞ (∆) as well. We denote the Lipschitz constant with respect to the first variable by lip 1 (K) for K ∈ Lip(∆).

Continuity
We establish continuity in two senses: the transform is continuous between suitable function spaces, and the image of a continuous function-and some L p functions-is continuous. Theorem 4. The transform I α K : L p → L q is well defined and continuous when α + 1/p < 1 + 1/q. In particular, this holds when p > 1/(1 − α), q < 1/α, or p = q. The norm of this mapping satisfies I α K L p →L q = O(sup ∆ |K|). Proof. We denote by inequalities involving constants independent of f. If K vanishes identically, the result is obvious. Otherwise we divide K by sup ∆ |K| so that we may assume |K| 1; the last part of the claim follows from this scaling.
We have This is the desired continuity estimate. □ M V de Hoop and J Ilmavirta Inverse Problems 33 (2017) 124003 Theorem 5. If K is continuous and p ∈ (1, ∞] is such that αp/( p − 1) < 1, then I α K : L p → C is continuous. Proof. It suffices to show that for every f ∈ L p the function I α K f is continuous. The space of continuous functions is a closed subspace of L ∞ and theorem 4 gives continuity in the L ∞ norm.
We denote g = I α K f . Let us first prove that g is continuous at 1. It is a natural interpretation of the definition of the integral transform that g(1) = 0. We can then use Hölder's inequality to find which establishes continuity at x = 1.
By definition we have Making the change of variable y Therefore We will show that each of the three terms tends to zero as x → x.
Let us define f x : ; notice that f x = f . We extend f by zero to a larger interval J = [− 1 1−a , 1] to make each f x well defined and we also extend f x to J by zero where the original formula would take the argument of f outside J. Let us show Fix any ε > 0. By density of continuous functions, there is a continuous function h ∈ L p (J) so that f − h L p (J) < (1 − a)ε/3. Let us define h x in terms of h like f x was defined in terms of f. Since these scaled and shifted functions differ from the original ones only by an affine change of variable, we get The function h is uniformly continuous, so there is This proves the desired auxiliary result: This implies that the functions f x have a uniform bound on their L p norm for x sufficiently close to x. Using Hölder's inequality like above, we find Since this bound is independent of x , the first term in (9) tends to zero as x → x. Essentially the same estimate together with uniform continuity of K shows that the second term tends to zero as well. The same estimate works for the third term too, since f x → f in L p ([0, 1]). □ We will also need a Lipschitz version of the previous theorem.
Theorem 6. If f and K are continuous and α ∈ [0, 1), then I α K f is continuous. If α > 0, and f and K are additionally Lipschitz, then I α K f is (1 − α)-Hölder everywhere and locally Lipschitz on [0, 1). In particular, the local Lipschitz constant of The estimate for the local Lipschitz constant implies, in particular, that the Lipschitz constant of the restriction Proof of theorem 6. We denote g = I α K f and prove the first claim first. By definition we have (15) for any x ∈ [0, 1). Take any x ∈ [0, 1). Making the change of variable y = 1−x 1−x y + x−x 1−x in the above integral gives Since f and K are uniformly continuous, it is easy to see that g(x ) → g(x) as x → x with a Hölder modulus of continuity. We have thus shown that g is continuous on [0, 1), so it remains to show that g(x) → 0 as x → 1. But this is elementary, as f and K are bounded.
The second claim follows easily from expression (16) for g(x ). □ The second claim of theorem 6 cannot be improved significantly:
To study injectivity of the Abel transform I α K , we define the integral which will turn out to be of great use. A simple application of Fubini's theorem shows that We will show injectivity of I α K under some conditions by showing that I 0 whenever x < y. This gives the following result; see [4] for an early use of such identities in ray tomography.

Proposition 7.
The mapping I α 1 : L p → L q is injective for any α ∈ (0, 1) whenever p and q are such that it is well defined. We have the inversion formula which holds almost everywhere. The formula holds pointwise for continuous f.

M V de Hoop and J Ilmavirta Inverse Problems 33 (2017) 124003
Proof. For f ∈ L 1 we have by (18) and (19) from which the result follows. The required version of the fundamental theorem of calculus for L 1 functions follows from Lebesgue's differentiation theorem. □ The general transform I α K can be inverted in a similar manner if one only finds such a kernel L : gives this result; the corresponding integral transform I 1/2 K is closely related to the Radon transform in the plane [4,16]. Here T k is the kth Chebyshev polynomial. The apparent singularity at zero is not an issue for injectivity; theorem 12 works so that any neighborhood of zero can be easily omitted.
For general K and α it is difficult to find a suitable L, so we approach the problem in a different way. The downside of the method below is that it does not yield an explicit inversion formula like (20). The result of proposition 7 can be slightly generalized as the next proposition demonstrates.

Proof. A simple calculation shows that
Let A be the space of absolutely continuous real functions on [0, 1] which vanish at 1. It is a classical result (see e.g. [38, chapter 3, section 3.2]) that I 0 1 (L 1 ) = A, so that we may equip A with the norm that makes I 0 1 : L 1 → A an isometry. For η ∈ [0, 1) we will use the subscript η in a function space to indicate restriction to [η, 1] instead of [0, 1]. Thus we write L p η = L p ([η, 1]) and A η for absolutely continuous functions on [η, 1] vanishing at 1. We write ∆ η = {(x, y); η x y 1}.

Proof. It follows immediately from
The and This estimate remains true if η is increased. We also have the estimate ϕ α η c α for any η.
by absolute continuity of the Lebesgue integral. Once we establish absolute continuity, this shows that g(1) = 0. For absolute continuity we have the estimate Let ε > 0; we wish to find δ > 0 so that whenever the intervals (a i , b i ) are disjoint and their union has measure at most δ. But this follows now immediately from the estimate (27), integrability of f and absolute continuity of the Lebesgue integral. We have now proven that Since g ∈ A, its derivative g exists almost everywhere and is in for almost all x.
Thus for any η ∈ [0, 1) This shows that ϕ α 2M as claimed. Since all variables in Δ take values in ∆ η only, we may replace supremum over Δ with supremum over ∆ η .
For a lower bound on the operator norm we use the constant kernel K ≡ 1. We have Proof. Since K cannot change sign at the diagonal, we may assume m : Since L(L 1 , L 1 ) is a Banach algebra, we may use the Neumann series to invert Dϕ α (K) = E − F: This concludes the proof. □ ; 0 x 1}, and Lipschitz continuous in some neighborhood of the diagonal.
One can write an inversion formula for theorem 12 in terms of Neumann series and iteration as may be read in the proof.

Remark 14.
The restriction result of theorem 12 (if I α K f vanishes above r, so does f ) is related to support theorems for ray transforms. Indeed, this result combined with the analysis of [4] proves Helgason's famous support theorem. Helgason [12, theorem 4.2] gave a different proof. More generally, this result gives a support theorem for ray transforms on spherically symmetric manifolds; see section 6.
Proof of theorem 12. We will only prove the theorem for α > 0. If α = 0, we may con- Apart from removing this one step the proof is unchanged.
The kernel K is strictly positive and Lipschitz continuous in {(x, y) ∈ ∆; y − x ζ} for some ζ ∈ (0, 1]. We will only use the values of K in this strip near the diagonal, so we may replace K with a Lipschitz extension elsewhere without altering the result and assume that K is Lipschitz in Δ. Let η ∈ [0, 1) be the constant of lemma 11 related to K and α. We denote δ = 1 − η; by possibly slightly decreasing δ, we may assume that (1 − r)/δ =: n is an integer. Lemma 11 remains true if δ is decreased and η thus increased. The and by lemma 11 this implies that f η = 0. Therefore f vanishes almost everywhere in [1 − δ, 1]. If n = 1 (η = r), we are done, so we assume that n 2. We for all x r + δ. If additionally x 1 − 2δ, lemma 11 may again be used to see that g vanishes on [ We can carry on inductively n times (recall that 1 − nδ = r), and finally conclude that f vanishes on [r, 1]. □

Differentiability
So far in this section we have studied Abel transforms with very low regularity. In this subsection we will study how Abel transforms preserve differentiability.
In particular, we study functions of the form where φ is a regular function. The question is how much regularity f inherits from φ, and an answer is given by the following proposition. The proposition will not be used in this paper directly, but it is employed in the follow-up work [5] and is used in the proof of proposition 28.
Proposition 15. Suppose ϕ : ∆ → R is continuous and k times continuously differentiable in the interior of Δ. Provided that α ∈ (0, 1), the function f defined by (35) is continuous in [0,1] and k times continuously differentiable in (0, 1). Furthermore, Proof. Continuity follows from theorem 6, so we only prove differentiability. Fix x ∈ (0, 1) and let x ∈ (0, 1). We will eventually let x → x and show that the difference quotient has the correct limit. Let k = 1 first.
Adding tildes to (35), we have We change the integration variable from ỹ to y so that y 2 − x 2 =ỹ 2 −x 2 . We write ỹ instead of the radical y 2 +x 2 − x 2 for the sake of brevity and legibility. We find Simple calculations give the following approximations: Using these, we obtain This gives (36). The general case k > 1 follows easily by induction on k, using (36). The integrands in (35) and (36) (36) is also true in the trivial case α = 0.

Fourier series and discs in balls
This section is devoted to auxiliary results that facilitate connecting the Abel transforms to ray transforms.

Measurability and integrability on submanifolds
The lemma below is given in more generality than needed, but it causes no added effort for the proof. Let G n 2 denote the Grassmannian of two-dimensional subspaces of R n . We give the next lemma only for the Euclidean metric, but this causes no loss of generality in the subsequent proofs. For another version of this statement, see [36].
where c > 0 is a constant and H 2 is the Hausdorff measure of dimension two. In particular, f | A∩P ∈ L p (A ∩ P) for almost every two dimensional plane P ∈ G n 2 .

Proof. It suffices to show
for all f ∈ L 1 (A). The identity (41) then follows after replacing f with |f | p .
Let us consider the set and the projection π : E → A to the first coordinate. This is a smooth bundle over A and A can be viewed as a Riemannian manifold with boundary. The manifold E inherits a natural Riemannian metric from the tangent bundle of A.
Since E has locally a product structure and each fiber π −1 (x) has the same measure, we have π * f ∈ L 1 (E) and It will be easier to convert the integral over E to an integral involving the Grassmannian. Consider the mapping F : E → G n 2 defined so that F(x, v) is the unique plane containing both x and v. This is a smooth surjection with an everywhere surjective differential, so the smooth coarea formula gives where σ P is the measure on the level set F −1 (P) and NJF is the normal Jacobian of F.
We could calculate NJF(z) explicitly, but we can deduce its form without such calcul ations. We have π(F −1 (P)) = P ∩ A for every P ∈ G n 2 , and due to rotation symmetry we must have for every plane P ∈ G n 2 , where φ is some weight function. Combining equations (44)-(46), we find If we let f be a spherically symmetric function supported near radius r and observe how both sides of our identity scale as r varies, we find that ϕ(r) = r n−2 up to a multiplicative constant. This proves identity (42). □ Remark 18. Lemma 17 is otherwise valid if we use the ball instead of the annulus (corresponding to R = 0), but the last statement is false. That is, equation (41) is valid but the restrictions of an L p function need not be in L p . The reason is that functions can concentrate near the origin but this is counteracted by the weight |x| n−2 in (41). For example, the function f (x) = |x| −n+1/2 is in L 1 (B(0, 1)) but its restriction to any proper subspace is not integrable.

Remark 19.
It is not hard to modify the proof given above to show that a function in L p (S n−1 ) restricts to an L p function on almost every great circle in the sense of the natural measure on the Grassmannian G n 2 .

Angular Fourier series
In this subsection we show that an L 2 function on a spherically symmetric surface can be written as a convergent Fourier series. It suffices to prove the next lemma in Euclidean geometry.

Lemma 20. Let
where a k : [R, 1] → C are L 2 functions. This series converges to f in L 2 (A). The L 2 norm of f can be written as Proof. We define Since f ∈ L 2 (A), each a k (r) is well defined for almost all r ∈ [R, 1].
To study convergence, we define the truncated series and an error function for m 0.
Since f ∈ L 2 (A), we have f (r, ·) ∈ L 2 (S 1 ) for almost all r ∈ [R, 1]. Thus for almost all r the Fourier series for f (r, ·) converges in L 2 (S 1 ), and 0 g m (r) → 0 monotonically almost everywhere as m → ∞. Monotonicity implies, in particular, that g m g 0 pointwise. Since g 0 (r) = f (r, ·) 2 L 2 (S 1 ) , we know that g 0 ∈ L 1 ([R, 1]). Therefore we may use the dominated convergence theorem to conclude that The estimate g m g 0 ∈ L 1 shows that f m ∈ L 2 (A) for all m, as can be seen in the above calculation.

The Funk transform
The Funk transform takes a function on the sphere S n to a its integrals over geodesics. It is a classical result (see lemma 22 below) that the Funk transform of a smooth function vanishes if and only if the function is odd with respect to the antipodal reflection. We extend this result to distributions-the result is probably not new but we give an explicit proof as we have been unable to locate one in the literature.

Funk transforms on smooth functions
Let G n k be the Grassmannian of k-dimensional subspaces of R n . We will work with G n k = G n+1 k+1 instead so that S n ∩ P has dimension k for all P ∈G n k . For 0 < k < n, we define the Funk transform F n k : All integrals we encounter in this section come with a natural measure on a sphere or a Grassmannian, and we only take average integrals. We could include the cases k = 0 and k = n, but these would only correspond to taking the even part and the average of the function.
We will use the integral operators I n k,l : This transform is known as the Radon transform on the Grassmannian [31]. These operators satisfy I n k,k = id, I n k,l • I n m,k = I n m,l , and whenever the indices are such that the operators are defined.

Duality
The operators introduced above have natural dual operators with respect to L 2 inner products on the manifolds involved. These are (F n k ) * : and These will allow defining the transforms I n k,l and F n K on distributions by duality. On a compact manifold M, let us denote the duality pairing between C −∞ (M) and C ∞ (M) by ·, · M . We now define F n k : for all f ∈ C −∞ (S n ) and g ∈ C ∞ (G n k ). When restricted to smooth functions, this agrees with the original definition.
Similarly, we define I n k,l : for all f ∈ C −∞ (G n k ) and g ∈ C ∞ (G n l ).

Convolutions
The group SO(n + 1) acts naturally on S n and G n k . The integral transforms F n k and I n k,l and their duals are equivariant under this action. We will use this action to define convolutions on our function spaces.
For η ∈ C ∞ (SO(n + 1)) and f ∈ C −∞ (S n ), we define the convolution η * f ∈ C ∞ (S n ) by where dg is the normalized Haar measure on the group SO(n + 1). We define the convolution similarly if f ∈ C ∞ (G n k ). It is a straightforward calculation to observe that for any η ∈ C ∞ (SO(n + 1)) and f ∈ C ∞ (S n ). Alternatively one can define the convolution first for only smooth functions and then extend to the case f ∈ C −∞ (S n ) by duality. We denote by η the reflected version of a function η ∈ C ∞ (SO(n + 1)), given by η (g) = η(g −1 ). One can define the convolution of a distribution f ∈ C −∞ (S n ) and a smooth function η ∈ C ∞ (SO(n + 1)) so that for every h ∈ C ∞ (S n ). This agrees with the definition given above.

Characterization of the kernel
Now we have the tools needed to characterize the kernel of the Funk transform on distributions. The corresponding result for smooth functions is a classical one:

Lemma 22 ([12, theorem 1.7, section 1B, chapter III]).
For any two integers 0 < k < n the kernel of the Funk transform F n k on smooth functions is precisely the space of odd functions.

Theorem 23. For any two integers 0 < k < n the kernel of the Funk transform F n k on distributions is precisely the space of odd distributions.
Proof. It is easy to observe that any odd distribution is indeed in the kernel. It therefore suffices to show that an even distribution with vanishing Funk transform has to vanish.
Let f be an even distribution on S n so that F n k f = 0. Take any smooth functions η ∈ C ∞ (SO(n + 1)) and h ∈ C ∞ (G n k ). Combining equations (63) and (64) gives and hence (That is, the intertwining property (63) holds for distributions.) Since η * f is an even smooth functions, it follows from lemma 22 that η * f = 0. This holds for any η ∈ C ∞ (SO(n + 1)). If we take a sequence (η m ) m∈N of smooth functions on SO(n + 1) converging weakly to the delta function at the identity element, we have f = lim m→∞ (η m * f ) in the sense of distributions. Since each η m * f vanishes, we must have f = 0. □ Since I n k,n−1 and F n k take even distributions to even distributions and we have the identity I n k,n−1 • F n k = F n n−1 by lemma 21, it would have sufficed to show the desired result for k = n − 1 above.

Remark 24.
Let us briefly outline an alternative proof of theorem 23 (see [39]). By equation (63) the group SO(n + 1) intertwines with the Funk transform. Therefore the kernel is a closed SO(n + 1)-invariant subspace. This implies that SO(n + 1)-finite functions (functions contained in finite dimensional representations) are dense in it. All SO(n + 1)-finite distributions are in fact smooth, and smooth elements in the kernel are known to be odd. Therefore all distributions in the kernel are limits of odd functions and therefore odd. For more details, see [11].

Geodesics and broken rays in spherical symmetry
In this section we will first review the basics of geodesics in spherical symmetry and then calculate integrals over functions along these geodesics.

Basic facts about geodesics
If we intersect our manifold M ⊂ R n (defined in the introduction) with any linear subspace of R n , we get a totally geodesic submanifold satisfying the Herglotz condition. Therefore it suffices to study geodesics in the case n = 2.
In two dimensions it is convenient to use polar coordinates (r, θ). A geodesic is uniquely determined by the location of its tip, the point with smallest r. Let us denote the coordinates of this tip by (r 0 , θ 0 ) and call r 0 the radius of the geodesic. There are two conserved quantities: the squared speed c(r(t)) −2 [r (t) 2 + r(t) 2 θ (t) 2 ] = 1 and the angular momentum c(r(t)) −2 r(t) 2 θ (t) = r 0 /c(r 0 ). These two quantities are conserved also when c has a jump discontinuity-these conservation laws at a jump point give Snell's law.
These two conserved quantities allow one to carry out most calculations quite simply. Let us calculate the length of a geodesic with a tip at radius r. Suppose the geodesic is parametrizesd by arc length as γ : [−L, L] → M. Due to symmetry, the length is 2 L 0 dt. We change the variable of integration from time t ∈ (0, L) to radius s ∈ (r, 1). The Herglotz condition ensures that this change of variable is possible, and that L < ∞. Using the two conserved quantities, we find that After simplification, we find that the length is In fact, the Herglotz condition is equivalent with all maximal geodesics having finite length. Similarly the angular distance between the endpoints of the geodesic (measured on the universal cover, may exceed 2π) is The integrand differs from that of (68) by the term θ which can be found using conservation of angular momentum. The lengths and angles can be expressed in this simple form regardless of discontinuities in c as long as the Herglotz condition is satisfied. The total length or opening angle of a geodesic can be decomposed into parts between the jumps. For each part where c ∈ C 1,1 , one can use the above formulas with appropriate changes to the limits of integration. The total integral is the sum of these parts. This makes no difference in the formulas since the conserved quantities are conserved across the jumps.

Integrals of functions over geodesics
Suppose for now that c ∈ C 1,1 without any jumps.
We consider functions of the form f (r, θ) = a(r)e ikθ on M, where a ∈ L 2 ([R, 1]) and k ∈ Z. For convenience, we will assume that a is continuous in the calculations, but the conclusions will hold true for L 2 functions with obvious modifications. It follows from lemma 20 that any function in L 2 (A) can be written as a sum of such functions, so we lose no generality studying functions of the chosen form.
Consider a geodesic γ parametrized by [−T, T] t → (r(t), θ 0 + ω(t)) with r(0) = r 0 and ω(0) = 0. Using the results of the previous subsection to convert arc length integrals to radius integrals and the identity ω(t) = Here we have With this result injectivity of ray transforms can be reduced to injectivity of the integral transform A k .
To simplify notation, we define so that By the analogue to the generalized Abel transforms studied in [16], we call T k the Chebyshev functions.
To make the integral transform more tractable, we make another change of variable: we change from r to r = ρ(r) = r/c(r), and denote other variables similarly with tildes. The Herglotz condition makes ρ a diffeomorphism. Multiplicative constants in the function c are irrelevant, so we may assume that c(1) = 1 and so ρ(1) = 1. We first notice that the singularity at r = z is simpler in the new variables. We denotẽ We write ã = a • ρ −1 (this implies ã(r) = a(r)) and define the transforms Ã k so that Ã kã (x) = A k a(x). In this notatioñ

Attenuated integrals
In section 5.2 we calculated the geodesic ray transform of functions of the form f (r, θ) = a(r)e ikθ for k ∈ Z. Now we will calculate the attenuated geodesic ray transform for a spherically symmetric attenuation λ on A; this calculation will be used to explicitly express the attenuated geodesic ray transform of a general function in L 2 (A) in section 6 below. Let λ : [R, 1] → R be a continuous function, fix some geodesic γ = γ r0,θ0 , and define where Since Λ + (0) = Λ − (0) and r(−s) = r(s), we have and thus Following the steps of equation (70), we obtain where and

Abel-type integral transforms
We define the integral transform A k taking functions on [R, 1] to functions on [R, 1] bỹ The nonsingular part of the kernel of this integral transform, K k (r, x) =T k (r; x)H(r; x) is Lipschitz continuous in r and x and K k (x, x) > 0 for all x. (The function H is obviously Lipschitz, and T k can be shown to be, too, with a proof identical to that of lemma 26.) Therefore Ã k falls in the category of integral transforms studied in section 2. Analogously with the Abel transform A k defined in equation (75) we define the attenuated Abel transform Zero attenuation gives the previously defined Abel transform: A 0 k = A k . The Abel transforms in the following result correspond to any C 1,1 wave speed satisfying the Herglotz condition. Proof. We define r; x)). It suffices to show that F 2 is Lipschitz continuous and vanishes on the diagonal, since cosh(F) is a Lipschitz function of F 2 . We will again change variables from r and x to r = ρ(r) and x = ρ(x). If we define F by demanding F (r;x) = F(r; x) and λ = λ • ρ −1 , we have (see equation (79))

Furthermore, the integral transform
Since ρ is bi-Lipschitz, we need to show that F 2 is Lipschitz and vanishes on the diagonal when λ is Lipschitz. Let lip 1F (r;x) denote the local Lipschitz constant of F with respect to r at (r;x), lip 2F (r;x) similarly, and lipλ(x) the local Lipschitz constant of λ at x.
Elementary estimates yield for some constant C.
For fixed r , the function F (r; ·) is an Abel-type integral transform of λ in the sense of section 2. Therefore we have from theorem 6 that for some constant C.
Combining estimates (91)-(93) we get for i = 1, 2, and F 2 is Lipschitz. Estimate (92) shows that F 2 vanishes on the diagonal, and the proof is complete. □ Remark 27. Proposition 8 gives a simple inversion formula for some Abel transforms. The transforms Ã 0 and A 0 do not fall in this category, but a calculation verifies an inversion formula for A 0 . We do not know of such a simple formula for other k, except in the Euclidean case (see [16, lemma 105]). The inversion formula for A 0 is and it holds almost everywhere for f ∈ L 1 ([R, 1]) and everywhere for continuous f.

Broken rays
Due to spherical symmetry, studying broken rays is simple once one understands geodesics. Just like a non-radial geodesic, a non-radial broken ray is contained in a unique two-dimensional subspace of R n and it therefore suffices to consider dimension two. The geodesic segments that constitute a broken ray only differ by rotations of the plane. Each of them has the same radius (radial coordinate of the tip), so we may call it the radius of the broken ray. It follows from spherical symmetry that if two distinct rotations of a geodesic share an endpoint, then they together form a broken ray that satisfies the usual reflection condition: the angle of reflection equals the angle of incidence.
Not all broken rays are periodic. Due to rotation symmetry periodicity only depends on the radius. Using the angle α defined in (69), it is easy to see that the broken ray corresponding to r ∈ (R, 1) is periodic if and only if α(r) ∈ πQ.
To prove injectivity results for the periodic broken ray transform we need to have enough periodic broken rays. One such result is provided by the following proposition, which assumes the countable conjugacy condition. This proposition is the sole reason for making the assumption in the results presented in this paper. (R, 1).

Proposition 28. If a C 1,1 wave speed satisfies the Herglotz condition (definition 1) and the countable conjugacy condition (definition 2), then the set of radii corresponding to periodic broken rays is countable and dense in
The proof of the proposition is somewhat involved, so we will not present it here. We only mention that by proposition 15 the opening angle α(r) defined in (69) is C 1 and that conjugate points at the boundary correspond to zeros of the derivative α . A proof can be found in [5, lemma 4.5].

The x-ray transform
Let us recall the definitions of the x-ray transform and the attenuated x-ray transform on a manifold M with boundary. The x-ray transform of a function f : M → C is a function on the space of maximal geodesics. Evaluated at a maximal geodesic γ : The linear operator I is the x-ray transform. Here the midpoint of the geodesic is γ(T/2), not γ(0). Take any continuous function λ : M → R. The attenuated x-ray transform with attenuation λ is defined via Notice that if the attenuation vanishes, then the attenuated x-ray transform is just the x-ray transform: I 0 = I . In particular, any function f ∈ L 2 (M) is uniquely determined by its integrals over all geodesics.

Remark 30.
The theorem is stated in a setting where a small ball of (coordinate) radius R is removed from a spherically symmetric manifold. The theorem is also applicable in a whole ball without anything removed. One can artificially remove a small ball of radius R from the center of the manifold and use the theorem for all R > 0 to conclude that the attenuated x-ray transform is injective. In fact, it suffices that all regularity assumptions (on f, λ and c) are satisfied locally in the punctured ball-the functions may blow up at the origin.

Remark 31.
Theorem 29 can be seen as a support theorem. It follows immediately from the theorem that if a function on the spherically symmetric manifold integrates to zero over all geodesics that stay outside a ball centered at the (coordinate) origin, then the function must vanish outside that ball. In Euclidean geometry-which is a special case of the theorem-this is Helgason's famous support theorem [12, theorem 2.6].
Remark 32. Theorem 29 provides injectivity, but not stability. It is well known that a spherically symmetric metric satisfying the Herglotz condition can have conjugate points. In dimension two the existence of conjugate points immediately implies instability for the unweighted x-ray transform [25].

Proof of theorem 29. Fix a Lipschitz attenuation
satisfies I λ f = 0. We wish to show that f = 0; this implies injectivity.
If we intersect our manifold M with a two dimensional subspace of R n , we get a totally geodesic submanifold with radial symmetry. The restriction of a function f ∈ L 2 (M) is still L 2 for almost all of these submanifolds by lemma 17. If the theorem holds in dimension two, then f vanishes on almost all of these submanifolds and therefore on almost all of M. Therefore it suffices to prove the theorem in dimension two.
Suppose the result is true for C 1,1 wave speeds. Let R < a 1 < a 2 < · · · < a N < 1 be the points where c fails to be C 1,1 . Consider the submanifold M ⊂ M where the radial coordinate is restricted to (a N , 1]. On M the metric is C 1,1 and satisfies the assumptions. Applying the result on the manifold M shows that f | M = 0.
Then we may restrict our problem to M \ M , since f vanishes elsewhere. In the same way we can deduce that f must vanish also between radii a N−1 and a N . Continuing this layer stripping argument at discontinuities of c eventually shows that f must vanish in all of M.
Therefore it suffices to prove the theorem when c ∈ C 1,1 . It is again convenient to use polar coordinates. We may write f (r, θ) as a Fourier series in θ. By lemma 20 there is a sequence of functions a k ∈ L 2 (R, 1) so that where the series converges in L 2 (M). Let us denote f k (r, θ) = a k (r)e ikθ . For any (r, θ) ∈ M there are two geodesics with their tip at this point. They are reverses of each other. If these geodesics are denoted by γ ± (r, θ), we write Combining equations (84) and (88), we have See section 5.3 for details. It is easy to check that E λ (r) is bounded, and lemma 25 guarantees that the integral transforms A λ k : L 2 (R, 1) → L 2 (R, 1) are equicontinuous. Therefore I λ 0 : L 2 (M) → L 2 (M) is continuous and Since E λ (r) never vanishes and we assumed that I λ f = 0, we have A λ k a k = 0 for every k ∈ Z. The Abel-type integral transform A λ k is injective by lemma 25, so in fact a k = 0 for all k. This means that f indeed vanishes identically. □

The broken ray transform
Injectivity results for the broken ray transform in the Euclidean disc or ball were shown in [16]. In this section we generalize those results to radially symmetric manifolds satisfying the Herglotz condition and the countable conjugacy condition. The Euclidean metric satisfies both conditions. We will only give abridged versions of the proofs, focusing on the new ingredients due to a non-Euclidean metric. The missing details can be found in the Euclidean arguments presented in [16].
By the boundary ∂M of our spherically symmetric manifold we only mean the outer boundary where r = 1. Let us choose a subset E ⊂ ∂M . The set E is the set of tomography, where broken rays have their initial and final points, and the set ∂M \ E is where the broken rays reflect. All measurements are therefore done on the set of tomography, and the rest of the boundary is essentially a mirror.
The broken ray transform takes a function on M into a function on the set of broken rays. Injectivity of this transform depends on the set of tomography E. The bigger E is, the more data are available.

Remark 35.
It is an open problem whether the broken ray transform is injective on smooth functions in the Euclidean unit disc.
Proof of theorem 33. It suffices to prove the theorem with R > 0; the results for R = 0 follow from combining the results for all positive R. It is enough to prove both statements in dimension two.
(1) All broken rays are in fact periodic broken rays, since the initial and final points must coincide. The set of radii corresponding to such broken rays is dense by proposition 28. In fact, for any N ∈ N, the set of radii of broken rays with at least N reflections is dense. This can be seen by induction, as excluding one number of reflections does not affect the density in the absence of conjugate points. Fix any z ∈ (R, 1). There is a sequence of broken rays γ i so that the number of reflections on γ i tends to infinity as i increases, γ i meets E only at its endpoints and the smallest radial coordinate on γ i tends to z. Using the aforementioned density result, one can find a sequence of periodic broken rays whose radius converges to z and whose number of reflections increases without bound.
The average integral of f over γ i tends to where A 0 is the Abel-type transform defined in (75), 1 denotes the constant function, and a 0 is the radially symmetrized f (as in lemma 20). In essence, the idea is that in the limit i → ∞ the normalized integral over the broken ray γ i tends to a rotation symmetric measure corresponding to A 0 . For details of this calculation, see [16, section 2.3].
It is easy to check that A 0 1(z) > 0 for z < 1, so the data determine the function A 0 a 0 . The transform A 0 is injective by lemma 25, so the data determine a 0 . This is exactly the claim.
(2) Suppose f integrates to zero over all broken rays. We may write it as a Fourier series as in lemma 20, and the series converges uniformly (see [16, lemma 8]). By the first part of the theorem we already know that a 0 = 0.
Since f integrates to zero over all broken rays, so do its angular derivatives of all orders. We shall use angular derivatives of even orders, ∂ 2n θ f . We may choose our polar coordinates (r, θ) so that the point (1,0) is in E. Consider a broken ray γ whose initial point is at θ = −ω and final point at θ = ω for some small ω > 0. Let z γ denote the distance from γ to the origin in R 2 . Then the average of ∂ 2n θ f over the broken ray γ is Here S k (γ) is a coefficient depending on the broken ray defined in [16, section 4.1]. It has exactly the same form in the non-Euclidean setting.
To find formula (103), one first needs to recall that if f is written as an angular Fourier series then the integral of f over a geodesic with a tip at (r, θ) is The Euclidean version of this statement was given (implicitly) in [16, equation (37)] and the more general statement in radial symmetry follows form (73). Suppose then that the corresponding geodesic has opening angle α(r) ∈ πQ and the corresponding periodic broken ray has N reflections. The integral over the boken ray is then N l=1 If (r, θ + lα).
This sum over l can be done explicitly for each term e ikθ , and this leads to the coefficients S k (γ). This argument is the same in Euclidean geometry; the only difference is in the definition of the Abel transforms A k . Since the sum (103) vanishes for all n ∈ N, then S k (γ)A k a k (z γ ) = 0 for all k ∈ Z. The assumption of quasianalyticity is needed for this step. The coefficient S k (γ) is non-zero whenever the broken ray γ cannot be extended to a periodic one-or equivalently, whenever ω / ∈ πQ. The set of corresponding z γ is dense in (R, 1) by proposition 28 (the complement is countable). Thus each of the functions A k a k , k ∈ Z, vanishes in a dense set.
We then turn to lemma 25 for properties of Abel transforms. The functions a k are continuous, and so are A k a k . The transforms are injective, so all of the functions a k must vanish. This concludes the proof. □

The periodic broken ray transform
In this section we consider the periodic broken ray transform. The related inverse problem is to reconstruct a function on a manifold with boundary from its integrals over all periodic broken rays. It turns out that in spherical symmetry the transform does not contain enough information to recover the entire function. In dimension three or higher one can recover the even part but very little information about the odd part. In dimension two one can recover the circlewise average (like in theorem 33) but very little other information.
The periodic broken ray transform was previously known to be injective in the Euclidean square but non-injective in the Euclidean disc [19]. If Y : S n−1 → C is a finite sum of odd spherical harmonics, the space of such functions a ∈ C ∞ 0 (R, 1) that the function M rω → Y(ω)a(r) integrates to zero over all periodic broken rays has finite codimension in C ∞ 0 (R, 1). In particular, the kernel of the periodic broken ray transform is infinite dimensional (but contained in the space of odd functions if n 3).

Remark 37.
We may without loss of generality assume that R > 0 in the proof, and we do so without mention. This choice eliminates radial geodesics. Letting R → 0 easily gives the result for R = 0 once it has been proven for all R > 0.

Remark 38.
A version of this problem is considered in [5] in connection with spectral rigidity problems. For rigidity within a radially symmetric class, one only needs to consider radially symmetric functions, and for those we do indeed have injectivity for n 2.
The tools needed for this proof can be found in the subsections below. Proofs of different parts of the theorem are given in sections 8.2 and 8.3. Each part of the theorem is given a separate proof.

The planar average ray transform
Let G n 2 be the Grassmannian of two-dimensional subspaces of R n . Notice that every non-radial geodesic is contained in exactly one plane in G n 2 and that M ∩ P is a totally geodesic submanifold of M for any P ∈ G n 2 . These submanifolds inherit properties from M: the Herglotz condition, the countable conjugacy condition, absence of conjugate points, regularity, and also the negations of these conditions.
Let us now define the planar average ray transform P , an integral transform closely related to the periodic broken ray transform. The planar average ray transform Pf of a sufficiently regular function f : M → C is a function on G n 2 × (R, 1). If If (γ) denotes the integral of f over a geodesic γ, then Pf (P, r) is the average of If (γ) over all geodesics γ of radius r and contained in the plane P. If n = 2, the set G n 2 is a singleton and we consider Pf to be a function of radius only.
In dimension two, let f 0 denote the angular average of a function f : M → C, that is, In the lemmas below, A 0 is the Abel transform corresponding to the conformal factor c(r) defined by equation (75). Proof. Fix any r ∈ (R, 1). We have for any angle θ ∈ S 1 ; see equation (70). If we take the integral average over θ ∈ S 1 , the lefthand side becomes Pf (r) and the right-hand side becomes A 0 f 0 (r). □ Lemma 40. Assume the assumptions of theorem 36. Let n = 2 and fix an exponent p ∈ [1, ∞). For almost every r ∈ (R, 1) we have Pf (r) = A 0 f 0 (r). If p > 2, this holds for every r.
Proof. We wish to show that for almost every r ∈ (R, 1) (every r if p > 2) our function f ∈ L p (M) is integrable over almost every geodesic of radius r, and that Pf = A 0 f 0 . Let us define the operator F : 1) is clearly continuous, and so is A 0 : L p (R, 1) → L p (R, 1) (by lemma 25).
Since continuous functions are dense, we can take a sequence of functions f k ∈ C(M) tending to f in L p (M). Since F |f k | → F |f | in L p (R, 1), we can-after passing to a subsequenceassume that F |f k | (r) → F |f | (r) and f k (r) → f (r) for almost every r ∈ (R, 1). If p > 2, we have F |f k | (r) → F |f | (r) for every r.
Take any such r. By lemma 39 we have the bound − S1 γ r,θ Due to our choice of r, this is a bounded sequence in k. Since f k → f almost everywhere, we have by the dominated convergence theorem that exists and is finite. This implies that Pf (r) makes sense and actually Pf (r) = A 0 f 0 (r). □ In the planar average ray transform we may restrict the set of radii to a subset of (R, 1), as we will do next. The following theorem characterizes the kernel of the planar average ray transform in dimensions n 3 as the set of odd functions. Theorem 41. Let R ∈ [0, 1). Let M =B(0, 1) \B(0, R) ⊂ R n , n 3, be equipped with a metric corresponding to a radially symmetric C 1,1 wave speed that satisfies the Herglotz condition (definition 1). Fix a set E ⊂ (R, 1) and consider the planar average ray transform P restricted to radii in E.
(1) If E has full measure, then P is well defined and determines the even part on L p for any p 1. (2) If E is dense, then P is well defined and determines the even part on L p for any p > 2.
(3) The planar average ray transform P vanishes on all odd functions.
That is, the kernel of P is precisely the set of odd functions (under either of the assumptions on E and p given above).
Proof. The last conclusion follows immediately from symmetry considerations, so we focus on the first two.
Let f ∈ L p (M) and suppose Pf (P, r) = 0 for every r ∈ E and almost every P ∈ G n 2 . We wish to show that f has to be odd.
For almost every plane P ∈ G n 2 the function f P = f | P∩M belongs to L p (P ∩ M) due to lemma 17. The planar average ray transform of f P can be considered a function of radius only. By lemma 40 we have A 0 f P 0 (r) = Pf P (r) = Pf (P, r) = 0 for all r ∈ E. We have f P 0 ∈ L p (R, 1). Let us consider the two remaining cases of the theorem: (1) If p 1, we know that A 0 f P 0 ∈ L 1 (R, 1) (see lemma 25). Since E has full measure, we must have A 0 f P 0 = 0 (as an elment of L 1 ). (2) If p > 2, we know that A 0 f P 0 ∈ C(R, 1) (see lemma 25). Since E is dense, we have A 0 f P 0 = 0.
Either way, the function A 0 f P 0 vanishes. Now injectivity of the Abel transform-guaranteed by lemma 25-implies that f P 0 (r) = 0 for almost every r ∈ (R, 1). This means that for almost every r ∈ (R, 1) the function f | ∂B(0,r) is in the kernel of the Funk transform, identifying spheres of different radii by scaling. The function f is in L 1 for almost every sphere, so by theorem 23 (with k = 1) the function f must be odd on almost every sphere. This concludes the proof. □

Positive results
We are now ready to start proving theorem 36. We give the positive results here and the negative ones in the next subsection.
Proof of theorem 36, part 1. Consider f ∈ L p (M) with vanishing periodic broken ray transform. Let P ⊂ (R, 1) be the set of radii corresponding to periodic broken rays. By averaging over rotations, we observe that the periodic broken ray transform determines the planar average ray transform with radii restricted to the set P. The set P is dense by proposition 28.
It follows from lemma 40 that A 0 f 0 (r) = Pf (r) for every r ∈ (R, 1). Since Pf vanishes in the dense set P and A 0 f 0 is continuous by lemma 25, we know that A 0 f 0 (r) = 0 for every r ∈ (R, 1). It follows from injectivity of A 0 (due to lemma 25) that f 0 = 0. □ Proof of theorem 36, part 2. Suppose the periodic broken ray transform of f ∈ L p (M) is known. We begin as in the previous proof and observe that the planar average ray transform of f is thus determined for radii in P. Theorem 41 allows us to reconstruct the even part of f from this information since P is dense. □

Negative results
It now remains to prove the non-injectivity results contained in theorem 36. In both claims we additionally assume the finite conjugacy condition. It follows from this condition that preimages of singletons under the function α defined by (69) are finite. To see this, observe that by the mean value theorem α has to vanish between any two points with equal value of α, but the assumed condition implies that α has finitely many zeros. A similar statement is true in the countable case, see [5, lemma 4.5].
Proof of theorem 36, part 3. The angle α(r) was discussed in section 5.1, and it was observed in section 5.5 that a trajectory of radius r is periodic if and only if α(r) ∈ πQ. Let us now calculate the integrals of functions over these trajectories one Fourier component at a time.
Suppose f (r, θ) = a(r)e ikθ for some continuous function a. Due to equation (73) the integral of f over the geodesic parametrized by (r, θ) is e ikθ A k a(r).
We parametrize periodic trajectories like the geodesics that constitute them. This parametrization is redundant-each trajectory is described as many times as it has geodesic segmentsbut it does not matter. We say that a radius r ∈ (R, 1) has index m ∈ N if the corresponding billiard trajectory is periodic and has m reflections. The index m of a radius r is the smallest natural number so that mα(r) ∈ πN.
Let r have index m and take any θ ∈ S 1 . Then the integral of f over the trajectory given by Here A k is the Abel-type integral transform introduced in sections 5.2 and 5.4. If kα(r) ∈ πN, then the sum is simply m. Otherwise it is zero, since e 2ikmα(r) = 1. The condition kα(r) ∈ πN is equivalent with m dividing k, so the integral is Gf (r, θ) = mδ m|k e ikθ A k a(r), where we have denoted Notice that if k is kept fixed and non-zero, this integral can only be nonzero for a finite amount of indices m. By standard convention every integer divides zero.
For an integer m 1, let us denote by Q m ⊂ (R, 1) the set of radii of index m. Each set Q m is finite, but the union m∈N Q m is dense by proposition 28.
By assumption for some natural number K and some coefficients y k . The assumption of zero average translates to y 0 = 0.
For a ∈ C ∞ 0 (R, 1), let us denote f a (r, θ) = a(r)Y(θ). For a radius r of index m we have by the above considerations In particular, Gf a (r, θ) can only be nonzero if Therefore the space E = {a ∈ C ∞ 0 (R, 1); A k a(r) = 0 for all 1 k K and r ∈ Q K } (116) is a subspace of the space we set out to study and it suffices to show that E has finite codimension in C ∞ 0 (R, 1). But since Q K is finite, the elements of E only satisfy a finite number of linear conditions. Therefore the codimension of E is indeed finite. □ Proof of theorem 36, part 4. This follows from part 3. We only need to know that when an odd spherical harmonic is restricted to a two dimensional subspace (which after intersecting with the sphere means a great circle) is an odd trigonometric polynomial with order bounded by that of the spherical harmonic. To see this, one can write the spherical harmonic as a harmonic polynomial of the same order as the spherical harmonic; its restriction to a subspace is a polynomial of at most the same order, and the order of a polynomial as a spherical harmonic is bounded by the order of the polynomial. □