A comparison of models in inertial surveying

The output of an inertial measuring unit contains the combined effect of vehicle acceleration, accelerations due to the reference system, gravitational attraction, calibration and alignment errors, and measuring noise. The separation of individual effects from the total output is the task of modelling. Differences in the existing models are mainly due to the treatment of system errors and of the anomalous gravity field. The paper discusses the state space formulation as a general approach to represent errors of a dynamic system and describes its application to the modelling of an inertial system. Complications due to the effects of the anomalous gravity field are discussed and possible solutions are indicated. Finally, the models in the existing inertial survey systems are compared and implementation requirements for survey applications are outlined.

---

Sommaire Les mesures effectuées par l’unité de mesure inertielle continnent les effets de l’accélération du véhicule et de celle due au système de référence, de l’attraction gravitationnelle, des erreurs de calibrage et d’alignement, et finalement du bruit de fond. La séparation de chacun de ces effets sur les mesures est le sujet de la modélisation. Les différences dans les modèles existant sont principalement dues au traitement du système d’erreurs et des anomalies du champ de gravité. La communication discute de la formulation du vector d’erreurs comme approche générale pour la représentation d’erreurs dans un système dynamique et décrit son application dans la modélisation d’un système inertiel. Les complications dues aux effets des anomalies du champ de gravité sont discutées et des solutions possibles sont indiquées. Enfin, les modèles dans les présents systèmes inertiels d’arpentage sont comparés et les exigences pour les applications en arpentage sont soulignées.

---

---

Presented at the 2nd International Symposium on Inertial Technology for Surveying and Geodesy, Banff, Canada, June 1–5, 1981.     

DORIS and the Determination of the Earth’s Polar Motion

DORIS (Détermination d’Orbite et Radiopositionnement Intégrés par Satellite) is a system used for precise orbit determination (POD) and ground-station positioning. It has been implemented on-board various satellites: the SPOT (Système pour l’Observation de la Terre) remote sensing satellites SPOT-2, SPOT-3, SPOT-4, SPOT-5, TOPEX/Poseidon and more recently on its successors Jason-1 and ENVISAT. DORIS is also a terrestrial positioning system that has found many applications in geophysics and geodesy; in particular, it contributes to the realization of the International Terrestrial Reference Frame, ITRF2000 and the forthcoming ITRF2005. Although not its primary objective, DORIS can bring information on Earth orientation monitoring, mainly polar motion and length of day (LOD) variations that complement other astrogeodetic techniques. In this paper, we have analyzed various recent polar motion solutions derived from independent analysis centers using different software packages and applying various analysis strategies. Comparisons of these solutions to the International Earth Rotation and Reference Systems Service (IERS) C04 solution are performed. Depending on the solutions, the accuracy of DORIS polar components are in the range of 0.5–1 mas corresponding to a few centimeters on the Earth’s surface. This is approximately ten times larger than results derived from GPS, which are typically 0.06 mas in both components. This does not allow DORIS results to be taken into account in the IERS–EOP combinations. A gain in the precision could come from technical improvements to the DORIS system, in addition to improvement of the orbit, tropospheric, ionospheric and Earth gravity field modeling.     

Le projet gradio et la determination a haute resolution du geopotentiel

Résumé Une des techniques de détermination fine et globale du champ de gravitation terrestre U est la gradiométrie spatiale, dans laquelle on mesure à bord d'un satellite sur orbite basse certaines combinations linéaires des composantes du tenseur ∂² U/∂x_i ∂x_j dans des axes {x _i } liés au satellite. Un tel projet, appelé GRADIO, est actuellement à l'étude en France et pourrait aboutir à partir de 1990. Après avoir rappelé les objectifs scientifiques d'une telle mission, nous en donnons les spécifications—étayées par une série d'études analytiques; nous définissons ensuite le satellite porteur et ses caractères techniques, en insistant sur les points délicats de la faisabilité (facteurs d'échelle des micro-accéléromètres constituant l'appareil, connaissance de l'attitude...) et en présentant des idées de solution en cours d'approfondissement.

---

Summary Satellite gradiometry arises as one of the methods for improving our knowledge of the global Earth gravity field at high resolution: by means of micro-accelerometers on board a low orbiting spacecraft, linear combination of the gravity tensor components ∂² U/∂x_i ∂x_j are measured in a satellite-fixed reference frame {x _i }. Based on this technique, a project named GRADIO is presently under study in France and could fly in 1990 at the earliest. After the scientific objectives of that experiment have been reviewed, the measurement specifications are given as coming from various analytical studies. The platform and its characteristics are then defined: the critical realization problems (scale factors of the micro-accelerometers, spacecraft attitude control and restitution) are pointed out together with some ideas for their solution which are under analysis and require further study.

     

Regional adjustment of inertial gravity disturbance vector measurements by optimal two-dimensional smoothing

A two-dimensional signal processing algorithm is developed to obtain smoothed estimates of the gravity disturbance vector from vector measurements obtained by an inertial surveying system. The method differs from a conventional least squares regional adjustment of such measurements in that it accommodates a signal model in the smoothing process. Using principles from the physical theory of geodesy, it is shown that for a local region on the surface of the earth, an appropriate signal model is obtained by applying the two-dimensional Laplacian operator to a function representing the surface disturbance potential and equating the result to spatial white noise. The model of the vector measurement is the three orthogonal spatial derivatives of a three dimensional disturbance potential evaluated at the surface contaminated by additive white noise. The problem of simultaneous smoothing of all the gravity disturbance measurements from all survey traverses in the region is solved by representing the surface disturbance potential by a two-dimensional Karhunen-Loeve expansion that makes no specific reference to either the geometry or the ordering of the parameter space, thereby making no assumptions of causality, stationarity or isotropy. The problem of estimating the gravity anomaly and the two vertical deflection components reduces to estimating the Karhunen-Loeve coefficients which are uncorrelated and rapidly converging. Simulation results as well as smoothing of actual gravity disturbance vector measurements obtained by the U.S. Army Engineer Topographic Laboratories (USAETL) with the Rapid Geodetic Survey System (RGSS) at the White Sands Missile Range (WSMR) are presented in the paper. An analysis of these results shows that the optimal two-dimensional smoother obtains a performance benefit relative to conventional regional least squares by a factor of 2 and a benefit relative to single-traverse smoothed results by a factor of 4.

---

Sommaire Un algorithme de traitement du signal en deux dimensions est développé pour obtenir une estimation lissée du vecteur de la perturbation de la pesanteur à partir des mesures de vecteur obtenues avec un système d’arpentage inertiel. La méthode diffère d’un compensation régionale conventionnelle par moindres carrés de telles mesures, par le fait qu’elle contient un modèle du signal dans le processus de compensation. En s’appuyant sur les principes de la géodésie physique, il est montré que pour une région locale de la surface de la terre, un modèle approprié du signal est obtenu en appliquant l’opérateur à deux dimensions de Laplace à une fonction représentant le champ perturbateur à la surface de la terre et égalisant le résultat à un bruit blanc spatial. Le modèle du vecteur de mesures est défini par les trois dérivées spatiales de la fonction tridimensionnelle du potentiel perturbateur évaluées à la surface et contaminées par un bruit blanc. Le problème du lissage de toutes les mesures de gravité perturbatrice obtenues à partir des polygonales effectuées est résolu en représentant le potentiel perturbateur à la surface à l’aide d’un développement Karhunen-Loeve à deux dimensions qui ne fait aucunement référence à la géométrie ou à l’ordre des paramètres; ceci prévient toute dépendance spatiale des points adjacents. Le problème de l’estimation de l’anomalie de la gravité et des deux composantes de la déviation de la verticale se réduit à celle des coefficients Karhunen-Loeve qui sont non-corrélés et convergent rapidement. Les résultats de simulation aussi bien que le lissage des données du vecteur de perturbation de la pesanteur foumi par l’U.S. Army Engineer Topographics Labs (USAETL) sont présentés. L’analyse de ces résultats montre que le lissage optimal à deux dimensions améliore les résultats par un facteur 2 comparés aux résultats d’une compensation régionale par moindres carrés, et par un facteur 4 comparés aux résultats lissés d’une simple traverse.

---

---

Presented at the Second International Symposium on Inertial Technology for Surveying and Geodesy, Banff, Canada, June 1–5, 1981.     

The IGS-combined station coordinates, earth rotation parameters and apparent geocenter

The International GNSS Service (IGS) routinely generates a number of weekly, daily and sub-daily products. Station coordinates and velocities, earth rotation parameters (ERPs) and apparent geocenter are among these products generated weekly by the IGS Reference Frame Coordinator. They have been determined since 1999 by combining independent estimates from at least seven IGS Analysis Centers (ACs). Two Global Network Associate Analysis Centers (GNAACs) also provide independent combinations using the same AC weekly solutions and they are currently used to quality control the IGS combination. The combined solutions are aligned to an IGS realization (IGS05) of the ITRF2005 using a carefully selected set of the IGS Reference Frame (RF) stations (nominally 132). During the combination process, the contributing solutions are compared and outliers are removed to ensure a high level of consistency of the estimated parameters. The ACs and the weekly combined solution are consistent at the 1–2 and 3–4 mm levels for the horizontal and vertical components. Similarly, the excess Length of Day (LOD), the pole positions and pole rates are consistent at the 10μs, 0.03–0.05 mas and 0.10–0.20 mas/day levels, respectively. The consistency of the apparent geocenter estimate is about 5 mm in the X and Y components and 10 mm in the Z component. Comparison of the IGS-combined ERP estimates with the IERS Bulletin A suggests a small bias of the order of ?0.04 mas and + 0.05 mas (both ±0.05 mas) in the x and y components.     

Superconducting tensor gravity gradiometer

A sensitive gravity gradiometer can provide much needed gravity data of the earth and improve the accuracy of inertial navigation. A complete determination of all five independent components of the gravity gradient tensor is possible in principle by rotating a single in-line component gradiometer. In order to avoid dynamically induced noise sources arising from rotation, a hard-mounted assembly of component gradiometers may nevertheless be advantageous in an ultra-sensitive device. Superconductivity and other properties of materials at low temperatures can be used to obtain a sensitive, low-drift, gravity gradiometer. By differencing the outputs of accelerometer pairs using superconducting circuits, it is possible to construct a non-rotating tensor gravity gradiometer. Additional superconducting circuits can be provided to determine the linear and angular acceleration vectors. A three-axis in-line gravity gradiometer (a “vector” gradiometer) is being developed for satellite geodesy. A two-dimensional spring concept for a three-axis cross component gradiometer is discussed. The superconducting tensor gravity gradiometer constitutes a complete package of inertial navigation instruments with angular and linear acceleration readouts as well as gravity signals. Accuracy of inertial navigation could be improved by use of such a gradiometeraided inertial navigator.

---

Sommaire Un gradiomètre gravimétrique sensible peut foumir une grande partie des données gravimétriques terrestres et augmenter la précision de la navigation inertielle. Une détermination complète des cinq composantes indépendantes du tenseur de gradient de gravité est possible en principe par la rotation d’un seul axe d’un gradiomètre. Dans le but d’éviter les sources de bruit provenant de la rotation, un assemblage complexe de gradiomètres peut être néanmoins avantageux dans un appareil ultra-sensible. La superconductivité et d’autres propriétés des matériaux à basse température peuvent être utilisées pour obtenir un gradiomètre gravimétrique sensible et à faible dérive. Il est possible, en différentiant les sorties de paires d’accéléromètres utilisant des circuits superconducteurs, de construire un gradimètre non-rotatif qui mesure le tenseur de gravité. Des circuits superconducteurs additionnels peuvent être ajoutés pour déterminer les vecteurs d’accélération linéaires et angulaires. Un gradiomètre gravimétrique à trois exes est présentement développé pour la géodésie spatiale. Le concept d’un ressort à deux dimensions comme composante d’un gradiomètre à trois axes est discuté. Le gradiomètre de mesure du tenseur de gravité à l’aide de circuits superconducteurs constitue l’ensemble complet pour un instrument de navigation inertielle avec lecture des accélérations angulaires et linéaires et des données de gravité. La précision de la navigation inertielle pourrait être améliorée par l’emploi d’un tel gradiomètre avec les systèmes inertiels.

---

---

Presented at the Second International Symposium on Inertial Technology of Surveying and Geodesy, Banff, Canada, June 1–5, 1981.     

Long geodesics on the ellipsoid

Summary This article examines the practical application of formulae for computing long lines on the ellipsoid. The main aim is to eliminate the successive approximation generally required. For the inverse problem, this is achieved by the method ofE. M. Sodano, Army Map Service, U.S.A. An adaptation of a method produced byG. T. McCaw is used for the direct problem. Results are given of five practical examples, including two which extend halfway round the world. Construction of further special tables is recommended to simplify the computations required by a problem which has an ever increasing application.

---

Zusammenfassung Der Verfasser—ein bekannter Fachmann auf diesem Gebiet—untersucht kritisch die von verschiedenen Geod?ten angewandten oder vorgeschlagenen Verfahren zur Berechnung sehr langer geod?tischer Linien.

---

Resumen El autor, especialista bien conocido en la cuestión, examina y discute los differentes métodos utilizados o preconizados por diferentes geodestas para el cálculo de lineas geodésicas de gran longitud en la superficie terrestre.

---

Résumé L'auteur, spécialiste bien connu de la question, examine et discute les différentes méthodes utilisées ou préconisées par différents Géodésiens, pour le calcul des lignes géodésiques de grande longueur à la surface de la terre.

---

Riassunto L'Autore, specialista ben noto dell'argomento, esamina e discute i diversi metodi utilizzati o proposti dai differenti Geodeti per il calcolo delle geodetiche lunghe sulla superficie terrestre.

     

Contribution au probleme de la compensation des reseaux de triangulation

Résumé Dans cette étude on expose une adaptation de la méthode des directions —pour la compensation d’un réseau de triangulation—aux possibilités de l’ordinateur électronique. On présente les formules par lesquelles on calcule les coefficients des inconnues dans les équations de condition (équations aux angles et équations aux c?tés) et on dresse leur matrice. Ensuite on traite de la formation de la matrice des coefficients des équations normales et son inversion, qui fournit les quantités corrélatives de Lagrange. Enfin, après avoir déterminé les corrections à apporter aux directions observées, on calcule l’erreur moyenne quadratique d’une observation isolée et l’erreur moyenne de chacune des directions compensées.

---

Summary In this study, an adaptation to the computers’ possibilities, regarding the method of directions, is disclosed in order to realise the adjustment of triangulation nets. The computing formulae of the coefficients of the unknown quantities in the condition equations and the creation of their matrix are given herebelow. The treatment leading towards the construction of the matrix of normal equations and its inversion, that furnishes Lagrange’s quantities, follows. After the computation of the corrections applied on the observed directions, the mean error of a single observation and the mean error of every adjusted direction are determined. The study is executed in such a way that it can be included in only one programme for automatic computation.

     

The deflection of the vertical in the Western and Central Mediterranean area

The deflection of the vertical at 23 stations was determined by gravimetrical methods and compared to the difference of astronomic observations and geodetic values of the European adjustment. The results show a fair agreement between the two independent sets of values, indicating that doubts about the relative reliability of submarine gravity observations are not justified. A systematic discrepancy in the meridian components of the deflection may reflect on theHayford ellipsoid as used in the European adjustment.

---

Zusammenfassung Auf 23 Stationen wurde die Lotabweichung gravimetrisch bestimmt und den astro-geod?tischen Werten der europ?ischen Ausgleichung gegenübergestellt. Die beiden Ergebnisse stimmen sehr gut überein. Damit werden die Bedenken hinf?llig, die sich gegen den Wert von Messungen in Unterseebooten richteten. Systematische Abweichungen in den meridionalen Komponenten beruhen m?glicherweise auf der art der Anwendung des Hayfordschen Ellipsoids in der Ausgleichung des Zentraleurop?ischen Netzes.

---

Resumen La desviación de la vertical en 23 estaciones ha sido determinada por el método gravimétrico y compara con los resultados astrogeodéticos de la compensación europea. Existe muy buen acuerdo entre los dos conjuntos de datos, que demuestra que las dudas emitidas en cuanto al valor de las observaciones en submarino son injustificadas. Una discordancia sistemática en las componentes meridianas podria ser debida al empleo del elipsoide deHayford en la compensación de la red europea.

---

Résumé La déviation de la verticale en 23 stations a été déterminée par la méthode gravimétrique et comparée aux résultats astrogéodésiques de la compensation européenne. Il y a très bon accord entre les 2 ensembles de données, ce qui montre que les doutes émis quant à la valeur des observations en sous-marins sont injustifiés. Une discordance systématique dans les composantes méridiennes pourrait être due à la fa?on dont on a employé l'ellipso?de deHayford dans la compensation du réseau central européen.

---

Sommario E'stata determinata col metodo gravimetrico la deviazione della verticale in 23 stazioni e comparata ai risultati astro-geodetici della compensazione europea. Vi é molto buon accordo tra i due insiemi di dati, ciò che dimostra che i dubbi emessi circa il valore delle osservazioni in sottomarino sono ingiustificati. Una discordanza sistematica tra le componenti maridiane potrebbe essere dovuta all'impiego nella compensazione della rete central europea dell'ellisoide diHayford.

     

Azimuth determination by equatorial stars

Summary Difficulties in the determination of precise azimuth from observations on Polaris at high latitudes suggested an investigation of azimuth determination from other stars. A proposed method for determining azimuth by observing the time of crossing of stars in the vertical plane through the reference object is discussed theoretically, and results of a field test of the method are given.

---

Zusammenfassung Genaue Azimutbestimmungen mit Hilfe des Polarsterns sind in h?heren Breiten schwierig. Der Verfasser untersucht daher die Azimutbestimmung mit Hilfe anderer Sterne. Der Aufsatz er?rtert die Theorie eines Messverfahrens, bei dem die Durchgangszeiten von Sternen durch die Referenz-Lotebene beobachtet werden. Ergebnisse von Versuchsmessungen werden mitgeteilt.

---

Resumen La dificultad de efectuar determinaciones precisas de azimut con la estrella polar, han conducido al autor a estudiar esta determinación con otras estrellas. El artículo discute la teoría de un método de determinación por observación del instante de pasos de estrellas por el plano vertical de la referencia y presenta resultados experimentales.

---

Résumé La difficulté des déterminations précises d'azimuth sur l'étoile polaire, aux hautes latitudes, ont amené l'auteur à étudier cette détermination sur d'autres étoiles. L'article discute la théorie d'une méthode de détermination par observation de l'instant de passage des étoiles dans le plan vertical de la référence et présente des résultats expérimentaux.

---

Riassunto La difficoltà di effettuare determinazioni precise d'azimut con la Polare, ad alte latitudini, ha indotto l'Autore a studiare tale determinazione con altre stelle. L'articolo tratta la teoria di un metodo di determinazione con osservazione dell'istante di passaggio delle stelle nel piano verticale di riferimento e presenta alcuni risultati sperimentali.

     

Evaluation et representation des erreurs sur les deformations d'un reseau geodesique: Utilisation de la methode de Monte Carlo

Résumé Le problème de la comparison des données obtenues à différentes époques sur un réseau geódésique pour en déduire les déformations subies par ce réseau au cours du temps est discuté, et conduit à proposer un mode de représentation simultané des déformations et des erreurs, utilisant les tenseurs de déformations classiques avec une méthode de Monte Carlo permettant de simuler un grand nombre de séries de mesures.

---

The repetition of geodetic measurements is employed in many cases to monitor ground deformations, and is widely used for civil engineering, geodynamics and slope instabilities studies. Generally the networks used are small and the main problem to solve is how to detect stable zones, i.e. where the signal/noise ratio of the deformation measurements are below unity. This is commonly done using statistical tests. We present here a new tool, the simultaneous display of the deformations and their related errors, using a classical Monte Carlo methodology when computing strain tensors. Such a representation provides to non-specialists in a quit efficient way the significance level of the deformation patterns. Some examples show that even when statistical test show a poor confidence level, this methodology may give anyway in some situations the pattern of the deformation (but—obviously—not its amplitude).

     

The generalization of the theorem of Jacobi

The present study deals with the inner mechanism of the method of least squares, and thus constitutes an extension of previous work, [2] and [3]. It covers an investigation as to whether there are corresponding relationships between the results of partial adjustments on the one hand, and a total adjustment of the same data on the other, as expressed by the theorem ofJacobi. It is shown that such relationships do, in fact, exist: the results of a total adjustment are identical (or proportional) to the weighted means based on the corresponding values found by the partial adjustments, with the determinants of the latter’s systems of normal equations as weights. This applies to the quantities determined by the adjustment and their accuracy, as well as to their internal correlation and the sum of squares of the residuals.

---

Zusammenfassung Der vorliegende Aufsatz besch?ftigt sich mit dem Mechanismus der Methode der kleinsten Quadrate and stellt somit eine Weiterführung der früheren Arbeiten [2] und [3] dar. Diesmal Wird untersucht, ob entsprechende Zusammenh?nge, welche durch den Satz vonJacobi zum Ausdruck kommen, zwischen den Ergebnissen von Teilausgleichungen einerseits und der Gesamtausgleichung desselben Beobachtungsmateriales andererseits, bestehen. Es wird nachgewiesen, dass dies tats?chlich der Fall ist: Die Ergebnisse der Gesamtausgleichung sind identisch mit (oder proportional zu) den gewogenen Mitteln der korrespondierenden Werte der Teilausgleichungen, wobei als Gewichte die Determinanten der Normalgleichunssysteme der letzteren auftreten. Dies trifft sowohl für die Werte der durch die Ausgleichung zu bestimmenden Gr?ssen, für die Genauigkeit und der internen Korrelation derselben, als auch für die Fehlerquadratsumme zu.

---

Résumé Le présent article traite du mécanisme interne de la méthode des moindres carrés et constitue donc une suite aux travaux précédents [2] et [3]. Les recherches portèrent cette fois-ci sur l’existence de relations correspondantes entre les résultats des compensations partielles (compensations ne comprenant qu’une partie des observations), d’une part, et une compensation totale de la même série d’observations, d’autre part, comme elle est exprimée par le théorème deJacobi. Il est démontré ici que de telles relations existent effectivement. Le résultat d’une compensation totale est identique (ou proportionnel) aux moyennes arithmétiques générales fondées sur des valeurs correspondantes résultant des compensations partielles et des déterminants des équations normales de ces dernières comme poids. Ceci est valable aussi bien pour la grandeur des quantités qui sont déterminées lors de la compensation que pour la précision et la corrélation interne de ces mêmes quantités, ainsi que pour la somme des corrections quadratiques.

     

A non-iterative and non-singular perturbation solution for transforming Cartesian to geodetic coordinates

The Cartesian-to-Geodetic coordinate transformation is re-cast as a minimization algorithm for the height of the Satellite above the reference Earth surface. Optimal necessary conditions are obtained that fix the satellite ground track vector components in the Earth surface. The introduction of an artificial perturbation variable yields a rapidly converging second order power series solution. The initial starting guess for the solution provides 3–4 digits of precision. Convergence of the perturbation series expansion is accelerated by replacing the series solution with a Padé approximation. For satellites with heights < 30,000 km the second-order expansions yields ~mm satellite geodetic height errors and ~10⁻¹² rad errors for the geodetic latitude. No quartic or cubic solutions are required: the algorithm is both non-iterative and non-singular. Only two square root and two arc-tan calculations are required for the entire transformation. The proposed algorithm has been measured to be ~41% faster than the celebrated Bowring method. Several numerical examples are provided to demonstrate the effectiveness of the new algorithm.     

Letters from Carlos Ibáñez e Ibáñez de Ibero to Aimé Laussedat: new sources for the history of nineteenth century geodesy

Recently found letters written by Carlos Ibáñez e Ibáñez de Ibero (1825–1891), first president of the International Geodetic Association, to geodesist–photogrammetrist Aimé Laussedat (1819–1907) are discussed. The contents of these epistolary conversations reflect their personal preoccupation with solving the geodetic problems of the time (1857–1889) and their desire to develop and propagate geodetic knowledge across all international boundaries.     

Neue Analyse der Chandler-Bewegung

Summary TheChandler movement is essentially an oscillation having an assumed period of the order of time of a double revolution of the moon’s nodes. The variation in its period is only apparent and relies upon a sudden change of phase which occurred at the beginning of 1928. If attention is paid to the residuals, a secondary oscillation with an 11-year period can be detected. This explains theWitting perturbations, and is possibly connected with the sun-spot cycle.

---

Resumen El movimientoChandler es en su esencia un balanceo que tiene una duración admitida del órden de la doble revolución de los nodos de la Luna. La variabilidad de su periodo no es más que simulada, debido a un salto de fase situada en el comienzo del a?o 1928. Utilizando los residuos, se puede demostrar un balanceo suplementario de 11 a?os, que explica lasperturbacionésWitting y que está posiblemente relacionado con el periodo de las manchas solares.

---

Résumé La période deChandler correspond, dans son essence, à un balancement ayant une durée présumée de l’ordre de la double révolution des nœuds de la Lune. La variabilité de sa période n’est que dissimulée et ceci grace à un saut de phase situé au commencement de l’année 1928. En utilisant les résidus, on peut démontrer un balancement supplémentaire de 11 ans, qui explique les perturbationsWitting, et qui est en liaison possible avec la période des taches solaires.

---

Sommario Il movimento diChandler è essenzialmente un’oscillazione avente presumibilmente un periodo doppio della rivoluzione dei nodi lunari. La variabilità di questo periodo è soltanto apparente, in consequenza di un salto di fase situato al principio del 1928. Utilizzando i residui, si può mostrare un’oscillazione supplementare del periodo di 11 anni, che spiega la perturbazione diWitting e che sta forse in relazione con il periodo delle macchie solari.

---

---

Communication présentée à l’Assemblée Générale de Bruxelles     

Analysis of secular polar motion and continental drift

The secular latitude variations of the five ILS stations of Mizusawa, Kitab, Carloforte, Gaithersburg and Ukiah were analyzed taking into account the recent continental drift theory. Using Le Pichon's 1968 reconstruction, the rate of rotation was computed from the astronomical data, fixing the pole of rotation by Le Pichon's determination. The most reasonable solution was obtained considering Mizusawa, Kitab and Carloforte lying on the Eurasia plate, the two American stations as one on the American plate (Gaithersburg) and the other on the North—East Pacific plate (Ukiah). The resulting relative rate between the Euro-American plates is found to be 0".0028/year and between the American—Pacific plates 0".0032/ years, or about 1°,3/10⁶ years and in excellent agreement with the plate tectonic theory. Luxembourg Meeting of the “Journées Luxembourgeoises de Géodynamique”, 1972.