Charakterystyka produktów radarowych - Instytut Meteorologii i ...
Charakterystyka produktów radarowych - Instytut Meteorologii i ...
Charakterystyka produktów radarowych - Instytut Meteorologii i ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Irena Tuszyńska<br />
<strong>Charakterystyka</strong> produktów<br />
<strong>radarowych</strong><br />
<strong>Instytut</strong> <strong>Meteorologii</strong> i Gospodarki Wodnej<br />
Państwowy <strong>Instytut</strong> Badawczy<br />
Warszawa marzec 2011
Autor opracowania<br />
mgr inż. Irena TUSZYŃSKA<br />
Irena.Tuszynska@imgw.pl<br />
Ośrodek Teledetekcji Naziemnej<br />
<strong>Instytut</strong> <strong>Meteorologii</strong> i Gospodarki Wodnej<br />
Państwowy <strong>Instytut</strong> Badawczy<br />
01-673 Warszawa<br />
ul. Podleśna 61<br />
Opracowanie prezentacji w HTML<br />
mgr Paweł Rychlewski<br />
Pawel.Rychlewski@imgw.pl<br />
Ośrodek Teledetekcji Naziemnej<br />
<strong>Instytut</strong> <strong>Meteorologii</strong> i Gospodarki Wodnej<br />
Państwowy <strong>Instytut</strong> Badawczy<br />
01-673 Warszawa<br />
ul. Podleśna 61<br />
Recenzent opracowania<br />
mgr Zdzisław Dziewit<br />
Zdzislaw.Dziewit@imgw.pl<br />
Ośrodek Teledetekcji Naziemnej<br />
<strong>Instytut</strong> <strong>Meteorologii</strong> i Gospodarki Wodnej<br />
Państwowy <strong>Instytut</strong> Badawczy<br />
01-673 Warszawa<br />
ul. Podleśna 61<br />
2
Spis treści<br />
1 Metodyka pozyskiwania informacji radarowej i jej przetwarzanie 5<br />
1.1 Radar meteorologiczny 5<br />
1.2 Radary z pojedynczą i podwójną polaryzacją 8<br />
1.3 Pozyskiwanie i przetwarzanie informacji radarowej 9<br />
1.4 Typy produktów generowanych przez system radarowy 12<br />
2 Podstawowe produkty radarowe 21<br />
2.1 Przekrój stożkowy - PPI 21<br />
2.2 Przekroje poziome typu CAPPI i PCAPPI 24<br />
2.3 Odbiciowość średnia w zadanej warstwie - LMR 25<br />
2.4 Odbiciowość maksymalna - MAX i CMAX 26<br />
2.5 Wysokość echa radarowego - EHT 28<br />
2.5.1 Wierzchołek echa radarowego 28<br />
2.5.2 Wysokość maksymalnej odbiciowości 29<br />
2.5.3 Wysokość podstawy echa radarowego 29<br />
2.6 Przekrój pionowy – RHI 31<br />
2.7 Przekrój pionowy przez dwa punkty - VCUT i przekrój wzdłuż łamanej –<br />
MLVCU<br />
32<br />
3 Produkty hydrologiczne 38<br />
3.1 Natężenie opadu - SRI 36<br />
3.2 Suma opadu - PAC 35<br />
3.3 Wodność scałkowana w pionie - VIL 36<br />
3.4 Akumulacja opadu w podzlewniach - RSA 37<br />
4 Produkty wiatrowe 40<br />
4.1 Produkty wiatrowe według geometrii ich tworzenia 40<br />
4.1.1 Przekroje stożkowe rozkładu prędkości radialnej i szerokości widma prędkości<br />
radialnych typu PPI<br />
40<br />
4.1.2 Rozkład prędkości radialnej i rozkład szerokości widma prędkości radialnych<br />
na stałej wysokości typu PCAPPI<br />
4.1.3 Rozkład prędkości i szerokości widma prędkości radialnych typu CMAX<br />
i MAX<br />
42<br />
43<br />
3
4.1.4 Rozkład prędkości radialnych i szerokości widma prędkości radialnych na<br />
przekrojach typu VCUT i MLVCUT<br />
46<br />
4.2 Pionowy profil wiatru - VVP 47<br />
4.3 Wiatr poziomy - HWIND i technika wiatru jednorodnego - UWT 49<br />
4.4 Prędkość wiatru radialnego w funkcji azymutu - VAD 51<br />
4.5 Produkty uskoków wiatru typu SHEAR 52<br />
4.5.1 Produkty SHEAR i ich kombinacje 52<br />
4.5.2 Poziomy gradient wiatru – HSHEAR 54<br />
4.5.3 Pionowy gradient wiatru – VSHEAR 55<br />
4..6 Turbulencje - LTB 56<br />
5 Produkty analizy groźnych zjawisk 58<br />
5.1 Wskaźnik groźnych zjawisk - SWI 58<br />
5.2 Prawdopodobieństwo wystąpienia gradu - ZHAIL 62<br />
6 Produkty prognostyczne 63<br />
6.1 Śledzenie komórek opadowych - RTR 63<br />
6.2 Śledzenie komórek burzowych - CTR 65<br />
7 Rozpoznawanie zjawisk w oparciu o produkty radarowe 66<br />
7.1 Wykorzystanie pomiarów z podwójną polaryzacją 66<br />
7.2 Rozpoznanie typu opadu – CAPPI(ET) i RHI(ET) 68<br />
7.2 Możliwości rozpoznawania zjawisk w oparciu o produkty z radarów<br />
pracujących z pojedynczą polaryzacją fali<br />
69<br />
8 Zbiorcza mapa radarowa 74<br />
4
1. Metodyka pozyskiwania informacji radarowej i jej przetwarzanie<br />
1.1 Radar meteorologiczny<br />
RADAR (RAdio Detection And Ranging) - urządzenie, którego zasada działania jest<br />
oparta na wykorzystaniu zdolności do odbijania przez obserwowane obiekty wysyłanego w<br />
ich kierunku promieniowania elektromagnetycznego. Radar znany był już w latach<br />
trzydziestych XX wieku. Dzisiaj jest wynikiem osiągnięć wielu lat badań. Pierwsze radary<br />
służyły tylko do zdalnego wykrywania samolotów i okrętów. Około lat pięćdziesiątych radary<br />
zaczęli wykorzystywać do swoich celów także meteorolodzy, tj. do wykrywania oraz<br />
śledzenia burz i opadów. Pojawiła się nowa dziedzina wiedzy: meteorologia radarowa<br />
zajmująca się badaniem różnic między odbiciami sygnałów <strong>radarowych</strong> pozyskanych od<br />
różnego rodzaju chmur i opadów.<br />
Wysłana przez radar wiązka fali elektromagnetycznej w momencie zetknięcia się z<br />
cząsteczkami chmurowymi wzbudza w nich elektryczne i magnetyczne dipole. One z kolei<br />
wypromieniowują w przestrzeń falę, która powraca do anteny radaru tworząc sygnał<br />
użyteczny, tzw. echo radarowe. Największy wkład w sygnał radarowy od obiektów<br />
meteorologicznych wnoszą duże cząsteczki – cząsteczki opadowe. Jednak tylko drobna część<br />
energii przechodzącej przez chmurę opadową jest przez nią odbijana. Większa cześć energii<br />
przechodzi dalej. Dzięki temu możemy obserwować chmury położone jedna za drugą.<br />
Podstawową cechą, wyróżniającą radar meteorologiczny od pozostałych radarów jest<br />
charakterystyka kierunkowa, czyli kształt<br />
emitowanej wiązki fal elektromagnetycznych i<br />
sposób przeszukiwania przestrzeni. Radary<br />
meteorologiczne emitują wiązkę "stożkową", o<br />
rozbieżności mniejszej lub równej 1° w pionie i w<br />
poziomie. Wiązka taka charakteryzuje się dużą<br />
koncentracją energii. Jej przekrój ma kształt<br />
kołowy, który jest tym większy im większa jest<br />
odległość od radaru. Dla danego ustawienia kątowego wiązki radarowej oświetlane są kolejne<br />
„próbki obiektu meteorologicznego”, w miarę jak impuls sondujący „przebija się” przez<br />
obiekt. Jest to wielką zaletą radarów meteorologicznych. Nie obserwuje on tylko granic<br />
chmur (jak aparat fotograficzny), ale prześwietla je działając podobnie jak rentgen czy USG.<br />
Pozyskane informacje z kolejnych próbek powracają do anteny w postaci echa radarowego<br />
(zwrotnej fali elektromagnetycznej). Ich rozmiary zależą od rozpiętości wiązki radarowej.<br />
5
Uwzględniając właściwości materiału kropli wody badacze tego zjawiska znaleźli zależność<br />
opisującą powierzchnię skuteczną kropli wody od jej średnicy. Rezultatem było<br />
sformułowanie zależności mocy echa radarowego P r od średnic kropli i ich liczby<br />
w jednostce objętości.<br />
P<br />
C<br />
N<br />
6<br />
r<br />
D<br />
2 i<br />
R i1<br />
gdzie: P r - moc echa radarowego [mW],<br />
D i - średnica kropli, R – odległość sondowanej próbki od radaru,<br />
C - współczynnik zależny od parametrów radaru i stanu skupienia hydrometeorów,<br />
N - liczba kropli w jednostce objętości chmury.<br />
Suma występująca w tej zależności jest nazywana odbiciowością i oznaczana symbolem Z<br />
czyli:<br />
Z <br />
N<br />
D i<br />
i1<br />
6 C<br />
[mm 6 /m 3 ] P r<br />
Z<br />
2<br />
R<br />
Wartości odbiciowości Z zmieniają się w bardzo szerokich granicach dlatego też najczęściej<br />
stosuje się tutaj skalę logarytmiczną, a jednostkę odbiciowości oznaczono symbolem dBZ<br />
(decybele ze stosunku Z do 1 mm 6 /m 3 ) i definiuje jako:<br />
1 dBZ = 10logZ<br />
a w meteorologii radarowej nazywamy ją odbiciowością radarową.<br />
W przypadku radarów meteorologicznych wykrywane są przestrzenne struktury<br />
meteorologiczne, które rozciągając się w przestrzeni wypełniają kolejne komórki (próbki)<br />
rozróżnialności w kolejnych wiązach skanowania radowego. Nie są one jednorodne tj. nie<br />
składają się tylko z kropli wody lecz są najczęściej mieszaniną hydrometeorów (kropli<br />
deszczu, płatków śniegu, krupu śnieżnej, gradu itp.). Pozyskana od nich wartość odbiciowości<br />
uzależniona jest jednocześnie od typu hydrometeorów, ich rozmiarów i rozkładu w<br />
zlokalizowanej strukturze (w jednostce objętości).<br />
Atmosfera od początku swego istnienia jest w nieustannym ruchu, pionowym jak<br />
i poziomym. Poziomy ruch powietrza to wiatr. Znajomość kierunku i prędkości wiatru to<br />
jedna z fundamentalnych informacji meteorologicznych. Standardowe pomiary kierunku i<br />
prędkość wiatru, jak wszystkie standardowe pomiary w meteorologii, są ograniczone do<br />
miejsca, w którym się je wykonuje. Natomiast radar nie jest w stanie wykryć ruchu (wiatru) w<br />
powietrzu. Jednak wiatr wywiera wpływ na obiekty znajdujące się bezpośrednio w powietrzu<br />
lub wchodzące z nim w interakcje. Obiektami tymi są, między innymi chmury, a ich<br />
6
przemieszczanie się to efekt działania wiatru. Technika radarowa obserwacji wiatru polega na<br />
śledzeniu zmian położenia hydrometeorów, będących w ciągłym ruch, w stosunku do<br />
obserwatora (radaru) jak i względem siebie, który jest wynikiem niejednorodnego rozkładu<br />
prędkości wiatru w obszarze danej struktury opadowej.<br />
Teledetekcyjne metody wyznaczania parametrów<br />
pola wiatr bazują głównie na wyznaczeniu<br />
dopplerowskiego przesunięcia częstotliwości<br />
emitowanego sygnału powodowanego ruchem<br />
cząstek chmurowych.<br />
Częstotliwość odbieranego sygnału od<br />
poruszających się struktur jest większa niż<br />
częstotliwość sygnału wysyłanego w ich<br />
kierunku w przypadku gdy obiekt zbliża się do radaru i mniejsza, gdy oddala się. Jest to<br />
typowy przykład występowania tzw. efektu Dopplera w odniesieniu do fal<br />
elektromagnetycznych. Zmianę częstotliwości fali elektromagnetycznej w kierunku od i do<br />
radaru nazywamy częstotliwością dopplerowską lub częstotliwością Dopplera.<br />
Zmiana częstotliwości fali elektromagnetycznej wywołanej<br />
zjawiskiem Dopplera<br />
Źródło rysunku: http://pl.wikipedia.org/wiki/Radar<br />
Wszystkie urządzenia oddanej do użytku w maju<br />
2004 roku sieci POLRAD to dopplerowskie radary meteorologiczne. Poza pomiarem<br />
odbiciowości radarowej dokonują także pomiaru składowej radialnej wiatru (ruchu cząstek<br />
opadowych) i spektrum widma prędkości radialnych w strukturze opadowej.<br />
Pomiar prędkości dopplerowskich dokonywany jest w odbiorniku radarowym. Prędkość<br />
przemieszczania się hydrometeorów, prędkość radialna wykorzystywana jest w radarach<br />
dopplerowskich także do eliminacji ech nie opadowych poprzez zastosowanie w trakcie<br />
pomiaru tzw. filtrów dopplerowskich i jest najskuteczniejszym kryterium selekcji ech<br />
pożądanych od zakłóceń.<br />
Zmierzona przez radar dopplerowski prędkość przemieszczania cząstek<br />
hydrometeorów w całym obiekcie jest prędkością rzeczywistą, gdyż hydrometeory są<br />
7
unoszone bez opóźnień. Jednak, aby móc określić wszystkie częstotliwości dopplerowskie,<br />
zależne od prędkości cząstek hydrometeorów, pojawiają się dodatkowe wymagania dla radaru<br />
dopplerowskiego. W tym momencie chodzi o tzw. częstotliwość powtarzania impulsu (PRF),<br />
która powinna być dostatecznie wysoka, aby móc określić wszystkie występujące<br />
częstotliwości dopplerowskie, zależne od prędkości cząstek hydrometeorów.<br />
1.2 Radary z pojedynczą i podwójną polaryzacją<br />
Aktualnie w polskiej sieci radarowej POLRAD większość urządzeń to radary<br />
pracujące z pojedynczą polaryzacją fali (z falą spolaryzowaną poziomo). W miarę ich<br />
zużywania się wymieniane są na radary pracujące z falą o podwójnej polaryzacji. Sieć w<br />
ramach dalszego rozwoju będzie zagęszczana.<br />
Zastosowanie do pomiaru techniki o podwójnej polaryzacji fali elektromagnetycznej<br />
otworzyło nowe możliwości dla meteorologii radarowej jeśli chodzi o identyfikację<br />
i klasyfikację obiektów meteorologicznych. W praktyce sprowadza się to do emisji sygnału<br />
radarowego, spolaryzowanego poziomo i pionowo, na przemian lub jednocześnie z kanałów<br />
nadawczych, a następnie odbioru sygnału odbitego, przez dwa kanały, z których każdy<br />
dostosowany jest do swojej polaryzacji.<br />
Pozyskujemy wówczas zróżnicowaną informację odnośnie charakterystyk obiektu ze względu<br />
na polaryzację sygnału sondującego. Niesymetrycznie rozpraszające obiekty meteorologiczne<br />
powodują pojawienie się składnika depolaryzacyjnego zależnego od formy i rozmiaru<br />
orientacji i dielektrycznych właściwości cząstek rozpraszających (hydrometeorów).<br />
Klasyczny radar meteorologiczny wysyła wiązkę mikrofal o polaryzacji liniowej<br />
poziomej. Oznacza to, że wektor elektryczny fali (elektromagnetycznej) oscyluje w jednej<br />
płaszczyźnie (poziomej). Odbicie fali elektromagnetycznej od powierzchni dielektryka<br />
(wody) jest w rzeczywistości generowaniem fali w wyniku zmian orientacji dipoli<br />
elektrycznych wody, wymuszonych przez falę padającą. Dla opadów stworzono teorię<br />
wiążącą depolaryzację z położeniem głównych osi cząstek, które aproksymuje się małymi<br />
elipsoidami i różnicuje się sygnał odbity w zależności od polaryzacji fali padającej, gdyż<br />
radar wysyłając falę “widzi” hydrometeor o wymiarze takim, który jest dominujący w<br />
płaszczyźnie polaryzacji fali wysłanej w jego kierunku.<br />
Wszyscy zapewne wiemy, ze kształt swobodnie spadającej<br />
kropli zależy od jej wielkości. Przy czym najmniejsze kropelki<br />
są prawie idealnie okrągłe, zaś duże krople spadając mają<br />
tendencję do „spłaszczania”. Zmniejszenie „pionowej średnicy” kropli jest tym większe im<br />
8
większa objętość kropli i jej prędkość spadania. W przypadku omiatania wiązką radarową<br />
takich spłaszczonych kropel powracający sygnał zależy od płaszczyzny polaryzacji fali.<br />
Wykonując pomiary radarowe radarem o podwójnej polaryzacji dokonujemy jakby pomiaru<br />
obu średnic kropel wzdłuż dwóch płaszczyzn (pionowej i poziomej) wynikających z<br />
opromieniowania hydrometeorów w strukturze meteorologicznej, jednocześnie falą<br />
o pionowej i poziomej polaryzacji. Skutkiem tego pozyskuje się dwie wartości sygnału<br />
odbitego od tych samych kropli różniące się miedzy sobą, a<br />
następnie porównuje się sygnały nadawane w dwóch<br />
wzajemnie prostopadłych kierunkach z sygnałami zwrotnymi.<br />
W przypadku kropel, których średnice są większe od 1 mm<br />
pomiar falą spolaryzowaną poziomo daje większą wartość<br />
odbiciowości niż falą spolaryzowaną pionowo (odbiciowość Z jest sumą szóstych potęg<br />
średnic kropel) w opromieniowywanej objętości atmosfery. Porównując sygnały odbierane w<br />
dwóch wzajemnie prostopadłych kierunkach możemy ocenić wielkości kropelek<br />
dominujących. Pozyskanie takiej informacji skutkuje dokładniejszym przeliczaniem<br />
odbiciowości na natężenie opadu.<br />
W oparciu o wartości wysłanego i zwrotnego sygnału podwójnie spolaryzowanego, w<br />
meteorologii radarowej wykorzystuje się operacyjnie różnego rodzaju wartości (wskaźniki),<br />
pozyskiwane z analiz tych sygnałów.<br />
1.3 Pozyskiwanie i przetwarzanie informacji radarowej<br />
Podstawowym zadaniem radaru jest zbieranie danych, a następnie ich przetworzenie<br />
przez specjalistyczne oprogramowanie. Antenę radarową możemy sterować tak, aby zbierała<br />
informację tuż przy ziemi lub ustawić na dowolnym kącie w stosunku do powierzchni ziemi.<br />
Dzięki temu przegląd przestrzeni nad radarem można wykonywać w przekroju pionowym lub<br />
poziomym.<br />
Podstawowymi sposobami pozyskania informacji radarowej są przeszukiwania<br />
przestrzenne atmosfery. Wynikiem takiego przeszukiwania są surowe (źródłowe) dane we<br />
współrzędnych sferycznych zebrane z kolejnych elewacji dla:<br />
• dalekiego zasięgu – skanowanie klasyczne,<br />
• bliskiego zasięgu – skanowanie dopplerowskie<br />
oraz dodatkowe dane ze<br />
• skanowania pionowego wykonanego na zadanym azymucie.<br />
9
Pomiary skanowania klasycznego zoptymalizowane są pod kątem pomiaru<br />
odbiciowości Z [dBZ] do zasięgu 250 km wokół radaru i wysokości 20 km. Najniższy kąt<br />
elewacji uwzględniający rozpiętość wiązki radarowej to 0,5. Wszystkie pozostałe kąty<br />
dobrane są tak, aby zapewnić optymalne pokrycie skanowanej przestrzeni (około 90%<br />
objętości skanowanej przestrzeni).<br />
20<br />
23,8 o<br />
18,5 o<br />
14,1 o<br />
10,6 o<br />
7,7 o<br />
5,3 o<br />
3,4 o<br />
15<br />
2,4 o<br />
10<br />
1,4 o<br />
5<br />
0,5 o<br />
0<br />
-5<br />
Odległość w [km]<br />
-10<br />
0 50 100<br />
150 200 250<br />
Przykładowa geometria i czas skanowania klasycznego (narastająco) podany w [sekundach].<br />
Powyższy diagram zawiera informacje o zasięgu pomiaru radarowego w odległości<br />
i wysokości. Przedstawia krzywiznę ziemi (brązowa linia) i pokrycie obszaru poniżej<br />
minimalnej i maksymalnej elewacji wynikające z rozpiętości wiązki (jasno zielony kolor).<br />
Natomiast w trakcie skanowania dopplerowskiego zbierane są odpowiednio dane o:<br />
- odbiciowości radarowej; Z [dBZ],<br />
- prędkości radialnej; V [m/s],<br />
- szerokości spektralnej widma prędkości radialnych ; W [m/s]<br />
do zasięgu 125 km wokół radaru i do wysokości 20km. Najniższy kąt elewacji to 0,5 (tak<br />
samo jak dla skanowania klasycznego), a najwyższy 32,8. Parametry skanowania<br />
zoptymalizowane są pod kątem pomiaru prędkości radialnej.<br />
10
Przykładowa geometria<br />
skanowania dopplerowskiego<br />
W trakcie skanowania pionowego radarem dopplerowskim dane mogą być zbierane<br />
po:<br />
- odbiciowości radarowej; Z [dBZ],<br />
- prędkości radialnej; V [m/s],<br />
- szerokości spektralnej widma prędkości radialnych ; W [m/s]<br />
do określonego zasięgu w odległości i wysokości. O tym, jakie dane są zbierane, decydujemy<br />
w momencie definiowania skanu pionowego.<br />
Wstępne przetwarzanie odebranego sygnału odbywa się na stacji radarowej i jest<br />
realizowane przez radarowy procesor sygnału RSP (Radar Signal Processor). RSP wykonuje<br />
wstępną obróbkę danych pochodzących ze skanowania. Wynik tego procesu stanowi źródło<br />
danych dla dalszego przetwarzania w systemie radarowym. Generuje on podstawowe<br />
wielkości, takie jak: odbiciowość Z, prędkość radialna V, szerokość widmowa W, czy<br />
składowa podwójnej polaryzacji ZDR itd. Jednocześnie, RSP zapewnia wykonywanie<br />
wszystkich niezbędnych obliczeń takich jak: filtrowanie sygnału ech stałych, uśrednianie<br />
sygnału, demodulacja sygnału do składowych dopplerowskich, dyskryminację wg. wartości<br />
progowych i innych wielkości zgodnie z zaprogramowanymi przez użytkownika<br />
ustawieniami w harmonogramie zadań.<br />
Pozyskaną informację system radarowy przetwarza na bardzo dużą ilość różnorodnych<br />
produktów według specjalistycznych algorytmów zbudowanych nie tylko w oparciu o wiedzę<br />
z dziedziny meteorologii radarowej, meteorologii ogólnej, klimatologii, ale także<br />
z uwzględnieniem nauk ścisłych.<br />
Standardowo w systemie zdefiniowane są produkty podstawowe, które są niezbędne<br />
do wyrażenia aktualnego stanu atmosfery. Produkty te są archiwizowane - zapisane w bazie<br />
11
w postaci plików graficznych typu PNG i w wersji źródłowej typu XML. Pozostałe produkty<br />
generowane są w zależności od potrzeb wynikających z:<br />
- konieczności wprowadzania danych <strong>radarowych</strong> do innych systemów,<br />
- zadań statutowych IMGW,<br />
- indywidualnych potrzeb użytkowników,<br />
- prowadzonych aktualnie prac badawczych.<br />
Dodatkowe produkty pozyskiwane w ramach systemu radarowego podzielić można ze<br />
względu na możliwości jakie wynikają z:<br />
- definicji tych produktów,<br />
- zobrazowania różnych stanów atmosfery na określonych wysokościach,<br />
- zobrazowania stanu atmosfery w określonej warstwie, wartości średnich lub<br />
maksymalnych,<br />
- geometrii tworzenia produktów w oparciu o te same algorytmy produktów,<br />
- zmian wartości progowych parametrów produktów określanych w momencie<br />
definiowania produktów.<br />
Wszystkie produkty ze skanu klasycznego pozyskujemy w zasięgu do 200-250 km<br />
wokół radaru, a ze skanu dopplerowskiego, w większości (tam gdzie jest to możliwe) do<br />
zasięgu 100-125 km wokół radaru.<br />
1.4. Typy produktów generowanych przez system radarowy<br />
Podstawowym zadaniem radarowej sieci meteorologicznej jest generowanie<br />
produktów <strong>radarowych</strong>. Wytwarzanie produktów poprzedzone jest wieloma uzupełniającymi<br />
się zadaniami (co wynika z lektury poprzednich rozdziałów) zmierzającymi do utrzymania<br />
przepływu informacji od radaru do modułu je generującego, a następnie ich zobrazowanie lub<br />
przesłanie do użytkownika w postaci końcowych produktów.<br />
Biorąc pod uwagę technikę i sposób tworzenia produktów możemy wydzielić dwie główne<br />
grupy, tak zwane:<br />
- produkty pierwszego rzędu “first-level” - pozyskiwane z przetwarzania surowych danych<br />
(typu 3D) z wykorzystaniem określonych algorytmów i szczegółowej definicji produktu,<br />
- produkty drugiego rzędu “second-level“ - obliczane na podstawie produktów pierwszego<br />
rzędu np. do wygenerowania PAC (Sumy opadów) niezbędny jest produkt SRI (Natężenie<br />
opadu).<br />
Zasadniczo, w ramach wymienionych dwu grup produktów, generowane produkty możemy<br />
podzielić na pięć podstawowych grup (ich nazwy traktować należy jako umowne):<br />
12
- produkty podstawowe – klasyczne produkty radarowe,<br />
- produkty hydrologiczne,<br />
- produkty wiatrowe,<br />
- produkty prognostyczne,<br />
- produkty groźnych zjawisk.<br />
Produkty podstawowe (Standard Products)<br />
Są to podstawowe produkty, wyliczane ze skanu klasycznego. Są wyrazem<br />
różnorodnych technik, metod wyliczania, prezentacji stanów atmosfery. Wyrażone są, dla<br />
większości produktów, w jednostkach odbiciowości radarowej [dBZ]. Wszystkie produkty tej<br />
grupy są obrazem odbiciowości radarowej lub są z niej pozyskiwane. Obrazują stan atmosfery<br />
na różnych wysokościach lub w różnych jej warstwach z obszaru znajdującego się w zasięgu<br />
pomiaru radarowego. Punktem odniesienia jest średni poziom morza.<br />
Produkty hydrologiczne (Hydrological Products)<br />
W ramach tej grupy generowane są produkty, które obrazują opad wyrażony w [mm]<br />
przypadających na jednostkę czasu np. [mm/godz.] lub w [mm] dla produktów wyliczających<br />
sumy opadów. Wartości tych produktów podawane są w skali logarytmicznej w celu lepszego<br />
zobrazowania mniejszych wartości opadu.<br />
Należy jednak zaznaczyć, że radar nie mierzy bezpośrednio natężenia opadu lecz<br />
wielkości z nim stowarzyszone – odbiciowość radarową, odbiciowość różnicową itp..<br />
Kwestia przeliczenia jednej wielkości na drugą stanowi cały odrębny dział wiedzy<br />
z dziedziny meteorologii radarowej.<br />
Produkty wiatrowe (Wind Products)<br />
Z zebranych w trakcie skanowania dopplerowskiego danych generowane są różnego<br />
rodzaju mapy prezentujące: rozkład prędkości radialnej, wiatru w postaci tradycyjnych<br />
strzałek wiatru lub wykresów, uskoków wiatru w pionie i poziomie oraz oszacowania<br />
turbulencji w obiektach meteorologicznych. Do ich wygenerowania niezbędne jest posiadanie<br />
dobrej jakości danych o prędkości.<br />
Należy w tym momencie zaznaczyć, że między wiatrem dopplerowskim, a ruchem pola<br />
opadu jest istotna różnica. Zmierzony wiatr dopplerowski jest wiatrem rzeczywistym -<br />
kropelki są unoszone bez opóźnień. Natomiast przesunięcie obliczono np. według korelacji<br />
czasowej pól opadu dotyczy całego układu opadowego, który ma swoją wewnętrzną<br />
dynamikę.<br />
13
Należy pamiętać, aby wszystkie produkty wiatrowe, do wyliczania których<br />
wykorzystywana jest składowa radialna wiatru, analizować łącznie z innymi<br />
produktami radarowymi jak i innymi źródłami informacji meteorologicznej np. z mapą<br />
synoptyczną, dla uniknięcia, możliwej, fałszywej interpretacji.<br />
Produkty prognoz i groźnych zjawisk (Forecasting & Phenomena Detection Products)<br />
W przykadku tej grupy produktów analizie poddawane są zarówno dane<br />
o odbiciowości jak i o prędkości radialnej. Wynikiem tego są produkty groźnych zjawisk typu<br />
opad i wiatr. Analizy wiatrowe prowadzone są w maksymalnym zasięgu do 125 km od radaru<br />
(ograniczenie wynikające ze skanu dopplerowskiego) co skutkuje takim samym<br />
ograniczeniem dla analiz odbiciowości jeśli chodzi o zasięg produktów groźnych zjawisk.<br />
Prawie wszystkie produkty dotyczące groźnych zjawisk mogą tworzyć z produktami<br />
meteorologicznymi, hydrologicznymi i wiatrowymi różne kombinacje. Produkty te<br />
w połączeniu np. z potencjalną wodnością, gradem, turbulencją oraz strefami uskoku pola<br />
wiatru mogą uzupełniać informację odnośnie genezy i tendencji rozwoju zjawisk<br />
meteorologicznych.<br />
Pozyskane w procesie przetwarzania informacje meteorologiczne i hydrologiczne<br />
prezentowane są na podkładach mapy fizycznej przeglądanego przez dany radar obszaru,<br />
opisane przez legendę zawierającą skalę kolorów z przypisanymi jej przedziałami wartości<br />
i podstawowymi informacjami dotyczącymi zasięgu oraz rozdzielczości prezentowanego<br />
obszaru. Na każdą mapę można nanieść dodatkowe informacje dotyczące prezentowanego<br />
terenu.<br />
Podstawowy format wyjściowy produktów <strong>radarowych</strong> bazuje na XML (Extensible<br />
Markup Language), powszechnie uznanym standardzie umożliwiający wymianę danych<br />
między różnymi systemami. Możliwe są także inne formaty np. HDF5 (format przyjęty<br />
obecnie za standard wymiany danych w Europie). Wyświetlane na ekranie monitora produkty<br />
wykonane są w grafice rastrowej 8- bitowej (na 1 piksel przypada 256 kolorów), a w pliku<br />
zapisywane są najczęściej w formacie PNG (Portable Network Graphics), choć dostępne są<br />
też inne formaty.<br />
Wszystkie produkty są dostępne dla pojedynczego radaru, a część z nich składana jest<br />
w mapy zbiorcze prezentujące obszar całej Polski lub regionu z zasięgu pomiaru wybranych<br />
stacji <strong>radarowych</strong>.<br />
Każde zobrazowanie radarowe produktów zawiera: mapę danych, podstawowe<br />
informacje dotyczące techniki skanowania, wartości niektórych parametrów służących do<br />
14
odczytu informacji zawartej w mapie oraz paletę barw z dowiązanymi wartościami. Z prawej<br />
strony produktów zamieszczona jest legenda opisująca produkt. Możliwe jest także<br />
generowanie zobrazowań bez legendy. Jednak wówczas taki produkt zawierający tylko mapę<br />
dla potencjalnego użytkownika byłby mniej czytelny. Mapy takie generowane są do<br />
zastosowań specjalistycznych.<br />
W polu legendy każdego z produktów, cztery pierwsze parametry zazwyczaj powtarzają się.<br />
Dlatego nie zawsze będą powtarzane przy opisie kolejnych produktów.<br />
Przykładowy produkt.<br />
PDF File - nazwa katalogu systemowego, w którym składowany jest produkt, a<br />
jednocześnie nazwa systemowa produktu (nadawana w procesie zapisu<br />
definicji produktu, właściwa tylko dla danego produktu).<br />
Clutter Filter - określa jakiego filtru dopplerowskiego użyto do eliminacji ech anomalnych,<br />
Time sampling - czas próbkowania (dla ustalonych wartości PRF oraz prędkości anteny<br />
procesor sygnału przetwarza pożądaną ilość impulsów przypadającą na jeden<br />
promień). Parametr dobierany w procesie definiowania skanu. Jest to ilość<br />
próbek w azymucie do uśredniania sygnału (czas jaki radar poświęca na<br />
próbkowanie komórek w przestrzeni w kierunku azymutalnym np. jeśli<br />
antena obraca się wolno to w tym samym czasie więcej impulsów wyśle w<br />
danym kierunku, zwiększy się częstość pozyskiwania próbek, tym samym<br />
dokładność pomiaru). Im większa liczba, tym lepiej uśredniane dane.<br />
PRF - częstotliwość powtarzania impulsów.<br />
Range - zasięg produktu (zobrazowany).<br />
15
Resolution - rozdzielczość produktu (rozmiar jednego piksela na mapie dla wyliczonej<br />
wartości produktu w tym punkcie) wyznaczana z zależności :<br />
Rozdzielczość = 2* Range/Size<br />
gdzie: Size jest parametrem ustalanym w procesie definiowania produktu.<br />
Height - wysokość w atmosferze, dla której policzono produkt.<br />
Alg type - typ zastosowanego algorytmu przyjęty do wyliczania produktu<br />
CAPPI Range - zasięg danych wyliczanych (uśrednianych) wg. algorytmu CAPPI.<br />
Data - rodzaj danych (źródło ich pochodzenia).<br />
Produkty bardziej specjalistyczne zawierające dodatkowe parametry, właściwe tylko dla<br />
produktu danego typu, będą dodatkowo opisywane.<br />
Większość produktów <strong>radarowych</strong> generowana jest w systemie pracy operacyjnej.<br />
Są to produkty niezbędne do wyrażenia aktualnego stanu atmosfery.<br />
W ramach systemu mogą być wytwarzane także produkty dodatkowe, standardowo nie<br />
generowane, a wykorzystywane do uzupełnienia lub wzbogacenia pozyskanej informacji.<br />
Poniżej zaprezentowano zestawienie omawianych w opracowaniu produktów z podziałem na<br />
produkty standardowo generowane i możliwe dodatkowo do wygenerowania.<br />
16
Nazwa<br />
skrócona/ nazwa<br />
angielska<br />
PPI(dBZ)<br />
Plan Position Indicator<br />
Polska nazwa produktu<br />
Produkty podstawowe (Standard products)<br />
Przekrój stożkowy – rozkład<br />
odbiciowości wzdłuż<br />
określonej elewacji.<br />
Źródło danych/Parametr istotny<br />
dla utworzenia produktu<br />
Skan klasyczny: Z[dBZ]<br />
Elevation - kąt położenia anteny<br />
(względem poziomu).<br />
PCAPPI(dBZ)<br />
Pseudo Constant<br />
Altitude PPI<br />
LMR<br />
Layer Mean Reflectivity<br />
RHI(dBZ)<br />
Range Height Indicator<br />
- reflectivity<br />
RHI(V)<br />
Range Height Indicator<br />
- radial velocity<br />
RHI(W)<br />
Range Height Indicator<br />
- spectral width<br />
MAX(dBZ)<br />
Maximum Display<br />
CMAX(dBZ)<br />
Maximum Display<br />
EHT(height)<br />
Echo Top<br />
EHT(height)<br />
Height of Max.<br />
Reflectivity<br />
EHT(height)<br />
Echo Base<br />
EHT(height)<br />
Layer Thickness<br />
Przekrój poziomy - rozkład<br />
odbiciowości na określonej<br />
wysokości.<br />
Odbiciowość średnia w<br />
zadanej warstwie.<br />
Przekrój pionowy na zadanym<br />
azymucie po odbiciowości<br />
radarowej.<br />
Przekrój pionowy na zadanym<br />
azymucie po rozkładzie<br />
prędkości radialnych.<br />
Przekrój pionowy na zadanym<br />
azymucie po rozkładzie widma<br />
prędkości radialnych.<br />
Maksymalna odbiciowość w<br />
rzucie na trzy płaszczyzny.<br />
Rozkład maksymalnej<br />
odbiciowości w rzucie na<br />
jedną płaszczyznę<br />
Wysokość echa radarowego<br />
Wysokość maksymalnej<br />
odbiciowości.<br />
Wysokość podstawy echa.<br />
Grubość struktury opadowej.<br />
Skan klasyczny: Z[dBZ]<br />
Height - wysokość n.p.m., dla<br />
której tworzony jest obraz.<br />
Skan klasyczny: Z[dBZ]<br />
Top, Bottom - górna i dolna<br />
granica danych (grubość<br />
warstwy).<br />
Skan pionowy: Z[dBZ]<br />
Azimuth - azymut, wzdłuż<br />
którego wykonywany jest<br />
przekrój.<br />
Skan pionowy: V[m/s]<br />
Azimuth – azymut, wzdłuż<br />
którego wykonywany jest<br />
przekrój.<br />
Skan pionowy: W[m/s]<br />
Azimuth - azymut, wzdłuż<br />
którego wykonywany jest<br />
przekrój.<br />
Skan klasyczny: Z[dBZ]<br />
Top, Bottom - górna i dolna<br />
granica danych.<br />
Skan klasyczny i dopplerowski:<br />
Z[dBZ]<br />
Top, Bottom - górna i dolna<br />
granica danych.<br />
Skan klasyczny: Z[dBZ]<br />
Skan klasyczny: Z[dBZ]<br />
Skan klasyczny: Z[dBZ]<br />
Skan klasyczny: Z[dBZ]<br />
17
Produkty hydrologiczne (Hydrological Products)<br />
SRI (dBR)<br />
Surface Rainfall<br />
Intensity<br />
PAC(dBA))<br />
(Precipitation<br />
Accumulation)<br />
VIL(dBA)<br />
Vertically Integrated<br />
Liguid Water<br />
Natężenie opadu na określonej<br />
wysokości.<br />
Suma opadu za: godzinę<br />
Suma opadu za 6 godzin<br />
Dobowa suma opadu<br />
Wodność scałkowana w pionie<br />
Skan klasyczny: Z[dBZ]<br />
Height Above Ground -<br />
wysokość ponad gruntem.<br />
Produkt drugiego rzędu bazujący<br />
na SRI(dBR).<br />
Skan klasyczny: Z[dBZ]<br />
Top, Bottom - górna i dolna<br />
granica danych (grubość<br />
warstwy).<br />
RSA<br />
River Subcatchment<br />
Accumulation<br />
<br />
Akumulacja opadu w<br />
podzlewniach<br />
Produkty wiatrowe (Wind Products)<br />
Produkt drugiego rzędu bazujący<br />
na produktach skanu<br />
klasycznego.<br />
Produktem wejściowym może<br />
być: CAPPI(dBZ), SRI(dBR) lub<br />
PAC(dBA).<br />
PPI(V)<br />
Plan Position Indicator -<br />
radial velocity<br />
PPI(W) -<br />
Plan Position Indicator -<br />
spectral width<br />
PCAPPI(V)<br />
Pseudo Constant<br />
Altitude PPI(V)<br />
PCAPPI(W)<br />
Pseudo Constant<br />
Altitude PPI(W)<br />
Przekrój stożkowy rozkładu<br />
prędkości radialnej<br />
Przekrój stożkowy obrazujący<br />
szerokości widma prędkości<br />
radialnych<br />
Rozkład prędkości radialnej na<br />
stałej wysokości typu PCAPPI<br />
Rozkład szerokości widma<br />
prędkości radialnych na stałej<br />
wysokości typu PCAPPI<br />
Skan dopplerowski: V[m/s]<br />
Elevation - kąt położenia anteny<br />
(względem poziomu).<br />
Skan dopplerowski: W[m/s]<br />
Elevation - kąt położenia anteny<br />
(względem poziomu)<br />
Skan dopplerowski: V[m/s]<br />
Height - wysokość n.p.m. dla<br />
której tworzony jest obraz<br />
Skan dopplerowski: W[m/s]<br />
Height - wysokość n.p.m.<br />
CMAX(V)<br />
Radial Velocity<br />
Maximum<br />
Rozkład maksymalnych<br />
prędkości radialnych w rzucie<br />
na jedną płaszczyznę<br />
Skan dopplerowski: V[m/s]<br />
CMAX(W)<br />
Spectral Width<br />
Maximum<br />
Rozkład maksymalnych<br />
wartości szerokości widma<br />
prędkości radialnych w rzucie<br />
na jedna płaszczyznę<br />
Skan dopplerowski: W[m/s]<br />
18
MAX(V)<br />
Radial Velocity<br />
Maximum<br />
Rozkład maksymalnych<br />
wartości prędkości radialnych<br />
zobrazowany na trzech<br />
płaszczyznach<br />
Skan dopplerowski: V[m/s]<br />
MAX(W)<br />
Spectral Width<br />
Maximum<br />
Rozkład maksymalnych<br />
wartości szerokości widma<br />
prędkości radialnych<br />
zobrazowany na trzech<br />
płaszczyznach<br />
Skan dopplerowski: W[m/s]<br />
VVP<br />
Volume Velocity<br />
Processing<br />
Pionowy profil wiatru<br />
Skan dopplerowski: V[m/s]<br />
HWIND<br />
Horizontal Wind<br />
Technique<br />
<br />
Rozkład wiatru poziomego na<br />
określonej wysokości<br />
Skan dopplerowski: V[m/s]<br />
Level - wysokość warstwy<br />
Number Of Grid Cells - ustala<br />
wielkości obszaru, na który<br />
przypada jeden wektor wiatru<br />
UWT<br />
Uniform Wind<br />
Technique<br />
<br />
Technika wiatru jednorodnego<br />
- wiatr poziomy na określonej<br />
wysokości<br />
Skan dopplerowski: V[m/s]<br />
VAD<br />
Velocity Azimuth<br />
Display<br />
Prędkość wiatru radialnego w<br />
funkcji azymutu<br />
Skan dopplerowski: V[m/s]<br />
SHEAR<br />
Radial Shear<br />
SHEAR<br />
Azimuthal Shear<br />
SHEAR<br />
Elevation Shear<br />
Radialny gradient wiatru<br />
Azymutalny gradient wiatru<br />
Elewacyjny gradient wiatru<br />
Skan dopplerowski: V[m/s]<br />
Elevation - kąt położenia anteny<br />
(względem poziomu), z którego<br />
pobierane są dane do wyliczenia<br />
produktu<br />
VSHEAR<br />
Vertical Shear<br />
Pionowy uskok wiatru.<br />
Skan dopplerowski: V[m/s]<br />
Top, Bottom -- górna i dolna<br />
granica danych (grubość<br />
warstwy)<br />
HSHEAR<br />
Horizontal Shear<br />
Poziomy uskok wiatru.<br />
Skan dopplerowski: V[m/s]<br />
Level - wysokość warstwy<br />
n.p.m.<br />
LTB<br />
Layer Turbulence<br />
Turbulencje w określonej<br />
warstwie atmosfery.<br />
Skan dopplerowski: W[m/s]<br />
Top i Bottom - wysokość górnej<br />
i dolnej warstwy, z której<br />
liczony jest produkt.<br />
19
Produkty groźnych zjawisk (Phenomena Detection Products)<br />
SWI<br />
Severe Weather Indicato<br />
<br />
ZHAIL<br />
Zhail detection<br />
Wskaźnik groźnych zjawisk<br />
Prawdopodobieństwo<br />
wystąpienia gradu.<br />
Skan dopplerowski: V[m/s],<br />
W[m/s], Z[dBZ],<br />
Skan klasyczny: Z[dBZ]<br />
Wymaga określenia izotermy<br />
Produkt prognostyczne (Nowcasting Products)<br />
RTR<br />
Rain Tracking<br />
CTR<br />
Cell Centroid Tracking<br />
<br />
Śledzenie pól opadowych.<br />
Śledzenie komórek<br />
burzowych.<br />
Produkt drugiego rzędu bazujący<br />
na produktach skanu<br />
klasycznego po Z[dBZ]<br />
i dopplerowskiego po V[m/s].<br />
Produkt drugiego rzędu bazujący<br />
na produktach skanu<br />
klasycznego po Z[dBZ]<br />
i dopplerowskiego po V[m/s].<br />
Produkty rozpoznania typu opadu (Dual-pol Products)<br />
CAPPI(ET)<br />
Echo and Hydrometeor<br />
Classification<br />
RHI(ET)<br />
Echo and Hydrometeor<br />
Classification<br />
Rozpoznanie typu opadu.<br />
Rozpoznanie typu opadu dla<br />
przekroju pionowego.<br />
Skan klasyczny: Z[dBZ],<br />
ZDR[dB], PhiDP [stopnie], KDP<br />
[stopnie/km], PhoHV [-].<br />
Skan pionowy: Z[dBZ],<br />
ZDR[dB], PhiDP [stopnie], KDP<br />
[stopnie/km], PhoHV [-].<br />
Zestawienie generowanych i możliwych do wygenerowania produktów <strong>radarowych</strong><br />
prezentowanych w niniejszym opracowaniu (stan na marzec 2011 r.)<br />
- kolorem żółtym oznaczono produkty, które są generowane w trybie pracy<br />
operacyjnej,<br />
Symbolem - oznaczono produkty nakładkowe.<br />
20
2. Podstawowe produkty radarowe<br />
2.1 Przekrój stożkowy – PPI (Plan Position Indicator)<br />
Produkt PPI(dBZ) jest odwzorowaniem rozkładu odbiciowości radarowej w<br />
zlokalizowanej strukturze meteorologicznej, pozyskanej z jednego obrotu anteny, wyrażonej<br />
w jednostkach odbiciowości radarowej [dBZ].<br />
Przykładowa mapa produktu<br />
PPI(dBZ).<br />
Elevation - kąt położenia anteny (względem poziomu), dla prezentowanej mapy równy 0,5 0<br />
Aktualnie w systemie POLRAD generowane są cztery produkty typu PPI dla odbiciowości<br />
radarowej. Pochodzą one z czterech<br />
pierwszych kątów podniesienia anteny:<br />
- PPI_1 dla 0,5 stopnia,<br />
- PPI_2 dla 1,4 stopnia,<br />
- PPI_3 dla 2,4 stopnia,<br />
- PPI_4 dla 3,4 stopnia.<br />
Wytwarzane są one na potrzeby systemu<br />
analityczno-prognostycznego NIMROD.<br />
Stanowią źródło danych <strong>radarowych</strong> na<br />
wejściu do tego systemu.<br />
Geometria tworzenia produktu PPI<br />
21
Podsumowanie<br />
PPI(dBZ) jest najszybciej tworzonym produktem, ponieważ generowany jest z<br />
danych dostarczonych tylko przez jeden obrót anteny.<br />
W zależności od odległości od radaru, wyświetlane dane pochodzą z różnych<br />
wysokości, co powoduje problemy z interpretacją danych.<br />
W pobliżu radaru występują echa stałe.<br />
Przydatność<br />
produktu<br />
Klasyczny produkt PPI (dla niskich kątów elewacji) wykonywany jest w celu<br />
uzyskania ogólnego przeglądu sytuacji.<br />
Bardzo przydatny w razie konieczności zwiększenia częstotliwości<br />
wykonywania przeglądu stanu atmosfery, szczególnie w przypadkach<br />
występowania zjawisk ekstremalnych. Zwiększenie ilości pozyskiwanych w<br />
jednostce czasu informacji pozwala na bardziej dokładną ocenę rozwoju<br />
sytuacji i lepsze monitorowanie atmosfery w niewielkich odstępach<br />
czasowych.<br />
Wymaga to jednak zmiany harmonogramu operacyjnego. W przypadku<br />
posiadania radaru mobilnego byłaby możliwość konfigurowania zadań<br />
dodatkowych w przypadku wystąpienia groźnych zdarzeń meteorologicznych<br />
np. w okresie zagrożenia powodziowego.<br />
Wówczas możliwe byłoby pozyskiwanie informacji radarowej np. co 5 min,<br />
lub częściej, o sytuacji nad danym radarem.<br />
Dane takie można wówczas wykorzystać do badań nad strukturą zjawisk jak<br />
i na potrzeby klientów strategicznych.<br />
22
2.2 Przekroje poziome – CAPPI(dBZ) i PCAPPI(dBZ) (Pseudo Constant Altitude<br />
Plan Position Indicator)<br />
Produkty CAPPI(dBZ) i PCAPPI(dBZ) są rozkładem odbiciowości radarowej,<br />
w zlokalizowanej strukturze meteorologicznej, wyliczonej dla określonej wysokości n.p.m..<br />
Przykładowa mapa produktu<br />
PCAPPI(dBZ)<br />
Level - wysokość warstwy CAPPI lub PCAPPI n.p.m., dla której wyliczany jest produkt<br />
Alg. Type - do liczenia produktu wykorzystano algorytm typu PCAPPI<br />
CAPPI Range - zakres pomiaru typu CAPPI<br />
Geometria tworzenia produktu PCAPPI.<br />
Operacyjnie generowanym produktem w ramach tej podgrupy jest produkt<br />
PCAPPI(dBZ) – odbiciowość radarowa na zadanej wysokości. Produktu ten zawiera rozkłady<br />
średnich ważonych mocy echa odbitego (odbiciowości) od obiektów na wysokości określonej<br />
ponad średnim poziomem morza, w zasięgu określonym przez parametry definiowalne<br />
produktu. Do wyliczenia produktów PCAPPI uwzględniane jest całe terytorium objęte<br />
skanem.<br />
23
Produkt PCAPPI(dBZ) to podstawowy produkt zobrazowania stanu atmosfery na ustalonej<br />
wysokości. Jest wynikiem przetworzenia informacji o odbiciowości radarowej. Na produkt<br />
składają się wartości odbiciowości wyliczone dla określonej wysokości nad średnim<br />
poziomem morza wyrażone w jednostkach [dBZ]. Typowa rozdzielczość danych to 1x1 km,<br />
zaś obszar prezentowanych (przez pojedynczy radar) danych to 400x400km. W legendzie<br />
każdego produktu podawane są granice zasięgu CAPPI.<br />
Mapa odbiciowości radarowej dla<br />
wysokości 1km policzona wg.<br />
algorytmu typu CAPPI, z terminu<br />
04.06.2007, godz. 16:10<br />
Podsumowanie<br />
Na mapie produktu PCAPPI(dBZ) mamy dane z całego obszaru w zasięgu<br />
skanowania klasycznego lub dopplerowskiego.<br />
Ograniczony zasięg dla trybu CAPPI – brak danych wokół radaru i w dalszej<br />
odległości zasięgu.<br />
W przypadku produktu PCAPPI dane spoza obszaru pomiaru CAPPI<br />
pochodzą z: wyższych warstw atmosfery (dla dalszych odległości) lub<br />
niższych (w pobliżu radaru). Jest to wynikiem właściwości pomiaru<br />
radarowego.<br />
Przydatność<br />
produktu<br />
PCAPPI to podstawowy produkt radarowy:<br />
- najlepiej odzwierciedlający stan atmosfery na zadanej wysokości,<br />
- wyjściowy do tworzenia wielu innych produktów typu opadowego.<br />
W związku z powyższym bardzo przydatny zarówno do badań jak<br />
i wykorzystania przez potencjalnego odbiorcę; do zobrazowań stanu atmosfery<br />
na różnych wysokościach.<br />
24
2.3 Odbiciowość średnia w zadanej warstwie – LMR (Layer Mean Reflectivity)<br />
Produkt LMR obrazuje średnie wartości odbiciowości radarowej dla dowolnie<br />
określonej warstwy atmosfery wyrażone w jednostkach odbiciowości radarowej [dBZ].<br />
Przykład produktu typu LMR.<br />
Bottom - dolna wysokość n.p.m. brana do wyliczenia produktu<br />
Top - górna wysokość n.p.m. brana do wyliczenia produktu<br />
Aktualnie produkty typu LMR nie są generowane operacyjnie z uwagi na brak<br />
zapotrzebowania.<br />
Podsumowanie<br />
Brak ech stałych na produkcie w przypadku gdy dolna granica zostanie<br />
ustalona wyżej niż wysokość instalacji radaru, w przeciwnym wypadku<br />
pojawiają się echa stałe.<br />
Redukuje efekt silnych zmian w pionowym profilu odbiciowości szczególnie<br />
w przypadku struktur opadowych składających się z mieszaniny<br />
hydrometeorów typu: śnieg/mokry śnieg/woda, grad/mokry grad/woda<br />
Ograniczony zasięg produktu. Produkt nie jest polecany do wykorzystania<br />
jako jedyny lecz jako pomocniczy wraz z grupą innych produktów.<br />
Przydatność<br />
produktu<br />
Użyteczny produkt do identyfikacji pionowej struktury atmosfery, a zwłaszcza<br />
komórek burzowych i struktury składającej się z mieszaniny hydrometeorów,<br />
o ile generowany jest jako zestaw map odpowiadających kolejnym,<br />
stosunkowo cienkim warstwom atmosfery.<br />
Przydatny szczególnie do analizy skali zjawiska w interesującej nas warstwie<br />
atmosfery np. dla lotu samolotu na poziomie określonego korytarza<br />
lotniczego.<br />
25
2.4 Odbiciowość maksymalna – MAX(dBZ) i CMAX(dBZ) (Maximum Display)<br />
Produkty MAX(dBZ) i CMAX(dBZ) są maksymalnymi wartościami odbiciowości<br />
radarowej, pochodzącymi z obszaru zlokalizowanej struktury meteorologicznej,<br />
zobrazowanymi w postaci rzutu na trzy lub jedną płaszczyznę.<br />
Przykładowa mapa obrazująca<br />
rozkład ech <strong>radarowych</strong> na<br />
produkcie MAX dla sytuacji<br />
konwekcyjnej (charakterystyczna<br />
dla okresu letniego).<br />
Height - wartość dolnej i górnej wysokości n.p.m., w zakresie których wyliczano produkt<br />
Hor.Res. - rozdzielczość produktu w poziomie dla obrazu „widok z góry”<br />
Vert.Res. - rozdzielczość obrazu w pionie dla zobrazowań bocznych<br />
Obraz MAX pochodzi jakby z „prześwietlenia” chmury w każdym kierunku i obrazuje<br />
maksymalną wartość odbiciowości na każdym z kierunków. W tym celu dla każdego piksela<br />
na mapie pobierana jest, z dostępnych danych, pionowa lub pozioma kolumna, a z niej<br />
wybierana największa wartość. W efekcie tworzony jest potrójny obraz analizowanej sytuacji,<br />
na który składają się:<br />
● widok najwyższych zmierzonych wartości odbiciowości w rzucie<br />
pionowym z góry na powierzchnię ziemi - czyli mapa odbiciowości<br />
maksymalnej dla każdej pionowej kolumny (pole A na rysunku),<br />
● widok największych zmierzonych wartości odbiciowości w kierunku<br />
północ-południe przedstawiony jako pasek górny produktu obrazuje<br />
odbiciowości maksymalne dla każdej linii prostej skierowanej z<br />
południa na północ (pole B na rysunku),<br />
● widok największych zmierzonych wartości odbiciowości w kierunku wschód-zachód<br />
przedstawiony jako zobrazowanie po prawej stronie mapy pokazuje odbiciowości<br />
maksymalne wzdłuż linii prostych skierowanych z zachodu na wschód (pole C).<br />
26
Produkt CMAX(dBZ)<br />
Produkt typu CMAX jest obrazem maksymalnych wartości pomiarów <strong>radarowych</strong><br />
generowany tylko dla płaszczyzny poziomej (na rysunku pomocniczym pole A), bez<br />
zobrazowań bocznych. Dla każdej z kolumn o wyżej wymienionej podstawie mogą być<br />
wyświetlane największe wartości: odbiciowości, prędkości radialnej lub rozkładu widma<br />
prędkości radialnych w całej przeglądanej przez radar przestrzeni (patrz rozdz. 4.4).<br />
Produkt CMAX - rozkład<br />
maksymalnych wartość<br />
odbiciowości utworzony dla<br />
zakresu wysokości od 1 do 15 km.<br />
Podsumowanie<br />
Zobrazowanie stanu atmosfery w rzucie na trzy płaszczyzny daje wrażenie<br />
trójwymiarowego spojrzenia na sytuację meteorologiczną i pozwala na<br />
zobrazowanie rozkładu ech <strong>radarowych</strong> w atmosferze. Oddaje obraz typu<br />
struktury opadowej (warstwowa, konwekcyjna czy mieszana).<br />
Gdy dolną granicę danych ustalimy odpowiednio wyżej niż wysokość anteny<br />
to do analizy nie będą brane echa stałe.<br />
Produkt mniej użyteczny do analizy danych o prędkości; im dalej od radaru<br />
tym wyżej nad horyzontem biegnie wiązka, tym większe wartości prędkości.<br />
Nie polecany do określania obszarów opadu.<br />
Przydatność<br />
produktu<br />
Bardzo użyteczny produkt do szybkiej oceny sytuacji meteorologicznej oraz<br />
analizy odbiciowości radarowej, także dla dalszych odległości od radaru w<br />
przypadku identyfikacji komórek typu cumulnimbus. Istotny, gdy występuje<br />
konieczność szybkiego zorientowania się w strukturze pionowej i obszarach<br />
groźnych zjawisk zachodzących w zidentyfikowanym obiekcie<br />
meteorologicznym.<br />
Polecany do wykorzystania jako produkt uzupełniający obraz sytuacji<br />
meteorologicznej, szczególnie w połączeniu z produktem SRI (natężenie<br />
opadu). Pokazuje cały obszar atmosfery, w obrębie którego zidentyfikowano<br />
echo radarowe bez względu na to czy jest ono źródłem opadu czy nie.<br />
27
2.5 Wysokość echa radarowego - EHT (Echo Height)<br />
Na podstawie analizy danych o odbiciowości możliwe są do wygenerowania<br />
następujące wartości:<br />
- wysokości wierzchołków echa radarowego – EHT Echo Top ,<br />
- wysokości echa radarowego maksymalnej wartości odbiciowości lub natężenia opadu –<br />
EHT Height Max Z ,<br />
- wysokości podstawy echa EHT Echo Base (dla przestrzeni objętej pomiarem),<br />
- grubość warstwy echa EHT Thicknes (różnica pomiędzy najwyższą, a najniższą wysokością<br />
echa w danym punkcie pomiaru) co odpowiada obliczeniom: Echo Top minus Echo Base.<br />
2.5.1 Wierzchołki echa radarowego EHT (Echo Top)<br />
Produkt EHT (Echo Top) jest obrazem wysokości wierzchołków ech <strong>radarowych</strong><br />
w obrębie zlokalizowanej struktury meteorologicznej, wyrażonych w [km].<br />
Przykładowa mapa produktu<br />
typu EHT.<br />
Min. Z - wartość progowa odbiciowości lub natężenia opadu brana do wyliczenia produktu<br />
Data - typ produktu, który może przyjmować następujące określenia:<br />
Echo Top - wysokości wierzchołków zlokalizowanych struktur meteorologicznych<br />
(wysokość wierzchołków ech)<br />
Echo Base - podstawa ech <strong>radarowych</strong> (wysokość podstawy echa)<br />
Height of Max. Reflectivity - wysokość maksymalnej odbiciowości<br />
Layer Thickness – grubości zlokalizowanej strukury.<br />
28
EHT Echo Top wytwarzany operacyjnie jest obrazem wysokości wierzchołków echa<br />
radarowego w zasięgu do 200 - 250 km od radaru, z rozdzielczością 1x1 km (punkt<br />
zobrazowania). W procesie generowania produktu analizowane są wszystkie miejsca gdzie<br />
wykryto echa radarowe, a zobrazowane wartości pochodzą z największych wartości<br />
wysokości, na których je zlokalizowano.<br />
2.5.2 Wysokość maksymalnej odbiciowości - EHT Height Max Z<br />
Produkt EHT(Height MaxZ) jest obrazem wysokości, na których zlokalizowano<br />
największą wartość odbiciowości radarowej, wyrażonych w [km].<br />
Mapa EHT - Height<br />
Max Z<br />
dla sytuacji z godziny<br />
15:50 dnia 29 lipca<br />
2005, radar Legionowo.<br />
Dla poszczególnych punktów mapy prezentuje on wartości wysokości, dla których<br />
wyliczono największe wartości odbiciowości w kolumnach o podstawie 1x1km od<br />
powierzchni ziemi do maksymalnej wysokości możliwego pomiaru.<br />
2.5.3 Wysokości podstawy echa radarowego – EHT Echo Base<br />
Produkt EHT(EchoBase) jest obrazem podstawy echa radarowego wyrażonej w [km].<br />
29
Mapa EHT - Echo Base<br />
dla sytuacji z dnia 29 lipca<br />
2005 godz. 15:50 - radar<br />
Legionowo.<br />
Podsumowanie<br />
Mapa wysokości wierzchołków obiektów meteorologicznych to bardzo<br />
„praktyczny” produkt. Pokazuje ona, dla zdefiniowanego zasięgu, wszystkie<br />
miejsca gdzie wykryto echo radarowe.<br />
Przydatny do automatycznej identyfikacji echa i jego pionowego rozwoju –<br />
maksymalnej rozpiętości w pionie, niezbędnej do oceny sytuacji w przypadku<br />
rozróżnienia występowania struktur konwekcyjnych i struktur warstwowych.<br />
W przypadku opadów o typowych „zwartych” strukturach, rozbudowanych w<br />
pionie, na produkcie typu Echo Top występuje „koncentryczna struktura<br />
kołowa”, która daje niejednoznaczny obraz danego obiektu<br />
meteorologicznego.<br />
Złagodzenie „koncentrycznej struktury kołowej”, na mapie produktu, można<br />
zniwelować poprzez zwiększenie liczby kątów elewacji składających się na<br />
cykl obserwacji. Nieznacznie efekt ten można poprawić (tylko dla<br />
wizualizacji produktu) poprzez zwiększenie zakresu skali kolorów dla<br />
prezentowanych wysokości. Nie wyeliminuje to jednak braków tego produktu<br />
wynikających z techniki pomiaru radarowego (tylko 10 kątów podniesienia<br />
anteny).<br />
Mała przydatność operacyjna produktu EHT Echo Base i EHT Thickness<br />
Przydatność<br />
produktu<br />
Bardzo użyteczny produkt do identyfikacji pionowego zasięgu struktury<br />
opadowej, a zwłaszcza komórek burzowych. Przy jego pomocy można<br />
określić wartości wysokości wierzchołków echa radarowego, z czego<br />
możemy wnioskować o występowaniu niebezpiecznych zjawisk w danym<br />
obszarze.<br />
W lotnictwie bardzo przydatny do wyznaczania obszarów lotu, na których<br />
możliwe jest ominięcie atmosfery zajętej przez rozbudowane struktury<br />
opadowe.<br />
30
2.6 Przekrój pionowy - RHI (Range Height Indicator)<br />
Produkt RHI(dBZ) jest rozkładem pionowym odbiciowości radarowej pozyskanej<br />
w procesie skanowania pionowego na ściśle określonym azymucie.<br />
Przykładowe produkty typu RHI.<br />
Range - zasięg produktu<br />
Vert Res - pionowa rozdzielczość obrazu<br />
Hor Res - pozioma rozdzielczość obrazu<br />
Height - wysokość skanowania pionowego (wysokość powyżej średniego poziomu morza)<br />
Elevation - rozpiętość sektora pionowego skanowania atmosfery<br />
Podsumowanie<br />
Produkt RHI charakteryzuje się wysoką rozdzielczością danych, zależną od<br />
ustawień parametrów skanowania pionowego.<br />
Produkt RHI(V) jest przydatny do badań związanych ze szczegółową analizą<br />
danej struktury meteorologicznej.<br />
Przekroje typu RHI są możliwe jedynie wzdłuż promienia od radaru. Chociaż<br />
możliwe jest ustalenie dowolnego azymutu, to musi on być zadany w trakcie<br />
definiowania harmonogramu skanowania pionowego, co wyklucza<br />
dynamiczne dostosowywanie się do aktualnej sytuacji meteorologicznej.<br />
Obecnie RHI jest w mniejszym, lub większym stopniu produktem<br />
wychodzącym z użycia i zastępowanym przez produkty VCUT i LMVCUT,<br />
które mogą być na bieżąco generowane w trakcie pracy operacyjnej. z danych<br />
typu 3D.<br />
Przydatność<br />
produktu<br />
Produkt przydatny (zalecany) w przypadku konieczności monitoringu stanu<br />
atmosfery na wybranym azymucie w stosunku do lokalizacji radaru np.<br />
w strefach podejścia samolotów do lądowania oraz możliwy do<br />
wykorzystania w przypadku obserwacji i analiz szczególnie niebezpiecznych<br />
zjawisk meteorologicznych w celach badawczych.<br />
31
2.7 Przekrój pionowy przez dwa punkty - VCUT (Vertical Cut) i przekrój wzdłuż<br />
łamanej MLVCUT (Multiple-Line Vertical Cut)<br />
Produkty VCUT(dBZ) prezentują rozkład pionowy odbiciowości radarowej dla<br />
przekroju wyznaczonego wzdłuż zadanego odcinka, poprowadzonego na dowolnie<br />
wygenerowanym produkcie. Przy czym produkt VCUT będzie tego samego typu co produkt<br />
wejściowy.<br />
Przykładowy obraz przekroju<br />
pionowego wzdłuż zadanego<br />
odcinka po odbiciowości<br />
radarowej.<br />
Każdy wykonany przekrój<br />
pionowy może być w systemie<br />
zapisany w oddzielnym pliku<br />
jako niezależny produkt.<br />
Height – przedział wysokości do zobrazowania danych (parametr ustawiany w konfiguracji<br />
wyświetlania produktu)<br />
Disp Len – długość trasy wyznaczonej do wykonania przekroju zobrazowania; wartość<br />
wyznaczona przez punk start i stop<br />
Hor Res - pozioma rozdzielczość obrazu<br />
Ver Res - pionowa rozdzielczość obrazu<br />
Przekroje wzdłuż łamanej MLVCUT (Multiple-Line Vertical Cut)<br />
Produkty MLVCUT(dBZ) prezentują rozkład pionowy odbiciowości radarowej<br />
przekroju wyznaczonego wzdłuż zadanej linii łamanej, poprowadzonej na wygenerowanym<br />
produkcie obrazującym daną strukturę meteorologiczną.<br />
32
Obraz rozkładu ech <strong>radarowych</strong> na<br />
produkcie MLVCUT wzdłuż łamanej<br />
składającej się z kilku odcinków.<br />
Sumaryczna długość przekroju przez<br />
strukturę to 357 km.<br />
Podsumowanie<br />
Przekrój typu VCUT i MLVCUT można wykonywać wzdłuż dowolnego<br />
odcinka lub łamanej. Początek i koniec cięcia ustalany jest w trybie<br />
interaktywnej pracy z produktem.<br />
Pozwala na pozyskanie orientacji o rozkładzie pionowym struktur opadowych<br />
nad wybraną lokalizacją.<br />
Wygodniej jest posługiwać się produktem typu VCUT (w stosunku do RHI)<br />
gdyż pozwala dowolnie orientować cięcia pionowe, bez konieczności<br />
wykonywania dodatkowego skanu.<br />
Wykonany obraz przekroju pionowego można także zapisać w postaci pliku<br />
źródłowego jak i graficznego.<br />
Gorsza rozdzielczość pionowa niż w przypadku produktu RHI.<br />
Produkty możliwe do wygenerowania tylko w przypadku dostępu do<br />
specjalistycznego oprogramowania i do danych w układzie sferycznym (typu<br />
3D).<br />
Przydatność<br />
produktu<br />
Bardzo użyteczny produkt, zwłaszcza przy analizie pionowego rozwoju<br />
dowolnie wybranej struktury meteorologicznej, zjawiska meteorologicznego<br />
(nie tylko ekstremalnego).<br />
33
3. Produkty hydrologiczne<br />
3.1 Natężenie opadu - SRI (Surface Rainfall Intensity)<br />
Produkt SRI jest obrazem rozkładu natężenia opadu w [mm/h] na zdefiniowanej przez<br />
użytkownika wysokości nad powierzchnią gruntu. Wysokość powierzchni ziemi obliczana<br />
jest z map orograficznych, które są także wykorzystywane do określania regionów gdzie<br />
ustawiona przez użytkownika warstwa jest niedostępna dla radaru.<br />
Dane produktu SRI są prezentowane w<br />
skali logarytmicznej, a ich wartości są<br />
wyskalowane w jednostkach natężenia<br />
opadu [mm/h].<br />
1 [mm/h] = 1 [litr/m 2 ] w ciągu godziny<br />
dBZ 10log10<br />
Z<br />
SRI H - wysokość nad powierzchnią terenu<br />
ZR - wartości współczynników dla relacji: odbiciowość radarowa (Z) - opad (R):<br />
a - czynnik A relacji Z-R przyjęty do przeliczania odbiciowości na natężenie opadu<br />
b - wykładnik B relacji Z-R przyjęty do przeliczania odbiciowości na natężenie<br />
opadu<br />
Alg. Type - typ algorytmu w oparciu, o który wyliczone są wartości produktu poza obszarem<br />
uwzględniającym orografię terenu:<br />
SRI - dane wyliczane w oparciu o algorytm typu SRI (analogicznie jak CAPPI),<br />
Pseudo SRI - dane wyliczane w oparciu o algorytm PseudoSRI (analogicznie jak<br />
PCAPPI).<br />
Standardowy produkt SRI jest obrazem rozkładu natężenia opadu na zdefiniowanej<br />
przez użytkownika wysokości nad powierzchnią terenu w obszarze 200 km wokół radaru<br />
zgodnie z algorytmem PseudoSRI.<br />
W procesie wyliczania produktu, pozyskane przez radar informacje o odbiciowości radarowej<br />
w [dBZ] poddawane są konwersji, z uwzględnieniem empirycznie wyznaczonych parametrów<br />
a i b, do natężenia opadu R według opracowanej przez Marshala-Palmer’a zależności:<br />
Z a R<br />
b<br />
34
gdzie:<br />
Z - jest w [mm 6 /m 3 ], R - w [mm/h],<br />
a, b – to wartości wyznaczone empirycznie.<br />
Założeniem produktów hydrologicznych jest liczenie opadu w warstwie<br />
uwzględniającej orografię terenu tzn. wartość prezentowana prze każdy punkt mapy winna<br />
pochodzić z stałej wysokości nad poziomem terenu określonej parametrem SRI H.<br />
Geometria tworzenia produktu SRI<br />
(Źródło: dokumentacja systemu<br />
RAIMBOW).<br />
Wysokości te są obliczane z map orografii terenu. Jeśli w definicji produktu nie jest możliwe<br />
uwzględnienie mapy orografii terenu to wówczas produkt SRI liczony jest na stałej wysokości<br />
nad średnim poziomem morza, zgodnie z zasadami przyjętymi przez algorytm PCAPPI(dBZ).<br />
Podsumowanie<br />
W prostej formie określa oszacowany w oparciu o dane radarowe rozkład<br />
opadu nad danym obszarem, wyrażony w skali barw, którym<br />
przyporządkowano określone wartości natężenia opadu.<br />
Używany jako zestaw danych wejściowych do generowania wielu innych<br />
produktów hydrologicznych.<br />
Powinien informować o opadzie tuż przy powierzchni ziemi, jednak im<br />
większa odległość od radaru tym większy błąd jego oszacowania. Wynika to z<br />
geometrii skanowania i specyfiki pomiaru radarowego.<br />
Przydatność<br />
produktu<br />
Bardzo przydatny produkt, który w łatwej do interpretacji formie (bez<br />
konieczności wnikania w zagadnienia konwersji odbiciowości radarowej)<br />
podaje rozkład natężenia opadu. W związku z tym może posługiwać się nim<br />
każdy użytkownik.<br />
Analiza map z kilku kolejnych SRI może być pomocna do określenia<br />
tendencji (spadek/wzrost) opadu na najbliższe dwie, trzy godziny.<br />
35
3.2 Suma opadów - PAC (Precipitation Accumulation)<br />
Produkt PAC jest zsumowanym natężeniem opadu w zdefiniowanym okresie czasu,<br />
podany w [mm]. Najczęściej generowanym produktem w oparciu o pozyskane dane radarowe<br />
jest godzinna suma opadu odświeżana co 10 minut.<br />
1 [mm] opadu oznacza:<br />
1 litr wody na powierzchnię 1 [m 2 ]<br />
Przykładowy obraz rozkładu<br />
godzinnej sumy opadu.<br />
Start - termin rozpoczęcia sumowania<br />
Interval - przedział czasu, dla którego wykonywano sumowanie opadu – łączny czas<br />
sumowania w: dniach (d), godzinach (h) lub minutach (m)<br />
Num Prod. - ilość produktów brana do analizy w celu wyliczenia PAC<br />
Miss Time – okres braku pomiarów<br />
W ramach systemu radarowego liczone są operacyjnie następujące sumy opadów:- godzinna<br />
suma opadów, 6-cio godzinne sumy opadów, dobowa suma opadów.<br />
Podsumowanie<br />
Sumy opadu są obliczane na podstawie kolejnych 10-cio minutowych<br />
pomiarów <strong>radarowych</strong>. Wartości te nadają się do porównania z wartościami<br />
sum opadu zmierzonych za dany okres czasu przez deszczomierze.<br />
Przydatność<br />
produktu<br />
Z uwagi na nieciągły charakter danych wejściowych (co 10 min.) zachodzi<br />
potrzeba interpolacji czasowej.<br />
Jako danych wejściowych potrzebuje także produktu typu prognostycznego -<br />
RTR (Śledzenie Komórek Opadowych).<br />
Produkt bardzo użyteczny dla zastosowań hydrologicznych.<br />
Przydatny do określania rozkładu opadu na danym obszarze w kolejnych<br />
jednostkach czasu.<br />
Szczególnie przydatny w momencie śledzenia rozwoju sytuacji opadowej w<br />
przypadkach występowania powodzi.<br />
Pomocny przy szacowaniu wezbrań na zlewniach.<br />
36
3.3 Wodność scałkowana w pionie - VIL (Vertical Integrated Liquid)<br />
Produkt VIL określa chwilową, oszacowaną zawartości wody w zdefiniowanej przez<br />
użytkownika warstwie atmosfery, wyrażoną w jednostkach sum opadu [mm]. Obraz produktu<br />
jest jej odzwierciedleniem w każdym punkcie wyznaczonym przez rozdzielczość produktu.<br />
Przykładowy produkt<br />
typu VIL<br />
Height – przedział wysokości, z którego wyliczono produkt<br />
ZM-C, ZM-D - wartości współczynników do przeliczenia odbiciowości na wodność:<br />
C - współczynnik C relacji Z-M<br />
D - wykładnik D relacji Z-M<br />
Podsumowanie<br />
Umożliwia policzenie wodności w dowolnie zadanych przedziałach<br />
wysokości co można wykorzystać do jej określenia na zadanych wysokościach<br />
np. tras korytarzy lotniczych (np. od 7 do 9 km) w celu wydzielenia komórek<br />
konwekcyjnych na tle rozbudowanych struktur opadowych.<br />
Produkt zaniża wartości w pobliżu stacji radarowej i w dalszych odległościach<br />
od radaru (na skrajach zasięgu).<br />
Przydatność<br />
produktu<br />
Szczególnie użyteczny do wyznaczenia obszarów o zawartości dużej masy<br />
wody nad danym obszarem. Pozwala na oszacowanie opadu jaki spadnie w<br />
najbliższym czasie (za kilka minut).<br />
Szczególnie przydatny do zobrazowania obszarów zagrożonych wystąpieniem<br />
opadów ekstremalnych pochodzących z pionowo rozbudowanych struktur<br />
opadowych; silnych burz.<br />
37
3.4 Akumulacja opadu w podzlewniach - RSA (River Subcatchment Accumulation)<br />
RSA to produkt umożliwiający oszacowanie sum opadu na obszarze zdefiniowanej<br />
przez użytkownika zlewni dla zdefiniowanego przedziału czasu.<br />
Przykładowy produkt obrazujący obszar zlewni z wylistowanymi wartościami<br />
wygenerowanymi przez RSA dla danego terminu obserwacji.<br />
Alg type - typ algorytmu użyty do obliczania sumy opadu<br />
Start - termin początkowy sumowania opadu (miesiąc, dzień, godz., min., sek., rok)<br />
Interval - przedział czasowy sumowania opadu określony w: dniach, godzinach i minutach<br />
Numprod - liczba produktów wejściowych w przedziale sumy<br />
Miss Time - łączny czas brakujących danych w przedziale czasu sumy (dz., godz, min.)<br />
RSA dostarcza następujących informacji dla zdefiniowanych przez użytkownika obszarów<br />
(najczęściej podzlewni rzek):<br />
- dla każdego kroku czasowego produktu wejściowego, obliczane są uśrednione przestrzennie<br />
natężenie opadu R i odpływ Q,<br />
- dla przedziału czasowego zdefiniowanego przez użytkownika odpływ całkowity Q tot ,<br />
odpływ uśredniany po czasie Q avg , opad całkowity R tot i natężenie opadu uśredniane po<br />
czasie R avg .<br />
W momencie listowania wyliczonych wartości dla każdej ze zdefiniowanych zlewni w<br />
otwierającym się oknie dostajemy:<br />
Name - nazwę zlewni,<br />
Size - powierzchnię zlewni w [km 2 ],<br />
Total Rain<br />
Average Rain<br />
- całkowity opad R tot w [mm],<br />
- średnie natężenie opadu R avg w [mm/h],<br />
38
Total Runoff - całkowity odpływ Q tot w [m 3 ],<br />
Average Runoff - średni napływ Q avg w [m 3 /s].<br />
Dane dla tego samego<br />
terminu wyświetlone w<br />
układzie X-Y<br />
Wyświetlone są na nim<br />
szeregi czasowe odpływu Q<br />
odnośnie zlewni o nazwie:<br />
Sub_1<br />
Druga zlewnia zdefiniowana<br />
w tym produkcie to Sub_2<br />
Natomiast gdy w preferencjach produktu wybrana zostanie opcja wyświetlania<br />
produktu jako „diagram”, szeregi czasowe odpływu Q są pokazywane jako dwuwymiarowe<br />
diagramy w układzie X-Y dla każdej z podzlewni (wyświetlane w momencie rozwinięcia<br />
zakładki Layer na produkcie głównym).<br />
Podsumowanie<br />
W oparciu o ten produkt istnieje możliwość policzenia całkowitego zasilenia<br />
zlewni przez opad w funkcji czasu.<br />
Konieczność zdefiniowania zlewni.<br />
Przydatność<br />
produktu<br />
Wartości produktu zapisane w plikach ASCII, mogą zostać użyte jako dane<br />
wejściowe do modeli hydrologicznych.<br />
Mogą być wykorzystane do różnego rodzaju badań i analiz np. związanych z<br />
porównaniem ilości opadu (wody) pochodzącego z różnych struktur<br />
meteorologicznych, względem wezbrań na zlewniach określonych rzek.<br />
39
4. Produkty wiatrowe<br />
W trybie pracy operacyjnej do pozyskania danych do wyliczenia produktów<br />
wiatrowych stosowane jest specjalistyczne, odrębne skanowanie atmosfery nazywane skanem<br />
dopplerowskim (patrz rozdział 1.3). Jak już wcześniej opisano jest to skanowanie<br />
zoptymalizowane pod kątem pomiaru danych o prędkości, zorganizowane tak, aby otrzymać<br />
wyniki o jak najwyższej jakości.<br />
4.1 Produkty wiatrowe według geometrii ich tworzenia<br />
4.1.1 Przekroje stożkowe rozkładu prędkości radialnej i szerokości widma prędkości<br />
radialnych<br />
Przekrój stożkowy rozkładu prędkości radialnej PPI(V) - (Plan Position Indicator)<br />
Produkt PPI(V) jest odwzorowaniem rozkładu prędkości radialnych pozyskanych<br />
z pomiaru dopplerowskiego dla jednego obrotu anteny, w pełnym zakresie azymutu i na<br />
zadanej elewacji. Prędkości prezentowane są w [m/s].<br />
Przykładowa mapa rozkładu<br />
prędkości radialnych typu<br />
PPI(V) dla kąta elewacji<br />
równego 0,8 ; zobrazowanie dla<br />
standardowo wyświetlanego<br />
zakresu prędkości.<br />
Elevation – kąt podniesienia anteny<br />
Produkty PPI(V) mogą być pozyskiwane z pomiaru prędkości radialnych dla wszystkich<br />
elewacji (kątów podniesienia anteny) skanu dopplerowskiego (patrz rozdział 1.2). Nazwa<br />
została przyjęta ze względu na analogię do geometrii tworzenia produktu PPI(dBZ) (patrz<br />
rozdział 2.1).<br />
Przekrój stożkowy obrazujący szerokości widma prędkości radialnych - PPI(W) (Plan<br />
Position Indicator).<br />
40
Produkt PPI(W) jest odwzorowaniem szerokości widma prędkości radialnych<br />
pozyskanych z pomiaru dopplerowskiego dla jednego obrotu anteny, w pełnym zakresie<br />
azymutu i zadanej elewacji. Szerokości widma prędkości radialnych prezentowane są w [m/s].<br />
Zobrazowanie rozkładu szerokości<br />
widma prędkości radialnych<br />
PPI(W)dla elewacji równej 0,8<br />
w zakresie wartości :<br />
od 0,1 m/s do 10 m/s.<br />
Elevation – kąt podniesienia anteny<br />
Produkty PPI(W) mogą być pozyskiwane z pomiaru szerokości widma prędkości radialnych<br />
dla wszystkich elewacji (kątów podniesienia anteny) skanu dopplerowskiego (patrz rozdział<br />
1.2). Nazwa została przyjęta ze względu analogię do geometrii tworzenia produktu PPI(dBZ)<br />
(patrz rozdział 2.1).<br />
Podsumowanie<br />
PPI(V) i PPI(W) są najszybciej tworzonymi produktami typu wiatrowego<br />
ponieważ mogą być generowane z danych pochodzących z pojedynczego<br />
obrotu anteny.<br />
PPI(V) przydatny do pośredniego określania kierunku wiatru wewnątrz<br />
struktury meteorologicznej, a PPI(W) do szacowania turbulencji.<br />
Wyświetlane dane pochodzą z różnych wysokości w zależności od odległości<br />
od radaru, co powoduje problemy z interpretacją danych.<br />
Dla najwyższych kątów elewacji, wyświetlany jest właściwie pionowy profil<br />
prędkości i szerokości widma prędkości radialnych.<br />
Przydatność<br />
produktu<br />
Produkt PPI(V) i PPI(W) - dla niskich kątów elewacji - wykonywany jest w<br />
celu uzyskania ogólnego, szybkiego przeglądu sytuacji w zakresie<br />
prezentowanych przez nie wartości. Ze względu na trudności interpretacyjne<br />
przydatny raczej do systemów automatycznych, działających bez ingerencji<br />
człowieka.<br />
41
4.1.2 Rozkład prędkości radialnej i rozkład szerokości widma prędkości radialnych na<br />
stałej wysokości typu PCAPPI (Pseudo Constant Altitude Plan Position<br />
Indicator)<br />
Produkt PCAPPI(V) jest odwzorowaniem rozkładu prędkości radialnych na określonej<br />
wysokości n.p.m., wyliczonych w oparciu o dane pozyskane z całego skanu dopplerowskiego.<br />
Rozkład prędkości radialnych prezentowany jest w [m/s].<br />
a) b)<br />
Rozkłady prędkości radialnej dla wysokości: a) 1km i b) 2 km z tego samego terminu.<br />
Produkt PCAPPI(V) jest jednym z podstawowych produktów zobrazowania składowej<br />
radialnej wiatru w zlokalizowanych echach <strong>radarowych</strong>, względem stacji radarowej.<br />
Produkt PCAPPI(W) jest odwzorowaniem rozkładu szerokości widma prędkości<br />
radialnych na określonej wysokości n.p.m. Rozkład ten prezentowany jest w [m/s].<br />
Produkty rozkładu szerokości widma prędkości radialnych PCAPPI(W) dla wysokości<br />
1 km i 2 km z tego samego terminu obserwacji.<br />
Height - wysokość warstwy CAPPI n.p.m.<br />
Alg Type - typ zastosowanego algorytmu, opcjonalnie do wyboru :<br />
CAPPI - wg. algorytmu typu CAPPI<br />
Pseudo CAPPI - wg. algorytmu typu PCAPPI<br />
CAPPI Range - zakres CAPPI; podany jest zasięg CAPPI dla danej wysokości<br />
42
Podsumowanie<br />
W porównaniu do PPI(V) wyeliminowany jest wpływ skrętu wiatru z<br />
wysokością, na powstający obraz prędkości.<br />
Jako samodzielny produkt jest on trudny do interpretacji, szczególnie w<br />
przypadku występowania zjawisk ekstremalnych typu wiatrowego.<br />
Wymagana jest praktyka i doświadczenie, aby poprawnie go zinterpretować.<br />
Przydatność<br />
produktu<br />
Rozkład prędkości radialnej typu PCAPPI(V) wykorzystywany jest do<br />
identyfikowania zjawisk typu: wiry i trąby powietrzne, zjawisk dywergencji<br />
i konwergencji w obszarze danej struktury.<br />
Produkt CAPPI(W) może być przydatny do określania turbulencji w danej<br />
warstwie atmosfery.<br />
4.1.3 Rozkład prędkości i szerokości widma prędkości radialnych na produktach typu<br />
CMAX i MAX<br />
Produkty rozkładu prędkości radialnych i szerokości widma prędkości radialnych<br />
mogą być także generowane według geometrii produktów CMAX i MAX. Produktami tymi<br />
są wówczas niżej wymienione mapy.<br />
CMAX(V) - rozkład maksymalnych wartości prędkości radialnej dla określonego przedziału<br />
wysokości.<br />
CMAX(W) - rozkład maksymalnych wartości szerokości widma prędkości radialnych dla<br />
określonego przedziału wysokości.<br />
Rozkład maksymalnych wartości prędkości radialnych i spektrum rozkładu prędkości<br />
radialnych wygenerowane według geometrii tworzenia produktu CMAX.<br />
43
MAX(V) - rozkład maksymalnych wartości prędkości radialnej dla określonego przedziału<br />
wysokości w rzucie na trzy płaszczyzny.<br />
MAX(W) - rozkład maksymalnych wartości szerokości widma prędkości radialnych w rzucie<br />
na trzy płaszczyzny<br />
Warunkiem utworzenia opisanych w tym podrozdziale produktów jest ich zdefiniowanie w<br />
systemie radarowym. Operacyjnie produkty te nie są generowane.<br />
Rozkład maksymalnych wartości prędkości radialnych i maksymalnych wartości widma<br />
prędkości radialnych w zlokalizowanej strukturze.<br />
4.1.4 Rozkład prędkości radialnych i szerokości widma prędkości radialnych na<br />
produktach typu VCUT<br />
Rozkłady danych o prędkości radialnej i szerokości widma prędkości radialnych mogą<br />
być także tworzone wg. geometrii produktów typu VCUT lub MLVCUT. Produktami tymi są<br />
wówczas prezentowane niżej mapy.<br />
VCUT(V) - rozkład pionowy prędkości radialnych wzdłuż określonego (wyznaczonego<br />
odcinka) przekroju przez strukturę meteorologiczną.<br />
VCUT(W) – rozkład pionowy szerokości widma prędkości radialnych wzdłuż określonego<br />
(wyznaczonego odcinka) przekroju przez strukturę meteorologiczną.<br />
44
Przykładowe produkty:<br />
b) rozkładu prędkości radialnych,<br />
c) szerokości widma prędkości radialnych<br />
na przekrojach pionowych typu VCUT<br />
wyznaczonych wzdłuż odcinka<br />
poprowadzonego na produkcie rozkładu<br />
odbiciowości (rysunek a).<br />
a)<br />
b) c)<br />
MLVCUT(V) - rozkład pionowy prędkości radialnych wzdłuż określonego (wyznaczonego<br />
przez linie łamaną) przekroju przez strukturę meteorologiczną.<br />
MLVCUT(W) – rozkład pionowy szerokości widma prędkości radialnych wzdłuż<br />
określonego (wyznaczonego przez linię łamaną) przekroju przez strukturę<br />
meteorologiczną.<br />
Przekrój MLVCUT(W) w obrazie 3D wzdłuż łamanej wyznaczonej na produkcie CMAX(W)<br />
45
Przekrój MLVCUT(V)<br />
oraz jego obraz 3D<br />
wzdłuż łamanej<br />
wyznaczonej na<br />
produkcie CMAX(V).<br />
Warunkiem utworzenia wymienionych w tym podrozdziale produktów jest dostęp do danych<br />
pochodzących ze skanu dopplerowskiego zapisanych w układzie sferycznym.<br />
Podsumowanie<br />
Zobrazowanie stanu atmosfery w rzucie na trzy płaszczyzny w przypadku<br />
MAX(V) ułatwia interpretację danych na obrazie głównym, który bywa mało<br />
czytelny z uwagi na „przemieszanie wysokości”, z których pochodzi wartość<br />
naniesiona na mapę.<br />
Produkt typu MAX(W) – pozwala na zlokalizowanie obszarów o dużych<br />
wartościach szerokości widma prędkości radialnych, a VCUT(W) na dokładne<br />
ich zobrazowanie i przeanalizowanie w połączeniu z lokalizacjami<br />
(obszarami) o silnej turbulencji.<br />
W przypadku produktu MAX(V) na płaszczyźnie poziomej produktu<br />
zobrazowane będą wartości z różnych, choć najczęściej górnych warstw<br />
atmosfery.<br />
Nie polecany do określania obszarów występowania maksymalnych prędkości<br />
radialnych wiatru.<br />
Należy pamiętać, że obszar maksymalnej prędkości radialnej wcale nie musi<br />
odpowiadać obszarowi maksymalnej prędkości wiatru.<br />
Przydatność<br />
produktu<br />
Produkty wykorzystywane do analiz zdarzeń ekstremalnych typu<br />
meteorologicznego oraz ekspertyz wykonywanych w przypadku katastrof<br />
lotniczych, w celu wykluczenia lub potwierdzenia wpływu na nie czynników<br />
meteorologicznych.<br />
46
4.2 Pionowy profil wiatru – VVP (Volume Velocity Processing)<br />
Produkt VVP obrazuje rozkład wektorów wiatru wyliczonych do określonej<br />
wysokości nad stacją radarową. Prędkość i kierunek wiatru liczone są dla jednakowo<br />
odległych poziomów wysokości. Uzyskujemy w ten sposób przybliżony profil pionowy<br />
wiatru dla momentu obserwacji, który może być wyświetlony w postaci: tradycyjnych<br />
wektorów wiatru lub wykresu gdzie: oś x to prędkość wiatru w [m/s], a oś y określa wysokość<br />
w [km].<br />
Dla każdego skanu dopplerowskiego wyliczone wartości<br />
VVP mogą być wyświetlane na oddzielnych<br />
diagramach/wykresach. W przypadku, gdy wybrano<br />
prezentację w postaci strzałek z upierzeniem, odchylenie<br />
strzałki względem góry odpowiada odchyleniu kierunku<br />
wiatru względem północy.<br />
Przykładowy produkt typu VVP.<br />
Range – podaje minimalny i maksymalny zasięg odległości od radaru dla obszaru, z którego<br />
pobierane są dane do wyliczenia produktu:<br />
minimum – granica odległości od stacji radarowej do wewnętrznego boku obszaru danych<br />
maksimum – granica odległości od stacji radarowej do zewnętrznego boku obszaru danych<br />
Alg type - typ pierwszej regresji używanej do wyznaczania VVP, która może być<br />
prowadzona w oparciu o analizę danych o prędkości radialnej z wykorzystaniem:<br />
47
trójwymiarowego wektora wiatru (typ - 3D Wind Vector) oraz dodatkowo z<br />
wykorzystaniem pochodnych poziomych (typ - 3D Vector and Hor. Derivatives)<br />
i pochodnych pionowych (typ - Complete).<br />
Elevation – zakres kątów podniesienia anteny w elewacji dla procesu skanowania<br />
dopplerowskiego.<br />
Second reg – może przyjmować wartość: On lub Off. On oznacza wyliczanie produktu<br />
z uwzględnieniem drugiej regresji.<br />
Uaktywnienie tego okna oznacza proces drugiej regresji (ponowne liczenie<br />
wektorów wiatru w oparciu o te same dane) wykonywany w celu dalszej eliminacji<br />
danych odstających od przyjętych wartości progowych. Jest to kontynuacja analizy<br />
danych pod kątem poprawy jakości produktu.<br />
Podsumowanie<br />
Profile pionowe generowane są w oparciu o analizę wielu parametrów<br />
wyliczanych z rozkładu prędkości radialnych. Co skutkuje dobrą jakością<br />
produktu.<br />
Radar podaje wiatr tam gdzie wykryje echa, ale czasami nawet w przypadku<br />
braku hydrometeorów (gdy nie ma opadów) algorytm jest w stanie<br />
wygenerować profile wiatru do wysokości kilku kilometrów. Dzieje się tak ze<br />
względu na wysoką czułość „kanału prędkości” radaru, który jest w stanie<br />
zarejestrować echo od zanieczyszczeń (pyły, sole) znajdujących się w<br />
powietrzu.<br />
Jako informacje wejściowe do wygenerowania tego produktu potrzebne są<br />
dane o prędkości, które powinny być dobrej jakości.<br />
Wymaga oddzielnego skanu, zoptymalizowany pod kątem pomiaru prędkości.<br />
Z przyczyn teoretycznych (skomplikowanych algorytmów) konieczne jest<br />
sprawdzanie jakości danych (pierwsza i druga regresja).<br />
Przydatność<br />
produktu<br />
Niesie bardzo dużo informacji na temat pionowego profilu wiatru w miejscu<br />
lokalizacji radaru.<br />
Bardzo użyteczny produkt dla lotnictwa oraz pracy synoptyka (jako<br />
uzupełnienie profili pionowych wiatru generowanych przez stacje<br />
aerologiczne).<br />
48
4.3 Wiatr poziomy – HWIND (Horizontal Wind) i Technika wiatru jednorodnego<br />
UWT – (Uniform Wind Technique)<br />
Produkty HWIND i UWT prezentują rozkład wiatru poziomego na określonej<br />
wysokości: jego kierunek oraz prędkość z obszarów gdzie na danej wysokości zlokalizowano<br />
jakiekolwiek echo radarowe. Pole wiatru poziomego na produkcie HWIND i UWT<br />
wyświetlane jest jako obraz poziomych wektorów wiatru na stałej wysokości w punktach<br />
regularnej siatki, opisanych przez długość i szerokość geograficzną.<br />
Produkt typu HWIND(V) z<br />
wygenerowanymi wektorami<br />
wiatru poziomego dla<br />
wysokości 2 km.<br />
Range – zasięg obliczeń<br />
Height – wysokość warstwy n.p.m., dla której wyznaczono wektory wiatru<br />
Alg. Type – wybrany typ algorytmu zastosowany do wyliczania produktu:<br />
HWIND - oznacza, że do wyznaczania wektorów wiatru, wartości rozkładu<br />
prędkości radialnych uśredniano (interpolowano) wg. algorytmu typu CAPPI,<br />
PHWIND - PseudoHWIND - oznacza, że do wyznaczania wektorów wiatru,<br />
wartości rozkładu prędkości radialnych interpolowano i pobierano wg. założeń<br />
algorytmu typu PCAPPI (patrz rozdział: 1.4).<br />
Poziome wektory wiatru są obliczane dla punktów siatki (bramek) wybranych z wnętrza<br />
obszaru pokrytego skanem dopplerowskim. W trakcie definiowania produktu wyznaczana jest<br />
liczba bramek siatki w kierunkach północ-południe i wschód-zachód określająca w pewien<br />
sposób „rozdzielczość produktu graficznego” tj. wielkość obszaru, na który przypada jeden<br />
graficznie zobrazowany wektor wiatru. Liczba wektorów wiatru (komórek siatki) jest tak<br />
dobierana, aby nie przeładować obrazu nadmiarem informacji.<br />
Do wyliczenia poziomych wektorów wiatru na określonej wysokości wykorzystywane są<br />
wartości rozkładu prędkości radialnej uśrednianej (interpolowanej) według różnych zasad.<br />
W przypadku produktu:<br />
49
- HWIND(V) zgodnie z przyjętym algorytmem typu CAPPI(V) lub PCAPPI(V),<br />
- UWT(V) z zastosowaniem sześciopunktowego schemat interpolacji liniowej.<br />
Dla produktu UWT, wykorzystywana jest uproszczona metodyka interpretacji pól<br />
poziomych dla zmierzonej prędkości dopplerowskiej. Założeniem tej metody jest stałość pola<br />
wiatru w analizowanym obszarze tj. poziome wektory wiatru liczone są, przy założeniu<br />
istnienia jednorodnego pola wiatru w określonym obszarze, dla niewysokich kątów<br />
nachylenia anteny (do 5 stopni).<br />
Wiatr poziomy policzony według algorytmu HWIND (po lewej) i UWT (po prawej) dla<br />
sytuacji występowania wiatru jednorodnego.<br />
Podsumowanie<br />
Pomiar radarem dopplerowskim umożliwia wyznaczenie wektorów wiatru<br />
poziomego w maksymalnym zasięgu do 125 km wokół każdej ze stacji<br />
<strong>radarowych</strong>. Produkty HWIND i UWT najlepiej je przedstawiają<br />
W celu pozyskania produktu niezbędne są dane o dobrej jakości. Z tego też<br />
powodu wymagany jest oddzielny skan, zoptymalizowany pod kątem pomiaru<br />
danych o rozkładzie prędkości radialnej.<br />
Ze względu na stopień skomplikowania algorytmów, konieczne jest<br />
sprawdzanie jakości danych w kolejnych krokach generowania produktów.<br />
Przydatność<br />
produktu<br />
Bardzo przydatny produkt do obrazowania pola wiatru zarówno rozległej<br />
struktury jak i pojedynczych komórek.<br />
Zarówno produkt HWIND jak i UWT mogą być kombinowane z wieloma<br />
innymi produktami. Ułatwia to ich wzajemną interpretację np. jeżeli<br />
połączymy je z odbiciowością radarową lub rozkładem prędkości radialnej,<br />
mamy połączoną informację o określonym obiekcie meteorologicznym.<br />
50
4.4 Prędkość wiatru radialnego w funkcji azymutu – VAD (Velocity Azimuth Display)<br />
Produkt VAD jest średnią prędkością radialną w funkcji kąta azymutu, dla ustalonego<br />
kąta elewacji w zakresie określonego przedziału odległości ustalonego wzdłuż wiązki<br />
radarowej.<br />
Zobrazowanie średnich prędkości radialnych prezentowane jest na produkcie w formie<br />
wykresu w układzie x/y, gdzie: oś x to azymut, a oś y to prędkość radialna.<br />
Standardowo wyświetlanym produktem jest zobrazowanie pochodzące z pierwszego kąta<br />
poniesienia anteny.<br />
Przykład produktu VAD. Dla tego przypadku analizy prowadzone są dla elewacji 0.5 o w<br />
zasięgu od 20 do 25 km (Min.Range=20 km, plus zasięg analiz danych ( pierścień=5 km).<br />
Rozkład prędkości radialnych (w pełnym zakresie azymutów) wykazuje kształt sinusoidalny.<br />
Rysunek po prawej stronie jest rysunkiem pomocniczym do zobrazowania obszaru analiz<br />
danych (pole zielonego pierścienia).<br />
Min. Range – minimalna odległość od radaru, dla której prowadzone są wyliczenia<br />
Minimalny zasięg VAD w [km]<br />
Range gap - rozmiar pierścienia (grubość pierścienia), w obrębie którego analizowane są<br />
dane w [km]<br />
Podsumowanie<br />
Szybkie, wizualne sprawdzanie danych o rozkładzie prędkości radialnej w<br />
określonym obszarze. Krótki czas obliczania.<br />
Interpretacja niejednorodnego wiatru wymaga sporego doświadczenia.<br />
Przydatność<br />
produktu<br />
Produkt do użytku powszechnego raczej nie stosowany. Wykorzystywany do<br />
analiz w przypadku konieczności przeprowadzenia specjalistycznych<br />
ekspertyz.<br />
51
4.5 Produkty uskoków wiatru typu SHEAR<br />
4.5.1 Produkty typu SHEAR<br />
Produkty typu SHEAR generowane są z danych o rozkładzie prędkości radialnej v r<br />
[m/s] (przekrojów stożkowych typu PPI(V)) i wyliczane są jako różnica prędkości radialnych<br />
między sąsiednimi punktami podzielona przez „odległość” między nimi:<br />
SHEAR<br />
v<br />
l<br />
2 1<br />
[m/s/km]<br />
2<br />
v<br />
l<br />
1<br />
Możliwe są do wygenerowania następujące niezależne obrazy składowych dla gradientu:<br />
radialny, azymutalny i elewacyjny. W zależności od wyboru może być obliczany<br />
i prezentowany obraz wartości tylko dla jednej składowej gradientu. Są one jakby produktami<br />
pośrednimi (półproduktami) obliczanymi w celu wyliczenia produktów SHEAR<br />
kombinowanych. Jednostką SHEAR jest [m/s/km].<br />
Radial SHEAR<br />
vr<br />
RS <br />
r<br />
Obliczana jest składowa radialna gradientu<br />
z prędkości radialnej v r , na określonej<br />
powierzchni elewacji PPI.<br />
Azimuth SHEAR<br />
v<br />
AS <br />
r<br />
r<br />
a<br />
Obliczana jest wartość gradientu z prędkości<br />
radialnej v r w kierunku zmiany azymutu w<br />
trakcie obrotu anteny.<br />
52
Elevation SHEAR<br />
v<br />
ES <br />
r<br />
r<br />
e<br />
Dla obliczenia pochodnej rozpatrywana jest<br />
wartość prędkości radialnej z elewacji bieżącej<br />
i poprzedniej - niższej (jeśli istnieje) oraz<br />
następnej - wyższej elewacji (jeśli istnieje).<br />
gdzie :<br />
r, r a , r e - odpowiednio współrzędne wzdłuż kierunku: radialnego, azymutalnego<br />
i elewacyjnego,<br />
v r<br />
- pochodna z prędkości radialnej policzona przy pomocy regresji liniowej dla kilku<br />
r ...<br />
pomiarów v r w obrębie przedziałów, których wielkość ustalamy wzdłuż kierunków:<br />
r, r a , lub r e .<br />
Podsumowanie<br />
Składowa radialna gradientu (radial shear) wykorzystywana jest do<br />
wykrywania obszarów zbieżności i rozbieżności.<br />
Składowa azymutalna gradientu (azimuth shear) jest dobrym wskaźnikiem<br />
rotacji: zgodnej z ruchem wskazówek zegara (wartość ujemna) lub przeciwnej<br />
do ruchu wskazówek zegara (wartość dodatnia).<br />
Produkty generowane są dla powierzchni stożkowych (PPI) mają więc<br />
wszystkie wady produktów opartych o ten algorytm (patrz rozdz. 4.1).<br />
Wymagają wysokiej jakości danych o prędkości radialnej. Zaburzenia danych<br />
o prędkości mogą powodować generowanie silnych wartości gradientów<br />
(generowanie fałszywych alarmów).<br />
Przydatność<br />
produktu<br />
Przydatny do detekcji drobnoskalowych obszarów zbieżności/rozbieżności<br />
rotacji: cyrkulacji cyklonalnych i antycyklonalnych. Pamiętajmy jednak, że są<br />
także produkty bardziej zaawansowane jak np. SWI do wykrywania tych<br />
obszarów. Generują one mniej fałszywych alarmów gdyż biorą pod uwagę<br />
więcej charakterystyk pola wiatru.<br />
53
4.5.2 Poziomy gradient wiatru – HSHEAR (Horizontal Shear)<br />
Produkt HSHEAR jest obrazem zmiany prędkości radialnej na płaszczyźnie poziomej<br />
w zlokalizowanej strukturze meteorologicznej, na zadanej wysokości uwzględniającej<br />
krzywiznę ziemi. Wartości uskoku poziomego podawane są w jednostkach [m/s/km].<br />
Przykładowy produkt typu<br />
HSHEAR wygenerowane dla<br />
wysokości 2 km<br />
Produkt HSHEAR wyprowadzany jest z rozkładu prędkości radialnych na stałej wysokości<br />
nad poziomem morza uśrednianych zgodnie z algorytmem typu CAPPI (patrz rozdział 1.4).<br />
Następnie obliczany jako pierwiastek z sumy kwadratów składowych gradientu w kierunku x<br />
i y:<br />
vr<br />
HSHEAR <br />
x<br />
2<br />
<br />
<br />
<br />
vr<br />
<br />
y<br />
2<br />
<br />
<br />
<br />
na podstawie oszacowania zmiany prędkości wiatru w kierunku północ - południe (w<br />
kierunku osi x) oraz wschód - zachód (w kierunku osi y) w warstwie poziomej CAPPI(V).<br />
Podsumowanie<br />
Wynik produktu jest zbliżony do kombinacji dwóch produktów: gradientu<br />
radialnego (radial shear) oraz gradientu azymutalnego (azimuth shear).<br />
Jednak analiza w tym przypadku wykonywana jest dla określonej warstwy<br />
o stałej wysokości nad poziomem morza co umożliwia wykorzystanie jej do<br />
wyliczenia poziomego uskoku wiatru.<br />
54
Wymagają wysokiej jakości danych o prędkości radialnej.<br />
Mały zasięg produktu dla niższych warstw atmosfery (poniżej 2 km).<br />
Zaburzenia danych o prędkości mogą powodować generowanie silnych<br />
gradientów prędkości skutkiem czego powstawać może niebezpieczeństwo<br />
fałszywych alarmów.<br />
Przydatność<br />
produktu<br />
Wykorzystywany jest do obserwacji horyzontalnego gradientu wiatru co<br />
oznacza, że umożliwia wykrywanie poziomego uskoku wiatru w warstwie na<br />
stałej wysokości nad poziomem morza.<br />
Produkt przydatny do obserwacji zjawisk ekstremalnych typu wiatrowego.<br />
Bardzo przydatny dla lotnictwa do wspomagania decyzji wymagających<br />
znajomości wartości poziomego uskoku wiatru na określonych wysokościach<br />
np. tras lotniczych.<br />
4.5.3 Pionowy uskok wiatru – VSHEAR (Vertical Shear)<br />
Produkt VSHEAR jest obrazem zmiany radialnej prędkości wiatru pomiędzy dwoma<br />
ściśle określonymi warstwami atmosfery. Jednostką uskoku pionowego jest [m/s/km].<br />
Przykład produktu<br />
VSHEAR – uskok<br />
pionowy wiatru<br />
określony pomiędzy<br />
warstwami 1-3km.<br />
VSHEAR <br />
v<br />
top<br />
<br />
v<br />
bottom<br />
odleglosc<br />
gdzie:<br />
v<br />
top<br />
- prędkości w górnej warstwie,<br />
v<br />
botom<br />
- prędkości w dolnej warstwie,<br />
odleglosc - odległość między warstwami.<br />
Wartości VSHEAR są obliczane jako absolutna różnica prędkości radialnych<br />
pomiędzy sąsiadującymi w pionie komórkami z dwu warstw, wyznaczonych według<br />
algorytmu typu CAPPI(V), podzielona przez odległość tych komórek.<br />
55
Podsumowanie<br />
Produkt obliczany na podstawie danych o rozkładzie prędkości radialnych na<br />
określonych wysokościach w celu umożliwienia określania pionowego uskoku<br />
wiatru dla sąsiadujących w pionie komórek.<br />
Gdy warstwa, w obrębie której jest wyliczany uskok pionowy jest dość<br />
„cienka” to tym produkt jest lepszej jakości (lepiej odpowiada temu czego po<br />
jego nazwie oczekujemy).<br />
Wymaga wysokiej jakości danych o prędkości radialnej.<br />
Zaburzenia danych o prędkości mogą powodować generowanie silnych<br />
gradientów prędkości skutkiem czego powstawać może niebezpieczeństwo<br />
fałszywych alarmów.<br />
Ograniczony zasięg produktu: dla niższych warstw atmosfery bardzo mały<br />
zasięg.<br />
Przydatność<br />
produktu<br />
Wykorzystywany do detekcji gwałtownych zmian wiatru pomiędzy<br />
określonymi warstwami atmosfery.<br />
Produkt przydatny do obserwacji zjawisk ekstremalnych typu wiatrowego.<br />
Bardzo przydatny dla lotnictwa do wspomagania decyzji wymagających<br />
znajomości wartości pionowego uskoku wiatru na określonych wysokościach<br />
np. trasach korytarzy lotniczych.<br />
4.6 Turbulencje – LTB (Layer Turbulence)<br />
Produkt LTB jest obrazem turbulencji pomiędzy dwoma określonymi warstwami<br />
atmosfery. Jego wartości wyrażone są w [m/s].<br />
Pole rozkładu wartości<br />
turbulencji dla warstwy<br />
od 1 do 3 km<br />
wyświetlone dla zakresu<br />
wartości od 0.1 do 6 m/s.<br />
Height - przedział wysokości, dla których prowadzona jest analiza danych i generowany jest<br />
produkt, określana jest przez wysokość dolnej i górnej warstwy n.p.m.<br />
56
Podsumowanie<br />
Produkt daje możliwość analizy danych o rozkładzie szerokości widmowej (W)<br />
pomiędzy dwoma dowolnie wybranymi warstwami atmosfery typu CAPPI (o<br />
stałej wysokości nad poziomem morza) dzięki czemu może być wykorzystany<br />
do obserwacji turbulencji pomiędzy nimi.<br />
Wymaga wysokiej jakości danych o rozkładzie szerokości widma prędkości<br />
radialnych.<br />
Przydatność<br />
produktu<br />
Szczególnie przydatny dla lotnictwa do detekcji turbulencji na różnych<br />
wysokościach np. tras lotniczych.<br />
Wykorzystywany do analizy zjawisk ekstremalnych w celu określenia<br />
intensywności ruchu hydrometeorów w zlokalizowanej strukturze<br />
meteorologicznej.<br />
57
5. Produkty analizy groźnych zjawisk<br />
5.1 Wskaźnik groźnych zjawisk – SWI (Severe Weather Indicator)<br />
Produkt SWI jest produktem analitycznym, gdzie analizie poddawane są<br />
trójwymiarowe dane radarowe, pozyskane ze skanowania dopplerowskiego celem wykrycia<br />
obszarów potencjalnie zagrożonych wystąpieniem groźnych zjawisk. Produkt wyznacza<br />
(sygnalizuje wystąpienie):<br />
• obszary intensywnych opadów i rdzeni struktur burzowych,<br />
• obszary gdzie możliwe jest zlokalizowanie mezocyklonow i mezo-antycyklonów,<br />
• obszary rozbieżności i zbieżności (dywergencji, konwergencji),<br />
• obszary zjawisk typu microburst.<br />
W momencie zidentyfikowania wymienionych zjawisk są one zobrazowane na produkcie<br />
SWI w odpowiednim obszarze, w postaci specjalnych znaczników.<br />
Dla lepszego<br />
wyodrębnienia<br />
struktur burzowych<br />
zaznaczono je<br />
kolorem bordowym<br />
(na produkcie<br />
generowanym<br />
operacyjnie jest to<br />
kolor czerwony).<br />
Przykładowy produkt SWI z dnia 20080815 godz. 14:54 z rozwiniętą legendą typu Layer w<br />
celu pokazania zjawisk wygenerowanych na produkcie SWI.<br />
58
Zidentyfikowane zjawiska wyświetlane są w postaci symboli, wycentrowanych na punktach,<br />
które podane są jako ich współrzędne geograficzne. Średnica symbolu odzwierciedla wielkość<br />
zjawiska. W trakcie operacyjnej pracy systemu bardziej szczegółowe informacje dotyczące<br />
charakterystyk groźnych zjawisk mogą być, w każdej chwili, wyświetlane na ekranie<br />
monitora w okienku pomocniczym.<br />
Zjawiska możliwe do zdiagnozowania przez produkt SWI pozyskiwane są w kilku<br />
krokach. Wykrywanie ich rozpoczyna się od rozpoznania obszarów występowania struktur<br />
burzowych, następnie mezocyklonów i obszarów zbieżności/rozbieżności. Końcowym<br />
etapem analiz jest rozpoznanie zjawiska typu mikroburst. Rozpoznanie tych zjawisk,<br />
wykonywane jest niezależnie i prezentowane może być jako oddzielne produkty. W systemie<br />
POLRAD są jednak operacyjnie obrazowane na jednym sumarycznym produkcie (SWI).<br />
Zjawiska rozpoznawane przez produkt SWI wyświetlone jako niezależne produkty:<br />
intensywny opad SSA (Storm Structure Analyses), dywergencja/konwergencja VERG<br />
(Con/di-vergence detection)<br />
Wykrywanie obszarów komórek konwekcyjnych z intensywnym opadem i rdzeni<br />
struktur burzowych.<br />
System opadowy<br />
Rdzeń opadowy<br />
Każda ze zidentyfikowanych struktur systemu opadowego,<br />
zaznaczona jest na mapie produktu ciągłą linią (czerwoną)<br />
i opisana w pliku produktu.<br />
Żółtą linią zakreślone są obszary o najintensywniejszych<br />
opadach jako rdzeń danej komórki burzowej.<br />
59
W celu zlokalizowania trójwymiarowych, konwekcyjnych komórek burzowych,<br />
wszystkie dane o odbiciowości, przekraczające zadany próg, są łączone w trójwymiarowe<br />
segmenty celem określenia obszarów występowania komórek konwekcyjnych i ich rdzeni<br />
(centrów komórek konwekcyjnych o dużej odbiciowości), które mogą być źródłem<br />
intensywnego opadu deszczu lub gradu.<br />
Analiza struktur mezo- i antycyklonalnych<br />
Mezocyklon<br />
Anty-mezocyklon<br />
Zidentyfikowane mezocyklony wyświetlane są na mapie<br />
produktu SWI w postaci symboli. Średnica symbolu<br />
odzwierciedla rozmiar MC. Dla wyliczonych charakterystyk<br />
możliwe jest wyświetlenie mezocylonów o rotacji<br />
cyklonalnej jak i antycyklonalnej.<br />
W przypadku produktu SWI algorytm lokalizujący obszary wystąpienia mezocyklonów<br />
analizuje dane wiatrowe zebrane z kolejnych elewacji, w trakcie skanowania<br />
dopplerowskiego, pod kątem wykrywania wirów. W przypadku zidentyfikowania wirów<br />
analizowane są sąsiednie warstwy w poszukiwaniu spełnienia warunków wg. określonych<br />
wzorców mezocyklonów. Następnie wykonywane jest porównanie czy zidentyfikowane wiry<br />
spełniają założenia odpowiadające warunkom mezocyklonu. Decyzja w jakim obszarze jest<br />
mezocyklon, podejmowana jest na podstawie wielu parametrów, które należy dobierać<br />
indywidualnie dla każdej lokalizacji radaru.<br />
Wykrywanie obszarów zbieżności i rozbieżności<br />
Konwergencja - zbieżność<br />
Dywergencja - rozbieżność<br />
W celu wykrycia zbieżności/rozbieżności analizowane są<br />
trójwymiarowe dane o prędkości radialnej w układzie<br />
współrzędnych sferycznych. W trakcie analiz poszukiwane<br />
są obszary pojawienia się struktur, spełniających założenia<br />
zastosowanych w algorytmie wzorców, określających<br />
wystąpienie rozbieżności lub zbieżności.<br />
60
Wykrywanie zjawisk typu mikroburst<br />
Mikroburst<br />
Prekursor Mikroburst<br />
Zjawisko mikroburst charakteryzuje się powstaniem na małym<br />
obszarze silnej rozbieżności wiatru blisko powierzchni ziemi,<br />
najczęściej związane z burzami. W trakcie jego „poszukiwania”<br />
prowadzona jest szczegółowa analiza danych <strong>radarowych</strong> pod<br />
kątem odbiciowości (poszukiwanie pewnych cech<br />
charakterystycznych dla burz), sprawdzane są dane o prędkości w<br />
poszukiwaniu oznak rozbieżności/zbieżności.<br />
Podsumowanie<br />
Automatyczne wykrywanie i analiza groźnych zjawisk. Produkt SWI może<br />
być wyświetlany jako nakładka łączona z innymi produktami.<br />
Prowadzenie analiz rozpoznawania zjawisk możliwe tylko w zasięgu do 125<br />
km od radaru.<br />
Jakość danych o prędkości radialnej powinna być możliwie wysoka.<br />
Zaburzenia danych o prędkości skutkują generowaniem silnych gradientów<br />
prędkości, a w konsekwencji generowanie rozpoznawania zjawisk, które<br />
w rzeczywistości nie miały miejsca (fałszywe alarmy).<br />
Przydatność<br />
produktu<br />
Użyteczny produkt do zautomatyzowanego wykrywania i analizy groźnych<br />
zjawisk zarówno typu opadowego jak i wiatrowego.<br />
Bardzo przydatny do badań zjawisk ekstremalnych (wykonywania ich<br />
charakterystyk) oraz badań związanych z cyklem życia komórek burzowych.<br />
61
5.2 Prawdopodobieństwo wystąpienia gradu – ZHAIL (Hail Detection)<br />
Produkt ZHAIL określa procentowe prawdopodobieństwo wystąpienia gradu, na danym<br />
obszarze, bazując na danych o odbiciowości (pozyskanych z radarów pojedynczej polaryzacji<br />
fali) i danych synoptycznych.<br />
Przykład produktu obrazującego<br />
prawdopodobieństwo wystąpienia<br />
gradu ZHAIL<br />
MinZ - minimalna odbiciowość w [dBZ], przy której może pojawić się grad<br />
H0 - wartość izotermy zero w [km]<br />
H Hail - minimalna wielkość warstwy w [km] pomiędzy wysokością izotermy zero, a<br />
wysokością wierzchołków ech określonych wg. zdefiniowanego progu odbiciowości<br />
Podsumowanie<br />
Umożliwia wyznaczanie prawdopodobieństwa wystąpienia gradu nad<br />
określonym obszarem.<br />
Standardowy algorytm do wykrywania gradu analizuje jedynie maksymalną<br />
odbiciowość w kolumnie pionowej, co daje wysokie prawdopodobieństwo<br />
fałszywych alarmów. W przypadku produktu ZHAIL brane są pod uwagę<br />
zarówno wysokość warstwy zamarzania jak i grubość warstwy silnej<br />
odbiciowości ponad nią, które to czynniki mają istotny wpływ na algorytm.<br />
Przewartościowanie wysokości ech zredukowano poprzez rozważanie efektu<br />
rozszerzania się wiązki radarowej wraz ze wzrostem odległości od radaru.<br />
Konieczność wprowadzania dodatkowych parametrów meteorologicznych, nie<br />
zawsze możliwych do pozyskania z danego obszaru powoduje niemożność<br />
generowania wartości produktu o odpowiedniej jakości przez radary<br />
z pojedynczą polaryzacja fali elektromagnetycznej.<br />
Ograniczony zasięg produktu do ok. 125 km od radaru związany między<br />
innymi z tzw. horyzontem radarowym (wraz ze wzrostem odległości od radaru<br />
wartość prawdopodobieństwa wykrywania gradu maleje). W dalszych<br />
odległościach skutkuje większymi wartościami FAR.<br />
Przydatność<br />
produktu<br />
Użyteczny produkt dla automatycznego wykrywania gradu, szczególnie jeżeli<br />
dostępne są jedynie dane o odbiciowości z pojedynczej polaryzacji. Najlepszy,<br />
gdy znana jest aktualna wysokość warstwy zamarzania.<br />
62
6. Produkty prognostyczne<br />
6.1 Śledzenie pól opadowych - RTR (Rain Tracking)<br />
Produkt RTR automatycznie śledzi i prognozuje pola opadu zlokalizowanych struktur<br />
meteorologicznych. Jest wynikiem analiz danych o odbiciowości i prędkości radialnej. Celem<br />
tych analiz jest identyfikacja pól opadowych i określenie kierunku oraz prędkości<br />
przemieszczania się pól opadowych.<br />
Przykładowy produkt typu<br />
RTR. Prognoza opadu z<br />
wyprzedzeniem na 30 min.<br />
(opcja Data = +30:00<br />
min).<br />
Alg Type - wybór algorytmu śledzenia wg. algorytmu używanego dla CTR lub VVP<br />
Data - czas wyprzedzenia prognozy podawany na każdej opcjonalnie wyświetlanej mapie<br />
produktu<br />
Produkt RTR składa się z kilku różnych warstw danych - map wyjściowych. Na produkcie<br />
głównym mogą być one prezentowane (opcjonalnie) jako mapy obrazujące:<br />
- bieżące natężenie opadu - aktualny opad wyliczany na podstawie aktualnych danych typu<br />
CAPPI lub SRI,<br />
- średnie natężenie wcześniejszego opadu,<br />
- krótkoterminowa prognoza akumulacji opadu,<br />
- krótkoterminowa prognoza natężenia opadu.<br />
63
Data: + 30:00 min. (current + time step /<br />
bieżące + krok czasowy)<br />
Spodziewana pozycja i natężenie opadu w<br />
bliskiej przyszłości, jako krótkoterminowa<br />
prognoza natężenia opadu z wyprzedzeniem<br />
czasowym na:<br />
▪ 30 minut<br />
Data: Pac + 30:00 min. (PAC+forecast time/<br />
suma opadu + czas wyprzedzenia)<br />
Oczekiwana suma opadu dla określonego<br />
czasu wyprzedzenia, jako krótkoterminowa<br />
prognoza sumy opadu z wyprzedzeniem na:<br />
▪ 10 minut,<br />
▪ 20 minut,<br />
▪ 30 minut (na przykładzie).<br />
Podsumowanie<br />
Istnieje możliwość wykorzystania produktu RTR jako danych wejściowych<br />
dla produktu PAC. Takie rozwiązanie gwarantuje uzyskanie lepszych<br />
wyników akumulacji opadu ponieważ komórka opadowa jest interpolowana<br />
pomiędzy kolejnymi cyklami obserwacji.<br />
Rezultaty prognozy RTR zależą od jakości wyników analizy wektora<br />
przesunięcia.<br />
Przydatność<br />
produktu<br />
Użyteczny produkt do zautomatyzowanego śledzenia i prognozowania pól<br />
opadu.<br />
64
6.2 Śledzenie komórek burzowych - CTR (Cell Centroid Tracking)<br />
Produkt CTR identyfikuje w zlokalizowanej strukturze meteorologicznej komórki<br />
burzowe, śledzi je i prognozuje kierunek ich przemieszczania się. Procesy te zobrazowane są<br />
na produkcie w postaci naniesionych konturów komórek: starych, aktualnych,<br />
prognozowanych.<br />
Przykładowy produkt typu<br />
CTR z wyświetlonym<br />
oknem opisującym<br />
prognozowaną komórkę<br />
burzową. Prognoza z<br />
wyprzedzeniem na 20<br />
minut.<br />
Podsumowanie<br />
Automatycznie wykrywa i śledzi kontury i środki ciężkości komórek<br />
burzowych.<br />
Może być wyświetlany jako nakładka na inne produkty.<br />
Wyniki prognozy, podawanej przez produkt CTR, są najlepsze dla<br />
odizolowanych komórek konwekcyjnych. Parametry wykorzystywanego<br />
algorytmu do wyliczenia wartości produktu powinny być dostosowane do<br />
lokalnych warunków klimatycznych.<br />
Prognozowanie komórek jest zdeterminowane także wewnętrznymi procesami<br />
zachodzącymi w komórce: zmianami kształtu oraz ograniczonym czasem<br />
życia komórki burzowej, a to nie jest uwzględniane przez algorytm.<br />
Przydatność<br />
produktu<br />
Produkt szczególnie użyteczny dla ośrodków zajmujących się<br />
prognozowaniem ultra-krótkoterminowym.<br />
Przydatny do zautomatyzowanego wykrywania i śledzenia komórek<br />
konwekcyjnych o charakterze burzowym.<br />
65
7. Rozpoznawanie zjawisk w oparciu o produkty radarowe<br />
7.1 Wykorzystanie pomiarów z podwójną polaryzacją<br />
Pozyskane dane radarowe z pomiaru falami o podwójnej polaryzacji (patrz rozdział<br />
1.3) są wykorzystane do analizy mikrofizycznych właściwości hydrometeorów tworzących<br />
zlokalizowaną strukturę meteorologiczną. Jednym z efektów tych analiz jest możliwość<br />
rozpoznania typu opadu z nich powstającego oraz wydzielenie struktur nie<br />
meteorologicznych.<br />
W celu pozyskania danych o jak najlepszej jakości zebrane dane poddawane są<br />
procesom wstępnego przetwarzania (preprocessing), który zmierza do korekcji niedostatków<br />
pomiaru (niektórych błędów pomiaru) wynikających z możliwych (ograniczonych) na dzień<br />
dzisiejszy technik pomiaru radarowego. Jest to ważny krok w procesie eksploracji danych.<br />
Ma szczególne zastosowanie do analizy zgromadzonych danych surowych pochodzących z<br />
bezpośredniego pomiaru (danych we współrzędnych sferycznych). Przetwarzanie nie<br />
skorygowanych danych, może dawać często błędne wyniki, które są rezultatem różnego<br />
rodzaju czynników wpływających na pomiar radarowy. Aby otrzymać dane dobrej jakości,<br />
przed uruchomieniem procesu analizy danych zmierzającej do pozyskania określonych<br />
produktów, dane poddawane są procesom preprocessing’u. Wstępne przetwarzanie obejmuje<br />
„czyszczenie” danych, ich ekstrakcję i selekcję pomiarów odstających od wartości progowych<br />
według ściśle zadanych algorytmów. Bazując na danych surowych: poszukiwane jest jasne<br />
pasmo i wykonywana korekta na pionowy profil odbiciowości, prowadzona jest korekta na<br />
tłumienie przez komórki o silnej odbiciowości na trasie sygnału, usuwanie wartości odstające<br />
od zadanych itp.<br />
Porównując produkty radarowe wygenerowane w oparciu o korektę danych typu 3D<br />
z wykorzystaniem procesów preprocessing’u i bez korekty, możemy potwierdzić zalety tych<br />
procesów.<br />
66
Mapa produktu MAX policzona z pomiarów bez<br />
korekcji sygnału radarowego.<br />
Mapa produktu MAX z tego samego<br />
terminu, ale policzona z wykorzystaniem<br />
procesów preprocessing’u.<br />
Na produkcie wykonanym w oparciu o dane „nie skorygowane”, na azymucie ok. 280 -290<br />
stopni występuje tłumienie (pochłanianie) sygnału wywołane przesłanianiem wiązki<br />
radarowej przez komórki o silnej odbiciowości.<br />
Zobrazowanie centralnej części produktu MAX(dBZ) przed i po procesach preprocessing’u.<br />
Wygenerowany produkt w oparciu o dane poddane procesom preprocesing’u potwierdza<br />
korektę tłumienia (między innymi) sygnału radarowego wzdłuż trasy propagacji.<br />
7.2 Rozpoznanie typu opadu – CAPPI(ET) i RHI(ET) (Echo and Hydrometeor<br />
Classification)<br />
W momencie zainstalowania w Polsce pierwszego radaru o podwójnej polaryzacji<br />
(Pastewnik) uruchomiono liczenie produktów w oparciu o dane pozyskane z tego radaru.<br />
Jednym z celów pomiaru polarymetrycznego jest szukanie możliwości i dążenie do<br />
rozpoznania typu opadu (rodzaju hydrometeorów w zlokalizowanej strukturze). Produktem<br />
takim w przypadku sieci radarów meteorologicznych POLRAD jest produkt CAPPI(ET),<br />
który może być wykonywany z zastosowaniem algorytmów sieci amerykańskiej NEXRAD, z<br />
wykorzystaniem metodyk rozpoznania typu hydrometeorów opartych na algorytmach dla<br />
67
obszaru australijskiego albo algorytmu opracowanego przez badaczy francuskich. Algorytmy<br />
rozpoznania typu opadów różnią się zasadniczo. W przypadku algorytmu australijskiego daje<br />
on możliwość rozpoznania dużej gamy hydrometeorów, którą dla polskiej strefy klimatycznej<br />
należałoby zmodyfikować. W oparciu o dane z pomiaru falą spolaryzowaną poziomo i<br />
pionowo algorytm „australijski”, dla prezentowanego niżej terminu obserwacji, rozpoznaje<br />
następujące zjawiska:<br />
- mieszanina gradu z deszczem<br />
- duży mokry grad<br />
- mały mokry grad<br />
- mokra krupa śnieżna<br />
- sucha krupa śnieżna<br />
- mokry śnieg<br />
- silny suchy śnieg<br />
- słaby suchy śnieg<br />
- deszcz<br />
- mżawka<br />
- obiekty biologiczne<br />
- zakłócenia<br />
- echa niesklasyfikowane<br />
- brak danych<br />
Natomiast zjawiska rozpoznawane przez algorytm „francuski” to:<br />
- mieszanina gradu z deszczem<br />
- duże krople opadowe<br />
- silny opad<br />
- średni opad<br />
- słaby opad<br />
- obiekty biologiczne,<br />
- zakłócenia<br />
- echa nie sklasyfikowane<br />
- brak danych<br />
68
.<br />
a) Produkt MAX<br />
b) Produkt RHI po odbiciowości<br />
Przykładowy produkt pomiaru maksymalnych odbiciowości i produkt przekroju<br />
pionowego typu RHI wykonanego wzdłuż azymutu 270 stopni z uwzględnieniem procesów<br />
preprocessing’u.<br />
Zobrazowanie przekroju pionowego typu<br />
hydrometeorów RHI(ET) rozpoznanych w<br />
zlokalizowanej strukturze meteorologicznej<br />
wzdłuż wykonanego przekroju.<br />
7.3 Możliwości rozpoznawania zjawisk w oparciu o produkty z radarów pracujących<br />
z pojedynczą polaryzacją fali<br />
Opad<br />
Produkty radarowe umożliwiają bardzo dokładny monitoring opadu. Pozwalają na<br />
określanie momentu jego pojawienia się oraz tendencji zmian, w czasie rzeczywistym,<br />
zarówno w przypadku gwałtownych, nawalnych jak i słabych opadów towarzyszących każdej<br />
porze roku. Mamy tu do dyspozycji wiele, zróżnicowanych produktów.<br />
Podstawowym produktem obrazującym natężenie opadu jest SRI. W oparciu o ten produkt<br />
tworzona jest mapa PAC godzinnych sum opadów, kilkugodzinnych sum opadu i dobowej<br />
sumy opadów. Natomiast ilość opadu potencjalnego, który w danym momencie jest<br />
oszacowany w zlokalizowanym obiekcie można ocenić w oparciu o produkt VIL. Cały obszar<br />
69
zagrożony wystąpieniem opadu pokażą nam na produkty typu MAX(dBZ) lub CMAX(dBZ).<br />
Analizę rozkładu intensywności zjawiska opadowego, dla określonej wysokości n.p.m. nad<br />
danym obszarem można prowadzić na produkcie CAPPI(dBZ). W zależności od potrzeb<br />
informacja o obszarach intensywnych opadaów może być wydzielona na produkcie typu SWI<br />
wraz z wyróżnieniem rdzenia opadowego (obszaru o najintensywniejszym opadzie - core) w<br />
przypadku szczególnie niebezpiecznych struktur opadowych. Dodatkową zaletą tego produktu<br />
jest możliwość przybliżonego oszacowania objętości struktury opadowej w [m 3 ] i ilości wody<br />
w niej zawartej wyrażonej w [tonach], czyli tego co w najbliższym czasie może, pod postacią<br />
opadu, pojawić się na powierzchni ziemi.<br />
W przypadku meteorologii radarowej gorzej przedstawia się sytuacja jeśli chodzi<br />
o monitorowanie słabych opadów typu mżawka. Pomimo, ze najmniejsza średnica kropel<br />
jakie mogą być wykrywane przez radar, to około 0,1 mm, to ich identyfikacja przez radar<br />
ograniczona jest przez tzw. horyzont radarowy i efekt zwany „przestrzeleniem chmury”<br />
(Moszkowicz, Tuszyńska, 2003). Jeszcze gorzej przedstawia się sytuacja z rozpoznawaniem<br />
mgły (składają się na nią zbyt małe drobiny hydrometeorów zawieszone tuż nad powierzchnią<br />
ziemi) czy gołoledzi, która niekiedy może być wynikiem zamarzania przechłodzonej mżawki<br />
na powierzchni ziemi.<br />
Burze<br />
Radar nie jest w stanie wykryć wyładowania atmosferycznego, które jest<br />
podstawowym zjawiskiem towarzyszącym burzy. Jednakże istnieje możliwość pośredniego (z<br />
pewnym prawdopodobieństwem) rozpoznawania burzy.<br />
W przypadku burz istnieje możliwość monitorowania ich rozwoju w całym obszarze pomiaru<br />
radarowego (skanu klasycznego – opadowego) na produktach hydrologicznych takich jak:<br />
SRI, PAC, VIL informujących o intensywności opadu i jego rozwoju nad danym obszarem.<br />
Podstawową informację o potencjalnych możliwościach pojawienia się i cyklu życia struktur<br />
burzowych można pozyskać już z produktów, generowanych z odbiciowości radarowej, typu:<br />
CAPPI(dBZ), MAX(dBZ) i EHT. Wydzielenie rdzeni struktur burzowych, w momencie<br />
najintensywniejszego ich rozwoju (Evans, J. E., Weber, M. E., 2000), dla sytuacji rozwiniętej<br />
konwekcji, można wykonać wykorzystując do tego celu produkt SWI, a do wyznaczenia<br />
obszarów narażonych na wystąpienie gradu produkt ZHAIL, który określa<br />
prawdopodobieństwo wystąpienia gradu w obszarze zlokalizowanej struktury.<br />
Analizę szczególnie niebezpiecznych struktur burzowych (opadowych) można wzbogacić<br />
poprzez obserwacje ich rozwoju na produktach dodatkowych takich jak:<br />
- przekroje pionowe: VCUT i MLVCUT dla obszarów szczególnie zagrożonych,<br />
70
- obrazy odbiciowości radarowej CAPPI(dBZ) na dodatkowych (różnych) wysokościach,<br />
- CMAX(dBZ), przy pomocy którego możemy określić obszar o najintensywniejszej skali<br />
zjawiska,<br />
- MAX(dBZ), który zobrazuje nam dodatkowo intensywność zjawiska w całej rozciągłości<br />
pionowej lub dla zadanego przedziału wysokości.<br />
W przypadku burz szczególnie przydatna może okazać się informacja jaką daje produkt<br />
EHT Echo Top w połączeniu z analizą odbiciowości na wysokości w zakresie od 5 do 8 km.<br />
Do śledzenia i określania kierunku przemieszczania się struktur burzowych i ich prędkości<br />
służyć może produkt CTR lub MAX(V), CMAX(V) - dla warstwy najniższej możliwej do<br />
wyznaczenia z danych pochodzących ze skanowania dopplerowskiego.<br />
Grad<br />
Zjawisko to najczęściej towarzyszy obszarom występowania burz. Ze względu na jego<br />
dużą szkodliwość dla różnych gałęzi gospodarki i społeczeństwa bardzo cenna jest<br />
jakakolwiek możliwość monitorowania obszarów narażonych na jego wystąpienie. Do jego<br />
identyfikacji (z pewnymi ograniczeniami) można wykorzystać, wspomniany wyżej, produkt<br />
ZHAIL (prawdopodobieństwo wystąpienia gradu). Informację w nim zawartą można<br />
wspomóc łączoną analizą produktów typu: MAX(dBZ) i EHT.<br />
Wiatr<br />
Pomimo poważnego ograniczenia jakim jest pomiar wiatru tylko w obszarach stref<br />
opadowych w przypadku radarów dopplerowskich, które pracują w sieci POLRAD tkwią<br />
duże możliwości, co do określania obszarów silnych wiatrów (przekraczających 17 m/s) oraz<br />
porywów wiatru (powyżej 5 m/s w ciągu 10 minut). Poza określaniem siły wiatru istnieje<br />
możliwość określania jego kierunku w formie tradycyjnych meteorologicznych strzałek<br />
wiatru. Podstawowym produktem obrazującym rozkład poziomy wiatru na określonej<br />
wysokości n.p.m. jest produkt HWIND. Praktycznie powinien być wykorzystywany do<br />
monitorowania wszystkich zjawisk typu wiatrowego.<br />
Dywergencja, konwergencja<br />
Z rozwojem struktur konwekcyjnych związane są najczęściej intensywne procesy<br />
przemieszczania się hydrometeorów w ich obrębie. Do oceny intensywności tych procesów<br />
wykorzystać można (poza innymi produktami) przede wszystkim produkt SWI. Umożliwia<br />
on pozyskanie informacji o aktywności procesów atmosferycznych, które mogą towarzyszyć<br />
silnej konwekcji i strukturom burzowym. Pozwala na identyfikację obszarów: lokalnej<br />
dywergencji i konwergencji, wystąpienia zjawisk typu mikroburst oraz zwiastunów<br />
mikroburst, struktur anty- i cyklonalnych. Produkt ten, poza możliwością identyfikacji<br />
71
wymienionych zjawisk, generuje ich charakterystyki z podaniem wartości wielu parametrów<br />
meteorologicznych, określa lokalizację, oszacowuje wysokości, na których mogły wystąpić.<br />
Należy jednak zaznaczyć, że analizy prowadzone są tylko w obszarze struktury opadowej,<br />
jaka możliwa jest do zidentyfikowania przez pomiar dopplerowski.<br />
Wiry i trąby powietrzne<br />
Poza produktem SWI określającym obszary potencjalnie zagrożone wystąpieniem<br />
cyrkulacji cyklonalnych i antycyklonalnych, dla szczególnie ekstremalnych sytuacji, kiedy<br />
możliwe jest pojawienie się trąb powietrznych, analizę zjawisk/procesów można uzupełnić<br />
wykorzystując dodatkowo produkty dopplerowskie typu PCAPPI(V) i CMAX(V) - do pewnej<br />
wysokości (maksymalnie 5-7 km) oraz PCAPPI(W) lub MAX(W) (Kożuchowski, 1998).<br />
Ponadto w przypadku wygenerowania przez produkt SWI symboli struktur cyklonalnych<br />
połączenie tej informacji z wymienionymi produktami może ustrzec przed ewentualni<br />
błędami, wynikającymi z algorytmu rozpoznania struktur cyklonalnych lub błędnego pomiaru<br />
radarowego (J.Szturc, 2009).<br />
Zaadoptowany w sieci POLRAD algorytm detekcji mezocyklonów opracowany został<br />
w 1985 roku dla sieci radarowej NEXRAD (Instruction Manual –Rainbow ® , 2007)<br />
Należy tutaj zaznaczyć, że produkt SWI jest tzw. produktem nakładkowym co<br />
oznacza, że produkt SWI może być wyświetlany (narysowany) na danych wygenerowanych<br />
przez inne produkty. Na przykład nałożenie na produkt SRI lub VIL (w zasięgu do 125 km)<br />
produktu SWI, pozwoli na jednoczesną analizę intensywności zjawisk wiatrowych w<br />
powiązaniu z intensywnością opadu na danym obszarze.<br />
Jako, że trąby powietrzne są zjawiskami bardzo trudnymi do zidentyfikowania, z powodu ich<br />
krótkiego cyklu życia, konieczne jest jednoczesne prowadzenie intensywnego monitoringu<br />
całokształtu sytuacji meteorologicznej na tradycyjnych produktach <strong>radarowych</strong>: CAPPI(dBZ),<br />
SRI, MAX i specjalistycznych typu: CAPPI(V), MAX(W).<br />
Uskoki wiatru<br />
Do analizy obszarów wystąpienia uskoków wiatru wykorzystywane mogą być<br />
produkty typu SHEAR. Do określenia obszarów występowania poziomego uskoku wiatru<br />
wykorzystywany jest produkt HSHEAR, a pionowego VSHEAR. W przypadku oceny<br />
uskoków wiatru w warstwie przyziemnej problemem jest mały zasięg pomiarowy wynikający<br />
nie tylko z jednoznacznie określonego zasięgu dla skanu dopplerowskiego (125 km od<br />
radaru). Dodatkowymi ograniczeniami są tutaj:<br />
- ograniczenie pomiaru radarowego wynikające z techniki pomiaru i tzw. horyzontu<br />
radarowego,<br />
72
- możliwości prowadzenia analiz tylko w obszarze pomiaru wg. algorytmu typu CAPPI<br />
(Moszkowicz, Tuszyńska, 2003),<br />
- pomiar jedynie w obrębie stref opadowych (z pewnymi wyjątkami),<br />
- lokalizacja stacji radarowej (jej wysokość n.p.m.).<br />
Jedynym rozwiązaniem tych problemów jest gęstsza lokalizacja stacji <strong>radarowych</strong> w sieci<br />
POLRAD.<br />
Turbulencje<br />
Dane o turbulencji pozyskiwane są z radarowego pomiaru rozkładu szerokości widma<br />
prędkości radialnych. Generowany w oparciu o ten pomiar produkt LTB identyfikuje obszary<br />
turbulencji w określonej warstwie atmosfery (struktury opadowej). Na dzień dzisiejszy innej<br />
możliwości analizy wymienionego zjawiska w obszarze eksploatowanych systemów<br />
teledetekcji nie ma.<br />
Marznący opad, przymrozki<br />
Możliwości monitorowania przymrozków w oparciu o dane radarowe praktycznie nie<br />
ma. Natomiast w przypadku występowania opadów marznących, szczególnie w<br />
przejściowych porach roku: zima-wiosna, jesień-zima można do ich identyfikacji posłużyć się<br />
produktami radarowymi. Wymaga to jednak dużego doświadczenia i wiedzy z dziedziny<br />
meteorologii radarowej. Jest to zagadnienie wybiegające poza zakres niniejszego<br />
opracowania.<br />
Najprostszą metodą, w tym wypadku, jest analiza produktu SRI i temperatury przy gruncie<br />
(około 0 o C) w obszarze występowania opadów.<br />
Upały – susze<br />
Wprawdzie pomiar temperatury nie leży w zakresie możliwości radaru jednak brak<br />
ech <strong>radarowych</strong> (opadów) w tym wypadku jest bardzo istotną informacją.<br />
Odnośnie takich zjawisk jak np. susza długofalowa i towarzyszące jej upały, trudno jest<br />
określić moment, kiedy się skończą. W przypadku danych <strong>radarowych</strong> pojawienie się<br />
jakichkolwiek struktur opadowych świadczyć może o ich zakończeniu lub złagodzeniu<br />
skutków. Radar w tym wypadku jest narzędziem, na którym najszybciej można to zauważyć<br />
ponieważ produkty typu MAX(dBZ) i CMAX(dBZ) są analizą atmosfery generowaną co 10<br />
minut. W przypadku mapy zbiorczej typu COMPO_MAX (patrz rozdział następny) mamy<br />
możliwość ich monitorowania na obszarze całej Polski.<br />
73
8. Zbiorcza mapa radarowa<br />
Typy mapy zbiorczej - COMP_CAPPI (Composite Constant Altitude Plan Position<br />
Indicator)<br />
Radarowe mapy zbiorcze generowane są, przez różne kraje, w oparciu o różne typy<br />
produktów, tworzonych przez sieci radarowe tych krajów.<br />
Sieć POLRAD umożliwia składanie map zbiorczych z wielu produktów. Najczęściej<br />
wykorzystywane są produkty pochodzące ze skanowania klasycznego, czyli tworzone w<br />
zasięgu 200 – 250 km od radaru.<br />
Radar Świdwin Radar Gdańsk<br />
Radar Legionowo<br />
Radar Poznań<br />
Metodyka składania radarowej<br />
mapy zbiorczej.<br />
Radar Rzeszów<br />
Radar Pastewnik<br />
Radar Ramża<br />
Radar Brzuchania<br />
Na dzień dzisiejszy są to następujące produkty radarowe:<br />
- CAPPI (dBZ) - daje mapę zbiorczą COMPO_CAPPI; [dBZ],<br />
- MAX (dBZ) - daje mapę zbiorczą COMPO_MAX; [dBZ],<br />
- SRI (dBR) - daje mapę zbiorczą COMPO_SRI; [mm/h],<br />
- PAC (dBA) - daje mapę zbiorczą COMPO_PAC; [mm],<br />
- EHT - daje mapę zbiorczą COMPO_EHT; [km].<br />
Zbiorcza mapa radarowa może być także tworzona w oparciu o takie produkty jak np.:<br />
- VIL (dBA) w przypadku gdy produkt wejściowy generowany jest do zasięgu 200 km od<br />
radaru – COMPO_VIL; [mm],<br />
- PPI (dBZ) – COMPO_PPI; [dBZ]<br />
74
lub w oparciu o produkty wiatrowe typu:<br />
- HWIND, UWT – COMPO_HWIND, COMPO_UWT; [m/s],<br />
- CAPPI(V) – COMPO_CAPPI(V); [m/s]<br />
z obszaru pomiaru dopplerowskiego (zasięg produktów wiatrowych to tylko 125 km od<br />
radaru).<br />
W przypadku map generowanych z produktów wiatrowych jest to mapa z lukami informacji,<br />
wynikającymi ze słabego zagęszczenia sieci radarowej. Warunkiem koniecznym do tworzenia<br />
kompletnych map z produktów wiatrowych jest gęstsza sieć radarów, zapewniająca pokrycie<br />
obszaru Polski przez zasięg pomiaru skanowania dopplerowskiego.<br />
Typ produktu, z którego składana jest mapa zbiorcza decyduje o tym co ona<br />
prezentuje. Poniżej przykłady podstawowych map zbiorczych generowanych aktualnie w<br />
sieci POLRAD (stan na początek roku 2011).<br />
Radarowa mapa zbiorcza<br />
typu PCAPPI<br />
Zbiorcza mapa radarowa<br />
odbiciowości radarowej na<br />
określonej wysokości n.p.m. Mapa<br />
generowana w oparciu o produkty<br />
PCAPPI z wysokości 0,7 km n.p.m.<br />
Radarowa mapa zbiorcza<br />
typu CMAX<br />
Zbiorcza mapa radarowa<br />
maksymalnych<br />
wartości<br />
odbiciowości.<br />
Mapa generowana w oparciu<br />
o produkty CMAX maksymalnych<br />
wartości odbiciowości radarowej w<br />
rzucie na jedną płaszczyznę (w<br />
układzie x-y)<br />
75
Radarowa mapa zbiorcza<br />
typu SRI<br />
Zbiorcza mapa radarowa natężeń<br />
opadów na określonej wysokości<br />
nad poziomem terenu.<br />
Mapa tworzona w oparciu<br />
o produkty SRI wygenerowane dla<br />
wysokości 1 km nad powierzchnią<br />
gruntu.<br />
Radarowa mapa zbiorcza<br />
typu PAC<br />
Zbiorcza mapa radarowa sum<br />
opadów z określonej wysokości nad<br />
poziomem ziemi.<br />
Operacyjnie generowana jest mapa<br />
w oparciu o produkty godzinnych<br />
sum opadu PAC.<br />
Radarowa mapa zbiorcza<br />
typu EHT<br />
Zbiorcza mapa radarowa wysokości<br />
wierzchołków ech <strong>radarowych</strong>.<br />
Mapa generowana w oparciu<br />
o produkty EHT - wartości<br />
wysokości wierzchołków ech<br />
<strong>radarowych</strong>.<br />
76
W ramach systemu radarowego POLRAD w trybie pracy operacyjnej oprócz<br />
godzinnej sumy opadów dla każdego terminu obserwacji (co 10 min = 144 godzinne sumy<br />
opadów w ciągu doby), liczone są:<br />
- 6-godzinne sumy opadów, na które składają się informacje o opadzie z 36 obserwacji<br />
poprzedzających wyliczanie sumy dla następujących terminów: 06:00 UTC, 12:00 UTC,<br />
18:00 UTC i 00:00 UTC,<br />
- dobowa suma opadów dla terminu obejmującego dobę hydrologiczną (doba hydrologiczna w<br />
IMGW liczona jest od godziny 06:10 danego dnia do godziny 06:00 dnia następnego);<br />
tworzona raz na dobę o godzinie 06:00 UTC.<br />
Z niektórych produktów (np. SRI, PCAPPI) mogą być także policzone mapy dla różnych<br />
wysokości nad powierzchnią gruntu (np. SRI) lub nad poziomem morza (np. PCAPPI(dBZ)).<br />
Mapy zbiorcze są szczególnie przydatne dla każdego użytkownika, praktycznie bez<br />
względu na branżę jaką reprezentuje, gdyż dają podgląd na aktualna sytuację meteorologiczną<br />
dla całego obszaru Polski. Analizując, np. przy pomocy animacji, kolejne mapy zbiorcze<br />
tworzone w odstępach 10-cio minutowych, możemy wyrobić sobie również pogląd o<br />
przewidywanym rozwoju zdarzeń meteorologicznych.<br />
W początkowych latach pracy z systemem, testowane były różne algorytmy jeśli<br />
chodzi o wybór wartości w tym obszarze mapy, gdzie dane pokrywają się. Możliwe jest do<br />
zastosowania w tym obszarze : uśrednianie wartości, wybór wartości maksymalnych, wybór<br />
wartości najbliższej radaru itp.<br />
Obecnie powszechnie udostępniana mapa radarowa na www.pogodynka.pl jest mapą zbiorczą<br />
składaną w oparciu o analizę danych, której efektem jest wybór maksymalnej wartości z<br />
obszarów gdzie możliwe jest pozyskiwanie danych z kilku radarów.<br />
Należy w tym momencie zaznaczyć, że proces liczenia mapy zbiorczej wg zadanych<br />
algorytmów przebiega w systemie prawidłowo. Jednakże powstałe mapy zbiorcze mogą<br />
nieraz uwypuklać szereg problemów związanych z specyfiką pomiaru radarowego, dając<br />
wrażenie błędów łączenia danych lub błędnej informacji.<br />
77
Wykorzystana dokumentacja i literatura<br />
Atlas D., The estimation of cloud parameters by radar, J. Meteorol., 11, 1954 . 309 - 317<br />
Battan L.J., Pogoda, tł.: S. Moszkowicz, PWN, Warszawa,1979<br />
Battan L.J., Radar observation of the atmosphere, The UNIVERSITY Chicago PRESS,<br />
Chicago Illinois U.S.A., 1973<br />
Bringi V. N., Chandrasekar V., Polarimetric Doppler weather radar, Cambridge University<br />
Press, 2001, s. 64-65<br />
Campbell S.. D., Use of features aloft n the TWDR microburst recognition algorithm, 24 th<br />
Conf. Radar Meteorology, Tallahaussee, USA, 1889, s.167-170<br />
Doppler Weather Radar System 500C – Documentation – File 01-07, Gematronik, Maj 2002<br />
Doppler Weather Radar System 360AC– Ext.Vendor Manuals – RCP1,RCP2,T/R1, 1996<br />
Doppler Weather Radar System 1500C – Documentation –j 2002; Gematronik.<br />
Douglas R.H., Hail size distribution, Proc. 11h Weather radar Conf., 1964, s. 146 - 149<br />
Doviak R.J., Zrnic D.S., Doppler radar and weather observations, Academic Press, San Diego<br />
California, 1993<br />
Evans, J. E., Weber, M. E., Weather Radar Development and Application Programs, Lincoln<br />
Laboratory Journal, Volume 12, Number 2, 2000, s 367-382<br />
Instruction Manual – Rainbow ® 5, Part 2, Release 5.18.0,: User Guide, Products &<br />
Algorithms, File Format SELEX Sistemi Integrati, kwiecień 2007<br />
Instrukcja skrócona Ravis 1.3, wydanie dokumentu 1.8, , Gematronik, 2001<br />
Klejnowski R., Atlas pogody, wyd. PASCAL, Bielsko Biała, 2007, s. 63-68<br />
Kożuchowski K., Atmosfera, klimat, ekoklimat, PWN, Warszawa, 1998<br />
Korczyk P., Wpływ wody chmurowej na drobnoskalową turbulencję – model laboratoryjny,<br />
Praca doktorska, <strong>Instytut</strong> Podstawowych Problemów Techniki PAN, Warszawa, 2008<br />
Kossakowska-Cezak U., Podstawy meteorologii i klimatologii, wyd. Szkoła Wyższa<br />
Przymierza Rodzin, Warszawa, 2007, s. 84<br />
Krasowski H., Zasada działania radaru meteorologicznego, poradnik Pilota B-767<br />
Marshall J.S., Palmer W., The distribution of raindrops with size, J.Meteorol, 5, 1948, s.165 -<br />
166<br />
Michelson D., Einfalt T., Holleman I., Gejsten U., Friderich K., Haase G., Lindskog M.,<br />
Jurczyk A., Weather Radar Data Quality In Europe: Quality Control and<br />
Charakterization, COST 717 Working Dokument WDD_MC_200403_1<br />
78
Moszkowicz S., Tuszyńska I, 2003. Meteorologia radarowa, <strong>Instytut</strong> <strong>Meteorologii</strong> i<br />
Gospodarki Wodnej, Warszawa, 2003<br />
Pietrak S., Jasiński J., Kroszczyński K., Badanie pola wiatru dla obszarów morskich na<br />
podstawie danych teledetekcyjnych i modeli numerycznych, Zeszyty Naukowe AM,<br />
Szczecin, 13/2008, s.59-64<br />
Rainbow. Krótki przewodnik, AMS Gematronic, Neuss-Rosellen, Germany, 2005,<br />
Rainbow. Podręcznik użytkownika, AMS Gematronic, Neuss-Rosellen, Germany, 2005,<br />
Rainbow. Produkty i algorytmy, AMS Gematronic, Neuss-Rosellen, Germany, 2005,<br />
Rainbow. Format plików, AMS Gematronic, Neuss-Rosellen, Germany, 2005,<br />
Retallack B.J., Podstawy meteorologii, skrypt WMO 1984, wyd. IMGW 1991 r., s. 31-37<br />
Ronald E. Rinehart, Radar for meteorologists, Published by R.E. Rinehart in Grand Forks,<br />
N.D, 1991<br />
Szostak Ł., Radary meteorologiczne – interpretacja wskazań, seminarium dyplomowe,<br />
Politechnika Rzeszowska, 2001/2002<br />
Tuszyńska I., Wykorzystanie technologii VSAT w sieci radarów meteorologicznych, <strong>Instytut</strong><br />
<strong>Meteorologii</strong> i Gospodarki Wodnej, Warszawa, 2005<br />
Waldteufel P., Cobrian H., On the analysis of single Doppler data, J.Appl. Meteorol. 18,<br />
1979, s. 532-542<br />
Waldvogel A., Federer B., Grimm P., Criteria fort he detection of hail cells, J.Appl.<br />
Meteorol., 18, 1979, s. 521-525<br />
http://pl.wikipedia.org/wiki/Radar<br />
www.CIMMS.ou.edu, Oklahoma, 2005<br />
http://pl.wikipedia.org/wiki/Radar<br />
79