Исследование полей концентрации взвеси на северо-западном шельфе Черного моря при взмучивании донных осадков движущимся циклоном

Термогидродинамика океана

УДК 561.46:551.585.1

Д.В. Алексеев, В.А. Иванов, Е.В. Иванча, В.В. Фомин, Л.В. Черкесов

Исследование полей концентрации взвеси на северо-западном шельфе Черного моря при взмучивании донных осадков движущимся циклоном

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

Введение. Северо-западный шельф Черного моря испытывает значительные антропогенные нагрузки. Загрязняющие вещества поступают с речными и промышленными стоками, с атмосферными выпадениями и в результате техногенных аварий в акватории моря. Одним из способов самоочищения водных экосистем является переход вредных примесей из водной среды в состав донных осадков [1, 2]. Тяжелые вещества непосредственно вовлекаются в процесс гравитационного осаждения, тогда как легкие предварительно адсорбируются частицами минеральной или органической взвеси. Однако, в свою очередь, донные осадки могут оказаться источником вторичного загрязнения морской среды. В мелководных областях с интенсивной динамикой вод, которая формируется, в частности, под воздействием атмосферных аномалий типа циклонов [3, 4], придонные сдвиговые напряжения приводят к взмучиванию отложений, в результате чего загрязняющие вещества снова переходят во взвешенное состояние.

В работах [5 - 7] рассматривается эрозия донных осадков мезомасштаб-ными циклоническими и антициклоническими вихрями, а также меандри-рующими вдольбереговыми и приливными течениями. Помимо этого представляет интерес определение положения областей подъема взвеси, направления ее переноса и вертикального распределения взвешенного вещества в случае, когда динамика вод формируется под действием перемещающегося циклона. В данной работе такие исследования проводятся для северозападного шельфа Черного моря на основе численной сигма-координатной гидродинамической модели [8], адаптированной к региональным условиям. При задании потоков эрозии и седиментации применяется аппроксимация, хорошо зарекомендовавшая себя при моделировании мутности в эстуариях [7,9].

©Д.В. Алексеев, В.А. Иванов, Е.В. Иванча, В.В. Фомин, Л.В. Черкесов, 2007

1. Математическая постановка задачи. Введем декартову систему координат с осями х, у, г, направленными соответственно на восток, север и вертикально вверх. Для описания динамики морской воды воспользуемся нелинейными уравнениями движения однородной несжимаемой жидкости в приближении гидростатики [10]:

(1и __ ^ + _ д (1( р дх дх

2А,

м

ди дх

ду

'ду + дил дх ду

дг

К

м

ди дг

¿V 1 др д

— + /и + —— =— (к р ду дх

А

1м

гду + дх ду

+

ду

2 А *

2Амду

дг

К

м

ду_ дг

др

+ = о,

ог

ди ду д\у Л — + — + — = 0, дх ду дг

где {и, V, м>} - компоненты скорости вдоль осей х, у, - время;

(1/(11 = д/д1 + ид/дх + уд/ду + ч>д/д1; / - параметр Кориолиса; р - давление;

р- 1015 кг/м3 - плотность воды; g - ускорение силы тяжести; - коэффициент горизонтальной турбулентной вязкости, зависящий от горизонтального разрешения модели Дх, Ду и вычисляемый по формуле Смагоринского [11]:

Ам = 0 Д ДхДу [(Зм/дх)2 + (ду/дх + ди/ду)2¡2 + (ду/ду)2 ]*/2; Км - коэффициент вертикальной турбулентной вязкости, который, в соответствии с теорией Мейлора — Ямады [10], является функцией кинетической энергии турбулентности и макромасштаба турбулентности.

Полагая, что взвешенные донные осадки состоят из способных к сцеплению частиц одного типа, запишем для их концентрации С уравнение переноса — диффузии [7, 9]

¿С А

а(^С) д( А дС

дг

дх

дх

ду

г11

дС

Эу.

ог

К,

дС_ дг

Здесь и'., - скорость оседания взвеси, значения которой обычно лежат в пределах 1(Г7 - 10~3 м/с [9]; Ан и Ки - коэффициенты горизонтальной и вертикальной турбулентной диффузии, в дальнейшем принимаемые постоянными.

Предположение об однородности жидкости позволяет проинтегрировать уравнение (3) по вертикали от свободной поверхности г/(х,у^) до произвольного горизонта г. р~ ра+ gp(т]~z). Считая радиус циклонического образования равным Я, представим атмосферное давление ра в виде

Ра =1

- р0 соэ Ра'

ЯГ ~2Ё

+ ра, г<Я\ г >Л,

где ра - фоновое значение давления, р0 - максимальное отклонение от ра в атмосферном возмущении, г - расстояние от центра движущегося циклона до точки с координатами {х, у}.

На свободной поверхности (г = фс,у,ф должны быть выполнены кинематическое и динамические условия:

УУ

^тЛУ-^^ГМ О)

а также условие отсутствия потока взвеси [7]:

-п.С-Кг

дС дг

= 0.

В (7) ра — 1,2 кг/м3 - плотность воздуха; {)¥х, 1¥у} - компоненты вектора скорости ветра W; СА = к2 (14,56 - 1п[(й^2 + )/и^о212 ~ коэффициент поверхностного трения, рассчитываемый на основе соотношения, предложенного в [12], где 1¥0= 1 м/с, к- 0,4 - постоянная Кармана. В неподвижном циклоне ветер рассчитываем, умножая модуль вектора геоциклострофической скорости ^ на эмпирический коэффициент ¡л =0,7 и учитывая, что направление ветра отклоняется от касательных к изобарам на угол у- 20° (угол втока) против часовой стрелки [13]. Считая, что барическое образование движется поступательно со скоростью с, получаем выражение для скорости приводного ветра:

\¥ = Т(90° + - + с, К =-— + г 2

(Ра) г Г Ра

1/2

Здесь (ра)г - радиальный градиент атмосферного давления (6); вектор г направлен от центра циклона к точке, для которой вычисляется XV; ^соб а -биюЛ

- матрица поворота на угол а.

Т(а) =

К эта соб а у

На дне (2 - -Н(х,у)) равна нулю нормальная составляющая скорости, а касательные напряжения связаны со скоростью квадратичной зависимостью [10]:

дН дН . \ди dv\ „ ( 2 2 V/2 f >

В (10) Со - коэффициент донного трения, который находится по формуле: 2о лГ 9 гДе - расстояние от дна, го = 0,001 м - масштаб шероховатости донной поверхности.

Поток взмученных осадков со дна определяется разностью потоков эрозии ¥е и седиментации ^ [7]:

wsC-KH^ = Fe-Fd. dz

Следуя [7, 9], будем считать Ге и функциями придонных сдвиговых напряжений ть = рС£}(и2 + V2). Чтобы учесть эффект истощения взмучиваемого вещества, введем величину запаса осадков, находящихся в данный момент времени на дне В(х,у, (). Седиментация и эрозия соответственно увеличивают и уменьшают этот запас:

дВ-F F

Поток эрозии отличен от нуля, когда придонные сдвиговые напряжения превышают критическое значение тсс и В не равно нулю:

F. =

МР

Л

т

\ L се у

, если ть>тсе и В > 0; в других случаях,

где Mq - коэффициент, характеризующий интенсивность взмучивания и изменяющийся в пределах КГ6-! 0~2 кг/(м2-с) [9].

Поток седиментации определяется концентрацией у дна С/, и отличен от нуля, когда придонные сдвиговые напряжения меньше критического значения Ted'

Fä =

w,Cb

Tcd J

Tu < Г

cd'

0,

Tu >Tt

cd *

На боковых границах выполняются условия прилипания для скоростей и условие отсутствия потока примеси в направлении внешней нормали п:

{м,у} = 0, =

дп

В начальный момент времени циклон находится за пределами акватории, движения жидкости отсутствуют, а свободная поверхность горизонтальна = Распределение запаса осадков на дне при этом предполагается

горизонтально-однородным В(х,у, 0) = 0,3 кг/м2 [7].

2. Метод решения. В исходных уравнениях (1), (2), (4), (5) и граничных условиях (7), (8), (10), (11) перейдем от декартовой к ачсоординатной системе по формулам:

У=У, <г =

где В = Н+г/. С учетом (16) перепишем (1), (2), (4) и (5) следующим образом (звездочки, обозначающие новые независимые переменные и новые искомые функции, здесь и далее опускаются):

д(иВ) . Вдра дг/ д

——- 4- Аи - /у£> +--— + gD—L = —

дг р дх дх да

Р да

д(уР) д1

+

. _ Рдра „дц д(Кмду Ау + /иР +--— + §Р—- =- —---

Р ду

ду да\Рда

дт] [ д(иР) | 8(уР) ^до)=0 д1 дх ду да

д(СР) д1

+ АС

да

да

Кн дС Р да

+ РЕг

Для оператора переноса используется обозначение

К(р^д{иРср) | д(уРср) х д(охр) дх ду да

Величина со в (19) определяется соотношением

со - \v- ii

дО дт] а— + —- | — V дх дх

г до д

а— + — ду ду

дР дт] а дг дг

В рассматриваемой модели выражения и берутся в виде, пред-

ложенном в [10]:

дх

.. -.Эм] д

Г ¿Н> + ди4 дх ду

ИРу = — дх

Г д\> 5м 4 Эх Эу

+ -

У

<Эх I дх

^ д'

ду

" ду

С учетом (22) граничные условия (7), (8) и (10), (11) записываются так:

Выполнение этих условий требуется соответственно при сг= 0 и а- -1.

Далее, поле горизонтальных скоростей представляется в виде суммы двух составляющих: средней по глубине компоненты и отклонения от нее. После этого уравнения движения (17) - (19) расщепляются на две системы. Из первой определяются осредненные горизонтальные скорости и уровень моря. Из второй - отклонения горизонтальных скоростей от их средних значений и вертикальная скорость.

Численный алгоритм решения указанных систем по времени базируется на применении двухслойных схем. Вначале на интервале времени Д/ = 4 мин решается система для осредненных скоростей и уровня с дробным шагом Д/1 =А//ЛГ, где Аг=20 выбрано в соответствии с критерием устойчивости Куранта для баротропных волн. Затем на том же интервале времени, но с шагом Дг, решаются уравнения для отклонений горизонтальных скоростей от их средних значений, а далее определяется вертикальная скорость из уравнения неразрывности. Операторы переноса аппроксимируются с помощью схемы БирегЬее [14], относящейся к классу 7У£>-схем.

Уравнение переноса — диффузии (20) решается методом расщепления с шагом Дг. При расщеплении отдельно выделяется задача трехмерного переноса и горизонтальной диффузии и одномерная задача диффузии по вертикальной координате. Решение последней задачи выполняется методом прогонки с учетом граничных условий (24), (25).

4 С целью повышения пространственного разрешения в модели реализована процедура вложенных сеток. На первом этапе определяются гидродинамические поля и поле концентрации взвеси на грубой сетке (разрешение

Дх = 7,8 км, Ау - 7,0 км) для всего бассейна Черного моря. На втором этапе проводится расчет на более мелкой сетке (разрешение Ах = 1,6 км, Ау = 1,4 км) только для северо-западной части моря, имеющей одну жидкую границу, параллельную оси х и проходящую по широте м. Тарханкут. Начальные условия, условия на свободной поверхности, на дне и на твердых боковых границах те же, что и на первом этапе. А в качестве условий на жидкой границе используются интерполированные на мелкую сетку данные расчетов на крупной сетке, полученные для каждого расчетного момента времени и соответствующие широте жидкой границы. При этом количество уровней по вертикали т0 в обоих случаях равнялось 25. Для лучшего разрешения по вертикали вблизи дна и свободной поверхности эти уровни сгущались с помощью соотношения

сгт = -1 + (т -1 )ах - (ах - а2)$т\2тх{7п -1)]/зт(2лах),

где ах =1/(/и0 -1); яг = 0,02;/и = 1; 2;...; т0.

3. Анализ результатов расчетов. При проведении численных экспериментов значения критических напряжений для эрозии тсе = 0,13 Н/м2 и седиментации гС£/=0,1 Н/м2 брались из работы [7]. Коэффициент М0, характеризующий интенсивность взмучивания донных осадков, принимался равным 10~5 кг/(м2-с), что обеспечивало согласование по порядку величин полученных результатов с данными измерений концентрации взвеси в северозападной части Черного моря в 25-м и 41-м рейсах НИС «Профессор Водя-ницкий» и в 18-м рейсе НИС «Профессор Колесников» [15]. Особую проблему представлял подбор значений скорости оседания м>х и коэффициента вертикальной диффузии Кн. Предварительное численное интегрирование уравнения (5) с граничными условиями (8), (11), выполненное в предположении горизонтальной однородности полей скорости и концентрации взвеси, а также при отсутствии вертикального адвективного переноса, показало, что адекватного описания распределения мутности с глубиной удается достичь, если Кн > м>л тах{£>Дсг}, где Дет- расстояние между соседними сг-уровнями. Чтобы удовлетворить полученному условию, в расчетах использовались Кн = Ю-4 м2/с и \у5 = Ю~5 м/с. Выбранные значения и коэффициента М0 соответствуют пелито-алевритовым осадкам с типичным размером частиц в несколько микрометров [5]. Коэффициент горизонтальной диффузии Ан полагался равным 10 м2/с.

Атмосферное возмущение имело радиус Я = 500 км, перепад давления между центром и периферией =15 гПа, скорость перемещения с - 5 м/с, что соответствует типичным параметрам циклонов, наблюдаемых в Черноморском регионе. Рассмотрим случай, когда барическое образование пересекает акваторию моря с запада на восток в широтном направлении (рис. 1).

с.ш.

46,8°

46,4(

46,0°-

45,6'

30,0° 30,5° 31,0° 31,5° 32,0° 32,5° 33,0° 33,5° в.д.

Р и с. 1. Батиметрия

Траектория движения центра циклона показана штриховой линией, на которой отмечены его положения в некоторые моменты времени.

Пространственная локализация максимумов концентрации взвеси связана прежде всего с районами взмучивания донных осадков. Положение этих районов меняется в процессе перемещения атмосферного образования и определяется как интенсивностью нестационарных течений, так и их взаимодействием с особенностями морфометрии дна. За время нахождения циклона над акваторией северо-западного шельфа дважды (в 35 - 40 и 55 - 60 ч) происходит смена преобладающих направлений ветра. В эти временные интервалы имеет место уменьшение скоростей течений, и существующие области взмучивания постепенно исчезают. По мере установления нового направления ветра скорости возрастают и возникают новые области эрозии донных осадков, но уже в иных местах шельфа.

На этапе вхождения циклона на акваторию моря наиболее интенсивные течения, обусловленные юго-западными ветрами, возникают вдоль западного побережья. Здесь к 20 ч после начала движения барического образования придонные сдвиговые напряжения начинают превышать критическое для эрозии значение, но протяженные области взмучивания образуются только к 25 ч. Они располагаются у западного побережья южнее 46,2°с.ш. С течением времени данные области расширяются на восток, а их мористые границы к 30 ч практически повторяют контуры изобаты 20 м (рис. 2, а). На расстоянии 1 м от дна горизонтальное распределение взвешенного вещества с концентрацией более 0,1 г/м3 соответствует конфигурации и положению районов эрозии (рис. 2, б). При этом мутность воды увеличивается по мере приближения к мелководью, и концентрация взвеси достигает максимальных значений 10 ISSN 0233-7584. Мор гидрофиз. журн., 2007, N2 1

порядка 100 г/м3 в узкой вдольбереговой полосе, в которой запас донных отложений полностью истощается. Имеющий место перенос взвешенных донных осадков течениями на северо-восток вдоль западного побережья не оказывает влияния на ситуацию севернее 46,2°с.ш., где мутность по всей глубине остается равной нулю. На свободной поверхности области с концентрацией 0,1 г/м3 и более пока не прослеживаются.

б

Р и с. 2. Поля придонных напряжений (а) и концентрации взвеси (г/м3) на расстоянии 1 м от дна (б) в момент времени г = 30 ч

К 35 ч юго-западные ветры передней части циклона с наибольшей силой воздействуют на акваторию Каркинитского залива, вызывая там усиление ISSN 0233-7584. Мор. гидрофиз. журн2007, № 1 И

северо-восточных течений. У западного побережья потоки морской воды постепенно ослабевают и отклоняются к востоку, вследствие чего области интенсивной эрозии здесь существенно уменьшаются и к 40 ч полностью исчезают. Новые области взмучивания появляются сначала у западной оконечности Тендровской косы, а затем и у о. Джарылгач (рис. 3, а). К 45 ч во всей акватории северо-западного шельфа устанавливается циклоническая циркуляция. Поле мутности на расстоянии 1 м от дна отражает положение как новых областей эрозии, так и областей, в которых процесс эрозии к этому времени прекратился (рис. 3, б). В последнем случае концентрация взвеси значительно уменьшается по сравнению с 30 ч и не превышает 65 г/м3. У западного побережья вертикальные потоки донных осадков с концентрацией 0,1-1 г/м3 достигают свободной поверхности, образуя на ней локальные пятна (рис. 3, в).

К 50 ч над акваторией моря уже перемещается тыловая часть циклона с преобладающими северо-западными ветрами, однако основные черты циклонической циркуляции сохраняются до 5$ ч. При этом возле о. Джарылгач прекращается донная эрозия, область взмучивания в районе Тендровской косы существенно уменьшается в размерах и вновь образуются два локальных очага подъема донных осадков у западного побережья северо-западного шельфа Черного моря. В последнем районе адвективный перенос в придонном слое проявляется в вытягивании с течением времени на юго-восток области с содержанием взвеси более 0,1 г/м3.

После 60 ч формируется антициклоническая циркуляция, существующая до момента ухода атмосферного возмущения с северо-западного шельфа. На данном этапе образуются самые обширные за все время исследования области эрозии (рис. 4, а). Прежде всего это происходит в центральной части акватории с глубинами 20 - 40 м, а также в районе Одессы и в Днепро-Бугском лимане. Затем интенсивный переход донных отложений в водную среду прослеживается с юга вдоль Тендровской косы, охватывает Бакальскую банку и узкую восточную часть Каркинитского залива. Распределение взвешенного вещества на расстоянии 1 м от дна становится менее локализованным (рис. 4, б). Сначала объединяются области мутности у Тендровской косы и в центральной части акватории, далее к ним присоединяются области у северного побережья Каркинитского залива и в его восточной части, а позднее - и у западного побережья шельфа. Таким образом, локализованными остаются подъемы взвеси вблизи Одессы и в Днепро-Бугском лимане. В западном регионе процесс эрозии уже не имеет места. Локальные пятна, связанные с выходом мутной воды на поверхность, здесь теперь сливаются в одно. На поверхности у Тендровской косы и о. Джарылгач проявляются следы очагов эрозии, действовавших между 40 и 45 ч (рис. 4, в). Помимо этого становятся заметными области мутности над Одесской и Бакальской банками и в Днеп-ро-Бугском лимане, тогда как поверхность центральной части шельфа остается пока свободной от взвеси.

в

Р и с. 3. Поля придонных напряжений (а\ концентрации взвеси (г/м3) на расстоянии 1 м от дна (б) и на поверхности (в) в момент времени t = 45 ч

с.ш. 46,8°

46,4°

46,0°

45,6°

с.ш.

46,8°

46,4° 46,0° 45,6°

с.ш. 46,8°

46,4°

46,0°

45,6°

в.д.

в

Р и с. 4. Поля придонных напряжений (а), концентрации взвеси (г/м3) на расстоянии 1 м от дна (б) и на поверхности (в) в момент времени / = 70 ч

в

Сокращение площадей районов эрозии, а также ее интенсивности начинается после 70 ч. Следует отметить, что запас донных отложений в центре всех крупных областей взмучивания полностью истощается. На данном этапе перемещения циклона существенным оказывается адвективный перенос. Вдоль западного побережья потоки взвеси направлены на северо-восток, в центральной части акватории - на север, у северного побережья - на восток, в Каркинитском заливе - вдоль Тендровской косы. Однако вынос донных осадков за пределы существующих областей мутности невелик и обнаруживается только в непосредственной близости от них. Более заметно перераспределение взвеси внутри данных областей. Следствием этого является вытягивание изолиний концентрации в направлении движения потоков.

После 80 ч очаги эрозии исчезают, а величины придонных сдвиговых напряжений становятся меньше критического значения для седиментации. По мере ухода циклона и ослабления течений изменения пространственного распределения взвеси начинают определяться преимущественно процессами турбулентной диффузии и осаждения. За счет горизонтальной диффузии площади областей, содержащих взмученное вещество, расширяются, а его концентрация уменьшается (рис. 5, а). Благодаря вертикальной диффузии на свободной поверхности к 200 - 250 ч появляются следы самого глубоководного из очагов эрозии в центральной части акватории (рис. 5, б). Процесс седиментации приводит к новому распределению запаса осадков на дне (рис. 5, в). Размеры областей, где он за счет эрозии существенно истощен, постепенно уменьшаются. Вблизи участков границ этих областей, через которые осуществлялся вынос взмученного вещества водными потоками, запас начинает превышать первоначальное значение 0,3 кг/м2. Так, в центральной части акватории это происходит севернее района взмучивания, вблизи Одессы - восточнее, у Тендровской косы и о. Джарылгач - юго-восточнее. У западного побережья увеличение запаса донных отложений прослеживается на северо-востоке, а также у выдающихся на юго-восток участков границ областей эрозии.

С целью получения более детального представления о распределении взвешенных донных осадков в толще воды рассмотрим изменение с течением времени вертикальных профилей концентрации на трех станциях (рис. 1). В месте расположения первой из них у западного побережья придонные сдвиговые напряжения только между 25 и 30 ч превышают критическое для эрозии значение. Станция №2 у Тендровской косы сначала (с 35 до 40 ч) попадает на периферию области взмучивания, а с 60 до 80 ч - в центр вновь образовавшегося очага. В узкой части Каркинитского залива у станции №3 течения не обладают интенсивностью, достаточной для эрозии. Однако через этот район проходят потоки, переносящие взвесь вдоль южного берега о. Джарылгач на этапе выхода циклона с акватории северо-западного шельфа.

После начала взмучивания вертикальные профили концентрации на станциях №1 и 2 характеризуются возрастанием значений от нулевых до максимальных в тонком придонном слое воды. С течением времени толщина этого слоя на указанных станциях увеличивается до 2 - 3 м, а мутность у дна достигает 60 и 1,5 г/м3 соответственно (кривые 1 на рис. 6, а и б).

в

Р и с. 5. Поля концентрации взвеси (г/м3) на расстоянии 1 м от дна (а), на поверхности (б) и запас осадков на дне (в) в момент времени t = 250 ч

Z, м

1—I—I-1—I-1—I—I—г

О 10 20 30 40 50 С, г/м3 а

z, м

С, г/м3

Z, м

-4 4

-6Н

\3 \

\

-Т-1-1-1-1-1-1-1-г-

0 1 2 3 4 С, г/м3 в

Р и с. 6. Вертикальные профили концентрации взвеси: а - на станции № 1

После ослабления течений продолжающийся подъем взвеси и начавшийся процесс седиментации приводят к более плавному нарастанию концентрации с глубиной и уменьшению ее максимальных величин. Между 50 и 55 ч на станции №1, за счет поступления взвешенного вещества из возникшей в непосредственной близости области эрозии, значение мутности у дна увеличивается с 16 до 32 г/м3 (кривая 2 на рис. 6, а). На станции №2 повторное нарастание концентрации вблизи дна до 140 г/м3 связано с началом второго этапа эрозии. При этом вертикальные профили распределения взвеси снова характеризуются наличием придонных участков резкого увеличения значений. Интенсивное взмучивание на станции №2 приводит к истощению запаса отложений за 5 ч, в результате чего концентрация у дна уменьшается до 27 г/м3 (кривая 2 на рис. 6, б). На станцию №3 заметное поступление взвешенного вещества начинается после 75 ч. До 80 ч это происходит в приповерхностном слое, что обусловливает образование соответствующего максимума 1,3 г/м3 на вертикальном профиле концентрации (кривая 1 на рис. 6, в). К 85 ч рост содержания взвешенных донных осадков имеет место уже во всей толще воды (кривая 2 на рис. 6, в).

По мере ослабления влияния адвективного переноса совместный эффект осаждения и диффузии приводит к тому, что с течением времени вертикальные профили распределения взвеси на всех станциях характеризуются близким к линейному возрастанием мутности с глубиной. При этом на станциях №1 и 2 еще до 120 - 140 ч продолжается увеличение концентрации на поверхности до 6 и 4 г/м3 соответственно, тогда как на дне она уменьшается до

17 и 15 г/м3 (кривые 3 на рис. 6, а и б). На станции №3 мутность растет по всей глубине до 140 ч, достигая 2 г/м3 на поверхности и 4 г/м3 на дне (кривая 3 на рис. 6, в).

Заключение. Сформулируем кратко основные результаты работы. Перемещение циклона с запада на восток над акваторией северо-западного шельфа Черного моря сопровождается последовательным возникновением и исчезновением областей донной эрозии. Под влиянием передней части атмосферного возмущения такие области образуются у западного побережья на глубинах до 20 м. Действие центральной части циклона вызывает эрозию у западной оконечности Тендровской косы и у о. Джарылгач. На этапе выхода барического образования с северо-западного шельфа обширные области взмучивания возникают в центральной части акватории с глубинами 20 — 40 м, в районе Одессы, в Днепро-Бугском лимане, вдоль Тендровской косы и в восточной части Каркинитского залива.

Распределение взвеси на расстоянии 1 м от дна отражает положение действующих и переставших действовать областей эрозии. Влияние нестационарных течений приводит к перераспределению взмученного вещества, что проявляется в вытягивании изолиний концентрации в направлении переноса, а у западного побережья - в образовании выдающегося на юго-восток участка с содержанием взвешенных донных осадков более 0,1 г/м3.

Свободной поверхности взвесь достигает со значительным запаздыванием по времени, и запаздывание тем существеннее, чем в более глубоководных районах происходит эрозия. Так, у западного побережья потоки донных осадков с концентрацией выше 0,1 г/м3 выходят на поверхность через 25 ч после

начала взмучивания, а в центральной части акватории - более чем через 100 ч.

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Цыбань А.В., Мошаров С.А. Биогенная седиментация и ее роль в переносе и депонировании загрязняющих веществ в морских экосистемах // Метеорология и гидрология. -1995.-№ 11.-С. 63-71.

2. Русаков В.Ю., Шимкус К.М., Зернова В.В. и др. Количественный и вещественный состав осаждающегося материала на северо-востоке Черного моря // Океанология. - 2003. - 43, №3.-С. 459-468.

3. Иванов В.А., Коновалов А.В., Манилюк Ю.В. и др. Математическое моделирование сгон-но-нагонных колебаний в Черном море // Метеорология и гидрология. - 1999. - №1. -С. 56-63.

4. Иванов В.А., Коновалов А.В., Черкесов Л.В. Влияние циклонов на изменение уровенной поверхности Азовского и Черного морей // Там же. - 2003. - №4. - С. 73-80.

5. Шапиро Г.И., Акивис Т.М., Пыхов Н.В. и др. Перенос мелкодисперсного осадочного материала мезомасштабными течениями в шельфово-склоновой зоне моря // Океанология. - 2000. - 40. №3. - С. 333-339.

6. Анцыферов С.Ai., Акивис Т.М. Модель транспорта наносов, взвешенных приливным течением в береговой зоне моря // Там же. - 2002. - 42, №4. - С. 502-510.

7. Burchard H., Bolding К., Villareal M. R. Three-dimensional modeling of estuarine turbidity maxima in a tidal estuary // Océan Dynamics. - 2004. - 54, №2. - P. 250-265.

8. Алексеев Д.В., Иванча E.B., Иванов B.A. u др. Моделирование эволюции волновых полей в районе северо-западного шельфа Черного моря при прохождении циклона // Морской гидрофизический журнал. - 2005. - №1. - С. 42-54.

9. Yang Zh., Hamrick I.M. Variational inverse parameter estimation in a cohesive sediment transport model: An adjoint approach // J. Geophys. Res. - 2003. - 108, №C2. - P. 37-1-37-10.

10. Blumberg A. F., Me II or G.L. A description of three dimensional coastal océan circulation model in Three-Dimensional Coast Océan Models // Coas. Estuar. Sci. - 1987. -4. - P. 1-16.

11. Smagorinsky J. General circulation experiments vvith primitive équations. I. The basic experi-ment//Mon. Wea. Rev. - 1963.-91. - P. 99-164.

12. Hsu S.A. A mechanism for the increase of wind stress coefficient with wind speed over water surface: A parametric model // J. Phys. Oceanogr. - 1986. -16. - P. 144-150.

13. Григоркина P.Г., Фукс В.P. Воздействие тайфунов на океан. - JL: Гидрометеоиздат, 1986.-244 с.

14. Sweby P. К. High resolution schemes using flux limiters for hyperbolic conservation laws // J. Numer. Anal. - 1984.-21. - P. 995-1011.

15. Кукушкин А.С., Агафонов E.A., Бурлакова З.П. и др. Изменчивость прозрачности и содержание взвешенного вещества в поверхностном слое северо-западной части Черного моря // Океанология. - 2004. -44, №6. - С. 870-881.

Морской гидрофизический институт HAH Украины, Материал поступил

Севастополь в редакцию 12.07.05

Морское отделение Украинского научно-исследователь- После доработки 23.09.05 ского гидрометеорологического института, Севастополь

ABSTRACT Processes of bottom sediment roiling, transport — diffusion of suspension and its sec-ondary sedimentation on the Black Sea northwestern shelf during a cyclone are studied based on the numerical sigma-coordinate model. Bottom sediments are assumed to consist of the particles of one type. The regions of the most intensive suspension roiling, directions of its transport by fluid flows and also the vertical profiles of concentration of suspended bottom sediments in a number of regions are determined. The investigation reveals that distribution on the suspended matter at a distance of 1 m from the bottom reflects the position of active and inactive erosion centers.