Simulations of tracer experiments conducted with a three-dimensional primitive-equation hydrodynamic and transport model are used to understand the processes controlling the rate of mixing between two rivers (Ebro and Segre), with distinct physical and chemical properties, at their confluence, upstream of a meandering reservoir (Ribarroja reservoir). Mixing rates downstream of the confluence are subject to hourly scale oscillations, driven partly by changes in inflow densities and also as a result of turbulent eddies that develop within the shear layer between the confluent rivers and near a dead zone located downstream of the confluence. Even though density contrasts are low-at most O(10 21 ) kg m 23 difference among sources-and almost negligible from a dynamic point of view-compared with inertial forces-they are important for mixing. Mixing rates between the confluent streams under weakly buoyant conditions can be of up to 40% larger than those occurring under neutrally buoyant conditions. The buoyancy effects on mixing rates are interpreted as the result of changes in the contact area available for mixing (distortion of the mixing layer). For strong density contrasts, though, when the contact area between the streams becomes nearly horizontal, larger density differences between streams will lead to weaker mixing rates, as a result of the stabilizing effect of vertical density gradients.