When rivers collide, complex three-dimensional coherent flow structures are generated along the confluence’s mixing interface. These structures play important roles in mixing streamborne pollutants and suspended sediment and have considerable bearing on the morphology and habitat quality of the postconfluent reach. A particular structure of interest - streamwise orientated vortices (SOVs) - were first detected in numerical simulations to form in pairs, one on each side of the mixing interface rotating in the opposite sense of the other. Since, it has proven difficult to detect SOVs in situ with conventional pointwise velocimetry instrumentation. Despite the lack of clear evidence to confirm their existence, SOVs are nevertheless considered important drivers of mixing and sediment transport processes at confluences. Additionally, their causal mechanisms are also not fully known which hinders a complete conceptual understanding of these processes. To address these gaps, we analyze observations of strongly coherent SOVs filmed in aerial drone video of a mesoscale confluence with a stark turbidity contrast between its tributaries. Eddy-resolved modeling demonstrates the SOVs’ dynamics could only be accurately reproduced when a density difference (Δρ) was imposed between the tributaries (Δρ = 0.5 kg/m$^{3}$) – providing compelling evidence the observed SOVs are indeed a density-driven class of SOV. This work confirms that SOVs exist, expands understanding of their generative processes and highlights the important role of small density gradients (e.g., less than 0.5 kg/m3) on river confluence hydrodynamics.

A small gradient in the densities (Δρ) of two rivers was recently shown to develop coherent streamwise orientated vortices (SOVs) in the mixing interface of their confluence. We further investigate this phenomenon at the Coaticook and Massawippi confluence (Quebec, Canada) using eddy-resolved numerical modelling to examine how the magnitude and direction of Δρ; affect this secondary flow feature. Results show that a front from the denser channel always slides underneath the lighter channel independent of the direction of Δρ. When the fast tributary (Coaticook) is denser, coherent clockwise rotating density SOVs tend to form on the slow (Massawippi) side. However, when the slow Massawippi is denser by the same magnitude, anticlockwise secondary flow caused principally by shear induced interfacial instabilities develop on the fast Coaticook side. This shows the inertia of the tributary opposing the lateral propagation of the dense front shapes the secondary flow characteristics of the mixing interface. Moreover, in the absence of a density difference, anticlockwise SOVs are predicted by the model which correspond well to new aerial observations of anticlockwise SOVs at the site. A densimetric Froude number (Fd) convention accounting for the direction of Δρ is proposed to accurately convey the local inertial forces that oppose the lateral propagation of the dense front. Finally, a conceptual model of the mixing interface’s secondary flow structure over a spectrum of plausible Fd values is proposed. The Fd convention provides a flexible and consistent metric for use in future studies examining the effects of Δρ on river confluence hydrodynamics.