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

асимптотически, так как из формулы (2.3.25) следует, что

соi= Cik-+0 при &-Ч) для всех значений Um.

Численные расчеты диаграмм устойчивости, произведенные выше, находятся в удовлетворительном качественном соответ­ ствии с материалами наблюдений. Однако до последнего вре­ мени не было возможности провести последовательное коли­ чественное соответствие между результатами теоретических моделей и основными параметрами крупномасштабных возму­ щений в Гольфстриме и других системах течений. Лишь не­ давно обработка детальных съемок меандров Гольфстрима от м. Гаттерас до 60° з. д. (э. с. «Эксплорер», 1965—1966 гг.) позволила Д. В. Ханзену t(1970) провести такое сопоставле­ ние. Результаты приводятся ниже в таблице, заимствованной

* Дюксбери рассмотрел только нейтральные возмущения (устойчивые волны). Фазовые скорости вычислены для первых трех мод в диапазоне длин волн, представляющем интерес.

** Значения в скобках соответствуют вязкой модели для v = 5-107 см2/сек, величины без скобок вычислены при v = 0 .

Д. Ханзен высказывается определенно в пользу моделей бароклинной неустойчивости (1965, 1969) по сравнению с баротропной моделью Гаурвитца и Пановски и «квазибаротропной» моделью Липпса, в которых источником роста воз­ мущений является кинетическая энергия среднего движения.

Численные значения, приведенные в последней строке таб­ лицы, были получены э/с «Эксплорер» в 1965—1966 гг. во время детальных съемок района меандрирования Гольфстри­

ма ниже м. Гаттерас. Положение изотермы 15°С на глубине 200 м принималось в качестве индикатора стрежня Гольф­ стрима. На рис. 18, взятом из работы Д. Ханзена (1970), по­ казана типичная картина меандров для нескольких последо­ вательных моментов времени.

93

Как видно из таблицы, введение горизонтальной вязкости в модели автора (1965) существенно улучшает количествен­ ное совпадение теоретических оценок с данными наблюдений. Это объяснимо физически, так как. согласно данным наблю­ дений в рассматриваемом районе (Webster, 1961) горизон­ тальные напряжения Рейнольдса в большинстве случаев пере-

Рис. 18. Положение стрежня Гольфстрима по Д. В. Ханзену (1970) для трех последовательных периодов времени, э/с «Эксплорер», 1965 г.

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

Другая группа моделей (Warren, 1963; Robinson, 1967),

рассматривающая Гольфстрим как свободную инерционную струю, подверженную влиянию топографии дна и изменения параметра Кориолиса (сохранение потенциального вихря вдоль траектории), как отмечает Д. Ханзен, дает удовлетво­ рительное описание среднего пути Гольфстрима, но не в со­ стоянии объяснить фактически наблюдаемую структуру меаш Дров.

Все модели дают завышенный коэффициент временного возрастания доминирующей волны, так как, например в>;=

94

—с,к= 10-6сек-1 соответствует удв_оению амплитуды за чет­ веро суток, тогда как фактически время удвоения, по-види­ мому, не менее 10 дней. Это очевидный недостаток линейных моделей, так как при достаточно больших амплитудах ско­ рость роста должна существенно отклоняться от экспонен­ циального за счет нелинейных эффектов. Тем не менее коэф­

expi(kx—'(at) волновое число принимается для простоты веще­ ственным, однако М. Гастер (1962) показал, предполагая ана­

ki = &i(dioT/dkr) - 1. Эта формула позволила Ханзену оценить величину ki в моделях, где k предполагается вещественным.

Значения Ld=2nlkr, со,, ki, полученные автором в вязкой модели ,(1965), находятся в хорошем соответствии с данными наблюдений (величины, заключенные в скобки в табл. 1). Однако значения фазовой скорости, полученные Орлански (1969), который в более полной модели учел влияние топогра­ фии дна, лучше согласуются с данными наблюдений. Посколь­ ку в рассматриваемом районе относительные изменения глу­ бины невелики, более заметное влияние на фазовую скорость возмущений должна оказывать существенная разница между толщиной главного термоклина и слоя глубинных холодных вод. В этом смысле двуслойная модель с разными толщинами слоев может лучше соответствовать реальности, нежели мо­ дель с постоянным сдвигом скорости по вертикали. Кроме того, в некоторых районах океана внутреннее число Фруда (величина, обратная числу Ричардсона в двуслойной модели)

Двуслойная и двухуровенная модели

В этом параграфе мы следуем в основ­ ном работам (Тареев, 1967, 19686).

Как уже было указано, большие значения вертикального сдвига скорости должны приводить к интенсивному развитию

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

дениями, в частности на участке Гольфстрима от мыса Гаттерас до большой Ньюфаундлендской банки. Между тем на

95

участке Гольфстрима от Флориды до м. Гаттерас наблюдения не отмечают заметного развития крупномасштабных волновых возмущений или вихрей, несмотря, на очень большие (по сравнению с критическим) значения вертикального сдвига скорости. Гольфстрим, проходящий на этом участке над об­ ластью материкового плато, представляет собой интенсивную узкую струю (ширина менее 100 км), ограниченную со сто­ роны берега и со стороны океана сравнительно медленно дви­ жущимися водами (см. рис. 8), в которых изопикнические поверхности практически горизонтальны. Ввиду этого, для описания свойств устойчивости Гольфстрима на этом участке мы должны принять во внимание конечную ширину течения. Кроме того (это будет ясно из дальнейшего), характерные параметры Гольфстрима на этом участке не позволяют ис­ пользовать при описании поля возмущений обычное прибли­ жение квазигеострофичности.

Поскольку Н. Филлипсом (Phillips, 1951) было показано, что простая двуслойная модель бароклинной неустойчивости правильно описывает многие важные свойства непрерывных моделей, ниже будет рассмотрена двуслойная модель основ­ ного течения, однако поле возмущений не будет предпола­ гаться квазигеострофическим. Если в метеорологии использо­ вание двуслойных (или двухуровенных) моделей имеет в зна­ чительной мере условный характер и диктуется в основном соображениями простоты, то в океанографии в связи с суще­ ствованием океанского термоклина разделение океана на два слоя может быть проведено более определенным образом Г Учет фактического различия в толщине термоклина и под­ стилающих холодных вод позволит получить, в частности, лучше согласующиеся с наблюдениями значения фазовых ско­ ростей неустойчивых возмущений по сравнению с соответст­ вующими значениями, полученными в § 2.5 в двухуровенной модели с постоянным вертикальным сдвигом скорости.

Рассмотрим двуслойное течение, направленное вдоль оси л: ширины I (рис. 19), в которой поперечный наклон поверхно­ сти раздела определяется формулой Маргулеса:

РР /

Индексы 1, 2 здесь и далее относятся соответственно к ниж­ нему и верхнему слою, р — среднее значение потенциальной плотности. При учете изменения параметра Кориолиса / с1

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

96

широтой будем, как и выше, использовать приближение р- плоскости, считая, что ось у направлена на север (f= /o +

+$уУ-

2

Рис. 19. Двуслойное течение с разными глубинами слоев, ограниченное по широте

Уравнения возмущенного движения, линеаризированные относительно основного течения (2 .6.1 ), имеют вид:

предположено, что свободная поверхность верхнего слоя не­ подвижна (это предположение отфильтровывает баротропные длинные волны, не представляющие здесь интереса). Урав­ нение (2.6.4) является следствием уравнения гидростатики и непрерывности полного давления на поверхности раздела. Со­ ставляя из (2 .6.2) уравнение вихря п подставляя значение го­

ризонтальной дивергенции из уравнения неразрывности (2 .6.3), получим

1 Там, где нет дифференцирования по у, будем считать / = / 0.

i i a ' - р » ; = ег ’

98