How to simulate slope stability in saturated and unsaturated regions separated by wetting lines

I have already seen how to do the strength reduction for slope stability simulation in this example(https://www.opengeosys.org/docs/benchmarks/small-deformations/mechanics-slope-stability/), but I am not sure how to consider the effect of soil saturation in the example. Many thanks.