Massefunktion

Als Massefunktionen werden wichtige Parameter des Erdschwerefeldes bezeichnet, die durch harmonische Analyse des Schwerepotentials (Reihenentwicklung nach Kugelfunktionen) bestimmt werden können.

Der Begriff wurde vom Geodäten Karl Ledersteger um 1965 geprägt, der erstmals untersuchte, wie weit die wesentlichen Ergebnisse der Dynamischen Satellitengeodäsie mit theoretischen Gleichgewichtsfiguren des Erdinnern in Einklang gebracht werden können.

Physikalisch stellen die Massefunktionen jene Parameter des Schwerepotentials dar, mit denen man die wesentlichen Abweichungen des Erdschwerefeldes von dem einer mittleren Erdkugel bzw. eines theoretischen Erdellipsoids beschreiben kann. Mathematisch sind sie die harmonischen Koeffizienten einer Reihenentwicklung des Schwerepotentials, das an der Erdoberfläche durch längen- und breitenabhängige Kugelflächenfunktionen bis zu einem bestimmten Grad n durchgeführt wird. Diese zweidimensionale Lösung wird mittels geeigneter Verfahren der „Feldfortsetzung“ in den Bereich der Bahnhöhen der geodätischen Erdsatelliten transformiert.

Die größte Massefunktion des Erdkörpers ist der Schwerekoeffizient J2, der etwa 0,00108 beträgt und direkt mit der Erdabplattung zusammenhängt (21 km Unterschied zwischen mittlerem Äquator- und Polradius der Erde). Physikalisch ist J2 die Differenz der Trägheitsmomente um die Rotationsachse der Erde und um eine im Äquator liegende Achse.
Als eigentliche Massefunktionen werden jedoch erst die nachfolgenden Koeffizienten J3, J4 usw. bezeichnet, die wesentlich kleiner sind (Wirkung auf die Erdfigur nur wenige Millionstel bzw. weniger als 20 Meter). Sie beschreiben die mittlere vertikale Abweichung des Geoids vom globalen Erdellipsoid in bestimmten Breitenzonen. Dieses rotationssymmetrische Potentialmodell (Mittelung des Gepods entlang der Parallelkreise) erhält man durch harmonische Entwicklung nach zonalen Kugelfunktionen, die sich aus breitenabhängigen Legendre-Polynomen aufbauen.

Die zonalen Massefunktionen niedriger Ordnung lassen sich noch direkt deuten, was bei höheren Ordnungen nicht mehr gelingt. So bedeutet J3 eine Nord-Süd-Asymmetrie der Erdfigur, die einer unterschiedlichen Abplattung der beiden Hemisphären entspricht und etwa 15 Meter ausmacht. Sie wurde bereits durch kleine Bahnstörungen des ersten US-Satelliten Explorer 1 entdeckt und erhielt in den Medien den unglücklichen Spitznamen „Birnenform der Erde“, obwohl sie gegenüber der Erdabplattung selbst weniger als ein Tausendstel beträgt und vom Weltraum aus nicht zu beobachten wäre. Die Massefunktion 4.Ordnung (J4) kann als mittlere Abweichung der Erdmeridiane von der Ellipsenform gedeutet werden und macht unter 45° Breite ebenfalls rund 15 Meter aus (also wie J3 etwa 2 Millionstel des Erdradius). Der Koeffizient J5 ist fast Null, J6 macht nur wenige Meter aus, und die folgenden Terme lediglich einige Dezimeter.

Seit etwa 1970 wurde es möglich, die vereinfachte zonale Reihenentwicklung (J2, J3 … J14) des Geoids zu einer flächendeckenden, zweidimensionalen Analysemethode zu erweitern, in der auch die geografische Länge berücksichtigt ist. In einem solchen Erdmodell mit tesseralen (schachbrettartigen) Kugelflächenfunktionen entspricht J2 dem harmonischen Koeffizienten C(2,0), J3 entspricht C(3,0) usw. Die allgemeinen Terme der harmonischen Koeffizienten werden dann als C(n,m) und S(n,m) bezeichnet. Dabei stehen C und S für einen Cosinus- und Sinus-Term der geografischen Länge, die in der Ordnung m auftritt, während n den Grad der Breitenabhängigkeit darstellt.

In einer solchen verallgemeinerten Reihenentwicklung können sich daher breiten- und längenabhängige Wirkungen auf das Geoid zwanglos überlagern, sodass man sie ohne jegliche Vorgaben (also nur aus den vorliegenden Messungen) bestimmen kann. Dies ist heute aus den Bahnen von Erdsatelliten bis etwa Grad & Ordnung 70 möglich.

Um etwa 1975 begann man, auch höhergradige Koeffizienten des Schwerefeldes zu bestimmen, indem man die aus Satellitenbahnen abgeleiteten Terme C(n,m) mit terrestrisch bestimmten Schwerefeldparametern kombiniuerte. Seit einigen Jahren ist dies bis etwa zu Grad/Ordnung 720 möglich (n = 2, 3, 4 … 720, m = 0, 1, 2 … n),. Damit haben die Geowissenschaften ein mathematisches Verfahren zur Verfügung, mit dem das Erdschwerefeld -- je nach Bedarf -- durch bis zu 100.000 Parameter beschrieben wird. Damit lässt sich das Geoid weltweit durch rekursive (also einfach zu programmierende) Formeln mit einer Auflösung von etwa 50 km darstellen, und die Wirkung auf Satellitenbahnen auf mindestens 1 Meter. Feinere Details können durch lokale terrestrische Messungen angefügt werden, während bei geringerem Genauigkeitsbedarf die Reihenentwicklung an passenmder Stelle abgebrochen wird. So haben die kleinen GPS-Empfänger das Schwerefeld der Erde bis etwa zum Grad n=8 oder 10 einprogrammiert, was etwa 100 Massefunktionen und einer Höhengenauigkeit von ca. 5 Meter entspricht.

Literatur

• Karl Ledersteger: Astronomische und Physikalische Geodäsie (Erdmessung). In: JEK. Band V, Kapitel 6, 7 und 12, Metzler, Stuttgart 1968.
• Günter Seeber: Satellitengeodäsie. de Gruyter, um 1975 und 1995.
• Manfred Schneider: Himmelsmechanik. In 4 Bänden. Band I und III. München 1995 und 1999.
• K. Seitz: Ellipsoidisch-topografische Effekte im geodätischen Randwertproblem. Dissertation (Abstract).

Siehe auch

• Erdmessung, Physikalische Geodäsie, Potentialtheorie
• Schwerefeld, Niveaufläche, Randwertaufgabe, Störpotential
• Bahnstörung, Altimetrie, Gradiometrie