Океанология, 2022, T. 62, № 1, стр. 125-134

Абразия берега, сложенного рыхлым материалом

И. О. Леонтьев *

Институт океанологии им. П.П. Ширшова РАН
Москва, Россия

* E-mail: igor.leontiev@gmail.com

Поступила в редакцию 03.12.2020
После доработки 03.03.2021
Принята к публикации 08.04.2021

Полный текст (PDF)

Аннотация

Рассматривается абразия берегов, при разрушении которых образуются наносы волнового поля (песок, гравий, галька). Развитие таких берегов контролируется обратной связью между темпом абразии и объемом растущего пляжа. Предложенная эволюционная модель, приводит к выводу, что абразия должна экспоненциально затухать со временем. Длительность процесса может измеряться десятками, сотнями и тысячами лет в зависимости от доли пляжеобразующего материала, поступающего из клифа. Повышение относительного уровня моря увеличивает длительность процесса. Для прогноза годового объема абразии клифа применен подход Ларсона [15], принимающий во внимание разность высот пляжа и волнового заплеска. Расчетная формула содержит эмпирический коэффициент, оцененный на основе данных по берегам Камского водохранилища и Охотского моря (где важную роль в абразии клифа играет прилив). По расчетам, даже небольшой пляж вызывает существенное снижение объема абразии клифа. Абразия подводного склона определяется на основе показателей надводной абразии и концепции “глубины замыкания”.

Ключевые слова: объем абразии, отступание клифа, развитие пляжа, эволюция абразионного берега

ВВЕДЕНИЕ

Абразия – процесс разрушения берега штормовыми волнами, сопровождающийся формированием абразионной террасы (бенча) и абразионного клифа [2]. Термин “абразия” в широком смысле применяется не только к берегам, выработанным в коренных породах, но и к аккумулятивным по своему генезису побережьям, которые в силу изменений природной обстановки стали испытывать рецессию [8]. По отношению к последним употребляется также термин “размыв”. По сути, речь идет об одном и том же явлении, когда волны разрушают клиф, а его материал выносится на пляж и подводный склон.

Возможны два сценария, отражающие неоднозначную роль обломочного материала, образующегося в ходе абразии: 1) поступающие из клифа обломки бомбардируют его, увеличивая скорость абразии; 2) скапливающийся материал формирует пляж, прикрывающий основание клифа и замедляющий абразию [1, 10, 11]. Вероятность того или иного сценария зависит от свойств материала, из которого построен берег. Первый сценарий реализуется в условиях прочных скальных и полускальных пород, фрагменты которых обладают значительным ударным действием. При уменьшении объема таких обломков (например, за счет истирания и выноса) абразия затухает [1].

В настоящей работе рассматриваются берега, развивающиеся по второму сценарию, когда в теле клифа в той или иной мере присутствуют частицы наносов волнового поля. Имеются в виду рыхлые несцементированные отложения песчаного, песчано-илистого, песчано-гравийного, песчано-галечного или гравийно-галечного состава, супеси, слабо уплотненные суглинки с вкраплениями песчаных и более крупных частиц (породы V класса по классификации [8]). Берега, о которых идет речь, распространены не только в морях (например, в Охотском, где абразию испытывают пересыпи северо-восточного Сахалина, рис. 1), но и в крупных водохранилищах, где еще продолжается процесс выработки профиля равновесия. Поступающий из клифа материал пляжеобразующих фракций отчасти выносится на подводный склон, наращивая осадочный слой на коренном бенче, а отчасти остается на пляже. Скорость отступания клифа оказывается обратно пропорциональной объему накоплений у его подножья [11].

Рис. 1.

Абразионный берег северо-восточного Сахалина (фото В.В. Афанасьева).

Представления о дальнейшей эволюции абразионного берега неоднозначны. Согласно исследованию [1], возможно установление равновесия между поступлением материала из клифа и его уходом (главным образом, за счет истирания). В этом случае должна поддерживаться некоторая оптимальная скорость абразии, которая при неизменности уровня моря и волнового климата будет продолжаться неограниченное время. Вместе с тем, нельзя не учитывать, что по мере расширения абразионной террасы и возрастающей диссипации энергии волн, их воздействие на клиф будет уменьшаться. В результате абразия может затухать и в конце концов прекратиться [2].

Прежние количественные методы оценки абразии разрабатывались в основном в связи с прогнозированием волновой переработки берегов искусственных водохранилищ. Выработанные рекомендации опирались, главным образом, на инженерно-геологические и геоморфологические аналогии, а также на среднегодовые показатели энергии волн в данном бассейне [3, 4, 9]. При этом пляж не рассматривался в качестве активного фактора. Однако и в более современных подходах обратная связь между абразией и развитием пляжа не принимается в расчет в должной мере [1, 5, 12].

Настоящая работа имеет целью в какой-то мере восполнить этот пробел. Здесь предлагается способ расчета годового объема абразии, учитывающий параметры пляжа, а также строится модель эволюции абразионного берега в длительной перспективе, когда свойства пляжа со временем изменяются.

АБРАЗИЯ В МАСШТАБЕ ГОДА

Представим себе отступающий берег, сложенный слабо устойчивыми к размыву отложениями, на котором перед абразионным клифом выстой hcl сформировался пляж с возвышением zb (рис. 2а). Для определения объема абразии воспользуемся моделью Ларсона [15], описывающей размыв авандюны. В ходе этого процесса на фронтальном склоне песчаной дюны также образуется клиф, подобный клифу абразионного берега, выработанному в рыхлых породах. Отличие только в том, что первый из них сложен несвязными частицами, а материал второго включает долю частиц, в той или иной мере сцепленных друг с другом. Это отличие влияет на скорость разрушения клифа, но не на суть самого процесса. В обоих случаях можно допустить, следуя Ларсону, что вес смытого материала ΔP (в нашем случае, материала клифа) пропорционален силе F, с которой бор, движущийся по поверхности пляжа (по сути, фронт потока заплеска) ударяет в клиф, т.е. ΔP ~ F.

Рис. 2.

Схемы к расчету абразии: (а) – абразионный клиф с формирующимся пляжем; (б) – эволюция абразионного берега; (в) – отложенный слой наносов и изменения уровня моря. Обозначения в тексте.

Детальное описание модели приведено в работе [15], а здесь ограничимся кратким пересказом основных моментов. Сила воздействия бора представляется в форме $F = m\left( {\frac{{du}}{{dt}}} \right) \approx {{mu} \mathord{\left/ {\vphantom {{mu} T}} \right. \kern-0em} T}$, где m – масса бора, u – его скорость, t – время, T – период волн. Переходя от веса к объему размыва ΔV и используя ряд аппроксимаций в отношении массы бора m, можно прийти к соотношению

${{\Delta V} \mathord{\left/ {\vphantom {{\Delta V} {\Delta t}}} \right. \kern-0em} {\Delta t}} \sim {{{{u}^{4}}} \mathord{\left/ {\vphantom {{{{u}^{4}}} {{{g}^{2}}T,}}} \right. \kern-0em} {{{g}^{2}}T,}}$
где Δt – длительность волнового воздействия, g – ускорение силы тяжести. Скорость бора аппроксимируется как ${{u}^{2}} = 2g(R - {{z}_{b}})$, где R – высота волнового заплеска (рис. 2а). Далее, полученное соотношение для объема размыва ΔV представляется в дифференциальной форме и интегрируется по времени. Результатом является формула, которая в нашем случае выражает объем абразии Va на единицу длины берега (м3 м–1) за промежуток времени tw:
(1)
${{V}_{a}} = {{K}_{a}}{{(R + \eta - {{z}_{b}})}^{2}}\frac{{{{t}_{w}}}}{T},$
Здесь Ka – безразмерный коэффициент пропорциональности, зависящий от свойств материала клифа (для песчаных наносов Ka порядка 10–3 [15]), η – изменение среднего уровня моря за счет штормового нагона и прилива. Очевидно, абразия возможна только при условии (рис. 2а)
(2)
$R + \eta > {{z}_{b}}.$
Высота волнового заплеска R может рассчитываться по известной формуле Стокдон и др. [16], основанной на натурных измерениях, в которую, однако, рекомендуется включить зависимость от угла подхода волн. Ясно, что при фронтальном подходе волны генерируют более высокий заплеск, чем при распространении под значительным углом Θ0 относительно береговой нормали (индекс “0” относится к глубокой воде). Наряду с другими прибрежными феноменами, волновой заплеск, в конечном счете, должен контролироваться потоком волновой энергии к берегу, прямо пропорциональным произведению $H_{0}^{2}\cos {{\Theta }_{0}}$, где H0 – высота волн. Так как высота заплеска R пропорциональна $H_{0}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}$ [16], то допуская связь R с потоком энергии, можно ожидать зависимость R также и от величины ${{\cos }^{{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-0em} 4}}}}{{\Theta }_{0}}$. Соответствующие изменения R при углах подхода волн Θ0 < 45° не превышают 9% и практически мало заметны. Однако с приближением Θ0 к 90° высота заплеска будет быстро уменьшаться. Таким образом, с учетом высказанных соображений формулу Стокдон и др. можно представить в следующем виде:
(3)
$\begin{gathered} R = (0.385\beta + 0.55\sqrt {0.563{{\beta }^{2}} + 0.004} ) \times \\ \times \,\,\sqrt {{{H}_{{s0}}}{{L}_{0}}} {{\cos }^{{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-0em} 4}}}}{{\Theta }_{0}},\,\,\,\,{{L}_{0}} = \left( {\frac{g}{{2\pi }}} \right){{T}^{2}}, \\ \end{gathered} $
где β – уклон дна, Hs0, L0, и T – значительная высота волн, их длина и период пика спектра соответственно. Для пологих берегов с уклоном 0.01–0.02 множитель перед $\sqrt {{{H}_{{s0}}}{{L}_{0}}} $ принимает значение около 0.04, а при сравнительно больших уклонах (около 0.1) увеличивается до 0.08. В дальнейшем под величиной β подразумеваем средний уклон надводного пляжа, $\beta = {{{{z}_{b}}} \mathord{\left/ {\vphantom {{{{z}_{b}}} {{{l}_{b}}}}} \right. \kern-0em} {{{l}_{b}}}}$, где lb – его ширина (рис. 2а).

Главной составляющей штормового нагона является подъем уровня воды ηw, вызванный касательным напряжением ветра τw, которое должно уравновешиваться наклоном водной поверхности:

(4)
$\frac{{d{{\eta }_{w}}}}{{dx}} = \frac{{{{\tau }_{w}}}}{{\rho g(h + {{\eta }_{w}})}},\,\,\,\,{{\tau }_{w}} = {{k}_{w}}\rho {{W}^{2}}\cos {{\Theta }_{0}},$
где h и x – глубина и горизонтальное расстояние, ρ – плотность воды, ${{k}_{w}} = 3 \times {{10}^{{ - 6}}}$ ветровой коэффициент, W – скорость ветра. Уравнение (4) линеаризуется и его решение представляется в виде трансцендентного соотношения, из которого определяется ${{\eta }_{w}}(x)$ для заданного профиля глубин h(x) в области разгона ветра [6]. С увеличением скорости ветра протяженность области нагона и его высота у берега возрастают.

Допустим, нам известна годовая повторяемость (продолжительность) различных градаций высот волн (и соответствующих периодов) для основных направлений их подхода. Тогда годовой объем абразии VA можно оценить суммированием элементарных величин Vaij, подсчитанных для отдельных градаций высот волн (i) и для выделенных направлений ( j):

(5)
${{V}_{A}} = {{\Sigma }_{j}}{{\Sigma }_{i}}{{V}_{{aij}}}.$
Величина VA характеризует абразию надводной части берегового профиля с возвышением, равным высоте клифа hcl (рис. 2б). Допуская, что форма профиля сохраняется, и все его точки отступают с одинаковой скоростью, можно записать
(6)
${{V}_{A}} = {{h}_{{cl}}}{{X}_{A}},\,\,\,\,{{C}_{A}} = \frac{{{{X}_{A}}}}{{{{t}_{Y}}}} = \frac{{{{V}_{A}}}}{{{{h}_{{cl}}}{{t}_{Y}}}},$
где XA – годовое отступание берега и клифа, CA – соответствующая скорость рецессии, а промежуток времени tY равен 1 году.

Абразия и рецессия надводной части профиля сопровождается соответствующими изменениями подводного склона. Это схематически отражено на рис. 2б, где показаны положения берегового профиля в последовательные моменты времени t. Абразия подводного склона начинается на некоторой глубине ${{h}_{*}}$, где волновое воздействие становится достаточно интенсивным, чтобы деформировать дно. При отступании клифа на расстояние X со дна должен быть вынесен определенный объем вещества. При условии, что форма подводного профиля существенно не меняется, этот объем можно оценить как $\frac{1}{2}X{{h}_{*}}$ (рис. 2б). Соответственно объем подводной абразии ΩA и полный объем абразии за год VA + ΩA определятся как

(7)
${{\Omega }_{A}} = \frac{1}{2}{{X}_{A}}{{h}_{*}} = \frac{1}{2}{{V}_{A}}\frac{{{{h}_{*}}}}{{{{h}_{{cl}}}}},\,\,\,\,{{V}_{A}} + {{\Omega }_{A}} = {{V}_{A}}\left( {1 + \frac{1}{2}\frac{{{{h}_{*}}}}{{{{h}_{{cl}}}}}} \right),$
где принято во внимание соотношение (6).

Глубина замыкания (closure depth) ${{h}_{*}}$ для песчаного берега определяется через высоту волн Hs12h, превышение которой возможно не более 12 ч в году [14]: ${{h}_{*}} = {{a}_{*}}{{H}_{{s12h}}}$, где ${{a}_{*}} \approx $ 1.6. Можно допустить возможность аналогичного определения ${{h}_{*}}$ и для рассматриваемого типа абразионного берега. Правда, в случае связного донного материала коэффициент ${{a}_{*}}$ должен уменьшаться. С другой стороны, глубина ${{h}_{*}}$ едва ли может быть меньше глубины обрушения hB, где воздействие волн на дно максимально. Высота обрушающихся волн HB и глубина hB связаны между собой как HB ≈ 0.8h. Тогда, полагая ${{H}_{B}} \approx {{H}_{{s12h}}}$, определим диапазон возможных значений глубины замыкания следующим образом:

(8)
$1.25{{H}_{{s12h}}} < {{h}_{*}} < 1.6{{H}_{{s12h}}}.$
Например, при ${{H}_{{s12h}}}$ = 4 м глубина ${{h}_{*}}$ может быть в интервале от 5.0 до 6.4 м.

КАЛИБРОВКА ФОРМУЛЫ АБРАЗИИ

Для оценки коэффициента Ka в формуле (1) следует сравнить результаты расчетов с наблюдениями за отступанием абразионных берегов при тех или иных условиях. В данном случае для сравнения выбраны два береговых объекта, один из которых относится к южному берегу Камского водохранилища (район пос. Висим), а второй расположен на побережье Охотского моря в северо-восточной части Сахалина (район залива Пильтун). Профили берегов показаны на рис. 3.

Рис. 3.

Профили береговых участков, использованные для калибровки модели: (а) – Камское водохранилище (пос. Висим); (б) – Охотское море (северо-восточное побережье Сахалина в районе залива Пильтун).

Камское водохранилище. Согласно [13], коренной берег в районе пос. Висим сложен суглинками с супесчаными и песчаными прослойками. К клифу высотой около 3 м (рис. 3а) примыкает песчаный пляж шириной 1.5 м с высотой zb ≈ 0.1 м (соответствующий уклон пляжа β ≈ 0.07). Участок открыт для воздействия ветров С, СЗ и З румбов с разгонами от 11 до 18 км и средней глубиной по разгону около 10 м.

Геолого-геоморфологические условия исключают возможность естественного питания данного сегмента с соседних береговых участков. Поступающие продукты размыва коренного берега вовлекаются во вдольбереговое движение и удаляются из береговой зоны. Очередной год начинается практически с того же самого профиля, отступившего вглубь суши. Скорость отступания клифа составляет 4–5 м в год.

Волновой климат характеризует табл. 1 [13], где представлены параметры волн (Hs0, T, Θ0) при различных скоростях ветра W для главных направлений (С, СЗ и З), а также максимальная годовая продолжительность выделенных волновых ситуаций (tw). Кроме того, в таблицу добавлены высоты штормового нагона ηw и волнового заплеска R. Расчеты нагона основывались на условном профиле дна с глубинами 8–10 м, длиной, соответствующей разгону ветра и уклоном дна 0.05 в прибрежной части (на глубинах менее 8 м).

Таблица 1.  

Режим ветра, волнения и уровня воды у южного берега Камского водохранилища (пос. Висим)

Румб Θ0, град. W, м/с Hs0, м T, с tw, ч ηw, м R, м
С –38.6 5 0.34 2.52 129 0.12
7 0.53 3.07 12 0.02 0.19
9 0.72 3.48 9 0.02 0.25
14 1.15 4.18 3 0.07 0.38
СЗ   6.4 5 0.33 2.47 144 0.13
7 0.51 3.00 15 0.02 0.19
9 0.70 3.42 3 0.03 0.26
14 1.12 4.12 3 0.09 0.39
З 51.4 5 0.38 2.73 252 0.13
7 0.60 3.32 57 0.02 0.21
9 0.80 3.72 6 0.03 0.27
14 1.23 4.36 3 0.09 0.39
19 1.52 4.61 3 0.15 0.46

Судя по результатам проведенных вычислений на базе табл. 1, абразия становится заметной при скоростях ветра не меньше 7 м/с и высоте волн более 0.5 м. Расчетные скорости отступания клифа CA попадают в диапазон наблюдаемых значений при величине коэффициента Ka = 0.003. В этом случае CA = 4.3 м/год, а годовой объем абразии составляет VA = 12.9 м3 м–1.

Несмотря на малый размер, пляж в рассмотренном случае играет весьма важную роль, так как защищает клиф от воздействия относительно небольших волн (0.3–0.4 м высотой), которые имеют сравнительно большую годовую продолжительность (см. табл. 1). Как показывает расчет, при полном отсутствии пляжа действие только этих волн могло бы вызвать дополнительное отступание клифа на 12 м в год.

Охотское море. Данные по сахалинскому побережью любезно предоставлены зав. Лабораторией береговых геосистем Института морской геологии и геофизики ДВО РАН д.г.н. В.В. Афанасьевым. На рассматриваемом участке имеет место абразия голоценовой террасы, сложенной преимущественно песками, а также супесчано-глинистым материалом с вкраплениями гравия и гальки. Перед клифом высотой 5 м расположен песчаный пляж с максимальным возвышением zb = 1.8 м и средним уклоном β = 0.04 (рис. 3б). Берег испытывает воздействие ветров и волн СВ, В и ЮВ румбов, а также приливов (с амплитудой около 0.6 м). Клиф отступает со скоростью 1.5–2 м в год.

Ветро-волновые условия характеризует табл. 2, основанная на волновых данных, приведенных в [7]. Здесь также приведены расчетные высоты волнового заплеска R и ветрового нагона ηw. При расчетах нагона принималось во внимание, что область его развития при наиболее сильных ветрах примерно соответствует расстоянию до 100-метровой изобаты (в данном случае от 40 до 70 км в зависимости от направления). Также учитывались расстояния до 50 и 10-метровой изобат и характерное значение уклона дна 0.01 от глубины 10 м до берега.

Таблица 2.  

Режим ветра, волнения и уровня воды у северо-восточного побережья о. Сахалин

W, м/с Hs0, м T, с СВ
Θ0 = 35°
В
Θ0 = –10°
ЮВ
Θ0 = –55°
tw, ч ηw, м R, м tw, ч ηw, м R, м tw, ч ηw, м R, м
5 0.25 7.0 4 0.22 18 0.23 70 0.21
7 0.75 8.1 44 0.45 166 0.01 0.47 788 0.41
8 1.25 11.2 70 0.01 0.80 745 0.02 0.84 2681 0.01 0.73
9 1.75 11.3 140 0.02 0.96 1060 0.03 1.00 1656 0.02 0.88
10 2.25 10.0 158 0.03 0.96 210 0.045 1.01 210 0.03 0.88
11 2.75 10.0 79 0.04 1.06 79 0.07 1.11 70 0.05 0.97
13 3.25 10.3 35 0.05 1.19 44 0.09 1.25 26 0.07 1.09
15 3.75 10.6 18 0.06 1.31 26 0.12 1.38 18 0.09 1.21
17 4.25 11.0 7 0.09 1.45 18 0.15 1.52 10 0.12 1.33
19 4.75 11.4 4 0.12 1.59 11 0.19 1.67 4 0.15 1.46
21 5.25 11.9 3 0.15 1.75 10 0.23 1.83 4 0.18 1.60
23 5.75 12.3 1 0.19 1.89 6 0.27 1.98 3 0.21 1.73
25 6.25 12.6 1 0.23 2.02 4 0.31 2.11 1 0.24 1.85
27 6.75 12.2       4 0.35 2.13 1 0.28 1.86
29 7.25 12.5       2 0.40 2.26      
31 7.75 13.0       2 0.45 2.43      

При включении в расчет приливных колебаний допускался их синусоидальный характер, при котором среднее возвышение уровня $\bar {a}$ в той или иной фазе связано с амплитудой прилива am как $\bar {a} = \left( {\frac{2}{\pi }} \right){{a}_{m}}$. В рассматриваемом случае am ≈ 0.6 м, а значит, среднее возвышение $\bar {a}$ ≈ 0.4 м. Далее предполагалось, что время действия tw каждой волновой ситуации, отраженной в табл. 2, условно делится пополам между фазами высокого ($ + \bar {a}$) и низкого ($ - \bar {a}$) уровня. Это подразумевает, что результирующий годовой объем абразии VA может быть определен полу-суммой величин, подсчитанных при высоком ($V_{A}^{ + }$) и низком уровнях ($V_{A}^{ - }$), т.е. ${{V}_{A}} = \frac{1}{2}(V_{A}^{ + } + V_{A}^{ - })$.

Результаты расчетов показывают, что на фоне отлива волновой заплеск в сумме с нагоном едва превышает верхнюю отметку пляжа только при самых сильных штормах. Поэтому абразия берега связана в основном с положительной фазой прилива. При этом скорость ветра и высота волн должны составлять не менее 15 м/с и 3.75 м для восточного румба и 17 м/с и 4.25 м для остальных румбов. Иначе говоря, абразия реализуется только при достаточно сильных штормах. Расчетные скорости отступания клифа CA попадают в диапазон наблюдаемых значений при величине коэффициента Ka = 0.002. Соответствующие значения CA и VA в этом случае составляют 2.0 м год–1 и 10.0 м3 м–1.

Как видно из табл. 2, уровень энергетических воздействий на берега северо-восточного Сахалина весьма высок. Поэтому сравнительно небольшой темп абразии клифа в рассмотренном примере может поддерживаться только благодаря значительному объему пляжа.

Таким образом, по результатам проведенной калибровки можно сделать вывод, что для клифов, сложенных рыхлыми отложениями (суглинки, супеси, пески) коэффициент абразии Ka близок к значениям 0.002–0.003, которые могут послужить ориентиром для расчетных оценок в аналогичных условиях.

АБРАЗИЯ В ДОЛГОВРЕМЕННОМ МАСШТАБЕ

При рассмотрении долговременной эволюции абразионного берега будем предполагать, что она контролируется двумя основными факторами – развитием пляжа и изменениями относительного уровня моря.

Допустим, что материал клифа содержит некоторое количество твердых частиц, способных удержаться у берега и сформировать пляж, причем их доля составляет nb < 1. Определенная часть поступающей из клифа пляжеобразующей фракции уходит на подводный склон, другая часть выносится из береговой зоны (например, в результате истирания), а на самом пляже остается только часть kb < 1. Объем пляжа Vb с течением времени должен возрастать за счет поступления материала из клифа пропорционально темпу абразии VA/tY и, следовательно, можно записать

(9)
$\frac{{d{{V}_{b}}}}{{dt}} = {{k}_{b}}{{n}_{b}}\frac{{{{V}_{A}}}}{{{{t}_{Y}}}}.$
Правую часть (9) можно представить с учетом (6) как ${{k}_{b}}{{n}_{b}}{{h}_{{cl}}}\frac{{{{X}_{A}}}}{{{{t}_{Y}}}} = \frac{{\Delta {{h}_{b}}}}{{{{h}_{{cl}}}}}\frac{{{{V}_{A}}}}{{{{t}_{Y}}}}$, где величина
(10)
$\Delta {{h}_{b}} = {{k}_{b}}{{n}_{b}}{{h}_{{cl}}}$
характеризует условную среднюю толщину слоя наносов, отложенного за год в пределах длины годового отступания клифа XA (рис. 2в). Тогда (9) перепишется в виде
(11)
$\frac{{d{{V}_{b}}}}{{dt}} = \frac{{\Delta {{h}_{b}}}}{{{{h}_{{cl}}}}}\frac{{{{V}_{A}}}}{{{{t}_{Y}}}}.$
Теперь допустим, что уровень моря изменяется за год на величину Δζ. Если уровень повышается (Δζ > 0), то фактическая толщина слоя осадков Δhb уменьшается на Δζ, а при понижении – увеличивается на ту же величину (рис. 2в). Следовательно, в условиях изменения относительного уровня уравнение (11) трансформируется в
(12)
$\frac{{d{{V}_{b}}}}{{dt}} = \frac{{\Delta {{h}_{b}} - \Delta \zeta }}{{{{h}_{{cl}}}}}\frac{{{{V}_{A}}}}{{{{t}_{Y}}}}.$
Заметим, что при изменении уровня моря высота клифа также должна меняться (рис. 2в). Однако для упрощения задачи будем допускать, что изменение уровня за рассматриваемый период времени t мало по сравнению с высотой клифа ($\Delta \zeta \frac{t}{{{{t}_{Y}}}} \ll {{h}_{{cl}}}$), и величину hcl в первом приближении можно считать постоянной.

По мере роста пляжа объем абразии должен уменьшаться, что отражает соотношение

(13)
$\frac{{d{{V}_{A}}}}{{dt}} = - \frac{{d{{V}_{b}}}}{{dt}}$
или с учетом (12)
(14)
$\frac{{d{{V}_{A}}}}{{dt}} = - \frac{{\Delta {{h}_{b}} - \Delta \zeta }}{{{{h}_{{cl}}}}}\frac{{{{V}_{A}}}}{{{{t}_{Y}}}}.$
Решение данного уравнения записывается в виде
(15)
${{V}_{A}} = {{V}_{{A0}}}{{e}^{{ - \frac{{\Delta {{h}_{b}} - \Delta \zeta }}{{{{h}_{{cl}}}}}\frac{t}{{{{t}_{Y}}}}}}},$
где VA0 – годовой объем абразии в начальный момент t = 0.

Как видно, развитие процесса абразии определяется соотношением величин Δhb и Δζ. Если пляж растет быстрее, чем уровень моря (Δhb > Δζ), то абразия экспоненциально затухает со временем. При падении уровня (Δζ < 0) затухание убыстряется. Если нарастание пляжа компенсируется подъемом уровня (Δhb = Δζ), то темп абразии остается постоянным. Наконец, в случае очень быстрого повышения уровня (Δζ > Δhb) абразия будет возрастать.

Чтобы оценить вероятность отмеченных трендов, зададимся характерными значениями определяющих параметров. Пусть, например, высота клифа равна hcl = 5 м, а содержание пляжеобразующего материала в нем составляет nb = 0.1, причем на пляже остается лишь десятая часть поступающего за год объема, т.е. kb = 0.1. Тогда отложенный за год слой наносов составит Δhb = 0.05 м. В то же время годовые изменения относительного уровня моря, вызванные глобальными изменениями климата или тектоникой, оцениваются величинами порядка $\Delta \zeta = {{10}^{{ - 3}}}$ м, что на порядок меньше Δhb. Следовательно, фактор нарастания пляжа доминирует, обусловливая постепенное затухание процесса абразии. Сценарии стабилизации или увеличения абразии при повышении уровня могут реализоваться только при ничтожно малом поступлении пляжеобразующего материала из клифа. Такая ситуация возможна, но не типична для берегов, рассматриваемых в данной работе.

Для оценки длительности процесса абразии Ta можно задаться некоторым минимальным значением VA, например, 0.01VA0, и рассматривать его как практический предел, соответствующий окончанию процесса. Тогда при t = Ta из (15) следует

(16)
${{T}_{a}} = - \ln 0.01\frac{{{{h}_{{cl}}}}}{{\Delta {{h}_{b}} - \Delta \zeta }}{{t}_{Y}} = 4.6\frac{{{{h}_{{cl}}}}}{{\Delta {{h}_{b}} - \Delta \zeta }}{{t}_{Y}}.$
При постоянном уровне моря Ta = $ = 4.6\frac{{{{h}_{{cl}}}}}{{\Delta {{h}_{b}}}}{{t}_{Y}} = \frac{{4.6}}{{{{k}_{b}}{{n}_{b}}}}{{t}_{Y}}$. График данной зависимости показан на рис. 4а. При очень малом содержании наносов в теле клифа (${{n}_{b}} \ll 1$) и преобладающем выносе их с пляжа (${{k}_{b}} \ll 1$) абразия может продолжаться сотни и тысячи лет. Если же доля наносов, остающихся на пляже составляет порядка десятой части всего объема абразии (kbnb = 0.1), то время Ta сокращается до нескольких десятков лет.

Рис. 4.

Зависимости, следующие из эволюционной модели абразии: (а) – длительность процесса абразии в зависимости от доли пляжеобразующего материала kbnb, поступающей из клифа и остающейся на пляже; (б) – изменения объемов пляжа и абразии с течением времени.

Чтобы найти суммарный объем абразии VΣ за определенный период, следует проинтегрировать величину VA/tY по времени: ${{V}_{\Sigma }} = \int {({{{{V}_{A}}} \mathord{\left/ {\vphantom {{{{V}_{A}}} {{{t}_{Y}}}}} \right. \kern-0em} {{{t}_{Y}}}}){\kern 1pt} dt} $. Полагая, что VΣ = 0 в момент t = 0, получим из (15)

(17)
${{V}_{\Sigma }} = \frac{{{{h}_{{cl}}}}}{{\Delta {{h}_{b}} - \Delta \zeta }}{{V}_{{A0}}}\left( {1 - {{e}^{{ - \frac{{\Delta {{h}_{b}} - \Delta \zeta }}{{{{h}_{{cl}}}}}\frac{t}{{{{t}_{Y}}}}}}}} \right).$
Очевидно, предельный объем абразии равен ${{V}_{{\Sigma m}}} = \frac{{{{h}_{{cl}}}}}{{\Delta {{h}_{b}} - \Delta \zeta }}{{V}_{{A0}}}$. С учетом (16) его можно выразить через время абразии Ta:
(18)
${{V}_{{\Sigma m}}} = 0.22{{V}_{{A0}}}\frac{{{{T}_{a}}}}{{{{t}_{Y}}}}.$
Максимальное расстояние Xm, на которое может отступить клиф постоянной высоты hcl, составляет
(19)
${{X}_{m}} = {{{{V}_{{\Sigma m}}}} \mathord{\left/ {\vphantom {{{{V}_{{\Sigma m}}}} {{{h}_{{cl}}}.}}} \right. \kern-0em} {{{h}_{{cl}}}.}}$
Подстановка значения VA из формулы (15) в уравнение (12) и интегрирование последнего (при условии Vb = 0 в момент t = 0) позволяют определить изменения объема пляжа в ходе абразии:
(20)
${{V}_{b}} = \frac{{\Delta {{h}_{b}}}}{{\Delta {{h}_{b}} - \Delta \zeta }}{{V}_{{A0}}}\left( {1 - {{e}^{{ - \frac{{\Delta {{h}_{b}} - \Delta \zeta }}{{{{h}_{{cl}}}}}\frac{t}{{{{t}_{Y}}}}}}}} \right).$
Максимальный объем пляжа, как видно, стремится к значению ${{V}_{{bm}}} = \frac{{\Delta {{h}_{b}}}}{{\Delta {{h}_{b}} - \Delta \zeta }}{{V}_{{A0}}}$. В условиях постоянного уровня моря величина Vbm соответствует годовому объему абразии в начальный момент: Vbm = VA0.

Из сравнения (17) и (20) видно, что объем пляжа и суммарный объем абразии изменяются со временем одинаковым образом, и отношение этих величин остается постоянным: $\frac{{{{V}_{b}}}}{{{{V}_{\Sigma }}}} = \frac{{\Delta {{h}_{b}}}}{{{{h}_{{cl}}}}} = {{k}_{b}}{{n}_{b}}$. Зависимость объемов от времени показана на рис. 4б. Величины Vb и VΣ здесь отнесены к их максимальным значениям Vbm и VΣm, а время – к продолжительности процесса абразии Ta.

Проиллюстрируем полученные результаты числовым примером. Как и прежде, примем hcl = = 5 м, nb = 0.1 и kb = 0.1, что соответствует отложенной за год толщине слоя наносов Δhb = 0.05 м. Допустим, что относительный уровень моря повышается, и годовой прирост составляет Δζ = = 0.001 м. В этом случае абразия будет продолжаться в течение времени Ta = 469 лет. При начальном годовом объеме абразии VA0 = 10 м3 м–1 ее суммарный объем за время Ta составит ${{V}_{{\Sigma m}}} = {{10}^{3}}$ м3 м–1. Берег отступит на расстояние Xm = 200 м, а объем пляжа на завершающей стадии процесса приблизится к величине Vbm = 10 м3 м–1. При понижении уровня моря с той же скоростью (Δζ = –0.001 м/год) показатели Ta, VΣm, Xm и Vbm уменьшатся, но всего на 2%.

Теперь представим, что поступление наносов из клифа на пляж на порядок больше, чем в приведенном примере, т.е. kbnb = 0.1 и Δhb = 0.5 м. Тогда время абразии и ее суммарный объем составят Ta = 46 лет и ${{V}_{{\Sigma m}}} = {{10}^{2}}$ м3 м–1, т.е. станут на порядок меньше.

ЗАКЛЮЧЕНИЕ

Проведенное исследование нацелено на прогнозирование абразии берегов, при разрушении которых образуются наносы волнового поля (песок, гравий, галька). По мере их накопления формируется пляж, поглощающий энергию волн и уменьшающий воздействие на абразионный клиф. Ключевым фактором развития таких берегов является обратная связь между объемом абразии и объемом растущего пляжа. Предложенная эволюционная модель приводит к оптимистическому выводу, что процесс абразии должен экспоненциально затухать со временем. Длительность процесса зависит от концентрации пляжеобразующих наносов в теле клифа, а также от доли материала, остающейся на пляже в результате перераспределения на береговом склоне и безвозвратных потерь. При неизменных внешних условиях абразия берегов рассматриваемого типа может продолжаться от нескольких десятков до нескольких сотен лет. После этого клиф отмирает, а берег переходит в категорию стабильного аккумулятивного. Разумеется, возможны особые случаи, когда в силу свойств прибрежной динамики практически весь поступающий материал уносится с пляжа, и замедления абразии не происходит, как, например, на рассмотренном ранее береговом участке Камского водохранилища.

Повышение относительного уровня моря, связанное с тектоникой или глобальным изменением климата, увеличивает длительность процесса. При типичных скоростях изменений уровня, наблюдаемых в настоящее время, этот фактор не слишком заметен. Но он может оказаться определяющим в условиях, когда на какой-либо стадии процесса подпитка пляжа по тем или иным причинам уменьшится или прекратится.

Полученные соотношения позволяют оценить объем абразии клифа, объем пляжа и темп отступания берега в требуемый момент времени, если заданы параметры процесса и известен начальный годовой объем абразии клифа. Последний подсчитывается суммированием элементарных объемов абразии Va при характерных волновых ситуациях, действующих в течение года.

Для оценки величины Va в данной работе применен подход Ларсона с соавторами [15], согласно которому клиф разрушается под действием силы, возникающей при ударе прибойного потока и зависящей от разности высот заплеска и пляжа. Расчетная формула содержит эмпирический коэффициент Ka, для оценки которого использовались данные по участкам берегов Камского водохранилища и Охотского моря. В последнем случае выявлено значительное влияние прилива на величину абразии. Сделан вывод, что в условиях рыхлого материала клифа типа суглинков, супесей и песков значение Ka близко к 0.002–0.003. В дальнейшем планируется уточнить величину коэффициента на основе дополнительных тестов для более широкого спектра условий. Проведенные расчеты показывают, что даже небольшой пляж способен снизить темп абразии клифа в несколько раз.

В заключение отметим, что рецессия надводной части профиля сопровождается соответствующей абразией подводного склона, которая может быть оценена на основе показателей надводной абразии при использовании концепции “глубины замыкания”. Эта тема нуждается в дальнейшей проработке.

Благодарности. Автор благодарен заведующему Лабораторией береговых геосистем Ин-та морской геологии и геофизики ДВО РАН д.г.н. В.В. Афанасьеву за предоставленные данные о береговых участках северо-восточного Сахалина.

Источник финансирования. Работа выполнена в рамках государственного задания ИО РАН (тема № 0128-2021-0004) при частичной поддержке РФФИ (грант № 18-55-34002 Куба_т).

Список литературы

  1. Есин Н.В., Савин М.Т., Жиляев А.П. Абразионный процесс на морском берегу. Л.: Гидрометеоиздат, 1980. 200 с.

  2. Зенкович В.П. Основы учения о развитии морских берегов. М.: Изд-во АН СССР, 1962. 710 с.

  3. Золотарев Г.С. Инженерно-геологическое изучение береговых склонов водохранилищ и оценка их переработки // Труды лаборатории гидрогеологических проблем. М.: Изд-во АН СССР. 1955. Т. 7. С. 38–57.

  4. Качугин Е.Г. Рекомендации по изучению переработки берегов водохранилищ // Всесоюзный институт гидрогеологии и инженерной геологии (ВСЕГИНГЕО). М.: Госгеолтехиздат, 1959. 89 с.

  5. Леонтьев И.О. Количественный анализ эволюции профиля абразионного берега // Проблемы развития морских берегов. М.: Изд-во ИОАН, 1989. С. 30–37.

  6. Леонтьев И.О. Прибрежная динамика: волны, течения, потоки наносов. М.: ГЕОС, 2001. 272 с.

  7. Леонтьев И.О. Афанасьев В.В. Динамика лагунного берега северо-восточного Сахалина на примере системы Ныйского залива и косы Пластун // Океанология. 2016. Т. 56. № 4. С. 618–626. https://doi.org/10.7868/S0030157416030138

  8. Леонтьев О.К. Основы геоморфологии морских берегов. М.: Изд-во МГУ, 1961. 418 с.

  9. Попов Б.А. Расчет абразионного профиля равновесия подводных склонов крупных водохранилищ // Изв. АН СССР. Сер. геогр. 1969. № 3. С. 89–101.

  10. Сафьянов Г.А. Геоморфология морских берегов. М.: Изд-во МГУ, 1996. 406 с.

  11. Шуйский Ю.Д. Проблемы исследования баланса наносов в береговой зоне моря. Л.: Гидрометеоиздат, 1986. 240 с.

  12. Хабидов А.Ш., Леонтьев И.О., Марусин К.В. и др. Управление состоянием берегов водохранилищ. Новосибирск: Изд-во СО РАН, 2009. 239 с.

  13. Хабидов А.Ш., Леонтьев И.О., Марусин К.В. и др. Мониторинг береговой зоны внутренних водоемов России. Новосибирск: Изд-во СО РАН, 2012. 139 с.

  14. Birkemeier W.A. Field data on seaward limit of profile change // J. of Waterway, Port, Coastal and Ocean Engineering. 1985. V. 111. No. 3. P. 598–602.

  15. Larson M., Erikson L., Hanson H. An analytical model to predict dune erosion due to wave impact // Coastal Engineering. 2004. V. 51. P. 675–696.

  16. Stockdon H.F., Holman R.A., Howd P.A., Sallenger A.H. Empirical parameterization of setup, swash, and runup // Coastal Engineering. 2006. V. 53. P. 573–588.

Дополнительные материалы отсутствуют.