Inhaltsverzeichnis
Der Erdmond
Mondposition
In diesem Kapitel werden die geozentrisch ekliptikalen Koordinaten des Mondes berechnet und dabei zwei Theorien gegenüber gestellt:
- ILE (Brownsche Mondtheorie), entnommen aus Astronomie mit dem Personal Computer (O. Montenbruck & T. Pfleger)
- ELP2000 (Chaprontsche Theorie), entnommen aus Astronomical Algorithms (J. Meeus)
in niedriger, aber für die meisten Zwecke ausreichender Genauigkeit. Es gilt für $T$ die Epoche $J1900.0$.
$$\label{epoch} T = \frac{JDE - 2415020.0}{36525}\tag{1}$$
Mittlere Bahnelemente
| Tabelle 1 | |
|---|---|
| Sonne: | |
| $L=$ | $279\overset{\circ}{.}696678 +36000\overset{\circ}{.}0 \cdot T$ |
| $+0\overset{\circ}{.}768925 \cdot T +3\overset{\circ}{.}03 \cdot 10^{-4} \cdot T^2$ | |
| $M=$ | $358\overset{\circ}{.}475833 +35640\overset{\circ}{.}0 \cdot T$ |
| $+359\overset{\circ}{.}049750 \cdot T -1\overset{\circ}{.}5 \cdot 10^{-4} \cdot T^2$ | |
| Mond: | |
| $l=$ | $270\overset{\circ}{.}434164 +480960\overset{\circ}{.}0 \cdot T$ |
| $+307\overset{\circ}{.}883142 \cdot T -1\overset{\circ}{.}133 \cdot10^{-3} \cdot T^2$ | |
| $m=$ | $296\overset{\circ}{.}104608 +477000\overset{\circ}{.}0 \cdot T$ |
| $+198\overset{\circ}{.}849108 \cdot T +9\overset{\circ}{.}192 \cdot10^{-3} \cdot T^2$ | |
| $\Omega=$ | $259\overset{\circ}{.}183275 -1800\overset{\circ}{.}0 \cdot T$ |
| $-134\overset{\circ}{.}142008 \cdot T +2\overset{\circ}{.}078 \cdot10^{-3} \cdot T^2$ | |
| Delaunay Argumente: | |
| $F=l-\Omega=$ | $11\overset{\circ}{.}250889 +483120\overset{\circ}{.}0 \cdot T$ |
| $+82\overset{\circ}{.}0251 \cdot T -3\overset{\circ}{.}211 \cdot 10^{-3} \cdot T^2$ | |
| $D=l-L=$ | $350\overset{\circ}{.}737486 +444960\overset{\circ}{.}0 \cdot T$ |
| $+307\overset{\circ}{.}1142 \cdot T -1\overset{\circ}{.}436 \cdot 10^{-3} \cdot T^2$ | |
$L$ = mittlere Länge der Sonne
$M$ = mittlere Anomalie der Sonne
$l$ = mittlere Länge des Mondes
$m$ = mittlere Anomalie des Mondes
$\Omega$ = Länge des aufsteigenden Knotens des Mondes
$F$ = Argument der Breite
$D$ = Argument der Länge
Störungsterme
In der Länge $\lambda$ fallen beide Theorien in den Störtermen zusammen:
\[\begin{align} \Delta\lambda = & +22640'' \cdot \sin(m)\\ &+769'' \cdot \sin(2 m) \\ &+36'' \cdot \sin(3 m) \\ &+4587'' \cdot \sin(2 D - m)\\ &+2370'' \cdot \sin(2 D) \\ &-668'' \cdot \sin(M) \\ &-412'' \cdot \sin(2 F) \\ &+212'' \cdot \sin(2 \cdot (D - m)) \\ &+206'' \cdot \sin(2 D - M - m)\\ &+192'' \cdot \sin(2 D + m)\\ &+165'' \cdot \sin(2 D - M)\\ &+148'' \cdot \sin(m - M) \\ &-125'' \cdot \sin(D) \\ &-110'' \cdot \sin(m + M)\\ &+55'' \cdot \sin(2 \cdot (D - F))\\ &-45'' \cdot \sin(2 F + m)\\ &-40'' \cdot \sin(2 F - m)\\ &+38'' \cdot \sin(4 D - m) \end{align}\tag{2}\]
In der Breite $\beta$ gehört $\Delta\beta_1$ zur ILE und $\Delta\beta_2$ zur ELP2000. Die Brownsche Mondtheorie benötigt vorab noch einen Hilfsterm:
\[\begin{align} \Delta G =&\; \frac{1}{3600''}\cdot \big[\Delta\lambda\\ &+412''\cdot \sin(2\cdot F)\\ &+541''\cdot \sin(M)\big] \end{align}\tag{3}\]
Störterme in der Breite für die DE200: \[\begin{align} \Delta\beta_1 = & +18520'' \cdot \sin(F + \Delta G)\\ &-526'' \cdot \sin(F - 2 D)\\ &+44'' \cdot \sin(F - 2 D + m)\\ &-31'' \cdot \sin(F - 2 D - m)\\ &-23'' \cdot \sin(F - 2 D + M)\\ &+11'' \cdot \sin(F - 2 D - M)\\ &-25'' \cdot \sin(F - 2 m)\\ &+21'' \cdot \sin(F - m)\\ &+24'' \cdot \sin(F + M)\\ &-14'' \cdot \sin(m) \end{align}\tag{4}\]
Störterme in der Breite für die ELP2000: \[\begin{align} \Delta\beta_2 = & +18462'' \cdot \sin(F)\\ &+1010'' \cdot \sin(m + F)\\ &+1000'' \cdot \sin(m - F)\\ &+624'' \cdot \sin(2 D - F)\\ &+200'' \cdot \sin(2 D + F - m)\\ &+167'' \cdot \sin(2 D - F - m)\\ &+117'' \cdot \sin(2 D + F)\\ &+62'' \cdot \sin(2 m + F)\\ &+33'' \cdot \sin(2 D - F + m )\\ &+32'' \cdot \sin(2 m - F) \end{align}\tag{5}\]
Auch bei der Parallaxe fallen beide Theorien in dieser Größenordnung zusammen: \[\begin{align} \Delta\Pi = & +187'' \cdot \cos(m)\\ &+10'' \cdot \cos(2 m)\\ &+34'' \cdot \cos(2 D - m)\\ &+28'' \cdot\cos(2 D)\\ &+3'' \cdot \cos(2 D + m) \end{align}\tag{6}\]
Die Parallaxe ist der vom Mond aus gesehene Erdhalbmesser.
Es fehlt nur noch der Störungsterm für den Radius $r$:
\[\begin{align} \Delta r = & -20905\,\mathrm{km} \cdot\cos(m\\ &-570\,\mathrm{km} \cdot \cos(2 m)\\ &-3699\,\mathrm{km} \cdot \cos(2 D - m)\\ &-2956\,\mathrm{km} \cdot\cos(2 D)\\ &+246\,\mathrm{km} \cdot \cos(2 (m - D))\\ &-205\,\mathrm{km} \cdot \cos(M - 2 D)\\ &-171\,\mathrm{km} \cdot\cos(m + 2 D)\\ &-152\,\mathrm{km} \cdot \cos(m + M - 2 D) \end{align}\tag{7}\]
Auch beim Radius $r$ fallen beide Theorien zusammen. Damit hat man die geozentrisch ekliptikalen Koordinaten des Mondes vorliegen:
| Tabelle 2 | |
|---|---|
| Länge: | $\large\lambda = l + \frac{\Delta\lambda}{3600''}$ |
| Breite: | $\large\beta_{k}=\frac{\Delta\beta_k}{3600''}$ |
| Radiusvektor: | $\large\Delta = 385000 \ km + \Delta r$ |
| Parallaxe: | $\large\Pi = 0\overset{\circ}{.}95333 + \frac{\Delta\Pi}{3600''}$ |
Der Index $k = 1, 2$ bezeichnet die jeweilige Theorie.
Beispiel
Man berechne die geozentrischen ekliptikalen Koordinaten des Mondes für den 15.4.2023 um 22:15 mitteleuropäische Sommerzeit (MESZ) mit dem oben vorgestellen „schnellen“ Algorithmus!
Im Beispiel wurden alle Werte mittels JavaScript ausgewertet und alle Kommastellen stehen gelassen, um dem Leser beim Nachvollziehen der Berechnung diese Werte anzugeben. Es macht natürlich keinen Sinn, so viele Kommastellen anzugeben! Werte für Winkel werden in der Regel auf 6 Kommastellen gerundet. Für $JD$ genügen 5 Kommastellen, da $0\overset{d}{.}00001 = 0\overset{s}{.}864$ sind. Für die Größe von $T$ sollten jedoch alle verfügbaren Nachkommastellen mitgezogen werden!
Für den gegebenen Zeitpunkt wurde der Julianische Tag hier bereits ermittelt zu $JD=2460050.34375$. Im Jahr 2023 war der Wert von $\Delta T = 69^{s}$, diese müssen hinzugefügt werden, um die Position des Mondes in der gleichmäßigen Skala der dynamischen Zeit $TD$ zu erhalten. Ein Tag hat 86400 Sekunden, daher folgt
\(\begin{align} JDE &= 2460050.34375 + \frac{69^{s}}{86400\frac{s}{d}}\\ &= 2460050.344548611\\ &= 2460050.34455 \end{align}\)
Die julianischen Jahrhunderte werden in diesem Algorithmus bezüglich der Epoche $J1900$ berechnet!
\(\begin{align} T &= \frac{(2460050.34455 - 2415020.0)}{36525}\\ &= 1.2328636427104733 \end{align}\)
Für die Hauptwinkel $L$, $M$, $l$, $m$, $\Omega$, $F$ und $D$ erhält man hier mithilfe der Reduktions-Funktion:
\(\begin{align} L &= 44663\overset{\circ}{.}736255799195 = 23\overset{\circ}{.}736256\\ M &= 44740\overset{\circ}{.}39521390763 = 100\overset{\circ}{.}395214\\ l &= 593608\overset{\circ}{.}107971898 = 328\overset{\circ}{.}107972\\ m &= 588617\overset{\circ}{.}2299879402 = 17\overset{\circ}{.}229988\\ \Omega &= -2125\overset{\circ}{.}3469280403915 = 34\overset{\circ}{.}653072\\ F &= 595733\overset{\circ}{.}4548382952 = 293\overset{\circ}{.}454838\\ D &= 549304\overset{\circ}{.}3716951403 = 304\overset{\circ}{.}371695 \end{align}\)
Die Korrekturterme für $\Delta\lambda$ und $\Delta r$ sind für beide Theorien gleich zu berechnen, man erhält
$\Delta\lambda = +891\overset{''}{.}3721325847027$
$\Delta r = -16979.93271594326\;\textrm{km}$
Für die Korrekturterme in der Breite $\Delta\beta$ erhält man je nach der entsprechenden Theorie
$\Delta G = 0\overset{\circ}{.}31183834919807063$
$\Delta \beta_1 = -17304\overset{''}{.}290756891038$
bzw.
$\Delta \beta_2 = -17264\overset{''}{.}092183245077$
Für die Horizontalparallaxe ergibt sich die Korrektur
$\Delta\Pi = +155\overset{''}{.}33209544205963$
Die geozentrischen ekliptikalen Koordinaten $\lambda, \beta, \Delta$ des Mondes sind dann
\( \begin{align} \lambda &= 328\overset{\circ}{.}107972 + \frac{891\overset{''}{.}37}{3600\frac{''}{\circ}}\\ &= 328\overset{\circ}{.}355575\\ \beta_1 &= \frac{-17304\overset{''}{.}29}{3600\frac{''}{\circ}} = -4\overset{\circ}{.}806747\\ \beta_2 &= \frac{-17264\overset{''}{.}09}{3600\frac{''}{\circ}} = -4\overset{\circ}{.}795581\\ \Delta &= 385000 \ km + (-16979.9)\\ &= 368020.1\,\textrm{km} \end{align} \)
Schließlich erhält man noch die Horizontalparallaxe $\Pi$ des Mondes mit
\(\begin{align} \Pi &= 0\overset{\circ}{.}95333\\ &+ \frac{155\overset{\circ}{.}332095}{3600\frac{''}{\circ}}\\ &= 0\overset{\circ}{.}996478 \end{align}\)
Man vergleiche diese Daten mit dem Beispiel unter Mondposition nach Meeus:
| Tabelle Beispiel | ||
|---|---|---|
| Schnell | Meeus | Differenz |
| $\lambda = 328\overset{\circ}{.}355575$ | $\lambda = 328\overset{\circ}{.}387192$ | $-113\overset{''}{.}9$ |
| $\beta_1 = −4\overset{\circ}{.}806747$ | $\beta = −4\overset{\circ}{.}806013$ | $-2\overset{''}{.}6 $ |
| $\beta_2 = −4\overset{\circ}{.}795581$ | $+37\overset{''}{.}6 $ | |
| $\Delta = 368020.1\;\textrm{km}$ | $\Delta = 367995.8\;\textrm{km}$ | $+24.3\;\textrm{km} $ |
| $\Pi = 0\overset{\circ}{.}996478$ | $\pi = 0\overset{\circ}{.}9931058$ | $+12\overset{''}{.}1 $ |
Ausführliche Mondtheorien nach J. Meeus, O. Montenbruck und T. Pfleger
Die beiden ausführlichen Theorien mit deutlich höherer Genauigkeit haben aufgrund ihres großen Umfangs zwei Extraseiten erhalten. Sie sind hier und hier zu finden. Ferner gibt es ein Vergleich der beiden Theorien.
Physische Ephemeriden
Durchmesser
Der geozentrische Durchmesser des Mondes lautet:
\[\begin{align} &\sin\left(\frac{\varnothing}{2}\right) = 0.2724934056\cdot \sin(\Pi)\\ &\textrm{mit}\quad\dfrac{R_M}{R_E} = 0.2724934056 \end{align}\tag{8}\]
Und den topozentrischen Durchmesser ermittelt man mit: $$\tan\left(\frac{\varnothing}{2}\right) = \frac{R_M}{\Delta'}\tag{9}$$
$\Delta'$ ist der topozentrische Abstand. Die Ausgabe ist in Grad und sollte ggf. durch Multiplikation mit $60$ in Bogenminuten umgerechnet werden.
Elongation, Phasen- und Sonnenpositionswinkel
Der Sonnenpositionswinkel $\theta$ und die Elongation $\eta$ werden in diesem Abschnitt berechnet. $\alpha_0$, $\delta_0$ müssen die Koordinaten der Sonne und $\alpha$, $\delta$ die Koordinaten des Mondes sein. Ist die Elongation bekannt, so kann mit guter Genauigkeit (J. Meeus) die Relation für den Phasenwinkel $\varphi$ gelten: $$\cos(\varphi) = -\cos(\eta)\tag{10}$$
Dies ist möglich, weil $\sigma \ll 1$ gilt. Der Beleuchtungsgrad $k$ ist dann:
$$k = \frac{1}{2} \cdot \big(1 + \cos (\varphi)\big) \cdot 100\%\tag{11}$$
Es ist üblich, den Beleuchtungsgrad $k$ auf 2 Nachkommastellen genau anzugeben, z.B. $k = 0.69$. Astronomiesoftware hingegen gibt den Beleuchtungsgrad oft in Prozentwerten an.
Der genäherte Phasenwinkel
Nach J. Meeus Astronomical Algorithms kann man den Phasenwinkel $\varphi$ genähert mit der folgenden Beziehung berechnen:
Man berechnet die Delauney-Argumente $D, M$ und $m$ mit
\[\begin{align} D &=\,297\overset{\circ}{.}8501921\\ &+445267\overset{\circ}{.}1114034\cdot T\\ &-0\overset{\circ}{.}0018819\cdot T^2 \end{align}\tag{12}\]
\[\begin{align} M &=\,357\overset{\circ}{.}5291092\\ &+35999\overset{\circ}{.}0502909\cdot T\\ &-0\overset{\circ}{.}0001536\cdot T^2 \end{align}\tag{13}\]
\[\begin{align} m &=\,134\overset{\circ}{.}9633964\\ &+477198\overset{\circ}{.}8675055\cdot T\\ &+0\overset{\circ}{.}0087414\cdot T^2 \end{align}\tag{14}\]
Es genügt hier, nur Terme bis zur 2. Potenz in $T$ zu berücksichtigen. $T$ sind wieder die julianischen Jahrhunderte bezüglich der Epoche $J2000$, nämlich
$$T = \frac{JD - 2451545.0}{36525}\tag{15}$$
Daraus erhält man einen guten Näherungswert für den Phasenwinkel mit \[\begin{align} \varphi &=\,180^{\circ} - D\\ &-6\overset{\circ}{.}289\cdot\sin(m)\\ &+2\overset{\circ}{.}100\cdot\sin(M)\\ &-1\overset{\circ}{.}274\cdot\sin(2\cdot D - m)\\ &-0\overset{\circ}{.}658\cdot\sin(2\cdot D)\\ &-0\overset{\circ}{.}214\cdot\sin(2\cdot m)\\ &-0\overset{\circ}{.}110\cdot\sin(D) \end{align}\tag{16}\]
Beispiel
Man berechne den Phasenwinkel $\varphi$ und Beleuchtungsgrad $k$ des Mondes am 15.4.2023 um 22:15 mitteleuropäische Sommerzeit (MESZ)
Für den gegebenen Zeitpunkt wurde der Julianische Tag hier bereits ermittelt zu $JD = 2460050.34375$.
Damit folgt
\( \begin{align} T &= \frac{2460050.34375 - 2451545.0}{36525}\\ &= 0.23286362080767\\ D &= 103984\overset{\circ}{.}36187462\\ &= 304\overset{\circ}{.}36187462\\ M &= 8740\overset{\circ}{.}3982969718\\ &= 100\overset{\circ}{.}3982969718\\ m &= 111257\overset{\circ}{.}21999958\\ &= 17\overset{\circ}{.}21999958\\ \end{align} \)
Wie immer wurden große Winkel mithilfe der Reduktionsfunktion auf das Intervall [0°-360°] gebracht.
Der Phasenwinkel des Mondes beträgt daher
\( \begin{align} \varphi &=\,180^{\circ} - 304\overset{\circ}{.}362\\ &-6\overset{\circ}{.}289\cdot\sin(17\overset{\circ}{.}22)\\ &+2\overset{\circ}{.}100\cdot\sin(100\overset{\circ}{.}398)\\ &-1\overset{\circ}{.}274\cdot\sin(2\cdot 304\overset{\circ}{.}362 - 17\overset{\circ}{.}22)\\ &-0\overset{\circ}{.}658\cdot\sin(2\cdot 304\overset{\circ}{.}362)\\ &-0\overset{\circ}{.}214\cdot\sin(2\cdot 17\overset{\circ}{.}22)\\ &-0\overset{\circ}{.}110\cdot\sin(304\overset{\circ}{.}362)\\ &=-122\overset{\circ}{.}358 \end{align} \)
Der Beleuchtungsgrad $k$ beträgt damit
$k = \frac{1+cos(122.35)}{2}\cdot 100\% = 23.24\%$
Der Cosinus ist eine gerade Funktion, daher ist der Wert $\cos(-\varphi) = \cos(\varphi)$.
Die Astronomiesoftware GUIDE berechnet die exakten Werte zu $\varphi = 122\overset{\circ}{.}45$ und $k = 23.17\%$.
Helligkeit
Mit Hilfe des Phasenwinkels $\varphi$ ist es nun möglich, die Helligkeit $H$ des Mondes am Firmament zu berechnen:
\[\begin{align} H = -12\overset{m}{.}74 &+ 5\cdot \log_{10}(R\cdot\Delta) \\ &+ 2\overset{m}{.}825\cdot \left(\frac{\varphi}{100}\right) \\ &-0\overset{m}{.}51\cdot \left(\frac{\varphi}{100}\right)^2 \\ &+ 0\overset{m}{.}525\cdot \left(\frac{\varphi}{100}\right)^3 \\ &+ 0\overset{m}{.}2\cdot \left(\frac{\varphi}{100}\right)^4 \\ \end{align}\tag{17}\]
$\Delta$ ist die geozentrische Monddistanz und R der heliozentrische Erdabstand. Beide Variablen sind in AE gehalten: $a$ ist die große Halbachse des Mondes mit $a = 383397.7725$ km. Diese muss in $AE$ (⇒ Wichtige Konstanten) umgerechnet werden.
$$a = \frac{383397.7725\,km}{149597870.7\tfrac{km}{AE}} = 0.002562855813\,AE\tag{18}$$
Die Helligkeit des Erdtrabanten ist ein Mittelwert aus den Gleichungen, die im D.L. Harris (Photometry and Colormetry of Planets and Satellites) und im C.W. Allen (Astrophysical Quantities) zu finden sind.
Rotation
Der Mond hat eine gebundene Rotation, d.h., er zeigt uns immer dieselbe Seite. Dennoch gibt es Rotationselemente für unseren Erdtrabanten. $T$ und $d = JD - 2451545$ beziehen sich dabei auf die Epoche $J2000.0$. Siehe auch den Abschnitt über den Julianischer Tag JD.
| Tabelle 3 | ||||
|---|---|---|---|---|
| $\alpha_0=$ | $269\overset{\circ}{.}9949$ | $+0\overset{\circ}{.}0031\cdot T$ | $-3\overset{\circ}{.}8787\cdot \sin(E_1)$ | $-0\overset{\circ}{.}1204\cdot \sin(E_2)$ |
| $+0\overset{\circ}{.}0700\cdot \sin(E_3)$ | $-0\overset{\circ}{.}0172\cdot \sin(E_4)$ | $+0\overset{\circ}{.}0072\cdot \sin(E_6)$ | $-0\overset{\circ}{.}0052\cdot \sin(E_{10})$ | |
| $+0\overset{\circ}{.}0043\cdot \sin(E_{13})$ | ||||
| $\delta_0=$ | $66\overset{\circ}{.}5392$ | $+0\overset{\circ}{.}0130\cdot T$ | $+1\overset{\circ}{.}5419\cdot \cos(E_1)$ | $+0\overset{\circ}{.}0239\cdot \cos(E_2)$ |
| $-0\overset{\circ}{.}0278\cdot \cos(E_3)$ | $+0\overset{\circ}{.}0068\cdot \cos(E_4)$ | $-0\overset{\circ}{.}0029\cdot \cos(E_6)$ | $+0\overset{\circ}{.}0009\cdot \cos(E_7)$ | |
| $+0\overset{\circ}{.}0008\cdot \cos(E_{10})$ | $-0\overset{\circ}{.}0009\cdot \cos(E_{13})$ | |||
| $W=$ | $38\overset{\circ}{.}3213$ | $+13\overset{\circ}{.}17635815\cdot d$ | $-1\overset{\circ}{.}4\cdot 10^{-12}\cdot d^2$ | $+3\overset{\circ}{.}5610\cdot \sin(E_1)$ |
| $+0\overset{\circ}{.}1208\cdot \sin(E_2)$ | $-0\overset{\circ}{.}0642\cdot \sin(E_3)$ | $+0\overset{\circ}{.}0158\cdot \sin(E_4)$ | $+0\overset{\circ}{.}0252\cdot \sin(E_5)$ | |
| $-0\overset{\circ}{.}0066\cdot \sin(E_6)$ | $-0\overset{\circ}{.}0047\cdot \sin(E_7)$ | $-0\overset{\circ}{.}0046\cdot \sin(E_8)$ | $+0\overset{\circ}{.}0028\cdot \sin(E_9)$ | |
| $+0\overset{\circ}{.}0052\cdot \sin(E_{10})$ | $+0\overset{\circ}{.}0040\cdot \sin(E_{11})$ | $+0\overset{\circ}{.}0019\cdot \sin(E_{12})$ | $-0\overset{\circ}{.}0044\cdot \sin(E_{13})$ | |
| Tabelle 4 | |
|---|---|
| $E_1 = 125\overset{\circ}{.}045 - 0\overset{\circ}{.}0529921\cdot d$ | $E_8 = 276\overset{\circ}{.}617 + 0\overset{\circ}{.}3287146\cdot d$ |
| $E_2 = 250\overset{\circ}{.}089 - 0\overset{\circ}{.}1059842\cdot d$ | $E_9 = 34\overset{\circ}{.}226 + 1\overset{\circ}{.}7484877\cdot d$ |
| $E_3 = 260\overset{\circ}{.}008 + 13\overset{\circ}{.}0120009\cdot d$ | $E_{10} = 15\overset{\circ}{.}134 - 0\overset{\circ}{.}1589763\cdot d$ |
| $E_4 = 176\overset{\circ}{.}625 + 13\overset{\circ}{.}3407154\cdot d$ | $E_{11} = 119\overset{\circ}{.}743 + 0\overset{\circ}{.}0036096\cdot d$ |
| $E_5 = 357\overset{\circ}{.}529 + 0\overset{\circ}{.}9856003\cdot d$ | $E_{12} = 239\overset{\circ}{.}961 + 0\overset{\circ}{.}1643573\cdot d$ |
| $E_6 = 311\overset{\circ}{.}589 + 26\overset{\circ}{.}4057084\cdot d$ | $E_{13} = 25\overset{\circ}{.}053 + 12\overset{\circ}{.}9590088\cdot d$ |
| $E_7 = 134\overset{\circ}{.}963 + 13\overset{\circ}{.}0649930\cdot d$ | |
$\alpha_0$, $\delta_0$ sind wieder die äquatorialen Pole und $W$ ist der Nullmeridian des Erdtrabanten. Im Kapitel über die planetozentrischen Koordinaten können durch Einsetzen die selenozentrischen Koordinaten $A$, $D$ und den Positionswinkel $P$ der Rotationsachse des Mondes ermittelt werden. Die Libration ist dabei schon enthalten.
Die Rotationselemente des Erdtrabanten stammen aus dem Report der IAU/IAG Working Group on Cartographic Coordinates and Rotational Elements of the Planets and Satellites.
Libration
Die Libration besteht aus drei Phänomenen, die hier kurz abgehandelt werden.
Optische Libration
Die Rotationsachse des Mondes, die Achse der Ekliptik (also die Senkrechte zu dieser Ebene) und die Achse der Mondbahn (die Senkrechte zu seiner Ebene) liegen immer in derselben Ebene. Der Schnittpunkt der Äquatorialebene des Mondglobus mit der Ekliptik verläuft parallel zur Knotenlinie der Mondbahn. Der absteigende Knoten des Mondäquators fällt mit dem aufsteigenden Knoten der Umlaufbahn zusammen. Die Rotationsachse des Mondes ist zur Mondbahnebene derart geneigt, dass seine Pole temporär sichtbar werden. Dieses Phänomen wird als die optische Libration oder Libration in der selenographischen Breite genannt. Die Periode dieser Libration ist der drakonitische Monat.
Physikalische Libration
Diese Drehung entsteht vor allem durch die leicht elliptische Umlaufbahn des Mondes. Für die Libration ist die große Ungleichheit verantwortlich. Dieses Phänomen wird als die physikalische Libration oder Libration in der selenographischen Länge genannt. Die Periode dieser Libration ist der anomalistische Monat.
Diese Libration ergibt sich aus der Tatsache, dass sich die Mondkugel gleichmäßig dreht (konstante Winkelgeschwindigkeit), während die Geschwindigkeit des Mondes auf seiner Umlaufbahn um die Erde variabel ist. In der Nähe des Perigäums, wo der Abstand zur Erde am geringsten ist, ist die Umlaufgeschwindigkeit größer als in der Nähe des Apogäums, wo der Abstand zur Erde maximal ist.
Abb.3 zeigt die Bewegung des Mondes um die Erde. Die Durchmesser von Erde und Mond sowie die Abstände sind natürlich nicht maßstabsgetreu. In Position 1 befindet sich der Mond im Perigäum. Sei $m$ zu diesem Zeitpunkt der Mittelpunkt der Mondscheibe (von der Erde aus gesehen). Man kann sich vorstellen, dass sich am Punkt $m$ ein großer Mondberg befindet, der auf die Erde gerichtet ist. Nach einer Viertelumdrehung (das sind ca. $\frac{27\overset{d}{.}32}{4} = 6\overset{d}{.}83$ Tage) hat sich die Mondkugel um genau $90^{\circ}$ gedreht. Der Mond hat dann Position 2 erreicht. Diese Position ist nicht genau bei $90^{\circ}$, sondern etwa bei $96^{\circ}$ Grad weiter in der Umlaufbahn als Position 1, und Punkt $m$ scheint nun nach Osten (nach „links“) verschoben zu sein. Am westlichen Mondrand sieht man etwas über den Rand der „mittleren“ sichtbaren Mondhälfte hinaus.
Eine halbe Umdrehung nach dem Perigäum erreicht der Mond sein Apogäum und hat von der Erde aus gesehen über $180^{\circ}$ zurückgelegt (Position 3). Auch die Mondkugel hat sich nun um $180^{\circ}$ gedreht und der Punkt $m$ liegt somit von der Erde aus gesehen wieder genau im Zentrum der Mondscheibe.
Eine Viertelumdrehung nach dem Apogäum hat sich Punkt $m$ um weitere $90^{\circ}$ in Bezug auf die Achse der Mondkugel gedreht. Aber der Mond hat sich nur über ca. $84^{\circ}$ in Bezug auf die Erde bewegt, weil die Geschwindigkeit in diesem Teil seiner Umlaufbahn geringer ist als sein Mittelwert. Der Mond befindet sich jetzt in Position 4 und man sieht etwas über den mittleren östlichen Rand seiner Scheibe hinaus; Punkt $m$ ist nun etwas westlich („rechts“) der Scheibenmitte zu sehen.
Nach einer weiteren Viertelumdrehung ist der Mond mit erhöhter Geschwindigkeit ins Perigäum zurückgekehrt und befindet sich, vom Erdpunktmittelpunkt aus gesehen, wieder genau im Zentrum der sichtbaren Hälfte der Mondkugel. Der zuvor genannte Winkel von $96^{\circ}$ ist ein Mittelwert: Tatsächlich kann dieser Winkel zwischen $95^{\circ}$ und $98^{\circ}$ variieren, da der Mond neben seiner elliptischen Umlaufbahn auch durch die Anziehungskraft der Sonne gestört wird.
Topozentrische Libration
Die topozentrische Libration (auch tägliche oder parallaktische Libration genannt) entsteht zum Einen durch den unterschiedlichen Beobachtungsort (z.B. Nord- oder Südhalbkugel) und zum Anderen durch die Erdrotation. Ein Beobachter $B$ sieht den Mond bei Mondaufgang von einer etwas anderen Position als beim Monduntergang.
Der Berechnung der Libration ist ein eigenes Kapitel gewidmet.
Konstellationen
Mondalter
Das Mondalter $MA$ bezeichnet die Anzahl der Tage seit dem letzten Neumond.
$$MA = \mathrm{red}(l − L,360^{\circ})\tag{19}$$
mit der Reduktionsfunktion red(x;y). Die Berechnung stammt aus K. Hempe & J. Molt.
Perigäumspassage
Die Korrektur für die Perigäumspassage wird zum julianischen Tag dazuaddiert:
\[\begin{align} \Delta JDE_\pi &= 6793\overset{d}{.}476501\\ &+ 1\overset{d}{.}24002\cdot 10^{-4}\cdot T\\ &+ 2\overset{d}{.}2325\cdot 10^{-9}\cdot T^2\\ &- 1\overset{d}{.}3985\cdot 10^{-13}\cdot T^3 \end{align}\tag{20}\]
Knotendurchgang
Die Korrektur für den Knotendurchgang wird zum julianischen Tag dazuaddiert:
\[\begin{align} \Delta JDE_\Omega &= 3232\overset{d}{.}60542496\\ &+ 1\overset{d}{.}68939\cdot 10^{-4}\cdot T\\ &+ 2\overset{d}{.}9833\cdot 10^{-9}\cdot T^2\\ &- 1\overset{d}{.}8809\cdot 10^{-13}\cdot T^3 \end{align}\tag{21}\]
$T$ ist wieder in Jahrhunderten der Epoche $J2000.0$. $\Delta JDE_{\pi,\Omega}$ muss nur noch zum aktuellen $JD$ hinzuaddiert werden. Die beiden Ausdrücke stammen von Viktor Reijs.
Für eine ausführliche Darstellung wird auf ein eigenes Kapitel verwiesen.
Mondphasen
Dem Thema Mondphasen ist ein eigenes Kapitel gewidmet. Es geht um die Berechnung der Zeitpunkte von Vollmond, Neumond etc. Dies ist auch für die Finsternisse essentiell.