Radio-Okkultation

Radio-Okkultation, auch als Okkultationsmethode bezeichnet, ist eine Messtechnik zur Sondierung planetarer Atmosphären unter Benutzung phasentreuer Radiosignale, die sich durch die Atmosphäre von einem Sender zu einem Empfänger ausbreiten. Sowohl Sender als auch Empfänger befinden sich während der Messphase außerhalb der zu sondierenden Atmosphäre.

Abbildung 1: Zur Geometrie der Radio-Sondierungsmessung. Zum Zeitpunkt t₁ befindet sich die Raumsonde noch oberhalb der Atmosphäre und der Radiostrahl bleibt von der Atmosphäre unbeeinflusst. Zum Zeitpunkt t₂ durchläuft der Radiostrahl die Atmosphäre und wird gebeugt.

Die Durchführung der Messung erfordert eine spezielle Geometrie der Raumsonde zur Empfangsstation, wobei die Raumsonde während der Messung aus Sicht des Empfängers hinter dem Planeten verschwindet und somit in Okkultation geht.

Während der Radiostrahl vom Weltraum oberhalb der Atmosphäre zu einem Punkt in die Atmosphäre läuft, findet eine kontinuierliche Aufzeichnung der Beobachtungsdaten statt. Das sondierte Medium wirkt in charakteristischer Weise auf das Radiosignal und verändert dessen Phasenwinkel, Amplitude und Polarisation.

Die Signalbeeinflussung durch das Medium erzeugt einen zeitabhängigen Datensatz, der dem Höhenprofil des Brechungsindex entspricht. Dieses Profil der Neutral-Atmosphäre ist für Gasgemische proportional zur Dichte, woraus sich mit der hydrostatischen Grundgleichung und dem idealen Gasgesetz Höhenprofile von Druck und Temperatur der Neutral-Atmosphäre berechnen lassen. Aus der Proportionalität der Elektronendichte zum Höhenprofil des Brechungsindexes sind zusätzlich Aussagen über die Elektronendichte der Ionosphäre möglich.

Geschichte

Die Idee der stellaren Okkultation zur Erforschung planetarer Atmosphären, wie sie erstmals 1904 von Anton Pannekoek in den „Astronomischen Nachrichten“ veröffentlicht worden war, diente 1962 Von R. Eshleman von der Stanford-Universität in Kalifornien als Vorlage für das Radio-Okkultationsexperiment.

Gleichzeitig entwickelten D. L. Cain und Kollegen vom JPL Fehleranalysen zur Satellitennavigation mit den damals neuen masergesteuerten Antennensystemen. Hauptaugenmerk der Analyse waren Fehlerbeiträge der Erdatmosphäre und -ionosphäre auf das Ortungssignal. Cain und seine Kollegen bemerkten, dass die Atmosphäre eines anderen Planeten eine nicht zu vernachlässigende Fehlerquelle darstellt, woraus aber wiederum Aussagen über die Planetenatmosphäre gewonnen werden können. Der Unterschied zu Eshleman bestand in den auf der Gesamtlichtlaufzeit basierenden Beobachtungsmethoden. Während das Stanford-Team um Eshleman von einem Einwegmodus ausging, wobei ein Sender an Bord des Raumfahrzeugs das Referenzsignal generieren sollte, schlugen Cain und seine Mitarbeiter vom JPL ein Zweiwegverfahren vor, wobei mit Hilfe der damals neuen Masertechnologie ein an Genauigkeit nicht zu übertreffendes Referenzsignal an der Bodenstation generiert werden sollte, welches dann vom Raumfahrzeug empfangen und leicht modifiziert phasentreu wieder zurückgesandt werden sollte.

Da das Einwegeverfahren mit der damaligen Technik noch nicht genau genug war, entschied die NASA, das Zweiwegeverfahren während der Mariner-4-Mission erstmals einzusetzen.

Die Bedeutung der neuen Messmethode wurde in den frühen 1960er-Jahren bei dem Mariner-Projekt, bei der auch eine Mars-Landung vorgesehen war, deutlich. Eine solche Landung stellt hohe Anforderungen an das Design und an die Landungsstrategie des Raumfahrzeugs. Genaue Kenntnisse der Marsatmosphäre waren unabdingbar, und obwohl Informationen aus bodengestützter Spektroskopie vorhanden waren, gab es Zweifel an den damals gemessenen Werten aus dem Jahr 1950 von de Vaucouleurs. Wie sich später zeigen sollte, war Vaucouleurs' Annahme des atmosphärischen Drucks mit 100 hPa viel zu hoch. Mariner 4 startete im Jahr 1964 zum Mars, wobei die Frage nach dem Oberflächendruck weiterhin ungelöst war. Mit Hilfe der neuen Messtechnik konnte ein Druck an der Oberfläche des Mars gemessen werden, der etwa zwei Größenordnungen unter dem damaligen Richtwert lag.

Messprozess

Beim Radio-Okkultationsexperiment durchleuchtet eine hochstabile kontinuierliche Welle im Mikrometerbereich die zu untersuchende planetare Atmosphäre. Gesendet wird im sogenannten X-Band bei 8,4 GHz und im sogenannten S-Band bei 2,3 GHz. Bei genauer Kenntnis aller beteiligten Geschwindigkeitskomponenten können aufgrund der Phasenstabilität der Referenzsignale Frequenzverschiebungen dem Brechungsindexprofil der Atmosphäre zugeordnet werden. Hieraus lassen sich in einem ersten Schritt entsprechende Beugungswinkel und Strahlparameter ermitteln. Der Beugungswinkel kennzeichnet die Ablenkung, die ein Radiostrahl in der Atmosphäre unterliegt. Zu jedem Beugungswinkel gibt es genau einen Strahlparameter, der senkrecht auf den Strahlasymptoten steht (Abbildung 2). In einem weiteren Schritt lässt sich mittels einer Abel-Transformation der zugehörige Brechungsindex ermitteln.

Gemessen wird mit hoher Präzision das vom Raumfahrzeug entsandte Signal, nämlich die Verschiebung der empfangenen Frequenz, relativ zur Sendefrequenz zu jedem Zeitpunkt. Die Verschiebungen in den gemessenen Frequenzen rühren aus den Geschwindigkeiten des Senders und des Empfängers. Diesem klassischen Dopplereffekt ist eine zusätzliche Verschiebung durch das durchleuchtete Medium aufgeprägt. Um diesen Anteil der Frequenzverschiebung aus den gemessenen Frequenzdaten zu extrahieren, subtrahiert man von den gemessenen Daten die Modellfrequenz, die aus den bekannten Geschwindigkeiten des Senders und Empfängers ermittelt werden können und keine atmosphärischen Anteile enthält. Dieser Anteil der Frequenzverschiebung wird als Residuenfrequenz bezeichnet und enthält im Idealfall lediglich Anteile des sondierten Mediums. Die Modellfrequenz berücksichtigt Geschwindigkeitsanteile aus der Rotation, Nutation und Präzession der Erde, idealerweise auch aus den Gezeiteneffekten sowie Anteile aus spezieller und allgemeiner Relativitätstheorie. Des Weiteren sind zur Trennung der Frequenzänderung der Trägerwelle durch die Planetenatmosphäre von denen der Erdatmosphäre, die Effekte der Erdatmosphäre und Erdionosphäre durch Modelle und empirisch ermittelte Werte wie Temperatur, Druck und Luftfeuchtigkeit zu berücksichtigen.

Abbildung 2: Zur Dopplerverschiebung. Dargestellt sind außer dem Beugungswinkel α und den Strahlparametern a, der Geschwindigkeitsvektor der Bodenstation zum Zeitpunkt der Messung, der Geschwindigkeitsvektor der Raumsonde zum Übertragungszeitpunkt und die entsprechenden Einheitsvektoren der Strahlasymptoten. Zum kinematischen Effekt der Frequenzverschiebung tragen lediglich die gezeigten Geschwindigkeiten bei, welche Projektionen der Geschwindigkeitsvektoren in die Okkultationsebene sind.

Es ist bekannt, dass eine elektromagnetische Welle in einem Medium gebrochen und phasenverschoben wird. Gegenüber dem Vakuumssignal legt es einen längeren Weg zurück. Dies lässt sich durch folgende Formel ausdrücken:

Die atmosphärische Frequenzverschiebung entspricht der Zeitableitung von , so dass bei dem implizit zeitabhängigen Sondierungsexperiment eine atmosphärische Frequenzverschiebung auftritt:

mit

f ist die Sendefrequenz und c die Vakuumlichtgeschwindigkeit.

In Verbindung mit der Geometrie der Radiosondierungsmessung und den beteiligten Geschwindigkeitsgrößen lässt sich obige Gleichung für die Frequenzverschiebung mit den beteiligten kinematischen Größen aus Abbildung 2 in Verbindung setzen:

Die atmosphärische Frequenzverschiebung ist also der Beitrag durch Projektion der Geschwindigkeiten auf die Strahlasymptoten, minus den kinematischen Dopplerbeitrag entlang der Sichtlinie. Hat man aus einer Messung erhalten, so muss obige Gleichung für die Komponenten gelöst werden. Die Annahme sphärischer Symmetrie liefert mit dem Brechungsgesetz von Snellius an Kugelflächen (Benndorff-Satz oder Bouguer-Satz) eine zweite Gleichung. Mit den beiden Gleichungen lassen sich die unbekannten Strahlasymptoten bzw. und bestimmen, was zu den Beugungswinkel und Strahlparametern führt.

Koordinatentransformation

Abbildung 3: Zur Okkultationsebene

Abbildung 4: Zur Koordinatentransformation

Bei Annahme sphärischer Symmetrie lässt sich zeigen, dass sich der Radiostrahl stets in einer Ebene in Richtung steigender Brechungsindex bewegt. Diese Ebene, im folgenden Okkultationsebene genannt, ist das Bezugssystem für die Berechnung der Beugungswinkel. Eine solche Ebene zu einem bestimmten Zeitpunkt, wird durch die drei Punkte Bodenstation, Planetenzentrum und Satellitenposition definiert, wie in Abbildung 3 dargestellt ist. Im Allgemeinen liegt zu jedem Zeitpunkt der Messung eine neue Okkultationsebene vor, die zu jedem Messzeitpunkt berechnet werden muss.

Sämtliche Orbitdaten der Raumsonde liegen ebenso wie die Ephemeriden in einem bestimmten Koordinatensystem (KS) vor. Üblicherweise handelt es sich hierbei um die Ebene des mittleren Erdäquators, bezüglich des dynamischen Äquinoktiums der Epoche J2000. Um in das Koordinatensystem der Okkultationsebene zu wechseln sind zwei Koordinatentransformationen nötig. Erstens in das planetozentrische Koordinatensystem, welches analog zum obigen Erdäquator-System ist und sodann in die Okkultationsebene. Abbildung 4 gibt dazu eine Übersicht. Andere Vorgehensweisen sind möglich.

Die Berechnung lassen sich zum Beispiel mit der Ephemeridenbibliothek SPICE durchführen. Die Spice-Bibliothek wird von der NAIF-Gruppe (NASA) frei zur Verfügung gestellt und liegt im Quelltext vor. Die NAIF-Gruppe liefert zudem die entsprechenden Orbitdaten im Spice-Format.

Beugungswinkel und Strahlparameter

Abbildung 5: Zur Berechnung der Beugungswinkel

Die Okkultationsebene, wie sie in Abbildung 5 dargestellt ist, wird aufgespannt durch die beiden Achsen z und r. Die Komponenten der einzelnen Geschwindigkeitsvektoren lassen sich in Termen dieser Achsen angeben. Die Unbekannten Richtungen der Einheitsvektoren und der Strahlasymptoten sind spezifiziert durch die dargestellten Winkel zwischen Strahlasymptote an der Bodenstation und Sichtlinie, sowie durch den Winkel zwischen Strahlasymptote an der Raumsonde und Sichtlinie. Des Weiteren sind die Winkel zwischen den Koordinatenachsen und der Sichtlinie dargestellt; gibt den Winkel zwischen z-Achse und Sichtlinie an, ist der Winkel zwischen r-Achse und Sichtlinie.

Die beiden unbekannten Winkel und lassen sich durch ein Gleichungssystem iterativ lösen. Das Gleichungssystem ist ausführlich in dargestellt. Aufsummierung der beiden unbekannten Winkel ergibt gerade den Beugungswinkel . Der Strahlparameter ergibt sich dann ebenfalls aus einem der unbekannten Winkel und dem Snellius-Gesetz an Kugelflächen:

Wurden die Beugungswinkel und Strahlparameter für jeden Strahl bestimmt, so kann mittels einer Abel-Transformation der zu jedem Beugungswinkel zugehörige Brechungsindex berechnet werden.

Abel-Transformation und Refraktivität

Abbildung 6: Zur Abel-Transformation

Abbildung 7: Schematische Darstellung der Sondierungsmessung. Zu jedem Abtastwert t_i erfolgt für jede Schicht 'i' aus der Messung von Δ f_i die Bestimmung der zugehörigen Strahlparameter a_i und Beugungswinkel α_i, woraus durch inverse Abel-Transformation der Brechungsindex n_i bestimmt wird.

Unter einer Abel-Transformation versteht man die Lösung einer Abelschen Integralgleichung. Integralgleichungen sind dadurch gekennzeichnet, dass die unbekannte Funktion als Integrand in einem bestimmten Integral vorkommt. Bekannte Integralgleichungen sind zum Beispiel die Fourier-Transformation und die Hankel-Transformation. In der Theorie der Integralgleichungen wird die Abelsche Integralgleichung den schwach singulären Volterraschen Integralgleichungen der ersten Art zugeordnet, zu der Niels Henrik Abel 1826 eine Lösung fand.

Im Allgemeinen führen Problemstellungen, bei welchen die Rekonstruktion zweidimensionaler radialsymmetrischer Verteilungen f(r) aus ihren Projektionen g(y) gefordert wird, zur Abel-Transformation. Sie findet zum Beispiel Anwendung in der Plasmaphysik, Astrophysik oder der medizinischen Computertomographie. Der Zusammenhang zwischen der radialen Verteilung f(r) und ihren Projektionen g(y), führt zur Gleichung:

wobei die Abel-Transformation bezeichnet. Eine gebräuchliche Form der inversen Abel-Transformation ist gegeben durch:

Der gemessene Beugungswinkel entspricht der Abel-Transformation eines radialsymmetrischen Brechungsindex. Eine inverse Abel-Transformation des Winkels ergibt also den Brechungsindex angegeben werden kann. Fjelbo u. a. konnten zeigen, dass eine Atmosphäre mit dem Brechungsindexprofil einen Strahl mit Strahlparameter um folgenden Betrag krümmt:

An der Strahlperiapsis, dem Punkt der nächsten Annäherung des Strahls am Planeten, ist das Produkt von Brechungsindex und Radius gleich , wobei durch die Inverse Abel-Transformation bestimmt ist:

Der Brechungsindex ist für Luft eine Zahl nahe bei 1. Um den Umgang mit dieser Zahl handlicher zu gestalten, subtrahiert man eine 1 und multipliziert mit , dies ist die Refraktivität N

wobei der Brechungsindex eine Funktion der Höhe h ist. Die Refraktivität der Neutralatmosphäre ist das Vermögen des Mediums auf die Kraftwirkung der elektromagnetischen Welle zu reagieren, was in der elektromagnetischen Theorie als dielektrische Polarisation bezeichnet wird. Die Polarisation ist als Mittelung über viele Teilchen definiert. Daraus ergibt sich eine Proportionalität der Refraktivität zur Teilchendichte. Die Teilchendichte unterliegt der barometrischen Höhenverteilung, so dass die Refraktivität ebenfalls eine Funktion der Höhe ist. Die Änderung in der Refraktivität, die sich in einer bestimmten Höhe einstellt, ist abhängig von den thermodynamischen Grundgrößen Druck und Temperatur in dieser Höhe, ferner wirken ionisierte Teilchen direkt auf das Signal.

Der korrespondierende Radius der Strahlperiapsis ist . Er kennzeichnet den Punkt der nächsten Annäherung des Strahls am Planeten, dies ist in Abbildung 6 mit dargestellt. Somit liegen nach der Abel-Transformation der Beugungswinkel Profile der Refraktivität als Funktion der Höhe vor. Mit diesen Informationen lassen sich Höhenprofile der charakteristischen Merkmale der planetaren Atmosphäre sehr genau bestimmen.

Temperatur-, Druck- und Dichteprofile

Mit den Ergebnissen aus dem vorhergehenden Abschnitt lassen sich nun die fundamentalen Zustandsgrößen der Neutral-Atmosphäre berechnen. Aus der Proportionalität der Refraktivität zur Anzahldichte, dem idealen Gasgesetz und der hydrostatischen Grundgleichung lassen sich Druck, Temperatur und Dichte berechnen. Ferner lässt sich das Elektronendichteprofil der Ionosphäre berechnen.

Temperaturprofil der Marsatmosphäre aufgenommen mit Mars Express. Die operationelle und wissenschaftliche Leitung des Radio-Okkultationsexperiments von Mars Express (Teil des Mars Express Orbiter Radio Science-Experiments) liegt bei Dr. Martin Pätzold, vom Institut für Geophysik und Meteorologie (IGM) der Universität zu Köln.

Druckprofil der Marsatmosphäre aufgenommen mit Mars Express

Dichteprofil der Neutralatmosphäre des Mars aufgenommen mit Mars Express

• Anzahldichte:

• Hydrostatische Grundgleichung:

• Ideales Gasgesetz:

• Elektronendichte:

Beispiele

Titans Smogschicht

Einzelnachweise

1. G. L. Tyler, B. Ahmad: Radio Occultation. Stanford University California, 2002.  (Buchentwurf)
2. A. Pannekoek: Über die Erscheinungen, welche bei einer Sternbedeckung durch einen Planeten auftreten. In: Astronomische Nachrichten. Band 164, 1904, S. 5–10. 
3. V. R. Eshleman, G. Fjeldbo: The Bistatic Radar-Occultation Method for the Study of Planetary Atmospheres. In: Journal of Geophysical Research. Vol. 70 (1965), S. 3217.
4. Yunck u. a.: Special issue of Terrestrial, Atmospheric and Oceanic Science. In: History of GPS Sounding. 11 (2000), S. 1–20.
5. G. de Vaucouleurs: Physics of the Planet Mars. trans. Patrick Moore, Faber and Faber, London 1954.
6. D. D. Morabito, S. W. Asmar: Radio-Science Performance Analysis Software. In: TDA Progress Report. 42, JPL NASA (1995), S. 120ff.
7. B. Ahmad: Accuracy and resolution of atmospheric profiles obtained from radio occultation measurements. Scientific Report No. DPD811-1998-1, Stanford University, 1998.
8. S. Remus: Untersuchungen zur Durchführung von satellitengestützten Radio Science Experimenten im interplanetaren Raum. Dissertation. Universität der Bundeswehr, München 2004.
9. G. A. Hajj u. a.: A technical description of atmospheric sounding by GPS occultation. In: Journal of Atmospheric and Solar-Terrestrial Physics. Vol. 64, Issue 4 (2002), S. 451–469.
10. G. Fjeldbo, A. J. Kliore, V. R. Eshleman: The Neutral Atmosphere of Venus as Studied with the Mariner V Radio Occultation Experiments. In: Astron. J. Vol. 76 No. 2 (1971), S. 123–140.
11. Pretzler u. a.: Comparison of Different Methods of Abel Inversion Using Computer Simulated and Experimental Side-On Data. In: Zeitung für Naturforschung. Vol. 47a (1992), S. 995–970.
12. Jenkins: Variations in the 13 cm Opacity below the Main Cloud Layer in the Atmosphere of Venus Inferred from Pioneer-Venus Radio Occultation Studies 1978–1987. Dissertation. Georgia Institute of Technology, 1992
13. D.P. Hinson u. a.: Initial results from radio occultation measurements with Mars Global Surveyor. In: Journal of Geophysical Research. 104 E11 (1999), S. 26297–27012.
14. M. Pätzold u. a.: MaRS: Mars Express Orbiter Radio Science. MARS EXPRESS: The Scientific Payload herausgegeben von A. Wilson (2004), S. 141–163, ESA Publications Division

Weitere Literatur

• J. Selle: Planung und Simulation von Radio-Science-Experimenten interplanetarer Raumfahrt-Missionen. Dissertation. Universität der Bundeswehr, München 2005.
• J. Wickert: Das CHAMP-Radiookkultationsexperiment: Algorithmen, Prozessierungssystem und erste Ergebnisse. Dissertation. Karl-Franzens-Universität Graz, Scientific Technical Report STR 02/07 (PDF; 11,8 MB), Potsdam, GFZ (2002)

• B. Stanek: Raumfahrt Lexikon. Halwag Verlag, Bern 1983, ISBN 3-444-10288-7, S. 214.
• Axel von Engeln: Satellite based Temperature Profile Determination using Passive Microwave and Radio Occultation Instruments. Logos Verlag, Berlin 2000, ISBN 3-89722-453-4.

Weblinks

• NAIF-Gruppe