L’article ci-dessous est la suite de « Le nocturlabe » article paru dans le N°41 - septembre 2014 de MathemaTICE.
Remerciement : je tiens à exprimer toute ma gratitude à Aymeric Picaud pour son travail de correction en vue de rendre l’article plus attrayant.
Étude de l’influence de la latitude et de la précession des équinoxes sur la lecture de l’heure aux étoiles.
À Denis Savoie pour son étude très fine du nocturlabe
1) Introduction
Je dois l’écriture de cet article à l’une de mes lectures, un document produit parDenis SAVOIE, astronome, qui m’a donné l’aimable autorisation de reproduire ses calculs de trigonométrie sphérique afin d’évaluer l’erreur commise lors de l’utilisation du nocturlabe.
Le dossier précédent sur le nocturlabe nous avait appris à mesurer l’heure aux étoiles en admettant que l’étoile polaire est sur l’axe du monde. Nous nous proposons d’établir le décalage dû au fait que l’étoile polaire est plus ou moins éloignée de l’axe du monde.
Au milieu du seizième siècle elle était approximativement située à 3° du pôle nord céleste et en sera au plus proche en 2100.
Nous verrons aussi l’influence de la latitude sur la précision des mesures. L’instrument devient de moins en moins précis au fur et à mesure que l’on monte en latitude.
La trigonométrie sphérique sera sollicitée par le biais des formules de Gauss.
Le logiciel Stellarium nous permettra de connaître les cordonnées équatoriales des étoiles à l’époque que l’on choisira et le tableur EXCEL fera le reste.
Il est important de préciser que la lecture de cet article suppose connues du lecteur les notions fondamentales d’astronomie que j’ai développées dans le dossier précédent sur le nocturlabe.
Notre étude visera à calculer :
Dh= angle horaire de l’étoile de référence- valeur lue sur le nocturlabe.
Nous visualisons ici le nocturlabe présenté dans Le nocturlabe paru dans le N°41 - septembre 2014 de MathemaTICE.
Schéma 1
Ce nocturlabe était basé sur le fait que l’étoile polaire est assez proche du pôle nord céleste au sens où, si l’on désigne par P l’étoile polaire, N le pôle nord terrestre ou céleste et O le centre de la Terre alors l’angle $\widehat{PON}$ est très petit.
Cet angle était de 3° au milieu du seizième siècle.
On pourra le vérifier avec par exemple le goniomètre de Stellarium.
Notons qu’on peut choisir indifféremment le centre de la Terre ou l’observateur, compte tenu des distances vertigineuses qui nous séparent des étoiles comparées aux dimensions de la Terre.
Le fait que l’étoile polaire ne soit pas rigoureusement sur l’axe du monde entraîne une erreur de lecture sur l’angle horaire et par voie de conséquence sur l’heure :
Sur le schéma 1, l’angle $\widehat{TPA}$ correspond à l’heure et l’angle $\widehat{APB}$ à l’angle horaire.
Voici à présent le schéma simplifié du nocturlabe.
- p est le pôle nord ; o centre de la sphère.
- [et] est le manche du nocturlabe.
- z : l’observateur
- e : l’étoile polaire au centre du nocturlabe
- E : sa projection sur la sphère
- a : l’étoile de référence
- A : sa projection sur la sphère.
Schéma 2
Rappel : l’orientation du manche du nocturlabe [et] se fait en considérant visuellement le plan contenant le manche, de telle manière que ce plan soit perpendiculaire au plan horizontal de l’observateur et qu’il passe par l’observateur.
Ce plan n’est plus le plan méridien de l’observateur comme dans le N° 41 consacré au nocturlabe.
Nous lisons alors un angle $\widehat{aet}= \widehat{Hp}$ sur le nocturlabe, mais cet angle n’est plus l’angle horaire de l’étoile de référence A du fait que la droite en vert n’est plus l’intersection du plan méridien avec le plan du nocturlabe.
Notre étude visera à calculer :
Dh= angle horaire de l’étoile de référence- Hp
2) Exercices
Exercice (shéma2) : Montrer que si l’on appelle E’ la projection orthogonale de l’étoile polaire sur le plan horizontal de l’observateur z alors ce plan (Etz) est aussi le plan (EE’z).
Exercice : Montrer que dans le schéma 3 l’angle sphérique $\widehat{AoB}$ est égal à la mesure de l’arc A’B’.
Schéma 3
Exercice
Soit une sphère de centre O ; (D’’) une droite passant par O ; (S) un plan perpendiculaire à (D’’).
z, p, a trois points de la sphère. Z, P, A les intersections respectives des droites (Oz) , (Op) et (Oa) avec le plan (S).
Soit (D) la droite d’intersection du plan (OPZ) avec le plan (S).
(D’) la droite d’intersection du plan (OPA) avec le plan (S).
Il s’agit de montrer que l’angle $\widehat{zpa}$ du triangle sphérique est égal à l’angle $\widehat{ZPA}$.
Schéma 4
Application : l’angle $\widehat{AEt}$ du nocturlabe est égal à l’angle du triangle sphérique aez (shéma2).
Nous allons dans ce document essayer de mieux comprendre les valeurs affectées aux cotés des triangles sphériques qui ont servi à établir les formules en vue du calcul de Dh.
3) Rappels sur les notions fondamentales d’astronomie
Les notions de latitude, de hauteur d’un astre, d’angle horaire, de coordonnées équatoriales, de point vernal, ont été exposées dans l’article n° 41.Nous les revoyons ici en nous concentrant sur certaines portions d’arc afférentes à ces notions.
1) la latitude : La mesure de l’arc $\overset{\frown}{ms}$ est la latitude de l’observateur et la mesure de l’arc $\overset{\frown}{Pm}$ que l’on appelle la colatitude est égale à 90°-latitude.
m désigne l’observateur, le méridien de m coupe le cercle équateur en s.
La latitude est la mesure en degrés de l’arc $\overset{\frown}{ms}$ , cette mesure étant affectée du signe + dans l’hémisphère nord et – dans l’hémisphère sud pour les formules mais on peut rajouter à cette mesure N POUR l’hémisphère nord et S pour l’hémisphère sud.
2) déclinaison d’un astre : Dans ce cas la mesure de l’arc $\overset{\frown}{PE}$ est égale à 90° -déclinaison de l’astre.
Il s’agit ici de rappeler ce qui a été vu dans Mathematice n°41.
La mesure en degrés de $\overset{\frown}{ES}$ est la déclinaison de l’astre qui est ici représenté sur la sphère céleste en E. On adopte les mêmes conventions que pour la latitude.
3) hauteur d’un astre
- Schéma 1
- O est l’observateur centre de la sphère céleste locale
- E désigne une étoile
- Z le zénith
- Le plan en gris est le plan équatorial
Masquons la sphère céleste (en violet).
Nous obtenons :
- Schéma 2
- e’ est la projection de l’étoile sur la plan horizontal.
Par définition la hauteur de l’étoile est $\widehat{EOe’}$.
Le schéma montre alors que la hauteur est complémentaire de $\widehat{EOZ}$.
4) Les triangles sphériques qui nous seront utiles
Voir la signification des lettres au niveau du schéma 2 du paragraphe 1.La sphère représentée ci-dessous peut être vue comme la sphère terrestre. Dans ce cas Z représente l’observateur.
On peut la voir aussi comme la sphère céleste et dans ce cas Z représente le zénith.
Schéma 5 l’explication des valeurs des arcs est donnée au paragraphe 3.
φ représente la latitude de z
δ’ la déclinaison de l’étoile polaire
δ la déclinaison de l’étoile de référence ;
K la distance angulaire entre l’étoile polaire et l’étoile de référence
h la hauteur de l’étoile de référence et p la mesure en degrés de l’arc ez
À présent nous pourrons appliquer les formules de la trigonométrie sphérique.
5) Démonstration des formules utilisées dans le fichier Excel intitulé : calcul d’erreur
D’après un calcul de Denis SAVOIE.
a ) RAPPEL :
a : étoile repère par exemple Dubhé ;
e : étoile polaire ;
P : pôle nord.
(ea) : direction dans laquelle est vue l’étoile repère.
Z : observateur
Hp : angle entre le manche du nocturlabe et la direction (ea) mesuré de 0° à 360° dans le sens des aiguilles d’une montre à partir du manche.
Notations : (et) désignent les ascensions droites de l’étoile de référence (par exemple Dubhé) et de l’étoile polaire. Et leurs déclinaisons respectives.
b) Application de la formule de GAUSS ( iii )aux différents triangles sphériques du schéma 5 du paragraphe 4 qui par commodité sont représentés comme ci-dessous.
Formules de Gauss.
Triangle sphérique EPA
$ \cos K= \sin \delta \sin \delta’ + \cos \delta \cos \delta’ \cos \widehat{EPA} $ (1)
Triangle sphérique APZ
$\sin h= \sin \varphi \sin \delta + \cos \varphi \cos \delta \cos \widehat{APZ} $ (2)
Triangle sphérique EPZ
$\cos p= \sin \varphi \sin \delta’ + \cos \varphi \cos \delta’ \cos \widehat{EPZ} $ (3)
Triangle sphérique AEZ
$\sin h=\cos K \cos p +\sin K \sin p \cos \widehat{AEZ} $ (4)
D’où
$\cos \widehat{AEZ} = \dfrac{\sin h -\cos K \cos p }{\sin K \sin p }$ (5)
Compte tenu de la périodicité de la fonction cosinus, on aura dans tous les cas de figure possibles :
$\cos \widehat{EPA} = \cos (\alpha - \alpha’)$
avec α et α’ ascensions droites de l’étoile repère (par exemple Dubhé) et de l’étoile polaire.
$\cos \widehat{APZ}= cos H $
H désignant l’angle horaire de l’étoile repère, par exemple Dubhé.
$\cos \widehat{EPZ} = \cos(\alpha - \alpha’ + H)$
$\cos EZ = \cos Hp $
où Hp a été défini au rappel 1.
D’où, si l’on tient compte de ce que 0 < p < π
$\sin h = \sin \varphi \sin \delta + \cos \varphi \cos \delta \cos H$
$\cos K = \sin \delta \sin \delta’ + \cos \delta \cos \delta’ \cos (\alpha - \alpha’)$
$\cos p = \sin \varphi \sin \delta’ + \cos \varphi \cos \delta’ \cos (\alpha - \alpha’ + H)$
$\sin h = \cos K \cos p +\sin K \sin p \cos Hp $
$\cos Hp= \dfrac{\sin h -\cos K \cos p }{\sin K \sin p }= \dfrac{\sin h -\cos K \cos p }{\sin K \sqrt{1-(\cos p)^2 }}$
6) La précession des équinoxes
Voir le paragraphe 2 de l’article Position, mouvement et distance des étoiles
7) Conclusion
Voici les résultats obtenus à l’aide du fichier Excel
pour l’étoile de référence Dubhé à une latitude moyenne de 49° en 2010.
angle horaire de | Dh en min |
---|---|
l’étoile de référence | |
0,000 | 1,568496981 |
30,000 | 2,832066826 |
60,000 | 4,431863665 |
90,000 | 5,95881666 |
120,000 | 7,020441416 |
150,000 | 7,328865022 |
180,000 | 6,780716381 |
210,000 | 5,505138482 |
240,000 | 3,847601774 |
270,000 | 2,274077852 |
300,000 | 1,223626299 |
330,000 | 0,972921361 |
Voici les résultats obtenus à l’aide du fichier Excel pour l’étoile de référence Dubhé à une latitude moyenne de 49° en 1550
H deg | Dh min |
---|---|
0 | 1,568496981 |
30 | 11,85983404 |
60 | 18,87410676 |
90 | 24,17262312 |
120 | 26,32826029 |
150 | 24,33295085 |
180 | 18,23701829 |
210 | 9,687611361 |
240 | 1,544145872 |
270 | -3,54765508 |
300 | -4,381656206 |
330 | -1,260748701 |
L’erreur de lecture est en nette augmentation.
Le fichier Excel permettra de faire varier la latitude et d’observer combien le décalage devient important aux hautes latitudes et de constater que plus on se rapproche de 2100, plus la distance entre l’étoile polaire et le pôle nord céleste diminue, ce qui a pour conséquence de réduire le décalage Dh.
Le changement de date permet d’être confronté à la précession des équinoxes qui modifie la position de l’étoile polaire sur la sphère céleste.
En effet, une étoile polaire est une étoile se trouvant approximativement sur l’axe de rotation de la Terre.
Ce fichier Excel comporte trois onglets :
- un onglet relatif à l’étoile Dubhé en 2010
- Un onglet relatif à l’étoile Dubhé en 1550
- Un onglet relatif à une étoile circumpolaire de vitre choix à la date que vous aurez choisie sur Stellarium.
L’étoile Kochab est intéressante pour notre époque alors que vers 1550, l’erreur pouvait atteindre 40 min.
Exercice : le vérifier en utilisant le fichier Excel après avoir récupéré les coordonnées équatoriales sur STELLARIUM ou sur le logiciel REDSCHIFT.
Pour conclure, nous avons montré que la latitude, dès lors qu’elle est élevée (à partir de 70°), et la précession des équinoxes affectent de manière notable la précision de l’instrument, celui-ci devenant de plus en plus précis au fur et à mesure que l’on se rapprochera de l’an 2100, date à laquelle la polaire sera au plus prés du pôle.
Par le plus heureux des hasards, certaines étoiles circumpolaires n’ont pas affecté la précision du nocturlabe, ce qui n’était pas le cas de toutes les étoiles circumpolaires comme par exemple l’étoile Dubhé qui en 1550 était à l’origine d’un décalage de 24 min environ.
Rendons hommage à l’ingéniosité humaine qui a permis, avec un mécanisme aussi simple de réaliser la prouesse de lire l’heure dans les étoiles !