In the higher latitudes of the northern hemisphere, ice jam related flooding can result in millions of dollars of property damages, loss of human life and adverse impacts on ecology. Since ice-jam formation mechanism is stochastic and depends on numerous unpredictable hydraulic and river ice factors, ice-jam associated flood forecasting is a very challenging task. A stochastic modelling framework was developed to forecast real-time ice jam flood severity along the transborder (New Brunswick/Maine) Saint John River of North America during the spring breakup 2021. Modélisation environnementale communautaire—surface hydrology (MESH), a semi-distributed physically-based land-surface hydrological modelling system was used to acquire a 10-day flow forecast. A Monte-Carlo analysis (MOCA) framework was applied to simulate hundreds of possible ice-jam scenarios for the model domain from Fort Kent to Grand Falls using a hydrodynamic river ice model, RIVICE. First, a 10-day outlook was simulated to provide insight on the severity of ice jam flooding during spring breakup. Then, 3-day forecasts were modelled to provide longitudinal profiles of exceedance probabilities of ice jam flood staging along the river during the ice-cover breakup. Overall, results show that the stochastic approach performed well to estimate maximum probable ice-jam backwater level elevations for the spring 2021 breakup season.