Seasonal thermal storage systems can reduce the temporal mismatch between renewable energy availability and energy demand. Ice storage systems exhibit a non-linear behaviour in the heat exchange and storage processes, complicating the formulation of optimal design and operation problems. In this work, we propose a mixed-integer quadratically-constrained programming formulation, which minimizes the Levelized Cost of Energy for space heating and cooling, including sizing of a supporting PV array. The optimization was repeated for different storage volumes, finding the system optimal operation in each case –and thereby the optimal system sizing. The heating and cooling demands were computed from an archetypal office building, placed in three reference locations with cold and semi-arid, warm and humid continental, and temperate and humid continental climates. Results show that the optimal PV size decreases with growing ice storage volume, and an ice storage works best in a temperate continental climate, covering up to 47% of the cooling demand with a 250 m³ storage.