Оценка коэффициента вертикальной турбулентной диффузии по данным ctd/ladcp-измерений в северо-западной части Черного моря в мае 2004 года

УДК 551.46

А.Н. Морозов, Е.М. Лемешко

Оценка коэффициента вертикальной турбулентной диффузии по данным CTD/LADCP-измерений в северо-западной части Черного моря в мае 2004 года

Приводятся осредненные по ансамблю станций профили частоты плавучести, вертикальных сдвигов измеренных и геострофических скоростей течений, рассчитанных по данным CTD/LADCP-измерений, выполненных в северо-западной части Черного моря в мае 2004 г. Оценки параметров вертикального перемешивания получены с использованием модели HWF86. Детально обсуждается методика расчетов. Получено, что скорость диссипации турбулентной энергии изменяется в пределах 10-10 - 210-9 Вт/кг в слое 70 - 500 м и имеет хорошо выраженный максимум в слое основного пикноклина. Коэффициент вертикальной турбулентной диффузии имеет значение 10-6 м2/с в центре пикноклина и увеличивается с глубиной по экспоненциальному закону, достигая 10-5 м2/с на глубине ~500 м.

Ключевые слова: вертикальное перемешивание, LADCP, Черное море.

Вертикальное турбулентное перемешивание - одна из актуальных тем современной океанографии. Знание коэффициентов вертикального перемешивания, их пространственной и временной изменчивости необходимо для моделирования как крупномасштабной циркуляции, так и локальных динамических процессов. В настоящее время достигнут значительный прогресс в установлении взаимосвязи параметров турбулентности с мелкомасштабными характеристиками гидрофизических полей [1 - 3]. Большинство полученных зависимостей для оценки коэффициента вертикальной турбулентной диффузии по данным мелкомасштабных измерений основано на соотношении вертикальных сдвигов течений и частоты плавучести. В Черном море характеристики вертикальных сдвигов изучены мало. Цель настоящей работы - частично восполнить этот пробел на основе данных измерений профилей скорости течений, полученных с использованием погружаемого акустического до-плеровского профилемера течений (LADCP), и представить оценки параметров вертикального перемешивания с использованием данных мелкомасштабных CTD/LADCP-измерений.

В работе использованы данные, полученные в рейсе НИС Akademik Болгарской АН 9 - 18 мая 2004 г. в северо-западной части Черного моря [4]. CTD-измерения выполнялись с помощью зонда SeaBird SBE9+. Профили температуры, солености, плотности интерполировались на сетку 1 м по глубине. Профили течений получены с помощью LADCP на основе WHM300 производства RDI, США. Измерения течений проводились при установленном размере сегмента глубины ( B ) для первой группы станций 10 м (16 ст.), для второй - 4 м (12 ст.). Обработка LADCP-данных выполнена в соответст-

© А.Н. Морозов, Е.М. Лемешко, 2014

58

ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 1

вии с [5]. На рис. 1 приведено расположение станций (квадратики - первая группа, треугольники - вторая), серыми линиями схематично показана циркуляция вод в районе, сплошными линиями обозначены Основное Черноморское течение и антициклонический вихрь, пунктирные линии ограничивают области сонаправленного движения вод [4].

Р и с. 1. Схема циркуляции вод и положение станций: квадратики - первая группа, треугольники - вторая

Профили вертикальных сдвигов, полученные в результате LADCP-изме-рений, содержат как составляющие стабильных течений, так и составляющие, вызванные внутренними волнами и другими динамическими процессами, имеющими случайные фазы на разных глубинах в разное время. Измеренный профиль сдвига оказывается значительно «зашумлен» по сравнению с профилем, рассчитанным по геострофическим соотношениям, поэтому для получения осредненного профиля измеренных LADCP сдвигов использовались данные не менее 10 станций [6]. Осреднение выполнялось относительно глубин залегания центра основного пикноклина на каждой станции с последующим смещением полученного профиля на среднее по ансамблю значение глубины [7]. В качестве параметра осреднения выбран квадрат вертикального

сдвига течений (БИ2ЫГ)СР ) :

ХИ^^Р (2) = иК2) + У?(2),

0233-7584. Мор. гидрофиз. журн., 2014, № 1

59

где иг ,Уг - производные по глубине широтной и меридиональной компонент скорости течения соответственно. Аналогичное осреднение было выполнено для сдвигов, рассчитанных по геострофическим соотношениям (8к'2Е05 ) :

( V £

Ж

лЛ , +

Ах

2 (Ару ^

Ау

/

где ^ - ускорение свободного падения (9,82 м/с2); / - параметр Кориолиса (10-4 с-1); р - плотность; Арху - приращение плотности в направлении х,у на расстоянии Ах,Ау между соседними станциями ~17 км. Аналогичное осреднение сделано и для квадрата частоты плавучести (N2):

N2 = £р, Р

где рг - производная плотности по глубине. На рис. 2 приведены результирующие осредненные профили для двух групп станций. Результирующие профили для первой группы станций показаны на рис. 2, а: ^к1Лоср) - черная сплошная линия; () - штриховая линия; {N2) - серая линия ((.......) - оператор осреднения по ансамблю станций). На рис. 2, б приведены

аналогичные профили для второй группы станций. Сопоставление профилей средних квадратов геострофических и измеренных сдвигов показывает их относительную соизмеримость в слое основного пикноклина (100 - 150 м),

отношение )/{) составляет ~ 0,25. Ниже основного пикнокли-

на ^И2Е08 ) уменьшается с глубиной значительно быстрее, чем (5к^ЛГ)ср ), и

на глубине 400 м )/<) составляет менее 0,01. Наблюдаемое

различие в изменении измеренных и геострофических сдвигов с глубиной свидетельствует о том, что сдвиги в слоях моря ниже основного пикноклина определяются агеострофическими процессами. Приведенные на рис. 2, а, б профили проявляют аналогичную зависимость от глубины. Некоторые различия могут быть вызваны разными районами проведения измерений (см. рис. 1). Дополнительно различие значений (^к21ЛГ^ср ) для двух групп станций определяется присущим ЛОСр пространственным осреднением, которое упрощенно представляется передаточной функцией (Н ) следующего вида [8]:

н (к) = \*т{лВк )Лб

рВк

где к - вертикальное волновое число.

2

60

ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 1

100 -

200 -

od Я

S 300 ю гу ч

l-Ч

400 -

500 -

B=10 м

(Sh2GEOs).-"

<Sh2L

I—i 11 iniij—i 11 iniij—i 11111n|

10-8 10-7 ю-6

<Sh2G

s), <Sh2,

10-6 10-5 10-4 c-2

.), < N2), с-2

100-

S 200 —

е£ - J

Я * *

S 300 — \

ю

у _

ч г' ъ

1-4 400 — г i

в

500 -

а  111

B =4 м

<sh2GEOs)

<Sh2L

б

I-1 I I lllllj-1 I I lllllj-1 I I lllllj-1 I I lllllj-ГТП

10-8 10-7 10-6 10-5 10-4 с-2

p), < N2), с-2

<Sh2GEOS), < Sh2LADCP/

Р и с. 2. Профили квадратов частоты плавучести, вертикальных измеренных и геострофических сдвигов: а - для первой группы станций

В основу большинства моделей, связывающих мелкомасштабные характеристики гидрофизических полей с параметрами вертикального перемешивания, положено допущение о равенстве скорости передачи кинетической энергии внутренних волн по спектру в область мелких масштабов и скорости диссипации кинетической турбулентной энергии (e). В соответствии с [9] это допущение выполняется для многих районов Мирового океана. В настоящее время для оценки e по данным мелкомасштабных измерений наиболее распространено использование моделей MM81 [10] и HWF86 [11]. В рамках данной работы для расчета e применяется модель HWF86, которая с учетом поправок на отклонение спектра внутренних волн от канонического вида GM76 [12] приводит к следующему соотношению [3, 13]:

N2 {Sh2)2 . f cosh-1

e = e0—т—-'—г F (Rm)A--—J—,

0 No {ShlM76)2 f0cosh

гдее0 = 7,8 •10-10 Вт/кг, f0 = 7,29 • 10-5 с-1, N0 = 5,24 -10-3 с-1 - постоянные, соответствующие спектру внутренних волн GM76; (ShGM76^ - среднее значение квадрата сдвига, полученное интегрированием спектра GM76; F (Rw) -функция частотного наполнения [3]. Расчет коэффициента вертикальной диффузии Kp производится по соотношению Осборна [14]:

e

Kr= Г

N2

ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 1

61

где Г - коэффициент эффективности перемешивания, который принято полагать равным 0,2 [15].

Функция частотного наполнения в соответствии с [16] может быть представлена следующим соотношением:

F (Rw) =

1 +1/R

4/3

(^ V Rw — 1

1/2

где Яа - отношение кинетической (ЕК) и потенциальной (Ер) энергий во внутренней волне. В Черном море частотный спектр внутренних волн имеет ярко выраженный максимум вблизи локальной инерционной частоты [17]. На

(2 3 ^

горизонтах 35 - 450 м в окрестности / I — / ——/ \ концентрируется не менее 75% энергии внутренних волн (// — N), некоторое расширение диапазона в сторону частот, меньших / , вызвано наличием фоновой завихренности

3

поля течений [18]. На частотах, больших —/, наклон спектра составляет

(0~2. Для оценки характерного для Черного моря значения Ят использовано следующее соотношение [16]:

„ ЕК Sh2 Л, =

EP N h2 '

где h - деформация. По ансамблю станций получены средние спектры числа Фруда и деформации для слоя глубин 250 - 450 м, где N меняется плавно. Расчет деформации выполнялся аналогично [19] по соотношению

N2 - N2

hz = — , N2

где N2 - «равновесный» профиль квадрата частоты плавучести, представлялся полиномиальной зависимостью для каждой станции. Соответствующие спектры приведены на рис. 3, а: черная линия - спектр числа Фруда, откорректированный с учетом передаточной функции (ФРг Сог); штриховая - спектр числа Фруда без коррекции (Фрг); серая - спектр деформации (Фэйшп). Среднее по спектру значение Rw составляет ~12, соответственно в расчетах F(Rw) полагалась равной 0,34. На рис. 3, б представлены скорректированные спектры числа Фруда для право- (ФРг CW, черная линия) и лево- (ФРг CCW, серая линия) вращательных компонент. В области длин волн Äv > 50 м доминирует вращение по часовой стрелке, отношение компонент Фргот/ФргСот = = 2,5, т. е. большая часть волновой энергии, ~70%, распространяется сверху 62 ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 1

вниз. Это может быть вызвано распространением внутренних волн с частотой, близкой к инерционной, которые являются одним из существенных источников энергии турбулентности в толще вод моря [18, 20].

Р и с. 3. Спектры числа Фруда и деформации

Средние квадраты сдвигов, вызванных внутренними волнами, по данным LADCP определяются, как правило, через расчет спектральной плотности в слоях ~300 м, на которые разбивается профиль [13]. Имеющиеся данные не позволяют использовать такую методику, так как глубина зондирований ограничивалась 500 м. Тем не менее можно получить оценку сверху значений соответствующего квадрата сдвига по соотношению

<.Sh2 > = {SHLADCP > - <ShGEOS > - <ShNOISE > ,

где <ShNOISE > - шумовая составляющая измерений, определенная по параметру ошибки, который регистрируется ADCP. Учет пространственного осреднения производится при расчете <ShGM76 > :

ks

<ShGM 76 > = { Фсмм7б(к, <N2 >)H (k)dk,

0

где ks - вертикальное волновое число среза, полагаемое равным 0,1 м-1; Ф<зм 76 - соответствующая GM76 спектральная плотность [21].

ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 1 63

На рис. 4 приведены результирующие средние профили е (рис. 4, а) и Кр (рис. 4, б), черная линия соответствует ансамблю станций с В = 4 м, серая - с В = 10 м. Поведение е с глубиной выявляет наличие значительного максимума в слое основного пикноклина. В целом по глубине значения е изменяются в пределах 2 10-10 - 210-9 Вт/кг, что характерно для многих районов Мирового океана [3, 13]. В слое основного пикноклина Кр имеет значение 10-6 м2/с. Изменение коэффициента диффузии с глубиной (рис. 4, б) выявляет близкую к экспоненциальной зависимость Кр » 4 • 10-7 ехр(7 • 10-3 г), и на глубине ~500 м его значение достигает 10-5 м2/с.

100 —

200-

св К

К 300 ю

Ч 1-ч

400-

500 -

100 -

200 -

В = 10 м В = 4 м

03 К

К 300

ю у

4

1-4

400 -

500 -

а

В = 10 м В = 4 м

10-1

л—I—I—I—I I 11—

10-9

е, Вт/кг

I—I I 11111—

10-6

—I—I—I I 11111—

10-5

Кр , м2/с

б т-1

Р и с. 4. Профили скорости диссипации турбулентной кинетической энергии - а и коэффициента вертикальной турбулентной диффузии - б

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

- 1989 г. [22] - Кр = 3,1 - 4,4-10-4 м2/с на глубине 500 м;

- 1996 г. [23] - Кр = 3 - 5• 10-5 м2/с в постоянном галоклине,

Кр = 3 - 5-10-4 м2/с на глубине 300 м;

- 2003 г. [24] - Кр = 105 м2/с по всей толще моря для ММ81 и ИЦ^86;

- 2006 г. [25] - Кр = 2 - 8-10-5 м2/с в основном пикноклине;

- 2007 г. [26] - Кр = 4-10-6 - 5• 10-5 м2/с в верхнем 180-метровом слое;

- 2008 г. [7] - Кр = 4-10-6 - 2-10-5 м2/с в слое 100 - 450 м.

64

0233-7584. Мор. гидрофиз. журн., 2014, № 1

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

В настоящее время принято считать оценки, полученные на основе микромасштабных данных [25, 26], наиболее объективными. Метод, обсуждаемый в данной статье, является одним из возможных подходов к задаче и может быть полезен для оценки параметров перемешивания на больших глубинах, где проведение микромасштабных измерений затруднено. Полученные оценки параметров вертикального перемешивания носят предварительный характер в силу допущений, которые были сделаны при выводе формул, и требуют уточнения по данным других измерений. Тем не менее, учитывая широкое использование LADCP в комплексе с CTD, изложенная в работе методика оценки коэффициентов вертикальной диффузии позволяет получать массовый материал о пространственной структуре этого параметра [6, 13, 16, 27].

Работа выполнена в рамках национального проекта «Фундаментальная океанология» и при частичной поддержке международного проекта РЕЯБЕиБ-№7, контракт № 287600.

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

1. Gregg M.C. Scaling turbulent dissipation in the thermocline // J. Geophys. Res. - 1989. -94. - P. 9686 - 9698.

2. Wijesekera H., Padman L., Dillon T. et al. The application of internal-wave models to a region of strong mixing // J. Phys. Oceanogr. - 1993. - 23. - P. 269 - 286.

3. Polzin K.L., Toole J.M., Smith R.W. Finescale parameterizations of turbulent dissipation // Ibid. - 1995. - 25. - P. 306 - 328.

4. Лемешко Е.М., Морозов А.Н., Станичный С.В. и др. Вертикальная структура поля скорости течений в северо-западной части Черного моря по данным LADCP в мае 2004 г. // Морской гидрофизический журнал. - 2008. - № 6. - С. 25 - 37.

5. Морозов А.Н., Лемешко Е.М. Методические аспекты использования акустического доп-леровского измерителя течений (ADCP) в условиях Черного моря // Там же. - 2006. -№ 4. - C. 31 - 48.

6. Polzin K., Kunze E., Hummon J. et al. The finescale response of lowered ADCP velocity profiles // J. Atmosph. Oceanic Techn. - 2002. - 19. - P. 205 - 224.

7. Морозов А.Н., Лемешко Е.М. Вертикальное перемешивание в Черном море по данным CTD/LADCP-наблюдений // Системы контроля окружающей среды. - Севастополь: МГИ НАН Украины, 2008. - С. 266 - 268.

8. Alford M.H., Gregg M.C. Near-inertial mixing: Modulation of shear, strain and microstructure at low latituide // J. Geophys. Res. - 2001. - 106, № C8. - P. 16947 - 16968.

9. GreggM.C., Sanford T.B., Winkel D.P. Reduced mixing from the breaking of internal waves in equatorial waters // Nature. - 2003. - 402. - P. 513 - 515.

10. McComas C.H., Muller P. The dynamic balance of internal waves // J. Phys. Oceanogr. -1981. - 11. - P. 970 - 986.

ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 1

65

11. Henyey F.S., Wright J., Flatte S.M. Energy and action flow through the internal wave field: an eikonal approach // J. Geophys. Res. - 1986. - 91. - P. 8487 - 8495.

12. Cairns J.L., Williams G.O. Internal Waves Observations From a Midwater Float, 2 // Ibid. -1976. - 81, № 12. - P. 1943 - 1950.

13. Naveira Garabato A.C., Oliver K.I.C., Watson A.J. et al. Turbulent diapycnal mixing in the Nordic seas // Ibid. - 2004. - 109, C12010. - 9 p.

14. Osborn T.R. Estimates of the local rate of vertical diffusion from dissipation measurements // J. Phys. Oceanogr. - 1980. - 10. - P. 83 - 89.

15. Moum J.N. Efficiency of mixing in the main thermocline // J. Geophys. Res. - 1996. - 101. -P. 12057 - 12069.

16. FerI. Scaling turbulent dissipation in Arctic fjord // Deep-Sea Res. II. - 2006. - 53. - P. 77 - 95.

17. Морозов А.Н. Спектральные характеристики инерционных колебаний в Черном море // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: МГИ НАН Украины, 2001. - Вып. 2. - C. 61 - 69.

18. Kunze E. Near-Inertial Wave Propogation in Geostrophic Shear // J. Phys. Oceanogr. -1985. - 15. - P. 544 - 565.

19. Mauritzen C., Polzin K.L., McCarthey M.S. et al. Evidence in hydrography and density fine structure for enhanced vertical mixing over the Mid-Atlantic Ridge in the western Atlantic // J. Geophys. Res. - 2002. - 107, № C10, 3147. - 19 p.

20. Иванов В.А., Белокопытов В.Н. Океанография Черного моря. - Севастополь: НПЦ «ЭКОСИ-Гидрофизика», 2011. - 209 с.

21. GreggM.C., Kunze E. Shear and strain in Santa Monica Basin // J. Geophys. Res. - 1991. -96. - P. 16709 - 16719.

22. Богуславский С.Г., Иващенко И.К. Вертикальная мезоструктура глубинных вод Черного моря // Морской гидрофизический журнал. - 1989. - № 5. - C. 25 - 32.

23. Еремеев В.Н., Кушнир В.М. Слоистая структура течений и вертикальный обмен в Черном море // Океанология. - 1996. - 36, № 1. - C. 13 - 19.

24. Самодуров А.С., Иванов Л.И. Среднее вертикальное распределение скорости диссипации турбулентной энергии в Черном море. Сравнение с существующими моделями // Морской гидрофизический журнал. - 2003. - № 3. - C. 3 - 8.

25. Самодуров А.С., Чухарев А.М. Оценка интенсивности вертикального турбулентного обмена в Черном море по экспериментальным данным // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: МГИ НАН Украины, 2006. - C. 524 - 529.

26. Зацепин А.Г., Голенко Н.Н., Корж А.О. и др. Влияние динамики течений на гидрофизическую структуру вод и вертикальный обмен в деятельном слое Черного моря // Океанология. - 2007. - 47, № 3. - С. 327 - 339.

27. Cisewski B., Strass V.H. and Prandke H. Upper-ocean vertical mixing in the Antarctic Polar Front Zone // Deep-Sea Res. - 2005. - 52, Issues 9-10. - P. 1087 - 1108.

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

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

E-mail: anmorozov@yahoo.com После доработки 23.10.12

66

ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 1

АНОТАЦ1Я Наводиться усереднеш за ансамблем станцш профiлi частоти плавучостi, вертика-льних зсувiв вимiряних i геострофiчних швидкостей течш, розрахованих за даними CTD/LADCP-вимiрювань, виконаних в пiвнiчно-захiднiй частинi Чорного моря в травш 2004 р. Оцшки параметрiв вертикального перемiшування отримаш з використанням моделi ЕШЕ86. Детально обговорюеться методика розрахункiв. Отримано, що швидкiсть дисипацп турбулентно! енергп змiнюеться в межах 10-10 - 210-9 Вт/кг в шарi 70 - 500 м i мае добре виражений максимум в шарi основного ткноклину. Коефiцiент вертикально! турбулентно! дифузй мае значення 10-6 м2/с в центрi пiкноклину i збшьшуеться з глибиною за експоненцiальним законом, досягаючи 10-5 м2/с на глибиш ~500 м.

Ключовi слова: вертикальне перемшування, LADCP, Чорне море.

ABSTRACT Averaged over a set of stations, the profiles of buoyancy frequency and vertical shears of the currents' measured and geostrophic velocities calculated by the data of CTD/LADCP-measurements carried out in the northwestern Black Sea in May, 2004 are represented. The estimates of vertical mixing parameters are derived using the HWF86 model. The method of calculations is discussed in details. It is found that the turbulent energy dissipation rate varies within the interval 10-10 - 210-9 W/kg in the 70 -500m layer and has a well-pronounced maximum in the layer of main pycnocline. The coefficient of vertical turbulent diffusion is equal to 10-6 m2/s in the pycnocline center and increases with depth according to the exponential law achieving 10-5 m2/s on the depth ~500 m.

Keywords: vertical mixing, LADCP, Black Sea.

ISSN 0233-7584. Мор. гидрофиз. журн, 2014, № 1

67