Continuous permafrost environments exhibit complex thermo-hydraulic interactions, particularly where lake taliks form within fractured bedrock aquifers. In the context of mining environments, these unfrozen zones act as preferential pathways for groundwater flow and solute transport, influencing permafrost integrity, discharge of deep brine, and mine water management. Based on a conceptual model of the Meliadine mine site, Nunavut, this research applies a two-stage numerical modelling framework using HEATFLOW/SMOKER, a finite-element platform that solves the coupled equations for density-dependent groundwater flow, advective–conductive heat transport with phase change, and advective-dispersive brine transport. The first stage focuses on calibrating a 2D cross-sectional model derived from a detailed geological framework, incorporating Quaternary sediments, fractured bedrock, and mapped fault zones. Calibration targets include hydraulic heads from multilevel piezometers, thermistor-based temperature profiles, and depth-dependent salinities reaching up to 64 g/L at the model base. In a second phase, scenario-based simulations explore the drivers of thaw front migration, talik development and brine transport. The results provide process-level insights into cryo-hydrogeological feedbacks, supporting predictive assessments for groundwater management and permafrost thaw mitigation in Arctic mining contexts. Future simulation scenarios will examine the impacts of climate forcing and mine dewatering on talik evolution.