Czym są meteorologiczne prognozy wiązkowe (EPS – Ensemble Prediction System)?

Opracowanie: dr Andrzej Mazur, Centrum Modelowania Meteorologicznego IMGW-PIB

Niektórych procesów atmosferycznych nie jesteśmy w stanie opisać z nieskończoną dokładnością. Dlaczego? Ze względu na:

  • niepewność (błąd) pomiaru i jego propagację do warunków początkowych modeli;
  • brak wiedzy o chwilowych warunkach pogodowych na tyle dokładnej, aby błąd w długookresowych obliczeniach był niezauważalny;
  • skończoną dokładność obliczeń, parametryzacje zamiast dokładnego opisu zjawisk;
  • błędy metod numerycznych, które są coraz bardziej widoczne wraz z wydłużaniem horyzontu czasowego obliczeń;
  • „efekt motyla” – wrażliwość na warunki początkowe (wg słów E. Lorenza, autora koncepcji „efektu motyla”: „dowolny układ fizyczny, który zachowuje się nieokresowo, jest nieprzewidywalny”).

Model meteorologiczny – jak wszystkie „modele” – opiera się na przyjętych założeniach, uproszczeniach i naszej wiedzy o warunkach początkowych, brzegowych, lokalnych nieciągłościach. Jednak nasza „niewiedza” jest olbrzymia w porównaniu do naszej „wiedzy”. ogromnej większości przypadków jesteśmy w stanie opisać jakieś zjawisko jedynie do pewnego poziomu szczegółowości i dokładności, wykorzystując dodatkowo przy tym parametryzacje, co do zasady upraszczające tenże opis.

W chwili obecnej nie jesteśmy w stanie przygotować stuprocentowo idealnej prognozy stanu atmosfery, jednak możemy określić, do jakiego stopnia możemy prognozie ufać, a narzędziem pozwalającym określić ten poziom zaufania są prognozy wiązkowe (ensemble) – zespół prognoz startujących z otoczenia stanu początkowego w granicach błędu pomiaru, czyli określenia tego stanu początkowego.

Ensemble (wiązka) jest to zatem zbiór realizacji prognozy, z których to realizacji każda startuje z różnego stanu początkowego, i/lub używa innej fizyki i/lub numeryki modelu, lub wręcz różnych modeli, np.:

  • Numeryka modelu – różne schematy numeryczne (eulerowski, Cranka-Nicholsona, leap-frog, Rungego-Kutty…).
  • Różna fizyka modelu – różne schematy parametryzacyjne różnych zjawisk (konwekcyjne, glebowe, turbulencji, radiacji, mikrofizyka chmurowa…).
  • Stochastyczne zaburzanie warunków początkowych i/lub brzegowych.
  • Różne modele, przygotowujące prognozę dla tego samego obszaru, na ten sam horyzont czasowy.
  • … i inne metody…

Wynikiem (rezultatem obliczeń EPS, czyli prognoz po zespole statystycznym) są między innymi:

  • Średnia po zespole (ensemble mean), czyli średnia ze zbioru wszystkich realizacji prognozy.

Obraz zawierający tekst, mapa, atlas

Opis wygenerowany automatycznie Obraz zawierający tekst, mapa, zrzut ekranu, żółty

Opis wygenerowany automatycznie

  • Niepewność prognozy, określana przez wartość rozrzutu (spread), zwykle definiowanego jako odchylenie standardowe względem średniej po zespole lub prognozy deterministycznej.

Obraz zawierający tekst, zrzut ekranu, diagram, mapa

Opis wygenerowany automatycznie Obraz zawierający tekst, mapa, zrzut ekranu, Jaskrawoniebieski

Opis wygenerowany automatycznie

  • Różnica między wartością maksymalną a minimalną w zbiorze, określana dla danego zjawiska. W pewnym sensie jest to analogon rozrzutu, z tym wyjątkiem, że wyznaczana jest na podstawie maximum maximorum (absolutnego maksimum) i minimum minimorum (absolutnego minimum) po zespole.

Obraz zawierający tekst, zrzut ekranu, mapa, Jaskrawoniebieski

Opis wygenerowany automatycznie Obraz zawierający tekst, mapa, zrzut ekranu, Jaskrawoniebieski

Opis wygenerowany automatycznie

  • Określenie prawdopodobieństwa przekroczeń i/lub wystąpienia jakiegoś zjawiska, obliczanego poprzez określenia pewnego kryterium. Określając taki warunek, sprawdza się, ilu członków zespołu go spełnia (nOK). Prawdopodobieństwo (%) – w naturalnym rozumieniu – wyznaczane jest przez iloraz tej liczby i całkowitej liczby nALL członków zespołu (nOK/nALL*100%). Prawdopodobieństwo takie można określić, poprzez wprowadzenie odpowiedniego kryterium (np. temperatura wyższa niż 30°C, wiatr silniejszy niż 15 m/s itp.) praktycznie dla każdego prognozowanego elementu. Jest to przydatne do szybkiego podsumowania prawdopodobieństwa wystąpienia określonych warunków, takich jak to, czy temperatura powietrza, czy temperatura odczuwalna osiągną niebezpiecznie wysoki poziom lub staną się wystarczająco niskie, aby zagrozić roślinności lub zdrowiu, lub też, czy opady mogą spowodować powodzie.

Obraz zawierający tekst, mapa

Opis wygenerowany automatycznie Obraz zawierający tekst, zrzut ekranu, mapa

Opis wygenerowany automatycznie

  • Uogólnione meteorogramy (EPS-gramy) to rozszerzenie standardowych meteorogramów – określających zmienność w czasie jednej wartości (np. z prognozy deterministycznej temperatury) dla jednej lokalizacji (np. stacji meteorologicznej) – do zawierającego przebieg czasowy średniej wraz z rozrzutem czy wartościami maksymalnymi i minimalnymi.

Obraz zawierający Wykres, zrzut ekranu, linia

Opis wygenerowany automatycznie Obraz zawierający Wykres, zrzut ekranu, diagram, linia

Opis wygenerowany automatycznie

  • Mapy baryczne i inne pola generowane na podstawie wyników prognozy wiązkowej na poziomach ciśnieniowych

 

Podsumowując:

  • Prognozy wiązkowe to zbiór (zespół statystyczny, ensemble) realizacji prognozy, każda startuje z różnego stanu początkowego, lub/i używa innej parametryzacji procesów fizycznych lub innych schematów numerycznych modelu lub różnych modeli.
  • Średnia ze zbioru realizacji, czyli tzw. „średnia po zespole” ma lepszą „sprawdzalność” (tzw. forecast skill, czy też zgodność z pomiarami) niż pojedyncze prognozy To założenie i podstawa systemu prognoz wiązkowych (EPS – Ensemble Prediction System, czyli prognoz po zespole statystycznym).
  • Przyjmuje się również generalną zasadę, że o ile założenia nie stanowią inaczej, żadna realizacja nie jest bardziej prawdopodobna od innej.

Obecnie, jako wyniki operacyjnego systemu prognoz wiązkowych w IMGW-PIB, dostępne są poniższe pola :

Na podst. wyników na poziomach hybrydowych i powierzchni:

Na podstawie wyników na poziomach ciśnieniowych:

  • Ciśnienie npg.
  • CAPE-MU (Most Unstable Level)
  • Indeks potencjału wyładowań (LPI)
  • Ciśnienie na poziomie morza
  • Temperatura 2m npg. (nad poziomem gruntu)
  • Temperatura punktu rosy 2m npg.
  • Opad śniegu
  • Temperatura powierzchni
  • Wilgotność względna 2m npg.
  • Zachmurzenie całkowite
  • Zachmurzenie niskiego poziomu
  • Opad całkowity
  • Opad deszczu
  • Temperatura odczuwalna 2m npg. (w cieniu)
  • Indeks stabilności mgły (FSI)
  • Zasięg widzialności
  • Prędkość wiatru 10m npg.

Dla tych pól, obok średnich, są też liczone prawdopodobieństwa przekroczenia zadanych progów, absolutne maksima i minima po zespole, rozrzut, różnice maksimum-minimum

  • Wysokość geopotencjalna poziomu standardowego
  • Mapa baryczna 300-950 hPa, 500-950 hPa, 700-950 hPa, 500-1000 hPa, 700-1000 hPa, 850-1000 hPa
  • Wysokość izotermy -20°C, -10°C i 0°C temperatury
  • Wysokość izotermy -20°C, -10°C i 0°C ttzw. (temperatury termometru zwilżonego)
  • Gradient temperatury 700-950 hPa, 700-850 hPa, 500-950 hPa, 500-850 hPa, 300-950 hPa, 300-850 hPa
  • Wirowość
  • Dywergencja
  • Temperatura powietrza
  • Temperatura punktu rosy
  • Temperatura term. zwilżonego
  • Niedosyt punktu rosy
  • Wilgotność względna
  • Prędkość i kierunek wiatru
  • Wiatr pionowy
  • Temperatura ekwiwalentno-potencjalna

Dla tych pól liczone są tylko wartości średnie po zespole.

Wszystkie pola są prezentowane w horyzoncie czasowym 36 godzin, rozdzielczość czasowa 3 godziny.

Dodatkowo, w wybranych lokalizacjach, przygotowywane są EPS-gramy prędkości wiatru 10m npg., temperatury powietrza 2m npg., ciśnienia zredukowanego do poziomu morza, wilgotności względnej 2m npg. oraz sumy i natężenia opadu.

— UDOSTĘPNIJ —