Spotlight #41 – Geofizyka i sejsmologia: wnętrze Ziemi, fale sejsmiczne, tektonika płyt i dynamo geomagnetyczne

Geofizyka to fizyka obserwacyjna doprowadzona do skrajności: nie możemy przewiercić się do płaszcza (najgłębszy odwiert, Kolski SG-3, osiągnął zaledwie 12,2 km), więc niemal wszystko, co wiemy o wnętrzu Ziemi, pochodzi z teledetekcji: sejsmologii, grawimetrii, magnetyki i przepływu ciepła. Zgodność między tymi niezależnymi typami danych a Wstępnym Referencyjnym Modelem Ziemi (PREM) to jedno z wielkich osiągnięć nauki XX wieku.

1. Wnętrze Ziemi: model PREM

Wstępny Referencyjny Model Ziemi (Dziewonski i Anderson, 1981) to jednowymiarowy, sferycznie symetryczny model, dopasowany do globalnej średniej obserwacji sejsmologicznych. Mimo ogromnej poprzecznej niejednorodności rzeczywistej Ziemi (uwidocznionej przez tomografię sejsmiczną), PREM dostarcza dokładnych czasów przebiegu dla fal teleseismicznych i pozostaje standardowym punktem odniesienia.

Struktura warstw PREM:
  Warstwa            Głębokość (km) v_P (km/s)  v_S (km/s)  ρ (g/cm³)
  Skorupa górna      0–15          6,0          3,5         2,60
  Skorupa dolna      15–35         6,8          3,9         2,90
  Płaszcz (górny)    35–220        8,1          4,5         3,32
  Strefa przejściowa 410–660       9,5–10,3     5,1–5,6     3,7–3,9
  Płaszcz dolny      660–2891      10,0–13,7    5,6–7,3     4,4–5,6
  Jądro zewn. (ciek) 2891–5150     8,1–10,4        0        9,9–12,2
  Jądro wewn. (stałe)5150–6371    11,3          3,6         12,8–13,1

Główne nieciągłości:
  Nieciągłość Mohorovičicia (Moho): podstawa skorupy na ~35 km (kontynentalna), ~8 km (oceaniczna)
    Zaznaczona gwałtownym przyspieszeniem fal P z 6,8 → 8,1 km/s
  410 km: przejście fazowe oliwin → wadsleyit (egzotermiczne; nachylenie Clapeyrona +3 MPa/K)
  660 km: ringwoodyt → perowskit + ferroperyklaz (endotermiczne; −2 MPa/K)
    Większy skok gęstości; granica termochemiczna częściowo blokująca cyrkulację płaszcza
  Granica jądro-płaszcz (CMB, 2891 km):
    Fale S nie mogą propagować się w ciekłym jądrze zewnętrznym → strefa cienia fal S
    Strefy ultraniskiej prędkości (ULVZ): ogniska częściowego topnienia na CMB, δV_P ≈ −10%
  Nieciągłość Lehmann (5150 km): ciekłe → stałe jądro wewnętrzne
    Granica jądra wewnętrznego: fazy Pkikp kontra Pkjkp badają anizotropię jądra wewnętrznego

Profil ciśnienia i temperatury:
  Ciśnienie: 0 (powierzchnia) → 24 GPa (660 km) → 136 GPa (CMB) → 364 GPa (środek)
  Temperatura: 0°C (powierzchnia) → 1600°C (podstawa litosfery) → 2500°C (CMB)
               → skok do 3600°C (góra jądra zewnętrznego) → ~5400°C (granica jądra wewnętrznego)
               [porównywalne z fotosferą Słońca]
          

2. Fale sejsmiczne i czasy przebiegu

Sejsmiczne fale objętościowe dzielą się na dwie rodziny: fale P (podłużne, kompresyjne, pierwotne), które przemieszczają się przez płyny i ciała stałe, oraz fale S (poprzeczne, ścinające, wtórne), które nie mogą propagować się w ciekłym jądrze zewnętrznym. Fale powierzchniowe — Rayleigha i Love’a — są prowadzone wzdłuż wolnej powierzchni i niosą większość niszczycielskich skutków dużych trzęsień ziemi.

Definicje prędkości fal:
  v_P = √((K + 4G/3)/ρ)    [K = moduł ściśliwości, G = moduł ścinania, ρ = gęstość]
  v_S = √(G/ρ)
  v_P/v_S = √((K/G + 4/3)) ≥ √(4/3) ≈ 1,155 (ciało stałe); → ∞ w płynie (G=0)

Typowe prędkości płaszcza (PREM, głębokość 400 km):
  v_P ≈ 9,1 km/s,  v_S ≈ 5,0 km/s,  v_P/v_S ≈ 1,82
Jądro zewnętrzne (czyste ciekłe żelazo):
  v_P ≈ 8,1–10,4 km/s,  v_S = 0 (brak ścinania → ciecz)

Prawo Snella na zakrzywionych granicach (parametr promienia p = const wzdłuż promienia):
  p = sin(i)/v(z) = r sin(i)/v(r)    (r = odległość radialna)
  Promienie zakrzywiają się z powrotem ku górze, gdy v(z) rośnie z głębokością (typowy przypadek Ziemi)
  Głębokość zawrócenia promienia: p = r_turn/v(r_turn) → rozwiązywane dla r_turn przy danym p

Strefy cienia:
  Strefa cienia fal P: odległość kątowa 103°–142° od epicentrum
    Jądro zewnętrzne silnie załamuje fale P do wewnątrz → luka w bezpośrednich wstąpieniach P
    Faza PKP: fale P przechodzące przez jądro (docierają >142°)
    PKIKP: przenika jądro wewnętrzne; dociera 0°–180° z opóźnieniem czasu przebiegu
  Strefa cienia fal S: 103°–180°
    Fale S nie mogą wejść do ciekłego jądra zewnętrznego → nieobecne powyżej 103°
    Odkrycie ciekłego jądra zewnętrznego: Oldham (1906), Gutenberg (1914)

Fale powierzchniowe (dyspersyjne, głębokość próbkowania zależna od częstotliwości):
  Rayleigha: retrogradny eliptyczny ruch cząstek w płaszczyźnie sagitalnej (sprzężenie SV + P)
    Prędkość grupowa ~0,92 v_S(głębokość ≈ λ/3)
  Love’a: poziomy ruch poprzeczny; wymaga uwarstwienia prędkości dla prowadzenia
    Prędkość fazowa c_L: 2kH tan(k_L H) = (ρ₁v_L1)/(ρ₀v_S0)...  (relacja dyspersji warstwy)
  Tomografia fal powierzchniowych: okres 20–300 s → czułość na głębokość 30–600 km
          

3. Sejsmologia: magnituda, mechanizmy ogniskowe i statystyka

Ilościowe określenie wielkości trzęsienia ziemi wymaga oddzielenia amplitudy fali zarejestrowanej przez odbiornik (funkcji odległości, wzorca promieniowania, tłumienia i odpowiedzi instrumentu) od wewnętrznej siły źródła. Skala magnitudy momentowej osiąga to, wiążąc magnitudę z fizycznym momentem sejsmicznym M₀, który całkuje spadek naprężenia po powierzchni uskoku.

Lokalna magnituda Richtera:
  M_L = log₁₀(A/A₀)    [A = maksymalna amplituda przemieszczenia w µm na 100 km, Wood-Anderson]
  Nasyca się przy M_L ≈ 7 (moment sejsmiczny rośnie szybciej niż amplituda)

Moment sejsmiczny:
  M₀ = μ · Ā · D̄    [μ = moduł ścinania ~30 GPa, Ā = powierzchnia uskoku, D̄ = średnie przemieszczenie]
  Wyznaczany z długookresowych fal objętościowych lub statycznej deformacji GPS

Magnituda momentowa (Hanks i Kanamori, 1979):
  M_w = (2/3)(log₁₀ M₀ − 9,1)    [M₀ w N·m]
  Nie nasyca się; powiązana z uwolnioną energią
  Przykłady:
    Sumatra 2004: M_w = 9,1, M₀ ≈ 4×10²² N·m, powierzchnia uskoku ≈ 1300×200 km²
    Valdivia 1960: M_w = 9,5 (największe zmierzone instrumentalnie; Chile)
    Uderzenie Chicxulub: M_w ≈ 10,5–11 (szacunkowo)

Mechanizm ogniskowy (rozwiązanie płaszczyzny uskoku):
  Tensor momentu M_ij: symetryczny tensor drugiego rzędu opisujący pary sił
  Wartości własne: osie T (rozciąganie), N (zerowa), P (ściskanie)
  „Beach ball": stereograficzna projekcja dolnej półkuli polaryzacji pierwszych wstąpień
    Kwadrantowy wzorzec kompresji (najpierw w dół)/dylatacji (najpierw w górę)
  Typy uskoków:
    Normalny:        skrzydło zwisające w dół (granica rozbieżna; naprężenie rozciągające)
    Nasunięcie/odwrotny: skrzydło zwisające w górę (granica zbieżna; ściskające)
    Przesuwczy:       ruch poziomy (uskok transformujący)

Relacja Gutenberga-Richtera:
  log₁₀ N(M) = a − bM
  b ≈ 1,0 (globalnie; lokalnie w zakresie 0,6–1,3)
  Oznacza: na każde M7,0 przypada ~10 trzęsień M6,0 i ~100 trzęsień M5,0
  Interpretacja fizyczna: bezskalowa, fraktalna geometria uskoków

Zanik wstrząsów wtórnych Omoriego-Utsu:
  n(t) = K/(t + c)^p    [n = częstość wstrząsów wtórnych, K,c,p = stałe; p ≈ 1,0]
  Produktywność wstrząsów wtórnych: N_total ∝ 10^{α(M_main − M_min)}, α ≈ 1
          

4. Tektonika płyt i dynamika płaszcza

Teoria tektoniki płyt, ugruntowana w latach 60. XX wieku (Hess, Vine, Matthews, Wilson, McKenzie, Parker), zjednoczyła wcześniej odrębne dziedziny dryfu kontynentów, spreadingu dna oceanicznego, sejsmiczności i wulkanizmu w jednych, kinematycznych ramach. Silnik napędowy — konwekcja płaszcza — pozostaje aktywnym obszarem badań, ponieważ zależy od reologii skał krzemianowych pod ekstremalnym ciśnieniem i temperaturą.

Reologia płaszcza:
  Lepkość efektywna: η_eff ≈ 10^{21} Pa·s (górny płaszcz, poglacjalne wydźwiganie)
                      10^{22}–10^{23} Pa·s (dolny płaszcz)
  Pełzanie dyslokacyjne + dyfuzyjne: ε̇ = A σⁿ exp(−E_a/RT) (pełzanie potęgowe, n≈3,5)
  Przejście od sprężystego do lepkiego: t_Maxwell = η/G ≈ 10^{21}/70×10^9 ≈ 450 lat

Bezwymiarowe liczby konwekcji płaszcza:
  Rayleigh: Ra = ρ_0 g α ΔT d³/(κη)
    ρ₀ ≈ 4000 kg/m³; g = 10 m/s²; α ≈ 2×10⁻⁵ K⁻¹; ΔT ≈ 2500 K
    d ≈ 2900 km; κ ≈ 10⁻⁶ m²/s; η ≈ 3×10²¹ Pa·s
    Ra ≈ 10⁷ → intensywna konwekcja (próg ~Ra_c = 1100 dla warstwy płaskiej)
  Pokrywa sztywna kontra pokrywa ruchoma: zależy od granicy plastyczności zimnej litosfery

Spreading dna oceanicznego i anomalie magnetyczne:
  Tempo: grzbiet śródatlantycki ~2,5 cm/rok (wolny); Wschodniopacyficzne Wzniesienie ~15 cm/rok (szybki)
  Vine-Matthews-Morley (1963): nowa skorupa oceaniczna zapisuje kierunek pola geomagnetycznego
    → naprzemienne pasy normalnej/odwróconej polaryzacji równoległe do osi grzbietu
  Wiek skorupy oceanicznej: maks. ~200 mln lat (jura); starsza jest subdukowana
  Wiek dna oceanicznego = (odległość od grzbietu) / (tempo spreadingu)

Prędkości płyt GPS (ITRF2020 / NNR-MORVEL56):
  Para płyt                  Tempo (mm/rok)  Azymut
  Pacyficzna–Północnoamerykańska  ~50           NNZ
  Indoaustralijska–Euroazjatycka  ~45           NNW
  Afrykańska–Euroazjatycka        ~7            NNZ
  Północnoamerykańska–Euroazjatycka ~24         WNW (otwieranie Atlantyku)
  Nubijska–Somalijska             ~6            WNW (Ryft Wschodnioafrykański)

Cykl Wilsona (pełne otwieranie i zamykanie oceanu):
  1 – Ryftowanie kontynentalne (Ryft Wschodnioafrykański dziś: wczesny etap)
  2 – Proto-ocean (Morze Czerwone: wąski ocean)
  3 – Otwarty ocean (Atlantyk: dojrzałe pasywne krawędzie kontynentów)
  4 – Inicjacja subdukcji → łuk wyspowy (Nazca–Ameryka Południowa)
  5 – Zamknięcie oceanu + kolizja kontynentalna (Alpy, Himalaje, Appalachy)
  Czas trwania: ~500 mln lat pełny cykl
          

5. Pole grawitacyjne Ziemi i geoida

Geoida to powierzchnia ekwipotencjalna pola grawitacyjnego Ziemi, która (średnio) pokrywa się ze średnim poziomem morza. Odchylenia od elipsoidy odzwierciedlają anomalie gęstości w skorupie i płaszczu. Misja satelitarna GRACE (2002–2017) oraz jej następczyni GRACE-FO (od 2018) mierzą pole grawitacyjne do stopnia harmonicznej sferycznej 120 (rozdzielczość ±330 km), śledząc zmieniającą się odległość między dwoma pojazdami kosmicznymi.

Rozwinięcie pola grawitacyjnego:
  V(r,θ,λ) = (GM/r)∑_l ∑_m (R/r)^l P_l^m(cosθ) [C_lm cos(mλ) + S_lm sin(mλ)]
  EGM2008 (NGA): kompletny do l_max = 2159 (rozdzielczość ~10 km)
  C₂₀ = J₂ (spłaszczenie Ziemi): −1,0826 × 10⁻³  (sferoida spłaszczona)

Poprawki grawitacyjne (redukcja do powierzchni standardowej):
  Poprawka wolnego powietrza: +0,3086 μGal/m w górę (−3,086 μGal/m wysokość obserwacji)
  Poprawka Bouguera: −2πGρh ≈ −0,1119ρh mGal/m  (masa warstwy między obs. a geoidą)
    Pełna poprawka Bouguera: poprawka terenowa Tc usuwa lokalną nierówność topografii

  Anomalia wolnego powietrza: Δg_FA = g_obs − g_normal + FAC         [usuwa efekt wysokości]
  Anomalia Bouguera:  Δg_B  = Δg_FA − 2πGρh + Tc             [usuwa warstwę skorupy]

Izostazja (modele Airy'ego i Pratta):
  Airy (George Biddell Airy, 1855): stała gęstość ρ_c; korzeń góry unosi się
    H_root = h · ρ_c/(ρ_m − ρ_c) ≈ 7h   (h = wysokość góry, nadmiarowa gęstość zanurzona)
    Przewidywanie: ujemna anomalia Bouguera pod górami (obserwowana globalnie)
  Pratt (John Henry Pratt, 1855): zmienna gęstość σ(x), stała głębokość bazowa T_0
    σ(x) · (T_0 + h) = σ₀ · T_0   (kompensacja przez lateralną zmianę gęstości)

Zastosowania GRACE:
  Masa lodu Grenlandii: −270 Gt/rok (trend 2003–2022); Antarktyda: −150 Gt/rok
  Wyczerpywanie wód podziemnych: basen Indusu −24 mm/rok ekwiwalentnej wysokości wody
  Poglacjalne wydźwiganie: Fennoskandia +10 mm/rok wzrostu grawitacji
  Cykl hydrologiczny: sezonowe wahania wysokości geoidy 10–20 mm w dorzeczach Amazonki/Konga
          

6. Pole geomagnetyczne i teoria dynama

Pole magnetyczne Ziemi jest generowane przez turbulentną konwekcję w ciekłym żelaznym jądrze zewnętrznym. Pole jest głównie dipolowe (dipol osiowy ≈ 80% pola powierzchniowego), ale zmienia się w skalach czasu od milisekund (wariacja wiekowa wynikająca z elektromagnetycznie sprzężonej dynamiki jądra) po miliony lat (inwersje polarności). Zapisy paleomagnetyczne w bazalcie dna oceanicznego zachowują historię tych inwersji aż do jury.

Źródła pola geomagnetycznego:
  Pole jądra: 95% całości; generowane przez dynamo MHD; zmienia się w skali lata–mln lat
  Pole skorupowe (litosferyczne): magnetyzm szczątkowy; anomalie ±2000 nT; quasi-statyczne
  Pole zewnętrzne: prądy magnetosfery + jonosfery; wariacja godziny-dni; spokojne Sq ≈ 30 nT

Równanie indukcji MHD:
  ∂B/∂t = ∇×(u×B) − ∇×(η ∇×B)
  η = 1/(μ₀σ) = dyfuzyjność magnetyczna [m²/s]
  Bilans: magnetyczna liczba Reynoldsa Rm = u L/η
  Jądro zewnętrzne Ziemi: u ≈ 0,3 mm/s, L ≈ 2000 km, η ≈ 1 m²/s → Rm ≈ 600 >> 1
  Wysokie Rm: dyfuzja pomijalna → pole „zamrożone" w płynie (twierdzenie Alfvéna)

Model konwekcji kolumnowej Bussego (1975):
  Rotacja planetarna (Ω = 7,3×10⁻⁵ rad/s) → siła Coriolisa organizuje konwekcję w
  cylindryczne walce równoległe do osi obrotu
  Twierdzenie Taylora-Proudmana: w szybko obracających się płynach przepływ jest 2D, prostopadły do Ω
  Walce konwekcyjne generują przepływ helikalny → efekt α (dynamo pola średniego):
    ⟨u×B⟩ = αB + ... (Steenbeck, Krause, Rädler, 1966)
  Efekt Ω: rotacja różnicowa rozciąga pole poloidalne w toroidalne
  Razem: dynamo αΩ → podtrzymuje oba składniki wbrew dyfuzji

IGRF (Międzynarodowe Referencyjne Pole Geomagnetyczne):
  Reprezentacja sferyczno-harmoniczna pola wewnętrznego do stopnia/rzędu 13
  Aktualny moment dipolowy: 7,94×10²² A·m² (spadek o 6% od 1840 roku)
  Anomalia Południowoatlantycka (SAA): region słabego pola; strumień cząstek auroralnych na niskiej wysokości
  Geomagnetyczny biegun północny (2024): ~86°N, 133°W (dryfuje ku Syberii ~40 km/rok)

Paleomagnetyczne inwersje:
  Geomagnetyczna Skala Czasu Polarności (GPTS) z oceanicznych anomalii magnetycznych + biostratygrafii
  Obecnie: normalny chron Brunhes (od 780 tys. lat temu)
  Ostatnia inwersja: Matuyama-Brunhes, 781 tys. lat temu (czas trwania przejścia 10–20 tys. lat)
  Średnia częstość inwersji: 4–5 na mln lat (kenozoik); kredowy superchron normalny: 40 mln lat bez inwersji
  Zapowiedzi: eskursje, minima natężenia, wielobiegunowa struktura pola podczas przejścia
  Potencjalne skutki inwersji: 10-krotne rozszerzenie SAA, zorza na średnich szerokościach geograficznych
          

Wypróbuj te symulacje