Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Тареев, Б. А. Динамика бароклинных возмущений в океане

.pdf
Скачиваний:
6
Добавлен:
20.10.2023
Размер:
7.86 Mб
Скачать

ных приливных внутренних волн не существует. Собственные решения однородного уравнения (3.3.15), удовлетворяющие граничным условиям (3..1.4) и (3.2.2.), соответствуют рас­ смотренным выше свободным внутренним волнам. Вынуж­ денные приливные внутренние волны соответствуют частно­ му решению неоднородного уравнения (3.3.15). Чтобы пере­ менные в уравнении (3.3.15) разделились, надо предполо­ жить, что это решение имеет вид

 

 

w = ср (г)

 

.

(3.3.16)

Тогда для ф (z)

получим

 

 

 

d2ф

+ k20 Л^ ~ (00

=

«рЩК N2

(3.3.17)

1z2~

“ о — / 2 .

(«о — / 2) g

 

 

 

 

В (3.3.17)

можно положить

©0 =

0, что соответствует

гидростатическому приближению. Частное решение неодно­

родного

уравнения (3.3.17) проще

всего

найти, положив

const

(физические выводы, впрочем, не изменятся, если

задаться более сложным законом изменения N с глубиной).

В этом

случае

решение

неоднородного

уравнения

(3.3.17)

дается рядом

 

 

 

 

 

 

 

 

 

л=0

 

 

 

(3.3.18)

 

 

 

 

 

 

 

коэффициенты

которого

определены

соотношениями

 

 

 

Ап = -----------^

------------ .

(3.3.19)

 

 

 

 

- h ( h n ! H f

 

 

Здесь

B = k0a>oaoN2lg;

Сп — коэффициенты разложения

единицы

в ряд

по sin(hnlH)z,

которые

при малых

п, оче­

видно, имеют порядок единицы. Интересно оценить порядок величины w в средних широтах, вдали от резонанса. Харак­ терные порядки величин, входящих в (3.3.19), а0= 10~4 см/сек2;

N2= 10-2— 10-3 сек-2, g = 103 см/сек2, &0=Ю -8 см-1, ©о~ f = = 10~4 сек-1, я2/Я2=10-9 см-2 дают для амплитуды верти­

кальной

компоненты скорости внутренних приливных волн

ничтожно

малое значение: 10~5'— 10-6 см/сек, которое при­

мерно на два порядка меньше характерного значения w для баротропных приливов в открытом океане. Таким образом, можно сказать, что под действием приливообразующих сил океан ведет себя как баротропная жидкость. Значения Ат могут стать достаточно большими вблизи «резонансных» широт (где ©о = f). Как известно, для полусуточного и суточ­ ного лунного прилива эти широты приближенно равны 76°

150

2

и 30°. Дополнительный член koN2 в знаменателе (3.3.19) очень мал и лишь незначительно уменьшает численное зна­ чение резонансной широты. Простые расчеты с помощью резонансной кривой показывают, что область, где w увели­ чивается в 103 раз и достигает значения 10-2— 10_3 см/сек, очень узка (порядка десятых долей градуса), так что дис­ сипативные процессы, как бы малы они ни были, по-види­ мому, полностью гасят резонансные эффекты. Наблюдения над обычными поверхностными приливами показывают, что на «резонансных» широтах не имеется каких-либо особен­ ностей в поведении колебаний уровня, но есть заметное уси­ ление Приливных течений. Это усиление легко объяснить, если предположить, что движение происходит почти без

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

и •= —iа°?°—

v =

SsL— e£(ft„*-<o.O

(3.3 20)

f-CD2

 

Д-СО2

 

Поскольку движение

происходит под действием массо­

вых сил и без градиентов

давления, плотность

вообще не

входит в формулы

(3.3.20). Порядок и, v в средних широтах

1 см/сек, поэтому

резонансное

увеличение горизонтальной

скорости только в 10 раз приводит к значительным и впол­ не измеримым значениям горизонтальной скорости.

Выводы о невозможности возникновения заметных внут­ ренних волн приливного периода в открытом океане посто­ янной глубины по существу уже содержались в работе (Defant, 1950), где рассматривались простые двуслойные модели. Имеются более поздние работы, претендующие на теоретическое доказательство возможности существования таких волн, однако детальный анализ применяемых рассуж­ дений показывает, что эти доказательства основаны на недо­ разумениях. Так, например, в работе (Krauss, 1959) собст­ венным решениям произвольно приписывается частота при­ ливообразующей силы, тогда как эту частоту должны иметь вынужденные решения.

Следует, однако, заметить, что по наблюдениям в некото­ рых районах океана, во всяком случае вблизи побережий (Рид Дж., 1956), внутренние волны приливного периода, по-видимому, имеют место. Возможность возникновения внутренних волн приливного периода на береговом шельфе была показана Раттри (Rattray, 1960). Если шельф доста­ точно широк, то эти волны являются почти плоскими и могут распространяться на значительные расстояния внутрь океана без существенного уменьшения амплитуды. Другой возможный механизм'—образование стоячих приливных

151

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

туда, очевидно, будет уменьшаться, по крайней мере как

__i_

г 2 , где г — расстояние от источника. Следует заметить, что данные о фактических частотах внутренних волн, полу­ ченные по наблюдениям в точке, могут быть значительно

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

§ 3.4. Возникновение внутренних волн при обтекании неровностей дна. Приложение к морской геологии

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

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

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

Поскольку в глубинных слоях океана вертикальные гра­ диенты температуры очень малы, вместо плотности удобно

152

ввести

в рассмотрение

энтропию, что позволит

наиболее

естественным путем рассмотреть изотермический случай.

Ряд

задач такого

рода (см., например, А. А. Дородни­

цын, 1938,

1950)

был

ранее

рассмотрен

в метеорологии.

Наряду

с

этим

такая

задача

представляет

также

океано­

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

типа внутренних волн

в

абиссальных областях океана.

С другой стороны,

эта

задача позволяет объяснить неко­

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

Экспедиционные исследования последних лет, в частно­ сти работы «Витязя» (отчет 34-го рейса, 1961), обнаружили тесную зависимость мощности донных отложений от харак­

тера рельефа дна.

Так, например, для холмистого

рельефа

с высотой холмов,

на порядок меньшей глубины

океана,

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

Холмистый, приближенно синусоидальный рельеф дна первоначально вулканического происхождения, прикрытый впоследствии слоем осадков, довольно часто встречается в океане, в частности в тропических частях Атлантики и Ти­ хого океана. «Длина волны» таких периодически повторяю­ щихся гряд порядка 1 км, высота 100—300 м при глубине океана 3—6 км. При таких относительно малых изменениях глубины изменения горизонтальной составляющей скорости морского течения малы и не могут привести к неравномер­ ному смыванию осадков с вершин холмов. Однако, как бу­ дет показано, возникающее при обтекании таких гряд хол­ мов поле вертикальных составляющих скорости может привести к существенному перераспределению осадочных

153

частиц, равномерно выпадающих из поверхностных слоев океана.

В метеорологии задача обтекания неровностей подсти­ лающей поверхности бароклинным потоком была достаточ­ но подробно исследована многими авторами. При рассмот­ рении океанографического варианта этой задачи мы будем следовать в основном работам автора (Тареев, 1964, 19656), в которых применен метод А. А. Дородницына (1938, 1950), но изменим соответствующим образом постановку задачи. Расположим начало координат на свободной поверхности океана и направим ось г вертикально вниз. Пусть форма неровностей дна и все движение не зависят от у и постоян­ ная скорость невозмущенного течения и направлена вдоль оси х. Считая морскую воду неоднородной несжимаемой жидкостью и пренебрегая эффектом вязкости, который мо­ жет быть существенным только в тонком придонном погра­ ничном слое (и в данном случае не представляет интереса), запишем уравнения задачи (считая движение квазистационарным и пренебрегая параметром Кориолиса)

, ди . ди \ др

 

р

и

ow

 

W-

 

 

др

 

gp-

(3.4Л)

 

~д7

 

 

 

дг

 

 

 

 

 

дг

 

 

 

 

 

 

 

 

 

- ^ + * 1

= 0,

 

 

 

 

 

 

 

 

 

дх

дг

 

 

 

 

 

 

 

 

и

ds

 

ds

 

г.

р =

,

ч

 

 

 

 

— +

w—

= 0 ,

р (s).

 

 

 

 

 

дх

 

дг

 

 

 

 

 

 

 

Здесь

и, w — составляющие

скорости

вдоль

осей х, z ;

р — плотность;

s — удельная

энтропия

морской

воды;

р — давление; g

— ускорение силы тяжести.

 

 

Пусть теперь:

 

 

 

 

 

 

 

 

 

 

 

 

и =

и +

и'

(х, z),

р = р0 (2) +

Р' (х, г),

 

 

 

w =

w' (х,

z),

s =

s0 (г) -f- s' (х, z).

 

(3.4.2)

 

 

 

 

Р =

Ро (2) +

р' (х, z).

 

 

 

 

Здесь

величины

со

штрихами — возмущения

основного

течения,

вызванные

обтеканием

малой

неровности

дна

| = 1(х),

причем | ( —оо)=0,

так что при х->----оо эти возму­

щения обращаются в нуль. Учитывая уравнение гидроста­ тики из (3.4.1) и (3.4.2), пренебрегая величинами второго порядка малости, получим:

154

 

dU'

 

dp’

p0u

dw'

dp'

gP

Ро^ dx

 

дх

дх

dz

 

d£_

>dsq

 

du'

dw'

(3.4.3)

u

 

= 0.

dx

W

dz

 

dx

dz

Считая в соответствии с вышесказанным, что изменения плотности связаны с йзменениями энтропии, но не давления, запишем:

Р =

Фо

ps •

(3.4.4)

dsn

 

 

 

Из уравнения (3.4.3) с учетом (3.4.4) получим уравнение для w'\

d2w' .

d2w |

1

dp0

dw' , klw' = 0.

(3.4.5)

dx2

dz2

p0

dz

dz

 

Здесь обозначено:

*g = _ I _ Y J L \ J*l .

Po^2 \ ds0 ) p dz

Учитывая, что относительные изменения плотности в океане порядка 10~2— 10~3, и сравнивая порядок членов в уравнении (3.4.5), найдем, что величина третьего члена в (3.4.5) по крайней мере на два порядка меньше остальных и им можно пренебречь1. Таким образом, с высокой сте­ пенью точности уравнение (3.4.5) можно переписать в виде:

d2w' .

d W i , 2 /

п

/о л c\

— — +

— — + k0w

= 0 .

(3.4.6)

дх2

ду2

 

 

Используя известные термодинамические соотношения с учетом уравнения гидростатики

I

Фо \

Т0

/

гф 0

\

ds0

\

dso 1 р

ср

\

dT

/ р

dz

перепишем выражение для kl:

dTa

Ро“2 L V дТ0 ) p dz

cp

- +

_g_

(

дро

 

dT0

 

\

дТ0

 

T0

dz

Po

р

Л- Л

(. Фо

 

>

0.

(3.4.7)

Pocp V dT

/ p.

 

 

 

Здесь

70(г )— равновесное

значение

абсолютной

темпера­

туры:

изменения

солености морской

воды не принимаются

во внимание;

ср — теплоемкость при

постоянном давлении.

В

качестве

граничного

условия для w' на дне имеет

обычное кинематическое

соотношение

w = и ----,

линеа-

 

 

 

 

 

 

dx

 

1 Заметим, что это приближение опять-таки

в точности

соответст­

вует приближению

Буссинеска,

неоднократно

использовавшемуся выше.

155

ризируя которое и принимая во внимание малость £ по сравнению с Н — глубиной океана, получим:

w' = и — при zzzH.

(3.4.8)

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

хю' = 0 при г = 0,

(3.4.9)

кроме того, w' (—002),

w' (+ 002) — ограничено.

(3.4.10)

Как известно, в основной толще вод океана, за исклю­ чением приповерхностного слоя главного термоклина тол­ щиной порядка 500—1000 м, температура меняется крайне медленно и это изменение можно считать линейным. Поэто­ му, полагая в (3.4.7) dT0/dz=const и заменяя р0 его средним

значением, получим, что kl = const.

Положим (отбрасывая в дальнейшем штрих у а/), что

w = - ^ 7 U ^ r + 4>(x,z),

(3.4.11)

лах

тогда, учитывая (3.4.6) и неоднородное граничное условие (3.4.8), получим для <р неоднородное уравнение:

(Э2ф

Э2ср

&оф = ■ U ■

d%

+ k20 J L

(3.4.12)

5л:2

 

Н

\ dx3

 

dx

 

с граничными условиями:

 

 

 

 

Ф (х, 0) = ф (х, Н) = ф(— оо, z) — 0,

 

(3.4.13)

 

 

Ф (оо, z) —ограничено.

 

Удовлетворяя граничным условиям по z, положим

 

 

 

О

 

 

 

 

 

ф =

5 ] Ф» (*) sin Y„z

 

 

 

(3.4.14)

 

 

7 1 = 1

 

 

 

 

Подставляя

(3.4.14) в (3.4.12)

и разлагая

правую часть

в ряд по sinyn^, получим для ф„(х) уравнение:

 

 

42<Р/г

+ (kt— yl) Ф„ = исп

+ kl

 

(3.4.15)

dx3

 

 

dx3

'dx

 

156

где

2

(3.4.16)

с„ = (— 1)" пп

Интеграл этого уравнения, удовлетворяющий

условиям

(3.4.13) по х, будет

X

Фя = р р = = 2 [ Sin [К&о — у2(х — £)] х

(3.4.17)

при kl — у«> 0,

при у» —&о > 0■

(3.4.18)

Пусть теперь

 

 

 

| 0

при

—оо < х < 0

(3.4.19)

( g

sinp*

при 0 < х < [о о .

 

Как покажут приведенные

ниже численные

оценки для

интересующих нас значений параметров, слагаемыми, кото­

рые

даст формула (3.4.18),

молено в

первом приближении

пренебречь,

так как они образуют * знакопеременный ряд,

члены которого убывают во всяком

случае не

медленнее,

чем

1In.

Решения

вида

(3.4.18)

появляются,

когда

Упко >0,

т. е. (как

будет

видно ниже)

при

достаточно

больших п.

Поэтому,

учитывая появление

сильно

убываю­

щих экспоненциальных множителей и рассматривая движе­ ние для достаточно больших х, ограничимся рассмотрением выражения (3.4.17).

Подставляя (3.4.19) в (3.4.17), выполняя квадратуры и простые тригонометрические преобразования, получим

Учитывая (3.4.11). (3.4.14), (3.4.16), (3.4.19), получим

(3.4.21)

157

Причем N соответствует последнему значению п, при кото­

ром все еще выполняется неравенство

kt — у« >0.

Уместно написать еще одно решение для случая, когда

рельеф дна задан функцией:

 

 

 

 

 

 

 

 

 

 

0

при

 

х < 0

 

 

 

 

£0/2 при

х —0

(3.4.22)

 

 

 

1£

при

 

х > 0

 

Тогда, очевидно, £' (х) = £8),

где б(х)

— функция Ди­

рака. Используя основное свойство б-функции:

ъ

 

 

 

 

 

 

 

 

 

J / (X -

 

1) 6<л) dl = ( -

1)"/<"> (х),

а < х < Ь,

а также формулы (3.4.11), (3.4.17), (3.4.18), получим:

w = ы£0

— б (х) + У

- ^ п-

sin V

kl — ylx х

 

Н

n=i

V k 2 v2

 

 

 

 

 

г к0

Yo

 

 

X sin

-

 

c«Yre

■sin

 

Y„xsm y„2.

 

 

Yo

 

 

 

 

л=ЛГ+1

 

 

 

 

(3.4.23)

 

 

 

 

 

 

 

 

 

Эта формула имеет особенность при х=0, что объяснимо физически, так как в окрестности этой точки линейная тео­ рия неприменима. Первая сумма в (3.4.23) должна вычис- I

ляться д л ях > 0,

сходимость

бесконечной

суммы в

(3.4.23)

при |х |=^=0

обеспечивается

наличием

экспоненциальных

множителей.

Для

х > 0 можем вместо (3.4.23) написать:

 

N

с пУ2п

•sin Vm- Y«xsinY„2.

 

W =

 

(3.4.24)

 

 

П~\ V k i ­ l l

Теперь проведем численные оценки, основываясь на ре­ шении (3.4.21), достаточно близко соответствующем реаль­

ным условиям.

___ Отг

Пусть X = — = 2 км, т. е. 0 = 3,14-10—5 см-1. Как из-

Р

вестно, за исключением верхнего слоя главного термоклина (500—1000 м), температура в океане изменяется с глубиной очень слабо. Поэтому для оценки положим dT0fdz и возьмем согласно (3.4.7) минимальное значение k0:

kl = —-— (

) л Г И .

Рйи 2 \

дТ0 J p V ср

158

Полагая,

что и 10 см/сек и имея в

виду,

что в CGS

g=103,

ро=1,

(ф 0/(37’о)р—1,5-10-4,

Т0= 2,7-102 гр,

Ср= 4 -107

эрг/град, получим

~ 0,4• 10—4 см-1.

получим в

Если

глубина

океана

Я = 5 км = 5-105 см,

(3.4.21) 7V= 6.

 

 

 

 

 

Распределение вертикальных составляющих скорости,

вычисленное

по

формуле

(3.4.21), при

вышеприведенных

значениях

параметров задачи показано на

рис.

34, где чис-

Рис. 34. Распределение вертикаль­ ных составляющих скорости над синусоидальным дном

ленные значения для получения размерных величин долж­ ны быть умножены на w0 = u$^o.

Если £

порядка 100 м =

104

см, то, очевидно, w =

= 3,1 см/сек. Горизонтальные

и

вертикальные

расстояния

на рисунке

даны соответственно

в

безразмерных

величинах

§х и z/H. Как видно из рисунка, при вышеприведенных зна­ чениях параметров, порядок величины вертикальной состав­ ляющей скорости в основной толще океана — 1 см/сек. Это значение намного превышает порядки вертикальных скоро­ стей, обычно встречающиеся при рассмотрении осредненной,

крупномасштабной

океанической

циркуляции

(10_3—

10~4 см/сек). Однако

легко видеть,

что рассматриваемые

бароклинные волны не дают вклада в результирующий вер­ тикальный водообмен между достаточно удаленными точка­ ми по вертикали ввиду квазипериодического характера дви­ жения. В самом деле, мы здесь имеем типичный случай стационарных внутренних волн, в которых вертикальные колебания изотерм (порядка 10 м), соответствующие сме­ щениям частиц жидкости, не должны вызывать удивления. Хорошо известно, что скорости вертикального движения во внутренних волнах в общем намного больше вертикальных скоростей в средней крупномасштабной океанской циркуля­ ции. Используя данные вычислений, приведенные на рисун­ ке, легко оценить, что в данном случае величина вертикаль-

159

Соседние файлы в папке книги из ГПНТБ