Индексы неустойчивости. Методы прогноза конвективной облачности и явлений.

Метеоцентр.Азия - наш новый сайт с высокодетализированными прогнозами погоды по пунктам России и мира
Облегчённая версия Метеоклуба (для КПК)
Карта активных участников Метеоклуба (105 кБ)    Таблица дней рождения активных участников Метеоклуба
Клуб любителей метеорологии (группа ВКонтакте)

Индексы неустойчивости. Методы прогноза конвективной облачности и явлений.

Сейчас в Метеоклубе:
Участников - 4 [ Balin, LESS, roman ashirov, Mark2 ]
Максимальное одновременное количество посетителей: 4 [29 Июн 2017 15:32]
Гостей - 0 / Участников - 4

 - Начало - Ответить - Статистика - Регистрация - Поиск -
МЕТЕОКЛУБ : независимое сообщество любителей метеорологии (Европа и Азия) : ФОРУМ О ПОГОДЕ И ПРИРОДЕ / Метеорология: наука и практика / Индексы неустойчивости. Методы прогноза конвективной облачности и явлений.
<< . 1 . 2 . 3 . 4 . 5 . 6 . 7 . 8 . 9 . 10 ... 30 . 31 . >>
Автор Сообщение
Cumulonimbus incus
Участник
Письмо
Кишинёв, Молдова
# Дата: 18 Июн 2012 15:17


Помогите с формулой для индекса Фауста! Весь форум излазил, все руководства в Интернете посмотрел, а там всё таблицы да графики. А для ГрАДСа нужна формула.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 18 Июн 2012 21:09


Cumulonimbus incus

Формулы как таковой нет. При дефиците ноль градусов кривые на графике - влажные адиабаты. При дефиците больше нуля поправка пропорциональная дефициту.

Cumulonimbus incus
Участник
Письмо
Кишинёв, Молдова
# Дата: 18 Июн 2012 21:49


Corvus

Я понял, но сразу же возникает вопрос: как эти товарищи его считают?
http://portal.chmi.cz/files/portal/docs/meteo/oa/p tu_grafy.html

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 18 Июн 2012 23:09


Cumulonimbus incus

Считать можно, просто скрипт будет громоздкий.

Например, по слоям толщиной в 50 гПа посчитать влажноад. градиент и так от 850 подняться до 500 гПа. Потом ввести поправку на средний дефицит в слое 850-500 гПа.

TornadoF5
Участник
Письмо
Харьков, Украина. (Игорь)
# Дата: 19 Июн 2012 20:28


Интересная ситуация с индексами вышла. По прогнозу завтра в Одесской и Николаевской обл. К-инекс будет зашкаливать (около 37-40), при том, что в районе ХФ на Слобожанщине он какой-то вялый (около 30). Основная неустойчивость и осадки ожидаются на северо-востоке страны (в зоне вронта), зато К-индекс что-то вообще непоказательный.
Вопрос: с чем может быть связаны такие значения К-индекса, ведь гроз в тех местах завтра не будет? Обычно он хорошо коррелирует с другими параметрами и правильно отражает грозовые очаги, а сейчас всё наоборот. Ниже привожу прогност. карты К-индекса, ЛИ и АД на завтра по Вознесенску (Николаевская обл.) на срок 15z:



Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 19 Июн 2012 20:33


TornadoF5

Я сделал вывод (и реализовал на сайте Азия), что вероятность грозы, определённую по К, нужно обязательно корректировать по КФ (индексу Фатеева).

Как-то так:

if( $v['KF'] >= 20 ) $v['PK'] = $v['PK'] + 20;
if( $v['KF'] >= 10 && $v['KF'] < 20 ) $v['PK'] = $v['PK'] + 10;
else if( $v['KF'] < 0 ) $v['PK'] = $v['PK'] - 20;

_Ilya_
Участник
Письмо
г. Белгород
# Дата: 19 Июн 2012 20:34


TornadoF5

Да, смотрел этот индекс. Если ему верить, завтра будет дежа-вю вчерашнего, с позволения сказать, фронта.

TornadoF5
Участник
Письмо
Харьков, Украина. (Игорь)
# Дата: 19 Июн 2012 20:36


_Ilya_
Если ему верить, завтра будет дежа-вю вчерашнего

Не надо нам такого опять! :)

TornadoF5
Участник
Письмо
Харьков, Украина. (Игорь)
# Дата: 19 Июн 2012 20:38


Corvus
Как-то так:

if( $v['KF'] >= 20 ) $v['PK'] = $v['PK'] + 20;
if( $v['KF'] >= 10 && $v['KF'] < 20 ) $v['PK'] = $v['PK'] + 10;
else if( $v['KF'] < 0 ) $v['PK'] = $v['PK'] - 20;


Что это?

Sunspot
Участник
Письмо
Ульяновск
# Дата: 19 Июн 2012 21:16


Что это?
TornadoF5

Если вероятность грозы, вычисленная по индексу Фатеева больше или равна 20, тогда вероятность по К=К+20, если по KF в интервале 10..20, тогда по К=К+10, иначе если KF<0 то К=К-20. Как то так, наверное. Извините, если не так понял, Владимир сам лучше объяснит.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 19 Июн 2012 21:22


Как то так, наверное.
Sunspot

Да, почти так :)

Если индекс Фатеева больше или равен 20, тогда вероятность грозы по К увеличиваем на 20%, если индекс Фатеева в интервале 10..20, увеличиваем на 10%, иначе если индекс Фатеева<0 то вероятность грозы уменьшаем на 20%.

TornadoF5
Участник
Письмо
Харьков, Украина. (Игорь)
# Дата: 19 Июн 2012 21:28


Sunspot
Corvus
Спасибо за разъяснения! :)

Cumulonimbus incus
Участник
Письмо
Кишинёв, Молдова
# Дата: 19 Июн 2012 21:32 - Поправил: Cumulonimbus incus


TornadoF5

Есть ещё один показатель - средний дефицит точки росы в слое 850-500 гПа. По Фаусту, если он превышает 14 - гроз не будет. В "Руководстве по краткосрочным прогнозам погоды" 1986 г. вообще указаны значения 8...10 в качестве предельных.
Завтра он в указанных Вами областях будет равен 14...22, и лишь вечером местами упадёт до 12. Так что я бы дал по этому региону слабые ливневые осадки с возможными спорадическими грозовыми разрядами, учитывая, что индекс Пескова будет в районе 0...2.

spralex
Участник
Письмо
г. Конотоп, Сумская обл., Украина
# Дата: 20 Июн 2012 23:31


На известном ресурсе (модель WRF, тестируемая УкрНИГМИ) появились индексы MCAPE и MCIN.
Пример на субботу:


Cumulonimbus incus
Участник
Письмо
Кишинёв, Молдова
# Дата: 21 Июн 2012 21:51


Corvus

Да, скрипт будет сложный. Представляю, сколько времени ГрАДС будет его считать...

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 21 Июн 2012 22:25


Cumulonimbus incus

Попытка не пытка. Главное, грамотно написать цикл, который будет считать Т по слоям с учётом влажноадиабатического градиента.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 8 Июл 2012 12:41


Можно рассчитать F (индекс Фауста) приближённо.

1. Считаем ЭПТ850 (эквивалентно-потенциальная температура) по Т850, Тд850 и давлению 850 гПа.
2. Считаем ЭПТ500 по Т500, Тд500 и давлению 500 гПа.
3. Считаем псевдоЛИ = ЭПТ850 - ЭПТ500.
4. Рассчитываем поправку к псевдоЛИ в зависимости от среднего дефицита в слое 850-500 гПа:

if($d<6) $r = 0;
else if($d==6) $r = 2;
else if($d==7 || $d==8) $r = 3;
else if($d==9 || $d==10) $r = 4;
else if($d==11 || $d==12 || $d==13) $r = 5;
else if($d==14 || $d==15) $r = 6;
else if($d>=16) $r = 9;

5. Эту поправку прибавляем к псевдоЛИ.

Так я сделал расчёт F в модели НАЕФС.

bullterrier
Участник
Письмо
Пермь
# Дата: 8 Июл 2012 14:27


Corvus
Подскажите пожалуйста, верна ли эта формула для подсчета Lifted Index?

Ew = 6.1078 * exp([(9.5939 * Tw) - 307.004]/[(0.556 * Tw) + 219.522])
E = Ew - 0.35 * (Td - Tw)
Tp = -1 * {[ln(E/6.1078) * 219.522] + 307.004} / {[ln(E/6.1078) * 0.556] - 9.59539}
LI= Tc(500mb) - Tp(500mb)

Tc(500mb) - это температура на уровне 500 гПа.
Tp(500mb) - это температура частицы, поднятой от земли до уровня 500 гПа.
Td - точка росы.
Tw - приземная температура (?)

Просто хотелось бы сделать его расчет по данным модели WRF, которую гоняют у нас в ПГНИУ.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 8 Июл 2012 20:39


bullterrier

А можно ссылку на первоисточник?

Похоже, Tw - температура смоченного термометра, но не уверен.

bullterrier
Участник
Письмо
Пермь
# Дата: 8 Июл 2012 21:35 - Поправил: bullterrier


А можно ссылку на первоисточник?

Corvus
URL

Мне LI нужен для подсчета других индексов и прогноза диаметра града по методу Заводченкова. Кстати, этот метод очень себя хорошо показал по модели GFS, на картах все очень наглядно получается.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 9 Июл 2012 07:20


Мне LI нужен для подсчета других индексов и прогноза диаметра града по методу Заводченкова
bullterrier

Я понимаю, он много для чего нужен, включая индекс Пескова.

Жаль, что на настоящий первоисточник ссылки нет (а только на наш же форум), это бы многое прояснило. У меня опасение, что эти формулы не совсем те.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 9 Июл 2012 09:09 - Поправил: Corvus


Ну вот, как я и опасался, это смесь бульдога с носорогом.

Вот это расчёт точки росы по психрометру:

'Saturation Vapor Pressure Wet
Ew = 6.1078 * exp([(9.5939 * Tw) - 307.004]/[(0.556 * Tw) + 219.522])
'Saturation Vapor Pressure Dry
Es = 6.1078 * exp([(9.5939 * Td) - 307.004]/[(0.556 * Td) + 219.522])

E = Ew - 0.35 * (Td - Tw) 'Actual Vapor Pressure

Relative Humidity
RH = (E / Es) * 100

Dew Point
Tp = -1 * {[ln(E/6.1078) * 219.522] + 307.004} / {[ln(E/6.1078) * 0.556] - 9.59539}


http://www.aprweather.com/pages/calc.htm

Где Tw - Т смоченного термометра, Td - Т сухого термометра, Tp - точка росы.

А в формуле для ЛИ Тр - это Т поднимающейся частицы и вычисляется она совсем иначе.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 9 Июл 2012 09:21


Интересный индекс попался

DCI = T(850 mb) + Td(850 mb) - LI(sfc-500 mb)

http://www.crh.noaa.gov/lmk/soo/docu/indices.php

При его значении более 30 высока вероятность сильных гроз и ливней.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 9 Июл 2012 09:30


Расчёт температуры смоченного термометра (нужен в некоторых методах прогноза шквалов и других конвективных явлений):

Wet Bulb Temperature


Variables:
Tc = Temperature in Degrees C
RH = Relative Humidity in form 88 (not 0.88)

Optional Variable (for more accuracy):
P = Pressure
or
Constant (with up to 0.2° inaccuracy):
P = 985


Tdc = ((Tc - (14.55 + 0.114 * Tc) * (1 - (0.01 * RH)) - ((2.5 + 0.007 * Tc) *
(1 - (0.01 * RH))) ^ 3 - (15.9 + 0.117 * Tc) * (1 - (0.01 * RH)) ^ 14))

E = (6.11 * 10 ^ (7.5 * Tdc / (237.7 + Tdc)))

WBc = (((0.00066 * P) * Tc) + ((4098 * E) / ((Tdc + 237.7) ^ 2) * Tdc)) / ((0.00066 * P) + (4098 * E) / ((Tdc + 237.7) ^ 2))



http://www.aprweather.com/pages/calc.htm

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 9 Июл 2012 09:48 - Поправил: Corvus


Бангладешский сайт про прогнозирование торнадо, анализируются метеоусловия наблюдавшихся случаев

http://bangladeshtornadoes.org/

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 9 Июл 2012 10:41


Считаем псевдоЛИ = ЭПТ850 - ЭПТ500.


Точнее сказать, это будет не псевдоЛИ, а индекс Шоултера, т.к. исходный уровень 850 гПа.

При расчёте ЛИ в качестве исходного берётся тот уровень (включая и приземный), на котором ЭПТ максимальная.

bullterrier
Участник
Письмо
Пермь
# Дата: 9 Июл 2012 13:57


Corvus
Спасибо большое. Придется тогда без Lifted Index обойтись. А в WRF оказывается есть переменная, прогнозирующая град. Но ее надо тестировать, смотреть, как она работает. Хорошая модель, только вот переменных CAPE/LI/CIN не хватает для полного счастья.:)

bullterrier
Участник
Письмо
Пермь
# Дата: 9 Июл 2012 14:09


А вероятность грозы буду считать по методу Вайтинга: 20-0% 40-100%, и сверху накладывать индексы Фатеева и TQ. Для первого раза, по моему, будет достаточно.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 9 Июл 2012 16:46


Спасибо большое. Придется тогда без Lifted Index обойтись.
bullterrier

Ну почему, его можно рассчитать аналогично индексу Шоултера:

1. Считаем ЭПТ850 (эквивалентно-потенциальная температура) по Т850, Тд850 и давлению 850 гПа.
2. Считаем ЭПТ500 по Т500, Тд500 и давлению 500 гПа.
3. SI = ЭПТ850 - ЭПТ500.


Только ЭПТ брать не для 850 гПа, а максимальную из приземной и значений на всех уровнях от земли до 700 гПа.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 11 Июл 2012 14:03 - Поправил: Corvus


Кстати, тот индекс, который в хСССР принято называть индексом Вайтинга, на Западе называют К-индексом. Предложен он неким Джорджем в 1960 году в книге:

George, J. J., 1960: Weather and Forecasting for Aeronautics. Academic Press, 673 pp.

<< . 1 . 2 . 3 . 4 . 5 . 6 . 7 . 8 . 9 . 10 ... 30 . 31 . >>
Ваш ответ

          Отменить *Что это?

 » Логин  » Пароль 
 
 
Полезная информация:



Поддержка: miniBB forum software © 2001-2017