A finite element collocation method for variably saturated flows in porous media