Along with the development in arctic region, the icebreaking technologies are gradually becoming the focus. As one of the most powerful and effective way to breaking ice, especially in the ability to solve ice jams, the study of the behaviour of the sea and river ice under dynamic loads is an urgent subject of scientific research and it attracts extensive attention. In addition, the temperature change in the process of ice failure cannot be neglected since that temperature plays an important role in the mechanical properties of the ice. In this study, a fully coupled thermoelastic ordinary state-based Peridynamic model is employed to investigate fragmentation of ice cover subjected to an underwater explosion. Both the deformation effect on the thermal effects and the thermal effects on deformation are taken into consideration. The pressure shocks generated by the underwater explosion are applied to the bottom surface of the ice cover for non-uniform load distributions. Crack propagation paths are investigated, the damage is predicted and compared with existing experimental results. The corresponding temperature distributions are also examined. Furthermore, the ice failure mode in both the top surface and the bottom surface of the ice sheet is investigated.