Les mouvements de l'ecorce terrestre et leur observation geodesique

Conclusion Les exemples cités dans cet article et étudiés en Tchécoslovaquie montrent que sur la terre entière, se produisent des mouvements de l'écorce terrestre, horizontaux et verticaux, positifs et négatifs. En même temps on peut constater que ce sont les régions de tectonique jeune qui subissent ces mouvements et qu'ils s'y manifestent beaucoup plus intensivement (système slovaco—carpathien) que sur les masses de continent plus stables de l'écorce terrestre (Bohême et Moravie—Massif de Bohême) et que dans les régions séismiquement actives (bassin de Komarno en Slovaquie du Sud) leur direction peut changer. On peut supposer que ces mouvements provoquent certaines tensions dans les roches de l'écorce terrestre. Si une réserve suffisante d'énergie potentielle s'accumule en forme de tensions élastiques, il suffit d'une impulsion insignifiante (par ex. des dérangements aux lieux, fatigués tectoniquement, une altération du milieu minier par l'exploitation, etc..) pour transformer l'énergie potentielle accumulée en énergie cinétique. Ensuite se produisent des pressions des roches, des secousses et des dérangements différents de la structure des roches et par suite des mouvements secondaires. A. PELNAR conclut également que les secousses dans les mines de Pribram en Bohême peuvent être causées par des mouvements verticaux petits, mais bien impétueux des massifs, montagneux, moyennant quoi une énergie considérable se libère par un, écrasement des établissements miniers. Comme les mouvements mentionnés sont fonction du temps, les changements de tension et les déformations dans les roches et leurs conséquences se présentent comme une fonction du temps. L'importance des études de l'écorce terrestre est immense non seulement au point de vue scientifique (géologie, géodésie, géophysique), mais aussi au point de vue pratique (géologie d'ingénieur, industrie minière). L'étude détaillée des mouvements de l'écorce terrestre, de leurs causes, de l'énergétique, de l'action, de la connection et des conséquences de leurs marques permettra de résoudre un nombre de problèmes, parmi les autres même le problème des pressions et des couches montagneuses et aidera à résoudre maintes discussions à cet égard (voir par ex. la discussion entre le savant tchèque A. PELNAR et le spécialiste G. SPACKELLER sur la question de savoir si les secousses à Ostrava sont causées par les dernières manifestations du plissement varisque ou si elles sont provoquées par l'exploitation des mines ou toute autre cause inconnue. Les mouvements du groupe B joueront aussi un r?le important dans la solution de ces problèmes.     

Valeur de la Pesanteur pour le Point Fondamental de Varsovie

Résumé Le système provisoire des mesures gravimétriques en Pologne a été établi pour le point fondamental de Varsovie sur la valeur: g=981,2412cm·sec⁻². De l'analyse des résultats obtenus par le système de rattachements modernes de Varsovie à Potsdam dont les mesures ont été effectuées par A. Kwiatkowski, P. Lejay, Reicheneder, Weiken, l'auteur déduit dans le système de Potsdam comme valeur de pesanteur la plus probable à Varsovie g=981,2400±0,0002 cm·sec⁻². Cette valeur a été vérifiée au moyen de l'analyse des rattachements indirects de Varsovie à Potsdam par les points pendulaires de Kielce, Cracovie, Racibórz (tableau 2), de Ciechanów, Goldap, Susz et Szczytno (tableau 3). Les données prises en considération proviennent de mesures effectuées avec des gravimètres de N?rgaard. La concordance des résultats particuliers a été vérifiée dans les limites d'erreurs moyennes d'observations. Dans la dernière partie de cette étude ont été consignés les résultats des dernières mesures absolues de la pesanteur à Washington et à Teddington qui ont permis de constater entre le système absolu et celui de Potsdam une différence de 15 à 20 mgals. Comité National Polonais de Géodésie     

Probleme de l'implantation d'une grille sur une shphere

Résumé Dans une première partie, traitée de fa?on élémentaire, l'auteur, après avoir rappelé diverses solutions particulières, propose deux méthodes pour résoudre le problème de l'optimisation d'un réseau de points de grille sur une sphère, admettant toutes deux une symétrie équatoriale, des points équidistants sur des parallèles régulièrement espacés en latitude, une ligne de points équatoriaux et deux points polaires. Une critique permet de dégager la meilleure de ces deux solutions; cette dernière solution convient, quelque soit le nombre N de points, si N>5, mais conviendra d'autant mieux que N sera plus grand. Une application numérique est donnée pour un réseau de 1.000 points, avec la localisation de ces points sur une carte de l'hémisphère terrestre Nord. Dans une seconde partie, traitée de fa?on plus théorique, l'auteur précise la notion de distance interpoints, montre les limites de validité des solutions précédentes pour finalement proposer une troisième solution, la solution dite “en spirale”, qui lui para?t définitivement la meilleure.     

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.     

Campagne Diademe Resultats obtenus en geodesie semi-dynamique a l'aide de mesures Doppler

Conclusions Le traitement des mesures Doppler sur des périodes voisines de 1 jour permet, par correction différentielle utilisant un bon modèle théorique, de retrouver les coordonnées,X, Y des stations du triangle méditerranéen avec des dispersions respectives de50 mètres et15 mètres. Bien que la coordonnéeZ soit moins bien déterminée, dispersion pouvant atteindre200 mètres, on peut retrouver les dimensions du triangle. L'analyse spectrale de la matrice de covariance fait appara?tre comme la plus grande composante de l'ellipso?de d'erreur, la cote du centre de gravité au-dessus du plan équatorial. Ceci peut s'expliquer par un manque de diversité des passages (période courte), un réseau d'observation trop restreint et dans le même hémisphère (mauvaise détermination dui), des imperfections du modèle de perturbations (Harmoniques de résonances, frottement). L'élimination des paramètres internes et la résolution, sur les seuls déplacements significatifs, diminue la dispersion à une dizaine de mètres surX, Y etZ. Ces résultats ne constituent qu'une première étape d'un traitement systématique en cours, qui sera exposé au COSPAR 1969.     

Revision of the first order triangulation in the region disturbed by the Nankaido Earthquake of 1946

Zusammenfassung Nach dem Erdbeden von Nanka?do im Jahre 1946 hat man das Dreiecksnetz I. Ordnung in einem Gebiet von etwa 300 km Umkreis wiederbeobachtet, und zwar auf 47 Punkten insgeasamt 68 Winkel. Die letzte Triangulation fand in den Jahren 1885–96 statt.—Nach der Ausgleichung des Netzes wurden Punkverschiebungen bis zu 1,69 m festgestellt.

---

Resumen Después del temblor de tierra de 1946, en Nanka?do, se ha vuelto a observar la red de triangulación de 1_H orden en un radio de 300 kilómetros aproximadamente. Los puntos de nuevo observados son 47 y los triángulos 68. La triangulaciín precedente se remontaba a los a?os 1885–96. Después de compensar la red, se han deducido los desplazamientos experimentados por los puntos sobre el terreno, que llegaron a ser de 1,69 metros.

---

Résumé Après le tremblement de tèrre de 1946 à Nankaido, on a réobservé le réseau de triangulation de premier ordre dans un rayon de 300 km environ. Les points réobservés sont au nombre de 47, et les triangles en nombre de 68. La triangtlation précédente remontait aux années 1885–96. La compensation du réseau effectuée, on en a déduit les déplacements subis par les points sur le terrain, qui importent jusqu’à 1,69 m.

---

Sommario Dopo il terremoto di Nankaido del 1946, è stata riosservata la rete di prim’ordine entro un raggio di circa 300 km. La precedente triangolazione risaliva agli anni 1885–96. Effettuata la compensazione, il confronto fra l’anticae la nuova triangolazione ha rivelato spostamenti nei vertici ammontanti fino a m 1,69.

---

---

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

La representation du potentiel terrestre par masses ponctuelles

Résumé Un modèle de représentation du potential terrestre par 126 masses ponctuelles de profondeurs situées entre 1000 et 1500 km a été construit à partir de la Standard Earth II du S.A.O. et utilisé avec succès pour la représentation du géo?de, des anomalies de gravité, ainsi qu’en calcul d’orbites par correction différentielle utilisant des observations réelles de satellites artificiels. Les propriétés de décroissance rapide de ces fonctions sont mises en évidence, et leur utilisation envisagée à l’analyse scientifique de mesures altimétriques.     

Adjustment of European Triangulation

Zusammenfassung Bei der Ausgleichung des Südwest-Blockes und des Nordblockes der europ?ischen Triangulation war es nicht m?glich, die Bowie-Methode anzuwenden, weil die Grundlinien und die Laplaceschen Punkte sich im allgemeinen nicht in den Knotennetzen zwischen Meridianund Parallelkreisketten befinden. Mit Rücksicht hierauf gibt es in dem Schema nur einige wenige Netzteile, die zun?chst für sich ausgeglichen werden konnten; in der Hauptsache mu\te fast der ganze Block in einem Gu? ausgeglichen werden. Die Ausgleichung des südwesteurop?ischen Blockes erfolgte nach der Methode der bedingten Beobachtungen; die Gleichungen der verschiedenen Art wurden durch die Mathematiker des Coast and Geodetic Survey aufgestellt und sofort auf Lochkarten der International Business Machines (IB.M.) übertragen, soda? alle Berechnungen mit Hilfe der I.B.M.-Maschinen rein mechanisch durchgeführt werden konnten. Die Aufl?sung der 2348 Normalgleichungen des süwesteurop?ischen Blocks erfolgte nach der Methode von Doolittle. Die Rechnungen erfolgten auf 6 Dezimalstellen. Die Ausgleichung des Südwest-Blocks erforderte fast 18 Monate Arbeit. Die Ausgleichung des nordeurop?ischen Blocks wurde im Gegensatz hierzu nach vermittelnden Beobachtungen durchgeführt, wobei die vorhandenen Basis- und Laplacebedingungen besonders berücksichtigt werden mu?ten. Sie umfa?t auch den D?nisch-Norwegischen Anschlu?, der mit Hilfe der Hochzieltriangulation hergestellt wurde. Die bereits verliegende Ausgleichung des baltischen Ringes hat die Arbeit weitgebend erleirchtert. Der nordeurop?ische Block ergab ein System von 2475 Normalgleichungen, die in der sehr kurzen Zeit von 3 Monaten gel?st wurden, wobei an den Maschinen sechs Tage pro Woche und t?glich jeweils 16 Stunden gearbeitet wurde. Die Durchführung der Arbeit hat gezeigt, da? die Aufl?sung gro?er Normalgleichungssysteme keine besonderen Schwierigkeiten mehr bereitet, und da? es kaum vorteilhaft ist, in den schw?cheren Teilen des Netzes Ketten auszuscheiden. Dieser Fall trifft in Italien zu. Weiterhin zeigt die Arbeit 1) da? es wohl angebracht w?re, in einigen Teilen des Netzes die Beobachtungen zu wiederholen, 2) da? die astronomischen Fundamental-Meridiane verschiedener L?nder, auf die sich die L?ngen beziehen, oft nur ungenügend übereinstimmen.

---

Resumen En la compensación europea no fué posible adoptar el método deBowie, dado que las bases geodésicas y los puntos deLaplace no se encuentran en general en los cruces entre las cadenas meridianas y paralelas; por esta razón no hay en el esquema sino algunas raras figuras compensadas de antemano; todo el resto de la red está compensada en un solo bloque. La compensación del bloque de la Europa Sud-Oeste ha sido efectuada por el método de observaciones condicionadas; las ecuaciones de los diferentes tipos fueron escritas por los matemáticos del Coast and Geodetic Survey, y llevadas inmediatamente a cartas perforadas de la “International Business Machines” (I.B.M.); a partir de ese momento, todo el cálculo se prosiguió por via mecánica con la intervención de las máquinas I.B.M. El método empleado para la resolución de los sistemas normales, cuyo número de ecuaciones para ese bloque era de 2348, es el de Doolittle, empleando inicialmente 6 cifras decimales; el cálculo de ese bloque exigió aproximadamente 18 meses de trabajo. La compensación del bloque de la Europa Norte ha sido efectuado. al contrario, por el método de variación de coordenadas, teniendo tam bién en cuenta el enlace geodésico entre Noruega y Dinamarca, efectuado por el método de los cohetes. La compensación que ya existía del anillo báltico ha facilitado mucho la tarea. El método empleado en el cálculo ha sido el de las observaciones indirectas condicionadas, siendo regidas las condiciones por la concordancia de las bases existentes y de los azimutes de Laplace. El bloque de la Europa Norte ha conducido a un sistema de 2475 ecuaciones que han sido resueltas en el tiempo extremadamente corto de 3 meses, trabajando las máquinas 6 dias por semana, durante 16 horas por dia. La experiencia adquirida durante la ejecución del trabajo, ha demostrado que no hay dificultad alguna para resolver grandes sistemas de ecuaciones simultáneas, y que, por tanto, no representa gran ventaja extraer cadenas en las partes más delicadas de la red. Es el caso que se ha presentado en Italia. Sin embargo, se infiere del trabajo efectuado, que sería conveniente volver a observar la red en algunas zonas, y que con frecuencia no hay concordancia suficiente entre los meridianos astronómicos fundamentales de los diferentes paises, en los que se apoyan las determinaciones de longitud.

---

Résumé Il ne fut pas possible, dans la compensation européenne, d’adopter la méthode deBowie, étant donné que les bases géodésiques et les points deLaplace ne se trouvent en général pas aux croisements entre cha?nes méridiennes et parallèles; il n’y a dans le schéma que quelques rares figures rigides compensées à l’avance; tout le reste du réseau étant compensé en un seul bloc. La compensation du bloc de l’Europe Sud-Ouest a été effectuée par la méthode des observations conditionnées; les équations des différents types étaient écrites par les calculateurs duCoast and Geodetic Survey, et reportées immédiatement sur des cartes perforées de l’International Business Machines (I.B.M.); à partir de ce moment-là, tout le calcul se poursuivit par voie mécanique à l’aide des machines I.B.M. La méthode employée pour la résolution des systèmes normaux, dont le nombre d’équations était pour ce bloc de 2348, est celle de Doolittle, en employant initialement 6 chiffres décimaux; le calcul de ce bloc demanda à peu près 18 mois de travail. La compensation du bloc de l’Europe Nord a été effectuée au contraire par la méthode de variation de coordonnées, tenant compte aussi de la jonction géodésique entre Norvège et Danemark effectuée au moyen de la méthode des fusées. La compensation existante de l’anneau baltique a grandement facilité la tache. La méthode employée dans le calcul a été celle des observations indirectes conditionnées, les conditions étant dues à l’accord des bases existantes et des azimuths laplaciens. Le bloc de l’Europe Nord a conduit sur un système de 2475 équations qui ont été résolues dans le temps très court de 3 mois, en travaillant aux machines six jours par semaine, pendant 16 heures par jour. L’expérience acquise pendant l’exécution du travail a montré qu’il n’y a aucune difficulté à résoudre des grands systèmes d’équations simultanées, et que pourtant il n’y a pas grand avantage d’extraire des cha?nes dans les parties plus délicates du réseau. C’est le cas qui s’est présenté en Italie. Du travail effectué, il ressort encore qu’il serait convenable de réobserver le réseau en quelques zones; et qu’il n’y a pas souvent accord suffisant entre les méridiens astronomiques fondamentaux des différents pays. auxquels s’appuient les déterminations de longitude.

---

Sommario Nella compensazione europea non fu possibile adottare il metodo diBowie, poichè le basi geodetiche ed i punti diLaplace non si trovano in generale agli incroci fra catene parallele e meridiane; per questa ragione non c’è nella rete che qualche rara figura rigida compensata preventivamente; tutto il resto della rete è stato compensato di un solo getto. La compensazione del blocco dell’Europa Settentrionale è stata effettuata con il metodo delle osservazioni condizionate; le varie equazioni sono state scritte dai matematici delCoast and Geodetic Survey, e trasportate immediatamente su schede perforate della International Business Machines (I.B.M.); da questo momento il calcolo si è svolto meccanicamente mediante le macchine I.B.M. Il metodo impiegato per la risoluzione del sistema normale, il cui numero di equazioni era, per questo blocco, di 2348, è quello detto di Doolittle; si sono impiegate inizialmente 6 cifre decimali, ed il calcolo per tutto il blocco ha richiesto 18 mesi di lavoro circa. La compensazione del blocco dell’Europa Settentrionale è stata invece effettuata con il metodo per variazione di coordinate, tenendo anche conto del collegamento geodetico mediante razzi illuminanti esistente fra la Danimarca e la Norvegia. L’esistente compensazione dell’anello baltico ha grandemente facilitato il compito. Il metodo di calcolo è stato quello delle osservazioni indirette condizionate, le condizioni essendo dovute all’esistenza di basi misurate e degli azimut laplaciani. Il calcolo del blocco dell’Europa Settentrionale ha richiesto la risoluzione di un sistema di 2475 equazioni, che sono state risolte nel tempo brevissimo di 3 mesi, operando le macchine sei giorni per settimana, e durante 16 ore al giorno. L’esperienza acquistata durante l’esecuzione del lavoro, ha mostrato che non c’è alcuna grande difficoltà nel risolvere grandi sistemi di equazioni simultanee, e che pertanto non si realizza un grande vantaggio nell’estrarre catene nelle parti più delicate di una rete continua. E’il caso che si è presentato in Italia. Ancora appare dal lavoro effettuato, che sarebbe opportuno riosservare la rete in qualche parte; e che spesso non c’è sufficiente accordo fra i meridiani astronomici fondamentali dei vari Paesi, ai quali si appoggiano le determinazioni astronomiche.

---

---

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

Un nouvel aspect concernant la determination du champ gravifique exterieur de la terre

Résumé Utilisant l’équation intégrale de la Géodésie-Physique de Molodenski, les auteurs proposent une méthode de calcul pour le champ gravifique en un pointM de la surface de la Terre en fonction des mesures astronomo-géodésiques effectuées dans les régions environnantes. Observant que les méthodes actuelles négligent les zones situées dans le voisinage immédiat deM à cause de la singularité du noyau de l’équation intégrale en ce point, les auteurs tiennent compte aussi de ces zones ainsi que de la configuration approximative du terrain avoisinant.     

Etude comparative des differents resultats obtenus en geodesie spatiale sur le reseau europeen

Conclusions Les différentes méthodes de géodésie spatiale utilisées sur le réseau européen ont prouvé la possibilité d'établir un réseau primaire à grandes milles, entre1500 et2000 kilomètres, avec une précision de3 à5 mètres sur les cordes, et d'environ0″5 sur l'orientation, et peut-être mieux. Le dépouillement complet de l'expérienceRCP 133 pourra confirmer ces chiffres et permettra de vérifier en particulier le facteur d'échelle trouvé. Les méthodes géométriques semblent les plus adaptées à ce travail, mais il serait nécessaire d'atteindre la précision ultime des Lasers:30 cm et d'utiliser les photographies d'écho Laser, l'expérience préliminaire réalisée par l'ONERA étant plus qu'encourageante à ce sujet. Les méthodes semi-dynamiques ont prouvé qu'avec un bon modèle de potentiel, on pourrait retrouver, même à l'échelle d'une portion de continent, la position du centre de gravité de la terre avec une précision d'environ20 mètres. Là encore, on peut espérer mieux: amélioration du modèle, choix de la période de traitement, utilisation d'autres satellites, utilisation de mesures photographiques (pour fixer l'orientation) et Laser. Le réseau constitué par les observations sur Pageos, avec pour échelle la distance San Fernando—Haute Provence, sera déjà une première réalisation de réseaux primaires. Mais il serait probablement souhaitable, à l'occasion du lancement de GEOS-C, d'établir un réseau optimisé et avec le maximum de précision.     

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