2017-08-25T16:31:17

This study develops a general quantification framework for consistent intermodel and intercatchment comparison of the nutrient and pollutant mass loading from multiple sources in a catchment area to downstream surface and coastal waters. The framework accounts for the wide spectrum of different transport pathways and travel times through the subsurface (soil, groundwater, sediment) and the linked surface (streams, lakes, wetlands) water systems of a catchment. The account is based on key flow partitioning and mass delivery fractions, which can be quantified differently by different flow and transport and reaction models. The framework application is exemplified for two Swedish catchment cases with regard to the transport of phosphorus and of a generic attenuating solute. The results show essential differences in model quantifications of transport pathways and temporal spreading, with important implications for our understanding of cause and effect in the catchment-scale nutrient and pollutant loading to downstream waters.

This study has investigated and outlined the possible quantification and mapping of the distributions of advective solute travel times through hydrological catchments. These distributions are essential for understanding how local water flow and solute transport and attenuation processes affect the catchment-scale transport of solute, for instance with regard to biogeochemical cycling, contamination persistence and water quality. The spatial and statistical distributions of advective travel times have been quantified based on reported hydrological flow and mass-transport modeling results for two coastal Swedish catchments. The results show that the combined travel time distributions for the groundwater-stream network continuum in these catchments depend largely on the groundwater system and model representation, in particular regarding the spatial variability of groundwater hydraulic parameters (conductivity, porosity and gradient), and the possible contributions of slower/deeper groundwater flow components. Model assumptions about the spatial variability of groundwater hydraulic properties can thus greatly affect model results of catchment-scale solute spreading. The importance of advective travel time variability for the total mass delivery of naturally attenuated solute (tracer, nutrient, pollutant) from a catchment to its downstream water recipient depends on the product of catchment-average physical travel time and attenuation rate.

Abstract. This paper quantifies and maps the effects of coupled physical and biogeochemical variability on diffuse hydrological mass transport through and from catchments. It further develops a scenario analysis approach and investigates its applicability for handling uncertainties about both physical and biogeochemical variability and their different possible cross-correlation. The approach enables identification of conservative assumptions, uncertainty ranges, as well as pollutant/nutrient release locations and situations for which further investigations are most needed in order to reduce the most important uncertainty effects. The present scenario results provide different statistical and geographic distributions of advective travel times for diffuse hydrological mass transport. The geographic mapping can be used to identify potential hotspot areas with large mass loading to downstream surface and coastal waters, as well as their opposite, potential lowest-impact areas within the catchment. Results for alternative travel time distributions show that neglect or underestimation of the physical advection variability, and in particular of those transport pathways with much shorter than average advective solute travel times, can lead to substantial underestimation of pollutant and nutrient loads to downstream surface and coastal waters. This is particularly true for relatively high catchment-characteristic product of average attenuation rate and average advective travel time, for which mass delivery would be near zero under assumed transport homogeneity but can be orders of magnitude higher for variable transport conditions. A scenario of high advection variability, with a significant fraction of relatively short travel times, combined with a relevant average Correspondence to: K. Persson (klas.persson@natgeo.su.se) biogeochemical mass attenuation rate, emerges consistently from the present results as a generally reasonable, conservative assumption for estimating maximum diffuse mass loading, when the prevailing physical and biogeochemical variability and cross-correlation are uncertain.