Variations Séculaires des Orbites Planétaires

Die Variations Séculaires des Orbites Planétaires (VSOP, deutsch „Säkulare Variation der Planetenorbits“) ist eine Planetentheorie, eine Methode zur Berechnung der Positionen für die Planeten des Sonnensystems mit sehr hoher Genauigkeit. VSOP87 etwa erreicht für die inneren Planeten eine maximale Winkelabweichung von einer Bogensekunde über den Zeitraum von 2000 vor Christus bis 6000 nach Chr.

Sie wurde im Jahre 1982 von Pierre Bretagnon^[1], einem Mitglied des Bureau des Longitudes in Paris, veröffentlicht (VSOP82), und seither weiterentwickelt: VSOP87, VSOP2000 und VSOP2002. Sie gilt – mit den Erweiterungen – noch heute (2006) als Referenz für die numerische Modellierung der Dynamik des Sonnensystems.

VSOP82

Die VSOP82-Theorie^[2] ist eine Methode zur Berechnung der Planetenpositionen für die Planeten Merkur bis Neptun.

Die elliptische Bahn eines Planeten im Zweikörpersystem lässt sich durch Angabe seiner sechs Bahnelemente vergleichsweise einfach beschreiben. VSOP82 beschreibt die Bahnen der großen Planeten ebenfalls durch Angabe der Bahnelemente. Diese Bahnelemente sind jedoch zeitlich veränderlich. Sie beschreiben diejenige Keplerbahn, die sich der tatsächlichen Planetenbahn zum gegebenen Zeitpunkt bestmöglich anschmiegt (sog. oskulierende Bahn). Diese Bahnelemente werden über geeignete Potenzreihen berechnet.

VSOP87

Die VSOP87^[3] ist eine Weiterentwicklung der VSOP82-Theorie, die im Jahre 1987 von Bretagnon und G. Francou veröffentlicht wurde. Bei ihr besteht die Möglichtkeit, den Rechenaufwand durch Weglassen hinterer Terme zu vermindern – natürlich auf Kosten der Genauigkeit. Außerdem bietet sie direkte Berechnung der heliozentrischen Koordinaten.

Es stehen mehrere Varianten der Theorie zur Verfügung:

VSOP87:	enthält ähnlich wie die weniger genaue Vorgängerversion VSOP82 Reihenentwicklungen für die (veränderlichen) Bahnelemente der Planeten. Nach der Bestimmung der für den gewünschten Zeitpunkt gültigen Bahnelemente muss aus diesen anschließend mit den üblichen Methoden der Ephemeridenrechnung die Planetenposition errechnet werden.
VSOP87A:	enthält Reihenentwicklungen, welche unmittelbar die heliozentrischen kartesischen Koordinaten der Planeten für das Standardäquinoktium J2000.0 liefern.
VSOP87B:	Reihenentwicklungen der heliozentrischen sphärischen Koordinaten (ekliptikale Länge, ekliptikale Breite und Radiusvektor) der Planeten für J2000.0.
VSOP87C:	Reihenentwicklungen der heliozentrischen kartesischen Koordinaten für das Äquinoktium des Datums
VSOP87D:	Reihenentwicklungen der heliozentrischen sphärischen Koordinaten für das Äquinoktium des Datums
VSOP87E:	Reihenentwicklungen der baryzentrischen kartesischen Koordinaten für J2000.0.

Neben der Bequemlichkeit, unmittelbar die gewünschten Koordinaten zu liefern, bieten die Varianten A bis E auch den Vorteil, bei geringeren Genauigkeitsansprüchen die Berechnung der Reihen abbrechen zu können, sobald die gewünschte Genauigkeit erreicht ist. Bei Verwendung der VSOP87 selbst wäre es in diesem Fall schwierig zu bestimmen, mit welcher Genauigkeit die von dieser Variante gelieferten einzelnen Bahnelemente berechnet werden müssen, um letztlich die daraus folgenden Koordinaten mit der gewünschten Genauigkeit zu erhalten.

Die VSOP87A–E beruht auf einer Potenzreihenentwicklung im Argument der Zeit bis in die 5. Potenz, deren jeweiligen Faktoren durch eine Fourieranalyse aufgeschlüsselt sind. Diese ist in Tabellen absteigenden Beitrags erfasst, so dass anhand der Koeffizienten der Beitrag zum Gesamtfehler abgeschätzt werden kann.

VSOP2000

Seit einigen Jahren gibt es eine Aktualisierung, die VSOP2000^[4] von Xavier Moisson und Pierre Bretagnon, die um den Faktor 10–100 genauer als die Vorgängerversionen ist und Fehler von nurmehr einigen 0,1 mas für Merkur, Venus und Erde für das Intervall 1900–2000 aufweist.

VSOP2002

Bretagnons letzte Arbeit war die Implementierung relativistischer Effekte, und eine weitere Steigerung um den Faktor 10 – die VSOP2002 blieb aber unvollendet und zeigt Schwächen bei Uranus and Neptun.^[5]

Publikation

Obwohl die Konstruktionsmethoden der VSOP82 und VSOP87 sowie ihre Eigenschaften in der astronomischen Literatur beschrieben wurden ^[2][3], sind diese Theorien selbst in den Publikationen nicht enthalten. Sie konnten ursprünglich nur auf Magnetband bezogen werden, sind aber inzwischen über das Internet erhältlich. Für Anwendungen mit geringeren Genauigkeitsansprüchen sind in dem Buch „Astronomische Algorithmen“ von Jean Meeus^[6] oder vom Österreichischen Astronomischen Verein^[7] Auszüge dieser Listen periodischer Terme veröffentlicht worden.

Beispiel

Für das tropische Jahr ergibt sich

1. gemäß der VSOP 87:

   365.242189623 – T × 0.000061522 – T² × 0.0000000609 + T³ × 0.00000026525

2. gemäß der VSOP2000:

   365.2421905166 – T × 0.000061560 – T² × 0.0000000684 + T³ × 0.0000002630 + T⁴ × 0.0000000032

T in julianischen Jahrtausenden (1000 × 365,25 Tage bezüglich J2000.0, d. h. bzgl. T = JD – 2 451 545.0) sind.

Literatur

1. ↑ Pierre Bretagnon (1943-2002)
2. ↑ ^{a b} P. Bretagnon: Théorie du mouvement de l'ensemble des planètes. Solution VSOP82. In: Astronomy and Astrophysics. Nr. 114, 1982, S. 278–288 (pdf 1,2 MB, abgerufen am 19. Juni 2006).
3. ↑ ^{a b} P. Bretagnon, G. Francou: Planetary theories in rectangular and spherical variables. VSOP87 solutions. In: Astronomy and Astrophysics. Nr. 202, 1988, S. 309–315 (pdf 841 KB, abgerufen am 19. Juni 2006).
4. ↑ X. Moisson, P. Bretagnon: Analytical Planetary solution VSOP2000. In: Springer (Hrsg.): Celestial Mechanics and Dynamical Astronomy. 80, Nr. 3–4, Juli 2001, S. 205–213 (abstract, abgerufen am 18. November 2008).
5. ↑ A. Fienga , J.-L. Simon: Analytical and numerical studies of asteroid perturbations on solar system planet dynamics. In: Astronomy and Astrophysics. Nr. 429, 2005, S. 361–367, doi:10.1051/0004-6361:20048159 (c ESO 2004, Webdokument, pdf).
6. ↑ Jean Meeus: Astronomical Algorithms. ISBN 0-943396-35-2.
7. ↑ Hermann Mucke: Wandelgestirnörter. In: Mucke (Hrsg.): Moderne astronomische Phänomenologie. 20. Sternfreunde-Seminar, 1992/93. Zeiss Planetarium der Stadt Wien und Österreichischer Astronomischer Verein, 1992, 2. Berechnen des heliozentrischen Orts der großen Planeten Merkur bis Neptun – Die Planetentheorien VSOP82 und VSOP87, S. 1–23 (Weblink. Sternfreunde, abgerufen am 19. Juni 2006. ).

Weblinks

• Die VSOP87 auf dem FTP-Server des Institut de mécanique céleste et de calcul des éphémérides (IMCCE) (Stand: 5. April 2005)