Abstract. Higher-dimensional models of peatland development are required to analyse the influence of spatial heterogeneity and complex feedback mechanisms on peatland behaviour. However, the current models exclude the mechanical process that leads to uncertainties in simulating the spatial variability in the water table position, vegetation composition, and peat physical properties. Here, we propose MPeat2D, a peatland development model in two dimensions, which considers mechanical, ecological, and hydrological processes together with the essential feedback from spatial interactions. MPeat2D employs poroelasticity theory that couples fluid flow and solid deformation to model the influence of peat volume changes on peatland ecology and hydrology. To validate the poroelasticity formulation, the comparisons between numerical and analytical solutions of Mandel's problems for two-dimensional test cases are conducted. The application of MPeat2D is illustrated by simulating peatland growth over 5000 years above a flat and impermeable substrate with free-draining boundaries at the edges, using constant and variable climate. In both climatic scenarios, MPeat2D produces lateral variability in the water table depth, which results in the variation in the vegetation composition. Furthermore, the drop in the water table at the margin increases the compaction effect, leading to a higher value of bulk density and a lower value of active porosity and hydraulic conductivity. These spatial variations obtained from MPeat2D are consistent with the field observations, suggesting plausible outputs from the proposed model. By comparing the results of MPeat2D to a one-dimensional model and a two-dimensional model without the mechanical process, we argue that mechanical–ecohydrological feedbacks are important for analysing spatial heterogeneity, shape, carbon accumulation, and resilience of peatlands.