Hysteretic upscaled constitutive relationships for vertically integrated porous media flow.