Abstract: | This paper presents key challenges in modeling water quality processes of riparian ecosystems: How can the spatial and temporal extent of water and solute mixing in the riparian zone be modeled? What level of model complexity is justified? How can processes at the riparian scale be quantified? How can the impact of riparian ecosystems be determined at the watershed scale? Flexible models need to be introduced that can simulate varying levels of hillslope‐riparian mixing dictated by topography, upland and riparian depths, and moisture conditions. Model simulations need to account for storm event peak flow conditions when upland solute loadings may either bypass or overwhelm the riparian zone. Model complexity should be dictated by the level of detail in measured data. Model algorithms need to be developed using new macro‐scale and meso‐scale experiments that capture process dynamics at the hillslope or landscape scales. Monte Carlo simulations should be an integral part of model simulations and rigorous tests that go beyond simple time series, and point‐output comparisons need to be introduced. The impact of riparian zones on watershed‐scale water quality can be assessed by performing simulations for representative hillsloperiparian scenarios. |