Glacial quarrying stems from the fracturing of subglacial bedrock. Much evidence shows that subcritical crack propagation of bedrock is closely related to subglacial water pressure fluctuations. Here we employ a model that assesses the impact of subglacial water pressure fluctuation on cavity length and subcritical crack propagation, while analyzing the effect of a pre-existing crack location using a phase-field model (PFM). Our results indicate that the cavity length is reduced during diurnal fluctuations in water pressure. There are two patterns of subcritical crack propagation on the corner of the step. The first stems from a rapid drop in water pressure. The second occurs after the water pressure recovers from the fluctuation to the initial steady state. This pattern is a consequence of enhanced stress concentration on the step since the modeled cavity length exceeds its steady value and has higher efficiency in promoting subcritical crack propagation. Additionally, based on the PFM results, we speculate that the subcritical crack initiation and propagation happen on a broader scale, including the ice-bed contact region and its adjacent region. Our findings imply that the duration of subcritical crack propagation is short and typically ceases once the cavity length adjusts to reduced water pressure levels.