Математическое моделирование эволюции рельефа сбросового уступа на примере Святоносского поднятия (Байкальская впадина)
Бызов Л.М., Саньков В.А.
В работе представлены результаты комплексного математического моделирования позднекайнозойской эволюции рельефа сбросовых уступов на примере Святоносского поднятия (Байкальская впадина). Моделирование выполнялось в программах CHILD и BFS, включающих параметры тектонических, флювиальных и склоновых процессов. Реконструкции показали высокую степень корреляции с реальным рельефом, что подтверждает состоятельность предложенного нами сценария эволюции. Также были выявлены некоторые особенности развития рельефа территории.
Ключевые слова: комплексное математическое моделирование, эволюция рельефа, сбросовый уступ, морфотектоника, Байкальский рифт, Святоносское поднятие.
Numerical Modeling of Landscape Evolution of NormalFault Scarp by the Example of Svyatoy Nos Upland
L. M. Byzov, V. A. Sankov
Abstract. We present results of complex numerical modeling of Late-Cenozoic landscape evolution of normal-fault scarps by the example of Svyatoy Nos Upland (Baikal Basin). For modeling we used the programs CHILD and BFS (by Greg Tucker) including algorithms of tectonical, fluvial and hillslope processes. Models showed the good correlation with real relief. These results confirm consistency of evolution's scenario which we suggested. Also we defined some features of relief development of territory.
Keywords: complex numerical modeling, normal-fault scarp, CHILD, landscape evolution, Baikal rift, Svyatoy Nos Upland.
Рельеф земной поверхности можно рассматривать как отражение воздействия на нее внешних и внутренних процессов. На размеры и форму рельефа влияют тектонические движения, физическое и химическое выветривание, флювиальные и склоновые процессы, антропогенное воздействие и многие другие факторы. Таким образом, в целях получения максимально полного представления о геоморфологическом строении той или иной территории необходимо рассматривать целый комплекс процессов, участвующих в его формировании. Для решения этой задачи в геоморфологии и геодинамике применяется метод комплексного математического моделирования. Суть его заключается в следующем. Интересующий участок суши или, к примеру, берегового склона принимается в виде пространственной сетки (если речь идет о трехмерных моделях) или двухмерного профиля. Каждый элемент сетки/профиля меняется во времени и пространстве под воздействием рельефообразующих процессов, выраженных в виде математических формул. Методика позволяет реконструировать рельеф, зная количественные параметры формирующих процессов, и напротив, владея информацией о рельефе, исследовать условия, его создавшие.
Математическое моделирование нельзя назвать принципиально новым методом, истоки идеи уходят к работе выдающегося американского геолога Г. К. Гилберта, а большинство практических закономерностей и первые модели были получены во второй половине XX в. (Куллинг, Леопольд и Мэд- док, Анхерт и др.). Революционным явлением для моделирования стало активное внедрение в геоморфологические исследования компьютерных технологий, что спровоцировало появление новых, до того времени невозможных, инструментов. Работы над созданием новых программ и подходов не прекращаются по настоящее время, при этом приветствуется кооперация усилий. В этой связи следует упомянуть сообщество CSDMS, созданное на базе Университета Колорадо и объединяющее исследователей со всего мира. Большинство методов, разработанных участниками сообщества, можно, с разрешения авторов, использовать в научных целях. Инструменты математического моделирования применяются для решения различных геоморфологических и геотектонических задач - например, для реконструкции флю- виальных (работы М. Атталя, А. Виттакера, Д. Брауна) и склоновых (Дж. Роеринг, А. Мартин и М. Чёрч) процессов, сбросовых (К. Пети), складчато-надвиговых (С. Миллер) структур и т. д. Сегодня существует множество компьютерных программ, направленных на комплексное моделирование рельефа; некоторые из них созданы для решения конкретной задачи, другие отличаются универсальностью и многофункциональностью. Однако следует заметить, что, несмотря на стремительное развитие компьютерных технологий, степень объективности математического моделирования зависит прежде всего от объективности исходных данных. И в случаях, когда речь идет об эпохах продолжительностью в миллионы лет, добиться точности всех многочисленных параметров практически невозможно. Таким образом, задача математического моделирования в геоморфологии заключается не в стремлении реконструировать полный аналог реального рельефа, но воссоздать стиль его развития, определить ключевые параметры рельефообразующих процессов и количественно проверить состоятельность той или иной гипотезы.
Обзор методики
В своей работе мы использовали программы CHILD (Channel-Hillslope Integrated Landscape Development) и BFS (Bedrock Fault Scarp), разработанные Грегом Такером [8]. CHILD является одним из наиболее универсальных инструментов геоморфологических трехмерных реконструкций. Программа написана на языке C++ для платформы UNIX (также возможно ее использование в среде Windows через эмулятор Cygwin). Она создает трехмерное изображение земной поверхности для каждого момента времени (шаг задается пользователем). Собственными функциями визуализации программа не располагает, для этих целей применяются инструменты Surfer или MatLab. Принцип моделирования основан на методе конечных разностей. Управление программой осуществляется через вводный текстовый файл, содержащий свыше 150 параметров, которые прямо или опосредованно характеризуют свойства модели. Определение многих параметров требует сложных поэтапных вычислений, с привлечением многочисленных дополнительных данных. Для этих целей нами был разработан оригинальный калькулятор на базе Microsoft Excel. В целом структуру CHILD можно представить в виде нескольких основных групп параметров: 1) климатические характеристики; 2) свойства коренного и осадочного материалов; 3) русловая геометрия; 4) русловые эрозионно-аккумулятивные процессы; 5) склоновый крип; 6) тектонический режим; 7) общие параметры сетки. Суть работы CHILD, как и любой подобной программы, можно представить в виде упрощенного выражения:
математический рельеф геолог
где ?h ?p??t?- изменение высоты h узла p за время t, U - изменение высоты за счет тектоники, E f - изменение высоты за счет русловой эрозии, Ed - изменение высоты за счет склоновой диффузии. Каждый элемент данного уравнения может быть как положительным, так и отрицательным (например, тектоническое поднятие/опускание или эрозия/аккумуляция).
Несмотря на уникальную, в своем роде, многофункциональность, CHILD имеет ряд недостатков, потребовавших внесения некоторых корректив в ход исследования. В частности, при моделировании сбросового уступа невозможно изменить угол падении разлома, который по умолчанию составляет 90°. Таким образом, крылья сброса взаимно смещаются вдоль плоскости, перпендикулярной поверхности Земли, что может не соответствовать реальным условиям. Для решения этой проблемы была использована программа BFS, моделирующая изменения двухмерного профиля склона под воздействием эрозии и смещения по разлому с заданным углом падения. Использование BFS позволило скорректировать угол наклона тектогенных граней (фасет), а также рассчитать количественные значения склонового крипа.
Объект исследования
В качестве объекта моделирования было выбрано Святоносское поднятие, которое представляет собой составную (хоть и обособленную) часть междувпадинной перемычки, разделяющей Южно-Байкальскую и СевероБайкальскую впадины, наряду с Ушканьим и Ольхонским горстами и горстом подводного Академического хребта [1]. Поднятие составляют три горста, кулисообразно сменяющие друг друга по простиранию; с востока его обрамляет широкая тектоническая ступень. Ступень перекрыта маломощным осадочным чехлом, и в центральной части он поднимается выше уровня Байкала, формируя Чивыркуйский перешеек (Мягкая Карга). Святоносское поднятие сложено разновозрастными магматическими (граниты, гра- нодиориты, сиениты) и метаморфическими (кристаллосланцы, гнейсы, мраморы) породами и частично перекрыто современными отложениями. Рельеф п-ва Святой Нос представлен высокими гольцовыми горами, максимальная отметка (1878 м над уровнем моря) приходится на южную кулису, центральная кулиса достигает высоты 1708 м, северная - 1639 м. Главная водораздельная поверхность полуострова в центральной его части представляет собой сравнительно широкое плато, а к периферии приобретает черты острого гребня. Морфология отдельных сегментов поднятия существенно различается. Южная кулиса (большая ее часть) представлена крутыми уступами, треугольными фасетами и глубокими У-образными долинами временных водотоков; на крайнем юге, у м. Нижнее Изголовье, высота и протяженность склонов значительно уменьшается, фасеты пропадают, а долины приобретают черты ручейковых водотоков. Склоны северной и центральной кулис отличаются большей протяженностью и, как следствие, более сложной дренажной сетью. Единого фронта фасет здесь не прослеживается, и рельеф характеризуется высокой степенью расчлененности. В береговой зоне большей части полуострова распространены озерные террасы.
Непосредственным объектом моделирования явился восточный склон южной части Святоносского поднятия (рис. 1).
Рис. 1. Объект исследования: 1, а - обзорная карта Прибайкалья; 1, Ь - п-в Святой Нос (по данным БЯТМ); белыми линиями обозначены границы моделируемых объектов; 1, с и 1, ё - сегменты А и В соответственно (сплошными линиями показаны водоразделы, прерывистыми - русла временных водотоков)
В морфотектоническом отношении это сбросовый уступ, представленный рядом тектогенных граней (фасет и антифасет). Длина и высота склона, высота и уклон фасет на всем протяжении уступа неодинаковы - северная часть отличается большей высотой и длиной склона, а также глубиной вреза речных долин. Морфология тектонических уступов (в частности, треугольных фасет) может служить важным индикатором геодинамиче- ских процессов. Исследованием морфотектоники сбросовых уступов в Байкальской рифтовой системе (БРС) занимался Г. Ф. Уфимцев, подробно описавший структуру и процесс формирования тектогенных граней [3]. В своей работе мы подходим к данной проблеме с другой стороны, пытаясь дать количественную оценку эволюции сбросового уступа с помощью математического моделирования. Подобные эксперименты уже проводились. Большое влияние на наши исследования оказали работы К. Пети, посвященные моделированию эволюции рельефа Прибайкалья [5] и Провинции Бассейнов и Хребтов [4], а также исследования М. Атталя по сбросовым уступам Аппе- нин [7]. Однако целью нашего эксперимента является попытка получить комплексную картину позднекайнозойского развития ландшафта с учетом максимально возможного количества факторов.
Условия моделирования
Технические особенности CHILD таковы, что пространственная неоднородность морфологии и вещественного состава Святоносского поднятия (даже в пределах одного склона) не позволяет воссоздать эволюцию всей структуры в рамках одного эксперимента. В этом случае неизбежны ошибки в реконструкции стока, эрозионно-аккумулятивных процессов, склонового крипа, поэтому для успешного моделирования необходимо разбить дискретную структуру на относительно однородные сегменты. В нашем исследовании были выбраны два участка восточного склона южной части Святого Носа (см. рис. 1, b). Они различаются размерами, морфологией и составом слагающих пород, обладая при этом общими свойствами, характерными для сбросовых уступов. В ходе изучения был протестирован ряд количественных моделей эволюции этих сегментов. Результатом эксперимента явились две реконструкции, отражающие возможный сценарий развития объекта на неотектоническом этапе. Далее приведены ключевые вводные параметры для моделей, а также основные формулы и закономерности, используемые в эксперименте.
Модель A отражает свойства одноименного сегмента (см. рис. 1, с). В качестве исходной поверхности рассматривается прямоугольная ровная наклоненная (3°) сетка 3500x7000 м с размером ячейки 100x100 м (выбор обусловлен размером пикселя в моделях SRTM). Модель охватывает период с позднего плиоцена (3,5 млн лет) по настоящее время, что соответствует наиболее активной фазе неотектонического этапа в БРС [2]. До начала активизации поднятия территория долгое время находилась в условиях тектонического затишья, о чем косвенно свидетельствуют данные трекового датирования по апатитам, полученные нами для Баргузинского хребта [6]. Это объясняет небольшой уклон и сглаженность исходной поверхности, поскольку эрозионные процессы, при практическом отсутствии поднятия, привели к денудации склона. Скорость тектонического поднятия в модели составляет 0,4 мм/год [1].
Климатические условия (интенсивность выпадения осадков, продолжительность и частота событий) рассчитывались по современным данным метеостанции «Баргузин» с условием, что они соответствуют средним за рассматриваемый период. Предпосылкой для такого допущения стали па- леоклиматические реконструкции, согласно которым климат Прибайкалья варьировал от теплого субтропического (как аридного, так и влажного) до холодного климата тундры и ледниковых зон [2]. Таким образом, современный климат региона, характерный для таежно-степной зоны, можно условно рассматривать как среднестатистический.
В качестве слагающего материала в модели принимается среднезернистый гранит с размером отдельностей 5 см [7] и критическим сдвиговым напряжением 36 и 14 Па (для коренной породы и реголита соответственно). Формирование стока определяется правилом Хортона, согласно которому сток начинается без избыточного накопления, что характерно для склонов с тонким реголитовым слоем. Геометрия потока в модели подчиняется законам классической русловой геометрии Леопольда - Мэддока.
Эрозионно-аккумулятивные процессы управляются уравнениями эрозионной (2) и транспортной (3) способности:
где E - интенсивность эрозии (м/год), Ke - коэффициент эрозии (зависит от свойств породы), ??и ?cr?- касательное и критическое сдвиговое напряжение (Па), ? - показатель степени, зависящий от выбора единиц измерения:
математический рельеф геолог
где T - транспортная способность (м/год), Kt - коэффициент переноса (зависит от свойств породы), W - ширина потока (м). Склоновые эрозионно-аккумулятивные процессы, не связанные с формированием русла, определяются уравнением нелинейной диффузии (4):