Search

The accuracy of atomic clocks is becoming small enough to determine altitude differences at the centimeter level in the near future. Current research activities consist in developing and comparing clock frequencies, and in setting up time transfer techniques with the required precision. Today, surveyor still has to provide clock altitude differences by traditional methods to allow these comparisons. This article reports levelling carried out by IGN at Paris Observatory for two European scientific projects in which the SYRTE Laboratory is involved., L'exactitude des horloges optiques permet d'envisager dans le futur leur utilisation pour déterminer des différences d'altitudes à une précision centimétrique. Les développements actuels consistent à comparer les fréquences d'horloges en cours de développement. Ainsi, la mise en place de protocoles de comparaison constitue une première étape vers cet objectif. Aujourd'hui, c'est encore au géomètre, au moyen des techniques traditionnelles, de déterminer l'altitude des horloges pour permettre ces comparaisons et fournir une mesure de référence. Cet article décrit les opérations de nivellement effectuées par l'IGN à l'Observatoire de Paris dans le contexte de projets scientifiques européens impliquant le laboratoire SYRTE de l'Observatoire de Paris. Abstract (english) : The accuracy of atomic clocks is becoming small enough to determine altitude differences at the centimeter level in the near future. Current research activities consist in developing and comparing clock frequencies, and in setting up time transfer techniques with the required precision. Today, surveyor still has to provide clock altitude differences by traditional methods to allow these comparisons. This article reports levelling carried out by IGN at Paris Observatory for two European scientific projects in which the SYRTE Laboratory is involved.

International audience

International audience

International audience; Current methods to determine the geopotential are mainly based on indirect approaches using gravimetric, gradiometric and topographic data. Satellite missions (GRACE, GOCE) have contributed significantly to improve the knowledge of the Earth's gravity field with a spatial resolution of about 90 km, but it is not enough to access, for example, to the geoid variation in hilly regions. While airborne and ground-based gravimeters provide the high resolution, the problem of these technics is that the accuracy is hampered by the heterogeneous coverage of gravity data (ground and offshore). Recent technological advances in atomic clocks are opening new perspectives in the determination of the geopotential. To date, the best of them reach a stability of 1.6×10-18 (NIST, RIKEN Univ. Tokyo) in just 7 hours of integration, an accuracy of 2.0×10-18 (JILA). Using the relation of the relativistic gravitational redshift, this corresponds to a determination of geopotential differences at the 0.1 m²/s² level (or 1 cm in geoid height). In this context, the present work aims at evaluating the contribution of optical atomic clocks for the determination of the geopotential at high spatial resolution. To do that, we have studied a test area surrounding the Massif Central in the middle of southern of France. This region, consists in low mountain ranges and plateaus, is interesting because, the gravitational field strength varies greatly from place to place at high resolution due to the relief. Here, we present the synthetic tests methodology: generation of synthetic gravity and potential data, then estimation of the potential from these data using the least-squares collocation and assessment of the clocks contribution. We shall see how the coverage of the data points (realistic or not) can affect the results, and discuss how to quantify the trade-off between the noise level and the number of data points used.

National audience

Most of the orbits of artificial satellites around the Earth have relatively low eccentricities. The calculation of their trajectories is very well under control, either by means of numerical methods when it comes to focus on accuracy and comparing observations, or either through analytical or semi-analytical theories to optimize the speed of calculations. This second category is used, in particular, for computing many long-term trajectories that could help to ensure the security and safety of outer space activities. Furthermore, there is also satellites operating in highly elliptical orbits for which the computations of the trajectories can be significantly improved, especially in analytical theories. Due to the fact that such orbits cover a wide range of altitudes, the hierarchy of the perturbations acting on the satellite changes with the position on the orbit. At low altitude, the oblateness of the Earth (the so-called J2 effect) is the dominant perturbation while on high-altitude the lunisolar perturbation acceleration can reach or exceed the order of the J2 acceleration. Moreover, the traditional analytical theories of Celestial Mechanics are not well adapted to this particular dynamic. Analytical solutions are generally quite developed in Fourier series of the mean anomaly. These infinite series converge very slowly in the general case, but in the case of small eccentricities convergence is accelerated thanks to the d’Alembert characteristic which ensures that a Fourier coefficient of order n keeps a degree of cancellation n with respect to the eccentricity e. This argument is no longer effective for large eccentricities. As is often the case, the solutions are more degraded in the case where the Fourier coefficients are developed in power series of the eccentricity. However, there are also theories whose the developments are made with angular variables other than the mean anomaly and well suited to the considered problem: the true anomaly for the inner potential and the eccentric anomaly for the perturbation due to external body attraction (Moon or Sun, often called third body). However, in this latter case, the methods for developing theories always use another type of approximation. Indeed, the explicit time dependence of the Hamiltonian is neglected for solving the partial differential equations that give the generator of the change of variables.This PhD thesis is devoted to propose several tools to overcome these limitations. Firstly, we expand the third-body disturbing function using Fourier series in multiples of the satellite’s eccentric anomaly (instead of the mean anomaly). We then perform a normalization of the expanded Hamiltonian, which aims to eliminate all periodic terms. To this end, we apply a change of canonical variables based on time-dependent Lie transforms. The construction of the generating function of the change of variables requires solving a partial differential equations (PDE) with respect to the angular variables of the third body and the satellite. To our knowledge this PDE has no exact solution. An approximation is usually done at this step to solve it by neglecting the terms related to the third body. We show how this approximation can be avoided by providing an iterative method for solving the PDE. This amounts to carrying out a power series expansion of a small ratio of frequencies << 1. In addition, for overall consistency of the theory we had to modify the classical solution of the central potential, especially the J 2 effect. Finally we get a theory which allows to extrapolate the osculating motion (not just the mean motion) over long periods (many decades) efficiently and with very high accuracy even for highly elliptical orbits (e > 0.7).Moreover, since the traditional numerical integration methods are not very effective for highly elliptical orbits, even with adaptive variable-step size, we show the benefits of the so-called geometric integrators, especially the variational integrators. To this end, we present some high-order numerical schemes and we test their performance., La plupart des orbites de satellites artificiels autour de la Terre ont des excentricités relativement faibles. Le calcul de leurs trajectoires est très bien maitrisé, soit au moyen de méthodes numériques quand il s’agit de privilégier la précision et la comparaison à des observations, soit au moyen de théories analytiques ou semi-analytiques pour optimiser la vitesse des calculs. Cette seconde catégorie est en particulier utilisée pour le calcul de nombreuses trajectoires à long terme dans le cadre de la sécurité de l’espace. Par ailleurs, il existe également des satellites sur des orbites à forte excentricité pour lesquels le calcul des trajectoires peut encore être très nettement amélioré, en particulier en ce qui concerne les théories analytiques. Comme ces orbites couvrent un large spectre d’altitudes, la hiérarchie des perturbations agissant sur le satellite change selon sa position sur l’orbite ; à basse altitude, la perturbation dominante est due à l’aplatissement de la Terre (appelé l’effet du J2 ) tandis que pour les altitudes élevées, les perturbations lunisolaires peuvent être aussi importantes et même supérieures à la perturbation de J2 . Par ailleurs, les théories analytiques classiques ne sont pas adaptées pour étudier la dynamique de ces orbites particulières. Les solutions analytiques sont assez généralement développées en série de Fourier de l’anomalie moyenne. Ces séries infinies convergent assez lentement dans la cas général, mais dans le cas des faibles excentricités la convergence est accélérée grâce à la caractéristique de d’Alembert qui assure qu’un coefficient de Fourier d’ordre n possède un degré d’annulation n par rapport à l’excentricité e. Cet argument perd toute son efficacité pour les fortes excentricités. Les solutions sont encore plus dégradées dans le cas où les coefficients de Fourier sont développés en séries entières de l’excentricité, comme c’est souvent le cas. Néanmoins, il existe aussi des théories dont les développements sont réalisés avec des variables angulaires, autres que l’anomalie moyenne, plus adaptées au problème étudié : l’anomalie vraie pour le potentiel intérieur et l’anomalie excentrique pour la perturbation d’un corps extérieur (Lune ou Soleil, souvent appelé troisième corps). Cependant, dans ce dernier cas, les méthodes de développement des théories utilisent toujours un autre type d’approximation : la dépendance temporelle explicite de l’Hamiltonien est négligée dans la résolution de l’équation aux dérivées partielles qui engendre le générateur du changement de variables.Cette thèse est consacrée au développement d’outils permettant de surmonter ces limitations. Dans un premier temps, nous développons la fonction perturbatrice de troisième corps en utilisant des séries de Fourier en multiples de l’anomalie excentrique du satellite (à la place de l’anomalie moyenne). Nous procédons ensuite à une normalisation du Hamiltonien ainsi développé, dans le but d’éliminer tous les termes périodiques. Pour y parvenir, nous appliquons un changement de variables canoniques construit à l’aide des transformées de Lie dépendantes du temps. La construction de la fonction génératrice du changement de variables nécessite la résolution d’une équation aux dérivées partielles (EDP) par rapport aux variables angulaires du troisième corps et du satellite. A notre connaissance cette EDP n’a pas de solution exacte. Habituellement, à cette étape, on fait une approximation qui consiste à négliger les termes de l’EDP liés au troisième corps afin de la résoudre. Nous montrons comment cette approximation peut être évitée, en proposant une méthode de résolution de l’EDP itérative, qui revient à effectuer un développement en série de puissances d’un rapport de fréquences petit devant 1. De plus, pour une cohérence globale de la théorie nous avons dû reconsidérer la solution classique de potentiel central et en particulier du J2 . Finalement nous obtenons une théorie qui permet d’extrapoler le mouvement osculateur (et pas seulement le mouvement moyen) sur de longues durées (des dizaines d’années) de façon efficace et avec une excellente précision y-compris pour des orbites très excentriques (e > 0.7).Par ailleurs, les méthodes d’intégration numérique classiques (même à pas variable) étant peu efficaces pour ces problèmes à très forte excentricité, nous montrons les avantages que procurent les intégrateurs dits géométriques, et particulièrement les intégrateurs variationnels. A cette fin, nous présenterons des schémas numériques d’ordre élevé dont nous testons les performances.

National audience

International audience