Literature from the past two decades demonstrates the feasibility of utilizing borehole geothermal energy storage (BTES) system for the heating of buildings in the cold climate region like Northern Quebec. However, BTES systems would generate an increase in temperature in surrounding soil formations, which may induce ground deformation and result in unexpected accidents. This study investigates the effect of BTES systems of 50-year service life period on thermal consolidation of surrounding soil formations in northern Quebec, where BTES sites are covered by a large quantity of unconsolidated glacial tills. A fully coupled thermal-hydro-mechanical modeling is conducted using Abaqus, and cases with different configuration and operation temperature are modeled. Thermally induced pore pressure generation and dissipation processes are simulated and analyzed. Thermally induced deformation of soil formations is also addressed. The glacial till in the upper layer is prone to sustain strain softening and the glacial till in the lower layer is prone to sustain strain hardening during the periodical thermal operation. A BTES system with a large geometric scale or with a high operation temperature (60°C in this study) would impose significant influence on the glacial till formation, which is displayed by significant changes in pore water pressure and ground deformations.