W pracy dokonano weryfikacji wybranych wzorów empirycznych do szacowania kwantyli przepływów maksymalnych rocznych o określonym prawdopodobieństwie przewyższenia, w zlewniach południowej Polski. Dane do analizy, obejmujące wielolecie 1986-2015, pozyskano z Instytutu Meteorologii i Gospodarki Wodnej w Warszawie. Dane poddano statystycznej weryfikacji na jednorodność i niezależność, istotność trendu, występowanie elementów odstających oraz równości wariancji. Kwantyle przepływów maksymalnych rocznych o określonym prawdopodobieństwie przewyższenia oszacowano za pomocą rozkładu Pearsona typu III oraz wzorów empirycznych: Punzeta i obszarowego równania regresji. Weryfikacja formuł empirycznych wykazała znaczne różnice pomiędzy rezultatami uzyskanymi za pomocą analizowanych wzorów oraz metody statystycznej. Świadczą o tym średnie wartości błędów względnych oszacowania kwantyli przepływów maksymalnych rocznych, które dla wzoru Punzeta wynosiły 64% natomiast dla obszarowego równania regresji 62%.