“…where ρ b is the bulk density of the subsurface medium (mg m −1 ), n is porosity (), and K d is the equilibrium sorption distribution coefficient (m 3 mg −1 ). The x denotes the spatial variability of K d in the x, y, and z directions, which results in spatial variability in R. More robust and mechanistic models (e.g., non-linear isotherms, Surface Complexation Models) could be used to model sorption (e.g., [67,68]), but our goal was to address the lack of understanding of how simple sorption processes affect solute transport within the hyporheic zone. Such knowledge is critical in designing remediation strategies for solutes that undergo sorption (e.g., perfluoroalkyl and polyfluoroalkyl substances (PFAS), 1,4-dioxane, perchloroethylene, trichloroethylene, and trace metals) as they are transported through porous media, particularly within the heterogeneous riverbed.…”