  We exhibit a family of graphs that develop turning singularities (i.e. their
Lipschitz seminorm blows up and they cease to be a graph, passing from the
stable to the unstable regime) for the inhomogeneous, two-phase Muskat problem
where the permeability is given by a nonnegative step function. We study the
influence of different choices of the permeability and different boundary
conditions (both at infinity and considering finite/infinite depth) in the
development or prevention of singularities for short time. In the general case
(inhomogeneous, confined) we prove a bifurcation diagram concerning the
appearance or not of singularities when the depth of the medium and the
permeabilities change. The proofs are carried out using a combination of
classical analysis techniques and computer-assisted verification.
