Stanford University · Gastforscher · 2025 – 2026

Neue Methoden für Eisradar

In den untersten paar hundert Metern des antarktischen Eises sieht Radar fast nichts, deshalb wird die Bewegung dort meist extrapoliert. Ein Jahr Radardaten über einem subglazialen See zeigt, dass diese Zone dennoch eine kohärente Phase trägt, und zwei neue Schätzverfahren machen daraus Geschwindigkeiten bis zur Eisbasis.

OrtMercer Subglacial Lake, Westantarktis
Aufzeichnung312 Tage · 1'878 Radarmessungen
InstrumentApRES-FMCW-Radar, 200–400 MHz
Eisdicke1'094 m

Basierend auf Evidence for Coherent Phase Signal in the Antarctic Echo-Free Zone von H. M. Stählin, D. M. Schroeder und M. R. Siegfried.

01

Die Messung

Ein Radar auf der Eisoberfläche wiederholt fast ein Jahr lang dieselbe Messung. Alles hier beruht darauf, wie sich seine Echos verändern.

bewegt sich mit dem Eis · 226 m/Jahr ApRES 0 m 600 m 1'094 m Interne Schichten hell, durchgehend Echos Schichtloses Eis 20–30 dB schwächer, keine durchgehenden Schichten Subglazialer See heller, flacher Reflektor
  1. 01

    Alle vier Stunden eine Radarmessung

    Eine Frequenzrampe von 200 bis 400 MHz geht ins Eis. Die Laufzeit des Echos ergibt die Tiefe in 5-cm-Schritten, bis hinunter zum See auf 1'094 m.

  2. 02

    Amplitude und Phase in jeder Tiefe

    Jede Tiefe liefert eine Amplitude (Echostärke) und eine Phase (Lage im Wellenzyklus). 1 mm Bewegung verschiebt die Phase um 1.3°.

  3. 03

    312 Tage Wiederholungen

    1'878 Messungen nebeneinander ergeben das Radargramm. Bewegung nach unten lässt die Phase driften; seitliche Bewegung verwürfelt das Echomuster allmählich.

Schichten lassen sich leicht verfolgen. Das unterste Eis zeigt keine Schichten, nur schwache Echos kleiner Streuer. Zeichnet deren Phase die Bewegung trotzdem auf?

02

Die Eissäule

312 Tage Echos und drei daraus berechnete Profile auf einer gemeinsamen Tiefenachse. Mit der Maus über eine Tiefe fahren, um sie abzulesen.

Radargramm

Kontrast nach Tiefe ausgeglichen. Schwache Echos bleiben sichtbar.

Kohärenz |γ|

Phasenstabilität von Messung zu Messung, 0 bis 1

Vertikale Geschwindigkeit

CW-MLPRBlockbewegung

Horizontale Dekorrelation

MDI P(v)Median

Tiefe-
Zone-
Kohärenz-
Vertikal (CW-MLPR)-
Modell der Blockbewegung-
Horizontal (MDI-Median)-

Das schichtlose Eis ist 20 bis 30 dB schwächer als die Schichten und wirkt deshalb im echten Massstab leer. Ausgeglichen werden seine gesprenkelten Echos sichtbar. In der Amplitude gleichen sie Rauschen; die Phase unterscheidet sie davon.

Geschwindigkeit ist nach unten positiv. Die Linie der Blockbewegung ist nur oberhalb von 600 m angepasst, ihre Übereinstimmung im schichtlosen Eis ist also eine unabhängige Prüfung. Rauten: herkömmliche Schichtverfolgung. Horizontale Werte sind relativ zum Instrument, daher sind Verhältnisse zwischen Tiefen verlässlicher als absolute Werte.

03

Signal im schichtlosen Eis

Ist die Phase im schichtlosen Eis mehr als Rauschen? Schritt für Schritt durch die zwei Korrekturen, mit dem Rauschen unterhalb der Eisbasis als Referenz.

Tiefenfenster
Amplitude Amplitude des gewählten Tiefenfensters über 312 Tage
Tag 0Tag 312
Phase Phase des gewählten Tiefenfensters über 312 Tage
Tag 0Tag 312

04

Zwei Schätzverfahren

Beide arbeiten auf dem komplexen Signal ohne Peak-Verfolgung und wenden jeden nichtlinearen Schritt erst nach der zeitlichen Mittelung an, sodass sich das Rauschen zuerst aufhebt.

CW-MLPR · vertikale Geschwindigkeit

Coherence-Weighted Multi-Lag Phase Regression

Jede Tiefe wird mit sich selbst bei Verzögerungen von 1 bis 8 Messungen korreliert. Die Phase wächst linear mit der Verzögerung, mit einer Rate, die durch die Geschwindigkeit bestimmt ist. Eine gewichtete Gerade durch den Ursprung liefert so die Geschwindigkeit, ohne dass die 2π-Sprünge der Phase korrigiert werden müssen. Mit nur einer Verzögerung entspricht sie dem Standard-Schätzer für ApRES.

v̂ = (λ/4π) · Σ |R(ℓ)| Φ(ℓ) τℓ / Σ |R(ℓ)| τℓ²
Phase jeder Verzögerung (Grösse = Gewicht) CW-MLPR-Anpassung wahre Geschwindigkeit
Wahr-m/Jahr
CW-MLPR-RMSE -
2π-Sprünge korrigiert + Gerade-RMSE -

Simulierte Aufzeichnung mit 1'878 Messungen. Der klassische Weg (2π-Sprünge der Phase korrigieren, dann eine Gerade anpassen) ist bei sauberen Daten präzise, versagt aber unterhalb eines SNR von etwa 2. RMSE über 40 Durchläufe.

MDI · Verteilung der horizontalen Geschwindigkeit

Multi-band Decorrelation Inversion

Streuer, die seitlich durch den Strahl driften, verlieren ihre Kohärenz mit einer Rate, die von ihrer Geschwindigkeit abhängt. MDI misst dies in drei Teilbändern und invertiert sie gemeinsam zur vollständigen Geschwindigkeitsverteilung.

|γ(τ; λ)| = ∫ P(v) · exp(−4π²σθ² v² τ² / λ²) dv
Kohärenz nach Teilband
Bestimmte Verteilung
tiefes volles hohes Band wahre Populationen MDI

Vorwärtsmodell mit den Einstellungen der echten Daten (σθ = 0.22 rad, drei Teilbänder, Verzögerungen bis 104 Tage), invertiert mit der nichtnegativen Methode der kleinsten Quadrate (NNLS). Eine Anpassung mit nur einer Geschwindigkeit (rot) landet zwischen den Populationen.