[1] A model is developed which couples fully saturated porous compaction to the viscous-plastic deformation of the skeleton matrix. The Darcy fluid flow during compaction is described by an advection-diffusion equation for the excess pressure with two source/sink terms that depend on the mechanical compressibility and viscous compaction of the pore space, the latter representing the effect of pressure solution. The incompressible deformation of the composite medium is described by a force balance equation and its rheology can be viscous, plastic, or viscoplastic (Bingham material). For the plastic and viscoplastic cases, the coupling between the compacting and plastically deforming parts of the system is through the Drucker-Prager frictional-plastic yield criterion modified by Terzaghi's principle, so that the yield strength depends on the effective dynamical pressure. The coupled system is solved using a two-dimensional (2-D) finite element method. Two problems are solved to demonstrate the behavior of our theory. The first considers compaction of a uniform sediment layer. The numerical results agree with the predictions of the nondimensional control parameters and previously published results. The second problem concerns 2-D kinematic progradation of deltaic sediments. Substratum and delta sediments have the same compaction properties and a Bingham rheology during deviatoric deformation, such that the delta undergoes linear postyield viscous flow. For certain depositional regimes, overpressure is generated. When pore pressures approach critical values, yielding occurs and the delta front fails and becomes unstable, spreading gravitationally under its own weight. The flow velocity is limited to geological rates by the Bingham viscosity. For the range of parameter values considered, pressure solution is the most effective mechanism for generating nearlithostatic fluid pressures that lead to initial failure, and it appears that mechanical compaction hardly contributes to the fluid overpressure at this stage.Citation: Morency, C., R. S. Huismans, C. Beaumont, and P. Fullsack (2007), A numerical model for coupled fluid flow and matrix deformation with applications to disequilibrium compaction and delta stability,