Wie sieht der J₂₂-Moment der Erde aus?

Das Geopotenzialmodell von Wikipedia gibt die Form der verknüpften Potenzialkomponente an J 2 als

u = J 2 P 2 0 ( Sünde θ ) R 3 = J 2 1 R 3 1 2 ( 3 Sünde 2 θ 1 ) = J 2 1 R 5 1 2 ( 3 z 2 R 2 )

und die Beschleunigungskomponenten (Gradient dieses Potentials) sind:

A X = J 2 X | R | 7 ( 6 z 2 1.5 ( X 2 + j 2 ) )

A j = J 2 j | R | 7 ( 6 z 2 1.5 ( X 2 + j 2 ) )

A z = J 2 z | R | 7 ( 3 z 2 4.5 ( X 2 + j 2 ) ) .

Aus dem Geopotentialmodell und dieser Antwort .

Ich habe geplant u Und | A | unten für die X z Flugzeug, wo z ist die Rotationsachse der Erde. Die Verteilung in drei Dimensionen wird einfach um diese gedreht z Achse.

Frage: J 2 ist der dominierende Teil des Quadrupolmoments der Erde, aber es ist nur eine Komponente.

  1. Tut J 22 den Rest des Quadrupolmoments der Erde darstellen?
  2. Was macht J 22 aussehen? Wie ist seine Form und Symmetrie?

Potential und Beschleunigung aufgrund des J2 der Erde

import numpy as np
import matplotlib.pyplot as plt

Re = 6378137. # meters
J2 = 1.7555E+25 # m^5/s^2
hw = 400
x = np.linspace(-2.5*Re, 2.5*Re, 2*hw+1)
x, z = np.meshgrid(x, x)
y = np.zeros_like(x)
x2, y2, z2 = [w**2 for w in (x, y, z)]
r2 = x2 + y2 + z2
r = np.sqrt(r2)
r7 = r**7

u = J2 * r**-5 * 0.5 * (3*z2 - r**2)
ax = J2 * (x/r7) * (6*z2 - 1.5*(x2+y2)) 
ay = J2 * (y/r7) * (6*z2 - 1.5*(x2+y2)) 
az = J2 * (z/r7) * (3*z2 - 4.5*(x2+y2)) 

for thing in (u, ax, ay, az):
    thing[r<Re] = np.nan

axz = np.sqrt(ax**2 + az**2)

if True:
    plt.figure()
    plt.subplot(1, 2, 1)
    umax = np.nanmax(np.abs(u))
    plt.imshow(u, cmap='seismic', vmin=-umax, vmax=umax)
    plt.colorbar()
    plt.title('specific potential energy (J/kg)', fontsize=16)
    plt.plot([hw, hw], [0.3*hw, 1.7*hw], '-k')
    plt.xlim(0, 2*hw+1)
    plt.ylim(0, 2*hw+1)
    plt.subplot(1, 2, 2)
    plt.imshow(axz, cmap='inferno')
    plt.colorbar()
    plt.title('acceleration magnitude (m/s^2)', fontsize=16)
    plt.plot([hw, hw], [0.3*hw, 1.7*hw], '-k')
    plt.xlim(0, 2*hw+1)
    plt.ylim(0, 2*hw+1)
    plt.suptitle("Earth's geopotential's J2 component", fontsize=16)
    plt.show()
@gerrit Ich werde gebührend ermahnt! Ich habe dies das erste Mal i.stack.imgur.com/1XdmH.png aus dem Standard von Python 3's matplotlib erhalten, also habe ich schnell mein altes Python 2 verwendet, das standardmäßig auf jet. Ich erinnere mich, dass ich einen Tweet mit einem alten Haiku über Jet gelesen habe, aber ich kann es jetzt nicht finden. Ich habe "Don't use Jet!" in diesem und diesem Plot und entschuldigte sich dafür auch in einem von ihnen. Ich werde nach einer besseren cmap suchen und updaten ;-)
@gerrit irgendwie denke ich, dass es jetzt ablenkender aussieht, aber ich habe die Farbkarten aktualisiert.
Ich empfehle colorbrewer2 :)

Antworten (2)

Repräsentiert J 22 den Rest des Quadrupolmoments der Erde?

Mit einer gebräuchlicheren Schreibweise l Grad , m Ordnung gibt es die zonalen l 2 m 0 , tesseralen l 2 m 1 und sektoralen l 2 m 2 sphärischen Harmonischen, aber Sie kennen sie wahrscheinlich bereits.

Wie sieht J22 aus ? Wie ist seine Form und Symmetrie?

Hier ist, was ich für l 2 m 2 bekomme :

Geben Sie hier die Bildbeschreibung ein

Zum Vergleich, hier ist, was ich für Ihr Diagramm "Beschleunigungsgröße" (l 2 m 0 ) bekomme (es ist nicht klar, was 800 darstellt):

Geben Sie hier die Bildbeschreibung ein

Vielen Dank für Ihre Antwort! Um es noch einmal zu überprüfen, die Antwort auf "Does J 22 den Rest des Quadrupolmoments der Erde darstellen?" ist Nein? Es gibt ein J 21 sowie? Tut ICH Und J dasselbe darstellen?
Informell betrachtet betrachten wir in der Geodäsie einen Quadrupol als etwas, das auf einer Ebene wirkt (gravitativer, magnetischer, elektrischer Quadrupol). Um das Quadrupolmoment der Erde darzustellen, benötigen Sie daher eine unendliche Anzahl von sphärischen Harmonischen; also meine antwort ist "nein".

Wikipedia geht tatsächlich detailliert auf die Komponenten des Gravitationskraftfeldes der Erde ein.

In der Wikipedia-Referenz haben wir nicht J M , N für M > N > 0 . Vielmehr wird dieser Term in eine Cosinus-Komponente zerlegt C N M und eine Sinuskomponente S N M . Diese Tabelle enthält Ergebnisse von einem aus einer Reihe von Modellen, die als JGM-Modelle bezeichnet werden und vom Goddard Space Center entwickelt wurden:

Zonale Koeffizienten

  N 2 0,1082635854 D 02 3 0,2532435346 D 05 4 0,1619331205 D 05 5 0,2277161016 D 06 6 0,5396484906 D 06 7 0,3513684422 D 06 8 0,2025187152 D 06

Tesserale Koeffizienten [sic; Fälle mit N = M sind sektoral]

N M C S 2 1 0,3504890360 D 09 0,1635406077 D 08 2 2 0,1574536043 D 05 0,9038680729 D 06 3 1 0,2192798802 D 05 0,2680118938 D 06 3 2 0,3090160446 D 06 0,2114023978 D 06 3 3 0,1005588574 D 06 0,1972013239 D 06 4 1 0,5087253036 D 06 0,4494599352 D 06 4 2 0,7841223074 D 07 0,1481554569 D 06 4 3 0,5921574319 D 07 0,1201129183 D 07 4 4 0,3982395740 D 08 0,6525605810 D 08

Hier D bedeutet Multiplizieren mit einer Zehnerpotenz, deren Exponent also folgt D 05 bedeutet multiplizieren mit 10 5 . Die absoluten Werte der obigen Zahlen stellen die relative Stärke des Gravitationspotentials dar, das sich aus der sphärischen harmonischen Komponente ergibt, in einem Abstand, der gleich dem durchschnittlichen Erdradius ist, und mit der isotropen Komponente "sphärische Masse" (was wäre J 0 in dieser Nomenklatur) skaliert auf 1.0 .

Nehmen wir nun an, Sie bringen einen Satelliten in eine niedrige Erdumlaufbahn, in einem radialen Abstand, der im Wesentlichen gleich der oben verwendeten Referenz ist. Unter den Quadrupolkomponenten eindeutig J 2 ist am größten, mit C 2 1 Und S 2 1 viele Größenordnungen kleiner, und C 2 2 Und S 2 2 dazwischen fallen. Sie können den Quadrupol daher vertretbar nur mit dem beschreiben ( 2 , 0 ) ( J 2 ) Und ( 2 , 2 ) ( C 2 2 , S 2 2 ) Komponenten, wobei die sehr kleinen ignoriert werden ( 2 , 1 ) Stücke.

Aber das ist nicht die ganze Geschichte. Während die verbleibende Quadrupolkomponente in der Nähe der Erdoberfläche winzig ist, einige polare Komponenten höherer Ordnung ( N 3 ) sind nicht. Wenn Sie die einschließen ( 2 , 2 ) Komponenten des Quadrupols zur Modellierung Ihrer erdnahen Bahnbewegung sollten Sie auch mindestens die ähnlich großen einbeziehen ( 3 , 0 ) , ( 3 , 1 ) Und ( 4 , 0 ) Begriffe, um Ihre Genauigkeit gegenüber dem rotationssymmetrischen Quadrupol wirklich zu verbessern. Die Abweichungen von einem rituell symmetrischen Quadrupol sind komplexer als nur eine zusätzliche Quadrupolkomponente.

Der Grund dafür ist, dass die Abweichungen von einem rotationssymmetrischen Quadrupol durch die Verteilung von großräumigen geografischen Merkmalen wie Gebirgszügen, Becken und tektonischen Subduktionsgrenzen bestimmt werden. Solche Merkmale auf der Erde sind zu komplex, um von einer oder zwei sphärischen harmonischen Komponenten erfasst zu werden. Wenn Sie die geologische Komplexität der Erde erfassen müssen, um Ihren Satelliten zu modellieren, benötigen Sie dazu viele sphärische harmonische Terme.

Sie können die Dinge vereinfachen, wenn Sie stattdessen Ihren Satelliten in eine große Höhe bringen, sagen wir eine geostationäre Umlaufbahn. Die Komponenten mit größer N Werte verblassen schneller mit der Entfernung als solche mit kleineren Werten, sodass Sie diese je nach der von Ihnen benötigten Genauigkeit möglicherweise auslassen ( 3 , M ) Und ( 4 , 0 ) Komponenten. Die geografischen Warzen, die wir aus der Nähe sehen, werden durch eine weiter entfernte Betrachtung geglättet. Das ist, wo Sie am ehesten Erfolg haben werden J 2 , C 2 2 , Und S 2 2 .

Die Erde ist nicht der einzige Körper mit solchen Gravitationskomplexitäten. Jeder scheinbar hydrostatisch ausgeglichene, aber solide Körper wird wahrscheinlich seine eigenen geografischen Unebenheiten haben. Insbesondere Mars hat den großen Höhenunterschied zwischen Nord und Süd, das Hellas-Becken und Olympus Mons. Wenn wir einen Mars-Orbiter mit hoher Genauigkeit verfolgen wollen, sollten wir besser unsere sphärischen Harmonischen auffrischen.

Okay, um zu überprüfen, ob ich es verstanden habe, 1) J 22 ist nur eine Größe und als einzelne Zahl geschrieben, fehlen genügend Informationen, um die Variation mit dem Längengrad anzuzeigen. 2) dito für J 21 aber es ist sehr klein und selten wird von ihm gesprochen.
Uh-huh, das ist es.
@OscarLanzi Sie haben n = m Terme in Ihre Tabelle "Tesseral-Koeffizienten" eingefügt, aber wenn n = m ist, wird dieser Koeffizient als "sektoral" bezeichnet.
Die Tabelle ist aus dem Quellartikel zitiert. Siehe die Anmerkung, die ich hinzugefügt habe.