The South Saskatchewan River (SSR) is one of the most important river systems in Saskatchewan and, arguably, in Canada. Most of the Saskatchewan residents, industries, and powerplants depend on the SSR for their water requirements. An established 1D modelling approach was chosen and coupled with the Hydrologic Engineering Center’s River Analysis System (HEC-RAS). The WASP (Water Quality Analysis Simulation Program) stream transport module, TOXI, is coupled with flow routing for free-flow streams, ponded segments, and backwater reaches and is capable of calculating the flow of water, sediment, and dissolved constituents across branched and ponded segments. Copper and nickel were chosen as two metals with predominantly anthropogenic (agriculture, mining, and municipal and industrial waste management) and geogenic (natural weathering and erosion) sources, respectively. Analysis was carried out at ten different sites along the South Saskatchewan River, both upstream and downstream of the City of Saskatoon, in the years 2020 and 2021. Model performance was evaluated by comparing model predictions with concentrations of copper and nickel measured in a previously published study. The model performed well in estimating the concentrations of copper and nickel in water samples and worked reasonably well for sediment samples. The model underestimated the concentration values at certain segments in both water and sediment samples. In order to calibrate the model more accurately, extra diffusive contaminant loads were added. While several default parameter values had to be used due to the unavailability of primary historical data, our study demonstrates the predictive power of combining WASP—TOXI and HEC-RAS models for the prediction of contaminant loading. Future studies, including those on the impacts of global climate change on water quality on the Canadian prairies, will benefit from this proof-of-concept study.

The lower Athabasca River was used as a test case using total suspended sediment, chloride and vanadium as the model variables. Upstream model boundary conditions included water from the tributary Clearwater River (right stream tube) and the upper Athabasca River extending upstream of the tributary mouth (left stream tube). This model will be extended to include the Peace-Athabasca Delta (PAD), to determine the implications of mining outfall discharges on a large region of the Athabasca – PAD region. A novel, quasi-two-dimensional surface water-quality modelling approach is presented in which the model domain can be discretised in two dimensions, but a one-dimension solver can still be applied to capture water flow between the discretisation units (segments). The approach requires a river reach to be divided into two stream tubes, along the left and right river sides, with flows exchanging through the segments longitudinally and also laterally between adjacent segments along the two streams. The new method allows the transverse mixing of tributary and outfall water of different constituent concentrations to be simulated along the course of the river. Additional diffuse loading of dissolved vanadium could be determined from the model’s substance balance. A scenario was then simulated in which the transport and fate of vanadium in a floodplain lake and a secondary channel was determined. • Quasi-2D modelling approach proves to be viable for transverse mixing. • Quasi-2D approach allows secondary channels and side lakes to be modelled. • Quasi-2D approach is appropriate to scale up to entire lower Athabasca River reach. • The approach allowed a diffuse loading of dissolved vanadium to be quantified.

A quasi-two-dimensional (quasi-2D) modelling approach is introduced to mimic transverse mixing of an inflow into a river from one of its banks, either an industrial outfall or a tributary. The concentrations of determinands in the inflow vary greatly from those in the river, leading to very long mixing lengths in the river downstream of the inflow location. Ideally, a two-dimensional (2D) model would be used on a small scale to capture the mixing of the two flow streams. However, for large-scale applications of several hundreds of kilometres of river length, such an approach demands too many computational resources and too much computational time, especially if the application will at some point require ensemble input from climate-change scenario data. However, a one-dimensional (1D) model with variables varying in the longitudinal flow direction but averaged across the cross-sections is too simple of an approach to capture the lateral mixing between different flow streams within the river. Hence, a quasi-2D method is proposed in which a simplified 1D solver is still applied but the discretisation of the model setup can be carried out in such a way as to enable a 2D representation of the model domain. The quasi-2D model setup also allows secondary channels and side lakes in floodplains to be incorporated into the discretisation. To show proof-of-concept, the approach has been tested on a stretch of the lower Athabasca River in Canada flowing through the oil sands region between Fort McMurray and Fort MacKay. A dye tracer and suspended sediments are the constituents modelled in this test case.