Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
синоптика.pdf
Скачиваний:
1055
Добавлен:
13.03.2016
Размер:
13.13 Mб
Скачать

9. Основные характеристики полей метеорологических величин

81

iС изменениями вихря скорости тесно связаны изменения барического

поля во времени

iЗнаки изменений вихря и давления (геопотенциала) противоположны

9.9.3. Уравнение тенденции вихря скорости

Уравнения движения позволяют получить формулы для определения изменения

∂Ω

вихря скорости со временем t в данном пункте, что может быть использовано для

P

оценки изменения давления t .

Для свободной атмосферы имеем

du

 

= −g

H

+ lv

 

 

 

x

 

 

dt

 

 

 

 

 

 

 

 

 

dv

 

= −g H − lu

 

 

 

 

y

 

dt

 

 

 

Продифференцируем первое уравнение по Y, второе – по X и найдём разность между вторым и первым. В результате получим:

∂Ω

∂Ω

∂Ω

∂Ω

∂l

 

∂l

 

∂τ v

∂τ ∂u

(9.9.7)

 

 

 

 

 

 

 

 

 

 

 

 

t = −(u x + v y + τ p ) (Ω + l)D (u

x + v

y ) (

x p

y p )

 

Данное уравнение есть уравнение тенденции вихря скорости. Согласно уравне-

∂Ω

нию тенденции вихря скорости, локальное изменение завихренности t определяется

действием ряда факторов.

 

 

 

 

Первое слагаемое в уравнении

 

 

 

 

∂Ω

∂Ω

∂Ω

∂Ω

(9.9.8)

t 1 = −(u

x + v

y + τ

p )

 

∂Ω

характеризует зависимость t от адвекции вихря скорости. Его называют вихревой

составляющей уравнения тенденции вихря скорости. Второе слагаемое

∂Ω

(9.9.9)

t 2 = −(Ω + l)D

 

Н.А. Дашко Курс лекций по синоптической метеорологии

9. Основные характеристики полей метеорологических величин

82

∂Ω

характеризует зависимость t от горизонтальной дивергенции скорости, его называ-

ют дивергентной составляющей уравнения тенденции вихря скорости. Что касается третьего слагаемого,

∂Ω

∂l

∂l

(9.9.10)

 

 

 

 

 

t 3 = −(u

x + v

y ) ,

 

∂l

то если направить ось X вдоль параллели, тогда x =0 и, следовательно,

∂Ω

= −v

∂l

(9.9.11)

 

 

t 3

y

 

∂Ω

данное слагаемое будет характеризовать зависимость t от меридиональных смеще-

ний воздушных масс – меридиональная составляющая уравнения тенденции вихря скорости.

Четвёртое слагаемое

 

 

 

 

 

∂Ω

∂τ v

∂τ ∂u

(9.9.12)

t 4 = −(

 

 

 

y p )

x p

 

∂Ω

показывает зависимость t от горизонтального градиента вертикальной скорости и

изменения ветра с высотой. Таким образом:

∂Ω

= (

∂Ω)вихр + (

∂Ω)див + (

∂Ω)мерид + (

∂Ω)верт .

(9.9.13)

t

 

t

t

t

t

 

Произведём оценку слагаемых в уравнении.

Характерные синоптические масштабы:

Скорость ветра U – 10 м/с, характерный масштаб изменения скорости имеет тот же порядок, что и скорость,

Время T – сутки, 105 с,

Длина L – 1000 км, 106 м, Высота Z – 1-10 км, 103 -104 м,

 

v

 

u

 

 

U

 

10

 

 

Завихренность []=

 

=

=

= 105 ,

(c1 ) ,

x

y

L

106

 

 

 

 

 

 

 

 

∂Ω

=

U

=

10

= 10

10

,

(c2 ) b,

 

 

 

 

 

LT

106

105

 

t

 

 

 

 

 

 

Н.А. Дашко Курс лекций по синоптической метеорологии

9. Основные характеристики полей метеорологических величин

 

 

83

 

∂Ω

 

U

 

 

10

 

10

, (c2 ) ,…

 

∂Ω

= 1010 , (c2 ) ,

u

 

= U

 

= 10

×

 

= 10

v

 

LL

106106

 

x

 

 

 

 

 

 

 

x

 

[P]=

[

 

]

 

 

 

× 104 = 105 , (Па)

(одна атмосфера – 105

 

− ρgz = −1 × 10

 

 

 

 

 

 

 

 

 

 

P

 

105

 

 

 

Вертикальная составляющая скорости [τ]=

 

=

 

= 1 ÷ 10

1

,

105

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

Перенос вихря по вертикали

 

 

 

 

 

 

 

 

∂Ω

 

ρgz

 

∂Ω

∂Ω

= 1010 ÷ 10

11 ,

 

(c2 ) ,

 

 

τ

p

=

 

×

 

 

= t

 

 

 

t

ρgz

 

 

 

Па),

Па

( с ) ,

 

Параметр Кориолиса [l]= [2ω sin ϕ]= 104 ,

(c1 ) ,

 

[β]=

 

∂l

 

2ω cos

ϕ

=

104

= 1011 , (м1c1 ) ,

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

R 3

 

107

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂τ

 

 

 

 

6

,

 

(c1 ) (из уравнения неразрывности),

 

Дивергенция [D]=

p

= 10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[lD]= 104106 = 1010 , (c2 ), [D]= 105106 = 1011 , (c2 ),

 

 

∂τ

v

 

 

∂τ ∂u

 

1

Па

 

 

 

1

м

 

 

 

 

 

 

 

 

,

 

 

 

 

с

 

= 1011 ÷ 1012 ,

2 )

 

 

 

 

=

 

 

с

 

×

 

 

 

 

 

 

 

y p

 

 

 

 

 

 

 

 

 

 

 

x p

 

 

10

6 м

 

105 Па

 

 

 

 

 

 

 

∂l

 

 

 

 

∂l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

, v

 

 

 

= 10 × 1011 = 1010 , (c2 ) .

 

 

 

x

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, подставляя полученные оценки в уравнение тенденции вихря скорости,

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[1010 ][1010 , 1010 ,

 

1011 ],

[1010 ,

1011 ],

[1010 , 1010 ], [1012 , 1012 ]

 

Оценка слагаемых в уравнении тенденции вихря скорости показывает, что вели-

 

∂Ω

 

 

 

 

 

∂Ω

 

 

 

∂Ω

 

 

 

 

 

 

 

∂l

 

 

∂l

 

 

чины

t

,

 

(u

 

x

+ v

y ),

lD, u

 

,

v

 

имеют одинаковый порядок и несколько

 

 

x

y

больше остальных членов уравнения.

∂Ω

В первом слагаемом переносом вихря по вертикали τ p можно пренебречь,

тогда

Н.А. Дашко Курс лекций по синоптической метеорологии

9. Основные характеристики полей метеорологических величин

84

∂Ω

= −(u

∂Ω

+ v

∂Ω

 

t 1

x

y ) .

 

Это слагаемое больше других и вносит наибольший вклад в изменение вихря со временем.

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

∂Ω = −lD . t 2

Последним слагаемым часто пренебрегают. Однако, вклад его нельзя недооценивать в горных районах.

Окончательно имеем:

∂Ω

 

∂Ω

 

∂Ω

∂l

t

= −(u

x

+ v

y ) − lD v

 

.

y

Отметим, что в исходных уравнениях движения

(9.9.14)

u v

x , y примерно в 10 раз

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

9.9.4. Уравнение тенденции вихря скорости

внатуральной системе координат

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

Поскольку имеет место соотношение

Ω = gl 2H ,

то можем записать:

Н.А. Дашко Курс лекций по синоптической метеорологии

9. Основные характеристики полей метеорологических величин

85

 

∂Ω

=

g

 

 

( 2H)

 

t

l ∂t

 

 

 

или

 

g

2 H

 

∂Ω

=

(9.9.15)

t

l

 

 

 

 

t

 

Ранее мы качественно показали, что изменения вихря скорости связаны с изменениями барического поля во времени:

Для циклонической завихренности:

Ω > 0, 2P > 0, P < 0

Для антициклонической завихренности:

Ω < 0, 2P < 0, P > 0 .

Для установления данной связи рассмотрим поле геопотенциала.

Поскольку поле геопотенциала имеет волновую структуру, то для каждой элементарной волны в первом приближении имеем:

H = A sin αx × cosβy,

2H = 2H + 2H = −[Aα2 sin αx × Aβ2 cosβy] = ∂x2 y2

= −A(α2 + β2 )(sin αx × cosβy) = −mH.

Здесь A – амплитуда волны. Тогда

2H = −mH ,

t 2H = 2 Ht = −m Ht .

∂Ω

=

g

2 H

= −

g

m

H

,

(9.9.16)

t

l

 

t

 

t

 

l

 

 

Рассматривая вихревую составляющую уравнения тенденции вихря скорости,

∂Ω

= −(u

∂Ω

+ v

∂Ω

t вихр

x

y ) ,

где вместо u, v и подставим выражения в геострофическом приближении: ug = − gl Hy , vg = gl Hx , Ω = gl 2H

то получим:

Н.А. Дашко Курс лекций по синоптической метеорологии

9. Основные характеристики полей метеорологических величин

∂Ω

= −(

 

g2

H

 

 

2H +

t

 

l2

 

y x

вихр

 

 

 

 

 

 

или через оператор Якоби

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Ω

 

 

= −

 

g2

 

 

(H, 2H)

 

t вихр

 

 

l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Привлекая выражение (9.9.16), получим

 

 

g

 

 

 

H

 

 

 

 

g2

 

 

m

 

 

 

 

= −

 

(H,

 

l

t

l2

 

m

H

 

=

g

 

(H, 2H),

 

t

l

 

 

 

 

 

 

 

 

 

 

 

 

H

=

g

 

(H, 2H).

 

t

lm

 

 

 

 

 

 

 

 

 

 

 

 

 

g2

H

2H)

l2

x y

 

.

2H),

86

(9.9.17)

(9.9.18)

Осуществим переход к натуральным координатам, где ось X (OS) направлена по касательной к изогипсе, ось n направлена в сторону возрастающих значений H, т.е. совпадает с осью Y, но противоположно направлена.

Тогда

H

=

g

(H, 2H) =

 

g

[

H

(

2H

+

2H

)

H

(

2H

+

2H

)],

 

 

 

 

 

 

 

 

 

 

x2

 

y2

 

 

 

 

2

y2

t

lm

lm

x

 

y

 

 

y x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

=

 

H

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

=

 

 

g

[

H

(

2 H

+

2 H

)].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

lm

y x

x2

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В случае прямолинейных изогипс вторая производная

2H

,

характеризующая

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

= tgα,

 

α = 0 ), равна нулю.

изменение тангенса угла наклона изогипсы H к оси S ( x

 

2H

Вслучае криволинейных изогипс x2 0.

Для её определения рассмотрим кривизну явно заданной кривой y=f(x), для которой в дифференциальной геометрии получена формула

Н.А. Дашко Курс лекций по синоптической метеорологии

9. Основные характеристики полей метеорологических величин

 

 

 

 

 

87

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

κ =

 

 

 

 

 

dx2

 

 

 

 

 

 

 

,

 

 

 

 

 

 

(9.9.19)

 

 

 

 

 

 

 

 

 

 

dy

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

) 2 ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[1 + (

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2 y

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выразим производные

 

 

 

 

,

 

 

 

. Поскольку H(x, y)=const – неявно заданная

dx2

 

 

dx

кривая, т.е. мы не знаем вида функции y=f(x), найдём

 

 

dH(x,y)

и (при условии, что

 

 

 

dx

 

H(x, y)=const), т.е.

dH(x,y)

 

0. По правилу дифференцирования сложной функции

dx

 

 

dH(X,Y)

 

 

 

 

H(X,Y)

 

dx

 

 

 

H( X,Y)

 

dy

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

+

 

 

 

0

,

 

 

 

dx

 

 

 

 

 

x

 

dx

 

y

dx

 

отсюда,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

H(x,y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

x

.

 

 

 

 

 

 

(9.9.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

H(x,y)

 

 

 

 

 

 

 

 

 

 

y

Кроме того, необходимо найти вторую производную, повторно продифференцировав полученное выражение (9.9.20) с учётом того, что все частные производные в правой части зависят и от x и от y, при этом y=f(x):

d2y

 

H(x,y)

[

2H(x,y)

+

2H(x,y) dy

] +

H(x,y)

[

2H(x,y) 2H(x,y)

+

2H(x,y) dy

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

y

x2

xy dx

 

x

yx

x2

y2

dx

. (9.9.21)

dx2

 

 

 

 

 

 

 

(

H

(x,y)

)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

dy

В натуральной системе координат выражение (9.9.20) dx =0. С учётом этого для

(9.9.20) получим:

 

 

H

(x, y)

 

2H

(x, y)

d 2 y

 

 

×

 

 

= −

 

y

 

x

2

 

 

 

dx 2

 

(

H

(x, y)

)2

 

 

 

 

 

 

y

 

 

 

 

2H(x, y) 2H(x, y)

 

 

x2

 

 

x2

 

 

(9.9.22)

= −

H(x, y)

=

H(x, y)

.

 

(

 

)

 

(

 

)

 

 

y

n

 

dy

Подставляя выражения из (9.9.22) и (9.9.20) в (9.9.19), учитывая, что dx =0, по-

лучим

Н.А. Дашко Курс лекций по синоптической метеорологии

9. Основные характеристики полей метеорологических величин

88

 

2H

 

 

2H

 

Hss

 

 

dx2

 

 

 

ds2

 

,

κ=

 

 

 

=

 

 

 

=

 

 

 

H

 

 

H

Hn

 

(

 

 

)

 

(

 

)

 

 

 

 

 

 

n

n

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2H

=

2H

=

κ

H

= κHn ,

dx2

ds2

 

n

 

 

 

 

 

 

 

 

 

Тогда для вихревой составляющей в натуральной системе координат имеем:

 

H

 

 

g H

H 2H

 

 

 

 

 

=

 

n [

 

(κ n ) +

 

( n 2

)],

 

 

t

lm

s

s

 

H

 

g H ∂κ ∂H H 2H

 

 

=

 

n [ s n

+ κ

 

( n ) +

 

(

n 2

)],

t

lm

s

s

или, следуя общепринятым обозначениям для натуральной системы, получим уравне-

ние Бугаева:

H

=

g

Hn (κs Hn + κHns + Hnns )

(9.9.23)

t

lm

 

 

 

Анализ данного уравнения показывает, что при положительной адвекции вихря скорости

∂Ω

>0 –

H

<0,

t

t

при увеличении антициклонического вихря, т.е. при отрицательной адвекции вихря:

∂Ω

<0 –

H

>0.

t

t

Hn >0, поскольку n направлена в сторону возрастающих значений H.

Hns – изменение Hn вдоль потока S: Hns >0 при сходимости изогипс, Hns <0

при расходимости изогипс.

κ– кривизна изогипс: κ >0 при циклонической кривизне изогипс, κ <0 при антициклонической кривизне изогипс.

κs – изменение кривизны по потоку: κs >0 – при увеличении циклонической

кривизны изогипс или уменьшении антициклонической; κs <0 – при уменьшении ци-

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

Hnns – изменение градиента геопотенциала вдоль потока и и по нормали к изо-

гипсам (вклад данного слагаемого невелик, обычно не рассматривается).

Н.А. Дашко Курс лекций по синоптической метеорологии

9. Основные характеристики полей метеорологических величин

89

H

Поскольку Hn >0, знак t определяется знаком слагаемых в скобках.

Первое слагаемое κs Hn >0:

При κs >0 – при увеличении циклонической кривизны изогипс или уменьшении антициклонической, κs Hn <0

При κs <0 – при уменьшении циклонической кривизны изогипс или увеличении

антициклонической,

Второе слагаемое:

κHns >0

Если κ >0, т.е. при циклонической кривизне изогипс, Hns >0 при сходимости изогипс;

Если κ <0, т.е. при антициклонической кривизне изогипс; Hns <0 при расходи-

мости изогипс.

κHns <0

Если κ <0, т.е. при антициклонической кривизне изогипс; Hns >0 при сходимо-

сти изогипс;

Если κ >0, т.е. при циклонической кривизне изогипс, Hns <0 при расходимости изогипс.

iНаибольший вклад в изменение знака

H

вносит слагаемое κs Hn

 

t

 

iЕсли знаки противоположны, то учитывается именно знак κs Hn

iВторое слагаемое κHns вносит усиливающий или ослабляющий эффект в

H

изменение знака t .

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

Практический опыт синоптика показывает, что наиболее благоприятные условия для развития циклона складываются, когда приземный центр его располагается под передней частью высотной барической ложбины на АТ500, где, при наличии значительных

Н.А. Дашко Курс лекций по синоптической метеорологии