Hydrologic exchange is a critical mechanism that shapes hydrological and biogeochemical processes along a river corridor. Because of limitations in field accessibility, computational demand, and complexities of geomorphology and subsurface geology, full threeâdimensional modelling studies to quantify hydrologic exchange fluxes (HEFs) have been limited mostly to localâscale applications. At reach scales, although surface flow conditions and subsurface physical properties are wellâknown factors that modulate hydrologic exchanges, quantitative measures that can describe the effects of these factors on the strength and direction of such exchanges do not exist. To address this issue, we developed a oneâway coupled surface and subsurface water flow model using the commercial computational fluid dynamics (CFD) software STARâCCM+ and applied it to simulate HEFs in a 7âkm long reach along the main stem of the Columbia River in the United States. The model was validated against flow velocity measurements from an acoustic Doppler current profiler in the river, vertical HEFs estimated from a set of temperature profilers installed across the riverbed, and simulations from a reactive transport model. The validated model then was employed to systematically investigate how HEFs could be influenced by surface water fluid dynamics, subsurface structures, and hydrogeological properties. Our results suggest that reachâscale HEFs are dominated primarily by the thickness of the riverbed alluvium layer, and then by the alluvium permeability, the depth of the underlying impermeable layer, and the pressure boundary condition. Our results also elucidate the scale dependence of HEFs on fluid dynamics that can be captured only by threeâdimensional CFD models. That is, while the net HEFs over the entire 7âkm domain are not significantly influenced by surface water dynamics pressure, the dynamic pressure induced by fluid dynamics can lead to more than 15% in net HEFs for a river section of a few hundred metres.