Тектонічні плити та конвекція мантії

Будова Землі

Земля — не однорідне тверде тіло. Її надра шаруваті, мов цибулина, і кожен шар має власний склад, температуру й механічну поведінку:

Літосфера — кора плюс жорстка найверхніша частина мантії — лежить поверх астеносфери, частково розплавленої, механічно слабкої зони у верхній мантії на глибині приблизно 80–200 км. Саме на цьому м’якому шарі тектонічні плити «плавають» і ковзають.

Конвекція мантії

Хоч мантія тверда, у геологічних масштабах часу (мільйони років) вона поводиться як надзвичайно в’язка рідина. Її в’язкість становить приблизно 10²¹ Па·с — приблизно у 10²⁴ разів в’язкіша за воду, та повільна теплова конвекція все одно можлива.

Механізм — проста термодинаміка: гаряча речовина глибоко поблизу ядра менш щільна й піднімається; холодніша, щільніша речовина поблизу основи літосфери опускається. Так виникають конвективні комірки, що тягнуть за собою плити над ними.

Основні рівняння — це рівняння Стокса у поєднанні з рівнянням теплоперенесення. Для нестисливої течії в мантії:

Де u — швидкість, P — тиск, η — в’язкість, ρ — щільність, T — температура, а κ — коефіцієнт температуропровідності. Перші два рівняння означають, що мантія тече, щоб урівноважити тиск і в’язкі сили проти гравітації — інерцією на цих масштабах можна знехтувати. Третє рівняння означає, що температура переноситься течією й розсіюється теплопровідністю.

Безрозмірне число Релея визначає, чи відбувається конвекція:

Де α — коефіцієнт теплового розширення, ΔT — різниця температур поперек шару, а d — товщина шару. Конвекція починається, коли Ra > ~1000. Для мантії Ra ≈ 10⁷, тож за цим критерієм конвекція дуже інтенсивна.

Типи меж плит

17 основних тектонічних плит взаємодіють на трьох видах меж, кожен з яких породжує характерні геологічні форми:

Дивергентні межі

Плити розходяться. Розрив заповнюється висхідною магмою, яка застигає, утворюючи нову океанічну кору — процес, що зветься спредингом морського дна. Серединно-Атлантичний хребет — найдовше гірське пасмо у світі: 16 000 км вулканічного хребта, що тягнеться центром Атлантики. Характерні неглибокі землетруси й базальтовий вулканізм.

Конвергентні межі

Плити стикаються. Щільніша океанічна плита заштовхується під легшу континентальну плиту в процесі, що зветься субдукцією. Занурювана плита опускається в мантію, тягнучи за собою поверхневу плиту (тяга плити, див. нижче). Зони субдукції породжують глибокофокусні землетруси (до 700 км), вибухові острівно-дугові вулкани й акреційні призми зіскребленого осаду. Тихоокеанське «Вогняне кільце» позначає ланцюг зон субдукції.

Там, де стикаються дві континентальні плити, жодна не субдукується легко (обидві плавучі). Натомість кора зминається вгору у гірські пасма. Так утворилися Гімалаї, коли Індія зіткнулася з Азією ~40–50 мільйонів років тому.

Трансформні межі

Плити горизонтально ковзають одна повз одну. Напруження наростає, коли нерівності блокують розлом, і зрештою вивільняється у вигляді великих землетрусів, коли блокування зривається. Розлом Сан-Андреас у Каліфорнії — класичний приклад: Тихоокеанська плита рухається на північний захід відносно Північноамериканської плити зі швидкістю ~5 см/рік.

Рушійні сили

Що насправді рухає плити? Відповідь охоплює кілька конкурентних сил:

• Тяга плити (домінантна, ~90%): стара, холодна океанічна літосфера щільніша за мантію під нею. У зонах субдукції вона опускається під власною вагою, тягнучи решту плити за собою, мов скатертину. Це найсильніший окремий рушій.
• Поштовх від хребта (~10%): на серединно-океанічних хребтах гаряча нова кора піднята й гравітаційно сповзає вниз геть від хребта. Менший, але реальний внесок.
• Мантійне тертя: конвекція в астеносфері чинить в’язке тертя на основу плит. Сучасні дослідження вказують, що воно може як рухати, так і гальмувати рух плит залежно від геометрії конвективної комірки.
• Базальне всмоктування: висхідні мантійні плюми (гарячі точки) під континентами можуть створювати ефект усмоктування, що розриває плити, започатковуючи рифтогенез.

Як працює симуляція

Браузерна симуляція тектонічних плит використовує спрощений, але фізично обґрунтований підхід. Замість розв’язання повних рівнянь Нав’є–Стокса (обчислювально неприйнятних на реальних масштабах) вона використовує три пов’язані моделі:

1. Сегментація плит за Вороним

Поверхня кулі поділена на діаграму Вороного з N плит. Кожна плита — це опуклий багатокутник із трикутників сітки. Зерна Вороного рухаються по сфері відповідно до векторів швидкості, отриманих із поля мантійної конвекції знизу.

2. Конвекція мантії через 2-вимірну сітку

Під корою сітка 256 × 256 зберігає температуру й швидкість. На кожному кадрі:

1. Адвекція температури: T_new[i] = T_old[i − v*dt] (білінійна інтерполяція)
2. Дифузія температури: 5-точковий лапласіанів шаблон з κ = 0.01
3. Обчислення швидкості плавучості з градієнта температури за допомогою спрощеного розв’язувача Стокса

Результат — поле швидкостей, що рухає зерна плит.

3. Деформація меж

На кожному граничному ребрі між плитами вектор відносної швидкості визначає тип межі:

• Дивергентна: розрив відкривається → вставляються нові трикутники, забарвлені в червоний (гаряча кора)
• Конвергентна: плити перекриваються → трикутники вилучаються, поверхня вигинається вгору шляхом зсуву Y-координат вершин уздовж нормалі до межі
• Трансформна: вершини ковзають по дотичній, нерівність межі накопичується як «енергія землетрусу», що вивільняється сплеском частинок, коли поріг перевищено

Горотворення наближено передається підняттям висот вершин пропорційно швидкості конвергенції з часом із подальшим застосуванням ґаусового згладжування для злиття рельєфу. Результат візуально переконливий, хоч і не геологічно точний.

Алгоритм оновлення сітки

Кожні 60 кадрів симуляції сітка перетріангульовується за допомогою обмеженої тріангуляції Делоне, що запобігає накопиченню вироджених трикутників на межах плит. Ключова ідея в тому, що більшість вершин лишаються нерухомими — лише граничні вершини в межах однієї ширини плити від активного ребра потребують перетріангуляції:

function updateBoundaryMesh(plateA, plateB, velocity) {
  const normal   = boundary.normal;
  const relSpeed = velocity.dot(normal);          // + = конвергентна, - = дивергентна
  const dt       = FRAME_TIME * SIM_SPEED;

  if (relSpeed > 0) {
    // Конвергентна: піднімаємо рельєф уздовж межі
    boundary.vertices.forEach(v => {
      v.y += relSpeed * dt * MOUNTAIN_SCALE;
      v.y  = Math.min(v.y, MAX_HEIGHT);
    });
    smoothTerrain(plateA, 3);
    smoothTerrain(plateB, 3);
  } else {
    // Дивергентна: вставляємо нові вершини, призначаємо колір молодої кори
    const newVerts = interpolateBoundary(boundary, -relSpeed * dt);
    insertVertices(plateA, newVerts, HOT_CRUST_COLOR);
  }
}

«Мантію можна уявити як гігантську лава-лампу — нагріту знизу ядром і охолоджувану зверху поверхнею. Воскові краплі, що пливуть угору й униз, — це мантійні плюми та низхідні потоки, а зовнішня оболонка лампи — наші континенти, що повільно дрейфують.»

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

Інтерактивна симуляція дає змогу спостерігати, як формуються конвективні комірки, як плити дрейфують і стикаються, здіймаються гори й відкриваються рифти — усе в реальному часі. Ви можете змінювати в’язкість мантії, потік тепла з ядра й кількість плит та спостерігати цикл суперконтиненту в пришвидшеному відтворенні.