Файл: Тареев, Б. А. Динамика бароклинных возмущений в океане.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 130

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

Поскольку градиентно-вихревые волны являются простым следствием принципа сохранения потенциального вихря, то возвращающая сила в уравнениях волнового движения в линейном приближении пропорциональна градиенту завих­ ренности невозмущенного состояния, т. е. в простых случаях, 'рассмотренных выше (см. формулы 2.1.9—2.1.11), пропор­

циональна ——= (3. Если изменения глубины слоя сущест- dy

венны, то уравнение (2.1.1) должно быть заменено на

(2.1.34)

и если Н = Н(у), то в формулах для зональных градиентно­

вихревых волн р следует заменить на

р —

Н

—. М. Лои-

ге-Хиггинс (Longuet-Higgins, 1964,

 

dy

1965а)

дал простое

наглядное объяснение западного направления фазовой ско­ рости градиентно-вихревых волн относительно среды для слоя постоянной глубины на p-плоскости или на вращающей­ ся сферической Земле. Это правило легко может быть обоб­ щено: фазовая скорость градиентно-вихревых волн направ­ лена под прямым углом влево от направления градиента потенциального вихря в невозмущенном состоянии (при обычном определении знака вектора вихря). Таким образом, если глубина океана уменьшается к востоку, то из уравнения (2.1.34) легко найти, что должна существовать компонента фазовой скорости волн, направленная к северу. В случае непрерывно стратифицированной жидкости определение по­ тенциального вихря усложняется (см., например, А. М. Обу­ хов, 1962), но предыдущее правило остается в силе. В част­ ности, из формул (2.1.9), (2.1.10) следует, что в слое постоян­ ной глубины (3-волны, имеющие чисто меридиональное направление, невозможны, а величина зональной фазовой скорости не зависит от направления распространения волны и определяется абсолютной величиной волнового числа k (или длиной волны 2 лIk).

Используя формулу для групповой'зональной скорости двумерных бездивергентных волн Россби:

вытекающую из (2.1.9), Дж. Педлоски (Pedlosky, 1965а)

показал, что западная интенсификация течений может быть, по крайней мере, качественно объяснена с точки зрения дина­ мики волн Россби, если предположить, что возмущения ма­


лых зональных пространственных масштабов {kx~>ky), отра­ жающиеся от западной стенки, диссипируются под влиянием трения вблизи стенки.

Все вышеприведенные соотношения имекй1 существенно дисперсионный характер, т. е. для полной определенности формул следует задать волновое число или частоту. Из фор­ мулы для частоты бароклинных волн Россби (2.1.11) в дву­ слойной модели легко получить, что максимальное абсолют­ ное значение частоты соответствует

 

К = 0,

kx =

//

 

о д

 

=

L

 

 

 

 

 

 

 

Hi -р #2

 

 

 

 

 

I ®max I =

------------- ~f---------- - =

Р/2&о-

(2.1.36)

 

 

2 [§ ' ^

Я 2/( Я 1 +

Я 2)]1/2

 

 

 

 

Полагая полную глубину океана Hi + H2 = H=A км,

тол­

щину главного термоклина # i = 600 м, g = 10 м/сек2

 

 

Pl~ P2

= 2.10—3, так что

л /

g ' ——^

— = 3 — 4'м/сек,

f =

р

 

 

у

&

H l +

H 2

 

 

 

1

= 10"4

сек-1, |3= 10-13 сек-1 см

*,

найдем,

что минимальный

период колебании

Tm in =

 

равен примерно году,

соот-

-----

 

 

 

 

^шах

сот ах/&о~2

см/сек,

а длина

ветствующая фазовая скорость

волны La= 2jt/&oта250 км.

 

 

вычисления

при Н = 4 км

Наряду с этим

аналогичные

для баротропно_й_дивергентной волны Россби (2.1.10) дают

Tmin ~ 7 сут. VgH=200 м/сек, С= <втах/&о« 20 м/сек, L0 = = 1200 км. Как видим, в баротропной модели длина волны вполне сравнима с размерами океанов. Поэтому М. Лонге-

Хиггинс (Longuet-Higgins, 1966), М. Раттри и Р. Чарнелл

(Rattray, Charnell, 1966) и другие исследователи уделили значительное внимание исследованию собственных квазигеострофических колебаний в замкнутых бассейнах. _В этом

случае уравнение для амплитуды функции тока ф (х, у) сводится к хорошо известному уравнению Гельмгольца при нулевом граничном условии на твердых границах

(у2 + ^)Ф = 0,

*,*=

+ Л .

(2.1.37)

 

\ 2(0 J

g H

 

Мы не будем останавливаться на исследовании собственных колебаний такого вида, а также на вопросах отражения планетарных волн Россби от твердых границ, а скажем только несколько слов о вынужденных решениях.

Наиболее важной проблемой динамики градиентно-вихре­ вых волн в физической океанографии является вопрос о ме­

56


ханизме возбуждения этих волн. Г. Стоммел и Дж. Веронис (Stommel, 1956) рассмотрели действие периодического во времени и пространстве ветра на неограниченный двуслой­ ный океан. М. Лонге-Хиггинс (Longuet-Higgins, 19656)

исследовал вопросы возникновения планетарных волн на основе уравнения

Л_\JL

д'

1

1

,

(2.1.38)

-----

J

V =

rot т

gH ) dt

дх

pH

 

 

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

действующей на слой океана, толщины Н: —г ■, по-види-

РН

мому, слишком груба при описании неустановившихся дви­ жений и совсем не годится в области высоких частот изменения ветра, соответствующих периодам, меньшим поло­ вины маятниковых суток. Во всяком случае численные оценки величины возбуждаемых скоростей имеют весьма ориентиро­ вочный характер. Более подробный обзор некоторых послед­ них работ дан в статье А. С. Монина, В. М. Каменковича, Р. В. Озмидова, Б. А. Тареева (1966). Большие величины возбуждаемых скоростей течений могут быть получены в слу­ чаях резонанса (Pedlosky, 19656). Однако эти случаи кажутся слишком специфичными, чтобы ими объяснить ши­ рокое распространение неустановившихся течений значитель­ ной амплитуды. Более того, мы покажем, что собственные пространственные и временные масштабы систем движения в океане существенно отличаются от таковых в атмосфере, что в среднем уменьшает возможности резонансного взаимо­ действия.

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

Как мы увидим, бароклинные градиентно-вихревые воллы, которые лежат в основе этого механизма, будут сущест­ венно изменены вследствие взаимодействия с полем средней циркуляции. Во-первых, и это самое важное, их амплитуда будет возрастать по мере распространения, а скорости рас­

пространения

и периоды

изменятся в очень большой

степени.

 

 

57


§ 2.2. Постановка задачи о бароклинной неустойчивости в океане в квазигеострофическом приближении

. с учетом сил горизонтального турбулентного трения

Проблема устойчивости крупномасштаб­ ных волновых квазигеострофических возмущений возникла в метеорологии (Charney, 1947; Eady, 1949; Kuo, 1952, 1953; Fjortoft, 1950; и др.) как основа для объяснения физической природы циклогенеза, рассматриваемого в качестве некото­ рого вида гидродинамической неустойчивости (главным об­ разом, бароклинной неустойчивости).

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

Как правило, в

качестве основного движения

рассматрива­

лось зональное

геострофическое течение, не

зависящее от

широты,- с постоянным вертикальным сдвигом

(градиентом)

скорости. Поле возмущений предполагалось квазигеострофическим и также не зависящим от широты, влиянием турбу­ лентного трения и теплопроводности пренебрегалось. Кроме того, различные авторы в процессе решения использовали другие упрощения. Так, например, Р. Фьортофт (Fjortoft, 1950) и Дж. Холомбо (Holomboe, 1959) предполагали от­ сутствие вертикальной 'статической устойчивости (верти­ кальной стратификации), Э. Иди (Eady, 1949) и Р. Фьор­ тофт (Fjortoft, 1950) пренебрегали влиянием р-эффекта.

Во всех случаях задача устойчивости рассматривалась не как задача с начальными условиями (задача Коши), а рас­ сматривались решения вида элементарных волн, аналогичных

(2.1.3)

 

 

ф(х, у, z) — (у, z)elk{x

ct).

(2.2.1)

В .этом случае процедура метода

возмущений

приводит

(в предположении независимости решения от одной из про­ странственных координат у или z) к задаче на собственные значений для ©быкновенного несамосопряженного дифферен­

циального уравнения, где амплитудный множитель ср являет­ ся неизвестной функцией. Волновое число k по определению считается действительным, а фазовая скорость волн с (или частота со) играет роль собственного значения задачи и

может, вообще говоря, принимать комплексные значения, т. е.

с — сг -}“ ici.

(2.2.2)

58


Из выражений (2.2.1) и (2.2.2) видно, что в этом случае величина сг сохраняет смысл фазовой скорости волны, а ве­ личина an= Cik — коэффициент при мнимой части комплекс­ ной частоты может быть назван коэффициентом возрастания, так как он характеризует степень экспоненциального изме­ нения амплитуды волны по мере ее распространения. Очевид­ но, если а - > 0 , то амплитуда волновых возмущений экспонен­ циально нарастает и такие решения соответствуют неустой­

чивости. Если с ,< 0 — имеют место затухающие

колебания.

Случай действительных собственных

значений с* = 0 соответ­

ствует нейтральным колебаниям, т.

е. волнам,

распростра­

няющимся без изменения амплитуды. Все рассмотренные в предыдущем параграфе примеры соответствуют как раз таким волнам, т. е. все дисперсионные уравнения дают только действительные значения для с при любых действительных k. Это объясняется тем, что все рассмотренные в § 2.1 модели были слишком просты и не содержали никаких источников энергии для роста волновых возмущений и, следовательно, никакого механизма для развития возмущений. Если, нако­

нец, сг= сг = 0, но ц>фО в (2.2.1), то имеет

место установив­

шееся вторичное течение, наложенное на основной поток.

за

Хотя общая теория гидродинамической

устойчивости

последнее время приобрела достаточную

физическую

яс­

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

или антиклонический вихрь, и поэтому на поздней стадии циклогенеза возможно только качественное сопоставление.

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

59