Growing occurrence of extreme floods in the Czech Republic has attracted attention to the security of protective earthfill embankments along the rivers. A suddenly increased amount of water on the waterside slope of the embankment may have destroying or even catastrophic consequences. Predictions of seepage patterns through the earth body are usually done considering the saturated flow beneath the free water level only, neglecting the saturated-unsaturated character of the soil water dynamics within earthfill dams. The importance of water dynamics within an earth dam is known and may be addressed using numerical simulation models. In this study the solution based on transient simulation of seepage through protection levee using saturated-unsaturated theory is presented. Simulations were carried out by a two-dimensional numerical model based on Richards’ equation for water flow in porous medium. It has been shown that proposed approach is, with certain limitations, suitable for large scale engineering applications. and Vyšší pravděpodobnost výskytu extrémních klimatických jevů obrací pozornost k ochraně před následky, které tyto jevy způsobují. Zájem se soustředí na protipovodňové ochranné zemní hráze a jejich bezpečnost při povodních. Výpočet průsaku zemními hrázemi se často omezuje pouze na tu část hráze, která byla plně nasycená vodou, to znamená na plně nasycené proudění. Tento způsob modelování průsaku je dodnes považován za standardní, přestože je velmi limitující. Bez zahrnutí nenasycené části tělesa hráze je zanedbán vliv časově i prostorově proměnlivého pole vlhkostí (např. při infiltraci vody ze srážky) na polohu hladiny. V naší studii je simulováno proudění v tělese hráze s použitím numerického modelu, který umožňuje řešit proudění vody v proměnlivě nasyceném heterogenním pórovitém prostředí, s obecnými okrajovými podmínkami. Výsledky potvrdily, že přístup, který uvažuje proudění i v nenasycené části hráze, lze k řešení průsaků zemními tělesy úspěšně využít.