Phosphorus export dynamics and hydrobiogeochemical controls across gradients of scale, topography and human impact


Ali, G., Wilson, H., Elliott, J., Penner, A., Haque, A., Ross, C., Rabie, M. (2017). Phosphorus export dynamics and hydrobiogeochemical controls across gradients of scale, topography and human impact, 31(18), 3130-3145.

Plain language summary

A variety of mathematical models were applied to describe how concentration of phosphorus in stream water relates to flow rate in a wide range of watersheds in the Lake Winnipeg Basin. A variety of different patterns were observed, but most sites showed increases in concentration with increasing flow. More complex patterns were observed in high slope or landscapes with more wetlands. This indicates that important factors aren't described in existing models and new models are needed that describe the chemistry and hydrology of more complex agricultural watershed on the Prairies.


Copyright © 2017 John Wiley & Sons, Ltd. Concentration-discharge (c-Q) plots are routinely used as an integrated signal of watershed response to infer solute sources and travel pathways. However, the interpretation of c-Q data can be difficult unless these data are fitted using statistical models. Such models are frequently applied for geogenic solutes, but it is unclear to what extent they might aid in the investigation of nutrient export patterns, particularly for total dissolved phosphorus (TDP) which is a critical driver of downstream eutrophication problems. The goal of the present study was therefore to statistically model c-Q relations (where c is TDP concentrations) in a set of contrasting watersheds in the Northern Great Plains—ranging in size from 0.2 to 1000+ km2—to assess the controls of landscape properties on TDP transport dynamics. Six statistical models were fitted to c-Q data, notably (a) one linear model, (b) one model assuming that c-Q relations are driven by the mixing of end-member waters from different landscape locations (i.e., hydrograph separation), (c) one model relying on a biogeochemical stationarity hypothesis (i.e., power law), (d) one model hypothesizing that c-Q relations change as a function of the solute subsurface contact time (i.e., hyperbolic model), and (e) two models assuming that solute fluxes are mostly dependent on reaction rates (i.e., chemical models). Model performance ranged from mediocre (R2 < 0.2) to very good (R2 > 0.9), but the hydrograph separation model seemed most universal. No watershed was found to exhibit chemostatic behaviour, but many showed signs of dilution or enrichment behaviour. A tendency toward a multi-model fit and better model performance was observed for watersheds with moderate slope and higher effective drainage area. The relatively poor model performance obtained outside these conditions illustrates the likely importance of controls on TDP concentrations in the region that are independent of flow dynamics.

Publication date


Author profiles