Position du Soleil

Ces équations, tirées de l'Almanach astronomique[3] - [4], peuvent être utilisées pour calculer les coordonnées apparentes du Soleil, l'équinoxe moyen et l'écliptique de la date, avec une précision d'environ 0° 0,01' (soit 36 secondes d'arc), pour les dates comprises entre 1950 et 2050.

On commence par calculer n, le nombre de jours (positifs ou négatifs) depuis midi de Greenwich, heure terrestre, le 1^er janvier 2000 (J2000.0). Si on connaît la date julienne pour l'heure souhaitée,

${\displaystyle n=\mathrm {DJ} -2\,451\,545,0}$

La longitude moyenne du Soleil, corrigée de l'aberration de la lumière, est ce nombre de jours fois le déplacement diurne moyen, ajouté à la position du périhélie (longitude écliptique mesurée depuis le point vernal) :

${\displaystyle L=280,460^{\circ }+0,9856474^{\circ }n}$

L'anomalie moyenne du Soleil (en fait, de la Terre dans son orbite autour du Soleil, mais il est commode de prétendre que le Soleil tourne autour de la Terre), est :

${\displaystyle g=357,528^{\circ }+0,9856003^{\circ }n}$

Mettre ${\displaystyle L}$ et ${\displaystyle g}$ dans la plage de 0° à 360° en ajoutant ou en soustrayant des multiples de 360° selon les besoins.

Enfin, la longitude écliptique du Soleil est :

${\displaystyle \lambda =L+1,915^{\circ }\sin g+0,020^{\circ }\sin 2g}$

La latitude écliptique du Soleil est presque :

${\displaystyle \beta =0}$,

comme la latitude écliptique du Soleil ne dépasse jamais 0,00033°[5],

et la distance du Soleil à la Terre, en unités astronomiques, est :

${\displaystyle R=1,00014-0,01671\cos g-0,00014\cos 2g}$.

${\displaystyle \lambda }$, ${\displaystyle \beta }$ et ${\displaystyle R}$ caractérisent complètement la position du Soleil dans le système de coordonnées écliptiques. Ils peuvent être convertis dans le système de coordonnées équatoriales en calculant :

• l'obliquité de l'écliptique, ${\displaystyle \epsilon }$ ;
• l'ascension droite ${\displaystyle \alpha =\arctan(\cos \epsilon \tan \lambda )}$, où ${\displaystyle \alpha }$ est dans le même quadrant que ${\displaystyle \lambda }$. Pour obtenir l'ascension droite dans le bon quadrant , on utilise la fonction arc tangente à double argument (atan2(y,x) dans de nombreux langages informatiques) : ${\displaystyle \alpha =\arctan 2(\cos \epsilon \sin \lambda ,\cos \lambda )}$ ;
• et la déclinaison, ${\displaystyle \delta =\arcsin(\sin \epsilon \sin \lambda )}$.

En coordonnées équatoriales rectangulaires à droite (où le ${\displaystyle X}$ l'axe est dans la direction du point vernal, et le ${\displaystyle Y}$ l'axe est de 90° à l'est, dans le plan de l'équateur céleste, et le ${\displaystyle Z}$ l'axe est dirigé vers le pôle nord céleste[6]), en unités astronomiques :

${\displaystyle X=R\cos \lambda }$

${\displaystyle Y=R\cos \epsilon \sin \lambda }$

${\displaystyle Z=R\sin \epsilon \sin \lambda }$

Lorsque l'obliquité de l'écliptique n'est pas obtenue ailleurs, elle peut être approximée :

${\displaystyle \epsilon =23,439^{\circ }-0,0000004^{\circ }n}$

à utiliser avec les équations ci-dessus.

Le Soleil semble se déplacer vers le nord pendant le printemps nordique, en contact avec l'équateur céleste lors de l'équinoxe de mars. Sa déclinaison atteint un maximum égal à l'angle d'inclinaison axiale de la Terre (23,44°)[7] - [8] au solstice de juin, puis diminue jusqu'à atteindre son minimum (−23,44°) au solstice de décembre, lorsque sa valeur est le négatif de l'inclinaison axiale. Cette variation produit les saisons.

Un graphique linéaire de la déclinaison du Soleil au cours d'une année ressemble à une onde sinusoïdale d'une amplitude de 23,44°, mais un lobe de la vague est plusieurs jours plus long que l'autre, entre autres différences.

Les phénomènes suivants se produiraient si la Terre était une sphère parfaite, dans une orbite circulaire autour du Soleil, et si son axe était incliné à 90°, de sorte que l'axe lui-même est sur le plan orbital (similaire à Uranus). À une date de l'année, le Soleil serait directement au-dessus du pôle Nord, donc sa déclinaison serait de +90°. Au cours des prochains mois, le point subsolaire se déplacerait vers le pôle Sud à vitesse constante, traversant les cercles de latitude à un rythme constant, de sorte que la déclinaison solaire diminuerait linéairement avec le temps. Finalement, le Soleil serait directement au-dessus du pôle Sud, avec une déclinaison de –90° ; alors il commencerait à se déplacer vers le nord à une vitesse constante. Ainsi, le graphique de la déclinaison solaire, vu de cette Terre très inclinée, ressemblerait à une onde triangulaire plutôt qu'à une onde sinusoïdale, zigzaguant entre plus et moins 90°, avec des segments linéaires entre les maxima et les minima.

Si l'inclinaison axiale à 90° est diminuée, les valeurs absolues maximale et minimale de la déclinaison diminueraient, pour égaler l'inclinaison axiale. De plus, les formes des maxima et des minima sur le graphique deviendraient moins aiguës (« pointues »), étant incurvées pour ressembler aux maxima et aux minima d'une onde sinusoïdale. Cependant, même lorsque l'inclinaison axiale est égale à celle de la Terre réelle, les maxima et les minima restent plus aigus que ceux d'une onde sinusoïdale.

En réalité, l'orbite de la Terre est elliptique. La Terre se déplace plus rapidement autour du Soleil près du périhélie, début janvier, que près de l'aphélie, début juillet. Cela accélère les processus comme la variation de la déclinaison solaire en janvier qu'en juillet. Sur le graphique, cela rend les minima plus aigus que les maxima. De plus, comme le périhélie et l'aphélie ne se produisent pas aux dates exactes comme les solstices, les maxima et les minima sont légèrement asymétriques. Les taux de variation avant et après ne sont pas tout à fait égaux.

Le graphique de la déclinaison solaire apparente est donc différent à plusieurs égards d'une onde sinusoïdale. Le calculer avec précision implique une certaine complexité, comme indiqué ci-dessous.

Diagramme de la position du soleil dans le ciel.

La déclinaison du Soleil, δ _☉, est l'angle entre les rayons du Soleil et le plan de l'équateur terrestre. L'inclinaison axiale de la Terre (appelée obliquité de l'écliptique par les astronomes) est l'angle entre l'axe de la Terre et une ligne perpendiculaire à l'orbite de la Terre. L'inclinaison axiale de la Terre change lentement au cours des milliers d'années, mais sa valeur actuelle d'environ ε = 23° 26' est presque constante, de sorte que le changement de déclinaison solaire au cours d'une année est presque le même que l'année suivante.

Aux solstices, l'angle entre les rayons du Soleil et le plan de l'équateur terrestre atteint sa valeur maximale de 23° 26'. Par conséquent, δ _☉ = + 23° 26' au solstice d'été nord et δ _☉ = −23° 26' au solstice d'été sud.

Au moment de chaque équinoxe, le centre du Soleil semble passer par l'équateur céleste, et δ_☉ vaut 0°.

La déclinaison du Soleil à un moment donné est calculée par :

${\displaystyle \delta _{\odot }=\arcsin \left[\sin \left(-23,44^{\circ }\right)\cdot \sin \left(LE\right)\right]}$

où LÉ est la longitude écliptique (essentiellement, la position de la Terre sur son orbite). Étant donné que l'excentricité orbitale de la Terre est petite, son orbite peut être approximée comme un cercle qui provoque jusqu'à 1° d'erreur. L'approximation du cercle signifie que l'EL serait à 90° devant les solstices de l'orbite terrestre (aux équinoxes), de sorte que sin (EL) peut s'écrire sin (90 + NDS) = cos (NDS) où NDS est le nombre de jours après le solstice de décembre. En utilisant également l'approximation selon laquelle arcsin [sin (d) · cos (NDS)] est proche de d · cos (NDS), la formule fréquemment utilisée suivante est obtenue :

${\displaystyle \delta _{\odot }=-23,44^{\circ }\cdot \cos \left[{\frac {360^{\circ }}{365}}\cdot \left(N+10\right)\right]}$

où N est le jour de l'année commençant par N = 0 à minuit, heure universelle (UT) au début du 1^er janvier (c'est-à-dire la partie jours de la date ordinale -1). Le nombre 10, en (N + 10), est le nombre approximatif de jours après le solstice de décembre jusqu'au 1^er janvier. Cette équation surestime la déclinaison proche de l'équinoxe de septembre jusqu'à + 1,5°. L'approximation de la fonction sinusoïdale conduit en soi à une erreur allant jusqu'à 0,26° et a été déconseillée pour une utilisation dans les applications d'énergie solaire[2]. La formule de Spencer de 1971[9] (basée sur une série de Fourier) est également déconseillée pour avoir une erreur allant jusqu'à 0,28°[10]. Une erreur supplémentaire allant jusqu'à 0,5° peut se produire dans toutes les équations autour des équinoxes si on n'utilise pas de décimale lorsqu'on sélectionne N pour ajuster le temps après UT minuit pour le début de la journée. Ainsi, l'équation ci-dessus peut avoir jusqu'à 2,0° d'erreur, environ quatre fois la largeur angulaire du soleil, selon la façon dont elle est utilisée.

La déclinaison peut être calculée avec plus de précision en ne faisant pas les deux approximations, en utilisant les paramètres de l'orbite de la Terre pour estimer plus précisément EL[11] :

${\displaystyle \delta _{\odot }=\arcsin \left[\sin \left(-23,44^{\circ }\right)\cdot \cos \left({\frac {360^{\circ }}{365,24}}\left(N+10\right)+{\frac {360^{\circ }}{\pi }}\cdot 0,0167\sin \left({\frac {360^{\circ }}{365,24}}\left(N-2\right)\right)\right)\right]}$

qui peut être simplifié en évaluant les constantes pour :

${\displaystyle \delta _{\odot }=-\arcsin \left[0,39779\cos \left(0,98565^{\circ }\left(N+10\right)+1,914^{\circ }\sin \left(0,98565^{\circ }\left(N-2\right)\right)\right)\right]}$

N est le nombre de jours depuis minuit (temps universel) au début du 1^er janvier (c'est-à-dire la partie jours de la date ordinale -1) et peut inclure des décimales pour s'adapter aux heures locales plus tard ou plus tôt dans la journée. Le nombre 2, en (N-2), est le nombre approximatif de jours après le 1^er janvier jusqu'au périhélie de la Terre. Le nombre 0,0167 est la valeur actuelle de l'excentricité de l'orbite terrestre. L'excentricité varie très lentement dans le temps, mais pour des dates assez proches du présent, elle peut être considérée comme constante. Les plus grandes erreurs dans cette équation sont inférieures à ±0,2°, mais inférieures à ±0,03° pour une année donnée si le nombre 10 est ajusté à la hausse ou à la baisse en jours fractionnaires, comme déterminé par la distance du solstice de décembre de l'année précédente avant ou après midi le 22 décembre. Ces précisions sont comparées aux calculs avancés de la NOAA[12] - [13] qui sont basés sur l'algorithme Jean Meeus de 1999 qui est précis à 0,01° près[14].

(La formule ci-dessus est liée à un calcul raisonnablement simple et précis de l'équation du temps, qui est décrit ici.)

Des algorithmes plus complexes[15] - [16] corrigent les modifications de la longitude écliptique en utilisant des termes en plus de la correction d'excentricité de premier ordre ci-dessus. Ils corrigent également l'obliquité de 23,44° qui change très légèrement avec le temps. Les corrections peuvent également inclure les effets de la lune en décalant la position de la Terre par rapport au centre de l'orbite de la paire autour du Soleil. Après avoir obtenu la déclinaison par rapport au centre de la Terre, une autre correction de parallaxe est appliquée, qui dépend de la distance de l'observateur par rapport au centre de la Terre. Cette correction est inférieure à 0,0025°. L'erreur dans le calcul de la position du centre du Soleil peut être inférieure à 0,00015°. À titre de comparaison, la largeur du Soleil est d'environ 0,5°.

Les calculs de déclinaison décrits ci-dessus n'incluent pas les effets de la réfraction de la lumière dans l'atmosphère, ce qui fait que l'angle d'élévation apparent du Soleil vu par un observateur est plus élevé que l'angle d'élévation réel, en particulier à faible élévation du Soleil[2]. Par exemple, lorsque le Soleil est à une altitude de 10°, il semble être à 10,1°. La déclinaison du Soleil peut être utilisée, avec son ascension droite, pour calculer son azimut et aussi sa véritable élévation, qui peut ensuite être corrigée pour la réfraction pour donner sa position apparente[13] - [17].

• Solar Position Algorithm, sur le site du Renewable Resource Data Center du National Renewable Energy Laboratory.
• Sun Position Calculator, sur pveducation.org. Une calculatrice interactive montrant la trajectoire du soleil dans le ciel.
• NOAA Solar Calculator, sur le site de la division de surveillance mondiale du laboratoire de recherche sur les systèmes terrestres de la NOAA.
• Calculateur de déclinaison et de position du soleil de la NOAA
• Système HORIZONS, sur le site du JPL. Positions très précises des objets du système solaire basées sur les éphémérides de la série JPL DE.
• Éphémérides générales des corps du système solaire, sur le site de l' IMCCE. Positions d'objets du système solaire basées sur les éphémérides de la série INPOP.
• Position solaire dans le package R. Insol.