<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Span. J. Soil Sci.</journal-id>
<journal-title>Spanish Journal of Soil Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Span. J. Soil Sci.</abbrev-journal-title>
<issn pub-type="epub">2253-6574</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">11108</article-id>
<article-id pub-id-type="doi">10.3389/sjss.2023.11108</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Science archive</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Nitrate and Phosphorus Transport in a Galician River (NW Iberian Peninsula): Insights From Fourteen Years of Monitoring</article-title>
<alt-title alt-title-type="left-running-head">Rodr&#xed;guez-Blanco et al.</alt-title>
<alt-title alt-title-type="right-running-head">Nitrogen and Phosphorus Transport</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Rodr&#xed;guez-Blanco</surname>
<given-names>M. L.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2108864/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Taboada-Castro</surname>
<given-names>M. M.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Taboada-Castro</surname>
<given-names>M. T.</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2129220/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Physical Geography Area, History, Art and Geography Department</institution>, <institution>GEAAT Group</institution>, <institution>University of Vigo</institution>, <addr-line>Ourense</addr-line>, <country>Spain</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Area of Soil Science and Soil Chemistry</institution>, <institution>Higher Technical School of Agricultural Engineering</institution>, <institution>University of Valladolid</institution>, <addr-line>Palencia</addr-line>, <country>Spain</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Faculty of Sciences</institution>, <institution>University of A Coruna</institution>, <addr-line>A Coru&#xf1;a</addr-line>, <country>Spain</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/89410/overview">Avelino N&#xfa;&#xf1;ez-Delgado</ext-link>, University of Santiago de Compostela, Spain</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: M. L. Rodr&#xed;guez-Blanco, <email>maria.luz.rodriguez.blanco@uvigo.es</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>13</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>13</volume>
<elocation-id>11108</elocation-id>
<history>
<date date-type="received">
<day>06</day>
<month>12</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>05</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Rodr&#xed;guez-Blanco, Taboada-Castro and Taboada-Castro.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Rodr&#xed;guez-Blanco, Taboada-Castro and Taboada-Castro</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The long-term development of water quality metrics is critical to estimate the status of water resources and successful catchment management. This study looks at the temporal evolution of nitrate and phosphorus loads and concentrations in a stream draining rural catchment in Galicia (NW Iberian Peninsula) during 14 hydrological years. The concentrations of both nutrients are relatively low because of the rural nature of the research catchment and low-intensity agriculture. However, catchment nutrient reference levels were widely exceeded. The hydrology controls nutrient fluxes as flow is the main driver of nutrient transfer. For both nutrients there were no consistent trends in nutrient fluxes. The analysed nutrients, however, have been seen to behave differently. As a result, N seems to be decreasing, mainly due to decreased transit in the autumn and spring. Phosphorus, which is primarily transported as particulate matter, exhibits an increase in winter linked to a greater flow and sediment transfer.</p>
</abstract>
<kwd-group>
<kwd>long-term</kwd>
<kwd>nutrients</kwd>
<kwd>Galicia</kwd>
<kwd>rural catchment</kwd>
<kwd>land use</kwd>
</kwd-group>
<contract-num rid="cn002">PGIDIT05RAG10303PR</contract-num>
<contract-sponsor id="cn001">Ministerio de Educaci&#xf3;n y Ciencias<named-content content-type="fundref-id">10.13039/100019954</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Xunta de Galicia<named-content content-type="fundref-id">10.13039/501100010801</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>The natural quality of inland waters is related to the geomorphology, climate, and land use in the catchment. Thus, size and slope of the catchment, rainfall, temperatures, vegetation, erosion and soil structure play an important role in water quality. However, due to human activities, the high loads of sediments and nutrients transported by rivers have become one of the greatest problems causing the deterioration of the quality of surface waters on a global scale (<xref ref-type="bibr" rid="B4">Bennet et al., 2001</xref>; <xref ref-type="bibr" rid="B49">Worrall et al., 2016</xref>; <xref ref-type="bibr" rid="B22">Grizzetti et al., 2021</xref>; <xref ref-type="bibr" rid="B27">Kelly et al., 2022</xref>), and this explains why in recent years, determined efforts have been made both in the field of research and management and restoration.</p>
<p>Although water management and use has been one of the main concerns of the European Water Policy since the 1960s, it was only in the year 2000 that progress was made in the protection of aquatic resources with the publication of the Water Framework Directive (WFD; <xref ref-type="bibr" rid="B13">Directive, 2000</xref>) by the European Parliament and the Council of the European Union, in which a community action framework was established in the field of water policy. The WFD focuses on the protection of all water bodies present in the territory of the European Community. Its objectives are to prevent the deterioration of water resources, protect and improve their condition, promote the sustainable use of the resource, guarantee a reduction in pollution, prevent the effects of floods and droughts, and achieve at least good ecological status by the year 2015. Despite the progress made with the scope of improving the environmental quality of the various water sources in Europe, the last report by the European Environment Agency (<xref ref-type="bibr" rid="B19">European Environment Agency, 2018</xref>) indicated that 47% of European surface waters continue to fail to meet the minimum objective of good ecological status, with agricultural activities and the urbanisation of rural landscapes being the main threats to fluvial courses (<xref ref-type="bibr" rid="B19">European Environment Agency, 2018</xref>). Therefore, the frequent use of agrochemicals to meet production demands means that agricultural systems play an important role as a source of contamination (non-point source) of water bodies (<xref ref-type="bibr" rid="B39">Oenema et al., 2005</xref>; <xref ref-type="bibr" rid="B47">Sharpley et al., 2008</xref>; <xref ref-type="bibr" rid="B30">Lassaletta et al., 2009</xref>). Moreover, over the last few decades, the progressive implementation of modern agriculture and livestock methods, with the introduction of increasingly intensive tillage techniques, an increase in the size of plots, the establishment of monocultures, etc. has promoted runoff and soil erosion (<xref ref-type="bibr" rid="B40">Panagos et al., 2015</xref>). In general, chemicals from non-point sources are transported to nearby water bodies either in dissolved form or along with sediment during heavy rainfall and erosion events, even in gently sloping areas (<xref ref-type="bibr" rid="B47">Sharpley et al., 2008</xref>; <xref ref-type="bibr" rid="B16">Dupas et al., 2017</xref>).</p>
<p>The latest reports from the European Environment Agency noted an improvement in 20% of the stations, directly attributed to the implementation of measures to control diffuse (non-point) contamination and increased efficiency in wastewater treatment (<xref ref-type="bibr" rid="B10">Carvahlo et al., 2019</xref>). Thus, based on reported long-term data for nitrates in European waters, there was a falling trend for nitrate concentration in rivers, although this has levelled off since 2009. The reason for nitrate persistence is due to overuse of nitrogen, which can remain in the subsurface flow for decades in catchments with long time response and add to increased nitrate loads long after application practices have changed, or mitigation strategies have been implemented (<xref ref-type="bibr" rid="B20">Fovet et al., 2015</xref>; <xref ref-type="bibr" rid="B17">Ehrbardt et al., 2021</xref>). These and other provenances of nitrate create a heterogeneous patchwork of source areas across the landscape that can be activated and deactivated in response to changing hydrological conditions (<xref ref-type="bibr" rid="B15">Dupas et al., 2019</xref>). For phosphorus, a marked decline was reported due to improved wastewater treatment, reduction of P in detergents, as well as reduction in soil fertilization rates and implementation of measures to reduce soil erosion and connectivity to the receiving waters (Oenema eta l., 2005; <xref ref-type="bibr" rid="B34">Minaudo et al., 2015</xref>; <xref ref-type="bibr" rid="B35">Moatar et al., 2017</xref>). However, phosphorus may persist through soil enrichment, and therefore it may take a long time before any major water quality improvements can be detected in the recipient water bodies after decreasing fertilisation of agricultural fields (<xref ref-type="bibr" rid="B48">Withers and Lord, 2002</xref>; <xref ref-type="bibr" rid="B18">Ekholm et al., 2015</xref>).</p>
<p>Correct assessment of the quality of water resources in the catchment is essential to estimate the status of the resources and start management and planning, if required. In fact, the WFD urges members to carry out studies on trends of different contaminants (with trend understood to be the monotonic variation of the concentration in relation to time) and identify those water bodies exceeding established reference nutrient levels. This requires sufficiently frequent monitoring of the water quality over long periods of time, which is an arduous, costly, and time-consuming task. At present, several countries, including Spain, must improve their control programs to guarantee broad and consistent monitoring of water bodies, with appropriate coverage of all relevant quality indicators, since there are still significant deficiencies. However, the number of checkpoints has decreased in recent years as compared to the first planning cycles of the WFD. This includes the Galician region, popularly known as &#x201c;the country of a thousand rivers.&#x201d; Despite this, the scientific community did not pay much attention to the health of Galician rivers. To date, they have not undergone extensive monitoring because hydrological and hydrochemical studies based on catchment monitoring are relatively recent in the community. This underlines the need to continue long-term monitoring efforts for water quality assessments and to identify trends.</p>
<p>The objective of this paper was to explore the temporal evolution of nitrate and phosphorus loads and concentrations in a stream draining a typical rural catchment in Galicia (NW Iberian Peninsula), i.e., small population centres and low-intensity agriculture. The study is based on information collected over a period of 14&#xa0;years in which regular manual sampling was combined with automatic sampling during runoff episodes.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>Material and Methods</title>
<sec id="s2-1">
<title>Study Area</title>
<p>The study site is the Corbeira stream headwater catchment located in NW Spain (<xref ref-type="fig" rid="F1">Figure 1</xref>, gauge at 43&#xb0; 12&#x2032;58&#x2033;N, 8&#xb0; 13&#x2032;48&#x2033;W). The Corbeira is a right-bank tributary of the Mero River, the main river contributing to the Cecebre reservoir (water supply for the city of A Coru&#xf1;a and surroundings (about 450,000 inhabitants) discharging into the Ria do Burgo, in the vicinity of the city of A Coru&#xf1;a (<xref ref-type="fig" rid="F1">Figure 1</xref>). The total length of the stream is 10&#xa0;km, and the whole catchment area covers about 16&#xa0;km<sup>2</sup>, with altitudes ranging from 65&#xa0;m a.s.l. at the confluence with the Mero River to 470&#xa0;m a.s.l. at the highest part. The mean slope is 19%, although approximately 27% of the catchment is steep, i.e., slopes higher than 25%. From a geological point of view, it is located within the Ordenes Complex and is dominated by basic schists (<xref ref-type="bibr" rid="B24">IGME, 1981</xref>), with phyllites and biotitic schists the most common petrological types in the area, which confers a certain uniformity to the lithological substrate. Schists are made up of easily alterable minerals, such as biotite (sometimes chlorite), plagioclase, and amphiboles. The soils are Umbrisols and Cambisol (<xref ref-type="bibr" rid="B26">IUSS Working Group WRB, 2014</xref>) with a silt and silty loam texture developed. The organic matter content in the upper horizon varies from less than 5.3% in cultivated lands, to 10.5% in forest soils, highlighting the importance of the latter in carbon sequestration. Total N soil concentrations ranged between 0.26% on cultivated land and 0.47% on grassland, whereas forest soils showed a mean value of 0.39% (<xref ref-type="bibr" rid="B45">Rodr&#xed;guez-Blanco et al., 2015</xref>). The P Olsen content varied between 14.0&#xa0;mg&#xa0;kg<sup>&#x2212;1</sup> in the forest soils and 69.0&#xa0;mg&#xa0;kg<sup>&#x2212;1</sup> in cultivated fields, with the mean content of grassland 19.5&#xa0;mg&#xa0;kg<sup>&#x2212;1</sup> (<xref ref-type="bibr" rid="B44">Rodr&#xed;guez-Blanco et al., 2013</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Location of the study area. Changes in land use observed in the study area: <bold>(A)</bold> 2000 and <bold>(B)</bold> 2018 based on SIOSE.</p>
</caption>
<graphic xlink:href="sjss-13-11108-g001.tif"/>
</fig>
<p>The main land use is forest, which covers 65% of the catchment area and is mostly used for commercial plantations of eucalyptus and pine. In addition, there is, as in most Galician River catchments, a dense riparian vegetation basically consisting of <italic>Alnus glutinosa</italic>, <italic>Betula alba</italic> and <italic>Salix spp</italic>. The agricultural surface represents 30% of the catchment and is dominated by pasture and meadows (26%), the latter being humid areas mainly in flood zones. They are used to produce green, dry, or ensiled forage for consumption. Cropland comprises approximately 4% of the catchment area and is essentially used for forage crops, such as maize and winter cereals. Potatoes, vegetables, and fruit trees are cultivated in small kitchen gardens and orchards. The remaining 5% of the catchment&#x2019;s surface is occupied by impervious areas, mostly road infrastructures (roads and forest paths) distributed throughout the catchment, crossing the main channel and its tributaries at numerous points. This information is based on land use maps from the Spanish Land Cover/Land Use Information System (SIOSE), used to calculate the cover (%) of land use for the years 2004 and 2018. A slight decrease was observed in the agricultural area during this period, whereas forests increased due to the replacement of cultivated land by plantations of fast-growing tree species, mainly eucalyptus. There is also a decrease in stocking density (dominated by cattle) during the study period, going from 0.96 livestock units (LU) per ha of utilised agriculture area (UAA) in 2004 to 0.73 LU per ha of UAA in 2018.</p>
<p>The distribution of land uses in the catchment is characterized by the mosaic presented in agricultural areas (<xref ref-type="fig" rid="F2">Figure 2</xref>), which is generally bordered by meadows and, in some cases, by small stone walls, characteristic of smallholdings in Galicia. This distribution could be classified as a low-intensity agricultural mosaic (cropland and grassland), with a light input of inorganic fertilizer, small field sizes and low livestock density (<xref ref-type="bibr" rid="B14">Dou et al., 2021</xref>).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Illustrative photographs of the study area. <bold>(A)</bold> Landscape mosaic in which meadows, cultivated land and forest area; <bold>(B)</bold> runoff generated on cultivated land devoid of vegetation and channeled through the ditch; <bold>(C)</bold> slurry application on the ground at field capacity, and <bold>(D)</bold> riverside vegetation.</p>
</caption>
<graphic xlink:href="sjss-13-11108-g002.tif"/>
</fig>
<p>The study area has a marked rural nature (&#x3c;150 inhab. km<sup>&#x2212;2</sup>; Organization for Economic Co-operation and Development- <xref ref-type="bibr" rid="B38">OECD, 2007</xref>), characterized by a low population density, which has gone from 35 inhab. km<sup>&#x2212;2</sup> at the beginning of the study (2004) to 32 inhab. km-<sup>2</sup> in 2018 (rates calculated from data of Spanish Statistical Office-<xref ref-type="bibr" rid="B25">INE, 2019</xref>)- a 15% decrease in just 15&#xa0;years. In fact, the Ministry of Agriculture statistics place Galicia at the head of the most isolated communities in rural areas. Between the years 2010&#x2013;2018, the population in Galician villages fell by 17.9% (<xref ref-type="bibr" rid="B25">INE, 2019</xref>). Most of the population is distributed in scattered dwellings that lack a sewage system, so that domestic wastewater is stored in cesspools constituting small sources of contamination.</p>
<p>The climate is temperate oceanic, with an annual average temperature of 13.5&#xb0;C ranging from 8.9&#xb0;C in January and February (winter) to 18.8&#xb0;C in July and August (summer). The rainfall is abundant (1,074&#xa0;mm- period 1983/1984&#x2013;2017/2018) and occurs throughout the year, although it is concentrated in the autumn and winter, with October, November and December being the wettest months and July and August the driest. The hydrological regime is pluvial oceanic; which implies two alternate periods: a high-discharge period during the wet period (i.e., late autumn and winter) and a low-discharge period during the dry season (i.e., late spring and summer) (<xref ref-type="fig" rid="F3">Figure 3</xref>; <xref ref-type="bibr" rid="B42">Rodr&#xed;guez-Blanco et al., 2020</xref>). Mean stream discharge is 0.181&#xa0;m<sup>3</sup>s<sup>&#x2212;1</sup>, ranging from 0.08 to 0.273&#xa0;m<sup>3</sup>s<sup>&#x2212;1</sup> (<xref ref-type="bibr" rid="B42">Rodr&#xed;guez-Blanco et al., 2020</xref>). Baseflow contribution to the total streamflow (i.e., baseflow index), calculated according to <xref ref-type="bibr" rid="B23">Gustard et al. (1992)</xref>, varies between 0.61 and 0.78 with an average of 0.71 (unpublished data), indicating that 71% of long-term streamflow in the study area is originated from groundwater.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Box plot of seasonal distribution of discharge in the Corbeira catchment. The length of the box plot shows the interquartile range, comprising 50% of the values, and the straight horizontal line inside the box represents the median. (Data from <xref ref-type="bibr" rid="B42">Rodr&#xed;guez-Blanco et al., 2020</xref>). A: autumn, W: winter, Sp: spring, S: summer.</p>
</caption>
<graphic xlink:href="sjss-13-11108-g003.tif"/>
</fig>
</sec>
<sec id="s2-2">
<title>Field Instrumentation and Water Sampling</title>
<p>The data used in this study covers 14 hydrological years from October 2004 to September 2018. Monitoring within the catchment included rainfall and discharge measurements and automatic and manual sampling of stream water at the catchment outlet. Rainfall was monitored at three locations across the catchment at 10-min intervals using tipping bucket rain gauges (0.2&#xa0;mm resolution), followed by the Thiessen Polygon method (<xref ref-type="bibr" rid="B31">Linsely et al., 1949</xref>) to calculate the mean rainfall in the area. Stream water level was measured at the catchment outlet using a pressure transducer (ISCO 720), and discharge was calculated from the rating curve. Measured water levels at 1&#xa0;min intervals were aggregated to 10-min values. Water sampling was carried out under baseflow conditions (manual sampling every 2&#xa0;weeks) and during runoff events (automatic sampling, ISCO 6712-FS). Full details of the design of the monitoring and water sampling can be found in <xref ref-type="bibr" rid="B44">Rodr&#xed;guez-Blanco et al. (2013</xref>, <xref ref-type="bibr" rid="B45">2015)</xref>.</p>
</sec>
<sec id="s2-3">
<title>Laboratory Analyses</title>
<p>Water samples were immediately brought back to the laboratory and refrigerated (4&#xb0;C) until analysed for nitrate and total and dissolved phosphorus concentrations. After sample filtration (0.45&#xa0;&#x3bc;m), nitrate concentrations (expressed as mg NO<sub>3</sub>
<sup>&#x2212;</sup> L<sup>&#x2212;1</sup>) were determined by capillary electrophoresis. Total phosphorus (TP) (&#x3bc;g P L<sup>&#x2212;1</sup>) was determined in unfiltered water samples by acid digestion with ammonium persulfate (<xref ref-type="bibr" rid="B3">APHA, 1998</xref>) followed by colorimetric determination using the molybdate-ascorbic acid method developed by (<xref ref-type="bibr" rid="B36">Murphy and Riley, 1962</xref>). Following water sample filtration, dissolved phosphorus (DP) was determined by ICP-MS. Particulate phosphorus (PP) was calculated as the difference between TP and DP.</p>
</sec>
<sec id="s2-4">
<title>Data and Statistical Analysis</title>
<p>Since the main objective of this work was to detect trends of nitrate and phosphorus in long time series, the load (of these nutrients) for different periods of time (monthly, seasonal, or annual) was calculated by summing up the products from the mean concentration between two consecutive samples and cumulative discharge for each time interval between two water samples. The discharge-weighted mean concentration of each nutrient for each specific period (monthly, seasonal, or annual) was calculated as the nutrient load of each period divided by the total discharge of the period. For the seasonal trend, the year was divided into four seasons: autumn (October, November, and December), winter (January, February, and March), spring (April, May, and June) and summer (July, August, and September).</p>
<p>Once the nutrient concentrations were estimated, two different methodologies were used to calculate the trend of nitrate concentrations. The first refers to the approach used by the member states in the context of the Nitrates Directive, which is based on calculating the differences between the average concentration of nitrates taken from two consecutive periods of 4 years. Differences in the range between &#x2212;1 and 1&#xa0;mg&#xa0;L<sup>&#x2212;1</sup> indicates the absence of trend, whereas those higher than 1 and 5&#xa0;mg&#xa0;L<sup>&#x2212;1</sup> are rated as a positive and strongly upward trend. In the second approach, trend analysis was performed using the non-parametric Mann-Kendall test and the Sen method in the Maksesen program (developed by the Finnish Meteorological Office). Both tests are widely used for trend detection in hydro-meteorological records and water quality (<xref ref-type="bibr" rid="B7">Bouraoui and Malag&#xf3;, 2020</xref>) and have also been recommended by the World Meteorological Organisation for analysing the trends in hydrometeorological time series (<xref ref-type="bibr" rid="B29">Kumar et al., 2009</xref>). The presence of a monotonic upward or downward trend was tested with the Mann-Kendall and the magnitude (i.e., slope of a linear trend) was estimated with the Sen method.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec id="s3-1">
<title>Exploratory Analysis: Annual and Seasonal Nitrate and Phosphorus Loads and Concentrations</title>
<p>The average nitrate and phosphorus export during the study period was 19.04&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup> (7.36 and 30.07&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup>) and 0.13&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup> (0.08 and 0.22&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup>) for nitrate and phosphorus respectively (<xref ref-type="fig" rid="F4">Figure 4</xref>). In both cases, the highest values were registered in the hydrological year 2009/2010 and the lowest were measured in the hydrological year 2016/17 (<xref ref-type="fig" rid="F4">Figure 4</xref>), similar to discharge variations in the study period. This highlights the importance of the hydrological response in the nutrient export in the study catchment, because small streams can be very sensitive to changes in the precipitation regime and runoff episodes compared to large basins (<xref ref-type="bibr" rid="B6">Bol et al., 2018</xref>). The results showed that PP is the main fraction of P mobilized from the catchment, representing approximately 72% of the phosphorus exported, although the proportions vary between 39% in the hydrological year 2013/14 and 86% in 2017/18. This variability is probably related to the characteristics of the rainfall-runoff events, management activities within the catchment in relation to meteorological conditions, and the source of sediments delivered to the river channel. The P losses occur primarily through erosion and surface runoff as P is mobilized and transported during runoff events (<xref ref-type="bibr" rid="B45">Rodr&#xed;guez-Blanco et al., 2015</xref>).</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Temporal evolution of annual <bold>(A)</bold> rainfall, <bold>(B)</bold> discharge, <bold>(C)</bold> NO<sub>3</sub>
<sup>&#x2212;</sup> load, P load <bold>(D)</bold>, mean <bold>(E)</bold> NO<sub>3</sub>
<sup>&#x2212;</sup> concentration and <bold>(F)</bold> P mean concentrations during the study period. TP: total phosphorus, PP: particulate phosphorus, DP: dissolved phosphorus.</p>
</caption>
<graphic xlink:href="sjss-13-11108-g004.tif"/>
</fig>
<p>Mean annual concentrations ranged from 4.1 to 5.2&#xa0;mg&#xa0;L<sup>&#x2212;1</sup> for NO<sub>3</sub>
<sup>&#x2212;</sup>, and 17.2 and 35.1&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup> for TP (<xref ref-type="fig" rid="F4">Figure 4</xref>), reaching mean values of 4.8&#xa0;mg&#xa0;L<sup>&#x2212;1</sup> and 26.3&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>, respectively, which signifies a greater variability of phosphorus (CV &#x3d; 20%) concentrations than nitrate (CV &#x3d; 6%). The highest mean concentration of nitrate was recorded in the hydrological year 2006/07, while the lowest values were obtained in 2004/05. For total phosphorus, the highest mean concentrations were recorded in hydrological year 2005/06 while the lowest values, also for dissolved phosphorus, were recorded in 2007/08. Therefore, the lowest concentrations are usually recorded in the driest years and with lower flow rates (<xref ref-type="fig" rid="F4">Figure 4</xref>).</p>
<p>Transport throughout the year (<xref ref-type="fig" rid="F5">Figure 5</xref>) showed that the highest nitrate loads were obtained in winter (8.45&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup>), followed by spring (4.78&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup>) and autumn (4.54&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup>), with the lowest values recorded in summer (1.25&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup>). The pattern for phosphorus was similar, which also returned the highest records in winter (0.050&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup>), followed by autumn (0.030&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup>) and spring (0.031&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup>), with the lowest values observed in summer (0.006&#xa0;kg&#xa0;ha<sup>&#x2212;1</sup>) (<xref ref-type="fig" rid="F5">Figure 5</xref>). The seasonal distribution of the concentrations (<xref ref-type="fig" rid="F6">Figure 6</xref>) followed the same pattern as the loads, meaning that nitrate showed the following decreasing order: winter (5.5&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>) &#x3e; autumn (5.3&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>) &#x3e; spring (4.5&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>) &#x3e; summer (4.1&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>), while phosphorus displayed the sequence: autumn (TP: 32.75&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; PP: 22.91&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; DP: 10&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>) &#x3e; winter (TP: 33&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; PP: 23&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; DP: 10&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>) &#x3e; spring (TP: 24&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; PP: 15&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; DP: 8&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>) &#x3e; summer (TP: 16&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; PP: 9&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; DP: 7&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>). Phosphorus concentrations showed large differences between autumn and winter against summer.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Box plot with the load of <bold>(A)</bold> nitrate (NO<sub>3</sub>
<sup>&#x2212;</sup>), <bold>(B)</bold> total phosphorus (TP), <bold>(C)</bold> particulate phosphorus (PP) and <bold>(D)</bold> dissolved phosphorus (DP). The length of the box plot shows the interquartile range, comprising 50% of the values, and the straight horizontal line inside the box represents the median. A: autumn, W: winter, Sp: spring, S: summer.</p>
</caption>
<graphic xlink:href="sjss-13-11108-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Box plot with the concentration of <bold>(A)</bold> nitrate (NO<sub>3</sub>
<sup>&#x2212;</sup>), <bold>(B)</bold> total phosphorus (TP), <bold>(C)</bold> particulate phosphorus (PP) and <bold>(D)</bold> dissolved phosphorus (DP). The length of the box plot shows the interquartile range, comprising 50% of the values and the straight horizontal line inside the box represents the median. A: autumn, W: winter, Sp: spring, S: summer.</p>
</caption>
<graphic xlink:href="sjss-13-11108-g006.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>Trends in Nitrate and Phosphorus Loads and Concentrations</title>
<p>
<xref ref-type="table" rid="T1">Tables 1</xref>&#x2013;<xref ref-type="table" rid="T3">3</xref> present the results obtained with the Mann-Kendall test. Statistically relevant trends in annual values were not observed for both nutrients (<xref ref-type="table" rid="T1">Table 1</xref>). Nitrate seems to show a limited downward prospect. The same results were obtained using the Nitrate Directive approach, because the difference between the NO<sub>3</sub>
<sup>&#x2212;</sup> flow-weighted concentrations of two consecutive reporting periods ranged between &#x2212;1 and 1&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>, indicating the absence of a trend. On a seasonal scale, a meaningful decrease is only observed in autumn for the nitrate load (<xref ref-type="table" rid="T2">Table 2</xref>), mainly linked to the decrease in October, as well as an increase in TP concentrations in winter (<xref ref-type="table" rid="T3">Table 3</xref>). A negative trend is also seen in the load of NO<sub>3</sub>
<sup>&#x2212;</sup> and P in the month of September, while the concentrations show a positive pattern, which could be related to the lower discharge in this month and the enrichment of the flow due to less dilution.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Mann-Kendall trend test results for nutrient yield and concentrations.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Variable</th>
<th align="center">Mann-Kendall test</th>
<th align="center">Sen&#x2019;s slope</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">NO<sub>3</sub>
<sup>&#x2212;</sup> (kg ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup>)</td>
<td align="char" char=".">&#x2212;1.20</td>
<td align="char" char=".">&#x2212;0.485</td>
</tr>
<tr>
<td align="left">TP (kg ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup>)</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">0.001</td>
</tr>
<tr>
<td align="left">PP (kg ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup>)</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">0.001</td>
</tr>
<tr>
<td align="left">DP (kg ha<sup>&#x2212;1</sup> year<sup>&#x2212;1</sup>)</td>
<td align="char" char=".">&#x2212;0.11</td>
<td align="char" char=".">0.000</td>
</tr>
<tr>
<td align="left">NO<sub>3</sub>
<sup>&#x2212;</sup> (mg L<sup>&#x2212;1</sup>)</td>
<td align="char" char=".">&#x2212;0.99</td>
<td align="char" char=".">&#x2212;0.018</td>
</tr>
<tr>
<td align="left">TP (&#x3bc;g L<sup>&#x2212;1</sup>)</td>
<td align="char" char=".">0.55</td>
<td align="char" char=".">0.277</td>
</tr>
<tr>
<td align="left">PP (&#x3bc;g L<sup>&#x2212;1</sup>)</td>
<td align="char" char=".">0.55</td>
<td align="char" char=".">0.176</td>
</tr>
<tr>
<td align="left">DP (&#x3bc;g L<sup>&#x2212;1</sup>)</td>
<td align="char" char=".">0.88</td>
<td align="char" char=".">0.154</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>A positive Sen&#xb4;s slope reveals increasing values, whereas a negative slope indicates decreasing values.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Mann-Kendall trend test results for seasonal and monthly nutrient loads.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left"/>
<th colspan="2" align="center">
<italic>NO3</italic>
<sup>
<italic>-</italic>
</sup>
</th>
<th colspan="2" align="center">
<italic>TP</italic>
</th>
<th colspan="2" align="center">PP</th>
<th colspan="2" align="center">DP</th>
</tr>
<tr>
<th align="left">Variable</th>
<th colspan="2" align="center">Sen&#xb4;s slope</th>
<th colspan="2" align="center">Sen&#xb4;s slope</th>
<th colspan="2" align="center">Sen&#xb4;s slope</th>
<th colspan="2" align="center">Sen&#xb4;s slope</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Autumn</td>
<td align="char" char=".">&#x2212;1.86<sup>&#x2b;</sup>
</td>
<td align="char" char=".">&#x2212;0.400</td>
<td align="char" char=".">&#x2212;1.31</td>
<td align="char" char=".">0.003</td>
<td align="char" char=".">&#x2212;1.42</td>
<td align="char" char=".">&#x2212;0.002</td>
<td align="char" char=".">&#x2212;1.09</td>
<td align="char" char=".">0.000</td>
</tr>
<tr>
<td align="left">Winter</td>
<td align="char" char=".">0.77</td>
<td align="char" char=".">0.194</td>
<td align="char" char=".">0.99</td>
<td align="char" char=".">0.003</td>
<td align="char" char=".">1.53</td>
<td align="char" char=".">0.002</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">0.000</td>
</tr>
<tr>
<td align="left">Spring</td>
<td align="char" char=".">&#x2212;0.33</td>
<td align="char" char=".">&#x2212;0.040</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">0.000</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">0.000</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">0.000</td>
</tr>
<tr>
<td align="left">Summer</td>
<td align="char" char=".">&#x2212;1.09</td>
<td align="char" char=".">&#x2212;0.052</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">0.000</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">0.000</td>
<td align="char" char=".">&#x2212;1.64</td>
<td align="char" char=".">0.000</td>
</tr>
<tr>
<td align="left">October</td>
<td align="char" char=".">&#x2212;2.19&#x2a;</td>
<td align="char" char=".">&#x2212;114.638</td>
<td align="char" char=".">&#x2212;1.42</td>
<td align="char" char=".">&#x2212;0.582</td>
<td align="char" char=".">&#x2212;1.42</td>
<td align="char" char=".">&#x2212;0.540</td>
<td align="char" char=".">&#x2212;1.86</td>
<td align="char" char=".">&#x2212;0.136</td>
</tr>
<tr>
<td align="left">November</td>
<td align="char" char=".">&#x2212;1.53</td>
<td align="char" char=".">&#x2212;199.571</td>
<td align="char" char=".">&#x2212;1.64</td>
<td align="char" char=".">&#x2212;1.773</td>
<td align="char" char=".">&#x2212;1.86<sup>&#x2b;</sup>
</td>
<td align="char" char=".">&#x2212;1.376</td>
<td align="char" char=".">&#x2212;1.31</td>
<td align="char" char=".">&#x2212;0.381</td>
</tr>
<tr>
<td align="left">December</td>
<td align="char" char=".">&#x2212;1.64</td>
<td align="char" char=".">&#x2212;372.552</td>
<td align="char" char=".">&#x2212;1.09</td>
<td align="char" char=".">&#x2212;2.275</td>
<td align="char" char=".">&#x2212;1.09</td>
<td align="char" char=".">&#x2212;0.628</td>
<td align="char" char=".">&#x2212;0.77</td>
<td align="char" char=".">&#x2212;0.387</td>
</tr>
<tr>
<td align="left">January</td>
<td align="char" char=".">0.55</td>
<td align="char" char=".">72.032</td>
<td align="char" char=".">1.31</td>
<td align="char" char=".">1.558</td>
<td align="char" char=".">0.99</td>
<td align="char" char=".">1.253</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">0.078</td>
</tr>
<tr>
<td align="left">February</td>
<td align="char" char=".">0.88</td>
<td align="char" char=".">104.815</td>
<td align="char" char=".">1.31</td>
<td align="char" char=".">1.963</td>
<td align="char" char=".">1.53</td>
<td align="char" char=".">2.083</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">0.032</td>
</tr>
<tr>
<td align="left">March</td>
<td align="char" char=".">0.44</td>
<td align="char" char=".">80.935</td>
<td align="char" char=".">0.44</td>
<td align="char" char=".">0.486</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">0.224</td>
<td align="char" char=".">0.77</td>
<td align="char" char=".">0.240</td>
</tr>
<tr>
<td align="left">April</td>
<td align="char" char=".">&#x2212;0.22</td>
<td align="char" char=".">&#x2212;28.660</td>
<td align="char" char=".">0.99</td>
<td align="char" char=".">1.038</td>
<td align="char" char=".">0.55</td>
<td align="char" char=".">0.243</td>
<td align="char" char=".">0.99</td>
<td align="char" char=".">0.206</td>
</tr>
<tr>
<td align="left">May</td>
<td align="char" char=".">&#x2212;1.42</td>
<td align="char" char=".">&#x2212;63.505</td>
<td align="char" char=".">&#x2212;0.77</td>
<td align="char" char=".">&#x2212;0.349</td>
<td align="char" char=".">&#x2212;1.09</td>
<td align="char" char=".">&#x2212;0.289</td>
<td align="char" char=".">0.11</td>
<td align="char" char=".">0.007</td>
</tr>
<tr>
<td align="left">June</td>
<td align="char" char=".">&#x2212;0.88</td>
<td align="char" char=".">&#x2212;47.283</td>
<td align="char" char=".">&#x2212;0.44</td>
<td align="char" char=".">&#x2212;0.092</td>
<td align="char" char=".">&#x2212;0.88</td>
<td align="char" char=".">&#x2212;0.172</td>
<td align="char" char=".">&#x2212;0.44</td>
<td align="char" char=".">&#x2212;0.029</td>
</tr>
<tr>
<td align="left">July</td>
<td align="char" char=".">&#x2212;0.66</td>
<td align="char" char=".">&#x2212;32.810</td>
<td align="char" char=".">&#x2212;0.66</td>
<td align="char" char=".">&#x2212;0.156</td>
<td align="char" char=".">&#x2212;0.66</td>
<td align="char" char=".">&#x2212;0.072</td>
<td align="char" char=".">&#x2212;1.09</td>
<td align="char" char=".">&#x2212;0.072</td>
</tr>
<tr>
<td align="left">August</td>
<td align="char" char=".">&#x2212;1.09</td>
<td align="char" char=".">&#x2212;16.652</td>
<td align="char" char=".">&#x2212;0.99</td>
<td align="char" char=".">&#x2212;0.074</td>
<td align="char" char=".">&#x2212;0.33</td>
<td align="char" char=".">&#x2212;0.022</td>
<td align="char" char=".">&#x2212;1.97&#x2a;</td>
<td align="char" char=".">&#x2212;0.059</td>
</tr>
<tr>
<td align="left">September</td>
<td align="char" char=".">&#x2212;1.86<sup>&#x2b;</sup>
</td>
<td align="char" char=".">&#x2212;34.130</td>
<td align="char" char=".">&#x2212;1.86<sup>&#x2b;</sup>
</td>
<td align="char" char=".">&#x2212;0.114</td>
<td align="char" char=".">&#x2212;1.53</td>
<td align="char" char=".">&#x2212;0.044</td>
<td align="char" char=".">&#x2212;2.08&#x2a;</td>
<td align="char" char=".">&#x2212;0.056</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>A positive Sen&#xb4;s slope reveals increasing values, whereas a negative slope indicates decreasing values. Significance level: &#x2a;: 0.05; <sup>&#x2b;</sup>: 0.1.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Mann-Kendall trend test results for seasonal and monthly nutrient concentrations.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left"/>
<th colspan="2" align="center">
<italic>NO3-</italic>
</th>
<th colspan="2" align="center">
<italic>TP</italic>
</th>
<th colspan="2" align="center">PP</th>
<th colspan="2" align="center">DP</th>
</tr>
<tr>
<th align="left">Variable</th>
<th colspan="2" align="center">Sen&#xb4;s slope</th>
<th colspan="2" align="center">Sen&#xb4;s slope</th>
<th colspan="2" align="center">Sen&#xb4;s slope</th>
<th colspan="2" align="center">Sen&#xb4;s slope</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Autumn</td>
<td align="char" char=".">&#x2212;1.31</td>
<td align="char" char=".">&#x2212;0.080</td>
<td align="char" char=".">&#x2212;0.33</td>
<td align="char" char=".">&#x2212;0.135</td>
<td align="char" char=".">&#x2212;0.55</td>
<td align="char" char=".">&#x2212;0.342</td>
<td align="char" char=".">0.22</td>
<td align="char" char=".">0.093</td>
</tr>
<tr>
<td align="left">Winter</td>
<td align="char" char=".">0.11</td>
<td align="char" char=".">0.007</td>
<td align="char" char=".">2.30&#x2a;</td>
<td align="char" char=".">1.006</td>
<td align="char" char=".">1.53</td>
<td align="char" char=".">0.832</td>
<td align="char" char=".">0.77</td>
<td align="char" char=".">0.168</td>
</tr>
<tr>
<td align="left">Spring</td>
<td align="char" char=".">&#x2212;0.22</td>
<td align="char" char=".">&#x2212;0.003</td>
<td align="char" char=".">0.44</td>
<td align="char" char=".">0.228</td>
<td align="char" char=".">&#x2212;0.11</td>
<td align="char" char=".">&#x2212;0.023</td>
<td align="char" char=".">1.53</td>
<td align="char" char=".">0.150</td>
</tr>
<tr>
<td align="left">Summer</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">0.006</td>
<td align="char" char=".">1.42</td>
<td align="char" char=".">0.393</td>
<td align="char" char=".">1.31</td>
<td align="char" char=".">0.292</td>
<td align="char" char=".">&#x2212;0.22</td>
<td align="char" char=".">&#x2212;0.056</td>
</tr>
<tr>
<td align="left">October</td>
<td align="char" char=".">&#x2212;1.53</td>
<td align="char" char=".">&#x2212;0.088</td>
<td align="char" char=".">&#x2212;0.77</td>
<td align="char" char=".">&#x2212;0.855</td>
<td align="char" char=".">&#x2212;1.20</td>
<td align="char" char=".">&#x2212;1.135</td>
<td align="char" char=".">&#x2212;0.22</td>
<td align="char" char=".">&#x2212;0.137</td>
</tr>
<tr>
<td align="left">November</td>
<td align="char" char=".">&#x2212;1.75<sup>&#x2b;</sup>
</td>
<td align="char" char=".">&#x2212;0.124</td>
<td align="char" char=".">&#x2212;0.66</td>
<td align="char" char=".">&#x2212;0.873</td>
<td align="char" char=".">&#x2212;1.09</td>
<td align="char" char=".">&#x2212;1.239</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">&#x2212;0.024</td>
</tr>
<tr>
<td align="left">December</td>
<td align="char" char=".">&#x2212;0.66</td>
<td align="char" char=".">&#x2212;0.084</td>
<td align="char" char=".">0.77</td>
<td align="char" char=".">0.458</td>
<td align="char" char=".">0.11</td>
<td align="char" char=".">0.173</td>
<td align="char" char=".">1.31</td>
<td align="char" char=".">0.285</td>
</tr>
<tr>
<td align="left">January</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">0.011</td>
<td align="char" char=".">1.31</td>
<td align="char" char=".">1.659</td>
<td align="char" char=".">0.77</td>
<td align="char" char=".">0.872</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">&#x2212;0.012</td>
</tr>
<tr>
<td align="left">February</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">&#x2212;0.011</td>
<td align="char" char=".">1.20</td>
<td align="char" char=".">0.939</td>
<td align="char" char=".">0.99</td>
<td align="char" char=".">0.696</td>
<td align="char" char=".">0.66</td>
<td align="char" char=".">0.165</td>
</tr>
<tr>
<td align="left">March</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">0.002</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">&#x2212;0.076</td>
<td align="char" char=".">&#x2212;0.33</td>
<td align="char" char=".">&#x2212;0.269</td>
<td align="char" char=".">1.42</td>
<td align="char" char=".">0.244</td>
</tr>
<tr>
<td align="left">April</td>
<td align="char" char=".">&#x2212;1.20</td>
<td align="char" char=".">&#x2212;0.031</td>
<td align="char" char=".">0.11</td>
<td align="char" char=".">0.144</td>
<td align="char" char=".">0.11</td>
<td align="char" char=".">0.128</td>
<td align="char" char=".">1.31</td>
<td align="char" char=".">0.165</td>
</tr>
<tr>
<td align="left">May</td>
<td align="char" char=".">&#x2212;0.22</td>
<td align="char" char=".">&#x2212;0.007</td>
<td align="char" char=".">0.33</td>
<td align="char" char=".">0.129</td>
<td align="char" char=".">&#x2212;0.33</td>
<td align="char" char=".">&#x2212;0.320</td>
<td align="char" char=".">1.42</td>
<td align="char" char=".">0.251</td>
</tr>
<tr>
<td align="left">June</td>
<td align="char" char=".">&#x2212;0.22</td>
<td align="char" char=".">&#x2212;0.013</td>
<td align="char" char=".">0.22</td>
<td align="char" char=".">0.209</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">&#x2212;0.005</td>
<td align="char" char=".">1.75<sup>&#x2b;</sup>
</td>
<td align="char" char=".">0.136</td>
</tr>
<tr>
<td align="left">July</td>
<td align="char" char=".">0.55</td>
<td align="char" char=".">0.014</td>
<td align="char" char=".">0.22</td>
<td align="char" char=".">0.215</td>
<td align="char" char=".">0.00</td>
<td align="char" char=".">0.016</td>
<td align="char" char=".">&#x2212;0.22</td>
<td align="char" char=".">&#x2212;0.031</td>
</tr>
<tr>
<td align="left">August</td>
<td align="char" char=".">0.44</td>
<td align="char" char=".">0.033</td>
<td align="char" char=".">0.66</td>
<td align="char" char=".">0.208</td>
<td align="char" char=".">1.42</td>
<td align="char" char=".">0.284</td>
<td align="char" char=".">&#x2212;0.33</td>
<td align="char" char=".">&#x2212;0.064</td>
</tr>
<tr>
<td align="left">September</td>
<td align="char" char=".">0.44</td>
<td align="char" char=".">0.017</td>
<td align="char" char=".">1.31</td>
<td align="char" char=".">0.196</td>
<td align="char" char=".">0.88</td>
<td align="char" char=".">0.133</td>
<td align="char" char=".">0.77</td>
<td align="char" char=".">0.220</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>A positive Sen&#xb4;s slope reveals increasing values, whereas a negative slope indicates decreasing values. Significance level: &#x2a;: 0.05; <sup>&#x2b;</sup>: 0.1.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>NO<sub>3</sub>
<sup>&#x2212;</sup> and P loads were similar or lower compared with other studies carried out in Galicia (<xref ref-type="bibr" rid="B5">Bern&#xe1;rdez et al., 2013</xref>; <xref ref-type="bibr" rid="B46">Serrano et al., 2015</xref>) and suggest that agricultural activities and forest plantation management in the study area have only a slight impact on surface waters of the Corbeira catchment relative to other mixed land use catchments in Europe (<xref ref-type="bibr" rid="B20">Fovet et al., 2015</xref>; <xref ref-type="bibr" rid="B16">Dupas et al., 2017</xref>; <xref ref-type="bibr" rid="B7">Bouraoui and Malago, 2020</xref>; <xref ref-type="bibr" rid="B17">Ehrhardt et al., 2021</xref>). In fact, annual mean concentrations of both elements were well below the recommendations of WFD, the threshold proposed by <xref ref-type="bibr" rid="B8">Camargo and Alonso (2006)</xref> for protecting the most sensitive freshwater species and the threshold stablished by <xref ref-type="bibr" rid="B9">Carleton et al. (2009)</xref> to prevent nuisance level of periphyton and cyanobacteria (<xref ref-type="table" rid="T4">Table 4</xref>). When the mean nitrate and phosphorus concentrations obtained in this study are compared to the threshold values reported by <xref ref-type="bibr" rid="B41">Poikane et al. (2021)</xref>, based on their relationship with macrophytes and phytobenthos, it is observed that mean nutrient concentrations of this study are within the range of predicted nitrogen good-moderate class boundaries for low alkalinity lowland and upland rivers (<xref ref-type="table" rid="T4">Table 4</xref>). They also are in the boundary values proposed by <xref ref-type="bibr" rid="B12">de Vries et al. (2013)</xref> to protect aquatic ecosystems from eutrophication. However, they exceed the threshold reported by Stevenson et al. (2008). In this case, the values reported by this last author coincide with concentrations in the catchment reference site (TP: 15 &#xb1; 7&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup>; mean and standard deviation from 120 water sampled collected from 2004 to 2016, unpublished data); i.e., in the headwaters of the studied stream, which drains a forest area and is assumed to be free of nutrient inputs from human activities. So, protection of stream with this threshold was deemed impractical. This highlights the need to adapt the thresholds to individual water bodies to protect aquatic ecosystems from eutrophication.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Nutrient thresholds establish by different ecology-based strategies.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Reference</th>
<th align="center">Nutrient criteria</th>
<th align="center">Approach for establishing criteria</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">
<xref ref-type="bibr" rid="B8">Camargo and Alonso (2006)</xref>
</td>
<td align="center">0.5&#x2013;1.0&#xa0;mg&#xa0;TN L<sup>&#x2212;1</sup>
</td>
<td align="left">Protect aquatic ecosystems from eutrophication and acidification</td>
</tr>
<tr>
<td align="left">
<xref ref-type="bibr" rid="B9">Carleton et al. (2009)</xref>
</td>
<td align="center">100&#xa0;&#x3bc;g&#xa0;TP L<sup>&#x2212;1</sup>; 2.7&#xa0;mg&#xa0;TN&#xa0;L<sup>&#x2212;1</sup>
</td>
<td align="left">Prevent nuisance level of periphyton and cyanobacteria</td>
</tr>
<tr>
<td align="left">
<xref ref-type="bibr" rid="B41">Poikane et al. (2021)</xref>
</td>
<td align="center">1.63&#xa0;mg&#xa0;TN L<sup>&#x2212;1</sup> (.71&#x2013;4.19), 32&#xa0;&#x3bc;g SRP L<sup>&#x2212;1</sup> (18&#x2013;58)</td>
<td align="left">Good moderate class boundary low alkalinity lowland rivers</td>
</tr>
<tr>
<td align="left"/>
<td align="center">1.79&#xa0;mg&#xa0;N L<sup>&#x2212;1</sup> (1.02&#x2013;4.04), 40&#xa0;&#x3bc;g&#xa0;L<sup>&#x2212;1</sup> (19&#x2013;78)</td>
<td align="left">Good moderate class boundary low alkalinity upland rivers</td>
</tr>
<tr>
<td align="left">Stevenson et al. (2008)</td>
<td align="center">10&#x2013;12&#xa0;&#x3bc;g&#xa0;TP L<sup>&#x2212;1</sup>
</td>
<td align="left">Protect high-quality biological conditions in streams (Mid-Atlantic Highlands)</td>
</tr>
<tr>
<td align="left">
<xref ref-type="bibr" rid="B12">de Vries et al. (2013)</xref>
</td>
<td align="center">1.0&#x2013;2.5&#xa0;mg&#xa0;N L<sup>&#x2212;1</sup>
</td>
<td align="left">Protect aquatic ecosystems from eutrophication and acidification</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>TP, Total phosphorus; SRP, soluble phosphorus; TN, total nitrogen.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Nitrate concentration exceeds the concentration collected by <xref ref-type="bibr" rid="B32">Mac&#xed;as et al. (1991)</xref> as reference values in slightly polluted rivers (1.19&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>), the background concentration proposed by <xref ref-type="bibr" rid="B33">Meybeck (1982)</xref> for the major unpolluted rivers (0.44&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>), as well as the levels suggested for uncontaminated natural rivers in Galicia (1&#x2013;4&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>, <xref ref-type="bibr" rid="B2">Antelo and Arce, 1996</xref>) and the mean concentration in the catchment reference site (1.99 &#xb1; 0.54&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>; mean and standard deviation from 120 water sampled collected from 2004 to 2016, unpublished data). <xref ref-type="bibr" rid="B37">Naves et al. (2019)</xref> observed nitrate contamination in the groundwaters (concentrations between 2.54 and 47.1&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>) of the study area linked to an inadequate management of the manure in the field and even to an occasional discharge of slurry from local farms. However, the results of this study are far from what was observed by these authors, since nitrate concentrations in summer (mean &#x3d; 4.13&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>) -when discharge is dominated by groundwater-are well below the values reported by <xref ref-type="bibr" rid="B37">Naves et al. (2019)</xref> as a reference in the study area (1.2&#x2013;9.9&#xa0;mg&#xa0;L<sup>&#x2212;1</sup>).</p>
<p>Loads and the concentrations of both nutrients showed a marked seasonal pattern related to the hydrological behaviour of the catchment, that is, the highest values were recorded in winter and autumn, coinciding with the highest discharge and runoff events (<xref ref-type="bibr" rid="B43">Rodr&#xed;guez-Blanco et al., 2019</xref>), while the lowest were obtained in summer, coinciding with the period of low waters and the greatest influence from the baseflow. This behaviour suggests a greater nutrient mobilization during the wet period (October-March), which is consistent with a predominant diffuse source of nitrate and phosphorus (<xref ref-type="bibr" rid="B1">Abbott et al., 2018</xref>; <xref ref-type="bibr" rid="B15">Dupas et al., 2019</xref>; <xref ref-type="bibr" rid="B21">Garzon-Vidueira et al., 2020</xref>). Thus, wet periods increase the connectivity between nutrient sources and the stream network (<xref ref-type="bibr" rid="B44">Rodr&#xed;guez-Blanco et al., 2013</xref>; <xref ref-type="bibr" rid="B45">2015</xref>), promoting the transfer of nutrients. Corbeira is a transport-limited catchment, where the nutrient fluxes are controlled by the hydrology, as the flow is the main driver of nutrient transport, since it controls the connectivity between the source (e.g., soils) and the stream. Nitrate is mainly transported in subsurface flow (<xref ref-type="bibr" rid="B45">Rodr&#xed;guez-Blanco et al., 2015</xref>) while P (mainly associated with sediment) in surface runoff (<xref ref-type="bibr" rid="B44">Rodr&#xed;guez-Blanco et al., 2013</xref>).</p>
<p>The results showed the lack of consistent trends of N and P over time, although a decreasing trend appears for N loads, mainly linked to the reduction of nitrate transport during autumn and spring (<xref ref-type="table" rid="T2">Table 2</xref>), although it might also be linked to a decline in agricultural activities in the catchment, which leads to a decrease in the application of manure and fertilisers&#x2014;ideally in the autumn and spring. Forest plantation cutting cycles should also be considered since forest cut decrease nutrient uptake, making more nutrients available for leaching into water courses (Thompson et al., 2009). Forest roads, especially when intercepting the stream can be an important source of nutrients associated (mainly P) with sediments in plantation areas. Phosphorus points to an upward trend, mainly in winter, which could be related to the positive flow trend in the study area (<xref ref-type="bibr" rid="B42">Rodr&#xed;guez-Blanco et al., 2020</xref>). An increase in discharge levels translates into a greater flow velocity and movement of sediments, which is essential for the transport of P in the catchment since it is intrinsically related to the particulate fraction. A similar pattern was observed in another mixed land use catchment (<xref ref-type="bibr" rid="B34">Minaudo et al., 2015</xref>; <xref ref-type="bibr" rid="B6">Bol et al., 2018</xref>), demonstrating the significance of suspended sediment in the dynamics of P. The results also show an increase in the nutrient concentrations, despite the fact that discharge and nutrient load decrease in the summer, meaning that the flow is being enriched owing to lesser dilution. As has been reported in several other catchments with forestry and agricultural uses (e.g., <xref ref-type="bibr" rid="B1">Abbott et al., 2018</xref>), this could be indicative of surplus nutrients applied to the soil being available for export to water bodies by drainage or runoff.</p>
<p>Despite the decline in population and reduced agricultural activities in the catchment, only a limited improvement in nitrate concentrations was observed. This may be due to the potentially long transit time in the catchment, so reducing non-point source nutrients may take decades to affect the water chemistry of rivers (<xref ref-type="bibr" rid="B20">Fovet et al., 2015</xref>; <xref ref-type="bibr" rid="B16">Dupas et al., 2017</xref>). The increase in winter flow could also counteract the effect of the reduction in the number of people, so that no trend at annual scale will be detected.</p>
</sec>
<sec id="s5">
<title>Concluding Remarks</title>
<p>The concentrations of N and P measured in the catchment are relatively low, which shows that agricultural practices and the forest plantation management in the catchment do not jeopardise (for the moment) the water quality, although the catchment nutrient reference levels were widely exceeded. The variability of export fluxes is high, despite their &#x201c;natural nature,&#x201d; due to the fact that connectivity between the nutrient source area and the stream is dependent on flow conditions. A lack of consistent trends of N and P over time was observed, but the nutrients analysed behaved differently. Thus, N appeared to show a downward trend, mainly associated with declining transport in autumn and spring. Phosphorus showed an upward trend in winter linked to a greater flow and sediment transfer. As a result, any increase in erosion and runoff will almost certainly result in an increase phosphorus load in the study area.</p>
<p>Despite the decrease in population and abandonment of agricultural land in the catchment, the concentrations of nitrate do not show a sizable reduction, which could be indicative of surplus nutrients provided to the soil being available for export to water bodies <italic>via</italic> drainage or runoff. As a result, a decrease in inputs from various sources does not result in rapid improvement in riverine water quality.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data Availability Statement</title>
<p>The datasets presented in this article are not readily available because the data are part of a large database. Requests to access the datasets should be directed to <email>maria.luz.rodriguez.blanco@uvigo.es</email>.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>This research was carried out within the projects REN2003-08143, funded by the Spanish Ministry of Education and Science, and PGIDIT05RAG10303PR and 10MDS103031PR, financed by the Xunta of Galicia.</p>
</sec>
<ack>
<p>We would like to thank two reviewers for their helpful comments and suggestions.</p>
</ack>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abbott</surname>
<given-names>B. W.</given-names>
</name>
<name>
<surname>Moatar</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Gauthier</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Fovet</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Antoine</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Ragueneau</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Trends and Seasonality of River Nutrients in Agricultural Catchments: 18 Years of Weekly Citizen Science in France</article-title>. <source>Sci. Total Environ.</source> <volume>624</volume>, <fpage>845</fpage>&#x2013;<lpage>858</lpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2017.12.176</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Antelo</surname>
<given-names>J. M.</given-names>
</name>
<name>
<surname>Arce</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>1996</year>). &#x201c;<article-title>Caracter&#xed;sticas fisicoqu&#xed;micas das augas superficiais</article-title>,&#x201d; in <source>As Aguas de Galicia</source>. Editor <person-group person-group-type="editor">
<name>
<surname>D&#xed;az-Fierros</surname>
<given-names>F.</given-names>
</name>
</person-group> (<publisher-loc>Santiago de Compostela, Spain</publisher-loc>: <publisher-name>Consello da Cultura Galega</publisher-name>), <fpage>353</fpage>&#x2013;<lpage>446</lpage>.</citation>
</ref>
<ref id="B3">
<citation citation-type="book">
<collab>APHA</collab> (<year>1998</year>). <source>Standard Methods for the Examination of Water and Wastewater 20th Edition</source>. <publisher-loc>Washington, DC</publisher-loc>: <publisher-name>American Public Health Association, American Water Work Association, Water Environment Federation</publisher-name>.</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bennett</surname>
<given-names>E. M.</given-names>
</name>
<name>
<surname>Carpenter</surname>
<given-names>S. R.</given-names>
</name>
<name>
<surname>Caraco</surname>
<given-names>N. F.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Human Impact on Erodable Phosphorus and Eutrophication: A Global Perspective</article-title>. <source>Bioscience</source> <volume>51</volume> (<issue>3</issue>), <fpage>227</fpage>&#x2013;<lpage>234</lpage>. <pub-id pub-id-type="doi">10.1641/0006-3568(2001)051[0227:hioepa]2.0.co;2</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bern&#xe1;rdez</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Ospina-&#xc1;lvarez</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Caetano</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Prego</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Fluvial Contributions of Nutrient Salts, Dissolved Trace Elements and Organic Carbon to the Sea by Pristine Temperate Rivers (SW Europe)</article-title>. <source>Environ. Chem.</source> <volume>10</volume> (<issue>1</issue>), <fpage>42</fpage>&#x2013;<lpage>53</lpage>. <pub-id pub-id-type="doi">10.1071/en12123</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bol</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Gruau</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Mellander</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Dupas</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Bechmann</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Skarb&#xf8;vik</surname>
<given-names>E.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Challenges of Reducing Phosphorus Based Water Eutrophication in the Agricultural Landscapes of Northwest Europe</article-title>. <source>Front. Mar. Sci.</source> <volume>5</volume>, <fpage>276</fpage>. <pub-id pub-id-type="doi">10.3389/fmars.2018.00276</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bouraoui</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Malag&#xf3;</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Trend Analysis of Nitrate Concentration in Rivers in Southern France</article-title>. <source>Water</source> <volume>12</volume> (<issue>12</issue>), <fpage>3374</fpage>. <pub-id pub-id-type="doi">10.3390/w12123374</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Camargo</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Alonso</surname>
<given-names>&#xc1;.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Ecological and Toxicological Effects of Inorganic Nitrogen Pollution in Aquatic Ecosystems: a Global Assessment</article-title>. <source>Environ. Int.</source> <volume>32</volume> (<issue>6</issue>), <fpage>831</fpage>&#x2013;<lpage>849</lpage>. <pub-id pub-id-type="doi">10.1016/j.envint.2006.05.002</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Carleton</surname>
<given-names>J. N.</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>R. A.</given-names>
</name>
<name>
<surname>Clough</surname>
<given-names>J. S.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Ecosystem Modeling Applied to Nutrient Criteria Development in Rivers</article-title>. <source>Environ. Manage</source> <volume>44</volume> (<issue>3</issue>), <fpage>485</fpage>&#x2013;<lpage>492</lpage>. <pub-id pub-id-type="doi">10.1007/s00267-009-9344-2</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Carvalho</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Mackay</surname>
<given-names>E. B.</given-names>
</name>
<name>
<surname>Cardoso</surname>
<given-names>A. C.</given-names>
</name>
<name>
<surname>Baattrup-Pedersen</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Birk</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Blackstock</surname>
<given-names>K. L.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Protecting and Restoring Europe&#x27;s Waters: An Analysis of the Future Development Needs of the Water Framework Directive</article-title>. <source>Sci. Total Environ.</source> <volume>658</volume>, <fpage>1228</fpage>&#x2013;<lpage>1238</lpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2018.12.255</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>de Vries</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Kros</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Kroeze</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Seitzinger</surname>
<given-names>S. P.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Assessing Planetary and Regional Nitrogen Boundaries Related to Food Security and Adverse Environmental Impacts</article-title>. <source>Curr. Opin. Environ. Sustain.</source> <volume>5</volume> (<issue>3-4</issue>), <fpage>392</fpage>&#x2013;<lpage>402</lpage>. <pub-id pub-id-type="doi">10.1016/j.cosust.2013.07.004</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Directive</surname>
<given-names>W. F.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Water Framework Directive</article-title>. <source>J. Ref. OJL</source> <volume>327</volume>, <fpage>1</fpage>&#x2013;<lpage>73</lpage>.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cosentino</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Malek</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Maiorano</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Thuiller</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Verburg</surname>
<given-names>P. H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A New European Land Systems Representation Accounting for Landscape Characteristics</article-title>. <source>Landsc. Ecol.</source> <volume>36</volume> (<issue>8</issue>), <fpage>2215</fpage>&#x2013;<lpage>2234</lpage>. <pub-id pub-id-type="doi">10.1007/s10980-021-01227-5</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dupas</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Abbott</surname>
<given-names>B. W.</given-names>
</name>
<name>
<surname>Minaudo</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Fovet</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Distribution of Landscape Units within Catchments Influences Nutrient Export Dynamics</article-title>. <source>Front. Environ. Sci.</source> <volume>7</volume>, <fpage>43</fpage>. <pub-id pub-id-type="doi">10.3389/fenvs.2019.00043</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dupas</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Mellander</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Gascuel-Odoux</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Fovet</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>McAleer</surname>
<given-names>E. B.</given-names>
</name>
<name>
<surname>McDonald</surname>
<given-names>N. T.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>The Role of Mobilisation and Delivery Processes on Contrasting Dissolved Nitrogen and Phosphorus Exports in Groundwater Fed Catchments</article-title>. <source>Sci. Total Environ.</source> <volume>599</volume>, <fpage>1275</fpage>&#x2013;<lpage>1287</lpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2017.05.091</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ehrhardt</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ebeling</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Dupas</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Kumar</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Fleckenstein</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Musolff</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Nitrate Transport and Retention in Western European Catchments Are Shaped by Hydroclimate and Subsurface Properties</article-title>. <source>Water Resour. Res.</source> <volume>57</volume> (<issue>10</issue>), <fpage>e2020WR029469</fpage>. <pub-id pub-id-type="doi">10.1029/2020wr029469</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ekholm</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Rankinen</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Rita</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>R&#xe4;ike</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Sj&#xf6;blom</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Raateland</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Phosphorus and Nitrogen Fluxes Carried by 21 Finnish Agricultural Rivers in 1985&#x2013;2006</article-title>. <source>Environ. Monit. Assess.</source> <volume>187</volume> (<issue>4</issue>), <fpage>216</fpage>&#x2013;<lpage>217</lpage>. <pub-id pub-id-type="doi">10.1007/s10661-015-4417-6</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="book">
<collab>European Environment Agency</collab> (<year>2018</year>). <source>European Waters Assessment. EEA Report N&#x00B0; 7/2018</source>. <publisher-loc>Copenhage, Denmark</publisher-loc>: <publisher-name>EEA</publisher-name>.</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fovet</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Ruiz</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Faucheux</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Mol&#xe9;nat</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Sekhar</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Vert&#xe8;s</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Using Long Time Series of Agricultural-Derived Nitrates for Estimating Catchment Transit Times</article-title>. <source>J. Hydrology</source> <volume>522</volume>, <fpage>603</fpage>&#x2013;<lpage>617</lpage>. <pub-id pub-id-type="doi">10.1016/j.jhydrol.2015.01.030</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garzon-Vidueira</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Rial-Otero</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Garcia-Nocelo</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Rivas-Gonzalez</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Moure-Gonzalez</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Fompedri&#xf1;a-Roca</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Identification of Nitrates Origin in Limia River Basin and Pollution-Determinant Factors. Agric</article-title>. <source>Ecosyst. Environ.</source> <volume>290</volume>, <fpage>106775</fpage>.</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grizzetti</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Vigiak</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Udias</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Aloe</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Zanni</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bouraoui</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>How EU Policies Could Reduce Nutrient Pollution in European Inland and Coastal Waters</article-title>. <source>Glob. Environ. Change</source> <volume>69</volume>, <fpage>102281</fpage>. <pub-id pub-id-type="doi">10.1016/j.gloenvcha.2021.102281</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Gustard</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Bullock</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Dixon</surname>
<given-names>J. M.</given-names>
</name>
</person-group> (<year>1992</year>). <source>Low Flow Estimation in the United Kingdom</source>. <publisher-loc>England</publisher-loc>: <publisher-name>Institute of Hydrology</publisher-name>, <fpage>19</fpage>.</citation>
</ref>
<ref id="B24">
<citation citation-type="book">
<collab>IGME (Instituto Tecnol&#xf3;gico Geominero de Espa&#xf1;a)</collab> (<year>1981</year>). <source>Mapa Geol&#xf3;gico de Espa&#xf1;a, 1:50,000</source>. <publisher-loc>Betanzos Spain</publisher-loc>. <comment>Hoja 45</comment>.</citation>
</ref>
<ref id="B25">
<citation citation-type="book">
<collab>INE</collab> (<year>2019</year>). <source>Nomenclator: Poblaci&#xf3;n del padr&#xf3;n continuo por unidad poblacional</source>. <publisher-loc>Madrid</publisher-loc>: <publisher-name>Instituto Estad&#xed;stico de Espa&#xf1;a</publisher-name>.</citation>
</ref>
<ref id="B26">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Iuss Working Group Wrb</surname>
</name>
</person-group> (<year>2014</year>). in <source>International Soil Clasification System for Naming Soils and Creating Legends for Soil Maps</source> (<publisher-loc>Rome</publisher-loc>: <publisher-name>World Soil Resources Reports. N&#x00B0; 106FAO</publisher-name>).<source>World Reference Base for Soil Resources.</source>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kelly</surname>
<given-names>M. G.</given-names>
</name>
<name>
<surname>Phillips</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Teixeira</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>V&#xe1;rb&#xed;r&#xf3;</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Herrero</surname>
<given-names>F. S.</given-names>
</name>
<name>
<surname>Willby</surname>
<given-names>N. J.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Establishing Ecologically-Relevant Nutrient Thresholds: A Tool-Kit with Guidance on its Use</article-title>. <source>Sci. Total Environ.</source> <volume>807</volume>, <fpage>150977</fpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2021.150977</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumar</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Merwade</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Kam</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Thurner</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Streamflow Trends in Indiana: Effects of Long Term Persistence, Precipitation and Subsurface Drains</article-title>. <source>J. Hydroloy</source> <volume>374</volume> (<issue>1-2</issue>), <fpage>171</fpage>&#x2013;<lpage>183</lpage>. <pub-id pub-id-type="doi">10.1016/j.jhydrol.2009.06.012</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lassaletta</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Garc&#xed;a-G&#xf3;mez</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Gimeno</surname>
<given-names>B. S.</given-names>
</name>
<name>
<surname>Rovira</surname>
<given-names>J. V.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Agriculture-induced Increase in Nitrate Concentrations in Stream Waters of a Large Mediterranean Catchment over 25 Years (1981&#x2013;2005)</article-title>. <source>Sci. Total Environ.</source> <volume>407</volume> (<issue>23</issue>), <fpage>6034</fpage>&#x2013;<lpage>6043</lpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2009.08.002</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Linsely</surname>
<given-names>R. K.</given-names>
</name>
<name>
<surname>Kohler</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Paulhus</surname>
<given-names>J. C.</given-names>
</name>
</person-group> (<year>1949</year>). <source>Applied Hydrology</source>. <publisher-loc>New York</publisher-loc>: <publisher-name>McGraw-Hill Book Co</publisher-name>.</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mac&#xed;as</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>&#xc1;lvarez Rodr&#xed;guez</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Calvo de Anta</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Impactos de origen agrario y urbano en la cuenca del r&#xed;o Ulla</article-title>. <source>Ecolog&#xed;a</source> <volume>5</volume>, <fpage>73</fpage>&#x2013;<lpage>86</lpage>.</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Meybeck</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>1982</year>). <article-title>Carbon, Nitrogen, and Phosphorus Transport by World Rivers</article-title>. <source>Am. J. Sci.</source> <volume>282</volume> (<issue>4</issue>), <fpage>401</fpage>&#x2013;<lpage>450</lpage>. <pub-id pub-id-type="doi">10.2475/ajs.282.4.401</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Minaudo</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Meybeck</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Moatar</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Gassama</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Curie</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Eutrophication Mitigation in Rivers: 30 Years of Trends in Spatial and Seasonal Patterns of Biogeochemistry of the Loire River (1980&#x2013;2012)</article-title>. <source>Biogeosciences</source> <volume>12</volume> (<issue>8</issue>), <fpage>2549</fpage>&#x2013;<lpage>2563</lpage>. <pub-id pub-id-type="doi">10.5194/bg-12-2549-2015</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moatar</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Abbott</surname>
<given-names>B. W.</given-names>
</name>
<name>
<surname>Minaudo</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Curie</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Pinay</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Elemental Properties, Hydrology, and Biology Interact to Shape Concentration&#x2010;discharge Curves for Carbon, Nutrients, Sediment, and Major Ions</article-title>. <source>Water Resour. Res.</source> <volume>53</volume> (<issue>2</issue>), <fpage>1270</fpage>&#x2013;<lpage>1287</lpage>. <pub-id pub-id-type="doi">10.1002/2016wr019635</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Murphy</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Riley</surname>
<given-names>J. P.</given-names>
</name>
</person-group> (<year>1962</year>). <article-title>A Modified Single Solution Method for the Determination of Phosphate in Natural Waters</article-title>. <source>Anal. Chim. Acta</source> <volume>27</volume>, <fpage>31</fpage>&#x2013;<lpage>36</lpage>. <pub-id pub-id-type="doi">10.1016/s0003-2670(00)88444-5</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Naves</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Samper</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Mon</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Pisani</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Montenegro</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Carvalho</surname>
<given-names>J. M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Demonstrative Actions of Spring Restoration and Groundwater Protection in Rural Areas of Abegondo (Galicia, Spain)</article-title>. <source>Sustain. Water Resour. Manag.</source> <volume>5</volume> (<issue>1</issue>), <fpage>175</fpage>&#x2013;<lpage>186</lpage>. <pub-id pub-id-type="doi">10.1007/s40899-017-0169-5</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="book">
<collab>OECD</collab> (<year>2007</year>). <source>OECD Regions at a Glance</source>. <publisher-loc>Paris</publisher-loc>: <publisher-name>OECD Publishing</publisher-name>, <fpage>252</fpage>.</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oenema</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>van Liere</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Schoumans</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Effects of Lowering Nitrogen and Phosphorus Surpluses in Agriculture on the Quality of Groundwater and Surface Water in the Netherlands</article-title>. <source>J. Hydrology</source> <volume>304</volume> (<issue>1-4</issue>), <fpage>289</fpage>&#x2013;<lpage>301</lpage>. <pub-id pub-id-type="doi">10.1016/j.jhydrol.2004.07.044</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Panagos</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Borrelli</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Poesen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ballabio</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Lugato</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Meusburger</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>The New Assessment of Soil Loss by Water Erosion in Europe</article-title>. <source>Environ Sci Policy</source> <volume>54</volume>, <fpage>438</fpage>&#x2013;<lpage>447</lpage>. <pub-id pub-id-type="doi">10.1016/j.envsci.2015.08.012</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Poikane</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>V&#xe1;rb&#xed;r&#xf3;</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Kelly</surname>
<given-names>M. G.</given-names>
</name>
<name>
<surname>Birk</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Phillips</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Estimating River Nutrient Concentrations Consistent with Good Ecological Condition: More Stringent Nutrient Thresholds Needed</article-title>. <source>Ecol. Ind.</source> <volume>121</volume>, <fpage>107017</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2020.107017</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodr&#xed;guez-Blanco</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Taboada-Castro</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Taboada-Castro</surname>
<given-names>M. T.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>An Assessment of the Recent Evolution of the Streamflow in a Near-Natural System: A Case Study in the Headwaters of the Mero Basin (Galicia, Spain)</article-title>. <source>Hydrology</source> <volume>7</volume> (<issue>4</issue>), <fpage>97</fpage>. <pub-id pub-id-type="doi">10.3390/hydrology7040097</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodr&#xed;guez-Blanco</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Taboada-Castro</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Taboada-Castro</surname>
<given-names>M. T.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>An Overview of Patterns and Dynamics of Suspended Sediment Transport in an Agroforest Headwater System in Humid Climate: Results from a Long-Term Monitoring</article-title>. <source>Sci. Total Environ.</source> <volume>648</volume>, <fpage>33</fpage>&#x2013;<lpage>43</lpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2018.08.118</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodr&#xed;guez-Blanco</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Taboada-Castro</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Taboada-Castro</surname>
<given-names>M. T.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Phosphorus Transport into a Stream Draining from a Mixed Land Use Catchment in Galicia (NW Spain): Significance of Runoff Events</article-title>. <source>J. Hydrology</source> <volume>481</volume>, <fpage>12</fpage>&#x2013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.1016/j.jhydrol.2012.11.046</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodr&#xed;guez&#x2010;Blanco</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Taboada&#x2010;Castro</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Taboada&#x2010;Castro</surname>
<given-names>M. T.</given-names>
</name>
<name>
<surname>Oropeza&#x2010;Mota</surname>
<given-names>J. L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Relating Nitrogen Export Patterns from a Mixed Land Use Catchment in NW Spain with Rainfall and Streamflow</article-title>. <source>Hydrol. Process</source> <volume>29</volume> (<issue>12</issue>), <fpage>2720</fpage>&#x2013;<lpage>2730</lpage>. <pub-id pub-id-type="doi">10.1002/hyp.10388</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Serrano</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Palleiro</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Rodr&#xed;guez-Blanco</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Taboada-Castro</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Taboada-Castro</surname>
<given-names>M. T.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Preliminary Assessment of the Riverine Nitrogen Concentration and Load in an Agroforestry Basin in NW Spain</article-title>. <source>Commun. Soil Sci. Plant Anal.</source> <volume>46</volume>, <fpage>326</fpage>&#x2013;<lpage>331</lpage>. <pub-id pub-id-type="doi">10.1080/00103624.2014.989071</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sharpley</surname>
<given-names>A. N.</given-names>
</name>
<name>
<surname>Kleinman</surname>
<given-names>P. J.</given-names>
</name>
<name>
<surname>Heathwaite</surname>
<given-names>A. L.</given-names>
</name>
<name>
<surname>Gburek</surname>
<given-names>W. J.</given-names>
</name>
<name>
<surname>Folmar</surname>
<given-names>G. J.</given-names>
</name>
<name>
<surname>Schmidt</surname>
<given-names>J. P.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Phosphorus Loss from an Agricultural Watershed as a Function of Storm Size</article-title>. <source>J. Environ. Qual.</source> <volume>37</volume> (<issue>2</issue>), <fpage>362</fpage>&#x2013;<lpage>368</lpage>. <pub-id pub-id-type="doi">10.2134/jeq2007.0366</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Withers</surname>
<given-names>P. J.</given-names>
</name>
<name>
<surname>Lord</surname>
<given-names>E. I.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Agricultural Nutrient Inputs to Rivers and Groundwaters in the UK: Policy, Environmental Management and Research Needs</article-title>. <source>Sci. Total Environ.</source> <volume>282</volume>, <fpage>9</fpage>&#x2013;<lpage>24</lpage>. <pub-id pub-id-type="doi">10.1016/s0048-9697(01)00935-4</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Worrall</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Jarvie</surname>
<given-names>H. P.</given-names>
</name>
<name>
<surname>Howden</surname>
<given-names>N. J.</given-names>
</name>
<name>
<surname>Burt</surname>
<given-names>T. P.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The Fluvial Flux of Total Reactive and Total Phosphorus from the UK in the Context of a National Phosphorus Budget: Comparing UK River Fluxes with Phosphorus Trade Imports and Exports</article-title>. <source>Biogeochemistry</source> <volume>whit130</volume> (<issue>1</issue>), <fpage>31</fpage>&#x2013;<lpage>51</lpage>. <pub-id pub-id-type="doi">10.1007/s10533-016-0238-0</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>