Научная статья на тему 'Следы цунамигенного землетрясения во вращающемся стратифицированном океане'

Следы цунамигенного землетрясения во вращающемся стратифицированном океане Текст научной статьи по специальности «Физика»

CC BY
45
7
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЛИНЕЙНАЯ ТЕОРИЯ ДЛИННЫХ ВОЛН / LINEAR LONG-WAVE THEORY / ГЕНЕРАЦИЯ ЦУНАМИ / TSUNAMI GENERATION / КОСЕЙСМИЧЕСКАЯ ДЕФОРМАЦИЯ ДНА / CO-SEISMIC BOTTOM DEFORMATION / ОСТАТОЧНЫЕ ГИДРОДИНАМИЧЕСКИЕ ПОЛЯ / RESIDUAL HYDRODYNAMIC FIELDS / ГЕОСТРОФИЧЕСКИЙ ВИХРЬ / GEOSTROPHIC VORTEX / СТРАТИФИКАЦИЯ / РАДИУС ДЕФОРМАЦИИ РОССБИ / ROSSBY RADIUS OF DEFORMATION / TWO-LAYER STRATIFICATION

Аннотация научной статьи по физике, автор научной работы — Носов Михаил Александрович, Нурисламова Гульназ Нуровна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Следы цунамигенного землетрясения во вращающемся стратифицированном океане»

Следы цунамигенного землетрясения во вращающемся стратифицированном океане

М. А. Носовa, Г.Н. Нурисламоваb

Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра физики моря и вод суши. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.

E-mail: a nosov@phys.msu.ru, bgulnaz1205@yandex.ru

Статья поступила 02.07.2013, подписана в печать 06.09.2013.

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

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

УДК: 551.466. PACS: 91.30.Nw, 92.10.Iv.

Введение

Косейсмические (остаточные) деформации дна океана, сопровождающие сильные подводные землетрясения, приводят к вытеснению значительных объемов воды [1, 2]. Вытеснение воды служит причиной формирования не только волн цунами [3, 4], но и долгоживу-щих «следов» цунамигенного землетрясения в океане — потенциального и вихревого остаточных гидродинамических полей [5]. Потенциальное остаточное поле — горизонтальные смещения частиц воды — возникает в результате «растекания» вытесненного объема под действием силы тяжести. Из-за влияния на этот процесс вращения Земли образуется вихревое остаточное поле — геострофический вихрь.

На возможность формирования вихревых остаточных полей при цунамигенных землетрясениях обратили внимание еще в конце XX в. [3, 6-9]. Потенциальные остаточные поля оставались без внимания вплоть до выхода наших работ [5, 10, 11].

В работе [11] нами было получено аналитическое решение, описывающее потенциальное и вихревое остаточные поля, которые возникают во вращающемся однородном океане постоянной глубины в результате вытеснения воды цилиндрически симметричной деформацией дна. Анализ этого решения показал, что в типичных природных условиях доминирующим является не вихревое, а именно потенциальное поле. Соответствующие ему горизонтальные смещения частиц воды достаточно велики (~ 100 м), следовательно, они могут быть зарегистрированы in situ и использованы для восстановления структуры источника и, возможно, для предупреждения о цунами. Что касается вихревого поля, то оно обычно является настолько слабым, что регистрация in situ скорости течения (~ 0.01 м/с) или смещения свободной поверхности (~ 0.01 м) в геострофическом вихре представляется малореальной.

Кроме того, в работе [11] было установлено, что в рассматриваемой системе в принципе возможны ре-

жимы, когда доминирующим оказывается не потенциальное, а вихревое остаточное поле. Эти режимы интересны тем, что при их реализации основная доля энергии, передаваемая водному слою сейсмическим источником, переходит не к гравитационным волнам (цунами), а к геострофическому вихрю. Иными словами, начальное возвышение свободной поверхности, образующееся в очаге цунами в результате вытеснения воды остаточной деформацией дна, связывается в геострофическом вихре. Впрочем, для реализации таких экзотических режимов необходимо, чтобы горизонтальный размер источника (очага цунами) заметно превосходил баротропный радиус деформации Россби, типичное значение которого ~ 106 м [12]. Выполнение этого условия на нашей планете практически невозможно. Дело в том, что только при очень сильных землетрясениях размер очага цунами приближается к баротропному радиусу деформации Россби; в большинстве случаев размер очага составляет ~ 105 м [1, 3, 4].

Но в типичных земных условиях эффект связывания начального возвышения в очаге цунами может реали-зовываться для возмущений устойчивой стратификационной структуры океана, вызванных сейсмическими движениями дна. Действительно, в этом случае горизонтальный размер очага цунами следует сопоставлять не с баротропным, а с бароклинным радиусом деформации Россби (~ 104 м) [12], который, как правило, существенно уступает размеру очага цунами. Следовательно, вращение Земли должно оказывать сильное влияние на эволюцию возмущения стратификационной структуры, препятствуя его распаду на серию внутренних волн. В этой связи картина эволюции цунами в стратифицированном океане, описанная в классической работе [13] без учета вращения Земли, может оказаться не вполне корректной. Впрочем, как следует из заключительного раздела этой работы, автор, несомненно, осознавал важность «эффекта Кориолиса».

Возникновение возмущений стратификационной структуры океана при подводных землетрясениях связа-

но не только с вытеснением воды при деформации дна, но и с интенсификацией вертикального обмена в океане над очагом подводного землетрясения [4, 14, 15]. В любом случае долгоживущий след цунамигенного землетрясения, несущий в себе информацию об источнике цунами, несомненно, представляет как теоретический, так и практический интерес.

Настоящая работа является логическим продолжением нашего исследования [11], выполненного в предположении однородного океана. Здесь мы расширяем постановку задачи, предполагая, что океан является устойчиво стратифицированным (двухслойным). Цель работы — получение аналитического решения модельной задачи и выявление особенностей остаточных полей, обязанных наличию стратификации, при условиях, свойственных реальным очагам цунами.

1. Математическая модель

Будем рассматривать безграничную вдоль горизонтальной плоскости несжимаемую двухслойную жидкость на вращающейся Земле. Постановка задачи схематически представлена на рис. 1. Верхний слой имеет

Рис. 1. Постановка задачи

толщину Н1 и плотность р, а нижний слой — толщину Н2 и плотность р + 5р (5р >0). Сферичностью Земли пренебрежем. Начало прямоугольной системы координат расположим на свободной невозмущенной поверхности верхнего слоя. Ось г направим вертикально вверх, а оси х и у — на восток и на север соответственно. Для описания движений жидкости применим уравнения линейной теории длинных волн, записанные с учетом силы Кориолиса и предположения о малости перепада плотностей между слоями (5р/р с 1)

дщ . f

It _ -g Их + [U1'

dvi d£i

if = -g ~ду -fui'

^ _ -g (^ + ^^ + fv2,

dt g \дх p дх ) 1 2' dv2 (д^1 Sp д& \ f

ИТ _ -ЧЪ + - fu2,

ди1 дv1 \ h д^1 д£2 _ ^ дх ду J 1 дt дt '

дщ2 + \ H + дп _ о

дх ду J 2 дt дt '

(1)

(2)

(3)

где щ и vi — компоненты горизонтальном скорости течения i-го слоя (i _ 1,2) вдоль осей х и у соответственно, i — смещение поверхности i -го слоя от равновесного положения, п — смещение поверхности дна от исходного положения, g — ускорение силы тяжести, f — параметр Кориолиса (f _ const).

Пусть до землетрясения положение дна определяется формулой zb _ -H, и оба слоя находятся в состоянии покоя ui _ vi _ ii _ 0. После землетрясения дно перемещается в новое положение zb _ -H + п»(х, у), где п»(х,у) — остаточная деформация дна (|п»| СH), которая является причиной формирования остаточных гидродинамических полей в рассматриваемой системе.

Подход к решению задачи (1)-(3) аналогичен тому, который был использован в нашей работе [11]. Вначале поле скорости течения представляется как сумма потенциальной и вихревой компонент:

u _ + i дх ду

^=др - (4)

ду дх

где — потенциал, г — функция тока I-го слоя. Компоненты скорости течения, выраженные через потенциал и функцию тока посредством формул (4), подставляются в уравнения (1)-(3), которые затем интегрируются по времени от 0 до то. В итоге получаем систему стационарных дифференциальных уравнений, описывающую потенциальное и вихревое остаточные поля в двухслойной вращающейся жидкости:

Н1ДФ1 = - Н2ДФ2 = П( - (5)

Дфк = [ ДФ¡, (6)

§Д6( +1Д-Ф1Ж = 0, §Д6( + g—+ !Д-Ф2Ж = 0,

р

(7)

где — остаточное смещение I -й поверхности в геострофическом вихре; ф( — функция тока, описываем

ющая остаточное вихревое поле; = | (И — по-

0

тенциал смещений, по которому рассчитывается вектор остаточного смещения частиц воды в горизонтальном направлении = УФ;.

Исключая из уравнений (5)-(7) функции и , приходим к следующей системе уравнений относительно функций и :

_(1 + в)(б» - 6»),

^2Д<?2» - <?2» _ -

1

1 + в

(п» + в&»),

(8) (9)

где Я0 = о0/1 и Я1 = е1/{ — баротропный и бароклинный радиусы деформации Россби [12], с0 = v/g(HГTH2) — скорость длинных гравитационных волн в однородной жидкости, С1 = = V ^рН1Н2/ р(Н1 + Н2) — скорость длинных внутренних волн в двухслойной жидкости [12], в = Н2/Н1 — отношение толщин нижнего и верхнего слоев.

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

По заданной остаточной деформации дна пто из решения системы (8), (9) определяются функции £1то и £2то, по которым с использованием уравнений (5) и (6) рассчитываются все остальные искомые функции Ф и ф^:

V-1« = -g &

= - g ( 6 о + — 6

ф1=- g ^ g

(10)

о + -р<?2о^ , ф2 = -g + ~&2о

-р Р

(11)

Вообще говоря, формулы (10), (11) верны с точностью до функции, являющейся частным решением плоского уравнения Лапласа, например С0 + С1 х + С2у, где С0, С1 и С2 — константы интегрирования. В силу того что искомые величины (остаточное смещение частиц воды и скорость в геострофическом вихре) выражаются через производные 1-го порядка по пространству от функций Ф; и ^¡то, величина константы С0 не имеет значения. А из условия стремления решения к нулю на бесконечности следует С1 = С2 = 0. Для иного частного решения уравнения Лапласа, например С11п(х2 + у2)+С0, получаем, что в силу ограниченности решения в начале координат константа интегрирования С1 = 0, а величина С0, как уже отмечалось, не имеет значения. Кроме приведенных выше двух примеров, уравнение Лапласа, конечно, имеет и другие частные решения [16]. Но во всех случаях по указанным выше причинам константы интегрирования в этих решениях следует выбирать равными нулю.

Аналогично тому как это было сделано для случая однородного океана в работе [11], в качестве простой модели остаточной деформации дна будем рассматривать осесимметричное поднятие, описываемое функцией

П0(Г)= П0 тах[1 - 0(г - Я)], (12)

где п0 тах — амплитуда деформации, Я — радиус источника, 0 — функция Хевисайда.

Непосредственному аналитическому решению система (8), (9) не поддается даже в наиболее простом случае цилиндрической симметрии. Существенно упростить задачу можно, полагая £1то = 0 (приближение «твердой крышки»). Но приближение «твердой крышки» может сильно исказить результат. Конечно, в типичных условиях смещение свободной поверхности в остаточном геострофическом вихре составляет всего порядка 1% от величины остаточной деформации дна: 6то ~ 0.01пто [11]. Но дело в том, что коэффициент в может принимать достаточно большие значения, например в = 50 при Н1 = 100 м и Н2 = 5000 м. Следовательно, величины в£1то и пто, входящие в правую часть уравнения (9), вполне могут оказаться сопоставимыми. Поэтому пренебрегать смещением свободной поверхности было бы некорректно. Приведем и второй, не менее важный аргумент против использования приближения «твердой крышки». Полагая £1то = 0, мы автоматически исключаем из рассмотрения баротропный геострофический вихрь и тем самым кардинально меняем постановку задачи.

Вместо приближения «твердой крышки» мы воспользуемся иным, более точным приближением. Будем полагать, что учет слабой стратификации не может сильно

изменить смещение свободной поверхности в геострофическом вихре. Таким образом, смещение свободной поверхности может быть найдено из решения задачи об остаточных полях для однородного океана [11]. Напомним, что эта задача сводится к неоднородному уравнению Гельмгольца относительно смещения свободной поверхности в геострофическом вихре £то:

я0Д£то - £то = -пто, (13)

для которого удалось найти аналитическое решение в случае осесимметричной деформации дна вида (12).

Подтвердим предположение о слабом влиянии стратификации на смещение свободной поверхности, используя оценку, основанную на аппроксимации оператора Лапласа в уравнениях (8), (9) и (13) по теории размерности: Д ~Я-2. В результате несложных преобразований получаем, что относительное изменение амплитуды смещения поверхности, которое обязано учету стратификации, определяется следующей формулой:

бто - бто Я (Я + Я2(1 + в)) ()

£то ~ (я2 + я02)(Я2 + Я2)+я2я2в' ( )

Подставляя в формулу (14) типичные значения (Я = 105 м, Я0 = 106 м, Я1 = 104 м, в = 10), получаем (£то - ^1то)/^то ~ 0.01. Иными словами, в типичных условиях учет стратификации действительно слабо сказывается на амплитуде смещения поверхности.

В цилиндрических координатах уравнение (13) имеет вид

д2£то , 1

дг*2 + г* Or*

- = -Р'0п'о >

(15)

где г* = r/R — безразмерная пространственная переменная, р0 = R/R0 — безразмерный параметр, выражающий отношение радиуса источника и баротропного радиуса деформации Россби. Решение уравнения (15) для деформации дна, определяемой формулой (12), выражается через модифицированные функции Бесселя [11]

, , п /1 - ß0K\(ß0)k(P0r*), 0 < г* < 1,

?о(г ) = П0 max <

[p0K0(p0r*)Ii (М0), г * ^ 1,

(16)

где I — функция Инфельда, K — функция Макдо-нальда.

Для получения приближенного решения системы уравнений (8), (9) будем полагать, что = , где функция определяется выражением (16). Если функция нам известна, то уравнение (8) для

дальнейшего решения уже не требуется — оно выпадает из рассмотрения. Оставшееся второе уравнение системы (9) также запишем в цилиндрических безразмерных координатах:

1 + ß

(По + ß6o), (17)

д ^2то + 1 д&то 2С =

+ г* - ^1«2то = -

где = Я/Я1 — безразмерный параметр, выражающий отношение радиуса источника и бароклинного радиуса деформации Россби. Уравнение (17) фактически аналогично уравнению (15), отличие состоит только в значении коэффициентов и виде правой части. Решение уравнения (17) дается следующей формулой:

6то(г*) = 1+|(С1(г*) + в(2(г*)) , (18)

0(г*) =

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1 - — 1 К1(—ОМ—1Г*), VlKo(vlr*)h (—1),

0 < Г* < 1, г* ^ 1,

(19)

1 +

—0—1 (——1/0(^—0)^ (—0) + —0/0(^—1 №(—1))

2 2 — 1 — —0

С2(Г*) =

0 < г* <1,

—0—1 (—1^(—0)Кэ(г*—0) — —0 /1 (^1 )К0 (г*—1))

22 —1 — —0

г* ^ 1.

(20)

Рассчитав потенциал смещений Ф; и функцию тока ф( через функции посредством формул (10) и (11), вычислим остаточные смещения частиц воды в радиальном направлении и скорость вихревого те-

чения

5Ф,-

п. =_-

^гоо — о »

дг

У-ОО -

дфи

дг

(21)

В рамках принятого приближения (= ) в верхнем слое эти величины соответствуют случаю однородного океана

А«>(г*) =

П0 шах#

Н1 + Н2

х(г*),

/ п0 шахл / ! с,

У1((г *) = - 1 х(г *),

Н1 + Н2

(22) (23)

х(г*) = | /1 (—0Г*)К (—0), 0 < г * <1, (24)

\/1 (М0)К1(М0Г*), г* ^ 1.

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

П2((Г*)= П0 ШаХ

л

Н1 + Н2

[х(г*) + Дх(г*)],

У2((Г *) = -

П0 шахЛ/

Н1 + Н2

[х(г*) + Дх(г*)]

Дх(г *) =

5р Мг*) + в^2(г*))

р

(1 + в)

(25)

(26) (27)

( *) —2] М—1Г*)КК—1), 0 < г* <1, °"1(г ) = —о <

—2 1/1(м1)К1(^1Г*), г * ^ 1.

^2(Г* ) = -2^-2 — 1 —0

/1 (-0Г*)К1 (-0) - /1 (—1Г*)К1 (—1),

0 < г* <1, /1 (—0)К1 (—0Г*) - /1 (—1 )К1 (—1Г*), г* > 1.

2. Результаты и их обсуждение

Полученное в предыдущем разделе решение, описывающее остаточные поля в двухслойном вращающемся океане, зависит от трех безразмерных параметров — в, — 0 и —1. Эти параметры удобны для получения аналитического решения, но они не являются независимыми. В частности, существует следующая связь: —2 = —0(1 + в)2/(в^р/р). Для представления результатов удобнее использовать набор независимых параметров: в, —0 и 6р/р.

Свойственное реальному океану соотношение толщин слоев (Н1 ~ 102 м, Н2 ~ 103 м) определяет типичное значение параметра в ~ 10. Впрочем, величина в может меняться в очень широком диапазоне от вШ1п ~ 10-3 (тонкий придонный перемешанный слой) до вшах ~ 103 (тонкий верхний перемешанный слой).

В условиях нашей планеты величина —0 варьируются от — шш « 0 в экваториальной зоне до —шах ~ 1 (высокие широты, протяженный очаг, шельфовые глубины). Типичное значение этого параметра при / ~ 10-4 с-1, Л ~ 105 м, Н1 + Н2 ~ 103 м составляет —0 ~ 10-1.

Во всех расчетах, представленных в этом разделе, мы будем основываться на фиксированной величине 5р/р = 0.003 — это типичное значение относительного перепада плотности для океана [12]. Два оставшихся свободных параметра будем варьировать в пределах, перекрывающих всевозможные реальные значения: 0.001 < —0< 10, 10-3< в < 103.

На рис. 2 показана форма смещения свободной поверхности воды в остаточном геострофическом вихре £1((г*) и форма смещения поверхности раздела слоев £2((г*), которые рассчитаны по формулам (16) и (18) при типичных значениях параметров —0 = 0.1 и в = 10. В первую очередь обращает на себя внимание тот факт, что смещение поверхности раздела слоев по амплитуде на порядок превосходит смещение свободной поверхности, т. е. наблюдается эффект связывания возмущения устойчивой стратификации в геострофическом вихре.

^ 1со /п Отах ^гоэ/Л Отах

Сг/а + Р)

10" г!Я

Рис. 2. Смещение свободной поверхности воды и смещение поверхности раздела слоев в остаточном геострофическом вихре. Смещение поверхности раздела слоев, рассчитанное в приближении «твердой крышки» ((1/(1 + в) )• Расчет выполнен по формулам (16) и (18) при —0 = 0.1, в = 10, 5р/р = 0.003

Пунктирной линией на рис. 2 показана форма смещения поверхности раздела слоев, рассчитанная в приближении «твердой крышки» (т.е. при (2(г*) = 0). Видно, что приближение «твердой крышки» неплохо работает внутри области источника (г* <1), но вне источника оно дает принципиально иной результат.

Смещение свободной поверхности и смещение поверхности раздела во всех случаях достигают максимальных значений в центре источника. Полагая в формулах (16) и (18) г* = 0, можно рассчитать эти максимальные значения в зависимости от параметров —0 и в. Результаты расчетов показаны на рис. 3. Видно, что во всех случаях амплитуда смещения поверхности раздела

Рис. 3. Амплитуда смещения свободной поверхности воды тах и смещения поверхности раздела слоев тах в остаточном геострофическом вихре как функция параметра р0. Амплитуда смещения поверхности раздела слоев, рассчитанная в приближении «твердой крышки» (£1 тах/(1 + в) )• Расчет выполнен по формулам (16) и (18) при 5р/р = 0.003 и различных значениях параметра в (указаны на рисунке)

превосходит амплитуду смещения свободной поверхности. А при малых значениях параметра в (тонкий нижний слой) амплитуда смещения поверхности раздела может превосходить амплитуду смещения свободной поверхности на несколько порядков.

Пунктирными линиями на рис. 3 показаны зависимости, рассчитанные в приближении «твердой крышки». Неработоспособность этого приближения, особенно ярко проявляется при больших значениях параметров и в.

На рис. 4 представлен вид функций х(г*) и Дх(г*), которые определяют пространственную структуру остаточных полей — остаточного смещения частиц воды Ао(г*) и скорости течения в геострофическом вихре ^¡ж(г*). Расчет выполнен по формулам (24) и (27) при типичных значениях параметров р0 = 0.1 и в = 10. При всех значениях аргумента х(г*) >0 и Дх(г*) >0. То есть в соответствии с формулами (22) и (25) поднятие дна (п0тах >0) сопровождается смещением частиц воды в положительном направлении — от центра

зс.Д% 10°

10"1

10~2

10"3

10"4

10"5

Ю-1 10° 101 102 г/Я

Рис. 4. Функции х + Дх, X и Дх, определяющие пространственную структуру остаточных полей. Расчет выполнен по формулам (24) и (27) при р0 = 0.1, в = 10, 5р/р = 0.003

источника. Отрицательный знак в формулах (23) и (26) означает, что в Северном полушарии, где параметр Кориолиса [ >0, поднятие дна вызывает вихрь, вращающийся в отрицательном направлении (антициклонический). Опускание дна (п0тах < 0) приводит к обратному результату — смещению частиц к центру и к циклоническому вихрю. Хорошо видно, что добавка к полю скорости в нижнем слое Дх(г*), которая обязана влиянию стратификации, незначительно увеличивает остаточное смещение частиц воды и скорость вихревого течения вблизи границы источника (г* = 1).

Максимальные значения функций х(г*) и Дх(г*) всегда достигаются в точке г* = 1. Зная это, можно рассчитать максимальные значения функций х(г*) и Дх(г*) в зависимости от параметров р0 и в. Результаты расчетов показаны на рис. 5. Из рисунка видно, что при типичном значении параметра р0 = 0.1 и при больших значениях параметра в (в > 10) стратификация незначительно влияет на амплитуду остаточных смещений частиц и скорость вихревого движения в нижнем слое. Но в тех случаях, когда в < 10, пренебречь этим влиянием становится невозможно. Более того, для относительно тонких нижних слоев в < 0.1 вклад стратификации в остаточные поля может оказаться преобладающим.

л тах' АЗС тах

- ЗСтах

тах

р= 10 -.....

"р=10"-.....

р= 10°......

р= 101

Р=Ю2 >....;

р=10 1 1 ..........

10"3 10"2 10"1 10° |Л0

Рис. 5. Максимальные значения функций х и Дх в зависимости от параметра р0. Расчет выполнен по формулам (24) и (27) при 5р/р = 0.003 и различных значениях параметра в (указаны на рисунке)

Заключение

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

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

Работа выполнена при финансовой поддержке РФФИ (гранты 13-05-00337, 13-05-92100 и 12-0531422).

Список литературы

1. Bolshakova A.V., Nosov M.A. // Pure and Applied Geo-phys. 2011. 168. P. 2023.

2. Носов М.А., Колесов С.В., Левин Б.В. // Докл. РАН. 2011. 441, № 1. C. 108.

3. Пелиновский Е.Н. Гидродинамика волн цунами. Н.Новгород, 1996.

4. Levin B.W., Nosov M.A. Physics of Tsunamis. Springer, 2008.

5. Носов М. А., Мошенцева А. В., Левин Б.В. // Докл. РАН. 2011. 438, № 5. C. 694.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

6. Ингель Л.Х. // Доклады РАН. 1998. 362, № 4. C. 548.

7. Доценко С.Ф. // Изв. РАН. Физика атмосферы и океана. 1999. 35, № 5. С. 706.

8. Доценко С.Ф. // Прикладная гидромеханика. 2000. 2(74), № 4. С. 24.

9. Доценко С.Ф., Шокин Ю.И. // Вычисл. технологии. 2001. 6, № 1. C. 13.

10. Nosov M.A., Moshenceva A.V., Kolesov S.V. // Pure and Appl. Geophys. 2013. 170, N 9, 10. P. 1647. http://dx.doi.org/10.1007/s00024-012-0605-2.

11. Носов М.А., Нурисламова Г.Н. // Вестн. Моск. ун-та. Физ. Астрон. 2012. № 5. C. 44.

12. GillA.E. Atmosphere-ocean dynamics. N.Y.; L.; P., 1982.

13. Hammack J.L. // J. Phys. Oceanogr. 1980. 10, N 9. P. 1455.

14. Левин Б.В., Носов М.А., Павлов В.П., Рыкунов Л.Н. // Докл. РАН. 1998. 358, № 3. С. 399.

15. Nosov M.A. // Volcanology and Seismology. 1998. 19, N 3. P. 371.

16. Полянин А.Д. Справочник по линейным уравнениям математической физики. М., 2001.

Traces of a tsunamigenic earthquake in the rotating stratified ocean M.A. Nosova, G.N. Nurislamovab

Department of Marine and Inland Water Physics, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.

E-mail: anosov@phys.msu.ru, bgulnaz1205@yandeх.ru.

Equations for description of potential and vortex residual fields that accompany tsunami generation by co-seismic bottom deformation in the rotating two-layer stratified ocean are derived within the framework of linear long-wave theory. An approximate analytical solution to the respective axisymmetric problem is obtained. On the base of the analytical solution, under typical conditions for a tsunami source, particular features of the residual fields that arise due to stratification are examined.

Keywords: linear long-wave theory, tsunami generation, co-seismic bottom deformation, residual hydrodynamic fields, geostrophic vortex, two-layer stratification, Rossby radius of deformation. PACS: 91.30.Nw, 92.10.Iv. Received 2 July 2013.

English version: Moscow University Physics Bulletin 6(2013). Сведения об авторах

1. Носов Михаил Александрович — докт. физ.-мат. наук, профессор; тел.: (495) 939-36-98, e-mail: nosov@phys.msu.ru.

2. Нурисламова Гульназ Нуровна —студентка; тел.: (495) 939-36-98, e-mail: gulnaz1205@yandex.ru.

i Надоели баннеры? Вы всегда можете отключить рекламу.