Unstable solutions when solving Richards' equation. Need help!

refine mesh and use Picard iteration