Quelltext der Seite Besselsche Elemente

[[Datei:BesselianElementsForTotalEclipse3D.svg|mini|hochkant=1.5|Umrisse des [[Kernschatten|Kern-]] und [[Halbschatten]]s (grün) auf der Erdoberfläche (grau) und in der Fundamentalebene (rot) während einer totalen [[Sonnenfinsternis]]]]
Die '''Besselschen Elemente''' sind geometrische Größen, die [[Friedrich Wilhelm Bessel]] einführte, um die lokalen Gegebenheiten bei einer [[Sonnenfinsternis]] an einem Beobachtungsort auf der Erde zu beschreiben. Neben Sonnenfinsternissen kann das damit verbundene Prinzip auch bei [[Okkultation|Stern- oder Planetenbedeckungen]] durch den [[Mond]] sowie den [[Planetentransit|Transiten]] von [[Venus (Planet)|Venus]] und [[Merkur (Planet)|Merkur]] vor der Sonne verwendet werden. Die bei [[Mondfinsternis]]sen vorgenommenen Berechnungen ähneln der Berechnung der Besselschen Elemente, 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 1829 durch Bessel entwickelt und später von [[William Chauvenet]] verfeinert.

Die grundlegende Idee des Verfahrens ist, dass die Besselschen Elemente die Bewegung des Schattens wiedergeben, den der bedeckenden Himmelskörper –&nbsp;bei Sonnenfinsternissen ist dies der Mond&nbsp;– auf einer gedachten ''Fundamentalebene'' verursacht. Bei dieser handelt es sich um die geozentrische [[Normalebene]] der Schattenachse, in der der Erdmittelpunkt liegt und die senkrecht auf der Achse des Schattenkegels steht. Letzteres ist die Gerade, 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 geeignet gewählten 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 Werte für die ''Erdoberfläche'' errechnet, in dem die [[Schnittkurve]]n der Schattenkegel mit der Erdoberfläche bestimmt werden, 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 ==
[[Datei:Friedrich Wilhelm Bessel.jpeg|hochkant|mini|Friedrich Wilhelm Bessel (1784–1846)]]
[[Datei:William Chauvenet.jpg|mini|hochkant|William Chauvenet (1820–1870)]]
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 [[Astronomische Nachrichten|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, Februar 1824, S. 17–28, {{bibcode|1824AN......3...17R|full=}}</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, September 1828, S. 1–16, {{bibcode|1828AN......7....1B|full=}}</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, Januar 1829, S. 121–136, {{bibcode|1829AN......7..119.|full=}}</ref><ref>Friedrich Wilhelm Bessel: ''Beiträge zur Theorie der Finsternisse und den Berechnungs-Methoden derselben (Beschluß).'' In: ''Astronomische Nachrichten'', Nr. 152, 7, Februar 1829, S. 137–144, {{bibcode|1829AN......7..137B|full=}}</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 (Astronomie)|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 ([http://books.google.de/books?id=V1M_AAAAcAAJ books.google.de])</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 der amerikanische Astronom [[William Chauvenet]] hervorhob. Er 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.'' S. 17 f., Philadelphia/London 1904</ref>

Wenn auch 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: [http://eclipse.gsfc.nasa.gov/SEcat5/beselm.html ''Besselian Elements of Solar Eclipses''.] NASA</ref>

== Sonnenfinsternisse ==
Die an einem Ort auf der Erdoberfläche beobachtbare gegenseitigen Bedeckung von [[Gestirn]]en hängt von den Bahndaten des bedeckten entfernteren sowie des bedeckenden näheren ab. Diese Daten ([[Ephemeriden]]) werden üblicherweise als Winkel [[Rektaszension]] und [[Deklination (Astronomie)|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 [[Okkultation]]styp dar, da sowohl der bedeckte Körper –&nbsp;die Sonne&nbsp;– als auch der bedeckende Körper –&nbsp;der Mond&nbsp;– nicht zu vernachlässigende [[Sehwinkel]] haben. Zudem muss die [[Scheinbar (Astronomie)|scheinbare]] Bewegung der Sonne während der Bedeckung berücksichtigt werden.<ref name="Green18.9">Robin M. Green: ''Spherical Astronomy.'' S. 459 ff., siehe Literatur</ref>

=== Definition der Besselschen Elemente ===
[[Datei:BesselianElementsFundPlane.svg|mini|hochkant=1.3|Fundamentalebene mit den Besselschen Elementen <math>x</math>, <math>y</math>, <math>l_1</math> und <math>l_2</math>]]
Zunächst wird ein [[Kartesisches Koordinatensystem|rechtwinkliges Koordinatensystem]] eingeführt, das als ''fundamentales'' oder ''Besselsches Koordinatensystem'' bezeichnet wird. Dabei wird von der Schattenachse ausgegangen, der [[Verbindungsgerade]]n des Zentrums von Sonne und Mond. Die [[Parallel (Geometrie)|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 –&nbsp;die der Richtung der <math>z</math>-Achse entspricht&nbsp;– wird durch die [[Deklination (Astronomie)|Deklination]] <math>d</math> und den Ephemeriden[[stundenwinkel]] <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 Öffnungs[[winkel]] des Halb- bzw. Kernschatten[[Kegel (Geometrie)|kegels]].<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.'' S. 424 f., 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.'' S. 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 [http://www.praxelius.de/ ''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 (Astronomie)|Deklination]] <math>\delta</math> und [[Rektaszension]] <math>\alpha</math> sowie der Entfernung <math>r</math> lassen sich zunächst die [[Ortsvektor]]en 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 [[Horizontalparallaxe|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 Ephemeriden[[stundenwinkel]] <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 [[Royal Greenwich Observatory|Greenwich]] bezogenen [[Sternzeit]] entspricht.

Zur Umrechnung in das fundamentale Koordinatensystem werden die [[Einheitsvektor]]en <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>

[[Datei:BesselianElementsCalc1.svg|mini|hochkant=1.6|Ortsvektoren von Sonne und Mond sowie deren Differenzvektor, der der <math>z</math>-Achse entspricht]]
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>

[[Datei:BesselianElementsCalc2.svg|mini|hochkant=1.8|Bestimmung der Konuswinkel für Halb- und Kernschatten (<math>f_1</math> und <math>f_2</math>); aus Gründen der Übersichtlichkeit ist jeweils nur eine Tangente dargestellt]]
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 [[Parallelverschiebung|parallel verschoben]], so dass sie durch den Mondmittelpunkt gehen (siehe Abbildung rechts). [[Hypotenuse]] beider Dreiecke ist die Verbindungslinie des Sonnen- und Mondmittelpunkts, die Gegen[[kathete]]n 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 –&nbsp;also der Endpunkt&nbsp;– 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 (Zahl)|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: [http://eclipse.gsfc.nasa.gov/SEpubs/20080801/TP214149b.pdf ''Total Solar Eclipse of 2008 August 01.''] (PDF; 7,4&nbsp;MB) März 2007, S. 6</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 stündlichen Intervallen für den gesamten Finsternisverlauf tabellarisch angegeben.<ref>Jean Meeus, Carl Grosjean, Willy Vanderleen: ''Canon of Solar Eclipses.'' Pergamon Press, Oxford, 1966</ref> Zwischenwerte können [[Interpolation (Mathematik)|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 [[Terrestrische Zeit|Terrestrischer Zeit]] (TT), 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: ''[https://aie.de/eclipsecalculator/hellt2.pdf Helligkeitsverlauf bei Sonnenfinsternissen (II)] (PDF; 71&nbsp;kB)''</ref>

In der Praxis wird häufig auf die vom [[Goddard Space Flight Center]] der [[NASA]]<ref>[http://eclipse.gsfc.nasa.gov/eclipse.html 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.'' S. 450–453, siehe Literatur</ref><ref>''Astronomical Almanac 2005.'' S. A78ff., Stationery Office Books, 2003 ([http://books.google.de/books?id=5CNSCAfINYQC&pg=PT82 books.google.de])</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 ====

{| class="wikitable float-right" style="text-align:right"
|+Besselsche Elemente (Polynomkoeffizienten)<br />Referenzzeit <math>t_0</math> = 11. August 1999 11:00:00 TT
|-
! <math>n</math> || <math>x</math> || <math>y</math> || <math>d</math> || <math>l_1</math> || <math>l_2</math> || <math>\mu</math>
|-
| &nbsp;0&nbsp; ||  0,070042{{0}} || 0,502841{{0}} || 15,32734{{0}} || 0,542469{{0}} || −0,003650{{0}} || 343,68741{{0}}
|-
| &nbsp;1&nbsp; ||  0,5443035 || −0,1184929 || −0,012035 || 0,0001168 || 0,0001163 || 15,002982
|-
| &nbsp;2&nbsp; ||  −0,0000406 || −0,0001158 || −0,000003 || −0,0000117 || −0,0000116 ||
|-
| &nbsp;3&nbsp; ||  −0,0000081 || 0,0000017 ||  ||  ||  ||
|- style="text-align:center"
|colspan="7"| <math>\tan f_1</math> = 0,0046129;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<math>\tan f_2</math> = 0,0045900
|}
Die nebenstehende Tabelle enthält die Besselschen Elemente der [[Sonnenfinsternis vom 11. August 1999]] in polynomialer Form.<ref>Fred Espeneak (NASA): [http://eclipse.gsfc.nasa.gov/SEbeselm/SEbeselm1951/SE1999Aug11Tbeselm.html ''Besselian Elements for the Total Solar Eclipse of 1999 Aug 11'']</ref> Ziel sei es nun, für 12:34:03&nbsp;[[MESZ]] (entspricht 10:34:03&nbsp;UT) die Position des Kernschattens in der Fundamentalebene zu berechnen.

Zunächst ist die Differenz zur Referenzzeit (11:00:00 [[Terrestrische Zeit|TT]]) zu ermitteln. Hierbei ist noch die Differenz [[Delta T|<math>\Delta T</math>]] zwischen TT 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 [[Schloßplatz (Stuttgart)|Stuttgarter Schloßplatz]] ({{Coordinate|name=StuttgartSchlossplatz|text=DMS|NS=48/46/42.78/N |EW=9/10/47.68/E |type=landmark |region=DE-BW}}) zu diesem Zeitpunkt in der Totalitätszone lag. Hierzu werden die Koordinaten des Schloßplatzes 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 [[Ellipsoidische Koordinaten|geodätischen Koordinaten]] des Schloßplatzes (<math>\phi</math> = 48,77855° und <math>\lambda</math> = 9,17991°) einschließlich der [[Ellipsoidische Höhe|ellipsoidischen Höhe]] (<math>h</math> = 295&nbsp;m<ref group="A">Die ellipsoidische Höhe ist für Stuttgart etwa 49&nbsp;m größer als die [[Normalhöhennull|NHN]]- oder [[Normalnull|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 [[Rotationsellipsoid]]s 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>
[[Datei:BesselianElementsCalc4.svg|mini|hochkant=1.6|Zusammenhang zwischen geodätischer Breite <math>\phi</math> und geozentrischer Breite <math>\phi'</math>]]

:<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&nbsp;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 Schloßplatzes vom Erdmittelpunkt in Einheiten des Äquatorradius aus, <math>\phi'</math> ist der Winkel zwischen der Äquatorebene und dem vom Erdmittelpunkt zum Schloßplatz zeigenden [[Ortsvektor]].

[[Datei:BesselianElementsCalc3.svg|mini|hochkant=1.5|Geozentrische Winkel in der Äquatorebene: <math>a</math> ist die Rektaszension der Schattenachse, <math>\lambda</math> der Längengrad des Beobachters, <math>\theta_0</math> die auf Greenwich bezogene [[Sternzeit]], <math>\mu</math> der Stundenwinkel der Schattenachse (bereits um <math>\Delta T</math> korrigiert) und <math>\theta</math> der Stundenwinkel des Beobachters gegenüber der Schattenachse]]
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 [[Ephemeridentag]]s (entsprechend der Terrestrischer 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 Terrestrischer Zeit und Universal Time korrigiert werden, der <math>\Delta T</math> entspricht. Zur Berechnung der entsprechenden Winkelkorrektur ist die [[Siderischer Tag|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 Schloßplatzes 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><!-- Formel bei Altmann fehlerhaft, Subtraktion statt Addition -->

Der Radius der Schnittfläche des Kernschattenkegels in der durch den Schloßplatz 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 Schloßplatzes 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 Schloßplatzes 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 Schloßplatzes in dieser Ebene kleiner ist als der Radius des Schattenkegels, lag der Schloßplatz also zum gegebenen Zeitpunkt innerhalb des Kernschattens. Weil es zu dieser Zeit in Stuttgart regnete, war allerdings auf dem Schloßplatz keine Beobachtung der verfinsterten Sonne möglich.<ref>{{Webarchiv|url=http://www.von-zeit-zu-zeit.de/index.php?template=thema&theme_id=85 |wayback=20100703081124 |text=''Schüchterne Sofi – die Sonnenfinsternis in Stuttgart''. |archiv-bot=2019-05-02 14:41:40 InternetArchiveBot }} von-zeit-zu-zeit.de</ref>

Durch iteratives Durchführen dieser Berechnungen für einen Zeitraum lassen sich prinzipiell die [[Konjunktion (Astronomie)#Zeitpunkte|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 ===
[[Datei:BesselianElementsCalc5.svg|mini|hochkant=1.4|Schattenzylinder bei einer Sternbedeckung]]
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 [[Äquatoriales Koordinatensystem|ä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 (Geometrie)|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 häufig nicht eine volle Stunde, sondern der Zeitpunkt der [[Konjunktion (Astronomie)|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&nbsp;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 dann wie folgt festgelegt:<ref name="Seidelmann9.41" />
{| class="wikitable"
|-
|<math>\!\ T_0</math>|| Der Zeitpunkt der [[Konjunktion (Astronomie)|Konjunktion in Rektaszension]], angegeben in [[Universal Time]] (UT)
|-
|<math>\!\ H</math>|| Der [[Stundenwinkel]] des Sterns zum Zeitpunkt <math>T_0</math>
|-
|<math>\!\ Y</math>|| Der Wert für <math>y</math> zum Zeitpunkt <math>T_0</math>
|-
|<math>\!\ x', y'</math>|| Die zeitliche Änderung von <math>x</math> und <math>y</math> pro Stunde
|-
|<math>\!\ \alpha, \delta</math>|| Rektaszension und Deklination des Sterns ([[Sternort]]) und gleichzeitig Richtung der <math>z</math>-Achse
|}
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 (Planet)|Jupiter]] und [[Saturn (Planet)|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 [[#Sternbedeckungen durch den Mond|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.'' S. 565, siehe Literatur</ref>

== Transit der unteren Planeten ==
Beim Transit der [[Unterer und oberer Planet|unteren Planeten]] [[Venus (Planet)|Venus]] und [[Merkur (Planet)|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.'' S. 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.'' S. 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 (Mathematik)|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 [[Joseph-Louis Lagrange|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. Bei der Berechnung von Mondfinsternissen werden entsprechende Sehwinkel ([[Polarkoordinaten|polare Koordinaten]]) bestimmt, was der Ermittlung der Besselschen Elemente ([[kartesische Koordinaten]] in der Besselschen Fundamentalebene) ähnelt. Deshalb werden die für Mondfinsternisse benutzten Sehwinkel gelegentlich ''auch'' Besselsche Elemente genannt.<ref name="Seidelmann8.4">P. Kenneth Seidelmann: ''Explanatory Supplement to the Astronomical Almanac.'' Seite 467–470.</ref> Eine Fundamentalebene wird aber weder für die Erde noch für den Mond benutzt, und es handelt sich bei der Beschreibung von Mondfinsternissen in der Regel ausschließlich um polare Koordinaten.

Wie bei Sonnenfinsternissen bezieht sich das fundamentale Koordinatensystem auf die Schattenachse, die bei Mondfinsternissen aber immer durch den Erdmittelpunkt geht. Die Berechnung ähnelt der bei Sonnenfinsternissen. [[Rektaszension]] und [[Deklination (Astronomie)|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 [[Bogensekunde]]n 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>

[[Datei:BesselianElementsCalc6.svg|mini|hochkant=1.8|Geometrie des Kernschattens bei einer Mondfinsternis]]
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 [[Horizontalparallaxe|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>

[[Datei:BesselianElementsCalc7.svg|mini|Besselsche Elemente für die Kontaktzeiten einer Mondfinsternis]]
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> –&nbsp;die stündliche Änderungsraten für die korrespondierenden Größen&nbsp;– gelten als Besselsche Elemente einer Mondfinsternis. Sie werden für eine Referenzzeit angegeben, beispielsweise den Zeitpunkt der Mond[[Opposition (Astronomie)|opposition]]. 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]] –&nbsp;also markante Punkte der Mondoberfläche&nbsp;– 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 Ein- und Austrittszeitpunkte bestimmter Mondkrater in den bzw. aus dem Kernschatten zeigen die auf diese Weise berechneten Daten keine brauchbare Übereinstimmung mit der Realität. Dies liegt zum einen daran, dass die Erde aufgrund ihrer [[Erdabplattung|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> wobei der Faktor 1,02 die Vergrößerung des Erdschattens durch die Wirkung der Erdatmosphäre um 1/50 und der Faktor 0,998340 die Abplattung der Erde im [[Mittelwert]] zwischen Äquator- und Poldurchmesser kompensieren soll:
: <math>f_1 = 1{,}02 \left(0{,}998340\ \pi_m + \pi_s + r_s\right)</math>
: <math>f_2 = 1{,}02 \left(0{,}998340\ \pi_m + \pi_s - r_s\right)</math>

[[André Danjon]] wies 1951 darauf hin, dass zur Berücksichtigung der Wirkung der Erdatmosphäre die beiden Schattenkegel nicht um den gleichen relativen Betrag von 1/50 zu vergrößern sind, sondern vielmehr eine Vergrößerung um denselben absoluten Betrag den tatsächlichen geometrischen Verhältnissen entspricht.<ref name="Danjon">A. Danjon: ''Les éclipses de Lune par la pénombre en 1951.'' In: ''L'Astronomie'', 65, p. 51–53</ref> Danjon geht von einer 75 Kilometer hohen Schicht der Erdatmosphäre aus, die absorbierend wirkt, was einer Vergrößerung des Erdradius bzw. der Parallaxe des Mondes um 1/85 entspricht. Der Faktor 1,01 kombiniert diese Vergrößerung mit dem Faktor für die Erdabplattung:<ref group="A">Danjon setzt für die Abplattung der Erde den Wert von 1/297 an, als Mittelwert zwischen Äquator- und Erddurchmesser wird 0,5 x 1/297 = 1/594 zur Korrektur der Schattengröße verwendet. Damit ergibt sich die Schattenvergrößerung zu 1 + 1/85 - 1/594 ~ 1,01.</ref>

: <math>f_1 = 1{,}01\ \pi_m + \pi_s + r_s</math>
: <math>f_2 = 1{,}01\ \pi_m + \pi_s - r_s</math>

[[Mondfinsternis#Magnitude (oder Größe)|Finsternisgrößen]] für Kernschatten-Finsternisse, die nach der 1/50-Regel berechnet werden, sind im Vergleich zur Rechnung nach Danjon um etwa 0,005 zu groß, für Halbschatten-Finsternisse um rund 0,026.<ref name="MeeusMucke.XXIV">J. Meeus, H. Mucke: ''Canon of Lunar Eclipses -2002 to +2526'', Seite XXIV</ref>

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, 1990, S. 293–305, {{bibcode|1990JBAA..100..293S|full=PDF; 2,0&nbsp;MB}}</ref>

== Anmerkungen ==
<references group="A" />

== 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, [http://books.google.de/books?id=6JlzJHd15vQC books.google.de]
* [[Jean Meeus]]: ''Elements of Solar Eclipses 1951-2200.'' Willmann-Bell, Richmond 1989, ISBN 0-943396-21-2 (Mit Rechenverfahren und den Besselschen Elementen aller Sonnenfinserisse im Zeitraum 1951 bis 2200.)
* Jean Meeus: ''Transits.'' Willmann-Bell, Richmond 1989, ISBN 0-943396-25-5 (Mit Rechenverfahren und den Besselschen Elementen aller Merkur-Transits im Zeitraum 1600 bis 2300 und aller Venus-Transits -2000 bis 4000.)
* Jean Meeus: ''Astronomical Tables of the Sun, Moon and Planets'' 3rd edition. Willmann-Bell, Richmond 2015, ISBN 978-1-942675-03-7 (Mit Rechenverfahren und den Besselschen Elementen für Bedeckungen heller Sterne im Zeitraum 2010 bis 2040.)

== Einzelnachweise ==
<references />

{{Exzellent|8. Dezember 2010|82389808}}

[[Kategorie:Sonnenfinsternis]]
[[Kategorie:Mondfinsternis]]
[[Kategorie:Himmelsmechanik]]
[[Kategorie:Friedrich Wilhelm Bessel als Namensgeber]]