CC BY
33
УДК 550.42
С.Г. Богуславский, Е.Е. Совга, С.И. Казаков
Влияние поля солености на положение верхней границы зоны сероводорода в Черном море
На основе материалов наблюдений оценивается изменчивость составляющих водного бюджета моря. Новые данные о режиме течений в Босфоре позволили уточнить среднюю величину испарения с поверхности Черного моря, которая оказалась равной 292 км3/год.
С помощью математической модели исследуется влияние материкового стока и водообмена через Босфор на горизонт залегания галоклина, градиент солености в нем и величину коэффициента вертикального обмена, от которых зависит положение верхней границы зоны сероводорода.
Описывается механизм зимних заморов при сгонной циркуляции у Южного берега Крыма.
Горизонт залегания верхней границы зоны Н23 является сложной функцией ряда физических и химических процессов. Мы ограничимся изучением основных физических причин, определяющих тенденцию многолетней изменчивости глубины залегания этой границы и ее отклонений от среднего положения.
Отличительной особенностью Черного моря является слабый вертикальный обмен между поверхностной и глубинной водными массами, разделенными галоклином. Максимальная величина градиента солености & галоклина в центральной части моря наблюдается вблизи горизонта 65 м, а на периферии моря — вблизи 80 м.
Основным источником солей Черного моря является нижнебосфорское течение, а стоком — верхнебосфорское. Расход нижнего течения составляет 300 км3/год при средней солености 35,6 %о, а верхнего — 600 км3/год при солености 17,8 %о [1]. Оба эти течения переносят около 340 т/с солей. При этом водный и солевой бюджеты моря сбалансированы. Если принять расход босфорских течений за основу и дополнить их средними величинами стока рек 338 км3/год, осадков 238 км3/год и результирующего поступления вод из Азовского моря 16 км3/го^, то для замыкания водного баланса моря испарение должно быть равно 292 км3/год (слой в 690 мм/год), а не общепринятой величине 396 км3/год [2]. Ошибка объясняется тем, что ранее использовались заниженные величины расходов течений в Босфоре, полученные на основе единичных измерений Макарова и Мерца (нижнее течение 176 км3/год, верхнее - 371 км3/год), а также неточностью современной теории океанического испарения, в которой не учитывается вертикальное распределение влажности в тропосфере над водоемом. Для Аральского моря, окруженного преимущественно пустыней, эта методика дает удовлетворительные результаты. Но при расчете испарения с поверхности океанов или больших морей, окруженных горами, как Черное, получаются завышенные результаты.
Для выяснения влияния современных изменений климата с учетом воздействия антропогенных факторов на режим полей солености и сероводорода
© С.Г. Богуславский, Е.Е. Совга, С.И. Казаков, 2006
целесообразно рассмотреть палеогидрологию Черного моря в плейстоцене (эпоха оледенений) и голоцене, когда климат менялся в широких пределах.
На основе материалов глубоководного бурения дна Черного моря и пролива Босфор [3], сведений об уровне водоема, палеобиологических индикаторов, математического моделирования в работе [4] выполнена палеогидроло-гическая реконструкция Черного моря для эпохи карангатской трансгрессии (рисс - вюрмское межледниковье), новоэвксинской регрессии (вюрм) и голоцена.
В ледниковые периоды уровень моря понижался на 80 — 100 м и в Босфоре оставалось лишь течение из Черного в Средиземное море. С момента вырождения нижнебосфорского течения начинался процесс распреснения моря, который длился около 10 тыс. лет. В конце этого интервала времени море превращалось в пресноводный водоем. Вертикальный обмен зимой захватывал всю толщу вод, сероводород отсутствовал, донные осадки формировались в кислородной среде.
В межледниковые периоды уровень моря восстанавливался до отметок, близких к современной, восстанавливалось нижнебосфорское течение и начиналось осолонение вод моря. В течение около 5 тыс. лет формировались поля солености и сероводорода, близкие к современным [4].
Исходя из вышеизложенного представления, целесообразно рассмотреть реально возможные изменения составляющих водного бюджета моря в XXI в. с учетом антропогенных факторов и оценить их влияние на поля солености и сероводорода. Главной положительной составляющей водного бюджета является сток рек, равный в среднем 388 км3/год. За многолетний период он изменялся от 246 км3 в 1949 г. (73% от среднего) до 492 км3 в 1970 г. (145% от среднего). Несмотря на изъятие на хозяйственные нужды части стока и строительство водохранилищ, в многолетнем ходе суммарного стока отмечается положительный тренд. Иногда сток р. Дунай характеризуется колебаниями с периодами, кратными периодам солнечной активности 22 года и 11 лет.
В среднем на поверхность моря выпадает 238 км3/год осадков. Минимум 170 км3/год наблюдался в 1948 г. (72% от среднего), а максимум 322 км3/год — в 1981 г. (136% от среднего). В многолетнем ходе осадков отмечен положительный тренд с угловым коэффициентом 0,86 км3/год.
В многолетнем ходе испарения выявлен отрицательный тренд с угловым коэффициентом, близким к коэффициенту осадков, но с обратным знаком. Изменение испарения связано с уменьшением скорости ветра.
Возрастание стока рек и осадков, а также уменьшение испарения ведут к распреснению поверхностных вод до галоклина. По данным наблюдений в многоводные периоды понижение солености прибрежных вод достигало в среднем 0,05%о, а вод центральной части моря 0,005%о. В это время ниже галоклина наблюдалось повышение солености [2]. А в галоклине возрастал градиент солености.
В конце XX в. отмечалось некоторое уменьшение толщины верхнего квазиоднородного слоя, что связано также с уменьшением средней скорости ветра и повторяемости штормов [2]. В апреле 1987 г. экспедицией МГИ на
НИС «Академик Вернадский» в центре западной части моря верхняя граница зоны H2S была обнаружена на горизонте 67 м.
По мере приближения галоклина к поверхности моря усиливается влияние ветрового перемешивания на вертикальный профиль солености. В период нагрева вод сезонный термоклин способствует подъему верхней границы галоклина.
На основе вышеизложенного рассмотрим некоторые задачи влияния реально возможных изменений водного бюджета моря на поле солености и, следовательно, на изменение границы зоны H2S.
При этом можно ограничиться одномерной моделью, осредняя изопо-верхности солености по исследуемой акватории. Для этого случая уравнение диффузии запишем в виде
dS , d2S dS — = ——, (1) dt 2 dz2 dz
где S — соленость, t — время, w — вертикальная компонента средней скорости, z — глубина.
Во многих случаях коэффициент вертикального обмена kz определяется путем деления потока тепла или солей на соответствующий градиент. При этом учитывается и перенос, осуществляемый вертикальной компонентой скорости. Поэтому А.Г. Колесниковым было введено понятие коэффициента
суммарного переноса [5] к*:
, dS /d2S к7 - к7 — w—/ —. 2 2 dz / dz
Второй член правой части этой формулы имеет размерность коэффициента диффузии и учитывает перенос вертикальной компонентой скорости.
Далее величину к* мы будем обозначать к2.
Для условий Черного моря зависимость к2 от градиента плотности можно записать в виде
к =_h._ (2)
2 1 + fi(dS/dz) '
где ко — коэффициент суммарного переноса при адиабатическом распределении плотности, /? — постоянная для данного района величина.
Для учета влияния стока рек и водообмена через Босфор на поле солености и величину к2 рассмотрим следующие начально-краевые задачи [6]:
dS_=d_
dt dz
ds
и dz
(3)
/ = 0, S(0, z) =Az), (4)
2 = '' йщщ,'^- (6)
где З^О^) — начальное распределение солености; — функция, которая характеризует временную изменчивость расхода верхнебосфорского течения;
— величина, определяемая расходом нижнебосфорского течения; / — средняя глубина черноморской впадины.
1. Рассмотрим случай, когда в течение 30 лет сток рек сократится на 10%. Такой процесс можно описать уравнением (3) при условиях (4), (5) и (6), где у(г) есть функция, которая убывает и характеризует сокращение стока, а
— величина постоянная. Численное решение системы (3) — (6) для этих условий дает следующие результаты. Сокращение стока рек и соответствующее уменьшение верхнебосфорского течения обусловливает постепенное осоло-нение вод выше галоклина. Одновременно происходит его заглубление.
С помощью формулы (2) можно оценить изменение к2 в течение 30 исследуемых лет. Расчеты показывают, что величина к2 в галоклине увеличивается при этом от 0,3 до 0,9 см2/с при одновременном увеличении глубины минимума к2 от 70 до 250 м. В этих условиях увеличивается поток кислорода через галоклин и верхняя граница зоны Н28 соответственно заглубляется.
2. Влияние более существенного сокращения стока рек исследовалось в работах [4, 6]. Рассмотрим случай сокращения стока рек на одну треть в течение 20 лет с последующей 10-летней его стабилизацией на достигнутом минимальном уровне. Для этих условий изменение величины к2 и горизонта залегания его минимума приведено в таблице. В работах [4,6] оценивается также влияние нижнебосфорского течения на поле солености и величину к2 в течение длительного интервала времени.
Влияние изменения стока рек на изменение величины кг и горизонта залегания его минимума
Изменение стока, % от современного Начальный уровень Конечный уровень
л А, см /с Глубина г, м к: см2/с Глубина г, м
Падение на 10 0,3 70 0,9 250
Падение на 30 0,3 70 5,0 280
Рост на 10 0,3 70 0,2 60
3. Поскольку в последние 20 лет на фоне периодических изменений наблюдается систематический рост речного стока, то целесообразно воспользоваться начально-краевой задачей (3) — (6) для выяснения влияния роста материкового стока на поле солености и величину к2. Рассмотрим случай роста стока на 10% за 30 лет. В этом случае в краевом условии (5) величина у(/) должна равномерно увеличиваться за 30 лет на 10%, а величину будем считать постоянной. Численное решение задачи дает следующие результаты:
градиент солености в галоклине увеличивается; величина к2 уменьшается за 30 лет от 0,3 до 0,2 см2/с; горизонт залегания минимума к2 уменьшается от 70 до 60 м (табл.).
Дальнейший подъем квазистационарного галоклина глубоководной части моря выше горизонта 60 м ограничен ветровым перемешиванием. По этой причине положение верхней границы зоны H2S в ближайшие десятилетия в среднем изменится мало. Но может увеличиться градиент солености в галоклине и концентрация H2S под ним.
Отклонение галоклина от своего среднего положения, обусловленное влиянием ветра, приводит иногда к локальным заморам на шельфе вблизи материкового склона у Южного берега Крыма (ЮБК).
Летом, когда развит сезонный термоклин, при сгонных ветрах вода холодного промежуточного слоя с температурой 8 — 10 °С поднимается к поверхности. При этом поднимается и галоклин и тянет за собой сероводородную зону. Но поскольку для одновременного подъема термоклина и галоклина требуется слишком большая работа ветра, то обычно при летних сгонах на шельфе у ЮБК заморы не наблюдаются. Но в конце зимы, когда термоклин отсутствует, при западных штормовых ветрах галоклин и граница зоны H2S поднимаются на несколько десятков метров и захватывают шельф до изобаты 50 - 60 м. При этом рыбы, которые зимуют в углублениях шельфа, отгороженных от глубоководной части моря порогами, гибнут. Такой зимний замор отмечался в феврале 1954 г. к западу от м. Кикинеиз. При этом происходила массовая гибель султанки и скорены.
В связи с отмеченной выше тенденцией медленного роста материкового стока в Черное море и увеличением градиента солености в галоклине есть основания ожидать в XXI в. увеличения повторяемости зимних заморов на шельфе у ЮБК.
В заключение отметим, что в глубоководной части Черного моря под га-локлином в верхнем слое сероводородной зоны происходит частичная минерализация органического вещества. Высвобождающиеся при этом биогенные элементы частично переносятся в зону фотосинтеза и могут способствовать увеличению биологической продуктивности моря.
СПИСОК ЛИТЕРАТУРЫ
1. Ozsoy Е., LatifM.A., Tugrul S. ai al. Water and nutrient fluxes through the Bosphorus, and exchanges between the Mediterranean and Black Seas // Problems of The Black Sea. — Sevastopol: MHI UAS, 1992. — P. 54 - 67.
2. Практическая экология морских регионов. Черное море /Под ред. В.П. Кеонджяна. — Киев: Наук, думка, 1990. — 250 с.
3. Геологическая история Черного моря по результатам глубоководного бурения /Под ред. Ю.П. Непрочнова. — М: Наука, 1980. — 202 с.
4. Еремеев В.Н., Богуславский С.Г., Жоров В.А. Особенности палеогидрологии Черного моря в различные геологические эпохи // Морской гидрофизический журнал. — 1994. — №4. — С. 28 - 36.
5. Колесников А.Г., Иванова З.С., Богуславский С.Г. О влиянии устойчивости на вертикальный перенос в Атлантическом океане // Океанология. — 1961. — вып. 4. — С. 592 - 599.
6. Березовский A.A., Богуславский С.Г. Математическая модель расчета и прогноза поля солености Черного моря // Асимптотическое интегрирование нелинейных дифференциальных уравнений. — Киев: Ин-т математики HAH Украины, 1992. — С. 133 - 165.
ABSTRACT Variability of the components of the seawater budget is estimated based on observations. New data on the current regime in the Bosporus permit to specify a mean value of evaporation from the Black Sea surface - 292 km3/year.
The mathematical model provides an opportunity to study the influence of the continental run-off and the Bosporus water exchange upon the depth of the halocline location, its salinity and the value of the vertical exchange coefficient that condition the location of the hydrogen sulfide zone upper boundary.
The mechanism of winter fish kill under the surge circulation nearby the Crimean southern coast is described.
Морской гидрофизический институт HAH Украины, Севастополь
Материал поступил в редакцию 29. 11.04
