Mathematical models provide insight to analyse peatland behaviour. However, the omission of mechanical processes by the existing models leads to uncertainties in their outputs. We proposed a peatland growth model in 2D that incorporates mechanical, ecological, and hydrological factors, together with the effect of spatial heterogeneity on the peatland system. Our model might assist in understanding the complex interactions and the impact of climate change on the peatland carbon balance.