Integrated climate-chemical indicators of diffuse pollution from land to water

Management of agricultural diffuse pollution to water remains a challenge and is influenced by the complex interactions of rainfall-runoff pathways, soil and nutrient management, agricultural landscape heterogeneity and biogeochemical cycling in receiving water bodies. Amplified cycles of weather can also influence nutrient loss to water although they are less considered in policy reviews. Here, we present the development of climate-chemical indicators of diffuse pollution in highly monitored catchments in Western Europe. Specifically, we investigated the influences and relationships between weather processes amplified by the North Atlantic Oscillation during a sharp upward trend (2010–2016) and the patterns of diffuse nitrate and phosphorus pollution in rivers. On an annual scale, we found correlations between local catchment-scale nutrient concentrations in rivers and the influence of larger, oceanic-scale climate patterns defined by the intensity of the North Atlantic Oscillation. These influences were catchment-specific showing positive, negative or no correlation according to a typology. Upward trends in these decadal oscillations may override positive benefits of local management in some years or indicate greater benefits in other years. Developing integrated climate-chemical indicators into catchment monitoring indicators will provide a new and important contribution to water quality management objectives.

Here, as part of this challenge, we examine the possible links between inter-annual weather cycles and water quality 2 considering phosphorus (P) as dissolved reactive phosphorus (DRP) or total reactive phosphorus (TRP) and nitrogen (N) as nitrate in runoff from agricultural land, with regard to their importance for drinking water quality and risk of eutrophication 15 . We specifically investigate the interactions between the local catchment scale and larger, oceanic scale climate and weather patterns.
We acknowledge that hydrometeorological processes drive P and N diffuse pollution and are linked to mobilisation processes which, at the landscape scale, are related to land use and soil nutrient management and status 16 . Such complexity ordinarily prevents the decoupling of agricultural effects from climate effects 14 . We limit our study to agricultural catchments in Western Europe for two reasons. Firstly, European Union (EU) and European Economic Area (EEA) countries are subject to diffuse pollution management and reviews under the EU Water Framework Directive (WFD -OJEU, 2000), associated daughter directives and national legislations. This provides a common policy framework of options of measures and targets. Secondly, we integrate results from three western European seaboard countries that are subject to larger scale oceanic climate and weather patterns and have high temporal resolution water quality monitoring records of diffuse pollution data in agricultural catchments defined as long term research observatories (Table 1). Such high temporal resolution monitoring is required to capture the full dynamics of nutrient loss to water over the year, without being skewed to representative sampling events and periods (baseflow/event flow, summer/winter etc.) 17 and to cover the full range of agri-environmental conditions.
The overall aim of the study was to investigate the influences and relationships between weather processes amplified by the North Atlantic Oscillation (NAO) during a sharp upward trend (2010-2016) and the patterns of diffuse P and N losses to waters across highly monitored catchments in Ireland, Norway and Northwestern France. The objectives were to 1) identify a catchment typology based on characteristics prone to NAO amplifications of diffuse pollution, 2) identify the relationship between NAOi and nutrient concentration trends in streams draining agricultural land, and 3) highlight policy issues for recovery and source management, taking into account these climatic and weather effects.
While anthropogenic warming has influenced weather in Western Europe, for example as observed by an increased number of days with extreme winter rainfall in southern England 18 , weather patterns and trends in this region are largely influenced by large-scale climate systems over the North Atlantic 19 . The NAO is one of the most prominent systems, caused by differentials in the Icelandic low and Azores high pressure systems, and its intensity can be expressed with an index (NAOi) 20 . In Northwestern Europe, a positive phase in the NAOi is often associated with elevated air temperatures in summer and more pronounced rainfall with more frequent large rain events in winter than normal 20,21 . This variation between summer and winter can influence soil drying and wetting periods, and consequently, influence soil chemical lability at the same time as increasing the magnitude and frequency of runoff events. Multi-decadal variability in ocean circulation can be twice as large as the background anthropogenic warming 22 and the effect of weather cycles change may be regionally amplified 23 . Additionally the response of summer temperatures to anthropogenic warming in Europe is likely influenced by the same processes that cause variations in the position of the summer storm track 24 , i.e. large-scale climate systems over the North Atlantic.

Methods
Weather analysis. Regional  NAO intensified and the time when Europe has developed datasets for the evaluation of current measures for mitigating N and P loss to water. The influence of a rise in NAOi on weather was analysed by correlating annual average temperatures, annual total amounts of rainfall and the annual number of large rain events (exceeding 10 and 25 mm).
Hydrochemical analysis. Nutrient loads from catchments are highly influenced by and positively correlated to rainfall-to-runoff processes 25 . Nutrient concentrations in rivers are more subtly influenced by an integration of weather influences, including rainfall-to-runoff processes, rain (or wetness) periodicity 26 and temperature controlled microbial soil nutrient cycling 27 . We therefore compared five to seven hydrological years of annually averaged reactive P and nitrate-N concentrations (hydrological years, 1 st Oct -30 th Sep) in stream water (aggregated from sub-hourly monitoring and bi-weekly samples, n = 26 to 35,040) to the average NAOi for the same period and time step, provided by the NOAA Climate Prediction Centre (www.noaa.gov), for six catchments in Ireland, five in Norway and two in France ( Table 1). The number of analysed years reflects those with available data in each catchment during the period when a sharp rise in NAOi occurred (2010-2016). A simple linear regression analysis of reactive P (RP) and nitrate-N concentrations to NAOi and levels of significance was conducted in SigmaPlot 11.0 ( Table 2). A Cook's distance analysis was conducted to identify potential outliers and influencers. If outliers could be explained by an influential occurrence that was not related to the changed weather conditions, they were omitted from the regression analysis. The NAOi data were in many cases filtered by antecedent simple moving averages (SMA), equally weighting 3 to 5 antecedent years with the current year. This was to avoid temporal delay, such as the yearly differences in time lags of nutrient mobilisation and transfer to stream water in monitored rivers, as well as other influences such as the Gulfstream, which may hide the correlation between RP and nitrate N concentrations with NAOi.
Hydrochemistry and weather data. In the Irish catchments, stream water level was recorded on a 10-min basis using OTT Orpheus Mini vented-pressure instruments installed in stilling wells in all catchments outlets. River discharge was calculated via rating curves developed (in WISKI-SKED software) on Corbett flat-v non-standard weirs at the catchment outlets using the velocity-area method with OTT Acoustic Doppler Current meters. Unfiltered TRP concentrations were monitored concurrently in all catchments outlets using a Hach-Lange Sigmatax-Phosphax suite of instruments. It was assumed that TRP was approximately equivalent to DRP since the flow-weighted mean DRP was previously reported to account for 98-99% of the flow weighted mean TRP in these catchments 28 . The bankside P analysers 29,30 continuously measure TRP (3-4 measurements per hour) by colorimetry using the molybdate-antimony method. The measuring range for TRP is 0.010 mg l −1 to 5.000 mg l −1 and the detection limit is 0.010 mg l −1 . Total oxidised N (TON) was monitored in the outlet on a sub-hourly basis using Hach-Lange Nitratax SC-Plus UV instruments (0.1 -50 mg l −1 ). It was assumed that TON was equivalent to nitrate-N based on low nitrite-N concentrations 31 . Weather data was provided by Met Éireann, the Irish Meteorological Service (www.met.ie/). In the Norwegian catchments hourly average, maximum and minimum discharge were calculated based on water levels, recorded at the catchment outlets every 10 s using a pressure transducer connected to a Campbell data logger, and site-specific rating curves. Composite water samples were collected automatically in the outlets n/a n/a n/a n/a n/a n/a n/a n/a on a volume proportional basis in all catchments 32 . The composite water sample in the container represented the average concentration during the sampling period. The composite water samples were collected for analysis every 14 days, but more frequently during periods with high runoff conditions. Heating-lamps or -cables together with more frequent maintenance was used in winter to prevent ice formation in the monitoring station, thereby guaranteeing year round reliable discharge measurements. Filtered samples (0.45 µm) were used to determine DRP spectrophotometrically by the ammonium molybdate method 33 , with ascorbic acid as a reducing agent and detection limit 0.002 mg l −1 . Nitrate-N was determined according to standard methods (EPA-600/4-79-020) with a detection limit of 0.002 mg l −1 . Results below the detection limit were estimated to be 50% of detection limit. Weather data was provided by the Norwegian Meteorological Institute (www.eklima.no).
In the French Kervidy-Naizin catchment, the stream water level was recorded at the outlets of the catchments every minute using a Thalimèdes OTT float-operating sensor while in the Moulinet catchment a Starflow Ultrasonic Doppler Instrument was used to record stream water level and velocity every 10 minutes. River discharge was calculated via rating curve developed for Kervidy-Naizin and for Moulinet by multiplying the current velocity by the wetted area. To obtain consistent 10 min frequencies, the 1 min measurements in Kervidy-Naizin were subsampled every 10 min. Two monitoring strategies complemented each other for P concentrations acquisition in order to sample both base and storm flows: regular sampling every 1 to 6 days (manual in Kervidy-Naizin; automatic in Moulinet) at approximately the same time (17:00 local time), and high-frequency sampling during storm events with autosamplers (ISCO 6712 Full-Size Portable Sampler). The autosamplers were placed in the shade at the outlet and collected samples when the stream water level reached a threshold, at a frequency of one sample every 30 min for 12 h. Samples were filtered (0.45 µm cellulose acetate filter) directly on-site and sent to the laboratory within 1 week for determination of DRP colorimetrically by reaction with ammonium molybdate with a precision of ±0.004 mg l −1 . In Kervidy Naizin another aliquot of each stream water sample was filtered directly on site to 0.22 µm for determination of nitrate concentration by ionic chromatography (DIONEX DX 100) with a precision of ±2.5%. Weather data was provided by Meteo France (www.meteofrance.com/accueil).
Despite differences in monitoring strategies, all three were designed to provide highly accurate annual P and N concentrations based on sample collection across all flow ranges.

Results and Discussion
Times of change. During the last approximate decade, there has been a distinct upward trend in the NAOi ( Fig. 1) and this has occurred with other records of regional meteorological change in the study areas. This may have offset both hydrological and chemical baselines in Western European Rivers. In Ireland, particularly influenced by the Atlantic Ocean 22 , the 1981-2010 mean annual precipitation over the country was 5% higher and the mean annual air temperature 0.5 degree higher compared with the 1961-1990 average 34 . There were, however, regional variations in precipitation and evidence of an increase in the number of days with heavy rain (10 mm or more) 34 . In Norway, annual average precipitation and air temperatures have increased since 1900 35 . While precipitation increased more in spring and less in summer, temperature increased more in spring and less in winter 35 . In France, the recent period of positive NAOi concurred with an increase in winter rainfall, decrease in spring and autumn rainfall and increase in annual temperature 36 . These national trends were also observed in the western part of France with a more pronounced seasonal distribution 36 .
During the monitoring period, coinciding with the sharp rise of NAOi (2010-2016), there was a significant and positive correlation of annual mean air temperature to the NAOi (Fig. 2) in three national weather stations, one in each country (close to the catchments analysed in this study), with the largest change in the Norwegian station. The annual total rainfall was only significantly positively correlated to the NAOi in the Norwegian station while the number of large rain events correlated positively to NAOi in both the Irish and Norwegian stations (Fig. 2). In the Irish station the number of rain days exceeding 25 mm correlated to NAOi and in the Norwegian station the number of rain days exceeding 10 mm correlated to NAOi, while in the French station there was no correlation of number of rain days exceeding either 10 or 25 mm.
Warmer summers and more frequent large rain events in winter affect soil P and N lability and subsequent mobilisation and transfer processes. Average total P flux and concentrations normally increase with higher rainfall volume and intensity 37 and variations in weather have previously explained observed changes in both stream N fluxes and concentrations 14 . Such changes can impact negatively on aquatic ecology by causing enhanced algal The NAO and nutrient delivery. Warmer summers and increase in heavy winter rains, responding to the sharp rise in NAOi, are divided into two main theoretical categories to explain potential trends of P and N concentrations in rivers.
1. Warmer summers and more rain or more frequent winter storms, as potentially expressed by a positive and larger NAOi, can increase P concentration in streams of headwater catchments if: (i) P is mobilised (erosion 41 and soil biochemistry 42 ) and (ii) surface flow pathways are dominating (conditions controlling flow) 5 , or (iii) below ground transfer dominate but are not attenuated (subsurface flow controls) 43 , or (iv) soil P sources are directly connected to the stream (landscape and management). If P is retained or P sources are not hydrologically connected, larger NAOi will likely decrease P concentration due to dilution or delay the process. That would also be a likely scenario if the main source of P is from a persistent point source 8 . Nitrate is highly mobile and its concentration may increase by an increased below-ground transfer of nitrate-rich groundwater volumes towards the river, and by shorter residence times induced by more runoff (more shallow groundwater tables, higher hydraulic gradients and quicker pathways activated) and therefore less denitrification along the whole transfer continuum, or decrease by: (i) dilution of surface flow paths (if there is high N storage in the groundwater or if there are long residence times) and/or (ii) denitrification due to more saturation. 2. An enhanced summer mineralisation can increase P and N lability at the start of the wet period 44 . In this case P and N concentration will increase with higher surface runoff and leaching with the first rains after summer 45,46 (but not if there is high storage potential in the groundwater or long residence times causing attenuation) 47 .
Integrating this theoretical framework with high temporal resolution nutrient and discharge data composited into annual averages (hydrological years that fully incorporate storm event data), both N and P (as nitrate-N and Reactive Orthophosphate) concentrations were in many cases (P: 10/13; N: 8/12) correlated to NAOi (Fig. 3a- Table 2). The degree of correlation was viewed in terms of the two categories described above (mobilisation, transfer pathway, storage, attenuation, hydrological connectivity and dilution). The influence of changes in fertilizer application to the fields of the catchments was seen to be of minor importance for P and N concentrations during the analysed period and further dampened by the release of legacy P and N 48 . In some catchments, P and N concentrations were positively correlated to the NAOi (P: R 2 = 0.49-0.99 (n = 6) and N: R 2 = 0.28-0.87 (n = 6)) and in other catchments there were negative correlations (P: R 2 = 0.64-0.88 (n = 4) and N: R 2 = 0.45-0.91 (n = 2)), poor or no correlation (P: n = 4 and N: n = 6). Since the catchments response may present  potential time-lags of losses to stream water, NAOi was in those catchments expressed as an antecedent three to five-year moving averages to achieve a stronger correlation between annual NAOi and N and P concentration. Reasons for this time-lag could involve the dominating nutrient mobilisation processes, transfer pathways and catchment size. With this method, we highlight catchment 'memory' effects. This was particularly the case for the Norwegian catchments as explained by the strong influence of the Gulf stream 49 at which latitude correlates to the NAOi with a two-year delay 21 . The results first present RP then nitrate-N, and for each, considering positive and then negative correlation.
Positive correlation, RP. Two catchments, Ballycanew and Mørdre, both have poorly drained soils and P was likely mobilised through erosion and lost to surface water via surface runoff. Reactive P may be desorbed from the eroded soil particles in the stream bed, and in these catchments the annual average RP concentrations were positively and strongly correlated to NAOi (Fig. 3a: Ballycanew, R 2 = 0.97 (excluding an outlying year with highly elevated RP concentration in baseflow conditions, possibly due to point-based sources) and Mørdre (R 2 = 0.99) (excluding an outlying year with intensive rain after sowing and fertilizer application, resulting in the highest registered mean RP concentration during the monitoring period)). Both of these catchments have flashy hydrology (high ratio of storm flow to base flow magnitudes), steep slope gradients (particularly Mørdre) and highly erodible soils 50,51 resulting in a large proportion of RP desorbed from mobilised soil particles resulting from rainfall to runoff partitioning influences 5,50 and a positive correlation of runoff to NAOi. The P concentrations in the rivers of these catchments were therefore highly and positively related to the NAO influence.
The four catchments Timoleague, Creggduff, Kervidy-Naizin and Moulinet are dominated by well-drained soils and mostly groundwater-fed flow. The RP concentrations in these catchments appeared sensitive to weather shifts ( Fig. 3a and Table 2) due to leaching and transfer via groundwater 43 . It was likely that soil chemistry and enhanced mineralisation influenced the mobilisation of P in these catchments, and therefore also the correlation of RP concentration and NAOi. The Timoleague, Kervidy-Naizin and Moulinet catchments have soils with high soil P status, and are iron and carbon rich (with subsequent P sorption competition), favouring P mobilisation into soluble form and leaching to shallow groundwater 42,43,52 . Alternatively, RP is directly exported from highly saturated areas, consequently responding by an increased RP concentration with more rain as hydrological gradients are elevated (R 2 = 0.82, 0.74 and 0.49 respectively). However, since the increase in NAOi had influenced the mean air temperature but not the annual rainfall nor the frequency of high rainfall days in the French catchments ( Fig. 2) a strong positive correlation of RP to NAOi was not expected there. There were only weak positive correlations of runoff to NAOi in the two Irish catchments (data not shown). In the Cregduff catchment, a karst spring contribution zone of pure limestone, P was largely retained in the aluminium rich soil and calcium rich weathered bedrock 53 . However, while a large proportion of P was retained in the karst, an excess of P was transferred to the conduits via slow small fissure and matrix flow where it was settled, remobilised and transferred further by larger event flows 53 . In this catchment P loss was to some degree defined by rainfall to runoff partitioning influences with a positive and strong correlation of P concentration to the antecedent five-year moving average NAOi (R 2 = 0.79). However, the mobilisation and transfer time from source to delivery point may be long and more reflected by an antecedent average NAOi rather than annual fluxes, and the P concentrations in the rivers of these four catchments are more reflective of previously mobilised, attenuated and subsequently remobilised P -i.e. a 'memory' of the NAO influence.
Negative correlation, RP. Despite a positive correlation of runoff to NAOi (data not shown) due to an increased rainfall (Fig. 2), there was a significant negative correlation of RP concentration and NAOi in the Corduff, Skuterud and Skas-Heigre catchments (R 2 = 0.88, 0.64 and 0.83, respectively), all on poor to moderately drained soils ( Fig. 3b and Table 2). The reasons for this relationship were not entirely clear and were likely to be catchment-specific. For example, in the Corduff catchment the negative correlations may be due to combinations of iron and calcium rich soils and high soil P retention in conjunction with hydro-chemical disconnection caused by natural landscape features 54 that strengthened dilution. In Skuterud, the livestock production ceased in 2010, which may influence point source RP concentrations. Skas-Heigre are dominated by high livestock density and high soil P status, but were influenced by an official P-mitigation project from 2010. Additionally, the Skuterud and Skas-Heigre, are both equipped with sedimentation ponds to mitigate sediment and particulate P transfer and which may also, in part, act as a mode to disconnect RP transfers by slowing down the flow, allowing for retention and cause dilution. However, the hydrological conditions that increase RP transfer during extreme positive phases of NAOi may also act to overcome catchment dis-connectivity and such catchments will require further research with focus on specific periods. An alternative explanation could be that these two catchments have small but persistent point sources likely to be diluted by more rain, more frequent large rain events (Fig. 2) and therefore more runoff.
No correlation, RP. The Castledockerell, Dunleer, Time and Vasshaglona catchments, all on well-drained soils, had no correlation of RP concentration to NAOi (R 2 = 0.01, 0.19, 0.02 and 0.20 respectively, Fig. 3c). In both Castledockerell and Time this can be explained by high soil P retention caused by soils rich in Al with less P solubility 43,55 . Positive correlation, nitrate-N. An increased number of days where soils have high soil temperatures and optimal soil moisture may enhance stored (short and long-term) organic N mineralisation and nitrification 45,46 , and if this is followed by more frequent large rain events, nitrate is easily leached to groundwater and loss to the stream water is likely to occur. In such cases, if the stream is groundwater fed and the pathways are not attenuated, The three catchments, Castledockerell, Timoleague and Time, are all dominated by well-drained soils and with mostly groundwater-fed rivers. Under such conditions excessive nitrate could be leached to groundwater and transferred to the rivers. The annual average nitrate-N concentrations were positively correlated to NAOi in these catchments (Fig. 4a: R 2 = 0.87, 0.57 and 0.76 respectively). The strongest correlation was found in the Castledockerell catchment, which largely lacks denitrifying conditions and has relatively short time lags due to a thick (up to ca. 20 m), highly permeable layer of weathered slate bedrock 56 and a positive correlation of rainfall and runoff to NAOi. The Timoleague catchment, with its permeable sandstone geology and denitrifying conditions in the near stream zones 56,57 had a lower correlation of nitrate-N concentration to NAOi, even if the source pressures were stronger than in Castledockerell. The Time catchment possibly has a large proportion of shallow groundwater together with a network of tile drains transferring nitrate to the river that was largely influenced by the NAOi.
The nitrate-N concentration in Ballycanew catchment, with a flashy hydrology, had a weak positive correlation to NAOi (R 2 = 0.28). This catchment has a transfer via shallow perched groundwater to the rivers. Moreover, nitrate-N is transferred to the river via a line of springs, where the upland well drained soils meet the low land poorly drained soils. Another important transfer pathway in that catchment is an extensive network of tile drains and ditches 58 that could transfer nitrate-N relatively quickly with little attenuation, thus responding positively to changes in NAOi.
Negative correlation, nitrate-N. In two of the catchments, Skas-Heigre and Kervidy-Naizin, the nitrate-N concentration responded negatively and correlated well to the increased NAOi ( Fig. 4b: R 2 = 0.45 and R 2 0.91 respectively). In Skas-Heigre, with a positive correlation of rainfall and runoff to NAOi, one explanation for the negative correlation is the potential mitigating effect of the sedimentation pond, which slows down the stream flow allowing for denitrification and dilution during rainfall. Furthermore, a high groundwater table in this catchment may contribute to high denitrification rates. However, this will need further testing. The Kervidy-Naizin catchment has a large N storage capacity in the groundwater which may mask the effect of leaching and during subsequent winter rain events there will be higher water flux but lower nitrate-N concentrations 57 . This catchment has been found to lack a relationship between annual N surplus and stream nitrate-N concentration due to higher temperatures and high nitrate-N fluxes from groundwater to stream, and long residence times 59 . No correlation, nitrate-N. The Dunleer, Cregduff, Corduff, Vasshaglona, Skuterud and Mørdre catchments had no correlation of nitrate-N concentration to NAOi (R 2 from 0.00 to 0.23, Fig. 3c). In the Dunleer catchment perched shallow groundwater is the highest source for N transfer to the stream, and this is only hydrologically connected, via open ditches that are active in winter, when the groundwater table is high 60 . It was noted that seasonal N concentration amplitudes were greatest in Dunleer during wet autumns following dry summers (i.e. increased lability -not shown). However, the overall effect appears to be dilution during periods of increased rainfall frequency in this catchment and so the ditches acted as limited and seasonal 'on-off switches' , hence weakening the correlation to NAOi. In the karst Cregduff catchment nitrate-N was likely leached and transferred via slow micro fissures to the conduits, where up to 90% of the P transfer occurred during spring flow events 53 , and with a large potential for denitrification and weakening of the influence of NAO. The Corduff and Vasshaglona catchments have shallow groundwater tables with a potential for denitrification in the near stream zones. The Skuterud and Mørdre catchments are both intensively tile drained and runoff is transferred quickly on the soil surface and via macropores to tile drains with little contact to the soil matrix and hence nitrate-N concentrations may be diluted. Furthermore, Skuterud catchment has a potential for denitrification in the sediment pond and this may also weaken the correlation to NAOi. This has not been shown.

Systems science integration.
Such relationships, only found by the analysis of data covering the full hydrochemical discharge range, have implications for the interpretations of data in catchments showing, for example, decreases in nutrient source pressures or implementations of other catchment mitigation measures 61 . An increase in NAOi, increasing the average air temperatures and causing more rain or more days with large rain events, may override these positive benefits in some years or indicate greater benefits in other years -and this is likely to be catchment-specific. A consideration of dominating mobilisation and transfer processes, storage and hydrological connectivity should, therefore, be included as moderating factors in policy reviews of diffuse pollution management.
Using the oceanic scale NAOi and local scale water quality data is an example of a 'climate-chemical indicator' in the framework of integrating systems-science 1 . However, planning for diffuse pollution mitigation, when some years will indicate smaller or greater benefits against a baseline of measures, may still leave water quality and agricultural actors and stakeholders with a conundrum. That is: how to integrate these chemical-climate systems further with socio-political indicators of environmental improvement against a back drop of expectations and targets. Furthermore, another integration is the need to sustain agricultural production for economic and food security. Over-engineering mitigations by increasing catchment resilience in vulnerable catchments to account for positive phases of NAOi (or other weather amplifications) or accepting phases of failed or exceeded targets will leave certain elements of the overall integrated system wanting. Developing a deeper understanding of the economic consequences of, for example, over-engineering or optimizing all elements of the systems related to catchment-specific diffuse pollution may be a solution that is more readily understood and implemented 62,63 . At the very least, developing climate-chemical indicators into catchment vulnerability indicators will provide a new and important component to water quality management objectives in agricultural catchments.

Conclusions
Oceanic-scale climate patterns, characterized by the intensity of the North Atlantic Oscillation (NAOi), was found to influence nutrient loss to water to a varying degree in agricultural river catchments in the northwest of Europe (Ireland, Norway and North Western France). A sharp upward trend in NAOi, which has amplified annual weather patterns in recent time, may override or underride benefits of local management, in a specific way in some catchments and years. A new and important component to diffuse pollution and water quality management objectives is therefore to develop integrated climate-chemical indicators into catchment response indicators. While simple mean residence time is often mentioned, it is also needed to consider the hydro-biogeochemical sensitivity of the catchments to climate variations.