Der Wasseranteil, der dem Grundwasser zuströmt, wird mit dem Perkolationsansatz nach Glugla (1969) zeitlich und mengenmäßig verzögert. Die Retention nimmt mit zunehmender Profiltiefe zu. Die Profile werden in k Schichten unterteilt, dabei wird die Schichtdicke so festgelegt, dass jede Schicht innerhalb eines Tages durchsickert wird, d.h. die Schichtdicke (Mx) ist abhängig vom kF-Wert. Ebenso wie bei der Berechnung des Aufteilungsoperators werden auch hier wieder Anfangs- (WSA) und Endwassergehalt (WSE) für jeden Zeitschritt benötigt, wobei der Endwassergehalt einer Schicht dem Anfangswassergehalt dieser Schicht im Folgezeitschritt entspricht. Im ersten Zeitschritt wird der Anfangswassergehalt festgesetzt, indem der für das gesamte Profil geltende Anfangswassergehalt (Sanf) auf die Schichten aufgeteilt wird. In den k Schichten spielen sich dann immer wieder (je nach Zusickerungsmenge aus der darüber liegenden Schicht (RG2pot_(k-1)) die gleichen Verzögerungsvorgänge ab. Der Endwassergehalt der jeweiligen Schicht k wird bei einer Zusickerung (RG2pot_(k-1) > 0) wie folgt berechnet:

Gl. 13

mit

Gibt es keine Zusickerung aus der darüber liegenden Schicht (RG2pot_(k-1) = 0), berechnet sich der Endwassergehalt nur in Abhängigkeit vom Anfangswassergehalt und dem bodenspezifischen Parameter (C) nach GLUGLA (1969):

Gl. 14

Der Anteil des Versickerungswassers, der eine Schicht vertikal durchsickert (RG2) berechnet sich aus der Differenz zwischen Anfangs- und Endwasserhaushalt zuzüglich der aus der darüber liegenden Schicht zusickernden Wassermenge:

Gl. 15

Für die folgende Schicht k wird der so berechnete Wert für RG2 wieder als Zusickerung (RG2pot_(k-1)) verwendet. Das Sickerwasser, das die letzte Schicht durchsickert, ist dann letztendlich das Wasser, das die Grundwasseroberfläche erreicht (RG2=GWN).