Oscillation tomografy study of Earth's composition and density with atmospheric neutrinos

Knowledge of the composition of the Earth's interior is highly relevant to many geophysical and geochemical problems. Neutrino oscillations are modified in a non-trivial way by the matter effects and can provide valuable and unique information not only on the density but also on the chemical and isotopic composition of the deep regions of the planet. In this paper, we re-examine the possibility of performing an oscillation tomography of the Earth with atmospheric neutrinos and antineutrinos to obtain information on the composition and density of the outer core and the mantle, complementary to that obtained by geophysical methods. Particular attention is paid to the D$^{\prime \prime}$ layer just above the core-mantle boundary and to the water (hydrogen) content in the mantle transition zone. Our analysis is based on a Monte-Carlo simulation of the energy and azimuthal angle distribution of $\mu$-like events generated by neutrinos. Taking as reference a model of the Earth consisting of 55 concentric layers with constant densities determined from the PREM, we evaluate the effect on the number of events due to changes in the composition and density of the outer core and the mantle. To examine the capacity of a detector like ORCA to resolve such variations, we construct regions in planes of two of these quantities where the statistical significance of the discrepancies between the reference and the modified Earth are less than $1\sigma$. The variations are implemented in such a way that the constraint imposed by both the total mass of the Earth and its moment of inertia are verified.


Introduction
Aside from its intrinsic interest, a detailed description of the inner parts of the Earth is essential for a proper understanding of basic geological phenomena such as volcanology, earthquakes, plate tectonics, and mountain building [1,2]. According to current knowledge, the Earth's internal structure is stratified and consist of successive layers with different chemical, geological, and physical properties. Most of the information about the layers and the boundaries between them has been acquired by examining how seismic waves created by natural a e-mail: dolivo@nucleares.unam.mx b e-mail: fis_pp@ciencias.unam.mx c e-mail: ismaelromero_@hotmail.com d e-mail: sampayo@mdp.edu.ar earthquakes are refracted and reflected as they propagate through the Earth [3,4]. Valuable information has also been obtained from measurements of magnetic and gravitational fields, observations of the planet's moment of inertia and precessional motion, and physical, chemical, and mineralogical analyzes of meteorites and xenoliths. Based on chemical composition, three main layers have been identified within our planet: the crust, the mantle, and the core. The mantle and the core are further subdivided into two regions each.
The mantle surrounds the core and extends from beneath the thin crust to a depth of 2,900 km. It consists mainly of dense silicate rocks rich in iron and magnesium. Seismic wave velocities in the mantle show three discontinuities at depths of 410, 660 and 2,700 km [5]. The first two correspond to the edges of the transition zone between the upper and lower mantle, and are best explained by mineral phase transformations without compositional changes. The third discontinuity is typically a 2.5-3.0% increase for both S-and P-waves observed at the top of the D layer, a transition shell about 200-300 km thick, at the base of the mantle, which presents a variety of seismic anomalies and is presumably the source of the large mantle plumes [6,7]. This layer is not yet fully understood, many of its characteristic may be attributed to the discovered MgSiO 3 postperovskite phase [8], but compositional differences may also play an important role in addition to the phase transformation [9]. Regarding the core, it is well established that it is composed primarily of an iron-nickel alloy, with Ni/Fe ∼ 0.06. The inner core is solid, while the lower density and absence of S-wave propagation are indications of a liquid outer core. The density deficit in the outer core cannot be simply explained by a difference in the state, but it requires about 5-10 wt% (weight percent) of lower atomic weight elements to reduce its density and melting point [10,11].
Good estimates of the abundance and distribution of "light" elements in the core are essential to understanding the formation and evolution of the Earth, as well as how the core and mantle interact in the region around the core-mantle boundary (CMB) [9,12]. The various processes that occur in this interface are highly influenced by heterogeneous structures at or near it. The great disparity in density prevents direct convection movement through the CMB, and nearby regions control the transfer of heat and material [13]. These processes strongly affect the convection in the mantle, responsible for the plate motion and continental drift, and the more vigorous convective flow in the outer core that is believed to be at the origin of the Earth's magnetic field (i.e., the functioning of the geodynamo) [14,15]. Due to vigorous convection, the liquid part of the core is usually assumed to have a uniform composition. However, seismological evidence indicates the possible existence of a ∼ 300 km thick region on top of the outer core (E layer), which shows anomalous low seismic velocities. This region is likely to be less dense than the rest of the outer core, but simply increasing the concentration of light elements also produces higher speeds, in contradiction with observations. None of the mechanisms that have been considered are without complications [16] and further observations are needed before a satisfactory explanation of the origin and nature of the E' layer can be formulated.
The study of seismic wave propagation and normal mode oscillation is undoubtedly the most effective and reliable method to search the Earth's deep structures and process. Nevertheless, the impressive progress done has not been accompanied with a concomitant improvement in the precision of the density estimations [17]. They are done through an underdetermined inversion problem, performed in two steps: the spatial distribution of seismic wave velocities is first inferred from seismological data and then the density distribution is inferred from the seismic velocities using some empirical relation [18]. Such procedure allows the average density along a path to be estimated with an uncertainties of about 5% [17] for the mantle and presumably larger for the core. For example, the density jump at the inner-core boundary, which play an important role in the maintenance of the geodynamo, has been inferred to be of 0.82 with an error of more than 20% [19]. While the density distribution can be obtained from seismological remote sensing, the compositional structure of the Earth [20,21] has been much more difficult to determine. Thus, the compositions of the lower mantle and the core remain quite uncertain, despite significant advances in recent years. Since in situ sampling is impossible, estimations are done by comparing density and sound velocity data from seismological observations with those from laboratory experiments and theoretical calculations [22]. Due to technical limitations, it has been difficult to perform reliable experiments for molten samples under high pressure and temperature conditions and the available information is still insufficient to infer the core composition. The most likely light elements in the core are oxygen, silicon, carbon, nitrogen, sulfur, and hydrogen [9,23], but there is still no consensus on the nature and proportion of the components.
In addition to the metallic core, significant amounts of hydrogen can be incorporated within the mantle, in silicate minerals, melts, and hydrous fluids. This is done in a variety of chemical species (OH, H 2 , · · · ) generically referred to as "water" [24][25][26]. The abundance and distribution of water has influenced not only the evolution, dynamics state and thermal structure of the deep interior, but also the evolution of the crust and hydrosphere [27]. Even small amount of water can affect properties like melting temperature, rheological properties, electrical conductivity, and seismic velocities of the mantle. How water is transported into Earth's deep interior and how it is distributed are today open questions. Almost no constraints exists on the water content of the lower mantle. The existing laboratory data cannot be used to infer the water content from geophysical observations at these high pressure and temperature conditions. Low-velocity regions have been observed near the top of the lower mantle and directly above the CMB [25]. If the low-velocity is interpreted as caused by partial melting, then some water is likely be present in these regions [28]. Water could instead be stored in the lowest parts of the mantle in Al-postperovskite [29]. There is a wide consensus that the mantle transition zone (MTZ), at 410-660 km deep, is a potential reservoir of water because its main mineral constituents can store up to ∼ 3% wt water [27]. However, the amount of water contained there is poorly constrained. Some water-rich inclusions recently found in diamonds [30,31] suggest a wet MTZ, but it is not clear if they are representative of the typical water content of the deep mantle or reflect local conditions. Geophysical methods (electrical conductivity and seismic observations) provides constraints on the water distribution in a global scale [32]. The inferences with these approaches are not direct and, despite the effort involved, a wide range of values have been reported in the literature [33,34]. As noticed in Ref. [33], geophysical estimates have large uncertainties and the water content in the MTZ can be heterogeneous, having significant lateral variability as revealed by a recent novel approach [35].
From the comments in the previous paragraph it is apparent how useful it would be to have additional experimental techniques, unaffected by the same uncertainties, which could provide complementary and independent information about the deep interior of the Earth. A promising candidate is neutrino tomography [36]. The basic idea is that neutrino propagation within the Earth is affected by their interactions with the particles present in the terrestrial matter. The cross section for neutrino interactions increase with energy and, in the case of absorption tomography, the density profile can be reconstructed from the attenuation of the flux of very high energy ( 10 TeV) neutrinos passing through the Earth [37 -44]. Another option is oscillation tomography, which takes advantage of the matter effect on flavor oscillations [45][46][47] of lower energy (MeV to GeV) neutrinos. In a medium, the transition probabilities between active neutrinos depend on the number density of electrons n e along the neutrino trajectory, which is proportional to the product of the matter density ρ times the average ratio Z/A of the atomic number Z and the mass number A [36,[48][49][50][51][52][53][54][55][56][57][58].
In this paper, we re-examine the feasibility of studying the internal structure of the Earth using atmospheric neutrino oscillation tomography. We analyze the ability of a detector such as ORCA to resolve deviations both of density and composition with respect to a standard Earth modeled in terms of 55 concentric shells integrated in five main layers, corresponding to the inner and outer core, the lower and upper mantle, and the crust. Constant densities are assigned to the shells from the mean value of the PREM densities within each shell. In our scheme, unlike other work on the subject, the densities of the main layers can be modified in a manner consistent with the well-measured total mass and moment of inertia of the Earth. In contrast, Z/A is a function of the chemical and isotopic composition of the medium and is not subject to either of these constraints. Therefore, in principle, one could constrain the allowed values of the density and composition of the Earth's deeper regions by studying the effects that changes in these quantities have on the events produced by atmospheric neutrinos after traversing the Earth, holding the total mass and moment of inertia of the planet fixed. We focus on the possible application to obtain information on the composition and density of the lower mantle regions above the CMB, in particular, the D region, paying special attention to the content of light elements, more specifically hydrogen. With the exception of hydrogen ((Z/A) H = 1), all other light elements have an almost equal number of protons and neutrons and hence Z/A ∼ = 0.5 for them. Thus, the presence of a significant amount of hydrogen (water) would produce appreciable Z/A changes compared to those in dry regions. In this sense, information from neutrino oscillation tomography could give valuable information on the water content in the lower mantle and core. We also allow for some variation in the location of the boundary between the lower and upper mantle as an effective way to account for the transition region between these two layers.
The paper is organized as follows. In Sect. 2 we present a model of the Earth's structure with 55 shells integrated within 5 main layers. In Sect. 3 we briefly review the formalism of matter neutrino oscillations and describe the algorithm to calculate the transition probabilities. In Sect. 4 we determine the number of µ-like neutrino events in a detector such as ORCA and the effects that changes in the composition and density of the outer core and lower mantle have on this observable. The results and final comments are presented in Sect.5, where we carry out a Monte Carlo simulation of the number of µ-like events and apply it to test different composition models of the outer core and mantle, changing also the density.

Model of the Earth's structure
The observed lateral variations in the Earth's properties are much less pronounced than the vertical variations. Therefore, the internal structure of the Earth can be well approximated by one-dimensional spherically symmetric models of seismic velocities, attenuation, and density as a function of depth [59,60].The most widely used of these models for seismic tomography has been the Preliminary Reference Earth Model (PREM) [59]. This model represents the mean properties of the Earth as a function of the radial distance r and was designed to fit different data sets and some basic data of the planet (radius, mass, and moment of inertia). In this study, we use a spherical model consisting of 55 concentric shells, each with a constant density equal to the average value of the PREM densities in the shell. The set of shells is divided into five large layers demarcated by concentric spheres of different radii: inner core (IC), outer core (OC), lower mantle (M 1 ), upper mantle (M 2 ), and crust (C). In Table 1 we give the values for the standard composition and the corresponding radii.  The primary information on the Earth's density as a function of r comes from the total mass of the Earth M ⊕ = 5.9724 × 10 27 g and its mean moment of inertia about the polar axis I ⊕ = 8.025 × 10 44 gcm 2 [61]. From these two quantities and R ⊕ = 6371 km for the Earth's radius, one gets I ⊕ 0.33M ⊕ R 2 ⊕ that is noticeably smaller than the moment of inertia of a homogeneous sphere of the same radius (0.40M ⊕ R 2 ⊕ ). This corroborates that there must be a concentration of mass towards the center of the planet or, in other words, that the inner regions are denser than average. We require that our model satisfy both constraints: where M i and I i , i = IC, OC, M 1 , M 2 ,C, are the masses and moments of inertia of the major layer specified above, which are given by the sums of the contributions of the shells contained in each of these divisions, as indicated in Table 1.
To modify the densities of the outer core and the lower and upper mantle we multiply the densities of all shells within each of these layers by the respective rescaling factor, f OC , f M 1 , and f M 2 . This is done in such a way that neither M ⊕ nor I ⊕ change. Then, Equating Eqs. (1) and (2) we obtain the following homogeneous system of linear equations: where are the relative changes of the densities in the outer core, lower mantle, and upper mantle, respectively. Solving this system, we can express δ M 1 and δ M 2 as functions of δ OC : where The value of the radius R M 1 set the position of the boundary between the regions M 1 and M 2 and varying it we can change the number of shells within each of these layers. This, in turn, modifies the values of M M 1,2 and I I 1,2 and makes the quantities δ 1,2 dependent on R M 1 . Fig. 2 shows the relative changes in the densities of layers M 1 and M 2 as a function of the relative change in the density of the outer core, for three different positions of the boundary between M 1 and M 2 .

Atmospheric neutrino oscillations
Neutrino oscillations are a well-verified and widely studied phenomenon that has proven beyond any doubt that neutrinos are mixed massive particles. The current experimental and observational data set can be interpreted in terms of the minimal extension of the Standard Model, where the known flavor states |ν α (α = e, µ, τ) are linear combinations of the states where c i j = cos θ i j and s i j = sin θ i j . A value of δ different form 0 or π implies CP-violation in the leptonic sector of the theory 1 .
In addition to the mixing angles θ i j and the CP-violating phase δ , the oscillations between the three active neutrinos are parametrized by two squared-mass differences: Five of the parameters (θ 12 , θ 13 , θ 23 , ∆ m 2 21 , and |∆ m 2 31 |), have been determined with remarkable precision (∼ 1 − 5%) by global fits of the data from solar, atmospheric, reactor and long baseline experiments. The issues still pending are: the sign of ∆ m 2 31 , the θ 23 octant, and the determination orfthe phase δ . The sign of ∆ m 2 31 characterize the normal ordering (NO), with m 3 > m 1,2 , and the inverted ordering (IO), with m 3 < m 1,2 . The 3ν oscillation parameters shown in Table 2 are the mean of the best-fit values for the allowed ranges at 1σ of the global analyses performed by three groups [62][63][64].
Neutrinos are produced and detected as flavor eigenstates. Consider a neutrino ν α produced at time t 0 that propagates in vacuum. Due to slight mass differences, the phases of the mass eigenstate components of the original flavor state change at different rates and, due to this, the flavor content of the neutrino beam oscillates along the trajectory. When neutrinos propagate in a medium, the coherent forward scattering of neutrinos with electrons is different for ν e and ν µ,τ , resulting in different refraction indexes for the electron neutrino and the other flavors. As a consequence, neutrino oscillations can be significantly modified in matter compared to oscillations in vacuum and new resonance enhancement effects appear. These effects are sensitive to the density and composition of the medium and we will take advantage of this in order to examine the inner parts of our planet by means of the oscillations of atmospheric neutrinos in the Earth. Atmospheric neutrinos have played a very important role in the study and characterization of the phenomena of neutrino oscillations. They are generated around the Earth as decay products in hadronic showers that result from collisions of cosmic rays with nuclei in the upper atmosphere. This provide a continuous source of neutrinos spanning a very wide range of energies and travelled distances before detection. On this work, we concentrate on those with energies in the range of 1-10 GeV and different nadir angles.
Let a neutrino ν α that enters the solid terrestrial matter at time t 0 . At any time t > t 0 the state of the system |ψ(t) can be expressed as |ψ(t) =Û (t,t 0 )|ψ(t 0 ) , where |ψ(t 0 ) = |ν α andÛ (t,t 0 ) is the evolution operator. The probability of having a neutrino of flavor β inside the Earth, at a distance t − t 0 ( = c = 1) from the entry point, is where the probability amplitude U β α ( ) = ν β |Û ( )|ν α is an element of the 3 × 3 unitary matrix U ( ) representing the evolution operator in the flavor basis. Rather than solving for U ( ) directly, it proves to be more convenient to determine the evolution operator in the basis of the mass eigenstates and then transform it to the flavor basis using the relation withH Here The first term in Eq. (10) is the Hamiltonian that drives the flavor evolution in vacuum, while the second term accounts for the matter effects. For antineutrinos, the sign of the second term is reversed and the matrix Γ is replaced by its complex conjugate.
In normal matter (n, p, e) where n e ( ) is the electron number density and G F is the Fermi constant. The electron number density depends on both the matter density ρ( ) and the chemical and isotopic composition of the medium: where m u = 931.494 MeV is the atomic mass unit and Z/A = ∑ λ λ (Z/A) λ . The summation runs over all the chemical elements present in the medium and (Z/A) λ denotes the ratio between the atomic number Z λ and the atomic mass A λ of the element that contributes a fraction λ to the mass at a given position. The relevant quantities for us are the oscillation probabilities for ν µ (ν µ ) → ν µ.τ (ν µ,τ ) and ν e (ν e ) → ν µ,τ (ν µ,τ ), which are required to compute the µ-like events produced in a detector by "upward" atmospheric neutrinos after traveling a distance L = 2R ⊕ cos η through the Earth (see Eq. (18)). According to the previous considerations, because of the dependence of these probabilities on the potential V( ), atmospheric neutrinos are sensitive to changes in the density and composition of the traversed layers. In what follows, we examine the feasibility of applying such an effect to obtain meaningful information about the deepest part of the Earth. With this in mind, we calculateŨ (L) for the Earth modeled as a sphere made up of 55 concentric spherical shells with different (constant) densities, as described in Sec. (2).
The complete evolution operator can be expressed as the product of the evolution operators for the consecutive shells through which the neutrinos pass on their way to the detector: where L J are the distances in each shell, such that L = ∑ j L j and The effective potential inH j takes the fixed value V j = √ 2G F n j e , where n j e denotes the constant number density of electrons in shell j. Since the Hamiltonians for the different layers do not generally commute between them, the exponential factors in Eq. (14) must be in the prescribed order. Each of these factors can be evaluated by applying the Cayley-Hamilton theorem which allows us to convert the infinity series into a polynomial:Ũ j (L j ) = a j 0 I + a j 1H + a j 2H 2 , where the coefficients are functions of the eigenvalues ofH j (which coincide with those of H j ) [55,66].   The distances L j are functions of the nadir angle η. According to Fig. 3, they can be determined as L j = |z j+1 −z j |, wherẽ To exemplify the effect that changes in composition and density have on the probabilities of flavor oscillations, in Figs. 4 and 5 we show the probabilities for ν µ → ν µ and ν e → ν µ as functions of the neutrino energy, for two different angles of incidence, both for NO and IO. In the next section, we use these probabilities in the calculation of the µ-like events produced by atmospheric neutrinos arriving at the detector after passing through the internal regions of the Earth.

Neutrino events and test of Earth's composition
Let us consider a generic detector with a mass of water or ice containing a number n N of nucleons as a target. Such a detector, like the planned ORCA, will efficiently detect the Cherenkov radiation emitted along the path of µ ± produced by the interactions of ν µ and ν µ with the nucleons in the instrumented volume or near it. In addition, it will be posible to detect showers due to the e or τ produced by the associated neutrinos [54,67,68].
To generate a data set that simulates the actual observations, with the inevitable statistical fluctuations, we perform a numerical calculation of the events due to neutrinos. Our observable is the number of µ-like events N µ and in what follows we use it to test different hypothesis about the composition and density of the outer core and the lower part of the mantle. To identify the most sensitive regions, we first do a scan dividing the cone under the detector into a series of angular and energy bins and calculate the number of µ-like events within given angular and energy intervals: where T is the detection time. The fluxes of atmospheric neutrinos and antineutrinos were taken from Ref. [69], while the charged-current cross sections for the ν µ (ν µ )nucleon and ν τ (ν τ )-nucleon scatterings, in the considered range of neutrino energies (1 − 10 GeV), are given approximately by [70]: In Eq. (18), Br τ→µ = Brτ →μ 0.17 are the branching ratios of the decays τ → µν µ ν τ and τ →μν µντ . The dependence of N µ on the density and composition of the medium is incorporated through the oscillation probabilities, which are calculated from the expressions given in Sec. 3. It is understood that these probabilities are evaluated at L. The values of the oscillation parameters are those given in Table 2. To estimate the impact of current uncertainties on the oscillation parameters, in addition to the best fits, we consider their values at the extremes of the 1σ ranges [62][63][64]. Thus, the muon numbers were also calculated by evaluating the oscillation probabilities at those values of a given parameter and keeping the best fits for the rest. We find that errors in the oscillation parameters introduce only few percent variations in the number of muon events and have little effect on our analysis.
To determine the angular and energy intervals where N µ is more sensitive to changes in the compositions of either the outer core or the lower mantle we introduce the quantity NO and IO. From the figures, it is apparent that in the case of the outer core the most sensitive regions correspond to energies around 5 GeV and angles compatible with the shadow of the outer core and M 1 . On the other hand, for the lower mantle the angular region has to be increased and the energy shifts to slightly lower values. The deviations in the number of events are considerably more pronounced for NO, mainly because for IO the resonance effects in matter occur for antineutrinos, whose charged cross-section is about two times smaller than that of neutrinos.
To simulate the number of events observed by a km 3 detector as ORCA, we follow the same procedure implemented in [55], which we reproduce here for completeness. In order  to perform the Monte Carlo simulation, based on the previous results, we consider events with 4 GeV < E < 11 GeV and 10°< η < 60°. Both of these intervals are divided into 200 subintervals. Every pseudo experiment is made up by tossing, in each square bin of the grid, a number of Poisson distributed events with the mean value equal to N µ as given by Eq. (18). Thus, each of the n exp pseudo experiments consists of 200 × 200 numbers corresponding to events, one for each bin. This sample of events are then distributed in angle and energy. We call them the true events and suppose that they are distributed according to the probability distribution function (pdf) f i exp t (E, η), i exp = 1, · · · , n exp . For this function we take the normalized histograms constructed by means of the Monte Carlo simulation.
The fractional number of true events in bin (i, j) for the i exp -th experiment is given by the integral of f i exp t (E, η) over the energy and angle intervals of the bin (i, j). To obtain a realistic distribution of events we must allow for a limited resolution of the detector, both in energy and angle. The net effect is the redistribution of the true events ("smearing") in the grid bins, which is implemented by folding the true distribution with a resolution function S (E o , η o |E, η). We assume a Gaussian smearing and write with the detector characterized by the angular and energy resolutions ∆ η(E) = α η / E/GeV and ∆ E(E) = α E E, respectively. For the values of the parameters α η and α E we consider different situations according to the discussion in Ref. [52]. To keep our analysis as simple as possible we assume a detection efficiency of 100%. When convoluted with the true events the kernel in Eq. 21 gives us what we call the observed events. That is, S (E o , η o |E, η) represents the conditional pdf for the measured values to be (E o , η o ) if the true values were (E, η). Since the event is observed somewhere, this function is normalized such that In terms of the resolution function, the number of observed events O i exp m,n in the bin (m, n) for the i exp -th experiment is given by As an illustration, in Fig. 8 and Fig. 9 we show the pdf of the true and observed events as functions of the energy and nadir angle, for NO and IO, respectively, in the case of a standard Earth. The figures were made using the same (large) number of bins for both kind of events, but these numbers generally differ. In what follows, to test how well the hypothesis of different Earth compositions is in agreement with the standard Earth, we take the observed events to be distributed into five angular bins and nine energy bins. From Eq. 23, for each  bin (m, n), with m = 1, · · · , 9 and n = 1, · · · , 5, we determine the number of events O si exp α for a standard Earth and O i exp α (H, κ) for an alternative Earth with different composition and density. Here α = 1, · · · , 45 label the two-dimensional bins (m, n). Thus, for each of these bins we have a sample of the observed events and determinate the corresponding mean valuesŌ For Poisson distributed events, in terms of the likelihood function L we construct the negative log-likelihood ratio function as According to Wilks' theorem [71] the χ 2 λ distribution can be approximated by the χ 2 distribution and, from it, the goodness of the fit can be established. The statistical significance of the χ 2 test is given, as usual, by the p-value: where n dof is the number of degrees of freedom and f χ (w, n dof ) is the chi-square distribution. In our case, n dof = 9 × 5 − 2 = 43. In this way, we can examine the levels of discrepancy between the standard Earth and different hypothesis about the composition and density of the outer core and the lower mantle. This is done in the next section, where the results and final comments are presented.

Results and Final Comments
Our goal is to determine to what extent the presence of light elements, in particular hydrogen, in the outer core and mantle can produce measurable effects due to the modifications it introduces in the flavor transformations of atmospheric neutrinos. We also pay attention to how uncertainties in the densities of the deepest regions of the planet can complicate obtaining compositional information. As discussed in Section 2, we rescale the mantle and outer core densities so that the constraints imposed by the Earth's total mass and moment of inertia are satisfied. At the same time, we allow for some changes in the compositions of the outer core and/or the lower mantle M 1 .
As the quantities to be fitted, we take the compositions of the outer-core and the lower mantle M 1 and the densities of the outer core, M 1 and the rest of the mantle M 2 . To evaluate the effects that changes of two of these quantities have on our observable, we construct regions in several planes: (δ OC , (Z/A) OC ), (δ M 1 , (Z/A) M 1 ), and ((Z/A) M 1 , (Z/A) OC ), where the statistical significance, given by the p-value, for the discrepancy between the standard and the modified Earth is less than 1σ . In Fig. (10) and (11) we show these regions for the oscillation parameters given in Table 2, for 10 years operation of an 8 Mton detector as ORCA, with resolution parameters α η = 0.25 and α E = 0.2. As can be seen, the regions for NO are more restricted than for IO. We consider two different values for the radius of the M 1 layer: 5600 km and 5871 km, which correspond to the lower and upper border of the transition zone between the lower and upper mantle.  In the graphs of composition versus density, non-zero correlations are observed, indicating that increases in composition are compensated by decreases in density, and vice versa, such that the electron density in the shell does not change. This compensation is not complete since, for the Earth's moment of inertia and mass to remain fixed, variations in density must also occur in other layers traversed by neutrinos and, as a consequence, the allowed regions become closed. The 1σ confidence regions for combined measures of Z/A in the outer core and M 1 agree with those obtained by [53]. Suppose that the variation of Z/A in the outer core or mantle is associated only with the abundance of hydrogen, then, as show in Fig.13, Z/A = (1 − H )(Z/A) 0 + H (Z/A) H , where (Z/A) 0 is the value with no hydrogen and H is the fractional contribution that hydrogen makes to the total mass of the layer. In this case, the 1σ regions are compatible with 5-8 wt% hydrogen in the lower mantle and outer core, respectively, values too high in light of geophysical estimates.  As a general rule, for the non fitted parameters we have kept their PREM values. However, for the sake of completeness, we paid some attention to effects from uncertainties in the densities when fitting the compositions. Thus, we have jointly fitted (Z/A) oc , (Z/A) M 1 and δ oc for the outer core density. (Notice, that the densities of the layers M 1 and M 2 vary in order to maintain the values of the Earth's mass and moment of inertia, as discussed in Sec-tion 2). The 1σ volume obtained in this way is projected onto the composition-composition plane, giving rise to the shaded area in Fig. 12. The observed negative correlation is caused by the compensation effect in matter neutrino oscillations, with an increase in composition counterbalanced by a decrease in density, and vice versa, to keep the electron number density unchanged. In the case of the composition-density plots (Fig.10a, Fig.10b, Fig.11a, and Fig.11b), the uncertainty in the outer core density is incorporated in the lower mantle density, since they are related by Eq. (4).  Fig.10 c, but including the uncertainties in the outer-core densities.
As is evident from Figs. 10 and 11, an experiment like ORCA has a limited potential to reveal a non-standard composition and/or density of the Earth. This worsens at 2 and 3σ and for the inverted ordering. A question then arises: how much exposure and how much resolution are required to constraint a non-standard composition? As a partial answer, in Fig.14 we shown the1σ regions for different exposures (in unity of ten years of ORCA operation) and angular and energy resolutions equal to α. We see that with an exposition of thirty years and resolutions of 0.1 it is possible to constraint the hydrogen content to 1%.
Finally, we pay special attention to the D region. Since the thickness of this remote interface between the rocky mantle and the iron core is relatively thin, changes in density and composition have little effect on our observable. Fig.15 shows the 1σ regions for the detector size versus the composition D , for a thickness of 200 km and various values of the detector resolution. As can be seen, the allowed regions are not very restrictive for reasonable values of resolution and size of the detector. The ability to constrain the composition is significantly improved in the case of a slightly thicker layer (500 km).
In summary, we have studied the possibility of conducting an oscillation tomography of the Earth based on the matter effects on the flavor oscillations of atmospheric neutrinos propagating through the depths of the planet. Using the µ-like events in a generic large Cherenkov detector as physical observables and making a Monte Carlo simulation of the energy and azimuthal angle distribution of these event, we tested possible variants with respect to a reference geophysical model with the densities as given by PREM and different composition in the outer-core and lower mantle. Unlike previous studies, the procedure we followed in this work allowed us to simultaneously vary the composition and density. When one of these quantities was fixed in the value of the reference geophysical model, our results, shown in Figs. 10, are compatible with those obtained by others authors in an uncorrelated way [50,53]. As shown, the study of questions of the type examined here would benefit greatly from the application of phenomena normally studied in the realm of particle physics and the advent of new and improved neutrino telescopes, such as KM3NeT, ORCA, PINGO and HyperK [72][73][74][75].