Long-Term Concentrations and Loads of Four Dissolved Macronutrients from Two Agroforestry Catchments in NW Spain

: Understanding hydrological processes controlling stream chemistry and quantifying solute concentrations over time is crucial for estimating future alterations of water quality due to land use or climate change impacts, as well as for setting preventive or remedial actions. In the current study, soluble sodium (Na + ), potassium (K + ), calcium (Ca 2+ ), and magnesium (Mg 2+ ) concentrations and loads were assessed in waters from two small catchments located at different distances to the sea (9 km Valiñas and 30 km Abelar) in NW Spain from 2003 to 2016. Solute concentrations were determined using spectrometric techniques, while streamﬂow data were employed for estimating loads. Moreover, concentration-discharge relations were calculated for each solute and catchment. The average concentrations of soluble Na + , K + , Ca 2+ and Mg 2+ were, respectively, 16.5, 2.6, 7.9 and 4.1 mg L − 1 in Valiñas, and 8.2, 0.9, 2.9 and 3.4 mg L − 1 in Abelar, although variability among samplings was high. The four soluble ions showed a dilution pattern in Valiñas, whereas in Abelar Na + tended to a chemostatic behavior and K + and Ca 2+ were positively related to streamﬂow. In conclusion, the dominant processes controlling these relationships are local and depend on catchment characteristics such as land use (including slurry applications in Abelar), distance to the sea, and vegetation cover.


Introduction
Peri-urban catchments are characterized by a wide range of land uses, including agriculture, forestry and recreation, that impact on water quality [1,2]. Moreover, solute concentrations in water depend on catchment characteristics such as climatology, topography, land use, size and bedrock nature [3,4]. In this context, understanding the hydrological processes that control stream chemistry, as well as the sources of chemical inputs, is highly relevant for estimating how changing climate and future land use may affect the chemistry of water bodies [5]. In fact, the material transported with water provides information about residence time, biogeochemical transformations and water sources [6]. Moreover, quantifying solute concentrations over time can allow for the assessment of the anthropogenic impacts on catchments, and help set preventive or remedial actions [7].
Rock weathering and solute transport are linked to hydrology in watersheds and this coupling is reflected in the relationships between stream discharge and solute concentrations [8]. Therefore, the study of concentration-discharge (C-Q) relationships allows for a better understanding of hydrological processes and chemical sources [5,[9][10][11]. The C-Q relations are usually linear in a log(C)-log(Q) space [8], thus the empirical relations can be approximated by a power law, C = a × Q b , where a and b are fitting parameters [4,8,11]. In The Valiñas River is a small headwater stream with a catchment that is 36.3 km 2 in surface. Topographically, this catchment consists of steep uplands and contrasting lowlands or local small plains. The bedrock consists of granite (85%) and some outcrops of schists (15%) from the Ordes series in the Northeast of the catchment [20]. The soils within this catchment present loam and sandy-loam textures, with pH ranging from 4.1 to 6.1 [20]. This catchment is located at 9 km from the coast, approximately, and it has been transformed by human activities. The main land use is agriculture (55% of the catchment is cultivated) but forest is also present (45% of the surface). The steep upland slopes are well vegetated with woodland and shrubs. Agriculture has dominated the lowlands and remains relevant nowadays. In fact, the lower part of the Valiñas catchment is dominated by pastures, old orchards and scattered fields in patches where maize is grown. However, an increasing area is being occupied by housing and light industry. Small villages and buildings are not connected to a public sewer system and release sewage into the river, which negatively affects water quality. Surface runoff from village streets, roads and a highway also impact the water quality of this stream [20]. The climate is humid and temperate with an annual total rainfall amount of 1102 mm and an average temperature of 13.8 °C from 1981 to 2010. Figure 2 shows the monthly rainfall amounts for a 30-year period (1981-2010) and those from the studied years (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016). These data were recorded at the weather station of Alvedro (A Coruña), located 3 km from the catchment outlet, approximately. It can be assumed that 2004, 2007, 2011 and 2015 were extraordinarily dry because annual rainfall was, respectively, 29%, 31%, 34% and 32% lower than the 30-year average ( Figure 2). However, rainfall amounts in 2003,2006,2008,2009,2013 and 2016 were close to the average, representing the rainy conditions of the area. The monthly rainfall amounts showed the high temporal variability of rainfall in this region [21]. The Valiñas River is a small headwater stream with a catchment that is 36.3 km 2 in surface. Topographically, this catchment consists of steep uplands and contrasting lowlands or local small plains. The bedrock consists of granite (85%) and some outcrops of schists (15%) from the Ordes series in the Northeast of the catchment [20]. The soils within this catchment present loam and sandy-loam textures, with pH ranging from 4.1 to 6.1 [20]. This catchment is located at 9 km from the coast, approximately, and it has been transformed by human activities. The main land use is agriculture (55% of the catchment is cultivated) but forest is also present (45% of the surface). The steep upland slopes are well vegetated with woodland and shrubs. Agriculture has dominated the lowlands and remains relevant nowadays. In fact, the lower part of the Valiñas catchment is dominated by pastures, old orchards and scattered fields in patches where maize is grown. However, an increasing area is being occupied by housing and light industry. Small villages and buildings are not connected to a public sewer system and release sewage into the river, which negatively affects water quality. Surface runoff from village streets, roads and a highway also impact the water quality of this stream [20]. The climate is humid and temperate with an annual total rainfall amount of 1102 mm and an average temperature of 13.8 • C from 1981 to 2010. Figure 2 shows the monthly rainfall amounts for a 30-year period (1981-2010) and those from the studied years (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016). These data were recorded at the weather station of Alvedro (A Coruña), located 3 km from the catchment outlet, approximately. It can be assumed that 2004, 2007, 2011 and 2015 were extraordinarily dry because annual rainfall was, respectively, 29%, 31%, 34% and 32% lower than the 30-year average ( Figure 2). However, rainfall amounts in 2003,2006,2008,2009,2013 and 2016 were close to the average, representing the rainy conditions of the area. The monthly rainfall amounts showed the high temporal variability of rainfall in this region [21]. The Abelar catchment covers 10.7 ha (0.107 km 2 ) and is located at 30 km from the coast, approximately. Geologically, the catchment is characterized by a rocky substrate, dominated by metamorphic schists from the Ordes series, in which permeability is associated with fissuring and surface alteration. The low permeability of these rock formations means that the piezometric level tends to follow the shape of the topography [22]. The soil is loam to silt-loam, with depths ranging from 45 to 125 cm and an acid pH, around 5.1 [23]. An operational farm has been established in this catchment for several years and excessive doses of slurry were applied to the soils. Later, the farm was replaced by a Eucalyptus globulus stand (1270 trees ha −1 ) covering most of the catchment area in 1998, approximately 85%. Two narrow power-line corridors with maize and pasture (11% of the area) and a small riparian woodland (4% of the area) still remain [23]. No applications of slurry were performed over the Eucalyptus plantation. However, a small maize field, approximately 3% of the catchment surface, remained within the studied site and, occasionally, received applications of slurry, namely twice between 1998 and 2010, affecting water flow and element concentration at the catchment outlet. A comprehensive description of this site can be found in Rodríguez-Suárez et al. [22]. The climate is humid and temperate with an annual total rainfall amount of 1787 mm and an average temperature of 13.0 °C, for the 1981-2010 period. Figure 3 shows the monthly rainfall amounts for a 30-year period (1981-2010) and those from the studied years (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016). These data were recorded at a weather station located within the study site, whereas 30-year data came from the meteorological station of Lavacolla (Santiago de Compostela, A Coruña province), located 10 km from the catchment outlet, approximately. It can be assumed that 2007 was dry, while 2011 and 2013 were extraordinarily dry because annual rainfall was, respectively, 48%, 37% and 33% lower than the 30-year average ( Figure 3). However, rainfall amounts in 2008, 2010, 2014 and 2016 were close to the average, representing the rainy conditions of the area. The Abelar catchment covers 10.7 ha (0.107 km 2 ) and is located at 30 km from the coast, approximately. Geologically, the catchment is characterized by a rocky substrate, dominated by metamorphic schists from the Ordes series, in which permeability is associated with fissuring and surface alteration. The low permeability of these rock formations means that the piezometric level tends to follow the shape of the topography [22]. The soil is loam to silt-loam, with depths ranging from 45 to 125 cm and an acid pH, around 5.1 [23]. An operational farm has been established in this catchment for several years and excessive doses of slurry were applied to the soils. Later, the farm was replaced by a Eucalyptus globulus stand (1270 trees ha −1 ) covering most of the catchment area in 1998, approximately 85%. Two narrow power-line corridors with maize and pasture (11% of the area) and a small riparian woodland (4% of the area) still remain [23]. No applications of slurry were performed over the Eucalyptus plantation. However, a small maize field, approximately 3% of the catchment surface, remained within the studied site and, occasionally, received applications of slurry, namely twice between 1998 and 2010, affecting water flow and element concentration at the catchment outlet. A comprehensive description of this site can be found in Rodríguez-Suárez et al. [22]. The climate is humid and temperate with an annual total rainfall amount of 1787 mm and an average temperature of 13.0 • C, for the 1981-2010 period. Figure 3 shows the monthly rainfall amounts for a 30-year period (1981-2010) and those from the studied years (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016). These data were recorded at a weather station located within the study site, whereas 30-year data came from the meteorological station of Lavacolla (Santiago de Compostela, A Coruña province), located 10 km from the catchment outlet, approximately. It can be assumed that 2007 was dry, while 2011 and 2013 were extraordinarily dry because annual rainfall was, respectively, 48%, 37% and 33% lower than the 30-year average ( Figure 3). However, rainfall amounts in 2008, 2010, 2014 and 2016 were close to the average, representing the rainy conditions of the area.

Sampling and Laboratory Determinations
The study period extended from January 2003 to December 2016 for Valiñas, and from January 2007 to December 2016 for Abelar. Water samples were collected at the catchment outlets under base and storm flow conditions. The sampling strategy was a stratified point sampling involving more frequent collections when flow was high [24] since a continuous sampling, which would improve the characterization of the rising and falling of the hydrograph [25], was not possible. A total of 984 samples were taken at the outlet of the Valiñas catchment (70 ± 13 samples per year), whereas 391 samples were collected at the outlet of Abelar (39 ± 8 samples per year). The number of samples taken each individual year was a consequence of the frequency and intensity of rainfall.
Samples were collected using polyethylene terephthalate (PET) bottles and they were stored at 4 °C for preservation purposes. The concentrations of four macronutrients, namely sodium (Na + ), potassium (K + ), calcium (Ca 2+ ) and magnesium (Mg 2+ ), were determined directly on the samples after filtration through a 0.45 μm filter [26]. These concentrations were measured by inductively coupled argon plasma-mass spectrometry (ICP-MS) using a Jobin Yvon JY 50-P instrument (Jobin Yvon, Longjumeau, France).

Mass Balance Calculations
In order to estimate the mass balance of Na + , K + , Ca 2+ and Mg 2+ , discharge data (Q) at the catchment outlets are required. In the current study, Q data were not continuously recorded. Therefore, daily Q was estimated using Visual Balan v.1.0 [27]. This software computes the water balance at a daily timestep. For streamflow calculation, the model uses the Muskingum-Cunge method [28]. The inputs for the water balance considered by Visual Balan are the contributions from tributary streams and rainfall, whereas, the outputs considered are runoff, infiltration, evaporation and surface discharge. The only data that this model requires are the catchment characteristics and daily records of rainfall and reference evapotranspiration. In the current study, the weather stations mentioned when describing the two catchments considered provided the daily records of rainfall and reference evapotranspiration. Estimations from this model show a reasonable fit with direct Q measurements in the studied catchments [29]. Moreover, the results

Sampling and Laboratory Determinations
The study period extended from January 2003 to December 2016 for Valiñas, and from January 2007 to December 2016 for Abelar. Water samples were collected at the catchment outlets under base and storm flow conditions. The sampling strategy was a stratified point sampling involving more frequent collections when flow was high [24] since a continuous sampling, which would improve the characterization of the rising and falling of the hydrograph [25], was not possible. A total of 984 samples were taken at the outlet of the Valiñas catchment (70 ± 13 samples per year), whereas 391 samples were collected at the outlet of Abelar (39 ± 8 samples per year). The number of samples taken each individual year was a consequence of the frequency and intensity of rainfall.
Samples were collected using polyethylene terephthalate (PET) bottles and they were stored at 4 • C for preservation purposes. The concentrations of four macronutrients, namely sodium (Na + ), potassium (K + ), calcium (Ca 2+ ) and magnesium (Mg 2+ ), were determined directly on the samples after filtration through a 0.45 µm filter [26]. These concentrations were measured by inductively coupled argon plasma-mass spectrometry (ICP-MS) using a Jobin Yvon JY 50-P instrument (Jobin Yvon, Longjumeau, France).

Mass Balance Calculations
In order to estimate the mass balance of Na + , K + , Ca 2+ and Mg 2+ , discharge data (Q) at the catchment outlets are required. In the current study, Q data were not continuously recorded. Therefore, daily Q was estimated using Visual Balan v.1.0 [27]. This software computes the water balance at a daily timestep. For streamflow calculation, the model uses the Muskingum-Cunge method [28]. The inputs for the water balance considered by Visual Balan are the contributions from tributary streams and rainfall, whereas, the outputs considered are runoff, infiltration, evaporation and surface discharge. The only data that this model requires are the catchment characteristics and daily records of rainfall and reference evapotranspiration. In the current study, the weather stations mentioned when describing the two catchments considered provided the daily records of rainfall and reference evapotranspiration. Estimations from this model show a reasonable fit with direct Q measurements in the studied catchments [29]. Moreover, the results obtained are in agreement with measured discharges in similar catchments from the same region [30]. In order to build a continuous database for mass balance calculations, daily concentrations of the different solutes were interpolated from measured data and discharge was obtained from model estimates. Daily solute loads from each catchment were estimated using the following equation (adapted from [31]): where E is the solute load (kg ha −1 ), Q is the discharge (L s −1 ), C is the concentration of the solute (kg L −1 ), t is time (seconds) and A is the area of the catchment (ha). Then, monthly seasonal and annual averages were computed. In the case of the season, we consider the period from January to March, as winter, April to June as spring, July to September as summer, and October to December as winter.

Relations between Concentration and Discharge
For each solute and each catchment, the relation between solute concentration and discharge (C = a × Q b ) was assessed in order to evaluate solute behaviors between catchments and whether this behavior can be generalized. These two variables are expected to show a linear relation in a log-log scale, expressed by two regression parameters: a, which is the intercept and has the same units as the concentration, and b, which is the unit-less exponent representing the slope of the interpolating line [11]. This b exponent determines the behavior of the solute [8]. As suggested by Meybeck and Moatar [32], C-Q relations were assessed for high and low flows separately in order to obtain a finer classification of the solute behavior into different categories. The criteria for determining high and low flows was the median of the discharge daily values [4,8].

Statistical Analysis
The main descriptive statistics of the studied datasets were computed, including mean, median, minimum, maximum, skewness, kurtosis and coefficient of variation (CV). These are calculated for each solute in each catchment aiming at highlighting differences between catchments, which are a consequence of natural and anthropogenic factors. Relations among the different solutes and Q were assessed through the Pearson's r correlation coefficient. The Kolmogorov-Smirnoff test was used for comparing, on a yearly basis, the frequency distributions of the different solutes studied. The statistical significance of having a b exponent different from zero was assessed by means of Student's t test at 0.05 level of significance [11].

Time Series Analysis
The monthly averages of the concentrations of each solute in each catchment were subjected to time series analysis in order to assess their trend and seasonality. The autocorrelation function was employed for examining how a value depended on preceding values over a period of time. This function is represented with a correlogram. The slope of this correlogram provides the time of response of the hydrological system [33]. Time series analysis has been widely applied in hydrology to assess the inertia of the system and to make forecasting and generate synthetic series which are inputs for the analysis of complex water resources systems [34,35].
Both statistical and time-series analyses were carried out using the R statistical environment v.4.0.5 [36]. Table 1 shows the statistical summary of the concentrations of the four water-soluble ionic species for the two catchments over the whole study period. The mean concentrations of these solutes were significantly higher in Valiñas than in Abelar. However, in both catchments, median concentrations were relatively close to the mean concentrations for each of the nutrients analyzed. In absolute terms, Na + ranged from 8.90 to 27.20 mg L −1 in Valiñas and from 3.65 to 23.20 mg L −1 in Abelar; K + from 1.13 to 13.49 mg L −1 in Valiñas and from 0.21 to 104.00 mg L −1 in Abelar; Ca 2+ ranged from 4.47 to 20.70 mg L −1 in Valiñas and from 1.34 to 9.40 mg L −1 in Abelar; finally, Mg 2+ ranged from 2.15 to 7.30 mg L −1 in Valiñas and from 1.62 to 7.20 mg L −1 in Abelar. This wide range of concentrations caused skewness and kurtosis values to be far from those corresponding to a normal distribution ( Table 1). The coefficients of variation were relatively low for Na + and Mg 2+ in Valiñas (< 20%), whereas this coefficient was extremely high for K + concentrations in Abelar (Table 1).

Magnitude
Pearson's r correlation coefficients ( Table 2) showed that discharge was negatively and significantly correlated with the concentrations of the studied solutes in the Valiñas catchment. In the case of Abelar, Na + and Ca 2+ concentrations were significantly correlated with discharge; in contrast, K + and Mg 2+ concentrations did not show a significant correlation with discharge ( Table 2). The concentrations of the four macronutrients studied were significantly and positively correlated in both catchments (Table 2). Table 2. Pearson's r correlation coefficients among discharge and solute concentrations for the whole dataset from each of the two studied catchments (Valiñas and Abelar). In the case of Abelar, the correlations are shown for both the entire dataset and by removing a single outlier.

Monthly and Seasonal Evolution of Macronutrient Concentrations
In both sites, the highest discharge values occurred in autumn and winter, between November and February, while the lowest values occurred in summer, between July and September ( Figure 4). A noticeable difference in discharge between the two catchments was observed, with higher values in Valiñas.

Monthly and Seasonal Evolution of Macronutrient Concentrations
In both sites, the highest discharge values occurred in autumn and winter, bet November and February, while the lowest values occurred in summer, between Jul September (Figure 4). A noticeable difference in discharge between the two catchm was observed, with higher values in Valiñas. For the Valiñas catchment, the highest concentrations of Na + , K + , Ca 2+ and tended to occur in late summer or early autumn, between August and October (F 5a). In contrast, no particular patterns were observed for Na + , K + and Mg 2+ in Abela to the presence of outliers in the data series, while Ca 2+ concentrations behaved i opposite way as those from Valiñas (Figure 6a). When removing the extreme ou likely caused by an anthropogenic source, Na + and Mg 2+ tended to be constant alon year in Abelar, while K + and Ca 2+ tended to appear in lower concentrations d summer. For the Valiñas catchment, the highest concentrations of Na + , K + , Ca 2+ and Mg 2+ tended to occur in late summer or early autumn, between August and October (Figure 5a). In contrast, no particular patterns were observed for Na + , K + and Mg 2+ in Abelar due to the presence of outliers in the data series, while Ca 2+ concentrations behaved in the opposite way as those from Valiñas (Figure 6a). When removing the extreme outlier, likely caused by an anthropogenic source, Na + and Mg 2+ tended to be constant along the year in Abelar, while K + and Ca 2+ tended to appear in lower concentrations during summer.
Autocorrelation functions showed a seasonality for the concentrations of all macronutrients in Valiñas (Figure 5b), with four-month cycle fluctuations. This seasonality was also observed for Ca 2+ concentrations in Abelar but not for the rest of the solutes, which tended to behave rather constantly in the case of Na + and Mg 2+ , and rather erratically in the case of K + (Figure 6b    In the Valiñas catchments, the concentrations of the four macronutrients showed a higher magnitude in 2008 and 2012 but presented and marked seasonal tendency ( Figures S1-S4). In contrast, the concentrations of the four studied macronutrients in Abelar tended to decrease in the long term ( Figures S5-S8). Kolmogorov-Smirnoff tests indicated that the concentrations of these solutes differed from year to year, with 2008 and 2009 being the years that showed different concentrations with respect to the other years studied (Tables S1-S3).
Moreover, the concentrations of Na + , K + and Ca 2+ were significantly higher (almost double) in Valiñas than in Abelar on a yearly and seasonal basis. Within a single catchment, the concentrations of Na + , Ca 2+ and Mg 2+ showed significant differences among seasons in Valiñas, namely the concentrations were higher in summer and autumn than in spring and winter, while K + concentrations did not differ among seasons (Table 3). In the case of Abelar, Na + and K + concentrations did not differ among seasons, whereas those of Ca 2+ and Mg 2+ were lower in summer (Table 3).  Figure 7 shows plots of concentration versus discharge for the four studied macronutrients in Valiñas. Decreasing concentrations with discharge were clearly observed for sodium and magnesium, whereas potassium and calcium concentrations tended to reduce with increasing discharge to a lesser extent. The b exponents were −0.127, −0.122, −0.123 and −0.156 for Na + , K + , Ca 2+ and Mg 2+ respectively. Solute concentrations decreased up to, approximately, 2 m 3 s −1 of discharge, after which concentrations tended to level off (Figure 7).

Concentration-Discharge Relations
When plotted on a log-log space (Figure 8), the four solutes showed a dilution behavior with increasing discharge and no significant changes in the slope of the regression lines were observed when considering flows lower or greater than the median, except for the case of K + in which a reduction in the slope value was detected under high flow conditions. When plotted on a log-log space (Figure 8), the four solutes showed a dilution behavior with increasing discharge and no significant changes in the slope of the regression lines were observed when considering flows lower or greater than the median, except for the case of K + in which a reduction in the slope value was detected under high flow conditions. In the case of Abelar (Figure 9), the concentrations of Na + and Mg 2+ remained almost constant independently of discharge, whereas those of K + and Ca 2+ tended to increase with discharge ( Figure 9). When plotted on a log-log space, the concentrations of Na + remained practically stable when flow was lower than the median, whereas they tended to decrease when discharge values were greater than the median (Figure 10a). In contrast, In the case of Abelar (Figure 9), the concentrations of Na + and Mg 2+ remained almost constant independently of discharge, whereas those of K + and Ca 2+ tended to increase with discharge ( Figure 9). When plotted on a log-log space, the concentrations of Na + remained practically stable when flow was lower than the median, whereas they tended to decrease when discharge values were greater than the median (Figure 10a). In contrast, K + concentrations increased with discharge and the slope of this increase was much higher under high flow conditions (Figure 10b). The concentrations of Ca 2+ increased with the discharge and no significant change in the slope of this relation was observed when discharge was greater than the median (Figure 10c). In the case of Mg 2+ , the concentration slightly increased with discharge but decreased under high flow conditions (Figure 10d). Considering the whole dataset, half of the C-Q relationships were negative for low flows, one was flat and three were positive. Under high flow conditions, five out of eight C-Q relationships were negative, two were positive and the remaining one was flat. Moreover, Student's t tests revealed that b exponents were significantly different from 0 for all macronutrients in Valiñas and Abelar (p-values ranging from 3.2 × 10 −8 to 0.025), except for that of Na + in Abelar (p-value = 0.205). Considering the whole dataset, half of the C-Q relationships were negative for low flows, one was flat and three were positive. Under high flow conditions, five out of eight C-Q relationships were negative, two were positive and the remaining one was flat. Moreover, Student's t tests revealed that b exponents were significantly different from 0 for all macronutrients in Valiñas and Abelar (p-values ranging from 3.2 × 10 −8 to 0.025), except for that of Na + in Abelar (p-value = 0.205).

Loads of Macronutrients in Each Catchment
The loads of the four major macronutrients studied showed a great variability at the monthly, seasonal and annual scales. Moreover, loads increased in autumn and winter with respect to spring and summer losses in both study sites (Table 4). In general, macronutrient exports from the Valiñas catchment were higher than those from Abelar. Depending on the solute, season and year, the ratio of the loads from Valiñas to those of Abelar ranged from 0.7 to 2.4 for Na + , from 0.5 to 5.0 for K + , from 0.8 to 3.3 for Ca 2+ , and from 0.4 to 1.4 for Mg 2+ .
In Valiñas, the annual transfers of Na + ranged from 57.17 kg ha −1 in 2015 to 99.28 kg ha −1 in 2006. On average for the study period, Na + loads were 56.10 kg ha −1 during autumn and winter, whereas they were 21.05 kg ha −1 during spring and summer. Moreover, Na + losses during the six driest months were 27.64% of the annual losses, on average for the studied period. The annual loads of K + ranged from 8.

Loads of Macronutrients in Each Catchment
The loads of the four major macronutrients studied showed a great variability at the monthly, seasonal and annual scales. Moreover, loads increased in autumn and winter with respect to spring and summer losses in both study sites (Table 4). In general, macronutrient exports from the Valiñas catchment were higher than those from Abelar. Depending on the solute, season and year, the ratio of the loads from Valiñas to those of Abelar ranged from 0.7 to 2.4 for Na + , from 0.5 to 5.0 for K + , from 0.8 to 3.3 for Ca 2+ , and from 0.4 to 1.4 for Mg 2+ .
In Valiñas, the annual transfers of Na + ranged from 57.17 kg ha −1 in 2015 to 99.28 kg ha −1 in 2006. On average for the study period, Na + loads were 56.10 kg ha −1 during autumn and winter, whereas they were 21.05 kg ha −1 during spring and summer. Moreover, Na + losses during the six driest months were 27.64% of the annual losses, on average for the studied period. The annual loads of K + ranged from 8.  In Abelar, the annual loads of Na + ranged from 39.23 kg ha −1 in 2012 to 89.06 kg ha −1 in 2013. On average for the study period, Na + loads were 49.42 kg ha −1 during autumn and winter, whereas they were 13.84 kg ha −1 during spring and summer. Moreover, Na + losses during the six driest months were 21.88% of the annual losses, on average for the studied period. The annual loads of K + ranged from 3.41 kg ha −1 in 2012 to 20.81 kg ha −1 in 2008. On average for the 14 years of study, K + losses were 7.33 kg ha −1 during autumn and winter, and 1.16 kg ha −1 during spring and summer. The K + transfers during the six driest months were 13.65% of the annual losses. The annual loads of Ca 2+ ranged from 14.82 kg ha −1 in 2012 to 40.30 kg ha −1 in 2010. On average, Ca 2+ loads were 20.90 kg ha −1 during autumn and winter, and 5.07 kg ha −1 during spring and summer. The Ca 2+ transfers during the six driest months were 19.54% of the annual losses. Finally, the annual transfers of Mg 2+ ranged from 17.07 kg ha −1 in 2012 to 39.19 kg ha −1 in 2010. On average, for the study period, Mg 2+ loads were 21.29 kg ha −1 during autumn and winter, and 5.86 kg ha −1 during spring and summer. The Mg 2+ transfers during the six driest months were 21.59% of the annual losses.

Influence of Catchment Features on Macronutrient Concentrations and Loads
In the current study, the mean concentrations of Na + , K + and Mg 2+ in waters from both catchments were greater than the reference values for unpolluted rivers [1], as previously reported for streams in Galicia [16,19]. In contrast, the mean concentration of Ca 2+ in waters from both catchments was lower than the reference value for unpolluted rivers. For instance, Na + , K + and Mg 2+ mean concentrations in Valiñas were, respectively, 4.5, 21.8 and 1.3 times greater than the reference values, while Ca 2+ concentrations were 0.6 times that of the reference. In Abelar, Na + , K + and Mg 2+ mean concentrations were, respectively, 2.2, 7.6 and 1.1 times greater than their corresponding reference, respectively; however, Ca 2+ concentrations in Abelar were 20% those of the reference for unpolluted waters. Therefore, the concentrations observed in the current study can be considered high for three of the macronutrients studied and low for Ca 2+ . This can be explained by the lithology of the area, consisting of rocks with low calcium contents. Nevertheless, the magnitudes observed are in agreement with previous reports in the region. For instance, high concentrations of Na + (from 8.74 to 34.35 mg L −1 ) have been reported in Galician streams due to the oceanic influence and the bedrock materials [19,[37][38][39][40]. However, the Mg 2+ concentrations measured in the current study were greater than those reported for other Galician streams [38][39][40].
Large differences in the magnitude of the macronutrients studied were detected between the two catchments. On average, Na + concentrations in the waters from the Valiñas catchment were two-fold those from Abelar, whereas K + and Ca 2+ concentrations were almost 3 times greater in waters from Valiñas than from Abelar; the concentration of Mg 2+ was 20% greater in waters from Valiñas than from Abelar. The proximity to the sea is one of the main factors influencing the Na + contents in rainfall, fog and water bodies [41] since these elements can reach the system through atmospheric deposition [13]. Therefore, the greatest concentrations of this solute in water samples from Valiñas were likely caused by a shorter distance (9 km) to the Atlantic Ocean than the Abelar catchment (30 km). Nevertheless, the nature of the soils and the lithology, as well as land use, might have also affected the concentrations of macronutrients in waters from both catchments. A great difference between Valiñas and Abelar regarding soil origin and bedrock exists. There is mainly granite in Valiñas and schists in Abelar. This would have led to differences between catchments in the base levels of each solute in the streams as well as in the process of rock weathering [42].

Human Activity Effects on Macronutrient Concentrations
Apart from the main features of catchments, the impact of anthropic activities on the concentrations of solutes detected in surface waters can be relevant and, sometimes, this impact is evident due to changes in magnitude, seasonality and long-term trends [11]. For instance, Na + concentrations have been reported to increase over time due to human influence (application of deicing salt) in Swiss catchments [11]. In the current study, the Na + and K + highest concentrations in the Abelar catchment seem to be originated from slurry applications, when this catchment was devoted to crop production, which buffered the seasonal oscillations expected. When the land use in this catchment changed after afforestation, the concentrations of these solutes in the stream were lower and tended to decrease over time. In the case of the Valiñas catchment, the concentrations of these solutes did not have large variations over time or within a given year, suggesting that processes of fertilization and transfer to the stream are limited, except for some high concentrations of K + in spring 2003 and summer 2014 (Figure 5a), which may indicate fertilizer leaching to the stream.
In contrast, Ca 2+ and Mg 2+ are geogenic solutes [11,43] and their concentrations in the streams are related to the nature of the bedrock materials existing in each of the studied catchments. Unexpectedly, the concentrations of Ca 2+ observed in the waters from Abelar were lower than those found in Valiñas because the bedrock in Abelar consists of schists, which usually have higher contents of Ca 2+ than granite. Nevertheless, bedrock composition does not determine the composition of water. The concentration of Ca 2+ is controlled by other factors such as carbonate content, pH and cationic exchange processes. In the case of the Abelar catchment, the highest concentration of Ca 2+ in the water was detected after an episode of slurry application by the end of October 2008 (Figure 6a), in agreement with former studies in which high concentrations of this solute in the streams were related to agricultural activities [11,16,38,43]. After catchment afforestation, the concentrations of this solute in the water were reduced. The concentrations of Mg 2+ in waters from Abelar were lower than in waters from Valiñas and increased in magnitude due to slurry applications. After afforestation, the concentrations of this solute were reduced and kept constant over time.

Concentration-Discharge Relationships
The current study revealed that concentration-discharge relations for Na + , K + , Ca 2+ and Mg 2+ differed between two catchments with different sizes, geological materials, land uses and located at different distances from the sea. According to the criteria proposed by Herdon et al. [44], limiting chemostatic behavior as −0.1 < b < 0, the C-Q curves observed in the current study indicated that the concentrations of Na + , Mg 2+ and Ca 2+ in Abelar followed a chemostatic behavior, which is usually attributed to homogeneous and uniform distribution of elements in a given catchment or, alternatively, to temporal stability of flow paths [3,8,43]. In contrast, the concentrations of the four macronutrients in Valiñas showed a negative chemodynamic behavior that can be attributed to dilution during high flows [45][46][47]. Finally, K + concentrations in Abelar showed a positive chemodynamic behavior that is usually attributed to enhanced erosion during high flows or to a reconnection of surface or subsurface catchment components where a given solute is more abundant [46,48]. It must be noted that the threshold of b > 0.1 is arbitrary and the use of Student's t tests are preferred [49]. In the current study, we employed this methodology and, despite the low values, the exponent b was significant in all cases, except for Na + in the Abelar catchment. Therefore, according to this criterion, all solutes showed a chemodynamic behavior, except for Na + in Abelar.
In the current study, the approach by Meybeck and Moatar [32], consisting of splitting the hydrograph at the median flow and calculating separate log-C log-Q slopes was followed [4,32]. Although this approach has been criticized because of being arbitrary, Diamond and Cohen [50] tested several breaking points for the C-Q curves and most of them were centered on the median flow, supporting our choice. In this study, depending on the solute and the catchment, five combinations were observed: down-down, down-flat, flat-down, up-up and up-down. In Valiñas, C-Q types were down-down for Na + , Ca 2+ and Mg 2+ , whereas a down-flat behavior was observed for K + . In Abelar, C-Q types were flat-down for Na + , up-up for K + and Ca 2+ and an up-down pattern was observed for Mg 2+ . Therefore, in agreement with observations made in French catchments [4], we found that less than 13% of all C-Q slopes were flat across the hydrograph, contrasting with previous studies that suggested the predominance of chemostatic C-Q relationships for major ions [8,51]. This can be explained by the hydrograph segmentation technique, which avoids confounding patterns for low and high flows [32].
In the Valiñas catchment, the negative relationship between Na + concentration and discharge suggests that the sodium that reaches the soil through atmospheric deposition is diluted during the wet season, as reported for other catchments [5,11]. In contrast, Na + concentrations in the Abelar catchment showed a chemostatic behavior at base flow conditions [8], since they remained constant despite the changes in discharge. However, when considering flows greater than the median, a slight trend to decrease Na + concentration with increasing discharge was observed in Abelar. The differences between catchments on the C-Q relationships for this solute can be attributed to catchment characteristics and land use [5,11,43].
The other solutes studied (K + , Ca 2+ and Mg 2+ ) showed a dilution effect in the Valiñas catchment, except for Mg 2+ at high flow conditions, which remained chemostatic. This suggests a source-limited behavior. However, the increase of K + and Ca 2+ concentrations during high flow conditions in Abelar suggests a source of these nutrients from sites in the system not normally accessed during base flow [5]. A recent study showed, through numerical simulations, that the vertical distribution of solute generation and activation thresholds of solute export are more influential in shaping the solute behavior (C-Q relation) than the timing of solute input [49]. The depth at which the solute is supplied or generated plays a major role in defining the b exponent for highly (Na + ) and weakly (Ca 2+ , Mg 2+ , and K + ) diluting solutes. Sporadic input applications increase temporal variability and add uncertainty to the solute behavior, the response of which remains the most difficult to constrain [49].
Finally, it is important to highlight that our data consists of low-frequency measurements, which allow for a powerful description of general patterns, although they obscure short-term C-Q dynamics, when hydrology, biological activity, and antecedent distribution of solute sources likely play a major role [4,46,52,53].

Conclusions
The long-term series of water quality analyzed in the current study allowed for assessing the signature of catchment features and the impact of human activities on Na + , K + , Ca 2+ and Mg 2+ concentrations observed in two small Atlantic watersheds. Changes in the magnitude of these major soluble ions were clearly related to human activities in the smallest catchment (Abelar), overwhelming the natural climatic variability. In contrast, catchment characteristics, such as distance to the coast and bedrock composition, had a significant influence on the concentrations of the macronutrients studied. Nevertheless, trends in long-term dynamics were not observed for the considered solutes in any of the studied catchments.
The analysis of the empirical C-Q power laws allowed the classification of the solute behaviors, which depended on the catchment features. While the four macronutrients showed a dilution pattern in the catchment located at the shortest distance to the coast (Valiñas), Na + tended to a chemostatic behavior in Abelar. Surprisingly, K + and Ca 2+ were positively related to streamflow in Abelar, which can be explained by the land use within the catchment, resulting in an up behavior for biologically reactive solutes at low flows and hydrological exports at high flows.
Despite the coarse resolution of the water quality data set, along with the relatively small sample size, the obtained results reinforce and expand the current understanding of the biogeochemical response of rivers, demonstrating that long-term observations allow for identifying several aspects of anthropic activities and watershed features on the concentrations of macronutrients in streams.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/ 10.3390/hydrology8030096/s1, Table S1: p-values from the Kolmogorov-Smirnoff tests comparing the temporal (on a yearly basis) distribution of sodium and potassium concentrations in waters from the Valiñas catchment. Yellow indicates those differences that are not significant, Table S2: p-values from the Kolmogorov-Smirnoff tests comparing the temporal (on a yearly basis) distribution of calcium and magnesium concentrations in waters from the Valiñas catchment. Yellow indicates those differences that are not significant, Table S3: p-values from the Kolmogorov-Smirnoff tests comparing the temporal (on a yearly basis) distribution of sodium, potassium, calcium and magnesium concentrations in waters from the Abelar catchment. Yellow indicates those differences that are not significant, Figure  Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.