Als Zeitgleichung ($ZGL$) bezeichnet man die Differenz von wahrer und mittlerer Ortszeit. Die mittlere Ortszeit ist jene, die wir auf einer Uhr ablesen. Eine Sonnenuhr hingegen zeigt die wahre Ortszeit.

Abb. 1: Die wahre Bewegung der Sonne verursacht ein Analemma

Wenn die mittlere Sonne den Meridian eines Beobachters durchläuft, ist dort mittlerer Mittag. Wahrer Mittag ist jener Moment, in dem die wahre Sonne den Meridian durchläuft. Die $ZGL$ ist die Differenz zwischen der wahren Ortszeit $WOZ$ und mittleren Ortszeit $MOZ$. Mit anderen Worten, die $ZGL$ ist die Differenz zwischen den Stundenwinkeln der wahren Sonne und der mittleren Sonne.

$$\large E = WOZ - MOZ$$

Im englischen Sprachraum wird oft der Buchstabe $E$ für „Equation of time“ verwendet. Aufgrund der Exzentrizität ihrer Umlaufbahn und in viel geringerem Maße aufgrund der Störungen durch den Mond und die Planeten variiert die heliozentrische Länge $L$ der Erde nicht gleichmäßig. Daraus folgt, dass die Sonne die Ekliptik mit einer ungleichförmigen Geschwindigkeit durchläuft. Aus diesem Grund und auch dadurch, dass sich die Sonne in der Ekliptik und nicht entlang des Himmelsäquators bewegt, nimmt die Rektaszension $\alpha_{\odot}$ der Sonne nicht gleichmäßig zu.

• Man stelle sich eine erste fiktive Sonne vor, die sich mit konstanter Geschwindigkeit entlang der Ekliptik bewegt und mit der wahren Sonne am Perihel und Aphel zusammenfällt (Sonnennähe bzw. Sonnenferne).
• Weiters betrachtet man eine zweite fiktive Sonne, die sich mit konstanter Geschwindigkeit den Himmelsäquator entlang bewegt und mit der ersten fiktiven Sonne zu den Tag-Nacht-Gleichen zusammenfällt. Diese zweite fiktive Sonne ist die mittlere Sonne, und ihre Rektaszension $\alpha_{\odot}$ nimmt definitionsgemäß gleichmäßig zu.

Definiert man also diese mittlere Sonne wie beschrieben, ergibt sich ein Zeitunterschied, der durch folgende Beziehung beschrieben werden kann:

$$ E = {L_0} - 0.0057183 - {\alpha _ \odot } + \Delta \psi \cdot \cos \varepsilon$$

In dieser Gleichung ist die Konstante $0.0057183^{\circ}$ die Summe aus dem Mittelwert der Aberration in Länge $-20\overset{''}{.}49552$ und der Korrektur zur Reduktion auf das FK5-System (Fundamental-Katalog), $-0\overset{''}{.}09033$. Im Rahmen der Genauigkeit, die hier besprochen wird, kann letztere Korrektur auch weggelassen werden.

$\alpha_{\odot}$ ist die scheinbare Rektaszension der Sonne, berechnet unter Berücksichtigung der Aberration und der Nutation. Die Größe $\Delta \psi \cdot \cos \varepsilon$, wobei $\Delta \psi$ die Nutation in Länge und $\varepsilon$ die wahre Schiefe der Ekliptik ist, wird benötigt, um die scheinbare Rektaszension der Sonne auf das Äquinoktium des Datums zu beziehen, ebenso wie die mittlere Länge $L_0$.

In der gegebenen Formel sollten die Größen $L_0, \alpha_{\odot}, \Delta\psi$ in Grad ausgedrückt werden. Damit wird auch die $ZGL$ in Grad ausgedrückt. Sie kann dann durch Multiplikation mit $4\frac{m}{\circ}$ in Minuten umgewandelt werden.

$$24^{h} \hat{=}\; 360^{\circ} \; \Rightarrow \; 1^{h} \hat{=}\; 15^{\circ} \; \Rightarrow \; 4^{m} \hat{=}\; 1^{\circ}$$

Die oben angegebene Vorgehensweise entspricht der exakten Berechnung der $ZGL$. Sie setzt jedoch genaue Berechnungen der Erdkoordinaten voraus, z.B. nach der VSOP87 oder einer anderen Planetentheorie, was schnell aufwändig werden kann. Auf dieser Seite wird daher ein einfacherer Algorithmus verwendet, der schneller ist, jedoch immer noch sehr gute Ergebnisse liefert.

Der hier verwendete Algorithmus zur Berechnung der $ZGL$ wird von J. Meeus in Astronomical Algorithms angeführt. Sie stammt aus dem Textbook on Spherical Astronomy von W.M. Smart. Dieser Algorithms verwendet zur Berechnung eine Reihenentwicklung für die mittlere Anomalie $M$ und die mittlere Länge $L_0$ der Sonne, sowie die Werte der Ekliptikschiefe $\varepsilon$ und der Exzentrizität $\epsilon$ der Erdbahn. Der Wert der $ZGL$ wird bei Smart mit $E$ bezeichnet und in Radiant ausgegeben.

Ein gegebenes Datum wird zunächst in julianische Tage $JD$ umgerechnet. Mit

$$T = \frac{JD - 2451545.0}{36525}$$

berechnet man dann folgende Grundgrößen:

Damit ist die Zeitgleichung (in Radiant, Hempe & Molt) gegeben durch den Ausdruck mit $y = \tan^2(\frac{\varepsilon}{2})$

\[\begin{align} E^{rad} &= y \cdot \sin (2\cdot L_0) \\ &- 2 \cdot \epsilon \cdot \sin (M) \\ &+ 4 \cdot \epsilon \cdot y \cdot \sin (M) \cdot \cos (2\cdot L_0) \\ &- \frac{1}{2} \cdot y^2 \cdot \sin (4\cdot L_0) \\ &- \frac{5}{4} \cdot \epsilon^2 \cdot \sin (2\cdot M) \\ \end{align}\]

Die Winkelwerte für $\varepsilon, L_0, M$ werden hier in in Grad berechnet. Für die $ZGL$ achte man auf die Argumente innerhalb der Sinus- bzw. Cosinusfunktionen. Die meisten Programmiersprachen benötigen Bogenmaßwerte als Argumente der trigonometrischen Funktionen und müssen durch Multiplikation mit $\frac{\pi}{180}$ entsprechend konvertiert werden.

Die Umrechnung des Ergebnisses von $E$ in Grad erfolgt mittels Multiplikation mit $\tfrac{180}{\pi}$, danach wird mit $4$ multipliziert, um den Zeitwert in dezimalen Minuten zu erhalten. Dieser wird schließlich in Minuten und Sekunden umgewandelt, um die Ausgabe in $\pm\textrm{mm:ss}$ zu erhalten. Die Genauigkeit des hier vorgestellten Verfahrens liegt etwa bei $\pm 4^{s}$ für Daten im Zeitraum $-2000$ bis $+5000$.

Die Kurve, die den Verlauf der $ZGL$ im Laufe eines Jahres darstellt ist bekannt und in vielen Astronomiebüchern zu finden (Abb.2). Gegenwärtig (2024) hat die Kurve ein tiefes Minimum nahe dem 11. Februar, ein hohes Maximum nahe dem 3. November und ein sekundäres Maximum und Minimum um den 14. Mai bzw. am 26. Juli. Der Verlauf der $ZGL$ ändert sich jedoch im Laufe der Jahrhunderte allmählich, weil sich die Schiefe der Ekliptik, die Exzentrizität und die Länge des Perihels der Erdbahn langsam ändern.

• Roter Bereich: Die Zeitgleichung ist negativ, die Sonnenuhr geht gegenüber der mittleren Sonne nach.
• Grüner Bereich: Die Zeitgleichung ist positiv, die Sonnenuhr eilt der mittleren Sonne voraus.

Abb. 2: Zeitgleichung für das Jahr 2024

Die nachstehende Abb.3 zeigt den Verlauf der $ZGL$ von $-2000$ bis $+5000$. Auf der vertikalen Skala werden die Daten in Abständen von fünf Minuten angegeben. Man sieht zum Beispiel, dass das Minimum im Februar in ferner Zukunft weniger tief sein wird.

Abb. 3: Die Zeitgleichung im Laufe der Jahrtausende

Im Jahr 1246 n.Chr., als das Perihel der Erde mit der Wintersonnenwende zusammenfiel, war der Graph der $ZGL$ bezüglich der Null-Linie genau symmetrisch: Das Minimum des Februars war genau so tief wie das Maximum im November hoch war, und das kleinere Mai-Maximum war genau so hoch wie der tiefste Wert des Juli-Minimums.

Zwischen 1600 und 2100 n.Chr. variieren die Extremwerte der Zeitgleichung, wie in der unten stehenden Tabelle gezeigt wird. Dies sind mittlere Werte: Die Berechnung basiert auf einer ungestörten elliptischen Bewegung der Erde, und die Nutation wurde hier nicht berücksichtigt. Die Werte weichen daher von der exakten Berechnung etwas ab.

Tabelle: Haupt- und Nebenmaxima/minima