82 tysiące pomiarów Wisły nie oznacza 82 tysięcy niezależnych eksperymentów, lecz niemal nieprzerwany zapis codziennego poziomu rzeki w Warszawie obejmujący okres od 1 lutego 1799 roku do 31 października 2024 roku, informuje redakcja TopFlop. Zespół badaczy zebrał dokładnie 82 453 dobowe wartości i ujednolicił je tak, aby pomiary wykonywane w czasach zaborów można było analizować razem z wynikami współczesnych elektronicznych wodowskazów. Powstała w ten sposób 225-letnia seria jest jedną z najdłuższych tego typu baz hydrologicznych w Europie.
- 82 453 wartości nie powstały w jednym systemie pomiarowym
- Jak z tysięcy pomiarów policzyć średnią, która rzeczywiście coś mówi
- Trend nie jest linią narysowaną między pierwszym i ostatnim pomiarem
- Co w takim zbiorze jest anomalią
- Największa trudność: zanim policzy się anomalie, trzeba znaleźć anomalie pomiarowe
- 2374 dni bez oryginalnego pomiaru: co zrobiono z lukami
- Algorytm również może stworzyć pozorną anomalię
- Jak analityk przeszedłby od 82 tysięcy liczb do wniosku
- Poziom wody a przepływ: liczby wyglądają podobnie, lecz mierzą coś innego
- Co 225 lat danych rzeczywiście pozwala badać
- Czy 82 tysiące pomiarów potwierdza, że Wisła wysycha
- Pytania i odpowiedzi
Sama liczba rekordów nie daje jednak odpowiedzi na pytanie, czy Wisła rzeczywiście staje się coraz płytsza, kiedy występowały nietypowe epizody ani czy obecny poziom można uznać za wyjątkowy. Najpierw trzeba ustalić, czy wartości z różnych epok rzeczywiście oznaczają to samo. W przypadku warszawskiej Wisły zmieniały się lokalizacje wodowskazów, poziomy ich zera, stosowane jednostki, dokładność odczytów, godziny pomiaru, a nawet system kalendarzowy. Dopiero po korektach można liczyć średnie, analizować rozkład danych i szukać anomalii bez mieszania zmian rzeki ze zmianami sposobu jej mierzenia.
82 453 wartości nie powstały w jednym systemie pomiarowym
Najstarsze zachowane codzienne obserwacje poziomu Wisły w Warszawie pochodzą z 1799 roku. Pierwszy wodowskaz znajdował się przy moście Ponińskiego, w rejonie dzisiejszej ulicy Bednarskiej. W 1864 roku urządzenie przeniesiono na most Kierbedzia, w 1937 roku punkt pomiarowy przesunięto ponownie, w 1958 roku trafił on w okolice Portu Praskiego, a od początku roku hydrologicznego 2018 pomiary wróciły w pobliże historycznej lokalizacji na lewym brzegu.
Każde takie przeniesienie ma znaczenie statystyczne. Odczyt „200 cm” na jednym wodowskazie nie musi oznaczać identycznego położenia zwierciadła wody jak 200 cm zapisane sto lat wcześniej w innym miejscu i względem innego poziomu odniesienia.
Jeszcze większym problemem było tzw. zero wodowskazu. Do 1834 roku pierwszy warszawski wodowskaz nie miał stałego zera powiązanego z geodezyjnym punktem odniesienia. Później jego poziom wielokrotnie korygowano. W 1958 roku zero obniżono aż o 200 cm, aby uwzględnić obserwowane obniżanie dna warszawskiego odcinka Wisły.
Dlatego analiza szeregu bez wcześniejszej homogenizacji mogłaby stworzyć sztuczny „trend”, który w rzeczywistości wynikałby z technicznej zmiany skali.
„Daily water level observations from the Warsaw gauge date back to 1799.”
— autorzy pracy opublikowanej w „Earth System Science Data”.
Dla szerszego kontekstu geograficznego przydatne jest także wcześniejsze opracowanie o najdłuższych rzekach Polski, przebiegu Wisły i miastach położonych nad jej biegiem. Wisła ma 1047 km długości, a Warszawa znajduje się na jej środkowym odcinku, dlatego pomiar w stolicy opisuje tylko jeden punkt dużego systemu rzecznego, a nie całą rzekę.
Jak z tysięcy pomiarów policzyć średnią, która rzeczywiście coś mówi
Najprostsza średnia arytmetyczna powstaje przez dodanie wszystkich obserwacji i podzielenie wyniku przez ich liczbę. Przy 82 453 rekordach samo działanie jest banalne dla komputera. Problemem jest interpretacja.
Średnia obliczona ze wszystkich dni od 1799 do 2024 roku może być użyteczna jako punkt odniesienia, ale nie mówi, jak wyglądał typowy poziom Wisły w lipcu, zimą, podczas wezbrania czy w okresie niskich przepływów. Rzeka ma wyraźny cykl sezonowy, dlatego analityk zwykle dzieli dane na mniejsze grupy.
Można więc liczyć:
- średnią dla całego szeregu;
- średnie roczne;
- średnie miesięczne;
- średnie dla poszczególnych pór roku;
- średnią kroczącą, np. z 30 dni, roku lub dekady;
- średnią dla konkretnego okresu porównawczego, np. 1951–2000 albo 1991–2020.
Dopiero takie zestawienia pokazują, czy obserwowana różnica wynika z jednego ekstremalnego dnia, sezonu czy z dłuższej zmiany.
W hydrologii sama średnia też nie wystarcza. Trzeba zestawić ją z medianą, minimum, maksimum, percentylami oraz zmiennością. Jeśli przez większość roku poziom wynosi około 100–150 cm, lecz kilka wielkich wezbrań daje wartości wielokrotnie wyższe, średnia może zostać podniesiona przez stosunkowo niewielką liczbę ekstremów.
Mediana odpowiada wtedy na inne pytanie: jaka wartość znajduje się dokładnie w środku uporządkowanego zbioru obserwacji.
Autorzy bazy przypominają, że wcześniejsze rekonstrukcje Wisły koncentrowały się głównie na rocznych wartościach maksymalnych, średnich i minimalnych. Nowy zbiór zachowuje natomiast rozdzielczość dobową, co pozwala analizować nie tylko wartość szczytową, lecz również czas rozpoczęcia, długość i dynamikę całego epizodu hydrologicznego.
Trend nie jest linią narysowaną między pierwszym i ostatnim pomiarem
Najczęstszy błąd przy analizowaniu długich szeregów polega na porównaniu pierwszej i ostatniej wartości. Jeżeli w 1799 roku wodowskaz pokazywał więcej niż w 2024 roku, nie oznacza to automatycznie, że Wisła przez 225 lat systematycznie traciła wodę.
Trend musi uwzględniać tysiące obserwacji między tymi punktami. Można go wyznaczać regresją liniową, metodami nieparametrycznymi albo analizować oddzielnie w kolejnych okresach. Przy danych hydrologicznych bada się również sezonowość i autokorelację, ponieważ poziom dzisiejszy nie jest statystycznie niezależny od poziomu wczorajszego.
Jeszcze istotniejsze jest to, że stan wody nie jest tym samym co przepływ. Stan wody jest wysokością powierzchni względem zera konkretnego wodowskazu. Przepływ mówi, jaka objętość wody przechodzi przez przekrój rzeki w jednostce czasu.
Warszawskie koryto w ciągu dwóch stuleci było regulowane i przekształcane. Zmieniała się jego geometria, następowała erozja dna, budowano infrastrukturę i zabezpieczenia przeciwpowodziowe. Te procesy mogą powodować zmianę poziomu wody nawet wtedy, gdy przepływ nie zmienia się proporcjonalnie.
Badacze podkreślają to wprost: długookresowego trendu poziomu na wodowskazie nie należy automatycznie interpretować jako takiego samego trendu przepływu ani jako bezpośredniego dowodu wysychania klimatu.
„Long-term trends in recorded stage should therefore not be interpreted directly as equivalent trends in discharge.”
— autorzy bazy 225-letnich pomiarów Wisły.
To rozróżnienie ma praktyczne znaczenie. Informacja „poziom Wisły spadł o X cm” bez podania wodowskazu, okresu odniesienia i przepływu może prowadzić do błędnych wniosków.
O tym, dlaczego stan rzeki w Warszawie zależy również od procesów zachodzących daleko od stolicy, można przeczytać w analizie regionów i dopływów wpływających na ilość wody docierającej Wisłą do Warszawy. Opady w dorzeczu i reakcja dopływów potrafią zmieniać sytuację w stolicy z opóźnieniem.
Co w takim zbiorze jest anomalią
Anomalia nie musi oznaczać rekordu. W statystyce jest to obserwacja nietypowa względem przyjętego punktu odniesienia.
Jeżeli średni poziom dla konkretnego dnia roku wynosiłby przykładowo 150 cm, wartość 300 cm byłaby dużym dodatnim odchyleniem. Jednak ten sam poziom może nie być niezwykły w miesiącu, w którym regularnie występują wezbrania.
Dlatego anomalię można obliczać względem:
| Punkt odniesienia | Co pokazuje |
|---|---|
| średniej całego szeregu | czy dzień odbiega od ogólnego poziomu |
| średniej miesięcznej | czy wartość jest nietypowa dla danego miesiąca |
| średniej sezonowej | czy odbiega od typowego przebiegu pory roku |
| średniej wieloletniej | czy określony rok różni się od okresu bazowego |
| odchylenia standardowego | jak duże jest odchylenie w jednostkach zmienności |
| percentyli | czy obserwacja należy np. do 1% najwyższych lub najniższych wyników |
W praktyce bardzo użyteczny jest standaryzowany wynik typu z-score:
z = (wartość – średnia) / odchylenie standardowe.
Wynik bliski zeru oznacza wartość znajdującą się w pobliżu średniej. Z-score równy +2 wskazuje na poziom o dwa odchylenia standardowe wyższy od średniej, a –2 na odpowiednio niski.
Nie oznacza to jednak, że każda obserwacja przekraczająca dwa odchylenia jest błędem. Powódź może być autentyczną ekstremalną obserwacją. Błędny zapis z archiwum także może być ekstremalny. Statystyka wykrywa podejrzany punkt, ale hydrolog musi ustalić jego przyczynę.
Największa trudność: zanim policzy się anomalie, trzeba znaleźć anomalie pomiarowe
W historycznej serii warszawskiej Wisły stosowano kilka systemów jednostek: warszawskie, holenderskie, nowopolskie, angielskie i rosyjskie. Dopiero w okresie I wojny światowej pomiary konsekwentnie przeszły na system metryczny z dokładnością centymetrową.
Zmiana skali sama może wyglądać jak zmiana zachowania danych.
Zespół badający bazę zastosował między innymi dwa wskaźniki jakości: UDV, czyli liczbę unikalnych dobowych wartości w roku, oraz DWC — liczbę przypadków, gdy kolejnego dnia zapisano dokładnie tę samą wartość.
Jeżeli pomiar wykonywany jest dokładniej, oczekuje się większej liczby unikalnych wartości i mniejszej liczby dni bez zmiany.
Analiza wykazała dwa istotne punkty zmiany: około 1818 i 1915 roku. Dla UDV statystyki testu SNHT wyniosły odpowiednio 26,4 oraz 18,42, podczas gdy wskazany w badaniu próg 99-procentowego poziomu ufności dla szeregu tej długości wynosił 12,982.
Po 1915 roku średnia roczna liczba unikalnych wartości wzrosła do 162, przy odchyleniu standardowym 18. Dla okresu 1800–1915 wynosiła 103 przy odchyleniu 22.
Jednocześnie średnia liczba dni, w których wodowskaz dzień po dniu pokazywał identyczną wartość, spadła z 74 do 30. Różnice były statystycznie istotne na poziomie p<0,001.
To nie dowód na nagłą zmianę Wisły w 1915 roku. Jest to przede wszystkim sygnał zmiany jakości i rozdzielczości pomiarów.
„The second point, 1915, corresponds to the historical switch of the Warsaw gauge to the metric system.”
— wyjaśniają autorzy badania.
2374 dni bez oryginalnego pomiaru: co zrobiono z lukami
Pełny szereg liczy 82 453 dni, ale 2374 dni, czyli około 2,9% całej serii, pierwotnie nie miało dostępnej obserwacji.
Największe luki przypadły na dwa wyjątkowe okresy.
Pierwszy trwał od 1 grudnia 1815 do 13 lipca 1817 roku. Zator lodowy zniszczył wtedy most i znajdujący się przy nim wodowskaz. Druga duża luka obejmowała okres od 29 lipca 1914 do 16 października 1915 roku i była związana z I wojną światową.
Krótkie luki, do siedmiu dni, można było uzupełniać prostszymi metodami. Dłuższe rekonstruowano z wykorzystaniem obserwacji z wodowskazów położonych w Toruniu i na Cyplu Mątowskim.
Odległości są znaczne: Toruń leży około 214 km poniżej Warszawy, a Cypel Mątowski około 366 km. Pomiędzy Warszawą a tymi punktami do Wisły uchodzą kolejne rzeki, w tym Narew, dlatego relacja nie jest idealna.
Do rekonstrukcji zastosowano sieć neuronową LSTM, czyli model przygotowany do analizy sekwencji. Nie chodziło o prognozowanie przyszłego stanu Wisły. Algorytm miał odtworzyć brakujący dzień na podstawie znanych fragmentów przebiegu rzeki w innych punktach.
Model uwzględniał m.in. poziom wody na stacji referencyjnej, sezonowość, tempo zmian dzień do dnia oraz opóźnienie przejścia fali hydrologicznej między stacjami.
Do oceny jakości rekonstrukcji wykorzystano MAE, RMSE, NSE i KGE. Pierwszy wskaźnik mierzy przeciętną wielkość błędu, drugi silniej karze duże pomyłki, a NSE i KGE służą do sprawdzania, jak dobrze odtworzona została dynamika hydrologiczna.
Algorytm również może stworzyć pozorną anomalię
Rekonstrukcja nie jest równoznaczna z rzeczywistym historycznym pomiarem. Autorzy bazy wyraźnie oznaczyli okresy o większej niepewności.
Szczególnie trudne okazały się zimy 1804, 1805, 1809 i 1815 roku. Zjawiska lodowe mogły podnosić wodę w jednym miejscu, ale nie powodować identycznej reakcji kilkaset kilometrów dalej. W takiej sytuacji zależność między wodowskazami przestaje działać tak dobrze jak w normalnych warunkach.
W jednym z modeli opartych na Cyplu Mątowskim zimowy średni błąd bezwzględny MAE osiągał 43,26 cm, a RMSE 55,81 cm. Latem wartości były wyraźnie mniejsze: odpowiednio 23,50 i 32,63 cm.
Oznacza to, że anomalia widoczna w odtworzonych danych musi być sprawdzona pod kątem źródła rekonstrukcji.
Dobry zbiór danych nie ukrywa niepewności. Informuje, które wartości są bezpośrednimi obserwacjami, które zostały przeliczone, a które odtworzono.
Jak analityk przeszedłby od 82 tysięcy liczb do wniosku
Praktyczna analiza takiego zbioru mogłaby wyglądać następująco:
- sprawdzenie kompletności i typów danych;
- identyfikacja wartości rekonstruowanych;
- kontrola jednostek i punktów odniesienia;
- rozdzielenie obserwacji na lata, miesiące i sezony;
- wyliczenie średnich, median, minimów, maksimów i percentyli;
- obliczenie odchylenia standardowego oraz rozstępu międzykwartylowego;
- przygotowanie średnich kroczących;
- osobne sprawdzenie zmian technicznych wodowskazu;
- analiza trendu w wybranych stabilnych okresach;
- wykrywanie wartości nietypowych;
- porównanie anomalii z historycznymi powodziami, suszami, zjawiskami lodowymi i zmianami koryta.
Dopiero po wykonaniu tych etapów można odpowiedzialnie napisać, że pewien wzorzec jest trwały.
Samo znalezienie linii opadającej w arkuszu kalkulacyjnym nie wystarcza.
Poziom wody a przepływ: liczby wyglądają podobnie, lecz mierzą coś innego
To jeden z kluczowych punktów całej analizy. Wodowskaz podaje wysokość zwierciadła wody względem ustalonego zera. Przepływ wyrażany jest zwykle w metrach sześciennych na sekundę.
Dwie różne geometrie koryta mogą przy takim samym przepływie dawać inny poziom.
Dane o dobowych przepływach Wisły w Warszawie są cyfrowo dostępne w zasobach IMGW dla okresu od 1 listopada 1950 do 30 października 2024 roku — informowali autorzy podczas procesu recenzji pracy.
Współczesne dane hydrologiczne IMGW zawierają m.in. stan wody, czas pomiaru, przepływ, temperaturę wody oraz informacje o zjawiskach lodowych. Oficjalny interfejs danych można sprawdzić w publicznej bazie hydrologicznej IMGW.
Porównanie poziomu z przepływem jest szczególnie ważne przy analizie długoterminowej. Jeśli dno rzeki się obniżyło, ten sam przepływ może odpowiadać niższemu wskazaniu wodowskazu niż wcześniej.
Co 225 lat danych rzeczywiście pozwala badać
Największą wartością bazy nie jest jedna spektakularna liczba.
Jej znaczenie wynika z możliwości prześledzenia pełnej sekwencji zdarzeń — dzień po dniu przez ponad dwa stulecia.
Dzięki temu można badać:
- czas występowania wysokich i niskich stanów;
- sezonowość wezbrań;
- długość epizodów niskiej wody;
- szybkość narastania i opadania fali;
- zmiany częstotliwości ekstremalnych zdarzeń;
- okresy występowania zjawisk lodowych;
- relację hydrologii z pogodą, klimatem oraz zmianami użytkowania rzeki.
Seria obejmuje okres od schyłku małej epoki lodowej, przez industrializację i duże przekształcenia Warszawy, po współczesny klimat. Polska Akademia Nauk wskazuje, że jest to jedna z najdłuższych tego rodzaju serii pomiarowych w Europie.
„Powstała jedna z najdłuższych tego rodzaju serii pomiarowych w Europie.”
— Polska Akademia Nauk, 1 października 2026 roku.
Pełny artykuł naukowy wraz z metodologią opublikowano 26 sierpnia 2026 roku w „Earth System Science Data”. Zbiór jest dostępny w otwartym dostępie w repozytorium Zenodo.
W interpretacji współczesnych zmian przydatny jest również szerszy obraz systemu Wisły. Rzeka zbiera wodę z rozległej części kraju, a Warszawa reaguje zarówno na sytuację w górnej Wiśle, jak i na dopływy. Ten mechanizm szerzej pokazuje wcześniejszy materiał o tym, skąd bierze się woda w Wiśle i które regiony wpływają na jej poziom w Warszawie.
Czy 82 tysiące pomiarów potwierdza, że Wisła wysycha
Nie można wyciągnąć takiego wniosku wyłącznie z trendu wskazań warszawskiego wodowskazu.
Właśnie na tym polega najważniejsza lekcja płynąca z nowej bazy. Długą serię należy interpretować równocześnie z informacjami o zmianach koryta, położeniu wodowskazów, ich zerze, przepływach oraz metodach pomiaru.
Warszawski odcinek rzeki został silnie przekształcony. Koryto jest zawężone, a jego dno podlegało erozji. W 1958 roku zmiana poziomu zera wodowskazu o 200 cm była bezpośrednio związana z obserwowanym obniżaniem dna.
Dlatego stwierdzenie, że „poziom jest niższy niż 100 lat temu”, może być prawdziwe jako opis wskazania wodowskazu, ale nie odpowiada jeszcze na pytanie, dlaczego jest niższy.
Do odpowiedzi potrzebne są równolegle dane o przepływie, opadach, dopływach, regulacji rzeki i geometrii koryta.
Pytania i odpowiedzi
Ile dokładnie pomiarów znajduje się w nowej bazie Wisły?
Seria obejmuje 82 453 dobowe obserwacje poziomu Wisły w Warszawie od 1 lutego 1799 do 31 października 2024 roku.
Czy wszystkie 82 453 wartości są bezpośrednimi pomiarami?
Nie. W pierwotnych materiałach brakowało danych dla 2374 dni, czyli około 2,9% całego okresu. Luki uzupełniono między innymi przy użyciu modelu LSTM oraz interpolacji.
Co jest lepsze do analizy Wisły: średnia czy mediana?
Obie miary odpowiadają na inne pytania. Średnia uwzględnia wszystkie wartości, ale silnie reaguje na skrajne wezbrania. Mediana pokazuje wartość środkową i jest odporniejsza na ekstrema. W praktycznej analizie hydrologicznej najlepiej porównywać je razem z percentylami, minimum, maksimum i miarami zmienności.
Czy niski poziom na wodowskazie oznacza niski przepływ?
Nie zawsze. Stan wody zależy także od geometrii koryta i poziomu zera wodowskazu. Ten sam przepływ może w różnych warunkach korytowych odpowiadać innemu stanowi wody.
Jak najprościej znaleźć anomalię w danych?
Najpierw trzeba ustalić właściwy punkt odniesienia, np. średnią dla danego miesiąca. Następnie można policzyć różnicę względem średniej, z-score albo sprawdzić pozycję wartości w rozkładzie percentylowym. Dopiero później trzeba zweryfikować, czy nietypowy wynik jest rzeczywistym zdarzeniem hydrologicznym czy problemem pomiarowym.
Dlaczego seria rozpoczynająca się w 1799 roku jest tak cenna?
Ponieważ obejmuje ponad dwa stulecia i zachowuje rozdzielczość dobową. Pozwala badać nie tylko pojedyncze rekordy, lecz także częstotliwość, czas trwania, sezonowość i dynamikę zdarzeń hydrologicznych w różnych okresach historycznych.
82 tysiące pomiarów Wisły pokazuje przede wszystkim, że duży zbiór danych nie daje odpowiedzi sam z siebie. Liczby trzeba najpierw ujednolicić, sprawdzić ich źródło, oddzielić bezpośrednie obserwacje od rekonstrukcji i zrozumieć, co fizycznie reprezentuje każda kolumna. Dopiero później średnia zaczyna mieć sens, trend przestaje być przypadkową linią, a anomalia może zostać odróżniona od zmiany wodowskazu, błędu zapisu czy zjawiska lodowego.
Właśnie dlatego 225-letnia seria Wisły ma większą wartość niż pojedynczy rekord niskiej lub wysokiej wody: pozwala zobaczyć zdarzenie na tle dziesiątek tysięcy dni, a nie tylko porównać dwa wybrane momenty.
Warto przeczytać także nasz kolejny materiał, w którym szerzej wyjaśniamy podobny temat: CERN szuka miniaturowych czarnych dziur: jak mała może być czarna dziura i jak ją wykryć
