Besselsche Elemente

Die Besselschen Elemente sind geometrische Größen, die Friedrich Wilhelm Bessel eingeführt hat, um die lokalen Gegebenheiten bei einer Sonnenfinsternis an einem Beobachtungsort auf der Erde zu beschreiben. Neben Sonnenfinsternissen kann das damit verbundene Prinzip auch bei Stern- oder Planetenbedeckungen durch den Mond sowie den Transiten von Venus und Merkur vor der Sonne verwendet werden. Die bei Mondfinsternissen vorgenommene Berechnung ist ähnlich, wobei in diesem Fall der Schatten nicht auf die Erde, sondern auf den Mond fällt.

Bei Sonnenfinsternissen kann beispielsweise basierend auf den Besselschen Elementen die Bedeckungsdauer an einem bestimmten Ort ermittelt werden, oder es ist der Pfad bestimmbar, auf dem der Kernschatten des Mondes die Erdoberfläche überstreicht. Dieses Berechnungsverfahren wurde in der ersten Hälfte des 18. Jahrhunderts durch Friedrich Wilhelm Bessel entwickelt und später von William Chauvenet verfeinert.

Die grundlegende Idee des Verfahrens ist, dass die Besselschen Elemente die Bewegung des durch den bedeckenden Himmelskörper verursachten Schattens – bei Sonnenfinsternissen ist dies der Mondschatten – auf der geeignet gewählten, sogenannten Fundamentalebene wiedergeben. Dabei handelt es sich um die geozentrische Normalebene der Schattenachse. Damit ist gemeint, dass diese Ebene durch den Erdmittelpunkt geht und senkrecht auf der Achse des Schattenkegels steht, also der Geraden, die durch die Zentren des bedeckten und des bedeckenden Himmelskörpers geht.<ref name="MuckeMeeus1992">Hermann Mucke, Jean Meeus: Canon of solar eclipses: -2003 to +2526. Astronomisches Büro, Wien 1992, Seite XXXIII–LI</ref>

Zur Beschreibung der Bewegung des Schattens in dieser Ebene ist die Angabe vergleichsweise weniger Größen ausreichend – bei hinreichender Genauigkeit. Dies liegt nicht zuletzt daran, dass der Schatten während des gesamten Finsternisverlaufs in dieser Ebene immer kreisförmig ist und keiner perspektivischen Verzerrung unterliegt.<ref name="Altmann"/> In einem zweiten Schritt werden die Schattenkegel mit der Erdoberfläche geschnitten und somit die Gegebenheiten an der Erdoberfläche errechnet, wobei erst dann die annähernde Kugelform der Erde, die Erddrehung sowie die Lage und Höhe des Beobachtungsorts berücksichtigt werden müssen.

Geschichte

Dieses Verfahren, Stern- und Planetenbedeckungen sowie Sonnenfinsternisse zu beschreiben, wurde vom deutschen Wissenschaftler Friedrich Wilhelm Bessel in den 1820er Jahren ausgearbeitet. Die erste Arbeit Bessels zum Thema Sternbedeckungen findet sich in den Astronomischen Nachrichten Nr. 50 aus dem Jahre 1824, in der er einige Berechnungen auf Basis zuvor beobachteter Sternbedeckungen anstellte.<ref>Friedrich Wilhelm Bessel: Berechnung verschiedener Sternbedeckungen von den Herren Rosenberger, Strehlke und Klupsz. In: Astronomische Nachrichten Nr. 50, 3:17–28, Februar 1824</ref> Im Jahr 1829 veröffentlichte er eine verallgemeinernde Arbeit Ueber die Vorausberechnung der Sternbedeckungen in den Astronomischen Nachrichten Nr. 145.<ref>Friedrich Wilhelm Bessel: Ueber die Vorausberechnung der Sternbedeckungen. In: Astronomische Nachrichten Nr. 145, 7:1–16, September 1828</ref> Noch im selben Jahr entwickelte er die Idee weiter, indem er das Verfahren mit dem Ziel der Anwendung für Planetenbedeckungen und Sonnenfinsternisse verallgemeinerte.<ref name="Bessel182901">Friedrich Wilhelm Bessel: Beiträge zur Theorie der Finsternisse und den Berechnungs-Methoden derselben. In: Astronomische Nachrichten Nr. 151, 7:121–136, Januar 1829</ref><ref>Friedrich Wilhelm Bessel: Beiträge zur Theorie der Finsternisse und den Berechnungs-Methoden derselben (Beschluß). In: Astronomische Nachrichten Nr. 152, 7:137–144, Februar 1829</ref>

Bis zu diesem Zeitpunkt wurden zur Berechnung zwei unabhängige Verfahren mit unterschiedlichen Zielen verwendet. Das erste Verfahren diente der Bestimmung der Gegebenheiten, wie sie sich einem Beobachter an einem konkreten Ort darstellten. Die hierbei verwendete Methode ging bereits auf Johannes Kepler zurück und war später von Jérôme Lalande und Johann Gottlieb Friedrich von Bohnenberger weiterentwickelt worden. Das zweite Verfahren, das auf Joseph-Louis Lagrange zurückzuführen ist, diente der Berechnung des Zeitpunktes der Konjunktion. Da sich dieses Verfahren auf den Erdmittelpunkt bezog und keine Aussage über lokale Gegebenheiten auf der Erdoberfläche machen konnte, wurde es zur Berechnung von Finsternissen weniger häufig angewandt als das erste. Es vereinfachte jedoch viele andere astronomische Berechnungen. Bessels Ansatz bestand nun darin, Lagranges Verfahren so weiterzuentwickeln, dass damit auch die Berechnung der lokalen Gegebenheiten möglich wurde, womit er eine Kombination beider Verfahren erreichte.<ref name="Bessel182901"/>

Im zweiten Band seiner Astronomischen Untersuchungen veröffentlichte Bessel 1842 eine vier Abschnitte umfassende Abhandlung mit dem Titel Analyse der Finsternisse. Darin fasste er seine bisher veröffentlichten Arbeiten zu diesem Thema zusammen und rundete sie durch einige Ergänzungen ab.<ref name="Bessel1842">Friedrich Wilhelm Bessel: Astronomische Untersuchungen. Band 2, Königsberg 1842 (online)</ref> Diese Veröffentlichung diente als Grundlage für viele Astronomen, die sich später mit diesem Thema auseinandersetzten. Peter Andreas Hansen verwendete in seinem 1858 veröffentlichten Werk Theorie der Sonnenfinsternisse und verwandter Erscheinungen abweichend von Bessel die Schnittgerade der Ekliptik mit der Fundamentalebene als <math>x</math>-Achse. Bessels Variante, die Verwendung der Äquatorebene statt der Ekliptik, besaß jedoch einige Vorteile, wie 1863 William Chauvenet hervorhob. Chauvenet, ein amerikanischer Astronom, folgte in seinem Manual of Spherical and Practical Astronomy größtenteils dem Verfahren Bessels, entwickelte aber für einige Teilprobleme eigene Lösungsansätze. Chauvenets Darstellung war daraufhin die Basis für viele weitere Entwicklungen auf diesem Gebiet.<ref>Roberdeau Buchanan: The Mathematical Theory of Eclipses According to Chauvenet's Transformation of Bessel's Method. Seite 17f, Philadelphia/London 1904</ref>

Auch in der heutigen Zeit, in der die Berechnungen von Finsternissen nicht mehr manuell, sondern elektronisch erfolgen, haben die Besselschen Elemente ihre Bedeutung nicht verloren. Im Gegenteil, sie stellen das Bindeglied zwischen den Berechnungen des Zeitpunkts des Auftretens einer Finsternis sowie den Berechnungen der lokalen Gegebenheiten dar. Viele Computerprogramme sind auf eine der beiden Berechnungen spezialisiert, wobei die Besselschen Elemente sozusagen als Schnittstelle fungieren.<ref>Fred Espeneak (NASA): Besselian Elements of Solar Eclipses</ref>

Sonnenfinsternisse

Die an einem Ort auf der Erdoberfläche beobachtbare gegenseitigen Bedeckung von Gestirnen hängt von den Bahndaten des bedeckten entfernteren sowie des bedeckenden näheren ab. Diese Daten (Ephemeriden) werden üblicherweise als Winkel Rektaszension und Deklination angegeben. Diese Winkel beziehen sich als geozentrische Koordinaten auf den Erdmittelpunkt, so dass aus ihnen nicht direkt die Gestirnsbedeckung zu entnehmen ist, die an einem bestimmten Punkt der Erdoberfläche zu beobachten ist.

Um eine Bedeckung in einem Punkt auf der Erdoberfläche zu beschreiben, müssen die aus Tafeln entnommenen oder anderweitig bekannten Bahndaten der beiden Himmelskörper umgerechnet werden. Die Besselschen Elemente dienen der Beschreibung des Verlaufs sowie der Größe des Kern- und Halbschattens in der Fundamentalebene. Es ist einerseits nicht schwierig, den Verlauf des Schattens in dieser Ebene ausgehend von den Bahndaten der Himmelskörper zu beschreiben, zum anderen ist auch eine recht einfache Umrechnung auf einen Beobachtungspunkt möglich. Für letztere Umrechnung enthalten die Besselschen Elemente auch Angaben, wie die Fundamentalebene gegenüber dem Nullmeridian und der Äquatorebene verdreht ist.

Die Bedeckung der Sonne durch den Mond stellt im Hinblick auf die Beschreibung der Gegebenheiten auf der Erde den kompliziertesten Okkultationstyp dar, da sowohl der bedeckte Körper – die Sonne – als auch der bedeckende Körper – der Mond – nicht zu vernachlässigende Sehwinkel haben. Zudem muss die scheinbare Bewegung der Sonne während der Bedeckung berücksichtigt werden.<ref name="Green18.9">Robin M. Green: Spherical astronomy. Seite 459ff, siehe Literatur</ref>

Definition der Besselschen Elemente

Zunächst wird ein rechtwinkliges Koordinatensystem eingeführt, das als fundamentales oder Besselsches Koordinatensystem bezeichnet wird. Dabei wird von der Schattenachse ausgegangen, der Verbindungsgeraden des Zentrums von Sonne und Mond. Die Parallele der Schattenachse, die durch den Erdmittelpunkt geht, stellt die <math>z</math>-Achse des Besselschen fundamentalen Koordinatensystems dar und folgt ständig dem Schatten, das Koordinatensystem dreht sich also mit der Richtung der Schattenachse. Die Fundamentalebene steht im Erdmittelpunkt senkrecht auf dieser Achse. In der Fundamentalebene wird die Position und Größe des Kern- und Halbschattens mittels der <math>x</math>- und <math>y</math>-Koordinate beschrieben. Die <math>x</math>-Achse ist dabei die Schnittgerade der Fundamentalebene mit der Äquatorebene und weist nach Osten, die <math>y</math>-Achse weist nach Norden.

Die ersten beiden Größen der Besselschen Elemente sind die Koordinaten <math>x</math> und <math>y</math> des Schnittpunkts der Schattenachse mit der Fundamentalebene. Die Richtung der Schattenachse – die der Richtung der <math>z</math>-Achse entspricht – wird durch die Deklination <math>d</math> und den Ephemeridenstundenwinkel <math>\mu</math> angegeben. Der Radius des Halbschattenkegels in der Fundamentalebene wird durch <math>l_1</math> beschrieben, der des Kernschattenkegels durch <math>l_2</math>. <math>l_2</math> ist dabei für eine totale Finsternis negativ, für eine ringförmige positiv.<ref name="MuckeMeeus1992"/> Die Werte <math>x</math>, <math>y</math>, <math>l_1</math> und <math>l_2</math> werden in der Regel in Einheiten des Äquatorradius der Erde angegeben.

Neben diesen sechs Größen, die sich im Verlauf der Finsternis ändern, gibt es noch zwei weitere Größen, die als konstant betrachtet werden können: Die Größen <math>\tan f_1</math> und <math>\tan f_2</math> definieren die halben Öffnungswinkel des Halb- bzw. Kernschattenkegels.<ref name="MuckeMeeus1992"/>

Berechnung der Besselschen Elemente

Die für Sonnenfinsternisse verwendeten Besselschen Elemente gehen aus vom zeitlichen Verlauf der geozentrischen Positionen von Sonne und Mond, die über deren Ephemeriden verfügbar sind..<ref group="A">Typischerweise werden vom Jet Propulsion Laboratory veröffentlichte Ephemeriden (DE200/LE200) als Grundlage verwendet. Diese Ephemeriden beziehen sich auf die Massezentren der Himmelskörper. Für Finsternisse ist jedoch das Zentrum der Mond-, Planeten- oder Sonnenscheibe maßgeblich. Dies wirkt sich störend im Falle des Mondes aus, dessen Massezentrum etwa zwei Kilometer näher bei der Erde liegt als sein geometrisches Zentrum. Die Größe der dadurch verursachten Abweichung zeigt deshalb einen Zusammenhang zur Libration. Falls bei der Berechnung Besselscher Elemente solche Korrekturen der Koordinaten des Mondes vorgenommen wurden, wird dies angegeben (<math>\Delta \beta</math> und <math>\Delta \lambda</math>).</ref><ref name="Seidelmann8.12">P. Kenneth Seidelmann: Explanatory Supplement to the Astronomical Almanac. Seite 424f, siehe Literatur</ref> Eine Möglichkeit zur Berechnung des Auftretens von Sonnenfinsternissen ist, die Positionen von Sonne und Mond sofort in das fundamentale Koordinatensystem umzurechnen.<ref name="Seidelmann8.32">P. Kenneth Seidelmann: Explanatory Supplement to the Astronomical Almanac. Seite 435–441, siehe Literatur</ref> Dann kann recht leicht ermittelt werden, ob und wann die Schattenachse die Fundamentalebene innerhalb des Erdglobus durchstößt – was bedeutet, dass sich eine zentrale, also totale oder ringförmige Finsternis ereignet.

Es gibt andere Möglichkeiten, das Auftreten von Sonnenfinsternissen zu berechnen, beispielsweise über die Finsternis-Limite. Aber auch in diesem Fall müssen die Positionen von Sonne und Mond für den Finsternisverlauf in das fundamentale Koordinatensystem umgerechnet werden, um auf Basis der Besselschen Elemente lokale Gegebenheiten an jedem Ort der Erde berechnen zu können.<ref name="Green18.6"/><ref>Otto Praxl Sonnen- und Mondfinsternisse.</ref>

Basierend auf den geozentrischen Koordinaten und den Entfernungen von Sonne und Mond können die Besselschen Elemente für einen bestimmten Zeitpunkt berechnet werden. Aus Deklination <math>\delta</math> und Rektaszension <math>\alpha</math> sowie der Entfernung <math>r</math> lassen sich zunächst die Ortsvektoren von Sonne und Mond wie folgt bestimmen:<ref name="Green18.6"/>

<math>\mathbf{r}_s = r_s \begin{pmatrix}\cos \alpha _s \cos \delta _s \\ \sin \alpha _s \cos \delta _s \\ \sin \delta _s \end{pmatrix}</math>

<math>\mathbf{r}_m = r_m \begin{pmatrix}\cos \alpha _m \cos \delta _m \\ \sin \alpha _m \cos \delta _m \\ \sin \delta _m \end{pmatrix}</math>

Als Einheit für die Entfernungen dient üblicherweise der Äquatorradius der Erde. In der Literatur wird die Entfernung häufig durch die Parallaxe ausgedrückt, die den Ephemeridentafeln entnommen werden kann. Da die Parallaxe <math>P</math> sich auf den Erdradius als Basis bezieht, kann die Entfernung <math>r</math> in Einheiten des Äquatorradius durch <math>r = 1 / \sin P</math> berechnet werden.

Im Folgenden werden als erste der Besselschen Elemente die Deklination <math>d</math> und der Ephemeridenstundenwinkel <math>\mu</math> berechnet, also die Äquatorialkoordinaten der Richtung der Schattenachse. Statt des Stundenwinkels wird hierbei zunächst die Rektaszension <math>a</math> berechnet, aus dieser kann der Stundenwinkel mittels der Formel <math>\mu=\theta -a</math> ermittelt werden, wobei <math>\theta</math> der auf Greenwich bezogenen Sternzeit entspricht.

Zur Umrechnung in das fundamentale Koordinatensystem werden die Einheitsvektoren <math>\mathbf{i}</math>, <math>\mathbf{j}</math> und <math>\mathbf{k}</math>, die in Richtung der Koordinatenachsen dieses Koordinatensystems zeigen, mittels der beiden Größen <math>d</math> und <math>a</math> ausgedrückt:<ref name="Green18.6"/>

<math>

\mathbf{i} = \begin{pmatrix}- \sin a \\ \cos a \\ 0 \end{pmatrix} ; \quad \quad \mathbf{j} = \begin{pmatrix}- \cos a \sin d \\ - \sin a \sin d \\ \cos d \end{pmatrix} ; \quad \quad \mathbf{k} = \begin{pmatrix}\cos d \cos a \\ \cos d \sin a \\ \sin d \end{pmatrix} </math>

Da die Richtung der <math>z</math>-Achse der Differenz der Ortsvektoren vom Erdmittelpunkt zu Sonne und Mond entspricht, lässt sich der Einheitsvektor in Richtung der <math>z</math>-Achse auch wie folgt ausdrücken:

<math>\mathbf{k} = \frac{\mathbf{r}_s - \mathbf{r}_m}{|\mathbf{r}_s - \mathbf{r}_m|}</math>

Durch Gleichsetzen der beiden Darstellungen von <math>\mathbf{k}</math> lassen sich nun <math>d</math> und <math>a</math> und somit alle Einheitsvektoren des fundamentalen Koordinatensystems bestimmen.

Unter Verwendung dieser Einheitsvektoren können nun die Koordinaten von Sonne und Mond in diesem Koordinatensystem bestimmt werden. Aufgrund der Definition der Fundamentalebene sind die <math>x</math>- und <math>y</math>-Koordinaten von Sonne und Mond identisch. Diese stellen gleichzeitig den Schnittpunkt der Schattenachse mit der Fundamentalebene dar und sind die nächsten ermittelten Besselschen Elemente. Weiterhin wird die <math>z</math>-Koordinate des Mondes bestimmt, da diese für die Berechnung der Schattenradien benötigt wird.<ref name="Green18.6"/>

<math>x = \mathbf{r}_m \cdot \mathbf{i} = r_m \cos \delta _m \sin \left( \alpha _m - a \right)</math>

<math>y = \mathbf{r}_m \cdot \mathbf{j} = r_m \left( \sin \delta _m \cos d - \cos \delta _m \sin d \cos \left( \alpha _m - a \right)\right)</math>

<math>z_m = \mathbf{r}_m \cdot \mathbf{k} = r_m \left( \sin \delta _m \sin d + \cos \delta _m \cos d \cos \left( \alpha _m - a \right)\right)</math>

Die Winkel zwischen der Schattenachse und den Tangenten an Sonne und Mond, die die Kegelmäntel des Halb- und Kernschattens bilden, können mittels eines Hilfsdreiecks ermittelt werden. Dabei werden die Tangenten parallel verschoben, so dass sie durch den Mondmittelpunkt gehen (siehe Abbildung rechts). Hypotenuse beider Dreiecke ist die Verbindungslinie des Sonnen- und Mondmittelpunkts, die Gegenkatheten der gesuchten Winkel bilden die auf den parallel verschobenen Tangenten rechtwinklig stehenden Strecken durch den Sonnenmittelpunkt. In diesen rechtwinkligen Dreiecken ist jeweils die Länge zweier Seiten bekannt, zum einen die Entfernung zwischen Sonne und Mond, zum anderen die Länge der Gegenkathete, die beim Halbschatten der Summe aus Sonnen- und Mondradius entspricht, beim Kernschatten der Differenz dieser beiden Größen. Somit gilt:<ref name="Green18.6"/>

<math>\sin f_1 = \frac{\rho_s + \rho_m}{|\mathbf{r}_s - \mathbf{r}_m|}</math>

<math>\sin f_2 = \frac{\rho_s - \rho_m}{|\mathbf{r}_s - \mathbf{r}_m|}</math>

Um die letzten beiden noch fehlenden Besselschen Elemente <math>l_1</math> und <math>l_2</math> zu errechnen, die Radien von Halb- und Kernschatten in der Fundamentalebene, wird der Abstand der Schnittpunkte der Tangenten mit der Schattenachse von der Fundamentalebene benötigt. Für den Halbschatten liegt dieser mit <math>V_1</math> bezeichnete Punkt auf der Schattenachse zwischen Sonne und Mond und stellt die Spitze des Halbschattenkegels dar. Der Schnittpunkt <math>V_2</math> liegt ebenfalls auf der Schattenachse und ist die Spitze – also der Endpunkt – des Kernschattens. Dabei gilt:<ref name="Green18.6"/>

<math>z_{V_1} = z_m + \frac {\rho_m}{\sin f_1}</math>

<math>z_{V_2} = z_m - \frac {\rho_m}{\sin f_2}</math>

Mittels dieser Abstände der Punkte <math>V_1</math> und <math>V_2</math> von der Fundamentalebene lassen sich die Radien der Schattenkegel in dieser Ebene wie folgt ermitteln:<ref name="Seidelmann8.32"/>

<math>l_1 = z_{V_1} \tan f_1</math>

<math>l_2 = z_{V_2} \tan f_2</math>

Wenn die Kegelspitze des Kernschattens vom Mond aus gesehen hinter die Fundamentalebene fällt, also eine totale Sonnenfinsternis vorliegt, ist <math>z_{V_2}</math> negativ, im anderen Fall positiv, was bei einer ringförmigen Sonnenfinsternis der Fall ist. Entsprechend der Konvention wird auch das Vorzeichen des Kernschattenradius <math>l_2</math> so gewählt, dass dieser im Falle einer totalen Sichtbarkeit negativ angegeben wird, bei ringförmiger Sichtbarkeit hingegen positiv. Die Größen <math>z_{V_1}</math> und <math>l_1</math> sind immer positiv.

Zur Berechnung wird ein Mondradius gewählt, der eine Mittelung der Unregelmäßigkeiten des Mondrandes darstellt (<math>\rho_m = 0{,}2725076</math>). Da aber die Totalität einer Finsternis nicht vorliegt, solange durch das tiefste Mondtal scheinende Sonnenstrahlen den Beobachtungsort noch erreichen, wird zur Berechnung der Totalitätszone und -dauer auch ein zweiter, kleinerer Wert (<math>\rho_m = 0{,}272281</math>) benutzt.<ref>Fred Espenak, Jay Anderson: Total Solar Eclipse of 2008 August 01. Seite 6, März 2007 (online; PDF; 7,8 MB)</ref>

Veröffentlichung der Besselschen Elemente

Die Besselschen Elemente sind zeitabhängig. Um eine Bedeckung zu beschreiben, müssen sie daher für einen Zeitraum angegeben werden, der beispielsweise zur vollständigen Beschreibung einer Sonnenfinsternis mehrere Stunden umfasst.

Es gibt verschiedene Varianten der Veröffentlichung der Besselschen Elemente einer Sonnenfinsternis. In manchen Fällen werden die Werte aller nicht als konstant anzusehenden Elemente (also <math>x</math>, <math>y</math>, <math>d</math>, <math>\mu</math>, <math>l_1</math> und <math>l_2</math>) in Zehn-Minuten-Intervallen für den gesamten Finsternisverlauf tabellarisch angegeben. Zwischenwerte können interpoliert werden.

Eine andere Variante ist, die Besselschen Elemente für eine Referenzzeit (<math>t_0</math>) anzugeben, beispielsweise die dem Maximum nächstliegende volle Stunde in Terrestrischer Dynamischer Zeit (TDT), und zusätzlich die stündlichen Änderungen für alle nicht als konstant anzusehenden Elemente. Dies ermöglicht die Berechnung der Werte für andere Zeitpunkte des Finsternisverlaufs als lineare Funktion der Zeit.<ref name="MuckeMeeus1992"/>

Die Angabe in polynomialer Form ermöglicht eine etwas genauere Näherung gegenüber der linearen Interpolation. Dabei werden für die veränderlichen Größen zusätzlich zum Wert zum Zeitpunkt <math>t_0</math> bis zu drei Polynomkoeffizienten angegeben. Die Berechnung des Werts zu einer bestimmten Zeit erfolgt dann in folgender Form:

<math>\sum_{n=0}^{n_\mathrm{max}} a_n t^n</math>

Dabei entspricht <math>a</math> einer der veränderlichen Größen, <math>t</math> ist die Differenz zur Zeit <math>t_0</math> in Stunden.<ref name="Altmann">Michael Altmann: Helligkeitsverlauf bei Sonnenfinsternissen (II) (PDF; 71 kB)</ref>

In der Praxis wird häufig auf die vom Goddard Space Flight Center der NASA<ref>NASA Eclipse Website</ref> in polynomialer Form veröffentlichten Besselschen Elemente zurückgegriffen. Bei den im Astronomical Almanac veröffentlichten Besselschen Elementen wird der Wert für <math>d</math> aus praktischen Gründen bereits unter Anwendung der in Frage kommenden Winkelfunktionen (Sinus sowie Kosinus) angegeben, zudem die Größen <math>d'</math>, <math>\mu'</math>, die für den gesamten Finsternisverlauf näherungsweise als konstant anzusehenden stündlichen Änderungen der Größen <math>\mu</math> und <math>d</math>.<ref name="Green18.6">Robin M. Green: Spherical astronomy. Seite 450–453, siehe Literatur</ref><ref>Astronomical Almanac 2005. Seite A78ff, Stationery Office Books, 2003 (online)</ref>

Beispiel der Anwendung der Besselschen Elemente

In folgendem Beispiel werden zunächst die Besselschen Elemente für einen vorgegebenen Zeitpunkt berechnet, womit Position und Größe des Kern- und Halbschattenkegels in der Fundamentalebene zu diesem Zeitpunkt bekannt sind. Für praktische Anwendungen muss anschließend untersucht werden, wie Punkte an der Erdoberfläche relativ zu diesen Schattenkegeln liegen. Alle hierzu erforderlichen Größen sind durch die Geometrie der Erde vorgegeben. Im Beispiel wird untersucht, ob ein gegebener Ort innerhalb des Kernschattenkegels liegt.

Ermittlung der Besselschen Elemente für einen bestimmten Zeitpunkt

Die nebenstehende Tabelle enthält die Besselschen Elemente der Sonnenfinsternis vom 11. August 1999 in polynomialer Form.<ref>Fred Espeneak (NASA): Besselian Elements for the Total Solar Eclipse of 1999 Aug 11</ref> Ziel sei es nun, für 12:34:03 MESZ (entspricht 10:34:03 UT) die Position des Kernschattens in der Fundamentalebene zu berechnen.

Zunächst ist die Differenz zur Referenzzeit (11:00:00 TDT) zu ermitteln. Hierbei ist noch die Differenz <math>\Delta T</math> zwischen TDT und Universal Time (UT) zu berücksichtigen, die zum Zeitpunkt der Finsternis 63,7 Sekunden betrug:

<math>t = 10{,}5675 - 11{,}0 + \frac{63{,}7}{3600} = -0{,}414805556</math>

Die Koordinaten des Schnittpunkts der Schattenachse mit der Fundamentalebene für die gewünschte Zeit errechnen sich wie folgt:

<math>x = 0{,}070042 + 0{,}5443035\ t - 0{,}0000406\ t^2 - 0{,}0000081\ t^3 = -0{,}155744523</math>

<math>y = 0{,}502841 - 0{,}1184929\ t - 0{,}0001158\ t^2 + 0{,}0000017\ t^3 = 0{,}551972467</math>

Analog errechnen sich Deklination und Stundenwinkel (die in der Tabelle fehlenden Werte für <math>n=2</math> oder <math>n=3</math> sind mit 0 anzusetzen):

<math>d = 15{,}32734 - 0{,}012035\ t - 0{,}000003\ t^2 = 15{,}33233167</math>

<math>\mu = 343{,}687410 + 15{,}002982\ t = 337{,}4640897</math>

Ebenso lässt sich nun der Radius des Kernschattens in der Fundamentalebene für diesen Zeitpunkt berechnen:

<math>l_2 = -0{,}0036500 + 0{,}0001163\ t - 0{,}0000116\ t^2 = -0{,}003700238</math>

Der Halbschattenradius kann auf die gleiche Weise berechnet werden, er wird allerdings für die folgende Berechnung nicht benötigt.

Prüfung, ob ein gegebener Punkt zu dieser Zeit in der Totalitätszone liegt

Im ersten Schritt wurden die Besselschen Elemente der Finsternis vom 11. August 1999 für 12:34:03 MESZ berechnet. Nun soll überprüft werden, ob der Stuttgarter Schlossplatz (48° 46′ 43″ N, 9° 10′ 48″ O48.778559.1799111111111{{#coordinates:48,77855|9,1799111111111|

   |dim=
   |globe=
   |name=StuttgartSchlossplatz
   |region=DE-BW
   |type=landmark
  }}) zu diesem Zeitpunkt in der Totalitätszone lag. Hierzu werden die Koordinaten des Schlossplatzes in das fundamentale Koordinatensystem umgerechnet. Sind diese Koordinaten bestimmt, kann leicht ermittelt werden, ob dieser Punkt innerhalb des Schattenkegels liegt, da die Schattenachse ja per Definition senkrecht auf der Fundamentalebene steht. 

Zunächst sind hierzu die gegebenen geodätischen Koordinaten des Schlossplatzes (<math>\phi</math> = 48,77855° und <math>\lambda</math> = 9,17991°) einschließlich der ellipsoidischen Höhe (<math>h</math> = 295 m<ref group="A">Die ellipsoidische Höhe ist für Stuttgart etwa 49 m größer als die NHN- oder NN-Höhe.</ref>) in geozentrische Kugelkoordinaten (<math>\phi'</math> und <math>\rho</math>) umzurechnen, wobei die Länge <math>\lambda</math> unverändert bleibt. Hierfür werden die numerische Exzentrizität <math>e</math> des Rotationsellipsoids der Erde und zwei weitere, daraus abgeleitete, breitenabhängige Hilfsgrößen verwendet:<ref name="Seidelmann8.33">P. Kenneth Seidelmann: Explanatory Supplement to the Astronomical Almanac. Seite 441–446, siehe Literatur</ref>

<math>e^2 = 1 - \frac{b^2}{a^2} = 0{,}006694380</math>

Dabei ist <math>a</math> der Äquatorradius und <math>b</math> der Polradius.

<math>C = \frac{1}{\sqrt{1 - e^2\sin^2\phi}} = 1{,}001899093</math>

<math>S = \left( 1 - e^2\right) C = 0{,}995191999</math>

Mit dem Äquatorradius <math>a</math> = 6.378.137 m lassen sich die geozentrischen Koordinaten wie folgt berechnen:<ref name="Altmann"/>

<math>\phi^\prime = \arctan \frac{\left( a S + h \right) \tan \phi}{a C + h}= 48{,}5877227^\circ</math>

<math>\rho = \frac{\left( a C + h\right) \cos \phi}{a \cos \phi^\prime}= 0{,}998156295</math>

Dabei drückt <math>\rho</math> den Abstand des Schlossplatzes vom Erdmittelpunkt in Einheiten des Äquatorradius aus, <math>\phi'</math> ist der Winkel zwischen der Äquatorebene und dem vom Erdmittelpunkt zum Schlossplatz zeigenden Ortsvektor.

Als Hilfsgröße wird nun der Stundenwinkel <math>\theta</math> des Beobachtungsorts gegenüber der <math>z</math>-Achse des fundamentalen Koordinatensystems ermittelt. Dabei ist zu beachten, dass bei den Besselschen Elementen der Stundenwinkel <math>\mu</math> unter Annahme eines Ephemeridentags (entsprechend der Dynamischen Zeit, früher: Ephemeridenzeit) berechnet wird. Da aber die tatsächliche Erdrotation nicht ganz regelmäßig ist, muss <math>\mu</math> zunächst um den Zeitunterschied zwischen Dynamischer Zeit und Universal Time korrigiert werden, der <math>\Delta T</math> entspricht. Zur Berechnung der entsprechenden Winkelkorrektur ist die siderische Taglänge maßgeblich, der Unterschied zur synodischen Taglänge (Sonnentag) wird durch den Faktor 1,002738 berücksichtigt.<ref name="MuckeMeeus1992"/><ref group="A">Es ist zu beachten, dass diese Korrektur von <math>\mu</math> bereits in der tabulierten Darstellung der Besselschen Elemente enthalten sein kann, wie beispielsweise im Astronomical Almanac.</ref> Eigentlich muss die geographische Länge des Beobachters (<math>\lambda</math>) von <math>\mu</math> abgezogen werden, da beide Winkel aber in entgegengesetzter Richtung gemessen werden, ist es eine Addition.<ref name="Seidelmann8.33"/>

<math>\theta = \left( \mu - 1{,}002738 \cdot \frac{360^\circ}{86400 \mathrm{s}} \Delta T\right) + \lambda = 346{,}3778563^\circ</math>

Damit lassen sich die kartesischen Koordinaten <math>\xi</math>, <math>\eta</math> und <math>\zeta</math> des Schlossplatzes im fundamentalen Koordinatensystem wie folgt ermitteln, wobei die Neigung der Fundamentalebene gegenüber dem geodätischen Koordinatensystem durch die Deklination <math>d</math> berücksichtigt wird:<ref name="Altmann"/>

<math>\xi = \rho \cos \phi^\prime \sin \theta = -0{,}155501299</math>

<math>\eta = \rho \sin \phi^\prime \cos d - \rho \cos \phi^\prime \sin d \cos \theta = 0{,}552271870</math>

<math>\zeta = \rho \sin \phi^\prime \sin d + \rho \cos \phi^\prime \cos d \cos \theta = 0{,}816780948</math>

Der Radius der Schnittfläche des Kernschattenkegels in der durch den Schlossplatz gehenden, zur Fundamentalebene parallelen Ebene liegt näher an Sonne und Mond und ist deshalb etwas größer als der Kernschattenradius in der Fundamentalebene. Er lässt sich auf Basis des in den Besselschen Elementen angegebenen Konuswinkels des Schattenkegels (<math>\tan f_2</math>) und des Abstands des Schlossplatzes von der Fundamentalebene (<math>\zeta</math>) berechnen. Dabei ist zu beachten, dass der Kernschattenradius bei einer totalen Finsternis per Definition negativ angegeben wird.<ref name="Altmann"/>

<math>l_{2,\zeta} = l_2 - \zeta\ \tan f_2 = -0{,}007449262</math>

Der Abstand <math>r</math> des Schlossplatzes von der Schattenachse in derselben Ebene lässt sich wie folgt ermitteln:

<math>r = \sqrt{\left(\xi-x\right)^2 + \left(\eta-y\right)^2} = 0{,}000385746 < 0{,}007449262 </math>

Da der Abstand des Schlossplatzes in dieser Ebene kleiner ist als der Radius des Schattenkegels, lag der Schlossplatz also zum gegebenen Zeitpunkt innerhalb des Kernschattens. Weil es zu dieser Zeit in Stuttgart regnete, war allerdings auf dem Schlossplatz keine Beobachtung der verfinsterten Sonne möglich.<ref>www.von-zeit-zu-zeit.de: Schüchterne Sofi - die Sonnenfinsternis in Stuttgart</ref>

Durch iteratives Durchführen dieser Berechnungen für einen Zeitraum lassen sich prinzipiell die Kontaktzeiten an einem bestimmten Ort ermitteln. Es gibt aber auch direkte Verfahren, um die Kontaktzeiten zu berechnen.<ref name="MuckeMeeus1992"/>

Weitere Gestirnsbedeckungen durch den Mond

Sternbedeckungen durch den Mond

Bei Sternbedeckungen kann die Berechnung der Besselschen Elemente gegenüber Sonnenfinsternissen stark vereinfacht werden, da es ausreichend genau ist, den bedeckten Himmelskörper als unendlich weit entfernt anzusehen. Diese Annahme ermöglicht es, die Lichtstrahlen des entfernten Objekts, die das Erde-Mond-System erreichen, als parallel zu betrachteten. Damit ergibt sich, dass die Richtung der Schattenachse, also die <math>z</math>-Achse des Besselschen fundamentalen Koordinatensystems, während des gesamten Verlaufs der Bedeckung immer genau in Richtung des Sterns zeigt und damit durch die äquatorialen Koordinaten des Sterns von vornherein gegeben ist.

Eine weitere Vereinfachung gegenüber einer Sonnenfinsternis besteht darin, dass kein Kern- und Halbschattenkegel beschrieben werden muss, sondern dass es ausreicht, den „Schatten“ als senkrecht auf der Fundamentalebene stehenden Zylinder aufzufassen. Der Radius dieses Zylinders entspricht dem Mondradius, der 0,2725 des Äquatorradius der Erde entspricht. Die Angabe von variablen Schattenradien sowie Öffnungswinkeln erübrigt sich damit.<ref name="Green18.9"/><ref name="Seidelmann9.41">P. Kenneth Seidelmann: Explanatory Supplement to the Astronomical Almanac. Seite 494–497, siehe Literatur</ref>

Die Fundamentalebene wird analog zu den Sonnenfinsternissen gewählt, also die durch den Erdmittelpunkt gehende Normalebene dieser Schattenachse. Die Schnittlinie der Fundamentalebene mit der Äquatorebene ist die <math>x</math>-Achse und zeigt nach Osten, senkrecht auf dieser steht im Erdmittelpunkt die <math>y</math>-Achse und zeigt nach Norden. Wie bei Sonnenfinsternissen erfolgen alle Angaben in diesem Koordinatensystem in Einheiten des Äquatorradius.<ref name="Green18.9"/>

Anders als bei der Sonnenfinsternis wird als Bezugszeitpunkt <math>T_0</math> für die Besselschen Elemente nicht eine volle Stunde, sondern der Zeitpunkt der Konjunktion in Rektaszension gewählt, also der Zeitpunkt, zu dem Stern und Mond dieselbe Rektaszension aufweisen. Zu diesem Zeitpunkt hat die <math>x</math>-Koordinate der Zylinderachse den Wert 0, so dass in Tabellen nur noch die <math>y</math>-Koordinate der Zylinderachse in der Fundamentalebene angegeben wird. Die Besselschen Elemente einer Sternbedeckung werden wie folgt festgelegt:<ref name="Seidelmann9.41"/>

Für Prognosenberechnungen ist es ausreichend, <math>x'</math> und <math>y'</math> während des gesamten Verlaufs der Bedeckung als konstant zu betrachten.<ref name="Green18.9"/>

Bedeckung der Planeten durch den Mond

Das Verfahren der Besselschen Elemente lässt sich auf beliebige Gestirnsbedeckungen anwenden, wenn beide Gestirne hinreichend genau kugelförmig sind. Es sind lediglich die Position und Größe der Sonne durch die des betreffenden Planeten zu ersetzen. Als Ausnahmen gab Bessel 1842 lediglich die Planeten Jupiter und Saturn an, da deren Abweichung von der Kugelgestalt damals messbar war.<ref name="Bessel1842"/> Um das Sichtbarkeitsgebiet für Bedeckungen von Planeten durch den Mond vorherzusagen, kann dasselbe vereinfachte Verfahren wie bei Sternbedeckungen angewandt werden (siehe oben).<ref name="Seidelmann9.41"/>

Sollen jedoch die Kontaktzeiten genau bestimmt werden, ist eine gegebenenfalls vorhandene Abweichung des Planeten von der Kugelform zu berücksichtigen und auch, welcher Teil der Planetenscheibe zum Zeitpunkt der Bedeckung von der Sonne angestrahlt wird. Dieses Verfahren wurde 1865 von Chauvenet beschrieben, da Bessels Verfahren für die zwischenzeitlich präziser gewordenen Beobachtungsmethoden nicht mehr genau genug war. Dabei wird der von der Sonne beschienene und von der Erde sichtbare Teil des Planeten direkt betrachtet und nicht in eine Fundamentalebene abgebildet.<ref name="Chauvenet.A.347">William Chauvenet: A manual of spherical and practical astronomy. Seite 565, siehe Literatur</ref>

Transit der unteren Planeten

Beim Transit der unteren Planeten Venus und Merkur vor der Sonne ist der bedeckende Himmelskörper der Planet. Dieser kann die Sonne niemals vollständig bedecken, denn der Kernschatten ist viel zu kurz, um auf die Erde zu fallen. Auch für diese astronomischen Ereignisse werden Besselsche Elemente zur Berechnung der lokalen Gegebenheiten verwendet.<ref>Jean Meeus: Transits. Willmann-Bell, 1989, ISBN 0-943396-26-3</ref> Es kann dabei genau dasselbe Berechnungsverfahren wie bei Sonnenfinsternissen verwendet werden, der Planet übernimmt dabei die Rolle des Mondes.<ref name="Chauvenet.A.357">William Chauvenet: A manual of spherical and practical astronomy. Seite 593–598, siehe Literatur</ref>

Da die Entfernung der unteren Planeten von der Erde wesentlich größer ist als die des Mondes, besteht bei Transiten die Möglichkeit einer vereinfachten Berechnung der Zeitpunkte des Ein- und Austritts der Planetenscheibe vor der Sonne.<ref name="Seidelmann8.5">P. Kenneth Seidelmann: Explanatory Supplement to the Astronomical Almanac. Seite 471, siehe Literatur</ref> Dieses Verfahren kommt ohne die Umrechnung der Ephemeriden in das Besselsche fundamentale Koordinatensystem aus. Dabei macht man sich zu Nutze, dass die quadrierte oder zu höherer Potenz erhobene Parallaxe der Planeten so klein wird, dass sie vernachlässigt werden kann. Ausgehend von den auf den Erdmittelpunkt bezogenen Kontaktzeiten können auf diese Weise die entsprechenden Zeitpunkte an jedem Punkt der Erde berechnet werden. Das Prinzip dieser vereinfachten Berechnung geht auf Lagrange zurück und wurde von William Chauvenet verbessert, indem er die Erdabplattung berücksichtigte.<ref name="Chauvenet.A.357"/>

Mondfinsternisse

Bei einer Mondfinsternis befindet sich ein irdischer Beobachter auf dem Himmelskörper, der den Schatten wirft. Somit sieht man von allen Orten auf der Erde genau denselben Finsternisverlauf, vorausgesetzt, der Mond ist sichtbar. Bessels Verfahren eignet sich auch zur Berechnung von Mondfinsternissen, obwohl die Berechnung der lokalen Gegebenheiten auf der Erde in diesem Fall nicht erforderlich ist.<ref name="Seidelmann8.4">P. Kenneth Seidelmann: Explanatory Supplement to the Astronomical Almanac. Seite 467–470, siehe Literatur</ref>

Wie bei Sonnenfinsternissen bezieht sich das fundamentale Koordinatensystem auf die Schattenachse, die bei Mondfinsternissen immer durch den Erdmittelpunkt geht. Die Berechnung ähnelt der bei Sonnenfinsternissen. Rektaszension und Deklination der Schattenachse ergeben sich in diesem Fall direkt aus den entsprechenden Werten der Sonne, die <math>z</math>-Achse zeigt dabei aber von der Sonne weg. Somit gilt:<ref name="Seidelmann8.4"/>

<math>

a = \alpha _s + 180^\circ ; \quad \quad d = - \delta _s </math>

Auf die gleiche Weise wie bei Sonnenfinsternissen kann der geozentrische Ortsvektor des Mondes in das fundamentale System umgerechnet werden. Über dessen <math>x</math>- und <math>y</math>-Komponente kann die Lage des Mondmittelpunkts in Bezug zur Schattenachse ermittelt werden. Da alle benutzen Winkel ihren Scheitelpunkt im Erdmittelpunkt haben, werden im Gegensatz zu Sonnenfinsternissen für die Umrechnung keine Längenangaben benötigt. Die Koordinaten <math>x</math> und <math>y</math> beziehen sich auf die Einheitskugel. Die daraus abgeleiteten Winkel werden in Bogensekunden angegeben. Die dabei verwendeten Formeln entsprechen bis auf die fehlende Einheitenumrechnung den bei der Sonnenfinsternis verwendeten.<ref name="Seidelmann8.4"/>

<math>x = \cos \delta _m \sin \left( \alpha _m - a \right)</math>

<math>y = \sin \delta _m \cos d - \cos \delta _m \sin d \cos \left( \alpha _m - a \right)</math>

Daraus lässt sich der Winkelabstand des Mondmittelpunkts von der Schattenachse berechnen:<ref name="Seidelmann8.4"/>

<math>m = \sqrt{x^2 + y^2 }</math>

Die Größe der Radien von Halb- und Kernschatten werden ebenfalls als geozentrische Sehwinkel angegeben. Die Größen <math>f_1</math> und <math>f_2</math> beschreiben hierbei den Sehwinkel der Schattenradien in der Mondumlaufbahn. In nebenstehender Abbildung deutet die gestrichelte Linie die Mondumlaufbahn an. Der Winkel <math>\pi_m</math> ist der Sehwinkel des Erdradius vom Mond aus gesehen und entspricht somit der Parallaxe des Mondes. Da dieser Winkel ein Außenwinkel des Dreiecks <math>\Delta EM'V_2</math> ist, gilt für den Sehwinkel <math>f_2</math> des Kernschattens in der Mondumlaufbahn

<math>f_2 = \pi_m - v_2\!\,</math>,

wobei <math>v_2</math> der halbe Öffnungswinkel des Kernschattenkegels ist.<ref name="Green18.2">Robin M. Green: Spherical astronomy. Seite 441f, siehe Literatur</ref> Analog kann über das Dreieck <math>\Delta ES'V_2</math> eine weitere Winkelbeziehung hergeleitet werden: Der Außenwinkel <math>r_s</math> entspricht dem Sehwinkel des Sonnenradius von der Erde, der Winkel <math>\pi_s</math> der geozentrischen Parallaxe der Sonne. Somit gilt:

<math>v_2 = r_s - \pi_s\!\,</math>

Aus beiden Winkelbeziehungen lässt sich nun durch Eliminierung des Konuswinkels <math>v_2</math> der gesuchte Winkel ermitteln:

<math>f_2 = \pi_s + \pi_m - r_s\!\,</math>

In analoger Weise kann auch der geozentrische Sehwinkel des Halbschattenradius im Mondorbit ermittelt werden. Für diesen ergibt sich folgende Beziehung:<ref name="Green18.2"/>

<math>f_1 = \pi_s + \pi_m + r_s\!\,</math>

Um die Ermittlung der Kontaktzeiten der Finsternis zu unterstützen, werden aus den Größen des Kern- und Halbschattens und dem Mondradius drei weitere Hilfsgrößen abgeleitet. Dies sind die Sehwinkel für den Abstand des Mondmittelpunkts von der Schattenachse während eines bestimmten Kontakts, die aus den Sehwinkeln der Schattenradien und dem Sehwinkel <math>r_m</math> des Mondradius berechnet werden:<ref name="Seidelmann8.4"/>

<math>L_1 = f_1 + r_m\!\,;</math> Ein- und Austritt des Mondes für den Halbschatten

<math>L_2 = f_2 + r_m\!\,;</math> Ein- und Austritt des Mondes für den Kernschatten

<math>L_3 = f_2 - r_m\!\,;</math> Beginn und Ende der totalen Finsternis

Die Größen <math>x\!\,</math>, <math>y\!\,</math>, <math>m\!\,</math>, <math>L_1\!\,</math>, <math>L_2\!\,</math>, <math>L_3\!\,</math>, <math>f_1\!\,</math> und <math>f_2\!\,</math> sowie <math>\dot x</math>, <math>\dot y</math> und <math>\dot m</math> – die stündliche Änderungsraten für die korrespondierenden Größen – gelten als Besselsche Elemente einer Mondfinsternis. Sie werden für eine Referenzzeit angegeben, beispielsweise den Zeitpunkt der Mondopposition. Es gibt allerdings im Gegensatz zu Sonnenfinsternissen keine allgemein anerkannte Art und Weise der Angabe der Kenngrößen.<ref name="Seidelmann8.4"/>

Die bisher dargestellten Berechnungen verwendeten nur Winkel zur Schattenachse und kamen ohne Definition der Fundamentalebene aus. Wenn berechnet werden soll, wann bestimmte Mondkrater – also markante Punkte der Mondoberfläche – in den Kernschatten ein- oder austreten, ist dies möglich, wenn man die Fundamentalebene so wählt, dass sie durch den Mondmittelpunkt geht – in ähnlicher Weise wie für Punkte der Erdoberfläche bei Sonnenfinsternissen.

Bei Überprüfung der berechneten Kontaktzeiten und insbesondere der Kernschatten-Ein- und Austrittszeitpunkte bestimmter Mondkrater zeigen die auf diese Weise berechneten Daten keine brauchbare Übereinstimmung mit der Realität. Dies liegt zum einen daran, dass die Erde aufgrund ihrer Abplattung keinen ausreichend kreisförmigen Schatten wirft. Zum zweiten liegt es an der Erdatmosphäre, durch die sich die Schattenkegel vergrößern. Um diese Effekte zu kompensieren, ist es üblich, in die Formeln zur Berechnung der Größe des Halb- und Kernschattenkegels zwei Korrekturfaktoren einzuführen:<ref name="Seidelmann8.23">P. Kenneth Seidelmann: Explanatory Supplement to the Astronomical Almanac. Seite 428–431, siehe Literatur</ref>

<math>f_1 = 1{,}02 \left(\pi_s + 0{,}998340\ \pi_m + r_s\right)</math>

<math>f_2 = 1{,}02 \left(\pi_s + 0{,}998340\ \pi_m - r_s\right)</math>

Aber auch auf diese Weise berechnete Daten zeigen noch keine besonders präzise Übereinstimmung mit der Realität. Dies wird vor allem darauf zurückgeführt, dass die Abplattung der Erdatmosphäre noch deutlich größer ist als die der Erdoberfläche. Es wird versucht, anhand der Beobachtungsdaten verschiedener Mondfinsternisse ein genaueres Korrekturverfahren zu entwickeln.<ref name="Soulsby">Byron W. Soulsby: Improved lunar eclipse ephemerides. In: Journal of the British Astronomical Association. 100: 293–305, 1990 (online; PDF; 2,0 MB)</ref>

Literatur

• P. Kenneth Seidelmann (Hrsg.): Explanatory Supplement to the Astronomical Almanac. University Science Books, Sausalito 2006, ISBN 1-891389-45-9
• Robin M. Green: Spherical astronomy. Cambridge University Press, Cambridge 1985, ISBN 0-521-23988-5
• William Chauvenet: A manual of spherical and practical astronomy. J. B. Lippincott & Co, Philadelphia 1863 (online)

Einzelnachweise

<references />

Anmerkungen

<references group="A" />