Search

As automated vehicles (AVs) become increasingly popular, the question arises as to how cyclists will interact with such vehicles. This study investigated (1) whether cyclists spontaneously notice if a vehicle is driverless, (2) how well they perform a driver-detection task when explicitly instructed, and (3) how they carry out such tasks. Using a Wizard-of-Oz method, 37 participants cycled a designated route and encountered an AV multiple times in two experimental sessions. In Session 1, participants cycled the route uninstructed, while in Session 2, they were instructed to verbally report whether they detected the presence or absence of a driver. Additionally, we recorded the participants' gaze behaviour with eye-tracking and their responses in post-session interviews. The interviews revealed that 30% of the cyclists spontaneously mentioned the absence of a driver (Session 1), and when instructed (Session 2), they detected the absence and presence of the driver with 93% accuracy. The eye-tracking data showed that cyclists looked more frequently and longer at the vehicle in Session 2 compared to Session 1. Furthermore, participants exhibited intermittent sampling of the vehicle, and they looked in front of the vehicle when it was far away and towards the windshield region when it was closer. The post-session interviews also indicated that participants were curious, felt safe, and reported a need to receive information about the AV's driving state. In conclusion, cyclists can detect the absence of a driver in the AV, and this detection may influence their perceptions of safety. Further research is needed to explore these findings in real-world traffic conditions.

The Tarantula Nebula in the Large Magellanic Cloud is known for its high star formation activity. At its center lies the young massive star cluster R136, providing a significant amount of the energy that makes the nebula shine so brightly at many wavelengths. Recently, young massive star clusters have been suggested to also efficiently produce high-energy cosmic rays, potentially beyond PeV energies. Here, we report the detection of very-high-energy $\gamma$-ray emission from the direction of R136 with the High Energy Stereoscopic System, achieved through a multicomponent, likelihood-based modeling of the data. This supports the hypothesis that R136 is indeed a very powerful cosmic-ray accelerator. Moreover, from the same analysis, we provide an updated measurement of the $\gamma$-ray emission from 30 Dor C, the only superbubble detected at TeV energies presently. The $\gamma$-ray luminosity above $0.5\,\mathrm{TeV}$ of both sources is $(2-3)\times 10^{35}\,\mathrm{erg}\,\mathrm{s}^{-1}$. This exceeds by more than a factor of 2 the luminosity of HESS J1646$-$458, which is associated with the most massive young star cluster in the Milky Way, Westerlund 1. Furthermore, the $\gamma$-ray emission from each source is extended with a significance of $>3\sigma$ and a Gaussian width of about $30\,\mathrm{pc}$. For 30 Dor C, a connection between the $\gamma$-ray emission and the nonthermal X-ray emission appears likely. Different interpretations of the $\gamma$-ray signal from R136 are discussed., Comment: 10+11 pages, 4+6 figures. Corresponding authors: L. Mohrmann, N. Komin

ULTRASAT (ULtra-violet TRansient Astronomy SATellite) is a wide-angle space telescope that will perform a deep time-resolved all-sky survey in the near-ultraviolet (NUV) spectrum. The science objectives are the detection of counterparts to short-lived transient astronomical events such as gravitational wave sources and supernovae. The mission is led by the Weizmann Institute of Science and is planned for launch in 2026 in collaboration with the Israeli Space Agency and NASA. DESY will provide the UV camera, composed by the detector assembly located in the telescope focal plane and the remote electronics unit. The camera is composed out of four back-metallized CMOS Image Sensors (CIS) manufactured in the 4T, dual gain Tower process. As part of the radiation qualification of the camera, Single Event Effect (SEE) testing has been performed by irradiating the sensor with heavy ions at the RADEF, Jyvaskyla facility. Preliminary results of both Single Event Upset (SEU) and Single Event Latch-up (SEL) occurrence rate in the sensor are presented. Additionally, an in-orbit SEE rate simulation has been performed in order to gain preliminary knowledge about the expected effect of SEE on the mission.

PSR B1259-63 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star. At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of ~100 hours of data taken from $t_p-24$ days to $t_p+127$ days around the system's 2021 periastron passage. We also present the timing and spectral analyses of the source. The VHE light curve in 2021 is consistent with the stacked light curve of all previous observations. Within the light curve, we report a VHE maximum at times coincident with the third X-ray peak first detected in the 2021 X-ray light curve. In the light curve -- although sparsely sampled in this time period -- we see no VHE enhancement during the second disc crossing. In addition, we see no correspondence to the 2021 GeV flare in the VHE light curve. The VHE spectrum obtained from the analysis of the 2021 dataset is best described by a power law of spectral index $\Gamma = 2.65 \pm 0.04_{\text{stat}}$ $\pm 0.04_{\text{sys}}$, a value consistent with the previous H.E.S.S. observations of the source. We report spectral variability with a difference of $\Delta \Gamma = 0.56 ~\pm~ 0.18_{\text{stat}}$ $~\pm~0.10_{\text{sys}}$ at 95% c.l., between sub-periods of the 2021 dataset. We also find a linear correlation between contemporaneous flux values of X-ray and TeV datasets, detected mainly after $t_p+25$ days, suggesting a change in the available energy for non-thermal radiation processes. We detect no significant correlation between GeV and TeV flux points, within the uncertainties of the measurements, from $\sim t_p-23$ days to $\sim t_p+126$ days. This suggests that the GeV and TeV emission originate from different electron populations., Comment: accepted to A&A

The spatial and spectral distributions of terahertz (THz) pulses emitted by two-color air plasmas are theoretically investigated for focused laser pulses and in the filamentation regime. We derive a so-called ''augmented'' conical emission model, which, similarly to the one originally proposed by You et al.\ [Phys.\ Rev.\ Lett.\ {\bf 109}, 183902 (2012)], involves phase matching between laser harmonics along the plasma channel, the plasma density and length, and the emitted frequency as key parameters. Our augmented model, however, treats envelope effects and accounts for transverse variations of both plasma- and Kerr-driven potential THz emitters. We highlight the importance of the characteristic spatio-spectral distributions of these two conversion mechanisms in the expression of the angular radiated power. The results of our model are successfully compared with data provided by a comprehensive, fully space and time-resolved unidirectional solver. Importantly, these numerical simulations clear up the effective plasma length along which THz emission develops, compared with the dephasing length along which the laser fundamental and second harmonic become out-of-phase. The impact of common optical aberrations, such as sphericity, astigmatism, and coma, on the THz generation is also investigated. Aberrations are shown to generally decrease the laser-to-THz conversion efficiency and potentially induce spatial asymmetries and narrowing in the THz spectra., Comment: 22 pages, 18 figures

The polarization states of terahertz (THz) radiation generated in a photo-ionized gas driven by strong two- or multi-frequency fields with locally controlled polarization are studied. We reveal a universal property of the resulting THz waveforms: the ellipticity of their polarization state increases linearly with the frequency. This ``linear chirp of ellipticity'' makes plasma-based THz generation unique among other THz sources. However, it also puts some constraints on the polarization properties of the generated THz radiation. We derive a general expression for the THz ellipticity and demonstrate how the polarization states of the generated THz waveforms can be manipulated and controlled by the polarization of the pump pulses.

We report on the development of thermal detectors based on large-size tellurium dioxide crystals (45x45x45 mm), containing tellurium enriched in $^{130}$Te to about 91%, for the CROSS double-beta decay experiment. A powder used for the crystals growth was additionally purified by the directional solidification method, resulting in the reduction of the concentration of impurities by a factor 10, to a few ppm of the total concentration of residual elements (the main impurity is Fe). The purest part of the ingot (the first ~200 mm, about 80% of the total length of the cylindrical part of the ingot) was determined by scanning segregation profiles of impurities and used for the $^{130}$TeO$_2$ powder production with no evidence of re-contamination. The crystal growth was verified with precursors produced from powder with natural Te isotopic composition, and two small-size (20x20x10 mm) samples were tested at a sea-level laboratory showing high bolometric and spectrometric performance together with acceptable $^{210}$Po content (below 10 mBq/kg). This growth method was then applied for the production of six large cubic $^{130}$TeO$_2$ crystals and 4 of them were taken randomly to be characterized at the Canfranc underground laboratory, in the CROSS-dedicated low-background cryogenic facility. Two $^{130}$TeO$_2$ samples were coated with a thin, $O$(100 nm), metal film in form of Al layer (on 4 sides) or AlPd grid (on a single side) to investigate the possibility to tag surface events by pulse-shape discrimination. Similarly to the small natural precursors, large-volume $^{130}$TeO$_2$ bolometers show high performance and even better internal purity ($^{210}$Po activity $\sim$ 1 mBq/kg, while activities of $^{228}$Th and $^{226}$Ra are below 0.01 mBq/kg), satisfying requirements for the CROSS and, potentially, next-generation experiments., Comment: Submitted to JINST; 22 pages, 11 figures, 4 tables

The CROSS experiment will search for neutrinoless double-beta decay using a specific mechanical structure to hold thermal detectors. The design of the structure was tuned to minimize the background contribution, keeping an optimal detector performance. A single module of the structure holds two scintillating bolometers (with a crystal size of 45x45x45 mm and a Ge slab facing the crystal's upper side) in the Cu frame, allowing for a modular construction of a large-scale array. Two designs are released: the initial $Thick$ version contains around 15% of Cu over the crystal mass (lithium molybdate, LMO), while this ratio is reduced to ~6% in a finer ($Slim$) design. Both designs were tested extensively at aboveground (IJCLab, France) and underground (LSC, Spain) laboratories. In particular, at LSC we used a pulse-tube-based CROSS facility to operate a 6-crystal array of LMOs enriched/depleted in $^{100}$Mo. The tested LMOs show high spectrometric performance in both designs; notably, the measured energy resolution is 5--7 keV FWHM at 2615 keV $\gamma$s, nearby the Q-value of $^{100}$Mo (3034 keV). Due to the absence of a reflective cavity around LMOs, a low scintillation signal is detected by Ge bolometers: ~0.3 keV (150 photons) for 1-MeV $\gamma$($\beta$) LMO-event. Despite that, an acceptable separation between $\alpha$ and $\gamma$($\beta$) events is achieved with most devices. The highest efficiency is reached with light detectors in the $Thick$ design thanks to a lower baseline noise width (0.05--0.09 keV RMS) when compared to that obtained in the $Slim$ version (0.10--0.35 keV RMS). Given the pivotal role of bolometric photodetectors for particle identification and random coincidences rejection, we will use the structure here described with upgraded light detectors, featuring thermal signal amplification via the Neganov-Trofimov-Luke effect, as also demonstrated in the present work., Comment: Submitted to JINST; 31 pages, 16 figures, 5 tables

The current experiments searching for neutrinoless double-$\beta$ ($0\nu\beta\beta$) decay also collect large statistics of Standard Model allowed two-neutrino double-$\beta$ ($2\nu\beta\beta$) decay events. These can be used to search for Beyond Standard Model (BSM) physics via $2\nu\beta\beta$ decay spectral distortions. $^{100}$Mo has a natural advantage due to its relatively short half-life, allowing higher $2\nu\beta\beta$ decay statistics at equal exposures compared to the other isotopes. We demonstrate the potential of the dual read-out bolometric technique exploiting a $^{100}$Mo exposure of 1.47 kg $\times$ y, acquired in the CUPID-Mo experiment at the Modane underground laboratory (France). We set limits on $0\nu\beta\beta$ decays with the emission of one or more Majorons, on $2\nu\beta\beta$ decay with Lorentz violation, and $2\nu\beta\beta$ decay with a sterile neutrino emission. In this analysis, we investigate the systematic uncertainty induced by modeling the $2\nu\beta\beta$ decay spectral shape parameterized through an improved model, an effect never considered before. This work motivates searches for BSM processes in the upcoming CUPID experiment, which will collect the largest amount of $2\nu\beta\beta$ decay events among the next-generation experiments.

HESS J1813$-$178 is a very-high-energy $\gamma$-ray source spatially coincident with the young and energetic pulsar PSR J1813$-$1749 and thought to be associated with its pulsar wind nebula (PWN). Recently, evidence for extended high-energy emission in the vicinity of the pulsar has been revealed in the Fermi Large Area Telescope (LAT) data. This motivates revisiting the HESS J1813$-$178 region, taking advantage of improved analysis methods and an extended data set. Using data taken by the High Energy Stereoscopic System (H.E.S.S.) experiment and the Fermi-LAT, we aim to describe the $\gamma$-ray emission in the region with a consistent model, to provide insights into its origin. We performed a likelihood-based analysis on 32 hours of H.E.S.S. data and 12 years of Fermi-LAT data and fit a spectro-morphological model to the combined datasets. These results allowed us to develop a physical model for the origin of the observed $\gamma$-ray emission in the region. In addition to the compact very-high-energy $\gamma$-ray emission centered on the pulsar, we find a significant yet previously undetected component along the Galactic plane. With Fermi-LAT data, we confirm extended high-energy emission consistent with the position and elongation of the extended emission observed with H.E.S.S. These results establish a consistent description of the emission in the region from GeV energies to several tens of TeV. This study suggests that HESS J1813$-$178 is associated with a $\gamma$-ray PWN powered by PSR J1813$-$1749. A possible origin of the extended emission component is inverse Compton emission from electrons and positrons that have escaped the confines of the pulsar and form a halo around the PWN., Comment: 13+5 pages, 13+11 figures. Accepted for publication in A&A. Corresponding authors: T.Wach, A.Mitchell, V.Joshi, P.Chamb\'ery

The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy $\gamma$ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula's $\gamma$-ray emission between 1 GeV and $\sim$100 TeV, that is, over five orders of magnitude in energy. Using the open-source software package Gammapy, we combined 11.4 yr of data from the Fermi Large Area Telescope and 80 h of High Energy Stereoscopic System (H.E.S.S.) data at the event level and provide a measurement of the spatial extension of the nebula and its energy spectrum. We find evidence for a shrinking of the nebula with increasing $\gamma$-ray energy. Furthermore, we fitted several phenomenological models to the measured data, finding that none of them can fully describe the spatial extension and the spectral energy distribution at the same time. Especially the extension measured at TeV energies appears too large when compared to the X-ray emission. Our measurements probe the structure of the magnetic field between the pulsar wind termination shock and the dust torus, and we conclude that the magnetic field strength decreases with increasing distance from the pulsar. We complement our study with a careful assessment of systematic uncertainties., Comment: 18+6 pages, 15+2 figures. Accepted for publication in A&A. Corresponding authors: M. Meyer, L. Mohrmann, T. Unbehaun. v2: after A&A language editing

Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin., Comment: 44 pages JCAP style (excluding author list and references), 19 figures; minor changes to match published version

We study the PSPACE-complete $k$-Canadian Traveller Problem, where a weighted graph $G=(V,E,\omega)$ with a source $s\in V$ and a target $t\in V$ are given. This problem also has a hidden input $E_* \subsetneq E$ of cardinality at most $k$ representing blocked edges. The objective is to travel from $s$ to $t$ with the minimum distance. At the beginning of the walk, the blockages $E_*$ are unknown: the traveller discovers that an edge is blocked when visiting one of its endpoints. Online algorithms, also called strategies, have been proposed for this problem and assessed with the competitive ratio, i.e. the ratio between the distance actually traversed by the traveller divided by the distance we would have traversed knowing the blockages in advance. Even though the optimal competitive ratio is $2k+1$ even on unit-weighted planar graphs of treewidth 2, we design a polynomial-time strategy achieving competitive ratio $9$ on unit-weighted outerplanar graphs. This value $9$ also stands as a lower bound for this family of graphs as we prove that, for any $\varepsilon > 0$, no strategy can achieve a competitive ratio $9-\varepsilon$. Finally, we show that it is not possible to achieve a constant competitive ratio (independent of $G$ and $k$) on weighted outerplanar graphs.

The radio galaxy M87 is a variable very-high energy (VHE) gamma-ray source, exhibiting three major flares reported in 2005, 2008, and 2010. Despite extensive studies, the origin of the VHE gamma-ray emission is yet to be understood. In this study, we investigate the VHE gamma-ray spectrum of M87 during states of high gamma-ray activity, utilizing 20.2$\,$ hours the H.E.S.S. observations. Our findings indicate a preference for a curved spectrum, characterized by a log-parabola model with extra-galactic background light (EBL) model above 0.3$\,$TeV at the 4$\sigma$ level, compared to a power-law spectrum with EBL. We investigate the degeneracy between the absorption feature and the EBL normalization and derive upper limits on EBL models mainly sensitive in the wavelength range 12.4$\,$$\mu$m - 40$\,$$\mu$m., Comment: 10 pages, 7 figures. Accepted for publication in A&A. Corresponding authors: Victor Barbosa Martins, Rahul Cecil, Iryna Lypova, Manuel Meyer, Perri Zilberman. Supplementary material: https://zenodo.org/records/10781524

BINGO is a project aiming to set the grounds for large-scale bolometric neutrinoless double-beta-decay experiments capable of investigating the effective Majorana neutrino mass at a few meV level. It focuses on developing innovative technologies (a detector assembly, cryogenic photodetectors and active veto) to achieve a very low background index, of the order of $10^{-5}$ counts/(keV kg yr) in the region of interest. The BINGO demonstrator, called MINI-BINGO, is designed to investigate the promising double-beta-decay isotopes $^{100}$Mo and $^{130}$Te and it will be composed of Li$_2$MoO$_4$ and TeO$_2$ crystals coupled to bolometric light detectors and surrounded by a Bi$_4$Ge$_3$O$_{12}$-based veto. This will allow us to reject a significant background in bolometers caused by surface contamination from $\alpha$-active radionuclides by means of light yield selection and to mitigate other sources of background, such as surface contamination from $\beta$-active radionuclides, external $\gamma$ radioactivity, and pile-up due to random coincidence of background events. This paper describes an R\&D program towards the BINGO goals, particularly focusing on the development of an innovative assembly designed to reduce the passive materials within the line of sight of the detectors, which is expected to be a dominant source of background in next-generation bolometric experiments. We present the performance of two prototype modules -- housing four cubic (4.5-cm side) Li$_2$MoO$_4$ crystals in total -- operated in the Canfranc underground laboratory in Spain within a facility developed for the CROSS double-beta-decay experiment., Comment: Submitted to Nucl. Instr. Meth. A; 28 pages, 11 figures, 4 tables

SS 433 is a microquasar, a stellar binary system with collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.), finding an energy-dependent shift in the apparent position of the gamma-ray emission of the parsec-scale jets. These observations trace the energetic electron population and indicate the gamma rays are produced by inverse-Compton scattering. Modelling of the energy-dependent gamma-ray morphology constrains the location of particle acceleration and requires an abrupt deceleration of the jet flow. We infer the presence of shocks on either side of the binary system at distances of 25 to 30 parsecs and conclude that self-collimation of the precessing jets forms the shocks, which then efficiently accelerate electrons., Comment: Submitted 20th Apr. 2023, published 25th January 2024 (accepted version)

We consider here the MultiBot problem for the scheduling and the resource parametrization of jobs related to the production or the transportation of different products inside a given time horizon. Those jobs must meet known in advance demands. The time horizon is divided into several discrete identical periods representing each the time needed to proceed a job. The objective is to find a parametrization and a schedule for the jobs in such a way they require as less resources as possible. Though this problem derived from the applicative context of reconfigurable robots, we focus here on fundamental issues. We show that the resulting strongly NP-hard Multibot problem may be handled in a greedy way with an approximation ratio of $\frac{4}{3}$., Comment: 26 pages, 4 figures

Shapley values are extensively used in explainable artificial intelligence (XAI) as a framework to explain predictions made by complex machine learning (ML) models. In this work, we focus on conditional Shapley values for predictive models fitted to tabular data and explain the prediction $f(\boldsymbol{x}^{*})$ for a single observation $\boldsymbol{x}^{*}$ at the time. Numerous Shapley value estimation methods have been proposed and empirically compared on an average basis in the XAI literature. However, less focus has been devoted to analyzing the precision of the Shapley value explanations on an individual basis. We extend our work in Olsen et al. (2023) by demonstrating and discussing that the explanations are systematically less precise for observations on the outer region of the training data distribution for all used estimation methods. This is expected from a statistical point of view, but to the best of our knowledge, it has not been systematically addressed in the Shapley value literature. This is crucial knowledge for Shapley values practitioners, who should be more careful in applying these observations' corresponding Shapley value explanations., Comment: Submitted to the 23rd European Young Statisticians Meeting (EYSM)

The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest for the faint signal of the primordial B-modes. A unique feature is the so-called "spectral imaging", i.e., the ability to recover the sky signal in several sub-bands within the physical band during data analysis. This feature provides an in-band spectral resolution of \Delta{\nu}/{\nu} \sim 0.04 that is unattainable by a traditional imager. This is a key tool for controlling the Galactic foregrounds contamination. In this paper, we describe the principles of bolometric interferometry, the current status of the QUBIC experiment and future prospects., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Gamma-ray observations have established energetic isolated pulsars as outstanding particle accelerators and antimatter factories in the Galaxy. There is, however, no consensus regarding the acceleration mechanisms and the radiative processes at play, nor the locations where these take place. The spectra of all observed gamma-ray pulsars to date show strong cutoffs or a break above energies of a few gigaelectronvolt (GeV). Using the H.E.S.S. array of Cherenkov telescopes, we discovered a novel radiation component emerging beyond this generic GeV cutoff in the Vela pulsar's broadband spectrum. The extension of gamma-ray pulsation energies up to at least 20 teraelectronvolts (TeV) shows that Vela pulsar can accelerate particles to Lorentz factors higher than $4\times10^7$. This is an order of magnitude larger than in the case of the Crab pulsar, the only other pulsar detected in the TeV energy range. Our results challenge the state-of-the-art models for high-energy emission of pulsars while providing a new probe, i.e. the energetic multi-TeV component, for constraining the acceleration and emission processes in their extreme energy limit., Comment: 38 pages, 6 figures. This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, Nat Astron (2023), and is available online at https://doi.org/10.1038/s41550-023-02052-3

Scalar-tensor theories with screening mechanisms come with non-linearities that make it difficult to study setups of complex geometry without resorting to numerical simulations. In this article, we use the $\textit{femtoscope}$ code that we introduced in a previous work in order to compute the fifth force arising in the chameleon model in the Earth orbit. We go beyond published works by introducing a departure from spherical symmetry $\unicode{x2014}$ embodied by a mountain on an otherwise spherical Earth $\unicode{x2014}$ as well as by implementing several atmospheric models, and quantify their combined effect on the chameleon field. Building on the numerical results thus obtained, we address the question of the detectability of a putative chameleon fifth force by means of space geodesy techniques and, for the first time, quantitatively assess the back-reaction created by the screening of a satellite itself. We find that although the fifth force has a supposedly measurable effect on the dynamics of an orbiting spacecraft, the imprecise knowledge of the mass distribution inside the Earth greatly curtails the constraining power of such space missions. Finally, we show how this degeneracy can be lifted when several measurements are performed at different altitudes., Comment: 36 pages, 20 figures. Accepted for publication in PRD