Як працюють вулкани — магма, тиск і фізика виверження

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

1. Утворення магми в мантії

Мантія — це переважно тверда порода — перидотит — що утримується твердою завдяки тиску попри температури 1 300–1 500 °C. Плавлення відбувається, коли цю стабільність порушує один із трьох механізмів:

Декомпресійне плавлення: на серединно-океанічних хребтах висхідна мантійна порода декомпресується, піднімаючись. Температура солідуса падає швидше, ніж порода охолоджується, що породжує розплав.
Флюїдне плавлення: у зонах субдукції вода, що виділяється з плити, яка занурюється, знижує солідус матеріалу мантійного клина зверху — вологий перидотит плавиться за значно нижчих температур, ніж сухий.
Теплоперенесення: мантійні плюми («гарячі точки») доставляють аномально гарячий матеріал із глибокої мантії, породжуючи розплав під внутрішньоплитовими гарячими точками, як-от Гаваї та Ісландія.

Частки часткового розплаву зазвичай становлять 1–30%. Вміст кремнезему в отриманій магмі визначає її в’язкість і тип виверження. Базальтові магми (низький SiO₂, ~50%) гарячі й рідкі; ріолітові магми (високий SiO₂, ~75%) в’язкі та насичені газом, схильні до вибухових вивержень.

η_basalt \approx 10²-10³ Pa\cdots (як мед за температури виверження ~1200 °C) η_rhyolite \approx 10⁶-10¹² Pa\cdots (майже тверда за поверхневих умов) В’язкість подвоюється приблизно з кожним зростанням SiO₂ на 10% і зростає експоненційно зі зниженням температури.

2. Підйом за рахунок плавучості

Магма (густина ≈ 2 500–2 700 кг/м³) менш густа за навколишній перидотит (≈ 3 300 кг/м³) і піднімається завдяки плавучості. Рушійний градієнт тиску такий:

ΔP/Δz = (ρ_rock - ρ_melt) \cdot g \approx (3300 - 2600) \cdot 9.8 \approx 6 860 Pa/m Для стовпа 10 км: ΔP \approx 68 MPa — достатньо, щоб розколоти кору.

Розплав збирається в кишені, потім мігрує крізь проникну породу за низької частки розплаву й зрештою — у більші канали чи дайки — вертикальні листоподібні тріщини, крізь які магма швидко піднімається. Поширення дайки рухається внутрішнім тиском магми, що перевищує в’язкість руйнування вмісної породи.

Піднімаючись, магма декомпресується, дозволяючи розчиненим газам (H₂O, CO₂, SO₂) виділятися й утворювати бульбашки — процес, званий везикуляцією. Це критичний крок, що відрізняє вибухові виверження від ефузивних.

Розчинність летких речовин: за 200 MPa (≈8 км глибини) базальтовий розплав розчиняє ~4 мас.% H₂O. Коли тиск падає під час підйому, розчинність води різко зменшується: за 10 MPa вона утримує лише ~0.5 мас.%. Виділення прискорюється експоненційно, коли магма наближається до поверхні.

3. Магматичні камери та накопичення тиску

У багатьох вулканічних системах висхідна магма накопичується в магматичній камері на середньокорових глибинах (5–20 км). Камера — це частково розплавлена зона, не просте рідке «озеро», а кристалічна каша з 40–70% твердих кристалів та силікатною розплавною міжкристалічною фазою.

Тиск у камері зростає, коли:

Нова магма проникає з глибини (поповнення)
Виділення летких речовин збільшує ефективний об’єм
Теплове розширення, коли камера нагрівається

Виверження запускається, коли надлишковий тиск ΔP перевищує межу міцності на розтяг породи, що лежить зверху (~10–40 MPa для типової кори). Сейсмічна томографія та вимірювання деформації ґрунту (InSAR) дозволяють вулканологам відстежувати роздування камери в реальному часі.

ΔP_critical \approx T_f + ρ_melt \cdot g \cdot h_chamber де T_f = в’язкість руйнування породи покрівлі (~10-30 MPa) Підняття ґрунту δ \approx (ΔP \cdot V_chamber) / (4π μ_crust d²) (модель точкового джерела Моги для поверхневої деформації)

4. Типи виверження

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

Гавайські виверження

Низьков’язка базальтова лава спокійно витікає у вигляді лавових фонтанів і потоків. Газ виходить поступово. Швидкості потоку можуть сягати 10 км/год на крутих схилах. Класичні приклади: Кілауеа (Гаваї), Стромболі (Італія).

Стромболіанські виверження

Періодичні помірні вибухи кожні кілька хвилин або годин. Великі газові кишені піднімаються каналом і лопаються на поверхні, викидаючи розжарені лавові бомби до 300 м. Характерні для базальтових та базальт-андезитових магм.

Вулканські виверження

Короткі потужні вибухи, коли затверділу пробку, що ущільнює жерло, вибиває газ під надлишковим тиском. Утворює густі попільні хмари та балістичні брили на дальність до кількох км.

Плінійські виверження

Найпотужніший тип — названий за описом виверження Везувію 79 року н. е. Плінієм Молодшим. Стійкий струмінь газу, попелу й пемзи утворює еруптивну колону, що може сягати 40–50 км у стратосферу. Швидкість виверженої маси може перевищувати 10⁸ кг/с.

Висота колони H \approx k \cdot Q^(1/4) (рівняння плюму Мортона-Тейлора-Тернера) де Q = масовий потік (кг/с), k \approx 0.24 км/(кг/с)^(1/4) Q = 10⁸ кг/с \to H \approx 0.24 \times (10⁸)^(0.25) \approx 24 км

Кальдероутворювальні «супервиверження»

Коли дуже велика магматична камера частково спорожнюється, покрівля може обвалитися всередину, утворюючи кальдерну западину (Єллоустон, Тоба, Кампі-Флегреї). Ці події вивільняють >1 000 км³ матеріалу (VEI ≥ 8) і можуть спричинити вулканічні зими.

5. Реологія потоку лави

Лава поводиться як бінгамівська пластична речовина з границею плинності τ_y — вона тече, лише коли напруга зсуву перевищує поріг плинності (через те, що кристали й везикули утворюють мережу). Понад τ_y вона дотримується ньютоноподібного співвідношення:

τ = τ_y + η_plastic \cdot (du/dy) (бінгамівське визначальне рівняння) Для канального потоку ширини W і глибини h на схилі θ: U_max = (ρg sin θ / 2η) \cdot (h - h_y)² де h_y = τ_y / (ρg sin θ) — товщина жорсткої пробки на поверхні.

Охолоджуючись, лава утворює тверду кірку навіть тоді, коли всередині вона лишається рухомою. Це створює:

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

Примітка щодо симуляції: SPH (згладжена гідродинаміка частинок) та LBM (метод ґраткового Больцмана) — основні чисельні методи для симуляції потоків лави, включно з неньютонівською реологією та тепловим зв’язком. Симуляція тектонічних плит використовує спрощену в’язку модель.

6. Пірокластичні потоки та тефра

Коли плінійська еруптивна колона обвалюється, або коли лавові куполи розколюються, суміш гарячого газу та фрагментованої породи (пірокласти) стає густішою за повітря й тече як гравітаційна течія — пірокластична густинна течія (PDC).

Швидкість: до 700 км/год на крутих схилах
Температура: 300–800 °C усередині
Дальність: 10–100 км від джерела
Смертоносність: відповідальна за більшість вулканічних смертей — викид у Помпеях 79 року н. е. досяг 500 °C за секунди

Тефра (вулканічні фрагменти в повітрі) класифікується за розміром зерна: вулканічні брили (>64 мм), лапілі (2–64 мм), вулканічний попіл (<2 мм). Дрібний попіл може лишатися завислим у стратосфері 1–3 роки, розсіюючи сонячне світло та знижуючи глобальні температури до 0.5°C (Пінатубо 1991 року охолодив Землю на ~0.4°C протягом 2 років).

Швидкість осідання за Стоксом для частинки тефри: v_s = (2r² \cdot (ρ_particle - ρ_air) \cdot g) / (9 η_air) Для попелу r = 25 μm: v_s \approx 0.03 м/с \to лишається в повітрі днями Для лапілі r = 1 мм: v_s \approx 50 м/с \to осідає за хвилини

7. Зв’язок із симуляцією

Симуляція тектонічних плит моделює великомасштабну мантійну конвекцію, що рухає плити й створює умови для вулканізму. Зони субдукції — де одна плита занурюється під іншу — є джерелом 75% вулканічних вивержень Землі.

Симуляція використовує спрощену модель в’язкої рідини на 2D-сферичній поверхні з залежною від температури в’язкістю. Коли плити сходяться, холодніший густіший матеріал занурюється, спричиняючи зворотний потік, що може виносити гарячі плюми до поверхні — візуальний аналог механізму флюїдного плавлення, описаного вище.

🌋 Відкрити симуляцію тектонічних плит →