Radarová meteorologie - nowcasting Milan Šálek, Petr Novák Český Hydrometeorologický ústav Nowcasting • aktuální stav počasí + velmi krátkodobá předpověď • detailní analýza současného stavu počasí spolu s předpovědí na několik hodin dopředu, která vychází z extrapolace trendu vývoje (Conway, 1998) • detailní předpověď na dobu 0-6 hodin s upřesněním předpovědi na dobu 6-12 hodin (COST-78) Nowcasting • využívá především informací o současném stavu počasí získávaných pomocí metod dálkové detekce (údaje z meteorologických radarů, družic a systémů detekce blesků), jejichž informace se mohou kombinovat s pozemními či aerologickými pozorováními a výstupy numerických modelů • přesnější předpovědi a lepší prostorová lokalizace než u delších předpovědí • krátký interval platnosti předpovědi • spolehlivost nowcastingu většinou rychle klesá s rostoucím časem • je nezbytné zajistit •častá aktualizace předpovědí •dostatečná rychlost výpočtu a bezodkladná distribuce uživatelům či následným systémům •vhodná forma prezentace výsledků (jednoduchost, přehlednost, geografická navigace, …) Nowcasting Nowcasting srážek a nebezpečných konvektivních jevů • využívá se především predikce radarového (popřípadě družicového) echa • základní přístupy: • extrapolace radarového (popř. družicového) echa na základě pohybového pole (získaného různými metodami) • „plošná“ extrapolace • extrapolace konvektivních buněk • expertní systémy - extrapolace s využitím empiricky získaných poznatků (např. koncepčních modely procesů tvorby srážek zejména konvektivních systémů) či výstupů z numerických modelů (popř. přízemních a aerologických pozorovaní), často se využívá následného statistického postprocessingu • numerické modelování Extrapolace konvektivních buněk (SCIT) • identifikace bouřkových jader v radarových snímcích a přiřazení odpovídavících jader na následujících snímcích. • vektor posunu pro extrapolaci echa je určen ze změny polohy centra (těžiště) buňky na následujících obrázcích • extrapolace se provádí metodou přímých trajektorií • je možné určit trend v zesilování/zeslabování jednotlivých jader (radarového echa) – tzv. growth-decay factor – aplikování tohoto faktoru na výsledné pole však nedává příliš dobré výsledky (u žádné z metod) • k jednotlivým identifikovaným buňkám lze přiřadit další parametry (množství bleskových výbojů, Echo Top, VIL, družicová IR teplota ) • není možné určit vznik nové oblačnosti (konv. buňky) • např.: SCIT “Storm Cell Identification and Tracking” algoritmus (WSR- 88D/NSSL), TITAN (NCAR) „Plošná“ extrapolace radarového echa • výpočet plošného pohybového pole, které je následně použito pro extrapolaci posledního dostupného radarového snímku • různé způsoby výpočtu pohybového pole: • COTREC - porovnání dvou následujích radarových snímků pomocí nějaké veličiny vyjadřující podobnost (střední absolutní chyba, korelační koeficient,...) + hlazení pomocí podmínek kontinuity, Mecklenburg 2000, Zgonc and Rakovec 1999 • NWP - vektor posunu pro extrapolaci echa je určen z výstupu numerického modelu předpovědi počasí (nejčastěji lokálního) jakožto pohybové pole v hladině řídícího proudění (3-5km ~ AT 700- 500 hPa), • extrapolace se provádí metodou zpětných trajektorií • je možné určit trend v zesilování/zeslabování radarového echa (COTERC) • není možné určit vznik nové oblačnosti (konv. buňky) Extrapolace radarového echa HYBRID • nalezení “obrazců” v radarovém snímku obdobně jako v metodě SCIT (malo- (bouřkové buňky) i velkoprostorových (srážkové systémy)) • vektor posunu pro extrapolaci echa je určen z polohy “obrazce” na předchozím snímku – poloha je určena obdobně jako u metody COTREC (stř. Absolutní chyba, korelační koeficient) • extrapolace se poté provádí metodou přímých nebo zpětných trajektorií (v závislosti na konkrétní metodě) • je možné určit trend v zesilování/zeslabování radarového echa • není možné určit vznik nové oblačnosti (konv. buňky) • např. Lakshmanan 2002, Growth and Decay Tracker (MIT/LL ) Expertní systémy • využití koncepčních modelů : GANDOLF, AUTONOWCASTER • metody dálkové detekce + NWP : NIMROD, INCA • statistický postprocessing, „fuzzy logic“ : UFA SAM • teoreticky lze určit trend v zesilování/zeslabování radarového echa i vznik nové oblačnosti (konv. buňky) – nutná podrobná analýza proudění, hranice výtoků z Cb • rozporuplné výsledky – především u předpovědí na krátký časový interval u konvekce a silné konvekce (SYDNEY 2000 – Field Demonstration Project) Nowcasting v ČHMÚ ● optimalizace radarového snímání ● zrychlení dostupnosti radarových dat ● extrapolační metody pro předpověď radarového echa ● distribuce a vizualizace dat koncovým uživatelům Prokládané snímání během 10ti minut jsou měřeny dva 5ti-minutové subscany (16 různých elevací) 5ti minutová data využíváno Řízením letového provozu (nově se zkouší i v ČHMÚ – odhady srážek i nowcasting) kompromis mezi častou aktualizací a dobrým vertikálním rozlišení speciální nízkohladinový scan pro dopplerovskéo využití – experimentálně na Skalkách Objemové měření od 2009 přechod z desetiminutové obnovy na prokládané snímání a následně na čistou pětiminutovou obnovu dat (12 různých elevací) 5ti minutová data – primárně pro Řízením letového provozu (nyní i v ČHMÚ – odhady srážek i nowcasting) kompromis mezi častou aktualizací a dobrým vertikálním rozlišení speciální nízkohladinový scan pro dopplerovské využití – experimentálně na Skalkách Zpracování a distribuce radarových dat – spolupráce s AČR – rychlé datové spojení na oba radary – umožnění přenosu objemových dat v reálném čase do centra – vývoj vlastního software pro plné zpracování a vizualizaci objemových dat – optimalizace pro rychlost zpracování – flexibilita v generování produktů Zpracování a distribuce radarových dat rychlý přenos dat a optimalizace zpracovatelských programů a skriptů (minimalizace časových prostojů při zpracování, zrychlení zpracovatelských algoritmů) výrazně zvýšila čerstvost dat pro koncové uživatele. radarová data jsou dostupná pro uživatele v první minutě po konci objemového měření (ve starém zpracování byla data dostupná v páté minutě) další „zrychlení“ dosaženo měřením od horních méně důležitých elevací ke spodním – od 2007 Zpracování radarových dat 1 km horizontální rozlišení – nový operativní standard umožňuje podrobnější zkoumání radarového echa (hlavně konvektivních bouří) než u dříve používaného 2km rozlišení Extrapolace radarového echa - ČHMÚ • založeno na 1x1km datech maximálních odrazivostí • 2 základní části • výpočet pohybového pole (pole pro následnou extrapolaci) • kritická část predikce radarového echa – konečný výsledek nejvíce závisí na přesnosti a stabilitě použitého pohybového pole • pohybové pole by mělo být dostatečně hladké, ale zároveň by mělo rozlišovat rozdíly v pohybech jednotlivých bouří • z tohoto důvodu jsou použity 2 rozdílné metody výpočtu • časová extrapolace radarového echa • předpokládá se časová konstantnost pohybového pole • extrapolace se provádí metodou zpětných trajektorií • není určován trend v zesilování/zeslabování rad. echa Předpověď radarového echa - COTREC • pohybové pole je určováno porovnáním 2 následujících radarových snímků za pomocí střední absolutní chyby jakožto kriteria podobnosti • radarové snímky jsou rozděleny do čtvercových oblastí – předchozí snímek je posouván podél osy X a Y porovnáván s posledním snímkem • 'posuv' s nejmenší chybou určí vektor pohybového pole • radarové snímky jsou porovnávány v několika úrovních (různé velikosti oblastí - snímky jsou nejdříve porovnávány celé, poté jsou rozděleny na 6 čtverců a poté každý z těchto čtverců je rozdělen na 25 menších čtverců (44x44km)) • zhlazování výsledného pole - je použita metoda SOR – pohybové pole je změněno tak aby splňovalo podmínku kontinuity proudění Předpověď radarového echa - COTREC t t-∆t vx= ∆x / ∆t` vy= ∆y / ∆t` Předpověď radarového echa - ALADIN • předpoklad: pohyb oblačnosti je řízen prouděním ve výškové hladině cca. 3-5km nad hladinou moře • geopotenciál v hladině 700hPa (předpovězený pomocí NWP LAM ALADIN) je interpolován do projekce radarového snímku (1x1km gnom.) a přepočítán na pohybové pole pomocí geostrofické aproximace • jednoduchá, ale robusní metoda • nelze s její pomocí určit vývoj oblačnosti (zesílení/zeslabení) Předpověď radarového echa – porovnání • předpovězené radarové snímky jsou porovnávány se skutečně naměřenými daty (pixel vs. pixel) pomocí •střední kvadratické chyby (RMSE), •střední absolutní chyby (MAE) •indexy kvality (critical success index (CSI), probability of detection, false alarm ratio, ...) pro předpověď jisté úrovně dBZ (12, 24, 36 dBZ) POD=YY/(YY+NY) ; FAR=YN/(YN+NN) ; CSI=YY/(YY+YN+NY) • CSI lépe vypovídá o kvalitě předpovědi polohy než RMSE nebo MAE • jednotlivé předpovědi jsou porovnávány mezi sebou i s perzistentní předpovědí (počáteční radarový snímek je nezměněn použit jako předpověď) • subjektivní porovnání – též velmi důležité Předpověď radarového echa – porovnání CSI lépe vypovídá o kvalitě předpovědi polohy než RMSE nebo MAE Předpověď radarového echa – porovnání • všechny předpovědní metody vykazují zlepšení radarové informace – jsou lepší než persistentní předpověď • míra poklesu kvality předpovědi s časem je obdobná pro všechny metody. V konvektivních situacích je předpověď použitelná zhruba do +30 minut (v některých případech do +50 minut – organizovaná konvekce), v stratiformních situacích může být předpověď využitelná až do +90 minut. • ve většině případů, metoda ALADIN je mírně horší než metoda COTREC, někdy však dává nejlepší výsledky (především v případě velkoprostorových stratiformních srážek a v případech, kdy radarové odrazy jsou na okraji radarového dosahu – není ovlivněna “artefakty” radarových měření (radarový horizont, bright-band)) • metoda ALADIN je méně hardwareově náročná než metoda COTREC • obě metody předpovídají přesun radarového echa nikoliv jeho vývoj (zesílení/zeslabení - ani vznik nových) Využití extrapolovaného radarového echa pro výpočet kvantitativní předpovědi srážek • přepočet extrapolovaných radarových snímků na intenzitu srážšk pomocí Marshall-Palmerova vztahu– stejné jako při odhadu srážek • integrace po stanovene časové období – stejné jako při odhadu srážek – v současnosti testovány 1h úhrny na +0-1h, +1-2h a +2-3h • výpočet průměrných srážek na povodí • testování v hydrolog. modelech – první operativní aplikace 2007 Quantitative Precipitation Forecats COTREC motion field applied to the last available PseudoCAPPI 2km composite image Prediction up to 3h with 5-minute step Extrapolated radar images converted into rainrate fields using standard Z-R relation Integration for 0-1h, 2-3h and 2-3h motion Calculation of mean precipitation over predefined catchments COTREC-based QPF can be for first 3 hour better than QPF from NWP model Since spring 2007 used as an operational input into hydrological model Hydrog Quantitative Precipitation Forecats Tests of COTREC QPF in hydrological model Hydrog Discharge [m3s-1] COTREC QPF + NWP ALADIN QPF ALADIN QPF 136 (alert) 19 h 21 h 197 (Q10) 21 h 22 h 243 (Q20) 21 h 22 h 320 (Q50) 22 h 23 h 390 (Q100) 22 h 24 h The time of prediction of important discharge values Porovnání 1h předpovědí srážek metodou COTREC s předpověďmi NWP modelu ALADIN QPF – Statistical evaluation IV. – IX. 2006 QPF – Statistical evaluation Recent Development and Outlook Improvement of calculation of radar-raingauge bias mean-field bias >> local bias improvements of QPE as well as COTREC QPF Improvement of COTREC QPF – enlarging of forecasting domain Inclusion of data from foreign radars (difficulty with timing) use of NWP motion field as a first guess Hydrological ensemble forecast based on different QPFs Recent Development and Outlook Improvement of QPF using statistical-advection model inspired by ADSTAT algorithm from US NWS predictors radar-raingauge QPE radar COTREC QPF ALADIN NWP model forecasts Meteosat MSG data ligthning data precipitation forecasts for 0-1h,1-2h, 2-3h - probabilistic as well as quantitative possibility to forecast also newly developed precipitation plan to be used operationally as a member of ensemble forecasts QPF – statistical-advection model 1h QPF - observation 3h QPF - observation Testování kvantitativního využití metody COTREC v hydrolog. modelu Hydrog PODHRADÍ - discharge forecast 11.8.2002 18 CEST 0 50 100 150 200 250 300 350 400 450 7.8.2002 0:00 9.8.2002 0:00 11.8.2002 0:00 13.8.2002 0:00 15.8.2002 0:00 17.8.2002 0:00 Q[m3 s-1 ] Real data ALADIN COTREC+ALADIN Q10 Q50 ALERT PODHRADÍ - discharge forecast 11.8.2002 24 CEST 0 50 100 150 200 250 300 350 400 450 7.8.2002 0:00 9.8.2002 0:00 11.8.2002 0:00 13.8.2002 0:00 15.8.2002 0:00 17.8.2002 0:00 Q[m3 s-1 ] Real data ALADIN COTREC+ALADIN Q10 Q50 ALERT Testování kvantitativního využití metody COTREC v hydrolog. modelu Hydrog Porovnání 1h předpovědí srážek metodou COTREC s předpověďmi NWP modelu ALADIN 1.4.- 30.9. 2006 0 0,1 0,2 0,3 0,4 0,5 0,6 Cot1 Cot2 Cot3 Aladin RMSE[mm] Porovnání 1h předpovědí srážek metodou COTREC s předpověďmi NWP modelu ALADIN CSI 1mm korelační koeficient Porovnání 1h předpovědí srážek metodou COTREC s předpověďmi NWP modelu ALADIN 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 4 5 6 7 8 9 měsíce RMSE[mm] Cot1 Cot2 Cot3 Aladin 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 4 5 6 7 8 9 měsíce CSI(práh1mm) Cot1 Cot2 Cot3 Aladin 0 0,1 0,2 0,3 0,4 0,5 0,6 Cot1 Cot2 Cot3 Aladin Per1 Per2 Per3 RMSE[mm] Porovnání 1h předpovědí srážek metodou COTREC s předpověďmi NWP modelu ALADIN a persistencí Problémy plošné extrapolace hlazení pohybového pole a využití zpětných trajektorií může způsobit selhání extrapolace pokud se nějaká konvektivní bouře pohybuje výrazně odlišným směrem než okolní buňky typický případ – supercely relativně řídký výskyt nicméně nebezpečné projevy počasí jeden z důvodů vývoje metody CELLTRACK CELLTRACK - identifikace buněk Rozdíl mezi „reálnou“ konvektivní buňkou (oblasti uspořádaných vzestupných (sestupných) pohybů vzduchu) a buňkou definovanou v nowcastingových metodách (tedy oblasti jisté vysoké odrazivosti) Různé metody využívají různá pravidla pro identifikace jader odrazivosti CELLTRACK - identifikace buněk Ve finání verzi používán jednoduchý práh odrazivosti o hodnotě 44 dBZ zkoušen též práh 36, 40, 44, 48 dBZ 44dBZ se jeví jako vhodný kompromis mezi identifikací i slabších buněk a neidentifikování více blízkých buněk jako jednoho jádra Testován i algoritmus s pohyblivým prahem odvozeným z metody TRACE3D lokální maximum odrazivosti – 10 dBZ Zvýšil se počet identifikovaných buněk (ze 4124 na 4406 => 93,6%), ale horší výsledky při jejich sledování. 3D identifikace – testováno ale s problematickými výsledky hledání souvislosti mezi jádry nalezenými na CAPPI hladinách jdoucích po sobě mnoho buněk nalezených pouze v jedné hladině souvisí zřejmě s prokládaným objemovým snímáním radarů a způsobem výpočtu CAPPI hladin - interpolace CELLTRACK - sledování jader Nejdříve se určují „shluky“ buňek, které spolu mohou souviset První odhad pohybu jader v předchozím termínu se provádí pomocí metody COTREC Pro každé takto posunuté jádro se hledají blízká skutečná jádra na následujícím snímku. Poloměr prohledávané oblasti závisí na rychlosti proudění (nicméně závislost úspěšnosti sledování jader na tomto poloměru není příliš silná); Prohledávání probíhá i „opačným“ směrem V dalším kroku se zpracovávají buňky v jednotlivých shlucích nejdříve se hledají nejpodobnější jádra (podobnost >0.85), která jsou dost blízko u sebe pak se hledají nejbližší buňky a zkoumá se jejich podobnost pokud je menší než 0,8 => štěpení, slučování CELLTRACK – extrapolace Extrapolace polohy jádra se provádí vektorem určeným z předchozí a současné pozice jádra (hledá se podobnost jader na po sobě jdoucích snímcích - lokální COTREC) Jednoduchá extrapolace (určená ze změny polohy těžiště jádra) není dostatečná zejména kvůli slučování a štěpení jader Pokud jádro nemá předchůdce, je použit průměr pohybových vektorů jader s alespoň jedním předchůdcem. Pokud neexistuje žádný předchůdce, předpověď se nekoná (typicky u prvního času se zachycením jádra). Předpověď je počítána na 10 - 90 minut s krokem po 10 minutách CELLTRACK - úspěšnost sledování jader odrazivosti Porovnání algoritmu sledování jader odrazivosti na následných snímcích s manuálním sledováním – pokud se přiřazení od sebe liší, je přiřazení provedené algoritmem označeno za chybné. hits misses wrong assignments CSI 44 dBZ 1018 76 59 0,88 TRACE3D 1035 102 99 0,84 CELLTRACK - úspěšnost sledování jader odrazivosti Porovnání předpovědi CELLTRACK s metodou COTREC CSI pro práh 44 dBZ Příprava operativního využití Experimentální zařazení metody CELLTRACK do vizualizačního software JSMeteoView Technické provedení: obrázek s transparentním pozadím vytvořen v jazyce C a přidán jako další vrstva do obrázku pomocí DHTML a JavaScriptu Budoucí vývoj Dokončení integrace do prohlížeče JSMeteoView Statistika časového vývoje vlastností identifikovaných buněk doba života velikost jader hodnota nejvyšší odrazivosti výška oblasti nejvyšší odrazivosti VIL, HAIL PROB, ECHO TOP snaha zahrnout i data ze sítě detekce blesků CELDN, a družice MSG Cíl -> uživatel by měl mít k dispozici přehledně předpovědí z metody CELLTRACK i COTREC (případně i aktuální a typický časový vývoj důležitých vlastností), které by mu měly pomoci při rozhodování CELLTRACK – charakteristiky konvektivních bouří Vizualizace dat je třeba dodat data uživatelům co nejrychleji a umožnit jim přesnou geografickou lokalizaci detekovaných jevů JSMeteoView – internetový prohlížeč radarových (a některých dalších) dat vývoj započal v 2001 – stále pokračuje univerzální zobrazení v moderních prohlížečích (Gecko-based browsers Mozilla/Netscape 6.x or Microsoft Internet Explorer 5.x/6.x) nezávislé na operačním systému (Windows, LINUX/UNIX, MAC OS) možnost zobrazení a kombinace i dalších meteorologických dat (Meteosat MSG, data detekce blesků, předpovědi NWP LAM ALADIN, SYNOP data) GIS funkce základní nástroj pro zobrazování dat distančních pozorování v ČHMÚ a též u některých externích uživatelů