![](/user_photo/_userpic.png)
книги из ГПНТБ / Тареев, Б. А. Динамика бароклинных возмущений в океане
.pdfных приливных внутренних волн не существует. Собственные решения однородного уравнения (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