Haulage networks are vital to underground mining operations as they constitute the arteries through which blasted ore is transported to surface.In the sublevel stoping method and its variations,haulage drifts are excavated in advance near the ore block that will be mined out.Numerical modeling is a technique that is frequently employed to assess the redistribution of mining-induced stresses,and to compare the impact of different stope sequence scenarios on haulage network stability.In this study,typical geological settings in the Canadian Shield were replicated in a numerical model with a steeplydipping tabular orebody striking EW.All other formations trended in the same direction except for two dykes on either side of the orebody with a WNW-ESE strike.Rock mass properties and in situ stress measurements from a case study mine were used to calibrate the model.Drifts and crosscuts were excavated in the footwall and two stope sequence scenarios-a diminishing pillar and a center-out one-were implemented in 24 mining stages.A combined volumetric-numerical analysis was conducted for two active levels by comparing the extent of unstable rock mass at each stage using shear,compressive,and tensile instability criteria.Comparisons were made between the orebody and the host rock,between the footwall and hanging wall,and between the two stope sequence scenarios.It was determined that in general,the center-out option provided a larger volume of instability with the shear criterion when compared to the diminishing pillar one (625,477 m3 compared to 586,774 m3 in the orebody;588 m3 compared to 403 m3 in the host rock).However,the reverse was true for tensile (134,298 m3 compared to 128,834 m3 in the orebody;91,347 m3 compared to 67,655 m3 in the host rock)instability where the diminishing pillar option had the more voluminous share.