The formulation and existence theory is presented for a system modeling diffusion of a slightly compressible fluid through a partially saturated poroelastic medium. Nonlinear effects of density, saturation, porosity and permeability variations with pressure are included, and the seepage surface is determined by a variational inequality on the boundary.