1 Introduction

Since the early twentieth century, rapid population growth, intensification of agriculture and extensive urbanisation have led to widespread phosphorus (P) enrichment of aquatic environments (Cordell et al. 2009) resulting in an array of detrimental economic (Pretty et al. 2003) and ecological impacts (Hilton et al. 2006; Jarive et al., 2018). As a naturally limiting nutrient of plant growth in waterbodies, P-enrichment fuels blooms of phytoplankton, periphyton and neuro-toxin secreting cyanobacteria colonies, which can dramatically lower species diversity and lead to a fundamental breakdown of ecosystem functioning (Smith et al. 1999).

Wastewater effluent discharged into rivers from sewage treatment works (STWs) represents one of the most important point sources of P pollution and is a major driver of freshwater eutrophication (Jarvie et al. 2006). STWs discharge effluent continuously with high-P concentrations and because the majority (70–90%) of the P is released as biologically available soluble reactive phosphorus (SRP), it carries increased risk of initiating deleterious changes in ecological communities (Demars et al. 2005; Edwards and Withers 2008). Consequently, under the EU Urban Wastewater Directive (91/271/EC), STWs serving a population of 10,000–100,000 people have water quality restrictions set at 2 mg P L−1 for effluent discharging into a surface watercourse, with this decreasing to 1 mg P L−1 for facilities serving > 100,000 people. However, smaller facilities (< 10,000 people) have no legal requirement to meet such targets and P concentrations in effluent discharges of up to 20 mg P L−1 have been recorded (Neal et al. 2005; Jarvie et al. 2006; Withers and Jarvie 2008). Previous water quality monitoring studies have therefore observed increases in dissolved P concentrations of > 600% downstream of STW outflows (House and Denison 2002; Demars and Harper 2005), with concentration orders of magnitude greater than the 40–120 μg P L−1 concentration recommended under the EU Water Framework Directive (2000/60/EC) to achieve ‘good’ ecological and chemical status.

Through various physical, chemical and biological processes, river systems have the internal capacity to cycle P between the sediments on the riverbed and the overlying water column and to transform P between organic, inorganic, particulate and dissolved forms (Jarvie et al. 2005). This cycling includes biotic and abiotic transformations related to P uptake/release, adsorption/desorption, precipitation/dissolution and advection/diffusion (Withers and Jarvie 2008). Importantly, this P cycling can create a self-regulating buffering system whereby riverbed sediments act to buffer SRP dissolved in the water column downstream of STWs (Wang and Li 2010). In this respect, the Equilibrium Phosphorus Concentration (EPC0) is a measure to determine whether riverbed sediments are acting as net sources or sinks of SRP. The EPC0 is defined as the concentration of SRP in solution which, when placed in contact with sediment, produces no net sorption or desorption of SRP over a period of 24 h (House et al. 1995; House and Denison 2000). When SRP > EPC0, the sediment has the potential to take up SRP from the water column and when SRP < EPC0, the sediment has the potential to release SRP into the water column (Palmer-Felgate et al. 2009). It is this differential between the EPC0 of the sediment and the SRP concentration of the bed sediment-water interface layer (defined to a depth of 0.1 m) that determines the uptake or release of SRP from riverbed sediments (Jarvie et al. 2005).

Whether riverbed sediment acts as a source or a sink for SRP is primarily related to the redox conditions and physicochemical properties of the sediment as this determines the bed sediment P-sorption behaviour (Palmer-Felgate et al. 2009). In particular, the particle size, mineralogy and the concentration of exchangeable phosphate already sorbed onto the sediment surface are key factors which influence the speciation of P in the sediment, as well as the potential for SRP exchange across the sediment-water interface (Stone et al. 1995; Zhang and Huang 2007).

Due to the continuous nature of STW effluent discharge, P concentrations downstream of STWs tend to display a highly seasonal pattern with higher concentrations during the summer low flows and lower concentrations during winter high flows due to dilution, thus causing P concentrations to peak during ecologically sensitive summer season when eutrophication risk is greatest (Cooper et al. 2002; Withers and Jarvie 2008; Bowes et al. 2012). The sorption and storage of SRP by riverbed sediments downstream of STWs can therefore play a vital role in reducing eutrophication risk during the ecologically sensitive low-flow summer period (Jarvie et al. 2006). However, because riverbed sediment chemistry is a major controlling factor on the chemistry of the overlying waterbody, deposition of P-rich fines can also act as a localised source of SRP if concentrations are greater than in the overlying waterbody or if sediments become entrained by biotic mixing and storm events (Jarvie et al. 2005; van der Perk et al. 2007). Increased P concentrations on the riverbed can also fuel excessive growth of macrophytes and benthic algae which are able to extract nutrients directly from the sediment rather than relying on dissolved SRP (Hilton et al. 2006). Understanding P cycling and the balance of these processes within ecosystems is therefore an important aspect of P transport which should be adequately represented within catchment management plans. Research into EPC0 values can assist by providing a clearer insight into abiotic cycling processes related to P uptake/release within fluvial systems.

The aim of this study was to investigate how the presence of STWs influences the EPC0 of riverbed sediments and dissolved SRP concentrations during the summer and early autumn low-flow periods of high-eutrophication risk in the River Wensum, UK—a river with multiple conservation designations suffering from elevated P concentrations and containing a large number of STWs. The main objectives were as follows:

  1. (i)

    Determine if riverbed sediments are effective at buffering dissolved SRP downstream of STWs;

  2. (ii)

    Contrast tributaries with and without a STW to assess where sediments are acting as net sources or sinks for P;

  3. (iii)

    Assess whether the ability of sediment to uptake excess P varies temporally.

To our knowledge, there have been no direct measurements of riverbed sediment EPC0 values in the River Wensum during the past decade and it is envisaged that this study will assist in providing an improved management plan for nutrient enrichment across the wider Wensum catchment.

2 Materials and methods

2.1 Study location

The River Wensum, Norfolk, UK, is a 78-km-length groundwater-dominated lowland (source = 75 m above sea level) calcareous river that drains an area of 660 km2 and has a mean annual discharge of 4.1 m3 s−1 near its outlet (Fig. 1). The river is designated a Site of Special Scientific Interest (SSSI) and a European Special Area of Conservation (SAC) in recognition of it being one of the best examples of a lowland calcareous river in the world due to the diversity of its internationally important flora and invertebrate fauna. However, 99.4% of the river habitat is considered to be in an unfavourable and declining state due to poor morphology (e.g. channelized, disconnected from the floodplain), sedimentation and eutrophication, whilst 27% of waterbodies across the catchment are at risk of failing P standards (Sear et al. 2006). There are 21 STWs across the catchment serving a population of c.57,000 people of which seven of these have tertiary P stripping technologies installed to reduce nutrient discharges (Evans 2012). Previous modelling studies using SAGIS and SEPARATE have estimated STWs are responsible for 26–47% of P in the River Wensum, compared with 29–40% derived from agriculture and 20–32% from other sources including urban runoff, storm tanks, atmospheric deposition and riverbank erosion (Natural England 2015). In addition, there are an estimated 1863 properties off mains sewage served by domestic septic tank systems, estimated to contribute a further 2–4% of the annual P load (Natural England 2015).

Fig. 1
figure 1

Distribution of sampling locations and sewage treatment works across the River Wensum catchment, Norfolk, UK

The catchment is underlain by Cretaceous White Chalk bedrock which is unconfined in the upper catchment and along sections of the river valley where the baseflow index (BFI) is 0.7–0.9. Over much of the rest of the catchment, the Chalk is confined by superficial deposits of Mid-Pleistocene diamicton glacial tills principally comprising chalky, flint-rich boulder clays of the Sheringham Cliffs (~ 0.2–10 m depth) and Lowestoft (~ 10–20 m depth) formations (Fig. S1, Electronic Supplementary Material). These are interspersed with layers of glaciofluvial and glaciolacustrine sands and gravels where the BFI is 0.5–0.7. Within the river valley, Holocene-age alluvium and river terrace deposits are present (Hiscock et al. 1996; Lewis 2014). Surface soils across the catchment range from low-permeability clay loams and sandy peats, to free draining sandy loams (Fig. S2, Electronic Supplementary Material). Arable agriculture (wheat, barley, sugar beet, oilseed rape) dominates land use (63%) with the remainder comprising 19% improved grassland, 9% mixed woodland, 5% unimproved grassland and 4% urban. The mean annual temperature is 10.1 °C and the mean annual rainfall total is 674 mm (1981–2010) (Met Office 2017).

2.2 Field campaign

River water and riverbed sediment samples were collected from ten tributary outlets across the catchment (Fig. 1) at monthly intervals between July and October 2016, such that 40 sediment and 40 water samples were collected in total during the field campaign. Of these ten locations, five were located downstream of STWs (sites 1, 6, 8, 14, 15) and five were on tributaries without a STW (sites 2, 9, 16, 20, 21). Sites influenced by STWs were situated > 500 m downstream of the facility to ensure the effluent discharge would be well mixed by the time it reached the sampling location. Riverbed sediments (~ 500 g) were collected from the upper 10 cm of the riverbed using a trowel and were sieved down to 2 mm on site to remove coarse sands and gravels which have limited P storage/release potential. River water was grab-sampled from the centre of the channel at ~ 0.5 × river depth in 1 L polypropylene bottles. All samples were returned to cold storage (4 °C) within 4 h to minimise biological degradation, and all analysis was completed within 7 days. Samples collected in July (1.50 m3 s−1), August (1.20 m3 s−1) and October (1.46 m3 s−1) were taken under baseflow conditions (based on the hydrograph of the site 8 gauging station) with limited rainfall (4, 7 and 14 mm, respectively) occurring in the preceding 5 days (Fig. 2). Higher discharges (1.83 m3 s−1) were observed during September reflecting higher rainfall totals (35 mm) in the 5 days before sampling.

Fig. 2
figure 2

Annual hydrograph for the River Wensum recorded at the Swanton Morley gauging station (sampling site 8) during the period 1969–2016 (CEH 2018). Precipitation record is for 2016. Vertical lines refer to sampling dates in 2016

2.3 Laboratory analysis

Riverbed sediment EPC0 values were derived following the methods of Jarvie et al. (2005). Approximately 2 g of each wet sediment sample was oven-dried (110 °C) for 18 h and weighed to determine water content. For each sediment subsample, wet sediment (equivalent to 0.5 g dry weight) was added into six 250 mL volume centrifuge tubes along with 200 mL of chilled CaCl2 solution (220 mg L−1). This solution was chosen to simulate similar calcium concentrations and conductivity to the calcareous waters of the River Wensum (~ 120 mg Ca L−1). The six centrifuge tubes were individually spiked with a phosphate standard (KH2PO4) to give concentrations of 0, 50, 100, 250, 500 and 1000 μg P L−1 and shaken in an incubated (10 °C) orbital shaker (150 rpm, 24 h) to ensure all sediment was in contact with the solution. After shaking, samples were centrifuged (6000 rpm, 5 min) and the supernatant extracted and filtered through a 0.45-μm syringe filter. Fifty milliliter of each filtered sample was added to separate 100 mL conical flasks to which 5 mL of mixed reagent (ammonium molybdate, sulphuric acid, ascorbic acid, potassium antimonyl tartrate) was added. After 30 min, 2 mL of each sample was transferred to cuvettes and analysed colorimetrically to determine SRP concentration using a spectrophotometer (885 nm). River water samples were similarly filtered through a 0.45-μm syringe filter and 50 mL added to 5 mL of mixed reagent prior to analysing colorimetrically for SRP concentration.

2.4 EPC0 calculations

The EPC0 was derived by first calculating the change in SRP concentration (ΔC in μg P L−1) of the solution after 24 h incubation with the riverbed sediment whilst correcting for any differences in sediment mass and solution volume (Eq. (1)):

$$ \Delta C=\frac{V}{m}\ {C}_0-{C}_{24} $$
(1)

where V is the total volume of the solution (L); m is the mass of dry sediment (g); C0 is the concentration of the SRP spikes (μg P L−1) at time zero and C24 is the SRP concentration of the solution after being shaken for 24 h (μg P L−1). ΔC was then plotted against the SRP concentration in solution after 24 h (C24), and the data points fitted with a non-linear least squares Freundlich isotherm (Eq. (2); House and Denison 2000), which yields the EPC0 where it intercepts the x-axis (Fig. 3):

$$ \Delta C={K}_{\mathrm{f}}{C_{24}}^{\frac{1}{n}} $$
(2)

where Kf is a constant related to the affinity of sorption and n is the Freundlich isotherm constant.

Fig. 3
figure 3

Example Freundlich isotherm plot for Wensum site 2 used to estimate the equilibrium phosphorus concentration (EPC0) of riverbed sediment

3 Results and discussion

3.1 Impact of STWs on dissolved P concentrations

Concentrations of dissolved SRP were high across the River Wensum catchment, with mean concentrations for all ten sampling locations exceeding the recommended 40–120 μg P L−1 standard for good chemical status (Table 1). The highest and lowest SRP concentrations were recorded at site 1 (1136 μg P L−1) and site 9 (46 μg P L−1), respectively, both in October 2016. The overall mean SRP concentration for all sites combined was 284 μg P L−1 (σ = 193 μg P L−1), but there were pronounced statistically significant differences (t test, p < 0.05) between sites depending on whether they were impacted by effluent discharges from STWs (Fig. 4). Mean SRP concentrations for individual sites downstream of STWs ranged from 278 to 647 μg P L−1, whilst they ranged from 124 to 247 μg P L−1 at sites without a STW, meaning that all sites impacted by STW outflows had higher SRP concentrations than unaffected sites. However, there was no significant relationship (p = 0.269; R2 = 0.02) between SRP concentration and the distance downstream from the STW outflows. The overall mean SRP for impacted sites (382 μg P L−1) was double that of unaffected sites (185 μg P L−1).

Table 1 River water-soluble reactive phosphorus (SRP) concentrations and riverbed sediment equilibrium phosphorus concentrations (EPC0) for ten sites across the River Wensum catchment with or without a sewage treatment works (STW) for the period of July–October 2016
Fig. 4
figure 4

Average a river water-soluble reactive phosphorus concentrations and b riverbed sediment EPC0 values for ten sampling sites on the River Wensum with and without sewage treatment works (STWs). Error bars represent one standard error

To put these values in context with the natural background P concentrations for this region, a previous study on the neighbouring, hydrochemically similar, River Bure recorded a mean concentration of just 4 μg SRP L−1 on a section of the river which flowed through an entirely wooded sub-catchment without point or diffuse source P inputs from STWs or agriculture (Moss et al. 1988). Additionally, groundwater analysed from 6, 12, 15 and 50 m depth boreholes drilled in the Blackwater Drain sub-catchment of the River Wensum had mean total P concentrations of < 10 μg SRP L−1 (Lewis 2014).

Although the largest two STWs (Dereham and Fakenham) in the River Wensum catchment have P stripping facilities installed, such technology is absent from the 14 smaller rural STWs as it is not considered cost-effective to install for works severing relatively small populations (Bewes et al. 2012). However, these results clearly demonstrate that STWs are significantly elevating SRP concentrations above the levels observed in unaffected tributaries and are adding to diffuse P inputs from the surrounding agricultural land. Without the installation of P removal technologies (e.g. chemical precipitation, enhanced biological removal) at more STWs, P enrichment is likely to be amplified over the coming decades by an increasing population resulting in greater sewage effluent discharge across the Wensum catchment.

3.2 Riverbed sediments as sources or sinks of P

The EPC0 values varied widely across the sites from a low of 18 μg P L−1 at site 9 in September 2016, to a high of 189 μg P L−1 at site 20 in August (Table 1). The overall mean EPC0 value for all sites combined was 81 μg P L−1 (σ = 54 μg P L−1), but as with SRP concentrations, there were statistically significant differences (t test, p < 0.05) between sampling locations depending on the presence or absence of STWs (Fig. 4). The mean EPC0 for individual sites downstream of STWs ranged from 74 to 147 μg P L−1, whilst sites without a STW ranged from 32 to 117 μg P L−1. The overall mean EPC0 for STW impacted sites (105 μg P L−1) was 70% higher than that of unaffected sites (62 μg P L−1), likely reflecting the deposition of P-rich fines on the riverbed downstream of STWs. There was no significant relationship (p = 0.06; R2 = 0.17) between EPC0 and the distance downstream from the STW outflows.

Note, for site 1, EPC0 values could not be calculated as all ΔC values were negative and the Freundlich isotherm never intercepted the x-axis. This, combined with the very high SRP concentrations recorded at this site, suggests that the sediments here were fully saturated with P and were likely releasing SRP into the water column even at high-spike concentrations (i.e. 1000 μg P L−1). Visual observations indicate this was likely due to the riverbed at site 1 being composed of a much higher proportion of coarse gravel and cobbles with low-P sorption capacity than observed at the other sampling sites (Fig. S3, Electronic Supplementary Material).

Plotting SRP against EPC0 (Fig. 5) reveals that riverbed sediments at all ten sites, regardless of STW influence, were almost always acting as a sink for SRP from the overlying water column, with SRP concentrations consistently greater than the sediment EPC0 values. The only exception to this was during October at site 9 when the EPC0 (56 μg P L−1) was marginally higher than the river water SRP concentration (46 μg P L−1). The residual concentration (i.e. SRP—EPC0) revealed sites without a STW were significantly (p < 0.05) closer to the 1:1 equilibrium line (mean residual = 124 μg P L−1; σ = 74 μg P L−1) than STW-impacted sites (mean residual = 211 μg P L−1; σ = 119 μg P L−1). This means that even though locations downstream of STWs had significantly higher SRP and EPC0 concentrations, the increase in EPC0 (+ 70%) relative to non-impacted sites was proportionally less than the increase in SRP (+ 106%), and as such the sediments here had greater capacity to absorb more SRP before they reached equilibrium (assuming the sediment was not already saturated with P). Similar results have previously been reported for the River Avon and River Wye, UK (Jarvie et al. 2005), and they reveal the sediments downstream of STWs were better able to buffer the river against increases in dissolved SRP released in sewage effluent.

Fig. 5
figure 5

Uptake and release potential for soluble reactive phosphorus of River Wensum bed sediments for sites a with or without STWs and b during different months. Points below the 1:1 equilibrium line indicate potential SRP uptake by riverbed sediments, whilst points above the equilibrium indicate potential SRP release into the overlying water column

The buffering potential is confirmed by a linear regression model of SRP concentration against EPC0 for all data points (Fig. 6a). This reveals a weak (R2 = 0.20) but significant (p < 0.05) positive correlation, with a gradient < 1 indicating that the riverbed sediments were not in equilibrium with SRP in the overlying water column and therefore had the capacity for sorption when exposed to higher dissolved SRP concentrations. This correlation strengthens further (R2 = 0.50) when considering only the mean EPC0 and SRP values for each site (Fig. 6b), with this positive relationship revealing that riverbed sediments generally have higher P concentrations at sites exposed to higher dissolved SRP concentrations downstream of STWs.

Fig. 6
figure 6

Correlation plots of soluble reactive phosphorus concentration against riverbed sediment EPC0 for a all data, b site means and c monthly means

The fact that riverbed sediments were acting as sinks for SRP across the catchment, regardless of whether the site is influenced by a STW, indicates there are significant P inputs from other sources elevating dissolved SRP in the water column. Given the dominance of intensive agricultural land across the catchment, this is most likely due to diffuse pollution from phosphate fertilisers, with a previous study in the River Wensum catchment revealing 1–2% of total annual P fertiliser application is lost into the river network (Outram et al. 2016).

Ultimately, the ability of sediments to sorb P is determined by the chemical composition of the catchment bedrock and superficial geology which whether to produce the soils and sediments that determine the natural background P buffering capacity (Walling et al. 2001; van der Perk et al. 2007; Ballantine et al. 2008). As such, it might be expected that EPC0 values in the Wensum catchment would be higher in areas where fine glacial till is the dominant superficial deposit and values would be lower in areas covered with coarse sands and gravels. This is because P sorption can only occur across the water-sediment interface and, as such, the specific surface area (SSA) of sediments strongly influences the rates of surface reactions and the ability of sediments to adsorb P (Evans et al. 2004; Rawlins et al. 2010). Sediments dominated by fine-grained particles have higher SSAs and tend to contain greater concentrations of clay minerals and metal oxyhydroxide complexes to which P has a high affinity, thus enhancing sorption rates (Hartikainen et al. 2010; Cooper et al. 2015). However, here EPC0 values displayed no obvious relationship with either the bedrock/superficial geology (Fig. S1, Electronic Supplementary Material) or the soil type (Fig. S2, Electronic Supplementary Material) present across the catchment, with instead the presence or absence of STWs remaining the dominant discriminator of EPC0 values.

3.3 Temporal P dynamics

Over the duration of the study, mean SRP concentrations for all sites were highest in September (363 μg P L−1; σ = 173 μg P L−1), followed by October (299 μg P L−1; σ = 309 μg P L−1), July (255 μg P L−1; σ = 82 μg P L−1) and August (248 μg P L−1; σ = 127 μg P L−1) (Table 1). September also yielded the highest mean SRP concentrations when considering sites with (496 μg P L−1) and without (229 μg P L−1) STW influence, thus indicating these September peaks were related to external climatic factors rather than fluctuations in effluent discharge. Heavy rainfall (35 mm) in the days preceding the September sampling resulted in higher river flows (14% above the 1969–2016 average for this date) than observed during July (9% below average), August (14% below average) and October (40% below average) when baseflow conditions prevailed. Under these higher September flows, dilutions in SRP concentration would typically be expected if the majority of P originated from continuous point source STW release. The concentration increases observed here therefore strongly suggests diffuse P from agriculture was a more important source of P input into the River Wensum at the time of the September sampling due to precipitation-induced soil erosion and flushing of P sources from across the wider catchment (Quinton et al. 2001). This is supported by previous studies which have shown up to 70% of annual riverine P load is transported from arable dominated sub-catchments of the River Wensum during storm events at mean rates of 0.1–0.2 kg P ha−1 (Outram et al. 2014; Outram et al. 2016). It is also likely that the elevated September discharge initiated the remobilisation of P-rich fines deposited on the riverbed during intervening periods of low flow when sediment and P retention times would have been greater.

Mean EPC0 values of the riverbed sediments were similarly highest during September at sites with (118 μg P L−1; σ = 51 μg P L−1) and without (70 μg P L−1; σ = 79 μg P L−1) STWs, revealing elevated levels of particulate P on the riverbed at some, but not all, sites at this time. This probably indicates the deposition of P-rich soil eroded from the surrounding agricultural land during this heavy rainfall event. Despite this increase, the mean residual concentration in September (239 μg P L−1) was significantly (p < 0.05) larger than for the other 3 months (136–140 μg P L−1), indicating greater potential to buffer the river against the higher dissolved SRP concentrations (Fig. 5b). The lowest EPC0 values were recorded in October for sites with (90 μg P L−1; σ = 58 μg P L−1) and without (43 μg P L−1; σ = 10 μg P L−1) STWs. The relationship between the average monthly SRP concentrations and EPC0 values reveals a strong (R2 = 0.91) positive correlation (Fig. 6c), albeit with limited data points, with an SRP:EPC0 ratio of 5:1 indicating sediment during all months was acting as a net sink for SRP.

3.4 Significance and implications

Despite being conducted over just a 4-month period, the results presented here clearly demonstrate riverbed sediments are acting as sinks for SRP across the River Wensum catchment, regardless of the presence or absence of STWs. However, although the EPC0 values of riverbed sediments have been shown to increase with an increase in SRP, the elevated concentrations of SRP downstream of STWs clearly demonstrate that the sediments have insufficient P sorption capacity to be able to fully buffer the river against effluent discharge. This is due to the combination of higher effluent contributions to river flow under baseflow conditions and sub-optimal contact time at the water-sediment interface for the river to reach equilibrium conditions with respect to P sorption (Jarvie et al. 2005). If P cycling cannot self-regulate the river system to lower dissolved SRP concentrations below ecologically damaging levels, then it is imperative that future P supply to the river is reduced to prevent eutrophication causing further ecosystem degradation.

Furthermore, sections of the River Wensum are currently being restored through actions including weir removal, reinitiating meanders, creating pool and riffle sequences, and gravel jetting to remove excess sediment, with the ultimate aim of regaining the original features of a lowland chalk river (e.g. faster heterogeneous flow with exposed gravel beds). If successful, these alterations will see a change from sediment-dominated to gravel-dominated riverbeds across the catchment. As gravel has reduced affinity to adsorb P than silt, the ability of bed sediments to buffer the river against high-SRP concentrations released by STWs could be diminished and P concentrations in the overlying water may increase (Stone et al. 1995; Jarvie et al. 2005). This means that changes to restore the original ecosystem of the River Wensum could in fact lead to more widespread eutrophication problems in the short-medium term unless accompanied by mitigation measures to reduce P loading that would enable the river to reach a new equilibrium at lower P levels.

4 Conclusions

Through the analysis of water and riverbed sediment samples collected across the River Wensum catchment, this research demonstrates that STWs significantly impact upon the P dynamics of this lowland calcareous river system. The mean dissolved SRP concentration recorded downstream of STWs was double that observed at sites unaffected by effluent discharges, whilst the mean EPC0 value of bed sediments was also 70% higher. Irrespective of STW influence and spatial and temporal variability in sample collection, riverbed sediments were almost always acting as sinks for SRP from the overlying water column, with SRP concentrations consistently greater than the sediment EPC0 values. This was particularly true downstream of STWs where sediments displayed greater potential to buffer the river against increases in dissolved SRP released in sewage effluent. With riverbed sediments consistently acting as SRP sinks, these results indicate significant P inputs from other sources, most likely diffuse fertiliser pollution from agriculture, particularly during higher flow flushing events. Despite riverbed sediments acting as sinks, elevated SRP concentrations downstream of STWs evidently demonstrate the sediments have insufficient P sorption capacity to fully buffer the river against effluent discharges. As such, SRP concentrations across the catchment continue to exceed recommended standard for good chemical status, emphasising the environmental requirement for advanced P-stripping technologies to be installed at smaller STWs across the catchment to minimise P loading.