Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона

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

Ausführliche Beschreibung

Gespeichert in:

Bibliographische Detailangaben
Datum:	2013
Hauptverfasser:	Михайлова, Э.Н., Шапиро, Н.Б.
Format:	Artikel
Sprache:	Russian
Veröffentlicht:	Морський гідрофізичний інститут НАН України 2013
Schriftenreihe:	Морской гидрофизический журнал
Schlagworte:	Математическое моделирование морских систем
Online Zugang:	https://nasplib.isofts.kiev.ua/handle/123456789/105100
Tags:	Tag hinzufügen Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:	Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона / Э.Н. Михайлова, Н.Б. Шапиро // Морской гидрофизический журнал. — 2013. — № 5. — С. 64-81. — Бібліогр.: 10 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	nasplib_isofts_kiev_ua-123456789-105100
record_format	dspace
spelling	nasplib_isofts_kiev_ua-123456789-1051002025-02-23T19:24:16Z Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона Михайлова, Э.Н. Шапиро, Н.Б. Математическое моделирование морских систем По результатам экспедиции выявлено, что особенностью субмаринной разгрузки в сентябре 2011 г. было существование четко выраженного промежуточного слоя пониженной солености около слоя резкого скачка температуры, которые наблюдались в открытой части моря вблизи и на удалении (до 60 м) от берега, причем на достаточно большом протяжении (несколько сотен метров) параллельно береговой линии. На основании проведенных численных экспериментов установлено, что главным фактором, влияющим на формирование такой структуры полей температуры и солености, является наличие источников пресной воды, расположенных глубже слоя скачка температуры в рассматриваемый период времени. Показано, что сравнительно тонкий слой пониженной солености и резкий термоклин могут быть обусловлены одновременным действием субмаринной разгрузки и нагонного ветра. За результатами експедиції виявлено, що особливістю субмаринного розвантаження у вересні 2011 р. було існування чітко вираженого проміжного шару зниженої солоності біля шару різкого стрибка температури, які спостерігалися у відкритій частині моря поблизу і на видаленні (до 60 м) від берега, причому на достатньо великому протязі (декілька сотень метрів) паралельно береговій лінії. На підставі проведених чисельних експериментів встановлено, що головним фактором, який впливає на формування такої структури полів температури та солоності, є наявність джерел прісної води, розташованих глибше за шар стрибка температури в даний період часу. Показано, що порівняно тонкий шар зниженої солоності та різкий термоклин можуть бути обумовлені одночасною дією субмаринного розвантаження та нагонного вітру. Results of the expedition show that the feature of submarine discharge in September, 2011 was existence of the well-pronounced intermediate layer with reduced salinity nearby the layer of sharp temperature drop. Both of them were observed in the open sea close to the coast and up to 60 km off it for rather a long space (a few hundreds meters) along the coastline. Based on the performed numerical experiments it is found that the basic factor influencing formation of such a structure of temperature and salinity fields consists in presence of freshwater sources located below the temperature drop layer during the period under consideration. It is shown that rather a thin layer with reduced salinity and sharp thermocline can be conditioned by simultaneous action of submarine discharge and surge wind. 2013 Article Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона / Э.Н. Михайлова, Н.Б. Шапиро // Морской гидрофизический журнал. — 2013. — № 5. — С. 64-81. — Бібліогр.: 10 назв. — рос. 0233-7584 https://nasplib.isofts.kiev.ua/handle/123456789/105100 551.465 ru Морской гидрофизический журнал application/pdf Морський гідрофізичний інститут НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	Russian
topic	Математическое моделирование морских систем Математическое моделирование морских систем
spellingShingle	Математическое моделирование морских систем Математическое моделирование морских систем Михайлова, Э.Н. Шапиро, Н.Б. Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона Морской гидрофизический журнал
description	По результатам экспедиции выявлено, что особенностью субмаринной разгрузки в сентябре 2011 г. было существование четко выраженного промежуточного слоя пониженной солености около слоя резкого скачка температуры, которые наблюдались в открытой части моря вблизи и на удалении (до 60 м) от берега, причем на достаточно большом протяжении (несколько сотен метров) параллельно береговой линии. На основании проведенных численных экспериментов установлено, что главным фактором, влияющим на формирование такой структуры полей температуры и солености, является наличие источников пресной воды, расположенных глубже слоя скачка температуры в рассматриваемый период времени. Показано, что сравнительно тонкий слой пониженной солености и резкий термоклин могут быть обусловлены одновременным действием субмаринной разгрузки и нагонного ветра.
format	Article
author	Михайлова, Э.Н. Шапиро, Н.Б.
author_facet	Михайлова, Э.Н. Шапиро, Н.Б.
author_sort	Михайлова, Э.Н.
title	Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона
title_short	Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона
title_full	Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона
title_fullStr	Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона
title_full_unstemmed	Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона
title_sort	моделирование особенностей субмаринной разгрузки пресных вод в районе м. айя в период нагона
publisher	Морський гідрофізичний інститут НАН України
publishDate	2013
topic_facet	Математическое моделирование морских систем
url	https://nasplib.isofts.kiev.ua/handle/123456789/105100
citation_txt	Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона / Э.Н. Михайлова, Н.Б. Шапиро // Морской гидрофизический журнал. — 2013. — № 5. — С. 64-81. — Бібліогр.: 10 назв. — рос.
series	Морской гидрофизический журнал
work_keys_str_mv	AT mihajlovaén modelirovanieosobennostejsubmarinnojrazgruzkipresnyhvodvrajonemajâvperiodnagona AT šapironb modelirovanieosobennostejsubmarinnojrazgruzkipresnyhvodvrajonemajâvperiodnagona
first_indexed	2025-11-24T15:45:42Z
last_indexed	2025-11-24T15:45:42Z
_version_	1849687165310599168
fulltext	© Э.Н. Михайлова, Н.Б. Шапиро, 2013 Математическое моделирование морских систем УДК 551.465 Э.Н. Михайлова, Н.Б. Шапиро Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона По результатам экспедиции выявлено, что особенностью субмаринной разгрузки в сентяб- ре 2011 г. было существование четко выраженного промежуточного слоя пониженной солено- сти около слоя резкого скачка температуры, которые наблюдались в открытой части моря вблизи и на удалении (до 60 м) от берега, причем на достаточно большом протяжении (не- сколько сотен метров) параллельно береговой линии. На основании проведенных численных экспериментов установлено, что главным факто- ром, влияющим на формирование такой структуры полей температуры и солености, является наличие источников пресной воды, расположенных глубже слоя скачка температуры в рас- сматриваемый период времени. Показано, что сравнительно тонкий слой пониженной солено- сти и резкий термоклин могут быть обусловлены одновременным действием субмаринной разгрузки и нагонного ветра. Ключевые слова: субмаринная разгрузка, негидростатическая двумерная модель, числен- ный эксперимент. Введение. Настоящая статья является логическим продолжением работ [1 − 6] по моделированию субмаринной разгрузки пресных вод в районе м. Айя на Южном берегу Крыма. Предыдущие работы были посвящены вос- произведению особенностей термохалинной структуры вод в самой большой карстовой полости во время гидрологических съемок в экспедициях в сен- тябре 2007 и 2008 гг. Эти съемки проводились в различных синоптических ситуациях, которые приводили к различной стратификации в открытой части моря (фоновой стратификации). Так, в 2007 г. имел место сгон, и температура воды в верхнем слое моря была достаточно низкой (около 14°C). В 2008 г. температура воды в море была сравнительно высокой (около 24°C) и близкой к типичной температуре для этого времени года. При этом температура вте- кающей через границы карстовой полости пресной воды по косвенным при- знакам (по температуре поверхности моря в полости) была выше фоновой температуры воды в море во время первой экспедиции и ниже фоновой тем- пературы воды в море во время второй экспедиции. В расчетах температура втекающей пресной воды принималась равной 15°C. Соленость втекающей воды принималась равной 0‰. Соленость в открытом море (фоновая соле- ность) во время проведения этих экспедиций увеличивалась с глубиной при- мерно с 17,8 до 18,2‰. Во время экспедиции в сентябре 2011 г. была другая синоптическая си- туация, по-видимому, имел место нагон, и фоновая стратификация имела со- ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 64 вершенно другой вид. А именно в поле температуры наблюдался резкий слой скачка: в верхнем 8-метровом слое температура практически не менялась с глубиной и равнялась 22°C, а ниже резко уменьшалась до 9°C и не менялась с глубиной до дна. Фоновая соленость при этом не менялась по глубине и со- ставляла 17,8‰. Заметим, что экспедиция проходила примерно в течение од- них суток, поэтому трудно судить о продолжительности такой ситуации. Наблюдения во время этой экспедиции, как и во время двух предыдущих, проводились не только в пределах наибольшей карстовой полости (грот № 1 по терминологии работ [1, 4]), но и в открытой части моря, прилегающей к системе гротов (карстовых полостей) на побережье м. Айя [7]. Эти наблюде- ния показали, что субмаринная разгрузка проявляется не только в наличии приповерхностного распресненного слоя в гроте № 1 и его ближайшей окре- стности, как было раньше, но и в появлении вне грота промежуточного слоя пониженной солености около слоя скачка температуры. Интересно, что этот подповерхностный слой наблюдался вблизи и вдали от берега (до 60 м), при- чем на достаточно большом протяжении (несколько сотен метров) парал- лельно береговой линии. Представляется, что, помимо источников пресной воды, действующих непосредственно в гроте № 1 и расположенных на сравнительно небольших глубинах (глубина карстовой полости не превышает 8 м), имеются источники пресной воды, расположенные вне гротов на большей глубине, по крайней мере глубже 9 м – глубины слоя скачка температуры. И именно поступающая из них пресная вода, перемешиваясь с морской водой, распространяется вблизи термоклина. Эта идея была в частном порядке высказана руководите- лем экспедиции А.И. Чепыженко. Данная работа как раз и посвящена моделированию субмаринной раз- грузки в случае, когда существует четко выраженный термоклин, а источник пресной воды находится на достаточно большой глубине ниже слоя скачка плотности. Следуя работе [5], моделирование будем проводить в рамках дву- мерной негидростатической модели с использованием приближения «твердой крышки». Выполнен цикл численных экспериментов для исследования механизмов формирования и эволюции особенностей стратификации в прибрежной зоне в районах, где существуют достаточно интенсивные глубоководные источники субмаринной разгрузки. Исследование носит качественный характер, по- скольку нет достоверных данных о реальных параметрах расходов и глуби- нах залегания источников подземных вод. Постановка задачи. Будем рассматривать движение в области 0 ≤ y ≤ B, 0 ≤ z ≤ H, предполагая, что все функции не зависят от координаты x. Здесь ось X направлена вдоль прямолинейного берега, ось Y – по нормали к берегу, ось Z – вертикально вниз, система координат – левая. При этом y = 0 – это берег, y = B – это открытая граница, z = 0 – это невозмущенная поверхность моря, z = H = const – это горизонтальное, для простоты, дно. Исходные уравнения запишем в приближении Буссинеска и в традици- онном приближении для силы Кориолиса, причем в дивергентном виде: ut + (vu)y + (wu)z – fv = Auzz + ALuyy, (1) ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 65 vt + (vv)y + (wv)z + fu = – py + Avzz + ALvyy, (2) wt + (vw)y + (ww)z = – pz + R + Awzz + ALwyy, (3) vy + wz = 0, (4) Tt + (vT)y + (wT)z = µ Tzz + µ LTyy, (5) St + (vS)y + (wS)z = µ Szz + µ LSyy, (6) R = R(T, S). (7) Здесь u, v, w − компоненты скорости течения вдоль осей x, y, z соответствен- но; T – температура; S – соленость; R = g ρ / ρ 0, где g – ускорение силы тяже- сти, ρ – плотность, ρ 0 = 1 г/см3; p − деленное на ρ 0 давление; f – параметр Кориолиса; (A, AL) и (µ , µ L) – коэффициенты вертикальной и горизонталь- ной турбулентной вязкости и диффузии соответственно. Параметр Кориолиса из-за малости масштабов считаем постоянным. Соотношение (7) – уравнение состояния для морской воды, плотность рассчитывается по формуле Мамаева ρ = 1,028154 – [(0,0735 + 0,00469 T) + (0,002 T – 0,802) (S – 35)]. (8) Рассмотрим постановку граничных и начальных условий. Поскольку на- блюдаемая фоновая структура вод является нетипичной и, по-видимому, кратковременной, то представляется разумным рассмотреть начальную ста- дию процесса перестройки термохалинной структуры вод под действием субмаринной разгрузки, не пытаясь получить стационарное решение. На поверхности моря, не учитывая потоки массы, тепла и соли и исполь- зуя приближение «твердой крышки», получим при z = 0 Auz = –τ x, Avz = –τ y, w = 0, Tz = Sz = 0, (9) где τ x, τ y – составляющие касательного напряжения трения ветра. На дне принимаем условия прилипания и отсутствия потоков тепла и соли: при z = H u = v = w = 0, Tz = Sz = 0. (10) На берегу задается втекание воды в слое Z1 ≤ z ≤ Z2 со скоростью v0 (с априори заданным полным потоком Q0), температурой T0 и соленостью S0. Таким образом, имеем при y = 0 и z ∉ [Z1, Z2] u = v = w = 0, Ty = Sy = 0, (11) при y = 0 и z ∈ [Z1, Z2] u = w = 0, v = v0, (vT – µ LTy) = v0T0, (vS – µ LSy) = v0S0. (12) На открытой границе ставятся условия свободного протекания для ско- рости течения и так называемые условия транспортивности для потоков теп- ла и соли: при y = B uy = vy = 0, w = 0, (13) (vT – µ LTy) = v+T + v−Tfon, (vS – µ LSy) = v+S + v−Sfon, ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 66 где v = v+ + v−; v+ = (v + v )/2 ≥ 0 – скорость вытекающего, v− = (v – v )/2 ≤ 0 – скорость втекающего в область течения. Условия транспортивности отражают тот факт, что в рассматриваемую область вода втекает с фоновыми температурой и соленостью, а вытекает со своей температурой и своей соленостью. Отметим, что условия для темпера- туры и солености на берегу (12) представляют собой частный случай условий (13). Условия свободного протекания для скорости позволяют рассчитать на открытой границе компоненты скорости течения на отдельных горизонтах. Эти условия представляются естественными с физической точки зрения, од- нако с математической точки зрения, как показывает опыт расчетов, они яв- ляются проблематичными, так как способствуют неустойчивости численной схемы. В начальный момент времени задаются поля температуры и солености, а именно фоновые поля, и предполагается отсутствие движения, т. е. при t = 0 u = v = w = 0, T = Tfon, S = Sfon. (14) Подчеркнем, что на открытой границе также задаются фоновые распре- деления температуры и солености, которые не меняются по времени. Использование приближения «твердой крышки» обусловлено, в первую очередь, тем, что оно обеспечивает сохранение количества воды в области. В этом случае не происходит накопления воды – сколько воды втекает из ис- точника, столько вытекает через открытую границу. Это следует из уравне- ния неразрывности. Интегрируя уравнение (4) по вертикали и учитывая ра- венство нулю вертикальной скорости на поверхности и дне моря, получим Vy = 0, где V = ∫ H vdz 0 – полный поток, причем V = const = Q0. Разумеется, уровень моря ζ не равен тождественно нулю и может быть вычислен апостериори из условия pz=0 = pa – ζρ 0=zg , (15) следующего из гидростатического соотношения, справедливого для тонкого верхнего слоя (0, ζ ), и граничного условия pz=ζ = pa, где pa – атмосферное давление. При этом нужно помнить, что уровень рассчитывается с точностью до константы, которая определяется, например, из условия сохранения объе- ма воды – равенства нулю среднего по площади уровня: dy B ∫ 0 ζ = 0. Для решения поставленной задачи используется уравнение вихря ω = = wy − vz, которое получается исключением давления с помощью перекрест- ного дифференцирования уравнений движения (2), (3) и с учетом уравнения неразрывности (4) и имеет вид ω t + vω y + wω z – fuz = Ry + Aω zz + ALω yy. (16) Уравнение (4), в свою очередь, позволяет ввести функцию тока ψ , так что ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 67 v = − ψ z, w = ψ y, (17) ω = ∆ψ = ψ yy + ψ zz. (18) Алгоритм решения задачи состоит в следующем. Сначала, зная распре- деления всех полей в момент времени t = n∆ t, из уравнений (5), (6) рассчи- тываются температура и соленость в следующий момент времени t = (n+1)∆ t. Затем с помощью уравнения состояния (7) вычисляется плот- ность R. Используя рассчитанные поля T, S, R, далее с помощью уравнения (16) рассчитывается вихрь ω . После этого решается уравнение Пуассона (18), определяется функция тока ψ и находятся компоненты скорости течения v, w. Наконец, из уравнения (1) вычисляется компонента скорости u и т. д. В принципе в каждый момент времени из уравнения (2) могут быть найдены давление p при z = 0 и уровень моря ζ . Граничные условия для вихря ω и функции тока ψ следуют из условий (9) – (13). Так, условия для вихря имеют вид: при z = 0 ω = τ y/A, (19) при z = H ω = – zv , (20) где при конечно-разностной аппроксимации производной vz предполагается линейная интерполяция условия v = 0 при z = H, причем v берется за n-й мо- мент времени. Аналогично используется конечно-разностная аппроксимация производ- ных, взятых за n-й момент времени, при задании вихря на берегу: при y = 0 ω = yw – zv . (21) На открытой границе ставится условие при y = B ω y = 0, (22) следующее из линейной аппроксимации условия w = 0, когда на границе и wyy = 0. Функция тока ψ на поверхности моря и на дне не меняется по y, по- скольку там равна нулю вертикальная скорость и имеет место равенство ψ y = w. Учитывая, что функция тока определяется с точностью до констан- ты, и задавая ее для определенности равной нулю на дне, условия для функ- ции тока примут вид: при z = H ψ = 0, (23) при z = 0 ψ = Q0, (24) при y = 0 ψ = – ∫ z vdz 0 , (25) при y = B ψ y = 0. (26) ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 68 Уравнения (5), (6), (16) и (1) решаются методом переменных направле- ний (продольно-поперечной прогонки). При этом используется двухслойная схема по времени с неявной аппроксимацией адвективных и диффузионных членов. Пространственная аппроксимация проводится с использованием бокс-метода на прямоугольной сетке типа C по терминологии Аракавы. Ад- вективные члены аппроксимируются схемой первого порядка точности, а именно направленными разностями. Полученные уравнения решаются мето- дом немонотонной прогонки. Уравнение Пуассона (18) аппроксимируется на 5-точечном шаблоне и решается методом верхней релаксации. Компоненты скорости u, v определяются в серединах граней бокса, па- раллельных оси z; компонента скорости w – в серединах граней, параллель- ных оси y; вихрь ω и функция тока ψ – в узлах бокса; температура T и соле- ность S – в центре бокса. Твердые границы области проходят через грани боксов, а именно: границы y = const – через грани, где рассчитываются ком- поненты скорости u, v, а границы z = const – через грани, где определяется компонента w; открытая граница y = B проходит через центры боксов. Как указывалось выше, в сентябре 2011 г. в фоновой стратификации вда- ли от берега наблюдалась довольно редкая ситуация, когда существовал очень резкий слой скачка температуры, а соленость при этом практически не менялась с глубиной. Представляется, что формирование такой стратифика- ции связано с нагоном вод, обусловленным действием ветров восточных румбов, имевших место как во время экспедиции, так, по-видимому, и перед экспедицией. Заметим, что береговая линия в рассматриваемом районе, где находятся все 8 гротов (карстовых полостей), направлена в целом с юго- востока на северо-запад (море – слева, берег – справа). Поэтому юго- восточный ветер, т. е. вдольбереговой ветер, дующий с юго-востока на севе- ро-запад, будет, в соответствии с теорией Экмана, приводить к нагону в по- верхностном слое теплых вод открытого моря к берегу. Для исследования особенностей субмаринной разгрузки в данной синоп- тической ситуации фоновую стратификацию, в соответствии с данными на- блюдений, зададим в следующем виде: Tfon = {22°C при 0 ≤ z ≤ 8 м; 9°C при 8 м < z ≤ H}, Sfon = 17,8‰. (27) Расчеты проводились для области с горизонтальным, для простоты, дном (H = 20 м, B = 200 м) на равномерной сетке с шагами ∆ y = 1 м, ∆ z = 10 см, шаг по времени ∆ t = 0,72 с. Коэффициенты турбулентной вязкости и диффу- зии принимались следующими: A = 10 см2/с, AL = 1000 см2/с, µ = 0,1 см2/с, µ L = 100 см2/с. Температура втекающей воды T0 = 15°C, соленость S0 = 0‰. Скорость втекания воды из источника, залегающего в метровом слое (Z1 = = 10 м, Z2 = 11 м), равна v0 = 0,5 см/с. К сожалению, нет достоверных данных о реальных значениях расходов и глубин залегания источников подземных вод, вследствие чего выбор указан- ных выше параметров не является однозначным и моделирование носит каче- ственный характер. Далее в работе обсуждаются возможные механизмы и моделируются сценарии развития ситуаций при разных соотношениях внеш- них воздействий. Результаты численных экспериментов, если говорить о ко- ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 69 личественном характере, оказались чувствительными к значениям перемен- ных, таких как скорость истекания источника, напряжение ветра, коэффици- енты турбулентной вязкости и диффузии. Выбор коэффициентов турбулент- ной вязкости и диффузии связан также с возможностями используемого ме- тода переменных направлений. Так, для выбранного пространственного раз- решения используемый метод оказывается неустойчивым при увеличении коэффициентов горизонтальной диффузии, и это несмотря на применение схем первого порядка точности при аппроксимации уравнений для темпера- туры и солености. Возможно, что на неустойчивость численной схемы влияет и принятое условие свободного протекания на открытой границе. Подчеркнем, что в двумерной модели, в отличие от трехмерной, понятие о расходе источника – относительно условное, поскольку в этой модели ис- точник представляет собой бесконечную щель. Поэтому задание скорости истечения v0 в двумерной модели не вполне адекватно ее заданию в трехмер- ной модели. Также отметим, что в описываемых ниже численных экспериментах роль процедуры «конвективного приспособления», как выявилось в расчетах, ока- залась несущественной. Как и ранее при исследовании субмаринной разгруз- ки в карстовой полости, когда пресные воды поднимались к поверхности мо- ря [5], главным фактором, влияющим на формирование наблюдавшейся трехмерной структуры полей температуры и солености, является образование вертикальной ячейки циркуляции. Эта ячейка циркуляции обусловлена за- вихренностью поля скорости, возникающей при втекании пресной воды в соленую воду (ω t ~ Ry) – так называемый бароклинный эффект [8]. Прежде чем переходить к описанию численных экспериментов, отметим, что имеет место полная аналогия субмаринной разгрузки с выпуском и рас- пространением в море сточных, фактически пресных, вод. Исследованию трансформации сточных вод у Южного берега Крыма посвящен цикл работ, выполненных в ИнБЮМе НАН Украины под руководством В.И. Заца [9, 10]. Было установлено, что для того чтобы сточные воды не выходили на поверх- ность моря в его прибрежной зоне, их надо выпускать на достаточно боль- шом удалении от берега и на достаточно большой глубине, по крайней мере глубже пикноклина. Пикноклин при этом является препятствием для подъема сточных вод в поверхностный слой моря. Кстати, эти рекомендации были учтены при создании глубоководных выпусков сточных вод в Ялте. Обсуждение результатов. Остановимся теперь на описании результатов нескольких численных экспериментов. Будем считать, что ось X направлена вдоль берега с юго-востока на северо-запад, ось Y – в открытое море, ось Z – вертикально вниз; начало координат находится на поверхности моря в точке пересечения с берегом, т. е. имеем левую систему координат. Поскольку, как уже говорилось, эксперименты носят качественный ха- рактер и имеют цель выявить роль различных факторов в механизме форми- рования промежуточного слоя пониженной солености, естественно было на- чать с простейшего случая. А именно с экспериментов, в которых не учиты- ваются действие ветра (τ x = τ y = 0) и вращение Земли (f = 0), а гидрологиче- ская структура вод формируется только за счет втекания пресных вод из глу- боководных подземных источников − основополагающего фактора в данном явлении. ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 70 0 50 100 150 Y,м 9 11 13 15 17 19 21 4 8 12 16 Z,м 0 50 100 150 Y,м 16 16.4 16.8 17.2 17.6 0 50 100 150 Y,м 11 11.4 11.8 12.2 12.6 13 13.4 T S R 0 50 100 150 Y,м -7 -5 -3 -1 1 3 4 8 12 16 Z,м 0 2 4 6 8Y,м -3 -2 -1 0 1 2 3 0 2 4 6 8Y,м -3 -2.6 -2.2 -1.8 -1.4 -1 -0.6 -0.2 v w Р и с. 1. Распределения температуры T (°C), солености S (‰) и плотности R (у. е.) − вверху, функции тока ψ (м2/с) и компонент скорости течения v, w (см/с) вблизи источника пресной воды − внизу (результаты расчета приведены на момент времени t = 4 ч; ветер и вращение Земли не учитываются; штриховыми кривыми указаны изолинии T = 12°C, S = 17‰, R = = 12 у. е., v = 0, w = 0) На рис. 1 показаны распределения полей температуры T, солености S, плотности R, функции тока ψ и горизонтальной v и вертикальной w компо- ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 71 нент скорости течения, которые сформировались через 4 ч действия источни- ка. В поле температуры хорошо виден термоклин, температура в котором резко меняется в слое от 6 до 7 м. В поле солености отмечается подповерхно- стный слой пониженной солености, который располагается около термоклина в слое (6 – 10 м), причем минимум солености, равный 16,2‰, достигается вблизи источника пресной воды на глубине примерно 7 м. В поле плотности наблюдается устойчиво стратифицированный по вертикали пикноклин, тол- щина которого близка к толщине слоя пониженной солености. В функции тока четко видны две ячейки вертикальной циркуляции, причем в верхнем слое моря происходит вращение против часовой стрелки, а в глубинном слое – по часовой стрелке. В результате в промежуточном слое 6 – 10 м малосоле- ная вода вытекает в открытое море, а из открытого моря в поверхностном и глубинном слоях более соленая вода, наоборот, подтекает к источнику. На рисунке также показаны детализированные распределения компонент скоро- сти течения v, w в 10-метровой окрестности от источника. Четко видны зоны интенсивного подъема вод (до −3 см/с), а также их оттока (до 2,4 см/с) и под- тока (до −2,6 см/с). На рис. 2 приведены распределения тех же полей, но сформировавшихся через 24 ч действия источника. Наблюдаются практически те же особенно- сти, только теперь термоклин, а также верхняя граница слоя пониженной со- лености и пикноклина поднялись ближе к поверхности моря. Теперь термо- клин находится на глубине примерно 4 м. Заметим, что нижние границы слоя пониженной солености и пикноклина остались практически на той же глуби- не, что и раньше, а именно на глубине залегания источника пресной воды. Понизилась при этом также соленость, ее минимальная величина уменьши- лась до 16‰. Таким образом, со временем промежуточный слой пониженной солености расширяется, распространяясь к поверхности, но не выходя за рас- сматриваемое время непосредственно к поверхности моря, и при этом все больше распресняется. Следующие численные эксперименты проводились с учетом вращения Земли, причем параметр Кориолиса принимался постоянным f = 10-4 с-1. Есте- ственно, что из-за малости масштабов влияние силы Кориолиса оказалось не очень существенным, разумеется, при выбранных значениях определяющих параметров. На рис. 3 показаны распределения полей температуры T, солено- сти S, плотности R и компонент скорости течения u, v и w, которые сформи- ровались через 36 ч действия источника при учете силы Кориолиса. При со- поставлении с рис. 1 и 2 видно, что по сравнению с расчетом на 24 ч и при f = 0 мало что изменилось. Разумеется, теперь появилось «поперечное» тече- ние (u ≠ 0), скорость которого не превышает скорости «продольного» тече- ния. Все дальнейшие расчеты проводились в приближении вращающейся плоскости. На следующем этапе была проанализирована эволюция со време- нем формирующейся термохалинной структуры вод. На рис. 4, а приведены вертикальные профили температуры T, солености S и плотности R на рас- стоянии 27 м от берега. Они показаны для различных моментов времени че- ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 72 рез каждые 12 ч, начиная с 12 и до 72 ч. Видно, что слой пониженной солено- сти постепенно расширяется, но так и не выходит на поверхность моря, при- чем мало изменяется после 36 ч. Также мало меняются после 36 ч профили температуры и плотности. 0 50 100 150 Y,м 9 11 13 15 17 19 21 4 8 12 16 Z,м 0 50 100 150 Y,м 15.8 16.2 16.6 17 17.4 17.8 0 50 100 150 Y,м 11 11.4 11.8 12.2 12.6 13 13.4 T S R 0 50 100 150 Y,м -6 -5 -4 -3 -2 -1 0 1 4 8 12 16 Z,м 0 2 4 6 8Y,м -3 -2 -1 0 1 2 0 2 4 6 8Y,м -3 -2.6 -2.2 -1.8 -1.4 -1 -0.6 -0.2 v w Р и с. 2. То же, что на рис. 1, только при t = 24 ч ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 73 0 50 100 150Y,м 9 11 13 15 17 19 21 0 50 100 150Y,м 15.8 16.2 16.6 17 17.4 17.8 0 50 100 150Y,м 11 11.4 11.8 12.2 12.6 13 13.4 T S R 0 50 100 150Y,м -1 -0.6 -0.2 0.2 0.6 1 0 50 100 150 Y,м -3 -2 -1 0 1 2 0 2 4 6 8 Y,м -3 -2.6 -2.2 -1.8 -1.4 -1 -0.6 -0.2 v wu 4 8 12 16 Z,м 4 8 12 16 Z,м Р и с. 3. Распределения температуры T (°C), солености S (‰) и плотности R (у. е.) − вверху, компонент скорости течения u, v, w (см/с) − внизу (вертикальная скорость показана только вблизи источника; результаты расчета приведены при t = 36 ч; ветер не учитывается; вращение Земли учитывается; штриховыми кривыми указаны изолинии T = 12°C, S = 17‰, R = 12 у. е., u = 0, v = 0, w = 0) ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 74 9 11 13 15 17 19 21 20 18 16 14 12 10 8 6 4 2 0 Z,м 16 16.6 17.2 17.8 11 11.6 12.2 12.8 13.4 S RT 12 12 12 12 72 72 72 9 11 13 15 17 19 21 20 18 16 14 12 10 8 6 4 2 0 Z,м 15.2 15.8 16.4 17 17.6 11 11.6 12.2 12.8 13.4 12 12 72 72 72 τx=τy=0 τx=0,5см2/с2 τy=0 T S R а б Р и с. 4. Вертикальные профили температуры T (°C), солености S (‰) и плотности R (у. е.) через каждые 12 ч, начиная с момента времени t = 12 ч и до t = 72 ч: а – ветер не учитывается, б – учитывается действие вдольберегового ветра ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 75 Как следует из проведенных расчетов, слой пониженной солености − ос- новной индикатор субмаринной разгрузки − со временем расширяется и под- нимается к поверхности, что не вполне соответствует данным наблюдений, согласно которым он как бы находится на одной глубине. Конечно, эти дан- ные получены в экспедиции, длившейся чуть более суток, при этом неизвест- ны внешние и внутренние факторы, обусловившие и поддерживающие суще- ствование такого слоя, поэтому следующие численные эксперименты были проведены, чтобы выявить возможные факторы и механизмы, способствую- щие его возникновению. Прежде всего, если вспомнить, что имеющая место во время экспедиции фоновая стратификация была сформирована, как нам представляется, в результате действия нагонного ветра, то имеет смысл рас- смотреть влияние действия такого ветра на результат субмаринной разгрузки. На рис. 4, б, как и на рис. 4, а, приведены вертикальные профили температу- ры T, солености S и плотности R для моментов времени с 12 и до 72 ч в слу- чае, когда на поверхности моря действует вдольбереговой ветер со скоростью около 5 м/с, т. е. с напряжением τ x = 0,5 см2/с2, τ y = 0. Хорошо видно, что учет ветра приводит к сжатию слоя пониженной солености, причем с ощути- мым понижением солености (почти как по данным измерений). Также как бы фиксируются на определенной глубине термоклин и пикноклин. Отметим, что при учете ветра кривые температуры T, солености S и плотности R со временем опускаются, а не поднимаются к поверхности моря, как это имеет место при отсутствии ветра. На рис. 5 приведены распределения полей температуры T, солености S, плотности R и компонент скорости течения u, v и w, которые сформировались через 36 ч действия источника и вдольберегового ветра. Видно, что по срав- нению с расчетом, когда действие ветра не учитывалось, теперь в самом деле получаем резкий термоклин практически на той же глубине, как и в открытом море. Хорошо видны тонкий слой сравнительно низкой солености и довольно резкий пикноклин. Отметим, что теперь формируется относительно сильное направленное вдоль берега течение (u > 0), очевидно вызванное ветром. Под- черкнем, что в поверхностном слое усиливаются нагонное течение (v < 0) и опускание вод (w > 0), вызывающие сжатие слоя пониженной солености. Если на рис. 3 и на рис. 5 вертикальная скорость была показана только вблизи берега, где она достигает максимальных величин, то на рис. 6 пред- ставлены распределения вертикальной скорости во всей области, причем как при отсутствии ветра, так и при его учете. Видно, что влияние ветра проявля- ется в основном в усилении подъема вод и расширении зоны подъема, а так- же в перестройке вертикальных движений в поверхностном слое моря. В качестве отступления отметим, что в численных экспериментах с более сильным ветром при фиксированном потоке воды из пресноводного источни- ка происходит вообще разрушение термоклина, пикноклина и слоя понижен- ной солености, образуются пятна распресненной воды, частично всплываю- щие на поверхность моря. Поскольку в наблюдениях такие явления не фик- сировались, анализ соответствующих расчетов, находящихся к тому же на ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 76 грани гидродинамической неустойчивости и неустойчивости используемой численной схемы, здесь проводиться не будет, это требует дополнительных исследований. 0 50 100 150Y,м 9 11 13 15 17 19 21 0 50 100 150Y,м 15.4 15.8 16.2 16.6 17 17.4 17.8 0 50 100 150Y,м 11 11.4 11.8 12.2 12.6 13 13.4 T S R 0 50 100 150Y,м 0 4 8 12 16 20 0 50 100 150 Y,м -2.5 -1.5 -0.5 0.5 1.5 0 2 4 6 8 Y,м -2.2 -1.8 -1.4 -1 -0.6 -0.2 0.2 v wu 4 8 12 16 Z,м 4 8 12 16 Z,м Р и с. 5. То же, что на рис. 3, только при учете действия вдольберегового ветра ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 77 5 10 15 Z,м 0 50 100 150 Y,м 5 10 15 Z,м (-2,6,-0,1)(-0,1,-0,01) (-0,01,0) (0) (0,0,01)(0,01,0,07) (-2,2,-0,1)(-0,1,-0,01) (-0,01,0)(0)(0,0,01)(0,01,0,25) Р и с. 6. Распределения вертикальной скорости течения w (см/с) на момент времени t = 36 ч без учета ветра (вверху) и при учете вдольберегового ветра (внизу) (в скобках соответствую- щего цвета указаны пределы изменения вертикальной скорости; штриховыми кривыми указа- ны изолинии w = 0) В заключение остановимся еще на двух численных экспериментах, ре- зультаты которых показаны на рис. 7. В этих расчетах не учитывалось дейст- вие ветра, а исследовалось влияние заглубления источника и усиления его мощности. ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 78 9 11 13 15 17 19 21 20 18 16 14 12 10 8 6 4 2 0 Z,м 17.1 17.3 17.5 17.7 11 11.6 12.2 12.8 13.4 6 6 6 36 36 36 T S R S RT 2 9 11 13 15 17 19 21 20 18 16 14 12 10 8 6 4 2 0 Z,м 14.8 15.6 16.4 17.2 11 11.6 12.2 12.8 13.4 2 2 12 12 12 а б Р и с. 7. Вертикальные профили температуры T (°C), солености S (‰) и плотности R (у. е.), полученные без учета ветра: а – через каждые 6 ч, начиная с момента времени t = 6 ч и до мо- мента t = 36 ч, при более глубоком положении источника пресной воды; б – через каждые 2 ч, начиная с момента времени t = 2 ч и до момента t = 12 ч, при более интенсивном источнике ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 79 На рис. 7, а показано, что заглубление источника пресной воды, который теперь располагается в слое 16 – 17 м, в основном сказывается на солености, увеличивая ее в промежуточном слое (до 17,1‰) и приводя к существенному растяжению этого слоя со временем, поскольку нижняя граница слоя остает- ся как бы привязанной к источнику. На рис. 7, б видно, что усиление мощности источника, когда v0 = 2 см/с, приводит к тому, что пресная вода, практически не перемешиваясь, быстро поднимается к поверхности моря. Буквально через 3 ч температура воды на поверхности моря достигает примерно 10°C, а соленость 15‰, и далее эти значения практически не меняются. Дополнительные эксперименты показа- ли, что действие более сильного (τ x = 1 см2/с2) вдольберегового ветра не мо- жет остановить выхода распресненных вод на поверхность моря. Заметим, что при отсутствии ветра и значении скорости v0 = 1 см/с рас- пресненная вода за те же 36 ч к поверхности моря не поднимается, так же как в варианте, когда v0 = 0,5 см/с. Это указывает на зависимость результатов расчета от интенсивности источника. Более того, можно сделать вывод о том, что наблюдаемая на полигоне термохалинная структура вод, а именно про- межуточный слой пониженной солености и резкий термоклин [7], может быть обусловлена сравнительно слабой субмаринной разгрузкой (c потоком V0 < 0,02 м2/с) и одновременным действием нагонного не очень сильного (со скоростью ≈5 м/с) ветра. Выводы. На основании проведенных численных экспериментов можно согласиться с частным мнением А.И. Чепыженко о том, что формирование наблюдавшейся трехмерной структуры полей температуры и солености в районе м. Айя в сентябре 2011 г. может быть связано с наличием источников пресной воды, расположенных вне гротов на сравнительно большой глубине. В данном случае расположенных глубже 9 м – глубины слоя скачка темпера- туры в рассматриваемый период времени. Можно предположить, что наблюдаемые на полигоне промежуточный слой пониженной солености и резкий термоклин обусловлены одновремен- ным действием субмаринной разгрузки и нагонного ветра. СПИСОК ЛИТЕРАТУРЫ 1. Иванов В.А., Прусов А.В., Рябцев Ю.Н., Шапиро Н.Б. Физические механизмы смешения морских вод с водами субмаринной разгрузки // Современные проблемы океанологии. − Севастополь: МГИ НАН Украины, 2009. − 90 с. 2. Иванов В.А., Прусов А.В., Чепыженко А.И. Субмаринная разгрузка подземных вод на м. Айя (Крым) // Доповiдi НАН Украïни. – 2010. – № 2. – С. 103 – 109. 3. Шапиро Н.Б. Моделирование трехмерной структуры гидрологических полей в карсто- вой полости под действием субмаринной разгрузки подземных вод // Морской гидро- физический журнал. – 2010. – № 1. – С. 46 – 62. 4. Шапиро Н.Б. Численная модель субмаринной разгрузки пресных вод в карстовой по- лости // Там же. – 2011. – № 5. – С. 66 – 84. 5. Михайлова Э.Н., Шапиро Н.Б. Роль гидростатического приближения при моделирова- нии субмаринной разгрузки пресных вод в карстовой полости // Там же. – 2012. – № 3. – С. 44 – 60. ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 80 6. Шапиро Н.Б. Моделирование трехмерной структуры гидрологических полей в карсто- вой полости на мысе Айя под действием субмаринной разгрузки // Экологическая безо- пасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. – Севастополь: МГИ НАН Украины, 2011. – Вып. 25. Т. 2. – С. 252 – 274. 7. Экспедиционные исследования МГИ НАН Украины «Субмаринная разгрузка пресных вод в районе м. Айя» в сентябре 2011 г. // Отчет. − Научные фонды МГИ НАН Украи- ны. − Севастополь, 2011. − 46 с. 8. Кочин Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика. Т. 1. – М.: ГИТТЛ, 1955. – 560 с. 9. Зац В.И., Немировский М.С., Андрющенко Б.Ф. и др. Опыт теоретического и экспери- ментального исследования проблемы глубоководного сброса сточных вод на примере района Ялты. – Киев: Наукова думка, 1973. – 274 с. 10. Гольдберг Г.А., Зац В.И., Ациховская Ж.М. и др. Моделирование процессов самоочище- ния вод шельфовой зоны. – Л.: Гидрометеоиздат, 1991. – 230 с. Морской гидрофизический институт НАН Украины, Материал поступил Севастополь в редакцию 03.07.12 E-mail: men_sh@mail.ru После доработки 04.09.12 АНОТАЦІЯ За результатами експедиції виявлено, що особливістю субмаринного розванта- ження у вересні 2011 р. було існування чітко вираженого проміжного шару зниженої солоності біля шару різкого стрибка температури, які спостерігалися у відкритій частині моря поблизу і на видаленні (до 60 м) від берега, причому на достатньо великому протязі (декілька сотень метрів) паралельно береговій лінії. На підставі проведених чисельних експериментів встановлено, що головним фактором, який впливає на формування такої структури полів температури та солоності, є наявність дже- рел прісної води, розташованих глибше за шар стрибка температури в даний період часу. По- казано, що порівняно тонкий шар зниженої солоності та різкий термоклин можуть бути обумо- влені одночасною дією субмаринного розвантаження та нагонного вітру. Ключові слова: субмаринне розвантаження, негідростатична двовимірна модель, чисель- ний експеримент. ABSTRACT Results of the expedition show that the feature of submarine discharge in September, 2011 was existence of the well-pronounced intermediate layer with reduced salinity nearby the layer of sharp temperature drop. Both of them were observed in the open sea close to the coast and up to 60 km off it for rather a long space (a few hundreds meters) along the coastline. Based on the performed numerical experiments it is found that the basic factor influencing for- mation of such a structure of temperature and salinity fields consists in presence of freshwater sources located below the temperature drop layer during the period under consideration. It is shown that rather a thin layer with reduced salinity and sharp thermocline can be conditioned by simultaneous action of submarine discharge and surge wind. Keywords: submarine discharge, non-hydrostatic two-dimensional model, numerical experiment. ISSN 0233-7584. Мор. гидрофиз. журн., 2013, № 5 81 mailto:men_sh@mail.ru

Моделирование особенностей субмаринной разгрузки пресных вод в районе м. Айя в период нагона

Institution

Ähnliche Einträge