See the description of the numerical scheme here.

TL;DR: We can set up an ODE system which describes the spread of water over the last physical layer timestep, and solve that with a single timestep of an implicit numerical method. The result can then be used to compute the amount of substance that came along with those volumes.