M. L. Rodríguez-Blanco
M. M. Taboada-Castro2
M. T. Taboada-Castro- 1Physical Geography Area, History, Art and Geography Department, GEAAT Group, University of Vigo, Ourense, Spain
- 2Area of Soil Science and Soil Chemistry, Higher Technical School of Agricultural Engineering, University of Valladolid, Palencia, Spain
- 3Faculty of Sciences, University of A Coruna, A Coruña, Spain
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.
Introduction
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 (Bennet et al., 2001; Worrall et al., 2016; Grizzetti et al., 2021; Kelly et al., 2022), and this explains why in recent years, determined efforts have been made both in the field of research and management and restoration.
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; Directive, 2000) 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 (European Environment Agency, 2018) 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 (European Environment Agency, 2018). 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 (Oenema et al., 2005; Sharpley et al., 2008; Lassaletta et al., 2009). 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 (Panagos et al., 2015). 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 (Sharpley et al., 2008; Dupas et al., 2017).
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 (Carvahlo et al., 2019). 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 (Fovet et al., 2015; Ehrbardt et al., 2021). 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 (Dupas et al., 2019). 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; Minaudo et al., 2015; Moatar et al., 2017). 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 (Withers and Lord, 2002; Ekholm et al., 2015).
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 “the country of a thousand rivers.” 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.
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 years in which regular manual sampling was combined with automatic sampling during runoff episodes.
Material and Methods
Study Area
The study site is the Corbeira stream headwater catchment located in NW Spain (Figure 1, gauge at 43° 12′58″N, 8° 13′48″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ña and surroundings (about 450,000 inhabitants) discharging into the Ria do Burgo, in the vicinity of the city of A Coruña (Figure 1). The total length of the stream is 10 km, and the whole catchment area covers about 16 km2, with altitudes ranging from 65 m a.s.l. at the confluence with the Mero River to 470 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 (IGME, 1981), 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 (IUSS Working Group WRB, 2014) 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% (Rodríguez-Blanco et al., 2015). The P Olsen content varied between 14.0 mg kg−1 in the forest soils and 69.0 mg kg−1 in cultivated fields, with the mean content of grassland 19.5 mg kg−1 (Rodríguez-Blanco et al., 2013).
FIGURE 1. Location of the study area. Changes in land use observed in the study area: (A) 2000 and (B) 2018 based on SIOSE.
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 Alnus glutinosa, Betula alba and Salix spp. 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’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.
The distribution of land uses in the catchment is characterized by the mosaic presented in agricultural areas (Figure 2), 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 (Dou et al., 2021).
FIGURE 2. Illustrative photographs of the study area. (A) Landscape mosaic in which meadows, cultivated land and forest area; (B) runoff generated on cultivated land devoid of vegetation and channeled through the ditch; (C) slurry application on the ground at field capacity, and (D) riverside vegetation.
The study area has a marked rural nature (<150 inhab. km−2; Organization for Economic Co-operation and Development- OECD, 2007), characterized by a low population density, which has gone from 35 inhab. km−2 at the beginning of the study (2004) to 32 inhab. km-2 in 2018 (rates calculated from data of Spanish Statistical Office-INE, 2019)- a 15% decrease in just 15 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–2018, the population in Galician villages fell by 17.9% (INE, 2019). 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.
The climate is temperate oceanic, with an annual average temperature of 13.5°C ranging from 8.9°C in January and February (winter) to 18.8°C in July and August (summer). The rainfall is abundant (1,074 mm- period 1983/1984–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) (Figure 3; Rodríguez-Blanco et al., 2020). Mean stream discharge is 0.181 m3s−1, ranging from 0.08 to 0.273 m3s−1 (Rodríguez-Blanco et al., 2020). Baseflow contribution to the total streamflow (i.e., baseflow index), calculated according to Gustard et al. (1992), 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.
FIGURE 3. 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 Rodríguez-Blanco et al., 2020). A: autumn, W: winter, Sp: spring, S: summer.
Field Instrumentation and Water Sampling
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 mm resolution), followed by the Thiessen Polygon method (Linsely et al., 1949) 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 min intervals were aggregated to 10-min values. Water sampling was carried out under baseflow conditions (manual sampling every 2 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 Rodríguez-Blanco et al. (2013, 2015).
Laboratory Analyses
Water samples were immediately brought back to the laboratory and refrigerated (4°C) until analysed for nitrate and total and dissolved phosphorus concentrations. After sample filtration (0.45 μm), nitrate concentrations (expressed as mg NO3− L−1) were determined by capillary electrophoresis. Total phosphorus (TP) (μg P L−1) was determined in unfiltered water samples by acid digestion with ammonium persulfate (APHA, 1998) followed by colorimetric determination using the molybdate-ascorbic acid method developed by (Murphy and Riley, 1962). Following water sample filtration, dissolved phosphorus (DP) was determined by ICP-MS. Particulate phosphorus (PP) was calculated as the difference between TP and DP.
Data and Statistical Analysis
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).
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 −1 and 1 mg L−1 indicates the absence of trend, whereas those higher than 1 and 5 mg L−1 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 (Bouraoui and Malagó, 2020) and have also been recommended by the World Meteorological Organisation for analysing the trends in hydrometeorological time series (Kumar et al., 2009). 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.
Results
Exploratory Analysis: Annual and Seasonal Nitrate and Phosphorus Loads and Concentrations
The average nitrate and phosphorus export during the study period was 19.04 kg ha−1 year−1 (7.36 and 30.07 kg ha−1 year−1) and 0.13 kg ha−1 year−1 (0.08 and 0.22 kg ha−1 year−1) for nitrate and phosphorus respectively (Figure 4). 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 (Figure 4), 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 (Bol et al., 2018). 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 (Rodríguez-Blanco et al., 2015).
FIGURE 4. Temporal evolution of annual (A) rainfall, (B) discharge, (C) NO3− load, P load (D), mean (E) NO3− concentration and (F) P mean concentrations during the study period. TP: total phosphorus, PP: particulate phosphorus, DP: dissolved phosphorus.
Mean annual concentrations ranged from 4.1 to 5.2 mg L−1 for NO3−, and 17.2 and 35.1 μg L−1 for TP (Figure 4), reaching mean values of 4.8 mg L−1 and 26.3 μg L−1, respectively, which signifies a greater variability of phosphorus (CV = 20%) concentrations than nitrate (CV = 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 (Figure 4).
Transport throughout the year (Figure 5) showed that the highest nitrate loads were obtained in winter (8.45 kg ha−1 year−1), followed by spring (4.78 kg ha−1) and autumn (4.54 kg ha−1), with the lowest values recorded in summer (1.25 kg ha−1). The pattern for phosphorus was similar, which also returned the highest records in winter (0.050 kg ha−1), followed by autumn (0.030 kg ha−1) and spring (0.031 kg ha−1), with the lowest values observed in summer (0.006 kg ha−1) (Figure 5). The seasonal distribution of the concentrations (Figure 6) followed the same pattern as the loads, meaning that nitrate showed the following decreasing order: winter (5.5 mg L−1) > autumn (5.3 mg L−1) > spring (4.5 mg L−1) > summer (4.1 mg L−1), while phosphorus displayed the sequence: autumn (TP: 32.75 μg L−1; PP: 22.91 μg L−1; DP: 10 μg L−1) > winter (TP: 33 μg L−1; PP: 23 μg L−1; DP: 10 μg L−1) > spring (TP: 24 μg L−1; PP: 15 μg L−1; DP: 8 μg L−1) > summer (TP: 16 μg L−1; PP: 9 μg L−1; DP: 7 μg L−1). Phosphorus concentrations showed large differences between autumn and winter against summer.
FIGURE 5. Box plot with the load of (A) nitrate (NO3−), (B) total phosphorus (TP), (C) particulate phosphorus (PP) and (D) 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.
FIGURE 6. Box plot with the concentration of (A) nitrate (NO3−), (B) total phosphorus (TP), (C) particulate phosphorus (PP) and (D) 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.
Trends in Nitrate and Phosphorus Loads and Concentrations
Tables 1–3 present the results obtained with the Mann-Kendall test. Statistically relevant trends in annual values were not observed for both nutrients (Table 1). Nitrate seems to show a limited downward prospect. The same results were obtained using the Nitrate Directive approach, because the difference between the NO3− flow-weighted concentrations of two consecutive reporting periods ranged between −1 and 1 mg L−1, indicating the absence of a trend. On a seasonal scale, a meaningful decrease is only observed in autumn for the nitrate load (Table 2), mainly linked to the decrease in October, as well as an increase in TP concentrations in winter (Table 3). A negative trend is also seen in the load of NO3− 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.
TABLE 1. Mann-Kendall trend test results for nutrient yield and concentrations.
TABLE 2. Mann-Kendall trend test results for seasonal and monthly nutrient loads.
TABLE 3. Mann-Kendall trend test results for seasonal and monthly nutrient concentrations.
Discussion
NO3− and P loads were similar or lower compared with other studies carried out in Galicia (Bernárdez et al., 2013; Serrano et al., 2015) 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 (Fovet et al., 2015; Dupas et al., 2017; Bouraoui and Malago, 2020; Ehrhardt et al., 2021). In fact, annual mean concentrations of both elements were well below the recommendations of WFD, the threshold proposed by Camargo and Alonso (2006) for protecting the most sensitive freshwater species and the threshold stablished by Carleton et al. (2009) to prevent nuisance level of periphyton and cyanobacteria (Table 4). When the mean nitrate and phosphorus concentrations obtained in this study are compared to the threshold values reported by Poikane et al. (2021), 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 (Table 4). They also are in the boundary values proposed by de Vries et al. (2013) 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 ± 7 μg L−1; 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.
TABLE 4. Nutrient thresholds establish by different ecology-based strategies.
Nitrate concentration exceeds the concentration collected by Macías et al. (1991) as reference values in slightly polluted rivers (1.19 mg L−1), the background concentration proposed by Meybeck (1982) for the major unpolluted rivers (0.44 mg L−1), as well as the levels suggested for uncontaminated natural rivers in Galicia (1–4 mg L−1, Antelo and Arce, 1996) and the mean concentration in the catchment reference site (1.99 ± 0.54 mg L−1; mean and standard deviation from 120 water sampled collected from 2004 to 2016, unpublished data). Naves et al. (2019) observed nitrate contamination in the groundwaters (concentrations between 2.54 and 47.1 mg L−1) 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 = 4.13 mg L−1) -when discharge is dominated by groundwater-are well below the values reported by Naves et al. (2019) as a reference in the study area (1.2–9.9 mg L−1).
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 (Rodríguez-Blanco et al., 2019), 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 (Abbott et al., 2018; Dupas et al., 2019; Garzon-Vidueira et al., 2020). Thus, wet periods increase the connectivity between nutrient sources and the stream network (Rodríguez-Blanco et al., 2013; 2015), 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 (Rodríguez-Blanco et al., 2015) while P (mainly associated with sediment) in surface runoff (Rodríguez-Blanco et al., 2013).
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 (Table 2), 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—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 (Rodríguez-Blanco et al., 2020). 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 (Minaudo et al., 2015; Bol et al., 2018), 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., Abbott et al., 2018), this could be indicative of surplus nutrients applied to the soil being available for export to water bodies by drainage or runoff.
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 (Fovet et al., 2015; Dupas et al., 2017). 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.
Concluding Remarks
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 “natural nature,” 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.
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 via drainage or runoff. As a result, a decrease in inputs from various sources does not result in rapid improvement in riverine water quality.
Data Availability Statement
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 maria.luz.rodriguez.blanco@uvigo.es.
Author Contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
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.
Acknowledgments
We would like to thank two reviewers for their helpful comments and suggestions.
Conflict of Interest
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.
References
Abbott, B. W., Moatar, F., Gauthier, O., Fovet, O., Antoine, V., and Ragueneau, O. (2018). Trends and Seasonality of River Nutrients in Agricultural Catchments: 18 Years of Weekly Citizen Science in France. Sci. Total Environ. 624, 845–858. doi:10.1016/j.scitotenv.2017.12.176
Antelo, J. M., and Arce, F. (1996). “Características fisicoquímicas das augas superficiais,” in As Aguas de Galicia. Editor F. Díaz-Fierros (Santiago de Compostela, Spain: Consello da Cultura Galega), 353–446.
APHA (1998). Standard Methods for the Examination of Water and Wastewater 20th Edition. Washington, DC: American Public Health Association, American Water Work Association, Water Environment Federation.
Bennett, E. M., Carpenter, S. R., and Caraco, N. F. (2001). Human Impact on Erodable Phosphorus and Eutrophication: A Global Perspective. Bioscience 51 (3), 227–234. doi:10.1641/0006-3568(2001)051[0227:hioepa]2.0.co;2
Bernárdez, P., Ospina-Álvarez, N., Caetano, M., and Prego, R. (2013). Fluvial Contributions of Nutrient Salts, Dissolved Trace Elements and Organic Carbon to the Sea by Pristine Temperate Rivers (SW Europe). Environ. Chem. 10 (1), 42–53. doi:10.1071/en12123
Bol, R., Gruau, G., Mellander, P., Dupas, R., Bechmann, M., Skarbøvik, E., et al. (2018). Challenges of Reducing Phosphorus Based Water Eutrophication in the Agricultural Landscapes of Northwest Europe. Front. Mar. Sci. 5, 276. doi:10.3389/fmars.2018.00276
Bouraoui, F., and Malagó, A. (2020). Trend Analysis of Nitrate Concentration in Rivers in Southern France. Water 12 (12), 3374. doi:10.3390/w12123374
Camargo, J. A., and Alonso, Á. (2006). Ecological and Toxicological Effects of Inorganic Nitrogen Pollution in Aquatic Ecosystems: a Global Assessment. Environ. Int. 32 (6), 831–849. doi:10.1016/j.envint.2006.05.002
Carleton, J. N., Park, R. A., and Clough, J. S. (2009). Ecosystem Modeling Applied to Nutrient Criteria Development in Rivers. Environ. Manage 44 (3), 485–492. doi:10.1007/s00267-009-9344-2
Carvalho, L., Mackay, E. B., Cardoso, A. C., Baattrup-Pedersen, A., Birk, S., Blackstock, K. L., et al. (2019). Protecting and Restoring Europe's Waters: An Analysis of the Future Development Needs of the Water Framework Directive. Sci. Total Environ. 658, 1228–1238. doi:10.1016/j.scitotenv.2018.12.255
de Vries, W., Kros, J., Kroeze, C., and Seitzinger, S. P. (2013). Assessing Planetary and Regional Nitrogen Boundaries Related to Food Security and Adverse Environmental Impacts. Curr. Opin. Environ. Sustain. 5 (3-4), 392–402. doi:10.1016/j.cosust.2013.07.004
Dou, Y., Cosentino, F., Malek, Z., Maiorano, L., Thuiller, W., and Verburg, P. H. (2021). A New European Land Systems Representation Accounting for Landscape Characteristics. Landsc. Ecol. 36 (8), 2215–2234. doi:10.1007/s10980-021-01227-5
Dupas, R., Abbott, B. W., Minaudo, C., and Fovet, O. (2019). Distribution of Landscape Units within Catchments Influences Nutrient Export Dynamics. Front. Environ. Sci. 7, 43. doi:10.3389/fenvs.2019.00043
Dupas, R., Mellander, P., Gascuel-Odoux, C., Fovet, O., McAleer, E. B., McDonald, N. T., et al. (2017). The Role of Mobilisation and Delivery Processes on Contrasting Dissolved Nitrogen and Phosphorus Exports in Groundwater Fed Catchments. Sci. Total Environ. 599, 1275–1287. doi:10.1016/j.scitotenv.2017.05.091
Ehrhardt, S., Ebeling, P., Dupas, R., Kumar, R., Fleckenstein, J. H., and Musolff, A. (2021). Nitrate Transport and Retention in Western European Catchments Are Shaped by Hydroclimate and Subsurface Properties. Water Resour. Res. 57 (10), e2020WR029469. doi:10.1029/2020wr029469
Ekholm, P., Rankinen, K., Rita, H., Räike, A., Sjöblom, H., Raateland, A., et al. (2015). Phosphorus and Nitrogen Fluxes Carried by 21 Finnish Agricultural Rivers in 1985–2006. Environ. Monit. Assess. 187 (4), 216–217. doi:10.1007/s10661-015-4417-6
European Environment Agency (2018). European Waters Assessment. EEA Report N° 7/2018. Copenhage, Denmark: EEA.
Fovet, O., Ruiz, L., Faucheux, M., Molénat, J., Sekhar, M., Vertès, F., et al. (2015). Using Long Time Series of Agricultural-Derived Nitrates for Estimating Catchment Transit Times. J. Hydrology 522, 603–617. doi:10.1016/j.jhydrol.2015.01.030
Garzon-Vidueira, R., Rial-Otero, R., Garcia-Nocelo, M. L., Rivas-Gonzalez, E., Moure-Gonzalez, D., Fompedriña-Roca, D., et al. (2020). Identification of Nitrates Origin in Limia River Basin and Pollution-Determinant Factors. Agric. Ecosyst. Environ. 290, 106775.
Grizzetti, B., Vigiak, O., Udias, A., Aloe, A., Zanni, M., Bouraoui, F., et al. (2021). How EU Policies Could Reduce Nutrient Pollution in European Inland and Coastal Waters. Glob. Environ. Change 69, 102281. doi:10.1016/j.gloenvcha.2021.102281
Gustard, A., Bullock, A., and Dixon, J. M. (1992). Low Flow Estimation in the United Kingdom. England: Institute of Hydrology, 19.
IGME (Instituto Tecnológico Geominero de España) (1981). Mapa Geológico de España, 1:50,000. Betanzos Spain. Hoja 45.
INE (2019). Nomenclator: Población del padrón continuo por unidad poblacional. Madrid: Instituto Estadístico de España.
Iuss Working Group Wrb, (2014). in International Soil Clasification System for Naming Soils and Creating Legends for Soil Maps (Rome: World Soil Resources Reports. N° 106FAO).World Reference Base for Soil Resources.
Kelly, M. G., Phillips, G., Teixeira, H., Várbíró, G., Herrero, F. S., Willby, N. J., et al. (2022). Establishing Ecologically-Relevant Nutrient Thresholds: A Tool-Kit with Guidance on its Use. Sci. Total Environ. 807, 150977. doi:10.1016/j.scitotenv.2021.150977
Kumar, S., Merwade, V., Kam, J., and Thurner, K. (2009). Streamflow Trends in Indiana: Effects of Long Term Persistence, Precipitation and Subsurface Drains. J. Hydroloy 374 (1-2), 171–183. doi:10.1016/j.jhydrol.2009.06.012
Lassaletta, L., García-Gómez, H., Gimeno, B. S., and Rovira, J. V. (2009). Agriculture-induced Increase in Nitrate Concentrations in Stream Waters of a Large Mediterranean Catchment over 25 Years (1981–2005). Sci. Total Environ. 407 (23), 6034–6043. doi:10.1016/j.scitotenv.2009.08.002
Linsely, R. K., Kohler, M. A., and Paulhus, J. C. (1949). Applied Hydrology. New York: McGraw-Hill Book Co.
Macías, F., Álvarez Rodríguez, E., and Calvo de Anta, R. (1991). Impactos de origen agrario y urbano en la cuenca del río Ulla. Ecología 5, 73–86.
Meybeck, M. (1982). Carbon, Nitrogen, and Phosphorus Transport by World Rivers. Am. J. Sci. 282 (4), 401–450. doi:10.2475/ajs.282.4.401
Minaudo, C., Meybeck, M., Moatar, F., Gassama, N., and Curie, F. (2015). Eutrophication Mitigation in Rivers: 30 Years of Trends in Spatial and Seasonal Patterns of Biogeochemistry of the Loire River (1980–2012). Biogeosciences 12 (8), 2549–2563. doi:10.5194/bg-12-2549-2015
Moatar, F., Abbott, B. W., Minaudo, C., Curie, F., and Pinay, G. (2017). Elemental Properties, Hydrology, and Biology Interact to Shape Concentration-discharge Curves for Carbon, Nutrients, Sediment, and Major Ions. Water Resour. Res. 53 (2), 1270–1287. doi:10.1002/2016wr019635
Murphy, J., and Riley, J. P. (1962). A Modified Single Solution Method for the Determination of Phosphate in Natural Waters. Anal. Chim. Acta 27, 31–36. doi:10.1016/s0003-2670(00)88444-5
Naves, A., Samper, J., Mon, A., Pisani, B., Montenegro, L., and Carvalho, J. M. (2019). Demonstrative Actions of Spring Restoration and Groundwater Protection in Rural Areas of Abegondo (Galicia, Spain). Sustain. Water Resour. Manag. 5 (1), 175–186. doi:10.1007/s40899-017-0169-5
Oenema, O., van Liere, L., and Schoumans, O. (2005). Effects of Lowering Nitrogen and Phosphorus Surpluses in Agriculture on the Quality of Groundwater and Surface Water in the Netherlands. J. Hydrology 304 (1-4), 289–301. doi:10.1016/j.jhydrol.2004.07.044
Panagos, P., Borrelli, P., Poesen, J., Ballabio, C., Lugato, E., Meusburger, K., et al. (2015). The New Assessment of Soil Loss by Water Erosion in Europe. Environ Sci Policy 54, 438–447. doi:10.1016/j.envsci.2015.08.012
Poikane, S., Várbíró, G., Kelly, M. G., Birk, S., and Phillips, G. (2021). Estimating River Nutrient Concentrations Consistent with Good Ecological Condition: More Stringent Nutrient Thresholds Needed. Ecol. Ind. 121, 107017. doi:10.1016/j.ecolind.2020.107017
Rodríguez-Blanco, M. L., Taboada-Castro, M. M., and Taboada-Castro, M. T. (2020). 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). Hydrology 7 (4), 97. doi:10.3390/hydrology7040097
Rodríguez-Blanco, M. L., Taboada-Castro, M. M., and Taboada-Castro, M. T. (2019). An Overview of Patterns and Dynamics of Suspended Sediment Transport in an Agroforest Headwater System in Humid Climate: Results from a Long-Term Monitoring. Sci. Total Environ. 648, 33–43. doi:10.1016/j.scitotenv.2018.08.118
Rodríguez-Blanco, M. L., Taboada-Castro, M. M., and Taboada-Castro, M. T. (2013). Phosphorus Transport into a Stream Draining from a Mixed Land Use Catchment in Galicia (NW Spain): Significance of Runoff Events. J. Hydrology 481, 12–21. doi:10.1016/j.jhydrol.2012.11.046
Rodríguez-Blanco, M. L., Taboada-Castro, M. M., Taboada-Castro, M. T., and Oropeza-Mota, J. L. (2015). Relating Nitrogen Export Patterns from a Mixed Land Use Catchment in NW Spain with Rainfall and Streamflow. Hydrol. Process 29 (12), 2720–2730. doi:10.1002/hyp.10388
Serrano, L., Palleiro, L., Rodríguez-Blanco, M. L., Taboada-Castro, M. M., and Taboada-Castro, M. T. (2015). Preliminary Assessment of the Riverine Nitrogen Concentration and Load in an Agroforestry Basin in NW Spain. Commun. Soil Sci. Plant Anal. 46, 326–331. doi:10.1080/00103624.2014.989071
Sharpley, A. N., Kleinman, P. J., Heathwaite, A. L., Gburek, W. J., Folmar, G. J., and Schmidt, J. P. (2008). Phosphorus Loss from an Agricultural Watershed as a Function of Storm Size. J. Environ. Qual. 37 (2), 362–368. doi:10.2134/jeq2007.0366
Withers, P. J., and Lord, E. I. (2002). Agricultural Nutrient Inputs to Rivers and Groundwaters in the UK: Policy, Environmental Management and Research Needs. Sci. Total Environ. 282, 9–24. doi:10.1016/s0048-9697(01)00935-4
Worrall, F., Jarvie, H. P., Howden, N. J., and Burt, T. P. (2016). 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. Biogeochemistry whit130 (1), 31–51. doi:10.1007/s10533-016-0238-0
Keywords: long-term, nutrients, Galicia, rural catchment, land use
Citation: Rodríguez-Blanco ML, Taboada-Castro MM and Taboada-Castro MT (2023) Nitrate and Phosphorus Transport in a Galician River (NW Iberian Peninsula): Insights From Fourteen Years of Monitoring. Span. J. Soil Sci. 13:11108. doi: 10.3389/sjss.2023.11108
Received: 06 December 2022; Accepted: 05 January 2023;
Published: 13 January 2023.
Edited by:
Avelino Núñez-Delgado, University of Santiago de Compostela, SpainCopyright © 2023 Rodríguez-Blanco, Taboada-Castro and Taboada-Castro. 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.
*Correspondence: M. L. Rodríguez-Blanco, maria.luz.rodriguez.blanco@uvigo.es