Dieses Kapitel befasst sich mit der Korrektur von Koordinaten, die durch spezielle himmelsmechanische, atmosphärische und relativistische Phänomene und Effekte auftreten. Nachstehend gilt für julianische Jahrhunderte $T$ bezüglich der Epoche $J2000$:

$$T = \frac{JDE - 2451545.0}{36525}\tag{1}$$

Abb. 1

In der Astronomie ist mit der Präzession eine Richtungsänderung der Erdachse gemeint, die durch die Anziehung des Mondes und der Sonne in Verbindung mit der Abplattung auftritt. Sie äußert sich durch das Fortschreiten des Frühlingspunkts entlang der Ekliptik, woraus sich auch die Bezeichnung Präzession (lateinisch für ‚Fortschreiten‘) herleitet. Die Präzessionsgleichungen beruhen auf O. Montenbruck. Zuerst werden einige Hilfswinkel benötigt; diese müssen vorab in Grad umgerechnet werden.

Abb. 2

Tabelle 1
Winkel	von	bis
$\Psi$	$1$	$(J2000)$
$\chi$	$1$	$(t)$
$90^\circ + z$	$(t)$	$3$
$90^\circ - \zeta$	$(J2000)$	$3$
$\Lambda$	$(t)$	$2$

$$\alpha' = \alpha + \zeta + z + \theta\cdot\tan(\delta)\cdot\sin(\alpha)\tag{2}$$ $$\delta' = \delta + \theta\cdot\cos(\alpha)\tag{3}$$

$$\lambda' = \lambda + p + \eta\cdot\tan(\beta)\cdot\cos(\lambda - \Pi)\tag{4}$$ $$\beta' = \beta - \eta\cdot\sin(\lambda - \Pi)\tag{5}$$

Bahnelemente:

$$i' = i - \eta\cdot\cos(\Pi - \Omega)\tag{6}$$ $$\Omega' = \Omega + p - \eta\cdot\frac{\sin(\Pi - \Omega)}{\tan(i)} \quad \forall i \neq 0^{\circ}\tag{7}$$ $$\omega' = \omega + \eta\cdot\frac{\sin(\Pi - \Omega)}{\sin(i)}\tag{8}$$

Abb. 3

Die astronomische Nutation (lat.: nutare ‚nicken‘) ist der schwankende Teil der Präzession der Erdachse im Raum unter dem Einfluss von Sonne und Mond. Die Stärke dieser Schwankungen wird durch die Präzession der Mondbahn verursacht, deren Knotenlinie mit einer Periode von $18.6$ Jahren die Erde umläuft. Die zur Berechnung der Nutation benötigten Mittelwerte lauten:

Es gilt:

$$\lambda' = \lambda + \frac{\Delta\lambda}{3600\tfrac{''}{^\circ}}\tag{9}$$ $$\varepsilon' = \varepsilon + \frac{\Delta\varepsilon}{3600\tfrac{''}{^\circ}}\tag{10}$$

Abb. 4

Die Aberration (lat. aberratio: Verirrung) beschreibt die scheinbare Position eines Himmelsobjekts durch einen bewegten Beobachter, auch bedingt durch die Endlichkeit der Lichtgeschwindigkeit $c$. Der Ort des Sterns oder Planeten erscheint verschoben und muss korrigiert werden.

Die Konstanten:

Nimmt man die Erdbahngeschwindigkeit $v_1 = 29.785 \frac{km}{s}$ und die Äquatorgeschwindigkeit $v_2 = 465.096 \frac{m}{s}$, so erhält man die 1. Aberrationskonstante $k_1$ bzw. 2. Aberrationskonstante $k_2$:

$$k_1 = \dfrac{20\overset{''}{.}4934472502}{3600''}\tag{11}$$ $$k_2 = \dfrac{0\overset{''}{.}3200138652125}{3600''}\tag{12}$$

Die tägliche Aberration:

Die tägliche Aberration ist die Abweichung durch die Erdrotation.

\[\begin{align} A_t &= k_2\cdot \frac{\cos(\beta_0) \cdot \cos(\theta - \alpha)}{15^h \cos(\delta)}\\ D_t &= k_2\cdot \cos(\beta_0) \cdot \sin(\theta - \alpha) \cdot \sin(\delta) \end{align}\tag{13}\]

Die neuen scheinbaren, äquatorialen Koordinaten bekommt man durch Addition: $$\alpha_t = \alpha_i - A_t \quad und \quad \delta_t = \delta_i - D_t\tag{14}$$

Die jährliche Aberration:

Die jährliche Aberration ist die Abweichung durch den Umlauf der Erde um die Sonne: \[\begin{align} A &= k_1\cdot \frac{\sin(L) \cdot \sin(\alpha) + \cos(L) \cdot \cos(\alpha) \cdot \cos(\varepsilon)}{15^h \cos(\delta)} \\\\ D &= k_1\cdot (\sin(\delta) \cdot \cos(\alpha) \cdot \sin(L)\\ &+ (\sin(\varepsilon) \cdot \cos(\delta)\\ &- \cos(\varepsilon) \cdot \sin(\delta) \cdot \sin(\alpha)) \cdot \cos(L)) \\\\ A_e &= k_1\cdot \epsilon\cdot \frac{\sin(\varpi) \cdot \sin(\alpha) + \cos(\varpi) \cdot \cos(\alpha) \cdot \cos(\varepsilon)}{15^h \cos(\delta)} \\\\ D_e &= k_1\cdot \epsilon (\sin(\delta) \cdot \cos(\alpha) \cdot \sin(\varpi)\\ &+ (\sin(\varepsilon) \cdot \cos(\delta)\\ &- \cos(\varepsilon) \sin(\delta) \cdot \sin(\alpha)) \cdot \cos(\varpi)) \end{align}\tag{15}\]

Die neuen scheinbaren, äquatorialen Koordinaten bekommt man wieder durch eine Addition:

$$\alpha_i = \alpha - (A + A_e) \quad und \quad \delta_i = \delta - (D + D_e)\tag{16}$$

Die Thematik der Abplattung wird aufgrund des umfangreichen Inhalts in einem eigenen Kapitel illustriert.

Abb. 5

Die Refraktion ist die Richtungsänderung eines Lichtstrahls auf seinen Weg durch die geschichtete Atmosphäre. Ursache ist der Anstieg des Brechungsindexes. Sie bewirkt eine Anhebung des Himmelsobjekts in der Höhe. Es werden die Formeln angegeben, mit denen man aus der »luftlosen«, wahren Höhe $h$ mittels einer Addition der Refraktion $R(h)$ die scheinbare Höhe $h'$ berechnet werden kann. J. Meeus stellt in seinen Astronomical Algorithms eine Zahlenwert-Gleichung von G.G. Bennet vor, die das berechnet:

$$R(h) = \frac{1}{\tan\left(h + \dfrac{7\overset{\circ}{.}31}{h + 4\overset{\circ}{.}4}\right)}\tag{17}$$

Dabei ist $h$ in Grad einzusetzen, das Ergebnis erhält man in Bogenminuten. Die Gleichung versagt jedoch im Zenit ($h = 90^\circ$):

• $R(90^{\circ}) = -0\overset{'}{.}0013515$, statt Null. Das wird durch Addition von $C = +0\overset{'}{.}0013515$ zur rechten Seite von $R$ bewerkstelligt
  (Für alle anderen Höhen gilt $C = 0'$.
• Laut Bennet ist die Gleichung für alle Höhen $h$ gültig mit einem maximalen Fehler von $0\overset{'}{.}07$, der ca. bei $h \approx 12^{\circ}$ auftritt.

Abb. 6: Refraktion R(h) nach Bennet

Auch für das umgekehrte Problem gibt es eine einfache Näherungsformel, wenn man aus der scheinbaren Höhe $h'$ die Refraktion ermitteln möchte:

$$R(h') = \frac{1.02}{\tan\left(h' + \dfrac{10\overset{\circ}{.}3}{h' + 5\overset{\circ}{.}11}\right)}\tag{18}$$

• Für $h' = 90^{\circ}$ ist auch hier die Refraktion $R(h')\ne 0$, was mit dem Korrekturfaktor $+0\overset{'}{.}0019279$ beseitigt werden kann.
• Die Formel stammt von Sæmundsson und ist konsistent zu Bennets Gleichung, mit einer Differenz von $4''$.

Der nächste Punkt ist eine weitere Korrektur $K$ (in Bogenminuten) des Refraktionsverlaufs für alle Höhen, bis auf den Zenit mit:

$$K = - 0.06\cdot \sin(14.7\cdot R(h) + 13.0)\tag{19}$$

Als Letztes wird noch der Druck $p$ in [hPa] und die Temperatur $T$ in [°C] berücksichtigt. Sind die obigen Berechnungen abgeschlossen, so muss noch der Wetterterm $W$ multipliziert werden:

$$W(p,T) = \frac{p}{1013.246 \ \textsf{hPa}} \cdot \frac{283.16 \ K}{273.16 \ K + T}\tag{20}$$

Es gilt:

$$R(h) = W(p,T)\cdot\big{(}R(h) + C + K\big{)}\tag{21}$$

Die nachfolgende Tabelle stammt aus G.D. Roth (Handbuch für Sternfreunde) und zeigt die zu erwartenden Ablenkungen $R$ bei verschiedenen Temperaturen $T$ für $p = 1013.246 \textsf{hPa}$.

Die Höhe $h$ wird dann korrigiert mit

$$h' = h + \frac{R(h)}{60'}\tag{22}$$

Das beobachtete Licht des Himmelsobjekts braucht zur Erde eine einige Minuten, bis zu einigen Stunden. Weil sich der Himmelskörper während dieser Zeit fortbewegt, befindet sich die beobachtete Position ein wenig hinter dem tatsächlich berechneten Ort zurück. Da die Lichtgeschwindigkeit $c$ erheblich größer ist als die Bewegung des beobachteten Objekts, liegt die Korrektur nur im Bereich von einigen Bogensekunden bis zu einer Bogenminute. Daher reicht hier die Näherung aus.

Abb. 7

Der Beobachter $E$ hat zum Zeitpunkt $t$ den Eindruck, Licht von $P$ zu empfangen. Er erhält jedoch tatsächlich ein Signal, welches zum Zeitpunkt $t − \tau$ von $P'$ ausgesandt wurde. In der Zeit $\tau$ hat sich das Objekt auf der (als gekrümmt dargestellten) Bahn weiterbewegt.

\[\begin{align} \tau &= \frac{1 AE}{1c\cdot 24^h\cdot 3600^s} \cdot \Delta \text{ Tage}\\ &= 0\overset{d}{.}005775518304 \cdot \Delta\\ &\approx \frac{\Delta}{173.14} \text{ Tage} \end{align}\tag{23}\]

Die Korrektur der Lichtlaufzeit $\tau$ liefert die astrometrischen Koordinaten.

Als Extinktion wird die Abschwächung des Lichts während seines Durchlaufs durch die Atmosphäre bezeichnet. Die folgenden Gleichungen stammen aus G.D. Roth (Handbuch für Sternfreunde). Der Luftmassenfaktor $M$ reduziert das einfallende Licht auf gleichlange Wege und bezeichnet gleichzeitig die Dicke der Luftmasse. Der Luftmassenfaktor M(z):

$$M(z) = \frac{1}{\cos(z - \Delta(z))}\tag{24}$$

mit

$$\Delta(z) \approx 0.05^{\circ}\cdot \frac{z}{93^{\circ} - z} \quad (0^{\circ} \leq z \leq 87^{\circ})\tag{25}$$

Im Zenit ist $M = 1$. Die Koeffizienten der Dunstextinktion $k_d$ und der Rayleigh - Streuung $k_r$ lauten:

$$k_r = \frac{0.00906\;p}{1013.25\;\textsf{hPa}}\cdot \left(\frac{\lambda}{1 \textsf{µm}}\right)^{-4}\tag{26}$$ $$k_d = 1.086\cdot \beta\cdot \left(\frac{\lambda}{1 \textsf{µm}}\right)^{-\alpha}\tag{27}$$

$\alpha$ ist hier der Wellenlängenexponent. Normalerweise ($\beta$ = 0.1) ist $\alpha = 1.3$, bei Staubstürmen ist $\alpha = 0.5$. Bei Wolken ist $\alpha = 0$ und $\beta$ sehr hoch. $\lambda$ ist die beobachtete Wellenlänge. Gelbes Licht: $\lambda = 550\;nm$. Im Zenit ist die Lichtausbeute $\mu (z) = 1$. Im Allgemeinen gilt deshalb $\beta\approx$ 0.1.

Es gilt $k_r$ + $k_d = 0.335252293177$. Das Licht eines Sterns im Zenit ist schon leicht geschwächt. Zur Ermittelung der Lichtausbeute wird deshalb der bereits angedeutete Faktor $\beta$ eingeführt, der den Zustand der Atmosphäre charakterisiert:

• $\beta = 0.05$: Himmel absolut klar.
• $\beta = 0.10$: Himmel leicht getrübt.
• $\beta = 0.20$: Himmel stark getrübt.

Als nächstes wird die optische Durchlässigkeit $T$ (Opazität) mit Hilfe von $M$ gebraucht. Sie gibt die Abschwächung des Lichts an.

$$m(z) = (k_r + k_d)\cdot M(z)\tag{28}$$ $$T(z) = 10^{-0.4 \cdot m(z)}\tag{29}$$

Die Extinktion $E$ ist dann die gewünschte Korrektur der Helligkeit des Himmelobjekts.

$$E(z) = (k_r + k_d + D(\lambda)\cdot t)\cdot (M(z) - 1)\tag{30}$$ $$\mu(z) = 10^{-0.4\cdot E(z)}\tag{31}$$

Der Term $D(\lambda)$ ist der Driftkoeffizient und kann je nach meteorologischen Bedingungen 0.001 $\tfrac{mag}{h}$ bis 0.1 $\tfrac{mag}{h}$ betragen. $E$ muss dann nur noch zur Größenklasse des Stern addiert werden. t ist die Zeit in Stunden, in der sich der Koeffizient D ändert. Die Extinktionstabelle stammt von der Sternwarte Höfingen.

Es muss eine bestimmte Reihenfolge eingehalten werden (Expl. Suppl.):

Die Lichtlaufzeit erfolgt als erstes, dann Aberration, Präzession und zum Schluss die Nutation. Bei der Parallaxe, Abplattung, Refraktion und Extinktion ist die Reihenfolge unwichtig.

Obwohl die Begriffe im Astronomie Lexikon erwähnt werden, sind sie hier nochmal zusammengefasst mit einer Tabelle zur besseren Übersicht wiedergegeben.

Die geometrischen Koordinaten eines Planeten ist die wahre Position eines Himmelskörpers zu einem gegebenen Zeitpunkt. Das heißt, dass die Aberration und Lichtlaufzeit nicht berücksichtigt wurde. Es handelt sich um die mittleren Bahnelemente und die heliozentrischen Koordinaten, die nur störungstheoretische Korrekturen erfahren haben. Sowie daraus durch Koordinatentransformation erzeugte geozentrische Koordinaten, sei es ekliptikal oder äquatorial.

Es handelt sich hier um die ekliptikalen oder äquatorialen Koordinaten, die durch die Lichtlaufzeit korrigiert wurden und deshalb in temporäre Sternkarten eingetragen werden können, wenn die Sternpositionen durch die Parallaxe und die Eigenbewegungen ebenfalls korrigiert wurden.

Das sind tatsächlich ursprünglich geo- oder astrometrische Koordinaten, die durch die Korrekturen wie Aberration, Refraktion und die Parallaxe berechnet werden und die wirkliche Position am Himmel wiedergeben, wie sie von einem Beobachter gesehen wird.

Werden Koordinaten mit mittleren Äquinoktium durch die Nutation korrigiert, so erhält man die wahren, ekliptikalen Koordinaten.