Спотлайт #41 – Геофізика та сейсмологія: надра Землі, сейсмічні хвилі, тектоніка плит та геомагнітне динамо

Геофізика — це спостережна фізика, доведена до крайності: ми не можемо пробурити мантію (найглибша свердловина, Кольська СГ-3, досягла лише 12,2 км), тому майже все, що ми знаємо про надра Землі, походить з дистанційного зондування: сейсмології, гравіметрії, магнетизму та теплового потоку. Узгодженість між цими незалежними типами даних та Попередньою еталонною моделлю Землі (PREM) — одне з найвизначніших досягнень науки XX століття.

1. Надра Землі: модель PREM

Попередня еталонна модель Землі (Дзєвонський і Андерсон, 1981) — це 1D сферично симетрична модель, що узгоджується з глобальним усередненням сейсмологічних спостережень. Незважаючи на величезну латеральну неоднорідність реальної Землі (яку показує сейсмічна томографія), PREM забезпечує точний час пробігу для телесейсмічних фаз і залишається стандартним еталоном.

Структура шарів PREM:
  Шар                Глибина (км)  v_P (км/с)  v_S (км/с)  ρ (г/см³)
  Верхня кора        0–15          6,0          3,5         2,60
  Нижня кора         15–35         6,8          3,9         2,90
  Мантія (верхня)    35–220        8,1          4,5         3,32
  Перехідна зона     410–660       9,5–10,3     5,1–5,6     3,7–3,9
  Нижня мантія       660–2891      10,0–13,7    5,6–7,3     4,4–5,6
  Зовнішнє ядро (рід)2891–5150     8,1–10,4        0        9,9–12,2
  Внутрішнє ядро(тв) 5150–6371    11,3          3,6         12,8–13,1

Основні розриви:
  Розрив Мохоровичича (Мохо): основа кори на ~35 км (континентальна), ~8 км (океанічна)
    Позначений різким прискоренням P-хвиль з 6,8 → 8,1 км/с
  410 км: фазовий перехід олівін → вадслеїт (екзотермічний; нахил Клапейрона +3 МПа/К)
  660 км: рінгвудит → перовскіт + феропериклаз (ендотермічний; −2 МПа/К)
    Більший стрибок густини; термохімічна межа частково блокує циркуляцію мантії
  Межа ядро-мантія (CMB, 2891 км):
    S-хвилі не можуть поширюватись у рідкому зовнішньому ядрі → тіньова зона S-хвиль
    Зони наднизької швидкості (ULVZ): ділянки часткового плавлення на CMB, δV_P ≈ −10%
  Розрив Леманн (5150 км): рідке → тверде внутрішнє ядро
    Межа внутрішнього ядра: фази Pkikp проти Pkjkp досліджують анізотропію внутрішнього ядра

Профіль тиску та температури:
  Тиск: 0 (поверхня) → 24 ГПа (660 км) → 136 ГПа (CMB) → 364 ГПа (центр)
  Температура: 0°C (поверхня) → 1600°C (основа літосфери) → 2500°C (CMB)
               → стрибок до 3600°C (верх зовнішнього ядра) → ~5400°C (межа внутрішнього ядра)
               [порівнянно з фотосферою Сонця]
          

2. Сейсмічні хвилі та час пробігу

Об'ємні сейсмічні хвилі поділяються на дві родини: P-хвилі (первинні, компресійні, поздовжні), що поширюються крізь рідини та тверді тіла, та S-хвилі (вторинні, зсувні, поперечні), що не можуть поширюватись у рідкому зовнішньому ядрі. Поверхневі хвилі — Релея та Лява — спрямовуються вздовж вільної поверхні і несуть більшість руйнівних наслідків від великих землетрусів.

Визначення швидкості хвиль:
  v_P = √((K + 4G/3)/ρ)    [K = модуль об'ємної пружності, G = модуль зсуву, ρ = густина]
  v_S = √(G/ρ)
  v_P/v_S = √((K/G + 4/3)) ≥ √(4/3) ≈ 1,155 (тверде тіло); → ∞ у рідині (G=0)

Типові швидкості мантії (PREM, глибина 400 км):
  v_P ≈ 9,1 км/с,  v_S ≈ 5,0 км/с,  v_P/v_S ≈ 1,82
Зовнішнє ядро (чисте рідке залізо):
  v_P ≈ 8,1–10,4 км/с,  v_S = 0 (немає зсуву → рідина)

Закон Снелліуса на викривлених межах (параметр променя p = const уздовж променя):
  p = sin(i)/v(z) = r sin(i)/v(r)    (r = радіальна відстань)
  Промені викривляються назад вгору, коли v(z) зростає з глибиною (звичайний випадок Землі)
  Глибина розвороту променя: p = r_turn/v(r_turn) → розв'язується для r_turn при заданому p

Тіньові зони:
  Тіньова зона P-хвиль: кутова відстань 103°–142° від епіцентру
    Зовнішнє ядро сильно заломлює P-хвилі всередину → розрив у прямих приходах P
    Фаза PKP: P-хвилі, що проходять крізь ядро (приходять >142°)
    PKIKP: проникає у внутрішнє ядро; приходить 0°–180° із затримкою часу пробігу
  Тіньова зона S-хвиль: 103°–180°
    S-хвилі не можуть увійти в рідке зовнішнє ядро → відсутні за межами 103°
    Відкриття рідкого зовнішнього ядра: Олдхем (1906), Гутенберг (1914)

Поверхневі хвилі (дисперсійні, залежна від частоти глибина зондування):
  Релея: ретроградний еліптичний рух частинок у сагітальній площині (зв'язок SV + P)
    Групова швидкість ~0,92 v_S(глибина ≈ λ/3)
  Лява: горизонтальний поперечний рух; потребує шаруватості швидкості для напрямку
    Фазова швидкість c_L: 2kH tan(k_L H) = (ρ₁v_L1)/(ρ₀v_S0)...  (дисперсійне співвідношення шару)
  Томографія поверхневих хвиль: період 20–300 с → чутливість до глибини 30–600 км
          

3. Сейсмологія: магнітуда, фокальні механізми та статистика

Кількісна оцінка розміру землетрусу вимагає відокремлення амплітуди хвилі, зареєстрованої приймачем (функції відстані, діаграми випромінювання, згасання та відгуку приладу), від власної сили джерела. Шкала моментної магнітуди досягає цього, прив'язуючи магнітуду до фізичного сейсмічного моменту M₀, який інтегрує падіння напруження по площі розриву.

Локальна магнітуда за Ріхтером:
  M_L = log₁₀(A/A₀)    [A = максимальна амплітуда зміщення в мкм на 100 км, Вуд-Андерсон]
  Насичується при M_L ≈ 7 (сейсмічний момент росте швидше за амплітуду)

Сейсмічний момент:
  M₀ = μ · Ā · D̄    [μ = модуль зсуву ~30 ГПа, Ā = площа розриву, D̄ = середнє зміщення]
  Виводиться з довгоперіодних об'ємних хвиль або статичної деформації GPS

Моментна магнітуда (Хенкс і Канаморі, 1979):
  M_w = (2/3)(log₁₀ M₀ − 9,1)    [M₀ у Н·м]
  Не насичується; прив'язана до вивільнення енергії
  Приклади:
    Суматра 2004: M_w = 9,1, M₀ ≈ 4×10²² Н·м, площа розлому ≈ 1300×200 км²
    Вальдівія 1960: M_w = 9,5 (найбільший інструментально зафіксований; Чилі)
    Удар Чиксулуб: M_w ≈ 10,5–11 (оцінка)

Фокальний механізм (розв'язок площини розлому):
  Тензор моменту M_ij: симетричний тензор другого порядку, що описує пари сил
  Власні значення: осі T (розтяг), N (нуль), P (стиснення)
  «Пляжний м'яч»: стереографічна проєкція нижньої півкулі полярностей першого руху
    Квадрантний шаблон стиснення (вниз-перше)/розрідження (вгору-перше)
  Типи розломів:
    Нормальний:      висяча стінка вниз (дивергентна межа; розтягуюче нормальне напруження)
    Насув/зворотний: висяча стінка вгору (конвергентна; стискаюче)
    Зсувний:         горизонтальний рух (трансформний розлом)

Співвідношення Гутенберга-Ріхтера:
  log₁₀ N(M) = a − bM
  b ≈ 1,0 (глобально; локально в діапазоні 0,6–1,3)
  Означає: на кожен M7,0 припадає ~10 M6,0 та ~100 M5,0 землетрусів
  Фізична інтерпретація: масштабно-незалежна фрактальна геометрія розломів

Спад афтершоків Омори-Уцу:
  n(t) = K/(t + c)^p    [n = частота афтершоків, K,c,p = константи; p ≈ 1,0]
  Продуктивність афтершоків: N_total ∝ 10^{α(M_main − M_min)}, α ≈ 1
          

4. Тектоніка плит і мантійна динаміка

Теорія тектоніки плит, консолідована в 1960-х роках (Гесс, Вайн, Метьюз, Вілсон, Мак-Кензі, Паркер), об'єднала раніше окремі галузі дрейфу континентів, спредингу морського дна, сейсмічності та вулканізму в єдиний кінематичний фреймворк. Двигун, що керує цим процесом — мантійна конвекція — залишається активною сферою досліджень, оскільки залежить від реології силікатних порід за екстремального тиску та температури.

Реологія мантії:
  Ефективна в'язкість: η_eff ≈ 10^{21} Па·с (верхня мантія, післяльодовиковий підйом)
                        10^{22}–10^{23} Па·с (нижня мантія)
  Дислокаційна + дифузійна повзучість: ε̇ = A σⁿ exp(−E_a/RT) (степенева повзучість, n≈3,5)
  Перехід від пружного до в'язкого: t_Maxwell = η/G ≈ 10^{21}/70×10^9 ≈ 450 років

Безрозмірні числа мантійної конвекції:
  Релея: Ra = ρ_0 g α ΔT d³/(κη)
    ρ₀ ≈ 4000 кг/м³; g = 10 м/с²; α ≈ 2×10⁻⁵ К⁻¹; ΔT ≈ 2500 К
    d ≈ 2900 км; κ ≈ 10⁻⁶ м²/с; η ≈ 3×10²¹ Па·с
    Ra ≈ 10⁷ → активна конвекція (початок ~Ra_c = 1100 для плоского шару)
  Застигла кришка проти рухомої кришки: залежить від границі текучості холодної літосфери

Спрединг морського дна та магнітні аномалії:
  Швидкість: середньоатлантичний хребет ~2,5 см/рік (повільний); Східнотихоокеанське підняття ~15 см/рік (швидкий)
  Вайн-Метьюз-Морлі (1963): нова океанічна кора набуває напрямку геомагнітного поля
    → чергування смуг нормальної/зворотної полярності, паралельних осі хребта
  Вік океанічної кори: макс. ~200 млн років (юра); старша субдукована
  Вік морського дна = (відстань від хребта) / (швидкість спредингу)

Швидкості плит GPS (ITRF2020 / NNR-MORVEL56):
  Пара плит                Швидкість (мм/рік)  Азимут
  Тихоокеанська–Північноамериканська  ~50           ПнПнЗ
  Індо-Австралійська–Євразійська      ~45           ПнПнС
  Африканська–Євразійська             ~7            ПнПнЗ
  Північноамериканська–Євразійська    ~24           СхПнС (розкриття Атлантики)
  Нубійська–Сомалійська               ~6            СхПнС (Східноафриканський рифт)

Цикл Вілсона (повне розкриття та закриття океану):
  1 – Континентальний рифтинг (Східноафриканський рифт сьогодні: рання стадія)
  2 – Прото-океан (Червоне море: вузький океан)
  3 – Відкритий океан (Атлантика: зрілі пасивні окраїни)
  4 – Ініціація субдукції → острівна дуга (Наска–Південна Америка)
  5 – Закриття океану + континентальне зіткнення (Альпи, Гімалаї, Аппалачі)
  Тривалість: ~500 млн років повний цикл
          

5. Гравітаційне поле Землі та геоїд

Геоїд — це еквіпотенціальна поверхня гравітаційного поля Землі, що (в середньому) збігається із середнім рівнем моря. Відхилення від еліпсоїда відображають аномалії густини в корі та мантії. Супутникова місія GRACE (2002–2017) та її наступниця GRACE-FO (з 2018) вимірюють гравітаційне поле до сферичного гармонічного ступеня 120 (розрізнення ±330 км), відстежуючи зміну відстані між двома апаратами.

Розклад гравітаційного поля:
  V(r,θ,λ) = (GM/r)∑_l ∑_m (R/r)^l P_l^m(cosθ) [C_lm cos(mλ) + S_lm sin(mλ)]
  EGM2008 (NGA): повний до l_max = 2159 (розрізнення ~10 км)
  C₂₀ = J₂ (сплюснутість Землі): −1,0826 × 10⁻³  (сплюснутий сфероїд)

Гравітаційні поправки (зведення до стандартної поверхні):
  Поправка вільного повітря: +0,3086 мкГал/м вгору (−3,086 мкГал/м висота спостереження)
  Поправка Буге: −2πGρh ≈ −0,1119ρh мГал/м  (маса шару між спост. та геоїдом)
    Повна поправка Буге: поправка на рельєф Tc усуває місцеву нерівність топографії

  Аномалія вільного повітря: Δg_FA = g_obs − g_normal + FAC         [усуває ефект висоти]
  Аномалія Буге:  Δg_B  = Δg_FA − 2πGρh + Tc             [усуває шар кори]

Ізостазія (моделі Ейрі проти Пратта):
  Ейрі (Джордж Біддел Ейрі, 1855): стала густина ρ_c; корінь гори плаває
    H_root = h · ρ_c/(ρ_m − ρ_c) ≈ 7h   (h = висота гори, надлишкова густина занурена)
    Передбачення: від'ємна аномалія Буге під горами (спостерігається глобально)
  Пратт (Джон Генрі Пратт, 1855): змінна густина σ(x), стала базова глибина T_0
    σ(x) · (T_0 + h) = σ₀ · T_0   (компенсація латеральною зміною густини)

Застосування GRACE:
  Маса льоду Гренландії: −270 Гт/рік (тренд 2003–2022); Антарктида: −150 Гт/рік
  Виснаження підземних вод: басейн Інду −24 мм/рік еквівалентної висоти води
  Післяльодовиковий підйом: Фенноскандія +10 мм/рік зростання гравітації
  Гідрологічний цикл: сезонні коливання висоти геоїда 10–20 мм у басейнах Амазонки/Конго
          

6. Геомагнітне поле та теорія динамо

Магнітне поле Землі створюється турбулентною конвекцією в рідкому залізному зовнішньому ядрі. Поле переважно дипольне (осьовий диполь ≈ 80% поверхневого поля), але змінюється в масштабах часу від мілісекунд (вікова варіація від електромагнітно пов'язаної динаміки ядра) до мільйонів років (інверсії полярності). Палеомагнітні записи в базальтах морського дна зберігають історію цих інверсій аж до юрського періоду.

Джерела геомагнітного поля:
  Поле ядра: 95% загального; створюється МГД-динамо; змінюється в масштабі роки–млн років
  Кірково-літосферне поле: залишковий магнетизм; аномалії ±2000 нТ; квазістатичне
  Зовнішнє поле: струми магнітосфери + іоносфери; варіація годинники-дні; спокійне Sq ≈ 30 нТ

Рівняння МГД-індукції:
  ∂B/∂t = ∇×(u×B) − ∇×(η ∇×B)
  η = 1/(μ₀σ) = магнітна дифузійність [м²/с]
  Баланс: магнітне число Рейнольдса Rm = u L/η
  Зовнішнє ядро Землі: u ≈ 0,3 мм/с, L ≈ 2000 км, η ≈ 1 м²/с → Rm ≈ 600 >> 1
  Високе Rm: дифузія незначна → поле «заморожене» в рідині (теорема Альфвена)

Модель стовпчастої конвекції Буссе (1975):
  Обертання планети (Ω = 7,3×10⁻⁵ рад/с) → Коріоліс організує конвекцію в
  циліндричні валки, вирівняні з віссю обертання
  Тейлор-Праудмен: у швидко обертових рідинах потік 2D, перпендикулярний до Ω
  Конвективні валки створюють спіральний потік → α-ефект (динамо середнього поля):
    ⟨u×B⟩ = αB + ... (Штеенбек, Краузе, Редлер, 1966)
  Ω-ефект: диференціальне обертання розтягує полоїдальне B у тороїдальне
  Разом: αΩ-динамо → підтримує обидва компоненти проти дифузії

IGRF (Міжнародне геомагнітне еталонне поле):
  Сферичне гармонічне представлення внутрішнього поля до порядку/ступеня 13
  Поточний дипольний момент: 7,94×10²² А·м² (зменшився на 6% з 1840 року)
  Південноатлантична аномалія (SAA): регіон слабкого поля; потік частинок авроральний на низькій висоті
  Геомагнітний Північний полюс (2024): ~86°Пн, 133°Зх (дрейфує до Сибіру ~40 км/рік)

Палеомагнітні інверсії:
  Шкала часу геомагнітної полярності (GPTS) з океанічних магнітних аномалій + біостратиграфії
  Поточний: нормальний хрон Брюн (з 780 тис. років тому)
  Остання інверсія: Матуяма-Брюн, 781 тис. років тому (тривалість переходу 10–20 тис. років)
  Середня частота інверсій: 4–5 за млн років (кайнозой); Крейдяний нормальний суперхрон: 40 млн років без інверсій
  Провісники: екскурси, мінімуми інтенсивності, мультиполярна структура поля під час переходу
  Потенційні наслідки інверсії: 10-кратне розширення SAA, полярне сяйво на середніх широтах
          

Спробуйте ці симуляції