LES PASSAGES DE MERCURE ET DE VÉNUS DEVANT LE SOLEIL
FICHES PÉDAGOGIQUES: FICHE N° 08

Afficher cette fiche au format Acrobat Reader (210Ko)

Méthodes simplifiées du calcul de la parallaxe : Exemples

I. Introduction

Dans cette fiche nous allons présenter un calcul simplifié de la parallaxe équatoriale moyenne du Soleil. Cette simplification du calcul se fait au prix dune contrainte importante pour les observations. Cette contrainte est la simultanéité des observations. Nous allons supposer que nous avons deux observations simultanées qui nous fournissent la distance entre les deux centres apparents de la planète Vénus devant le disque solaire. Nous indiquerons les approximations et simplifications que nous effectuons.

Nous donnerons également deux exemples dutilisation des formules simplifiées de la méthode de Delisle et de Halley permettant de calculer une première approximation de la parallaxe équatoriale moyenne de Soleil.

II Hypothèses

Soient deux lieux dobservations M1 et M2, suffisamment éloignés ; deux observateurs notent au même instant t la position du centre apparent de la planète Vénus devant le disque solaire, puis à laide de ces deux observations, ils déterminent la distance qui joint ces deux centres apparents de Vénus. La mesure de cette distance exprimée en rayon solaire, permet de calculer la parallaxe moyenne équatoriale du Soleil p0. Nous verrons que cette mesure est loin dêtre simple.

Figure 1. – Observation du passage de Vénus depuis deux lieux au même instant.

Soient O le centre de la Terre, C le centre du Soleil, V le centre de Vénus et V1 et V2 les centres de Vénus sur le disque solaire vus respectivement depuis les points M1 et M2. Notons D1 et D2 les angles CM1V et CM2V formés par les directions des droites joignant les deux points dobservations aux centres de Vénus et du Soleil et notons ps langle sous lequel on voit le segment M1M2 depuis le Soleil et pv langle sous lequel on voit le segment M1M2 depuis Vénus. Ces deux angles sont les parallaxes du Soleil et de Vénus vues depuis les lieux M1 et M2 (figure 1).

Si les deux points M1 et M2 sont quelconques sur la surface terrestre dans la zone de visibilité du passage, il ny a aucune raison pour que les quatre points M1, M2, V et C soient dans un même plan. Donc les droites M1C et M2V ne sont pas dans le même plan et ne se coupent pas. On ne peut donc pas appliquer les règles de la géométrie plane dans la figure1.

Et la relation suivante :

est fausse. Elle nest vraie que lorsque les quatre points sont coplanaires.

Par contre la différence des parallaxes est égale à la distance angulaire entre les deux centres apparents de Vénus (figure 2).

Figure 2. - Positions apparentes de Vénus sur le disque solaire.

On vérifie bien que cette différence est égale à D2 D1 lorsque les quatre points sont coplanaires cest-à-dire lorsque V1, V2 et C sont alignés.

La valeur que les observateurs vont mesurer est donc la distance Dp entre les centres apparents de Vénus et cest la relation qui va nous permettre de calculer les parallaxes.

Pour cela nous allons exprimer les deux parallaxes en fonction des distances entre le centre de la Terre et le centre des deux astres. Soient rV la distance entre le centre du Soleil et le centre de Vénus et rT la distance entre le centre de la Terre et le centre du Soleil, la distance Vénus-Terre est donc égale à rTrV. Pour exprimer cette parallaxe nous devons également connaître la projection d de la distance entre les deux points M1 M2 sur le plan normal à la direction Terre Soleil (figure 3).

Figure 3. - Parallaxe solaire relative aux points M1 et M2.

Comme le rayon terrestre et la distance entre les deux points sont petits par rapport aux distances Terre-Soleil et Terre-Vénus, les parallaxes sont données par les formules approchées suivantes :

           (1)

En réalité les parallaxes exactes sont données par :

On a donc la relation suivante :

           (2)

Et

ou encore

           (3)

La mesure nous donne la valeur Dp exprimée en diamètre solaire, on doit également mesurer le diamètre du Soleil, car si la distance Terre-Soleil est inconnue on ne peut pas la calculer.

Pour connaître la parallaxe solaire, il faut donc également connaître le rapport des distances Soleil-Terre et Soleil-Vénus. Or ce rapport peut être calculé grâce aux lois de Kepler.

III. Le calcul du rapport des distances au Soleil à laide des lois de Kepler

La première loi de Kepler nous dit que les planètes décrivent des orbites elliptiques autour du Soleil et que le Soleil occupe un des foyers de ces ellipses. A un instant donné le rayon vecteur rp joignant le centre du Soleil à une planète p se calcule à laide de la formule suivante :

           (4)

ap est le demi-grand axe de lellipse, ep est lexcentricité de lellipse et E est un angle appelé anomalie excentrique qui permet de placer la planète sur son orbite.

La troisième loi de Kepler fournit une relation entre les demi-grands axes des orbites et les périodes de révolution des planètes ; ainsi pour un même corps central toutes les orbites des planètes qui gravitent autour de ce corps central vérifient la relation suivante :

           (5)

Les lois de Kepler décrivent donc les orbites du système solaire à un facteur déchelle près. Lobservation des périodes de révolution des planètes nous donne les rapports des demi-grands axes, ainsi le rapport des demi-grands axes des orbites de Vénus et de la Terre est égal à :

           (6)

et à un instant t quelconque le rapport des rayons vecteurs est égal à

           (7)

Donc les lois de Kepler permettent de calculer le rapport des rayons vecteurs pour un instant t quelconque.

Notre mesure nous permet de calculer la valeur ps, il convient donc maintenant de passer de cette valeur à la valeur de la parallaxe équatoriale moyenne du Soleil p0 .

IV. Calcul de la parallaxe moyenne du Soleil

La parallaxe équatoriale moyenne du Soleil p0 est par définition langle sous lequel on voit le rayon équatorial de la Terre depuis le centre du Soleil lorsque le Soleil se trouve à une unité astronomique de la Terre.

On a donc la relation suivante :

           (8)

R étant le rayon équatorial terrestre et a lunité astronomique.

Léquation (1) nous donne la valeur de la parallaxe solaire ps en fonction de la distance rT Terre-Soleil et de la projection d de la distance entre les points dobservations sur le plan normal à la direction Terre-Soleil.

Il suffit dexprimer cette distance d en rayon terrestre et la distance Terre-Soleil en unité astronomique pour avoir une relation entre ps et p0.

           (9)

Il ne reste plus quà calculer le rapport d sur R.Le rapport a/rT nous est fourni par la première loi de Kepler (cf. formule 4). Si lon fait le produit vectoriel des deux vecteurs et on obtient :

           (10)

Or le produit de la longueur du premier vecteur par le sinus de langle entre les deux vecteurs est égal à la distance d. De même la longueur de est égale à la distance rT(figure 4).

Figure 4. - . Parallaxe solaire relative aux points M1 et M2.

La résolution de léquation (10) nous donne la valeur de d.

           (11)

Remarque : si la notion de produit vectoriel nest pas connue, on peut utiliser le produit scalaire des mêmes vecteurs, cela permet de calculer le cosinus de langle, puis son sinus par la relation : .

Description de ce calcul

Ce calcul sur les vecteurs demande de connaître les coordonnées cartésiennes des deux points M1 et M2 et du centre du Soleil C dans un repère orthonormé (O, x, y, z) centré au centre de la Terre. Nous allons utiliser le repère équatorial apparent géocentrique pour ce calcul.

Ce repère est défini par le plan de léquateur terrestre à linstant t de lobservation (plan Oxy) et par la direction du pôle céleste nord de laxe de rotation de la Terre (Oz). Dans ce repère on peut définir un système de coordonnées cartésiennes (x, y, z) et un système de coordonnées polaires (a, d, r) les deux angles portent le nom dascension droite et de déclinaison (figure 5). On passe dun système à lautre par les relations suivantes :

           (12)

et les relations inverses

           (13)

La direction de laxe Ox à linstant t est la direction de léquinoxe de printemps au même instant.

Les éphémérides (cest-à-dire les lois de Kepler) nous donnent les coordonnées équatoriales géocentriques du centre du Soleil (a, d) la distance nest pas connue mais cela na pas dimportance car on peut remplacer le vecteur par son vecteur unitaire dans léquation 11.

Le problème le plus complexe est la détermination des coordonnées cartésiennes des points M1 et M2 dans ce repère équatorial.

Figure 5. - . Coordonnées équatoriales géocentriques.

Les positions dun point de la surface terrestre sont données par sa latitude et sa longitude géographique, la latitude est donnée par rapport à léquateur terrestre cest donc une variable angulaire identique à la déclinaison, la longitude est donnée par rapport à un méridien origine (méridien de Greenwich) cest donc une variable angulaire identique à lascension droite, mais qui a une origine différente de celle des coordonnées équatoriales célestes. Il convient donc de connaître à chaque instant langle entre la direction de laxe Ox et la direction de la projection du méridien origine dans le plan de léquateur (cf. figure 5). Cet angle est lié à la rotation de la Terre sur elle-même, il porte le nom de temps sidéral au méridien de Greenwich et il varie de 360° en 23h 56m 4s (révolution sidérale de la Terre).

Il suffit donc de connaître le temps sidéral à Greenwich TG à 0h UT le jour du passage pour connaître le temps sidéral à Greenwich à linstant t puis le temps sidéral en tout point de la Terre de longitude l.

           (14)

On passe du temps sidéral à Greenwich au temps sidéral au lieu M de longitude l en ajoutant ou en retranchant cette longitude.

Attention, le temps sidéral augmente lorsque lon séloigne vers lest du méridien de Greenwich, il convient donc de bien faire attention à la convention de signe utilisée pour noter les longitudes.

Si les longitudes sont comptées négativement vers lest alors la relation liant le temps sidéral local au méridien du lieu de longitude l et le temps sidéral au méridien de Greenwich est la suivante :

           (15)

Attention, les deux angles doivent être exprimés avec la même unité (degrés ou heures).

Alors les coordonnées cartésiennes dun point M1 de coordonnées géographiques (j1, l1) à linstant t sont données par :

           (16)

La longueur du vecteur (son module) et ses coordonnées (X, Y, Z) sont donnés par :

           (17)

Le vecteur unitaire de la direction « centre de la Terre-Soleil » est donné par :

           (18)

Le produit vectoriel et son module sont alors :

           (19)

et finalement en utilisant la formule 11, on obtient :

           (20)

Et la parallaxe équatoriale moyenne p0 est donnée par daprès (9) :

           (21)

V. Application numérique

Nous allons prendre pour exemple lobservation Antananarivo (Madagascar) et à Helsinki (Finlande) à linstant t=8h 30min le 8 juin 2004.

Les coordonnées géographiques de Antananarivo sont les suivantes :

Latitude :18° 52' sud, longitude : 47° 30' est donc j1 = –18,866667° et l1 = –47,5°.

Les coordonnées géographiques dHelsinki sont les suivantes :

Latitude :60° 8' nord, longitude : 25° 3' est donc j2 = 60,133333° et l2 = –25,05°.

Les coordonnées équatoriales géocentriques du Soleil à 8h 30m UTC sont données par les éphémérides :

Ascension droite du Soleil as = 76°49' 36.493"

Déclinaison du Soleil ds = +22°53' 16.237"

Le temps sidéral à Greenwich à un instant t en UTC est donné par la formule suivante :

TG(t UTC) = 17h 6m 51,31s + 1,002737908 t

Donc le temps sidéral à Greenwich à 8h 30min est égal à :

TG = 17h 6m 51,31s + 8h 31m 23,78s = 25h 38m 15,09s = 1h 38m 15,09s

Il convient de le convertir en degrés avant de lutiliser pour calculer le temps sidéral local au deux lieux considérés.

TG = 1h 38m 15,09s = 24,562875°.

Doù on déduit le temps sidéral local à 8h 30m à Antananarivo :

Tl1 = 24,562875 – (–47,5°) = 72,062875°

Et le temps sidéral local à 8h 30m à Helsinki :

Tl2 = 24,562875 – (–25,05°) = 49,612875°

On en déduit les coordonnées cartésiennes équatoriales des deux villes :

Antananarivo :

Helsinki :

Les coordonnées du vecteur unitaire de la direction Terre-Soleil sont obtenues à laide de la formule (18) :

Le vecteur a pour coordonnées :

La formule (20) nous permet de calculer la valeur de d :

Les éphémérides nous donnent les rapports des rayons vecteurs ainsi que le rapport de la distance Terre-Soleil sur le demi-grand axe de lorbite terrestre à linstant considéré :

Il ne reste plus quà faire une hypothèse sur les valeurs mesurées, cest-à-dire sur Dp et sur le diamètre solaire :

Nous allons faire les hypothèses suivantes :

Ce qui donne pour valeur de Dp : 28,359"

La formule 3 nous donne la valeur de la parallaxe solaire :

Et la formule (21) nous donne la valeur de la parallaxe équatoriale moyenne :

La valeur que lon trouve est relativement proche de la réalité, mais elle repose uniquement sur la mesure de la distance des centres apparents de Vénus sur le disque solaire et la grandeur du diamètre solaire. La taille apparente du diamètre solaire peut être mesurée avec une bonne précision, par contre la mesure de la distance entre les centres apparents de Vénus nest pas évidente, sur un cliché photographique classique, le diamètre apparent est de lordre 20mm, la distance des centres est alors de 0,3mm et une précision de lordre du millième correspond à une mesure à 0,02mm près. On remarque également que la valeur Dp : 28,359" est inférieure au rayon apparent (28,884") de Vénus. Les deux disques apparents de Vénus sont donc partiellement superposés.

Dans les formulaires précédents on a occulté un certain nombre de difficultés pour simplifier le problème. Voici la liste des complications qui apparaissent si lon veut faire un calcul rigoureux :

  • En raison des perturbations mutuelles, les orbites desplanètes ne suivent pas les lois de Kepler (valables uniquement pour deuxcorps) mais des trajectoires plus complexes.
  • Ce nest pas la Terre qui a une orbite quasi-elliptique autourdu Soleil mais le barycentre du système Terre-Lune.
  • Suite au mouvement de laxe de rotation de la Terre(précession et nutation), lorigine Ox du repère équatorial nest pasfixe dans le temps.
  • La lumière se propageant avec une vitesse finie, les positionsdu Soleil et de Vénus à un instant t ne sont pas des positionsgéométriques, mais celles des deux corps au instants t – tp, tp représentant le temps mis par la lumière pourparcourir la distance entre chaque corps et la Terre. Comme ces distances nesont pas supposées connues, il convient de réitérer les calculs pour en tenircompte.
  • Nous avons supposé la Terre comme sphérique, en réalité elleest aplatie.

VI. Calcul de la parallaxe à partir des instants des contacts ou de la durée du passage.

Nous avons vu dans la fiche n°04b, quil existe deux formules simplifiées qui permettent un calcul direct de la parallaxe à partir de la comparaison des instants dun même contact en deux lieux distincts (méthode de Delisle) ou à partir de la comparaison de la durée des passages en deux lieux distincts (méthode de Halley).

Nous allons traiter simultanément ces deux aspects à partir de lexemple numérique précédent.

La parallaxe équatoriale moyenne solaire sobtient en comparant deux contacts identiques à laide de la formule simplifiée suivante (cf. formule 16 de la fiche n°06) :

           (22)

Si lon néglige les incertitudes et les erreurs alors cette formule devient :

           (23)

De même la parallaxe équatoriale moyenne solaire sobtient en comparant deux durées identiques à laide de la formule suivante (cf. formule 21 de la fiche n°06) :

           (24)

i et j sont des indices liés aux mêmes contacts : i = 1, j = 4 pour les contacts extérieurs et i = 2, j = 3 pour les contacts intérieurs.

Les coefficients A, B, C et le terme sont calculés pour chaque contact et sont données par le tableau suivant :

Description du contact A B C dD/dt
"/min
Premier contact extérieur (indice 1) 2,2606 -0,0194 1,0110 -3,0846
Premier contact intérieur (indice 2) 2,1970 0,2237 1,1206 -2,9394
Dernier contact intérieur (indice 3) -1,0929 -1,1376 1,9090 2,9391
Dernier contact extérieur (indice 4) -0,9799 -1,3390 1,8383 3,0842

VII. Exemples numériques

Nous allons reprendre lexemple des deux villes précédentes avec les hypothèses dobservations suivantes :

Ville n°1 : Antananarivo (j1 = –18,866667° et l1 = –47,5°)

Instant du premier contact intérieur observé (indice 2): t2 = 5h 35m 30s UTC.
Instant du dernier contact intérieur observé (indice 3): t3 = 11h 8m 4s UTC.
Durée du passage intérieurobservée: t3- t2 = 5h 32m 34s.

Ville n°2 : Helsinki (j2 = 60,133333° et l2 = –25,05°)

Instant du premier contact intérieur observé (indice 2): t2 = 5h 38m 38s UTC.
Instant du dernier contact intérieur observé (indice 3): t3 = 11h 2m 20s UTC.
Durée du passage intérieurobservée: t3- t2 = 5h 23m 42s.

Dans les formules (22) et (23) les facteurs des coefficients A, B, C sont identiques et peuvent être calculés séparément :

Calcul de la parallaxe à laide des premiers contacts :

Lécart des temps des premiers contacts intérieurs est de –3m 8s (–3,1333m), et lusage des valeurs des coefficients A2, B2, C2 et dans la formule (22) nous donnent la relation suivante :

Calcul de la parallaxe à laide des durées des passages intérieurs :

Lécart de durée des passages intérieurs est de 8m 52s (8,866m), et lusage des valeurs des coefficients A2, B2, C2, A3, B3, C3 et dans la formule (23) nous donnent la relation suivante :

Attention, cest la valeur et surtout son signe qui doit être utilisée.

On rappelle que ces méthodes ne sont pas exactes, et que lon doit utiliser des formules plus complexes pour réduire les observations.

Contact

Pour tous renseignements concernant ces pages, veuillez contacter ou le Service de renseignements de l'IMCCE.
Caution : This Website was created with the Ministère de l'Education Nationale, the CNRS and the CNES support. Any use of the data published on this website requires the IMCCE agreement..