Ice-rich permafrost thaw can create expanding thermokarst lakes as well as shrinking large wetlands. Such processes can have major biogeochemical implications and feedbacks to climate systems by altering the pathways and rates of permafrost carbon release. We developed a new model parameterization that allows a direct representation of surface water dynamics with subsidence. Our results show increased surface water fractions around western Siberian plains and northeastern territories of Canada.