Hi,
When using an unsaturated (Richards' equation) model, FEFLOW handles both the saturated and unsaturated zone in the same model. You can adjust the temporal parameters of the model in "Temporal&Control Data" Menu. You can also define time-vaying (power-) functions and use them for transient modelling of BC, material data, etc.