A 2D model for predicting clogging of a landfill leachate collection system and subsequent leachate surface position (mounding) is described. A transient finite element fluid flow model is combined with a reactive, multiple-species finite element transport model. The transport model considers biological growth and biodegradation, precipitation, and particle attachment and detachment. It uses a geometrical relationship to establish porosity from the computed thickness of the accumulated clog matter and a relationship between the porosity and hydraulic conductivity of elements in the system. The model represents the flow path within the drainage layer in profile. An iterative method is used to solve for the new hydraulic heads, surface and internal nodal positions, and redistributed clog properties (clog quantity, porosity, hydraulic conductivity) for each element and for each time step. The porosity (and consequently hydraulic conductivity) of the media can therefore change spatially and temporally. The mesh is regenerated automatically each time step (including the addition or subtraction of nodes) taking into account allowable element aspect ratios, the interfaces between differing hydrostratigraphic layers, and static point sources and openings. An integrated alternate solution for very thin mounds is included. The application of the model is demonstrated using a hypothetical field case.