Лето на 69-й параллели. Научная фантастика

1.1. Пространственно-временная структура вертикальных движений

Неотектонические движения в Антарктиде в подавляющем большинстве обусловлены гляциоизостатической корректировкой (GIA) — реакцией литосферы и мантии на изменение ледовой нагрузки после последнего ледникового максимума (21 тыс. лет назад). Сеть непрерывных GNSS-станций POLENET [1] (127 пунктов, 1995—2025) позволяет измерять вертикальные скорости с точностью до 0,5 мм/год. Результаты показывают ярко выраженную асимметрию: Восточно-Антарктический кратон (мощность литосферы до 250 км) поднимается со скоростями 0—3 мм/год, в то время как Западная Антарктида (литосфера 60—80 км, тепловой поток 85—100 мВт/м²) демонстрирует подъём до 15 мм/год в районе Земли Мэри Бэрд и Антарктического полуострова. Исключение — бассейн Уилкса (Восточная Антарктида), где наблюдается слабый прогиб (—6 мм/год), интерпретируемый как остаточная депрессия от более мощного льда в прошлом.

Для разделения GIA и тектонических вертикальных движений использованы три независимые модели: ICE-6G_C [2], IJ05_R2 [3] и W12a [4]. Расхождения моделей в глубинных районах Восточной Антарктиды (Купол C, Купол A) достигают 3—4 мм/год, что вносит систематическую ошибку в оценки горизонтального дрейфа (см. главу 2). Причина расхождений — неопределённость в истории дегляциации внутренних областей Восточной Антарктиды, где лед оставался стабильным последние 10 млн лет.

1.2. Горизонтальные деформации как индикатор внутриплитной активности

Долгое время Антарктиду считали асейсмичной и тектонически пассивной, однако современные сейсмические сети (GSN, AGAP, TAMNNET) зафиксировали сотни событий с магнитудами Mw 2—4,5, а также несколько событий Mw 5—6 в море Скотия. Пространственное распределение эпицентров приурочено к трём структурам:

Трансантарктические горы (ТАМ) — разломы Такер, Лип-Йир, Лантерман, Авиатор с правосдвиговой кинематикой. Скорости скольжения, оцененные по повторным GPS-измерениям, составляют 0,4—0,9 мм/год [41].

Западно-Антарктическая рифтовая система (WARS) — нормальные сбросы и бассейны растяжения (море Росса, бассейн Бирдмора). Здесь деформации сдвига достигают 10⁻⁸ год⁻¹ [42].

Антарктический полуостров и остров Десепшн — активный вулканизм и вулкано-тектоническая сейсмичность (последнее извержение 1970 г.) [27].

Наличие измеримых внутриплитных деформаций означает, что приближение «абсолютно жёсткой плиты» является только нулевым приближением. В интегральной модели (глава 3) это учтено через эффективную вязкость под кратоном и допущение о том, что на временах ≥50 млн лет деформации отдельных блоков могут быть усреднены.

1.3. Влияние ледовой нагрузки на момент инерции плиты

По данным Bedmap2 [15], современная ледовая нагрузка распределена крайне неравномерно: центр масс льда смещён относительно центра масс литосферной плиты на 320 км в направлении Трансантарктических гор. Это смещение создаёт асимметрию в тензоре инерции плиты, что приводит к появлению дополнительного гравитационного момента, когда плита вращается. Величина этого момента мала по сравнению с моментом от спрединговых сил (отношение ∼0,02), но, как показано в главе 3, она не равна нулю и должна учитываться в точных прогнозах.

2.1. Палеомагнитные реконструкции на интервале 50—0 млн лет

Для углублённого анализа на интервале ±50 млн лет были собраны и обработаны палеомагнитные данные из базы PALEOMAGIA [40] по Антарктиде за последние 50 млн лет. Всего использовано 28 определений виртуальных геомагнитных полюсов (VGP) для интервалов: 50—40, 40—30, 30—20, 20—10, 10—0 млн лет. Типичная погрешность определений составляет ±3—5° по широте. Получена кривая широтного дрейфа центра масс плиты: 50 млн лет назад — около 63° ю.ш., 30 млн лет назад — 68° ю.ш., 10 млн лет назад — 73° ю.ш., современное положение — 75,5° ю. ш. Эти данные служат основой для ретроспективной верификации модели в главе 4.

2.2. Современная кинематика по GNSS и глобальным моделям

Обработка данных 68 непрерывных GNSS-станций (1995—2025) выполнена в системе ITRF2020 с использованием программ GAMIT/GLOBK и GIPSY-OASIS. Горизонтальные скорости станций подчиняются вращательному полю с высокой когерентностью (среднее отклонение от модели жёсткого вращения составляет 0,8 мм/год, что близко к уровню шума). Получены следующие параметры в различных системах отсчёта:

NNR-MORVEL56 [6] (система без суммарного вращения, близкая к мантийной): полюс вращения (58,2° с.ш., 127,5° з.д.), угловая скорость ω = 0,219 ± 0,007°/млн лет.

ITRF2020 [19] (центр масс Земли, включает вращение относительно мантии): (59,1° с.ш., 126,8° з.д.), ω = 0,225 ± 0,006°/млн лет.

Горячие точки (модель HS3-NUVEL-1A): (57,0° с.ш., 129,2° з.д.), ω = 0,210 ± 0,010°/млн лет.

Различия между системами достигают 0,015°/млн лет (около 7% от ω), что важно для долгосрочных прогнозов. В данной монографии в качестве базовой используется система NNR-MORVEL56 как наиболее приближённая к мантийной системе отсчёта. Линейная скорость центра масс плиты (приблизительно 75,5° ю.ш., 120° в.д.) составляет 11,2 ± 0,4 мм/год, направление — азимут 343°(северо-северо-запад).

2.3. Дрейф полюса вращения: терминология и оценки на интервале ±50 млн лет

В монографии принято следующее разграничение: истинное полярное блуждание (TPW) — движение всей Земли относительно оси вращения; относительный дрейф полюса плиты — изменение ориентации вектора Ω в системе горячих точек. Для интервала ±50 млн лет скорость относительного дрейфа полюса Антарктической плиты составляет 0,09 ± 0,03°/млн лет (по [7]) и считается постоянной в первом приближении. За 50 млн лет полюс смещается на 4,5 ± 1,5°. Это значение используется как в ретроспективном, так и в прогнозном расчётах.

2.4. Вековые изменения угловой скорости: экстраполяция на ±50 млн лет

Аппроксимация палеомагнитных данных и реконструкций океанических изохрон [10] даёт угловую скорость как функцию времени: ω (t) = ω₀ * exp (-t/τ) с ω₀=0,219°/млн лет, τ=78±12 млн лет. Здесь t отсчитывается от современности в будущее, а в прошлое — t отрицательное, т.е. ω (t) = ω₀ * exp (-|t|/τ) для t <0 (симметричная экспонента). Однако данные показывают, что 50 млн лет назад ω была выше — около 0,34°/млн лет, что соответствует τ=78 млн лет. Таким образом, экспоненциальная модель с постоянным τ справедлива как в прошлое, так и в будущее, но с разной начальной амплитудой (в будущем — от современного значения, в прошлом — от более высокого). В расчётах используется единая формула: ω (t) = ω₀ * exp (-t/τ) для t≥0 (будущее) и ω (t) = ω₀ * exp (+t/τ) для t <0 (прошлое), где t отрицательно. Это обеспечивает непрерывность в t=0.

3.1. Полный вывод уравнений движения плиты с учётом неотектонического момента

Рассмотрим литосферную плиту как твёрдое тело с переменным тензором инерции I (t), которое движется в вязкой мантии. Уравнение вращательного движения в системе отсчёта, связанной с центром масс Земли (невращающейся, но с ускоренным движением центра масс — пренебрегаем), имеет вид:

d/dt (I·Ω) = M_ridge (t) + M_drag (t) + M_load (t) + M_TPW (t), (1)

где Ω — вектор мгновенной угловой скорости плиты; M_ridge — момент от спрединговых толкающих сил (ridge push), распределённых по границам плиты; M_drag — вязкий момент на подошве литосферы от мантийного течения; M_load — момент, создаваемый несимметричной ледовой нагрузкой (неотектоническая компонента); M_TPW — момент от истинного полярного блуждания (влияет на переориентацию плиты относительно оси вращения Земли, но не на относительное движение в мантии — в рамках модели NNR этим членом можно пренебречь, считая TPW уже включённым в M_drag через градиент мантийного потока).

Для квазистационарного движения (dΩ/dt мало, характерное время изменения Ω много больше периода релаксации астеносферы, ~10⁶ лет) левая часть уравнения (1) обращается в ноль, и мы получаем баланс моментов:

M_ridge + M_drag + M_load = 0. (2)

В приближении, что M_drag пропорционален Ω (линейное вязкое сопротивление), можно записать:

M_drag = -γ · Ω, (3)

где γ — тензор вязкого сопротивления (в изотропном приближении — скаляр). Тогда из (2) и (3):

Ω = (1/γ) (M_ridge + M_load). (4)

Теперь необходимо выразить M_load через параметры ледовой нагрузки. Ледовая нагрузка создаёт на поверхности литосферы вертикальное напряжение σ_zz = -ρ_ice g H (x,y) (H — мощность льда). Изостатический прогиб литосферы h_iso (x,y) определяется уравнением упругого изгиба пластины на вязком основании. Для длинноволновых компонент (λ> 100 км) прогиб приближённо равен: h_iso = — (ρ_ice/ρ_mantle) H, где ρ_mantle ≈ 3300 кг/м³. Тогда вертикальное смещение поверхности литосферы создаёт наклон, который приводит к появлению касательных напряжений на подошве литосферы. Через эти касательные напряжения и интеграл по площади получается момент сил:

M_load = ∫∫ (r — r_CM) × [-ν ρ_ice g ∇ h_iso] dS, (5)

где ν — эффективная вязкость астеносферы (в м²/с, после преобразования размерностей — коэффициент пропорциональности между градиентом давления и касательным напряжением). После подстановки h_iso = — (ρ_ice/ρ_mantle) H, получаем:

M_load = ∫∫ (r — r_CM) × [ν ρ_ice g (ρ_ice/ρ_mantle) ∇ H] dS = (ν ρ_ice² g / ρ_mantle) ∫∫ (r — r_CM) × ∇ H dS. (6)

Интеграл ∫ (r — r_CM) × ∇ H dS по площади плиты может быть преобразован в интеграл от момента инерции ледовой нагрузки. В случае, когда центр масс льда не совпадает с центром масс плиты, этот интеграл отличен от нуля. Для Антарктиды с её характерным смещением центра льда на 320 км, численная оценка даёт:

|M_load| ≈ (ν ρ_ice² g / ρ_mantle) * A * ΔR * (H_avg / L),

где A ≈ 14×10⁶ км² — площадь плиты, ΔR = 320 км, H_avg = 2,1 км, L = 500 км (характерный размер градиента). Принимая ν ≈ 10⁹ м²/с (что соответствует вязкости 10²¹ Па·с и плотности 3300 кг/м³), получаем |M_load| ≈ 5×10²⁵ Н·м. Для сравнения: M_ridge для Антарктиды оценивается в (1—2) ×10²⁷ Н·м. Таким образом, M_load составляет около 2—5% от M_ridge, что согласуется с порядком величины, полученным из GPS-данных (0,003°/млн лет против ω=0,219°/млн лет, то есть 1,4%, см. ниже).

3.2. Переход к угловой скорости: вклад неотектоники и дрейфа полюса

Из (4) и (6) получаем выражение для угловой скорости:

Ω = (1/γ) M_ridge + (1/γ) M_load ≡ Ω_ridge + Ω_load. (7)

Величина Ω_load может быть вычислена численно с использованием современных моделей ледовой нагрузки. Выполнены расчёты по формуле (6) с использованием грида 20×20 км по данным Bedmap2 [15]. Результат: |Ω_load| = 0,0031 ± 0,0005°/млн лет, направление — север-северо-запад, близкое к Ω_ridge. Отношение |Ω_load|/|Ω_ridge| = 0,014 ± 0,002, то есть около 1,4%.

Теперь необходимо добавить эффект дрейфа полюса вращения плиты в системе мантии. В системе отсчёта, где Ω = ω * e (e — единичный вектор полюса), изменение e (t) задаётся скоростью дрейфа v_pole = |de/dt|. По данным горячих точек [7], v_pole = 0,09 ± 0,03°/млн лет. Полный вектор угловой скорости как функция времени:

Ω (t) = ω (t) e_0 + ω (t) * (v_pole t) × e_0 + Ω_load (t), (8)

где ω (t) = ω_0 exp (-t/τ) — затухание угловой скорости, e_0 — начальное направление полюса. Интегрирование (8) по времени даёт угол поворота плиты и смещение любой точки.

3.3. Анализ чувствительности к вязкости и дегляциации

Главный источник неопределённости — вязкость мантии под Восточной Антарктидой, которая по разным данным варьирует от 1×10²¹ до 5×10²¹ Па·с [5, 24]. Это меняет Ω_load пропорционально ν (см. (6)), т.е. в 5 раз. Для базового расчёта принято ν = 2×10²¹ Па·с. При варьировании ν в указанном диапазоне Ω_load изменяется от 0,0015 до 0,0075°/млн лет, что на 50 млн лет даёт дополнительный вклад в смещение от ±30 км до ±120 км (включено в итоговую неопределённость). Дегляциация (уменьшение H) уменьшает Ω_load со временем; в прогнозе на 50 млн лет мощность льда уменьшается линейно до 0,3 км, что снижает Ω_load примерно в 7 раз. Этот эффект учтён в численном интегрировании.

4.1. Методология двунаправленного моделирования

Для исследования на интервале от —50 млн лет (50 млн лет назад) до +50 млн лет (50 млн лет вперёд) использованы следующие положения:

Уравнения движения (1) — (8) из главы 3 симметричны относительно направления времени при замене знака момента сил, но из-за необратимости дегляциации и диссипации в мантии решение не является строго симметричным. Для ретроспективного расчёта знак времени изменён на противоположный, а мощность льда задана по палеоклиматическим реконструкциям.

Параметры: ω (t) по экспоненциальному закону с τ=78 млн лет; дрейф полюса v_pole = 0,09°/млн лет постоянный; вклад неотектоники Ω_load (t) меняется пропорционально средней мощности льда H (t), которая для прошлого восстановлена по данным об изменении уровня моря и изотопов кислорода [51]: 50 млн лет назад — ледниковый щит отсутствовал (H≈0), 34 млн лет назад — начало формирования, современная мощность 2,1 км.

Для ретроспективного расчёта также учтено, что полюс вращения плиты в прошлом дрейфовал; начальное положение полюса 50 млн лет назад восстановлено обратным дрейфом.

4.2. Ретроспективная верификация: интервал 0—50 млн лет назад

Расчётное положение центра масс Антарктиды для моментов 10, 20, 30, 40 и 50 млн лет назад сравнивалось с палеомагнитными данными (28 точек, сгруппированных в 5 временных срезов). Модельные значения получены с точностью ±0,2°.

Для времени 10 млн лет назад палеомагнитная широта составляет 72,8° ю.ш. с погрешностью ±3,5°, модельная — 72,6° ю.ш., расхождение 0,2°. Для 20 млн лет назад: 70,1±4,0° и 70,6°(расхождение 0,5°). Для 30 млн лет назад: 67,5±4,2° и 68,0°(0,5°). Для 40 млн лет назад: 64,9±4,5° и 65,4°(0,5°). Для 50 млн лет назад: 62,6±5,0° и 63,1°(0,5°). Среднеквадратичное расхождение по широте — 0,46°(≈51 км), что лежит в пределах погрешности палеомагнитных определений (±3—5°). Без учёта дрейфа полюса (v_pole=0) расхождение составило бы 1,3°, с постоянной угловой скоростью — 2,1°. Таким образом, модель ПВС с дрейфом полюса успешно верифицирована на всём 50-миллионном интервале.

4.3. Прогноз на 50 млн лет вперёд (детерминированный и вероятностный)

Прогноз выполнен методом Монте-Карло (1000 реализаций) с варьированием следующих параметров (равномерное распределение): начальная угловая скорость ω₀ в диапазоне 0,212 — 0,226°/млн лет; постоянная затухания τ в 66 — 90 млн лет; скорость дрейфа полюса v_pole в 0,06 — 0,12°/млн лет; вязкость ν в 1×10²¹ — 5×10²¹ Па·с; мощность льда через 50 млн лет в 0,2 — 0,5 км.

Результаты (медианы и 90%-й интервал): смещение центра масс — 720 км (90% доверительный интервал: 620—820 км); азимут — 341°(северо-северо-запад); изменение широты центра масс — 6,5°(с 75,5° ю.ш. до 69,0° ю.ш.); изменение долготы — +8,5°(с 120° в.д. до 128,5° в.д.); положение Антарктического полуострова — 57,1° ю.ш. (интервал 55,5—58,7°). Сравнение с прогнозом без учёта дрейфа полюса (только Ω_ridge + затухание) даёт смещение 640 км, т.е. вклад дрейфа полюса составляет +80 км (около 12,5%).

4.4. Асимметрия геодинамического режима на интервале ±50 млн лет

Сравнение эволюции ключевых параметров в прошлом (от —50 до 0) и будущем (от 0 до +50) выявляет следующую асимметрию.

Угловая скорость ω (t): за 50 млн лет до настоящего ω уменьшилась с 0,34 до 0,219°/млн лет (падение на 36%). За 50 млн лет вперёд ω уменьшится с 0,219 до 0,168°/млн лет (падение на 23%). Причина — нелинейный экспоненциальный закон: скорость убывания пропорциональна текущему значению, поэтому в прошлом, когда ω была выше, абсолютное падение было больше.

Вклад неотектоники Ω_load/Ω_ridge: 50 млн лет назад ледникового щита практически не было (H≈0), поэтому Ω_load ≈ 0. В современную эпоху Ω_load/Ω_ridge = 1,4%. Через 50 млн лет, при H≈0,3 км, это отношение уменьшится до ≈0,2%. Таким образом, неотектонический вклад был нулевым, достиг максимума в настоящем и снова упадёт в будущем.

Положение полюса вращения: за 50 млн лет до настоящего полюс дрейфовал в направлении, обратном современному (азимут 132°), сместившись примерно на 4,5°. В будущем он продолжит дрейф в том же направлении (азимут 312°), также на 4,5°. Симметрия дрейфа полюса относительно настоящего времени соблюдается, так как скорость дрейфа принята постоянной.

Итог: режим не является симметричным из-за эволюции ледовой нагрузки и экспоненциального затухания ω. Прошлое характеризовалось быстрым уменьшением угловой скорости и нарастанием неотектонического вклада, будущее — более плавным уменьшением и затуханием этого вклада. Количественно асимметрия может быть выражена отношением падения ω в прошлом и будущем: (0,34—0,219) / (0,219—0,168) = 0,121/0,051 ≈ 2,37. То есть падение в прошлом было в 2,4 раза круче.

4.5. Климатические и гляциологические следствия на интервале ±50 млн лет

В прошлом, 50 млн лет назад, Антарктида находилась на широтах 62—63° ю.ш., ледниковый щит отсутствовал, континент был покрыт умеренными лесами (данные палинологии). Смещение к полюсу и похолодание привели к оледенению около 34 млн лет назад. В будущем, через 50 млн лет, северные части континента вернутся на широты, сходные с широтами 50-миллионнолетней давности, но климатическая система будет иной из-за более высокого содержания CO₂ (прогнозируемого). Это вызовет деградацию щита, но не его полное исчезновение.