Équation du temps

L’équation du temps est un paramètre utilisé en astronomie pour rendre compte du mouvement apparent relatif du Soleil par rapport au soleil moyen, lesquels peuvent différer l'un par rapport à l'autre de plus ou moins un quart d'heure environ. D'une année sur l'autre, la courbe d'évolution annuelle de ce paramètre se répète quasiment à l'identique. La connaissance de l'équation du temps donne le moyen de corriger à tout instant l'heure donnée par un cadran solaire pour trouver l'heure légale, d'écoulement uniforme. Autrefois, elle permettait de contrôler la marche d'une horloge, à écoulement théoriquement uniforme, par rapport aux indications d'un cadran solaire, notamment au moment du midi vrai, alors important socialement, moment repéré sur un cadran ou une méridienne.

Résultant des caractéristiques du mouvement de la Terre autour du Soleil, l'équation du temps peut se calculer très précisément. On en trouve des tables détaillées dans les éphémérides astronomiques^[1].

Remarque sur le mot « équation » : en astronomie ancienne, le terme « équation » désignait une correction ajoutée algébriquement à une valeur moyenne pour obtenir une valeur vraie. C'est une telle acception qui a survécu dans l'expression « équation du temps », et qui se retrouve aussi dans « équation du centre » ou « équation des équinoxes ». Il s'agit bien d'un paramètre, et non d'une équation au sens habituel du terme (égalité avec inconnues, comme c'est le cas d'une équation polynomiale ou d'une équation différentielle).

Définition

Par convention, l’équation du temps, à un instant donné, est la différence entre le temps solaire moyen et le temps solaire vrai^[2],[3].

• Le temps solaire moyen est basé sur le soleil moyen, défini comme un objet qui, tout au long de l'année, se déplacerait sur l'équateur à une vitesse constante, telle que la durée du jour solaire moyen soit de 24 heures exactement.
• Le temps solaire ou temps vrai est une mesure du temps basée sur le soleil vrai, tel que donné par un cadran solaire. En particulier, le midi solaire correspond à l'instant de la journée où le Soleil atteint son point le plus élevé dans le ciel.

Une valeur positive de l'équation du temps indique que le soleil vrai est en retard sur le soleil moyen, c'est-à-dire plus à l'est, et une valeur négative qu'il est en avance, c'est-à-dire plus à l'ouest. Par exemple, lorsque l'équation du temps vaut + 8 minutes, cela signifie qu'il est 12 h 08 du temps solaire moyen lorsque le cadran solaire indique midi vrai.

C'est du moins la convention de signe utilisée en France, où l’équation du temps est l'équation du temps vrai, c'est-à-dire ce qu'il faut ajouter au temps vrai pour obtenir le temps moyen. Dans certains pays, tels que le Royaume-Uni, les États-Unis ou la Belgique, l'équation du temps est souvent définie avec la convention de signe inverse : c'est l'équation du temps moyen, c'est-à-dire la quantité qu'il faut ajouter au temps moyen pour obtenir le temps vrai. Les deux variables, « équation du temps vrai » et « équation du temps moyen » ont des valeurs opposées.

Autre forme de la définition : l’équation du temps, à chaque instant, est la différence entre l'ascension droite du soleil moyen et celle du soleil vrai.

Évolution annuelle de l'équation du temps

Équation du temps (courbe rouge)

L'évolution de l’équation du temps sur une année complète est représentée par la courbe rouge sur la figure ci-contre. En première approximation, sa forme s'analyse comme résultant de la superposition de deux sinusoïdes :

• en bleu sur le diagramme : une sinusoïde de période égale à un an, d'amplitude égale à 7,66 minutes et s'annulant aux passages de la Terre aux apsides : périgée le 3 janvier et apogée début juillet. Cette composante reflète l'excentricité de l'orbite terrestre ;
• En vert sur le diagramme, une sinusoïde de période égale à une demi-année, d'amplitude 9,87 minutes et s'annulant aux solstices et aux équinoxes. Cette composante résulte de l'obliquité de l'écliptique sur l'équateur.

L'équation du temps, en rouge, s'annule quatre fois par an, vers le 15 avril, le 13 juin, le 1^er septembre et le 25 décembre. Son maximum, atteint vers le 11 février, vaut 14 min 15 s, et son minimum, atteint vers le 3 novembre, vaut − 16 min 25 s.

Analemme

Exemple d'analemme, pour Greenwich en 2006.
Latitude 51,4791° nord, colatitude = 38,5209°, longitude 0°^[4]

L'évolution annuelle de l'équation du temps, en un lieu donné, peut être visualisée à l'aide d'une courbe appelée analemme ou courbe en 8, définie comme suit : chaque point de cette courbe représente une position du soleil (vrai) lorsqu'il est 12 h pour le soleil moyen, c'est-à-dire lorsque ce dernier passe au centre du diagramme. Les axes sont les suivants, avec des échelles différentes, de façon à mieux mettre en évidence la légère asymétrie de la courbe :

• l'axe horizontal représente l'azimut en degrés (180° correspond au sud). L'équation du temps se lit le long de cet axe, donc comme l'écart horizontal par rapport à la ligne 180°. Avec la convention de signe adoptée, elle est positive à gauche de la ligne 180°. La correspondance entre l'angle et le temps est 360° = 24 h, donc 1° = 4 minutes ;
• l'axe vertical représente la hauteur du Soleil au-dessus de l'horizon, liée aux variations de sa déclinaison ;
• le soleil moyen, midi moyen, se trouve au centre du diagramme (azimut = 180°, hauteur = 90° - latitude du lieu).

Sur l'exemple ci-contre^[5], le premier jour de chaque mois est affiché en noir, et les positions des solstices et équinoxes sont affichées en vert. On lit par exemple :

• le 3 novembre, avance maximale du soleil vrai sur le soleil moyen, et l'équation du temps, qui est négative, vaut - 16 min 23 s ;
• le 12 février, le retard est maximal, et l'équation du temps, qui est positive, vaut + 14 min 20 s ;
• au solstice d'hiver, vers le 21 décembre, la hauteur du Soleil est minimale et vaut 15,08° (hauteur au solstice d'hiver = colatitude du lieu - obliquité de l'écliptique = 38,52° - 23,44°) ;
• au solstice d'été, vers le 21 juin, la hauteur est maximale et vaut 61,96° (hauteur au solstice d'été = colatitude du lieu + obliquité de l'écliptique = 38,52°}} + 23,44°) ;
• aux équinoxes, le Soleil passe dans le plan de l'équateur et a, à ce moment-là, la même hauteur que le soleil moyen, égale à la colatitude du lieu.

Certains cadrans solaires sont munis d'un analemme. Ils peuvent même donner directement le temps moyen, soit parce que les droites horaires sont transformées en courbes corrigées de l'équation du temps, soit parce que le gnomon a reçu une forme tenant compte de cette correction. Dans les deux cas, il faut tenir compte de la période de l'année ou disposer de deux cadrans.

Homonymie

Il ne faut pas confondre cet analemme avec la figure du même nom, qui lui est historiquement bien antérieure^[6], et qui servait à tracer des cadrans solaires ou établir géométriquement la hauteur du Soleil. Elle résultait de la projection de la sphère céleste sur le plan méridien.

Variation de cette évolution au cours du temps

La forme de la courbe « équation du temps », c'est-à-dire la valeur des extrema et les instants où on les observe, ainsi que les instants où la courbe s'annule, évoluent très lentement au cours des années pour au moins deux raisons :

1. la Terre dans son mouvement autour du Soleil subit l'influence des autres planètes du système solaire, ce qui entraîne une variation de l'excentricité de son orbite, ainsi qu'une lente rotation de la ligne joignant le périhélie à l'aphélie de l'orbite, appelée ligne des apsides ;
2. la Terre, dans sa rotation sur elle-même, subit l'influence du couple (Lune, Soleil), ce qui entraîne une variation de son obliquité en inclinaison et direction ; ces phénomènes sont connus et décrits sous le nom de nutation en longitude, nutation en obliquité et précession des équinoxes.

Ces évolutions provoquent notamment un glissement relatif des dates des passages aux apsides par rapport à celles des solstices et des équinoxes, qui sont fixes par construction de l'année tropique. Sur une durée de 70 siècles, de l'an - 2000 à + 5000, les extrema sont définis par le tableau suivant^[7] :

Année	Premier maximum	Premier minimum	Deuxième maximum	Deuxième minimum
- 2000	+ 18 min 33 s, 31 janvier	- 12 min 45 s, 20 mai	+ 2 min 06 s, 10 août	- 9 min 30 s, 26 octobre
- 1000	+ 18 min 18 s, 3 février	- 10 min 14 s, 21 mai	+ 2 min 06 s, 6 août	- 11 min 45 s, 27 octobre
0	+ 17 min 27 s, 6 février	- 7 min 44 s, 20 mai	+ 2 min 57 s, 1er août	- 13 min 45 s, 29 octobre
+ 1000	+ 16 min 04 s, 9 février	- 5 min 27 s, 18 mai	+ 4 min 30 s, 29 juillet	- 15 min 20 s, 1er novembre
+ 2000	+ 14 min 15 s, 11 février	- 3 min 41 s, 14 mai	+ 6 min 30 s, 26 juillet	- 16 min 25 s, 3 novembre
+ 3000	+ 12 min 08 s, 14 février	- 2 min 37 s, 10 mai	+ 8 min 41 s, 25 juillet	- 16 min 57 s, 6 novembre
+ 4000	+ 9 min 52 s, 15 février	- 2 min 24 s, 6 mai	+ 10 min 48 s, 25 juillet	- 16 min 54 s, 9 novembre
+ 5000	+ 7 min 38 s, 15 février	- 3 min 00 s, 3 mai	+ 12 min 38 s, 26 juillet	- 16 min 17 s, 12 novembre

Analyse intuitive de l'équation du temps

Dans cette section, on examine successivement deux effets séparément :

• Influence de l'ellipticité de l'orbite de la Terre. Pour cela, on supposera que l'obliquité est nulle, c'est-à-dire que l'axe de rotation de la Terre est perpendiculaire à l'écliptique.
• Influence de l'obliquité de la Terre. Pour cela, on supposera que l'ellipticité de l'orbite de la Terre est nulle, c'est-à-dire que l'orbite est circulaire.

Influence de l'ellipticité de l'orbite de la Terre

La seconde loi de Kepler (loi des aires) indique comment la vitesse de la Terre varie le long de son orbite : elle est maximale au périhélie (30,287 km/s vers le 3 janvier), minimale à l'aphélie six mois plus tard (29,291 km/s début juillet), et égale à sa valeur moyenne (vers début avril et début octobre).

D'octobre à mars (demi-orbite voisine du périhélie), la vitesse de la Terre sur son orbite est plus élevée que sa valeur moyenne. Donc, vu de la Terre, le soleil vrai va plus vite que le soleil moyen dans son mouvement annuel. Or le mouvement annuel en question est rétrograde ; donc, dans un repère local, le soleil vrai prend du retard par rapport au soleil moyen, et l'équation du temps (composante bleue) est croissante. D'avril à septembre, la situation est inversée, le soleil vrai rattrape son retard et l'équation du temps décroît. Et la composante bleue de l'équation du temps rend bien compte que c'est début avril et début octobre que se trouvent ses extrema.

En première approximation, l'équation du temps varie sinusoïdalement avec une période d'une année, s'annule au périhélie et à l'aphélie, et est extrémal entre ces deux points (courbe bleue de la figure équation du temps). L'expression de ce retard dû à l'ellipticité, exprimé en minutes, est le suivant :

ΔT_C(d) = 7,53.cos(B) + 1,5.sin(B) = 7,678.sin(B + 1,374).

Voir la définition de B(d) ci-dessous.

Influence de l'obliquité de la Terre

Ce schéma montre le plan de l'écliptique et celui de l'équateur céleste qui interceptent chacun un grand cercle sur une sphère virtuelle appelée sphère céleste.
L'intersection de ces deux plans définit l'axe vernal.
L'angle entre ces 2 plans est appelé l'obliquité.

Illustration des étapes nécessaires pour qu'un méridien terrestre donné retourne face au Soleil d'un jour au suivant (point de vue géocentrique)

On suppose ici que l'orbite de la Terre est circulaire. Même dans ce cas, le mouvement apparent du Soleil le long de l'équateur céleste n'est pas uniforme, par suite de l'inclinaison de l'axe de rotation de la Terre par rapport à son plan orbital.

La figure ci-contre montre les trois étapes du retour d'un méridien face au Soleil, en adoptant un point de vue géocentrique, c'est-à-dire la Terre étant fixe au centre de la figure et le Soleil orbitant autour de la Terre :

• une rotation complète (360°) de la Terre sur elle-même pour passer de 1 à 2 ;
• ce faisant, le Soleil a avancé sur son orbite autour de la Terre, et de ce fait la Terre montre ce même méridien non pas face au Soleil mais face aux étoiles lointaines, point 2 ;
• une rotation complémentaire de la Terre sur elle-même est alors nécessaire pour que le méridien soit à nouveau face au Soleil, point 3.

Le Soleil a avancé de façon régulière sur son orbite située dans le plan écliptique, alors que la rotation complémentaire de la Terre sur elle-même pour se remettre face au Soleil est mesurée dans le plan de l'équateur céleste. Il faut donc rapporter le mouvement du Soleil dans ce plan de l'équateur céleste pour apprécier le retard ou l'avance du temps solaire par rapport à une horloge régulière. C'est cette opération, appelée réduction à l’équateur, qui explique que le mouvement apparent du Soleil le long de l'équateur céleste n'est pas uniforme.

L'orbite étant supposée circulaire, le module du vecteur vitesse du Soleil est donc constant le long de son orbite. Une composante de ce vecteur est portée par l'axe vernal, l'autre est portée par un vecteur orthogonal à cet axe vernal et situé dans le plan écliptique. La première composante se projette sans modification sur le plan de l'équateur céleste, la seconde se projette avec un facteur de réduction égal au cosinus de l'obliquité. De façon intuitive, la somme de ces deux projections sur le plan de l'équateur céleste sera minimale sur l'axe vernal et maximale sur la quadrature de ce même axe. La variation de vitesse sera donc nulle en ces quatre points, et de même pour l'avance ou le retard.

En première approximation, il s'agit d'une sinusoïde de période une demi-année (courbe verte de la figure Équation du temps, cf. plus haut) qui s'annule quatre fois sur une année, en particulier à l'équinoxe de printemps. L'expression de ce retard dû à l'obliquité, exprimé en minutes, est le suivant :

ΔT_R(d) = − 9,87.sin(2B).

Voir la définition de B(d) ci-dessous.

Version simplifiée de l'équation temps

La somme des deux formules précédentes fournit une première approximation de l’équation du temps :

ΔT(d) = ΔT_C(d) + ΔT_R(d),

c'est-à-dire :

ΔT(d) = 7,678.sin(B + 1,374) − 9,87sin(2B)

avec : $B(d) = \frac{2\pi(d-81)}{365}$, exprimé en radians, dépend du numéro du jour de l'année :

d = 1 le 1^er janvier ; d = 81 à l'équinoxe de printemps.

Étude détaillée

Influence de l'ellipticité de l'orbite de la Terre

• Calcul de l'anomalie moyenne $M(d)=357,5291\ ^\circ+0,98560028(J-J_{2000}),$

avec J₂₀₀₀ = 2451545. J est le jour julien de la date considérée. En première approximation, (J − J₂₀₀₀) peut être remplacé par le numéro d du jour dans l'année (d = 1 le premier janvier).

• Contribution de l'ellipticité de la trajectoire : c'est l'équation du centre en radian $C(d)=\left(2e-\frac{1}{4}e^3\right)\sin(M)+\frac{5}{4}e^2\sin(2M)+\frac{13}{12}e^3\sin(3M)$ où e = 0,01671 représente l'excentricité de la trajectoire de la Terre autour du Soleil.

Application numérique :

$C(d)=1,9148\ ^\circ\sin(M)+0,0200\ ^\circ\sin(2M)+0,0003\ ^\circ\sin(3M)\,$.

Influence de l'obliquité de la Terre

• Calcul de la longitude écliptique $\lambda_s=280,4665\ ^\circ+0,98564736(J-J_{2000})+C,$.

la petite différence de période entre λ_s et M est due à la précession des équinoxes.

• Contribution de l'obliquité de la Terre : c'est la réduction à l'équateur (en radian) $R(d)=(-y^2)\sin(2\lambda_s) +\left(\frac{1}{2}y^4\right)\sin(4\lambda_s)-\left(\frac{1}{3}y^6\right)\sin(6\lambda_s) \,$ où $y=\tan(\epsilon/2)$ ; $\epsilon=\rm 23,43929^o$ représente l'inclinaison de l'axe de la Terre par rapport au plan de l'écliptique.

Application numérique :

$R(d)= -2,4680\ ^\circ\sin(2\lambda_s) + 0,0530\ ^\circ\sin(4\lambda_s) -0,0014\ ^\circ\sin(6\lambda_s)\,$.

Équation du temps

Équation du temps en degrés :

E(d) = C(d) + R(d), où C et R sont exprimés en degrés.

Équation du temps en minutes :

$\Delta T(d) = 4 \times E(d)$, où E est exprimé en degrés.

Explications et démonstration de la formule

La figure 1 montre la Terre T qui tourne sur elle-même et qui tourne autour du Soleil S en un an dans le plan de l'écliptique. La situation présentée correspond à l'automne. Le point P est le périhélie, atteint au début du mois de janvier. L'angle θ s'appelle anomalie vraie. L'axe γ, appelé axe vernal ou point vernal, est l'intersection du plan de l'écliptique avec le plan équatorial. Il sert d'origine pour mesurer la longitude écliptique λ_s.

Figure 1

La figure 2 représente la Terre dans un repère fixe par rapport aux étoiles. L'obliquité ε est l'angle entre le plan de l'écliptique et le plan de l'équateur.

Figure 2

Appelons t le temps qui s'écoule. Considérons un point $A\left(t\right)$ fixé sur Terre et positionné sur l'équateur. Il fait donc un tour en un jour sidéral de façon régulière.

Partant du centre de la Terre, le point $S\left(t\right)$ est situé en direction du Soleil. Il se situe donc sur le cercle de l'écliptique. Le point $S\left(t\right)$ fait un tour en une année sidérale.

Comme l'orbite terrestre est elliptique et d'après les lois de Kepler, $S\left(t\right)$ ne tourne pas de façon régulière. Considérons le méridien passant par $S\left(t\right)$ et appelons $B\left(t\right)$ l'intersection de ce méridien avec l'équateur. Remarquons qu'il est midi solaire au point A lorsque le point $A\left(t\right)$ traverse ce méridien (par exemple lorsque les points $A\left(t\right)$

et $B\left(t\right)$ coïncident). Remarquons aussi qu'un jour solaire vrai est la durée qui sépare deux croisements de $A\left(t\right)$ et $B\left(t\right)$. Plus généralement l'heure solaire vraie est l'angle entre $A\left(t\right)$ et $B\left(t\right)$ :

$H_{V}\left(t\right)=\widehat{BA}=\widehat{B\gamma}+\widehat{\gamma A}$

H_V est l'heure qu'indiquerait un cadran solaire.

Pour définir l'heure solaire moyenne, il faut se référer à des mouvements réguliers (moyennés). Nous avons vu que le point $A\left(t\right)$ a un mouvement régulier. Ce n'est pas le cas du point $B\left(t\right)$, ni même du point $S\left(t\right)$. À la place de $S\left(t\right)$ on considère un point virtuel $\tilde{S}\left(t\right)$ sur l'écliptique qui a un mouvement régulier et de même période (on verra que $\tilde{S}\left(t\right)$ est directement relié à l'anomalie moyenne $M\left(t\right)$).

Par conséquent l'heure solaire moyenne est :

$H_{M}\left(t\right)=\widehat{\tilde{S}\gamma}+\widehat{\gamma A}$.

Par définition l'équation du temps est la différence :

$E\left(t\right)=H_{M}-H_{V}=\widehat{\tilde{S}\gamma}+\widehat{\gamma B}$.

Or une relation de trigonométrie donne :

$\tan\left(\widehat{\gamma B}\right)=\cos\varepsilon\tan\left(\widehat{\gamma S}\right)$.

En effet, projetons à partir du centre de la Terre le triangle sphérique SγB sur le plan tangent à la Terre au point vernal γ. Il devient un triangle rectangle d'angle ε au sommet γ et de côté adjacent $a=\tan(\widehat{\gamma B})$ et d'hypothénuse $h=\tan(\widehat{\gamma S})$. On déduit la relation a = hcos ε.

On déduit :

$\tan\left(\widehat{\gamma\tilde{S}}+E\right)=\cos\varepsilon\tan\left(\widehat{\gamma S}\right)$

et l'expression de l'équation du temps :

$E\left(t\right)=-\widehat{\gamma\tilde{S}}+\arctan\left(\cos\varepsilon\tan\left(\widehat{\gamma S}\right)\right)$.

Remarques

• Dans cette dernière équation tout est connu. D'une part, d'après la figure 1, il apparaît que l'angle $\widehat{\gamma S}$ est relié à la longitude écliptique λ_s par $\widehat{S\gamma}+\lambda_{s}=(\widehat{\overrightarrow{TS},\gamma})+(\widehat{\gamma,\overrightarrow{ST}})=\pi$

et λ_s elle-même est reliée à l'anomalie vraie $\theta\left(t\right)$ par $\lambda_{s}=\theta+\overline{\omega}$ où $\overline{\omega}$ est la longitude du périhélie. Donc

$\begin{matrix} \widehat{\gamma S} & = & \theta\left(t\right)+\overline{\omega}-\pi & \simeq & \theta\left(t\right)+102.94\ ^\circ \end{matrix}$.

De même $\widehat{\gamma\tilde{S}}$ est relié à l'anomalie moyenne $M\left(t\right)$ par

$\widehat{\gamma\tilde{S}}=M\left(t\right)+\overline{\omega}-\pi$.

• Traditionnellement on décompose E(t) de la façon suivante : $E\left(t\right)=\underbrace{\left(\widehat{\gamma S}-\widehat{\gamma\tilde{S}}\right)}_{C}+\underbrace{\left(-\widehat{\gamma S}+\arctan\left(\cos\varepsilon\tan\left(\widehat{\gamma S}\right)\right)\right)}_{R}$.

Le premier terme, C, est appelé « contribution de l'ellipticité » ou équation du centre. On a : $C=\left(\widehat{\gamma S}-\widehat{\gamma\tilde{S}}\right)=\theta\left(t\right)-M\left(t\right)$. C est dû à l'ellipticité de l'orbite terrestre. Dans un modèle où la Terre aurait un mouvement circulaire et régulier, on aurait $\tilde{S}=S$ et seul le deuxième terme, R, appelé réduction à l'équateur et dû à l'obliquité ε, interviendrait. Remarquer que si on avait ε = 0 (le Soleil en permanence dans le plan de l'équateur), ce dernier terme serait nul.

• À l'aide d'un développement limité, et en posant $y=\tan\frac{\varepsilon}{2}\simeq0.20$, le terme de la réduction à l'équateur ci-dessus peut s'écrire

$\begin{align} R & = & \arctan\left(\cos\varepsilon\tan\left(\widehat{\gamma S}\right)\right)-\widehat{\gamma S} \\ & = & \arctan\left(\cos\varepsilon\tan\left(\lambda_{s}-\pi\right)\right)-\left(\lambda_{s}-\pi\right) \\ & = & -\arctan\left(\frac{\sin\left(2\lambda_{s}\right)y^{2}}{1+y^{2}\cos\left(2\lambda_{s}\right)}\right) \\ & \simeq & -y^{2}\sin\left(2\lambda_{s}\right)+\frac{1}{2}y^{4}\sin\left(4\lambda_{s}\right)-\frac{1}{3}y^{6}\sin\left(6\lambda_{s}\right)+\ldots \end{align}$.

Notes et références

1. ↑ Dans les Éphémérides Astronomiques publiées par la Société astronomique de France, la valeur de l'équation du temps est donnée, pour chaque jour de l'année à 0 h de Temps universel.
2. ↑ Voir la définition donnée par l'« Institut de mécanique céleste et de calcul des éphémérides » (Observatoire de Paris - Bureau des longitudes - CNRS) Temps vrai, temps moyen, équation du temps.
3. ↑ Voir la définition donnée par Laplace dans son livre Exposition du système du monde - Livre Premier, chapitre 3, fin du §3.
4. ↑ Schéma établi en utilisant les données de hauteur et d'azimut fournies par le site JPL Horizons
5. ↑ Un diagramme similaire se trouve sur le site [1]
6. ↑ L'analemme d'Anaximandre à Ptolémée
7. ↑ Source : J. Meeus et D. Savoie (cf. Bibliographie).

Voir aussi

Articles connexes

• Anomalie excentrique
• Anomalie vraie
• Anomalie moyenne
• Analemme
• Cadran solaire

Liens externes

Sur les autres projets Wikimedia :

• « Équation du temps », sur Wikimedia Commons (ressources multimédia)

• Les échelles de temps par l'Institut de mécanique céleste et de calcul des éphémérides
• L'équation du temps de Kepler
• (en) Article très complet sur l'équation du temps

Bibliographie

• Pierre-Simon Laplace, Exposition du système du monde
• Jean Meeus et Denis Savoie, L’Équation du temps, in revue L'Astronomie, vol. 109, juin 1995, pp. 188-193.