The ice generation is one of the challenges facing the methane hydrate depressurization, which, however, has not been fully addressed by existing numerical models for hydrate-bearing sediments (HBS). In this study, we develop a high-fidelity, fully coupled thermo-hydro-mechanical-chemical numerical model that incorporates the effect of ice. The model, developed using COMSOL, takes into account water–ice phase change, thermally...