Application of the RAPS Method for Determining the Dependence of Nitrate Concentration in Groundwater on the Amount of Precipitation

: Protecting groundwater from contamination is today’s most current environmental protection topic. What can man do in his environment to reduce the harmful impact of contamination on the environment, and thus the immediate effect on groundwater? Agricultural production is an ongoing source of groundwater contamination due to the increasingly frequent use of nitrates in fertilizers, which are washed out from the soil into groundwater due to precipitation. This paper investigates three wellﬁelds in the north of the Republic of Croatia near the town of Varaždin. With the application of the RAPS method, the dependence of nitrate concentration in groundwater on the amount of precipitation was established. The analysis results show the connection of the observed parameters, especially in the upper aquifer layer. In this layer, the coefﬁcients of correlation are greater than 0.80 at all locations, which shows a strong positive connection between the parameters. In the lower aquifer, the values of the coefﬁcients of correlation are lower, and the results mostly indicate a weak correlation. The obtained results will serve as a starting point for future studies, which will aim to precisely determine the factors that inﬂuence groundwater quality in the observed area.


Introduction
A high-quality water supply is one of the necessary prerequisites for sustainable development. Still, in the last few decades, due to increasingly intense anthropogenic influences, there has been a significant strain on the water supply system. Agricultural production is recognized worldwide as one of the primary sources of groundwater contamination, primarily due to nitrates, an important component of mineral fertilizers. Plant growth, development, and stable yields are only possible with sufficient nitrogen. However, the nitrates that the plant does not receive are easily and quickly washed out from the soil into the groundwater by rainwater. Nitrates introduced into the soil with mineral or organic fertilizers are easily soluble; they do not bind to the colloidal complex of the soil and can be washed into the surface or groundwater. Under favorable conditions, nitrates in groundwater can be accumulated for an extended period and transported over long distances, thereby risking the contamination of intake structures.
Over the past years, when the negative impacts of industrial and agricultural development and urban areas on the environment have been increasingly apparent, groundwater quality and its changes in space and time have become the subject of numerous scientific investigations. Madison and Brunett (1985), in a study on nitrates in groundwater, determined indicative limits of nitrate concentrations in groundwater based on which it is Earth 2023, 4 267 possible to identify anthropogenic influence [1]. The results of the study show that nitrate concentrations less than 1 mg/L NO 3 − can be attributed to natural nitrates in the aquifer, concentrations in the interval 1-13 mg/L NO 3 − mark the transitional area, i.e., concentrations that may or may not indicate anthropogenic influence; and nitrate concentrations above 13 mg/L NO 3 − prove anthropogenic influence. Addiscott et al. (1991) monitored the increase in nitrate concentration over time on a global scale. They concluded that the increase in nitrate concentration positively correlates with the increase in agricultural production [2]. A high concentration of nitrates in groundwater originating from agriculture is also reported by Obeidat et al. (2007), who stated that because nitrates are highly mobile and present in agricultural and industrial waste as well as household waste, they can be considered to be an indicator of groundwater contamination [3]. Nemčić-Jurec et al. (2007) indicated the impact of agriculture as a diffuse and point source of contamination with drinking water nitrates in northwestern Croatia [4]. Three different locations were investigated for potential contamination. At a place where no contamination was expected, which was used as a control group, a low nitrate concentration of 4.6 mg/L NO 3 − was found. In the other two locations, an area of intensive agriculture and a suburban area, it was shown that the concentration of nitrates was statistically significantly different and amounted to 28.7 and 26.5 mg/L NO 3 − , which indicates that the application of mineral and organic fertilizers and sewage are primary in sources of contamination. Recent research shows that a regression analysis of aggregated data can efficiently estimate an average contribution of point and diffuse contamination sources to the nitrate concentration in the groundwater at Varaždin wellfield [5]. In recent times, numerical models are often used to manage water resources. Šrajbek et al. (2022), using numerical simulations, located and quantified the sources of nitrate contamination in the influential area of the Varaždin wellfield, which was put into a state of rest due to the constant high concentration of nitrates [6].
The research area in this paper includes three wellfields located in the area of very intensive agricultural production. Earlier research has shown that agricultural production near the wellfields has a negative effect on groundwater quality, primarily by increasing the concentration of nitrates. As nitrates are washed out of the soil by precipitation and transported into the groundwater, this work aims to establish the connection between the concentration of nitrates in the groundwater at the observed wellfields and the amount of precipitation. Of course, groundwater quality depends on many factors, such as geology and soil properties, relief, and slopes determining surface runoff and seepage. Additionally, it depends on the quality of precipitation, sediments in the ground, and the duration of contact between the water and the rock. Given the above, it is clear that this question is extremely complex, and it demands detailed field and laboratory research. As such researches are very often long and expensive, the goal of this paper is to determine a connection between two basic parameters, namely, the amount of precipitation and the concentration of nitrates in groundwater. The results will serve as a basis for further research.

Research Area
The research area is located in the north of the Republic of Croatia and includes the space of three wellfields in the vicinity of the town of Varaždin (Figure 1). The groundwater quality is good, and the quality parameters are generally uniform at all observed wellfields, except for the concentration of nitrates, which varies significantly in time and space. Earlier research has shown that groundwater quality in the lower aquifer is considerably higher than in the upper aquifer [7,8] (Table 1). In addition, a statistically significant difference in nitrate concentration in aquifers at the observed wellfields was found [9]. In the regional water supply network, nitrate concentration was reduced by increasing pumping at wellfields where the values were below the minimum allowed concentration [10]. Due to the constantly high concentration of nitrates, the Varaždin wellfield has not been in operation since 2004. regional water supply network, nitrate concentration was reduced by increasing pumping at wellfields where the values were below the minimum allowed concentration [10]. Due to the constantly high concentration of nitrates, the Varaždin wellfield has not been in operation since 2004.

Geological and Hydrogeological Conditions of the Varaždin Aquifer
The geological structures consist of coarse-grained sediments of gravel and sand (middle and upper Pliocene age), whereas the underlying deposits are composed of silt and clay (upper Miocene age) [11].
The flow direction in the aquifer is parallel to the course of the Drava River, and its thickness increases from the west to the east. At the far west end, the thickness is 5 m, whereas in the town of Varaždin, it is about 75 m. In the eastern part, the thickness of the aquifer increases to its maximum value of 148 m and then decreases gradually. The aquifer is composed of two permeable gravel and sand layers divided mainly by a semipermeable silty and clay aquitard [12] (Figure 2). The aquifer's upper layer is formed of coarsegrained gravel and sand, whereas the lower part is mainly finer-grained gravel with sand. The groundwater flow is in the west-east direction, parallel to the Drava river.
The hydraulic conductivity values change slightly in the plan from 100 to 300 m/day, so it is higher in the west and lower in the east, resulting from the size drop of the granulometric fractions in the west-east direction. At the location of the Varaždin wellfield, the hydraulic conductivity values were measured to be in the range of 185-241 m/day for the upper and 92.6 m/day for the lower layer [13].

Geological and Hydrogeological Conditions of the Varaždin Aquifer
The geological structures consist of coarse-grained sediments of gravel and sand (middle and upper Pliocene age), whereas the underlying deposits are composed of silt and clay (upper Miocene age) [11].
The flow direction in the aquifer is parallel to the course of the Drava River, and its thickness increases from the west to the east. At the far west end, the thickness is 5 m, whereas in the town of Varaždin, it is about 75 m. In the eastern part, the thickness of the aquifer increases to its maximum value of 148 m and then decreases gradually. The aquifer is composed of two permeable gravel and sand layers divided mainly by a semipermeable silty and clay aquitard [12] (Figure 2). The aquifer's upper layer is formed of coarse-grained gravel and sand, whereas the lower part is mainly finer-grained gravel with sand. The groundwater flow is in the west-east direction, parallel to the Drava river.

Experimental Data
Groundwater at all three observed wellfields is drawn from the same alluvial aquifer from the upper and lower aquifers. As prescribed by the current ordinance [14], water was sampled at least once a week in active wells from 1993 to 2017. In the observed period  The hydraulic conductivity values change slightly in the plan from 100 to 300 m/day, so it is higher in the west and lower in the east, resulting from the size drop of the granulometric fractions in the west-east direction. At the location of the Varaždin wellfield, the hydraulic conductivity values were measured to be in the range of 185-241 m/day for the upper and 92.6 m/day for the lower layer [13].

Experimental Data
Groundwater at all three observed wellfields is drawn from the same alluvial aquifer from the upper and lower aquifers. As prescribed by the current ordinance [14], water was sampled at least once a week in active wells from 1993 to 2017. In the observed period, four wells from the upper and four from the lower aquifer were active at the Bartolovec wellfield. Eight wells from the upper and one well from the lower aquifer were active at the Varaždin wellfield. In comparison, three wells were active at the Vinokovščak wellfield, which, due to the wedging out of the poorly permeable layer at that location, captured water from a single aquifer.
Based on the chemical analysis of the groundwater quality at the three observed water wellfields, databases were formed, and further studies were conducted based on those databases. As water is pumped separately from both aquifers at the Bartolovec and Varaždin wellfields, analyses were carried out separately for each layer at these wellfields. The results related to the upper aquifer are labeled Bartolovec 1 and Varaždin 1, and the results related to the lower aquifer are marked Bartolovec 2 and Varaždin 2. At the Vinokovščak wellfield, the poorly permeable clay layer that divides the aquifer into two parts wedges out, and there is a unique aquifer layer at that location.
Based on the results of nitrate concentration measurements at all wells, at all three observed wellfields, in both aquifers separately, average monthly values of nitrate concentration in the period from 1993 to 2017 were calculated. Additionally, based on the daily results of measuring the amount of precipitation in the same period, the monthly precipitation amounts were calculated ( Figure 3).

Experimental Data
Groundwater at all three observed wellfields is drawn from the same alluvial aquifer from the upper and lower aquifers. As prescribed by the current ordinance [14], water was sampled at least once a week in active wells from 1993 to 2017. In the observed period, four wells from the upper and four from the lower aquifer were active at the Bartolovec wellfield. Eight wells from the upper and one well from the lower aquifer were active at the Varaždin wellfield. In comparison, three wells were active at the Vinokovščak wellfield, which, due to the wedging out of the poorly permeable layer at that location, captured water from a single aquifer.
Based on the chemical analysis of the groundwater quality at the three observed water wellfields, databases were formed, and further studies were conducted based on those databases. As water is pumped separately from both aquifers at the Bartolovec and Varaždin wellfields, analyses were carried out separately for each layer at these wellfields. The results related to the upper aquifer are labeled Bartolovec 1 and Varaždin 1, and the results related to the lower aquifer are marked Bartolovec 2 and Varaždin 2. At the Vinokovščak wellfield, the poorly permeable clay layer that divides the aquifer into two parts wedges out, and there is a unique aquifer layer at that location.
Based on the results of nitrate concentration measurements at all wells, at all three observed wellfields, in both aquifers separately, average monthly values of nitrate concentration in the period from 1993 to 2017 were calculated. Additionally, based on the daily results of measuring the amount of precipitation in the same period, the monthly precipitation amounts were calculated ( Figure 3).  Monthly amounts of precipitation represent the optimal time measure because they include the entire annual course of precipitation amounts. Therefore, the basic characteristic of the yearly revolution of monthly precipitation amounts in the observed area is a maximum in July, whereas the minor monthly precipitation amounts fall in January and February, which is a characteristic of the continental precipitation regime (Figure 4).
Monthly amounts of precipitation represent the optimal time measure because they include the entire annual course of precipitation amounts. Therefore, the basic characteristic of the yearly revolution of monthly precipitation amounts in the observed area is a maximum in July, whereas the minor monthly precipitation amounts fall in January and February, which is a characteristic of the continental precipitation regime (Figure 4).

RAPS Method and Correlation Intensity
As can be seen from Figure 2, large oscillations were observed in the monthly amounts of precipitation, which is understandable considering the nature of that parameter, i.e., the annual trend in the amounts of precipitation shown but also the deviations from the average monthly quantities in specific years. However, there were years when there was no precipitation at all in certain months as well as years when the amount of precipitation exceeded a month's average many times over. Therefore, the RAPS method described below was used for further analysis due to large data oscillations.
The RAPS (rescaled adjusted partial sums) method uses the sum deviation curve to analyze the observed quantities' time distribution. It enables overcoming small systematic and random changes, errors, and variability in the researched time series. The graphic display points to several sub-periods and trends, sudden jumps and declines in value, irregular fluctuations, periodicity, and others [15,16]. So far, this method has been used to analyze the series of connections of hydrological and meteorological indicators and determine hydrological regimes of watercourses [17][18][19][20]. It was also used to analyze the interconnection between wastewater quality indicators and temperature and precipitation [21], and this method was used to research nitrates' impact on groundwater resources [22].
The following equation defines the expression for the calculation of the RAPS method: whereas: Yi-the value of an individual sample i = 1,2,…, N, Y -average value of the observed sample, SY-the value of the standard deviation of the time series, N-the amount of data in the time series.
The RAPS method was used to analyze monthly precipitation and average monthly nitrate concentrations at all locations in 1993-2017, and an attempt was made to establish their interrelationship. For this purpose, correlation coefficients were calculated. Correlation describes the degree of connection between two variables, where the correlation coefficient takes a value between −1 and 1. The positive value of the correlation coefficient indicates the proportionality of the observed variables; that is, the growth of the value of both variables. A negative correlation coefficient indicates an increase in the value of one and a decrease in the value of the other variable. The closer the correlation coefficient

RAPS Method and Correlation Intensity
As can be seen from Figure 2, large oscillations were observed in the monthly amounts of precipitation, which is understandable considering the nature of that parameter, i.e., the annual trend in the amounts of precipitation shown but also the deviations from the average monthly quantities in specific years. However, there were years when there was no precipitation at all in certain months as well as years when the amount of precipitation exceeded a month's average many times over. Therefore, the RAPS method described below was used for further analysis due to large data oscillations.
The RAPS (rescaled adjusted partial sums) method uses the sum deviation curve to analyze the observed quantities' time distribution. It enables overcoming small systematic and random changes, errors, and variability in the researched time series. The graphic display points to several sub-periods and trends, sudden jumps and declines in value, irregular fluctuations, periodicity, and others [15,16]. So far, this method has been used to analyze the series of connections of hydrological and meteorological indicators and determine hydrological regimes of watercourses [17][18][19][20]. It was also used to analyze the interconnection between wastewater quality indicators and temperature and precipitation [21], and this method was used to research nitrates' impact on groundwater resources [22].
The following equation defines the expression for the calculation of the RAPS method: whereas: Y i -the value of an individual sample i = 1,2, . . . , N, Y-average value of the observed sample, S Y -the value of the standard deviation of the time series, N-the amount of data in the time series. The RAPS method was used to analyze monthly precipitation and average monthly nitrate concentrations at all locations in 1993-2017, and an attempt was made to establish their interrelationship. For this purpose, correlation coefficients were calculated. Correlation describes the degree of connection between two variables, where the correlation coefficient takes a value between −1 and 1. The positive value of the correlation coefficient indicates the proportionality of the observed variables; that is, the growth of the value of both variables. A negative correlation coefficient indicates an increase in the value of one and a decrease in the value of the other variable. The closer the correlation coefficient value is to 0, the weaker the connection, and the closer to 1, whether positive or negative, the stronger the relationship. The correlation's intensity is defined based on the correlation coefficients' value (

Results and Discussion
At the Bartolovec 1 wellfield, a connection between the observed data sets was observed, and three sub-periods can be discerned from the RAPS data time series diagram: the first from the beginning of the observed period to 1999, the second from 1999 to 2012, and the third after that year ( Figure 5). The value of the correlation coefficient for the entire period is 0.51. However, as the analysis of the time series of RAPS data indicated the existence of three sub-periods, the linear correlations were analyzed separately for each sub-period. The first sub-period is characterized by an upward trend, and the value of the correlation coefficient is 0.85, which indicates a strong connection ( Figure 6). In the second sub-period, the trend is downward, and the correlation coefficient value is 0.81, showing a strong relationship. The third sub-period is characterized by separating the RAPS data curves of the observed variables, which confirms the value of the correlation coefficient of −0.24. Based on the obtained results, it can be concluded that at this location, the amount of precipitation affects the nitrate concentration in the groundwater in the first and second sub-periods. At the same time, such a connection was absent in the third sub-period. tion coefficients' value (

Results and Discussion
At the Bartolovec 1 wellfield, a connection between the observed data sets was observed, and three sub-periods can be discerned from the RAPS data time series diagram: the first from the beginning of the observed period to 1999, the second from 1999 to 2012, and the third after that year ( Figure 5). The value of the correlation coefficient for the entire period is 0.51. However, as the analysis of the time series of RAPS data indicated the existence of three sub-periods, the linear correlations were analyzed separately for each subperiod. The first sub-period is characterized by an upward trend, and the value of the correlation coefficient is 0.85, which indicates a strong connection ( Figure 6). In the second sub-period, the trend is downward, and the correlation coefficient value is 0.81, showing a strong relationship. The third sub-period is characterized by separating the RAPS data curves of the observed variables, which confirms the value of the correlation coefficient of −0.24. Based on the obtained results, it can be concluded that at this location, the amount of precipitation affects the nitrate concentration in the groundwater in the first and second sub-periods. At the same time, such a connection was absent in the third sub-period.  Three sub-periods can also be discerned from the curve diagram of the time series of RAPS data at the Bartolovec 2 wellfield: the first until 2002, the second from 2002 to 2012, and the third after that year (Figure 7). The correlation analysis did not show a significant connection between the observed data sets in any sub-period, so the correlation coefficients in the first and second sub-periods are −0.35 and 0.42, indicating a weak correlation. In the third sub-period, the value of the correlation coefficient is −0.66, which means a negative correlation (Figure 8). The analysis results indicate that in this location, the amount of precipitation does not affect or weakly affect the concentration of nitrates in the groundwater. Three sub-periods can also be discerned from the curve diagram of the time series of RAPS data at the Bartolovec 2 wellfield: the first until 2002, the second from 2002 to 2012, and the third after that year (Figure 7). The correlation analysis did not show a significant connection between the observed data sets in any sub-period, so the correlation coefficients in the first and second sub-periods are −0.35 and 0.42, indicating a weak correlation. In the third sub-period, the value of the correlation coefficient is −0.66, which means a negative correlation (Figure 8). The analysis results indicate that in this location, the amount of precipitation does not affect or weakly affect the concentration of nitrates in the groundwater.   Three sub-periods can also be discerned from the curve diagram of the time series of RAPS data at the Bartolovec 2 wellfield: the first until 2002, the second from 2002 to 2012, and the third after that year (Figure 7). The correlation analysis did not show a significant connection between the observed data sets in any sub-period, so the correlation coefficients in the first and second sub-periods are −0.35 and 0.42, indicating a weak correlation. In the third sub-period, the value of the correlation coefficient is −0.66, which means a negative correlation (Figure 8). The analysis results indicate that in this location, the amount of precipitation does not affect or weakly affect the concentration of nitrates in the groundwater.  Since the Varaždin 1 wellfield was put into a state of rest in 2004, groundwater has not been pumped at that location. Following that, the number of samplings is practically negligible and the analysis is only possible until that year. The curves of the time series of RAPS data are well aligned during the entire observed period, and considering the trend, two sub-periods can be discerned: the first until 1999 and the second after that year ( Figure  9). The first sub-period is characterized by an upward trend, and the second by a downward trend. The values of the correlation coefficients in both sub-periods are high and amount to 0.87 for the first and 0.95 for the second sub-period ( Figure 10). Because the points on the dispersion diagram do not discern two data sets, the correlation coefficient for the entire period was calculated and amounted to 0.91. This analysis confirms the high Since the Varaždin 1 wellfield was put into a state of rest in 2004, groundwater has not been pumped at that location. Following that, the number of samplings is practically negligible and the analysis is only possible until that year. The curves of the time series of RAPS data are well aligned during the entire observed period, and considering the trend, two sub-periods can be discerned: the first until 1999 and the second after that year (Figure 9). The first sub-period is characterized by an upward trend, and the second by a downward trend. The values of the correlation coefficients in both sub-periods are high and amount to 0.87 for the first and 0.95 for the second sub-period ( Figure 10). Because the points on the dispersion diagram do not discern two data sets, the correlation coefficient for the entire period was calculated and amounted to 0.91. This analysis confirms the high correlation between the two variables, i.e., it shows that the amount of precipitation significantly affects the nitrate concentration value in the groundwater at the location of the Varaždin 1 wellfield.
not been pumped at that location. Following that, the number of samplings is practically negligible and the analysis is only possible until that year. The curves of the time series of RAPS data are well aligned during the entire observed period, and considering the trend, two sub-periods can be discerned: the first until 1999 and the second after that year ( Figure  9). The first sub-period is characterized by an upward trend, and the second by a downward trend. The values of the correlation coefficients in both sub-periods are high and amount to 0.87 for the first and 0.95 for the second sub-period ( Figure 10). Because the points on the dispersion diagram do not discern two data sets, the correlation coefficient for the entire period was calculated and amounted to 0.91. This analysis confirms the high correlation between the two variables, i.e., it shows that the amount of precipitation significantly affects the nitrate concentration value in the groundwater at the location of the Varaždin 1 wellfield.  The Varaždin 2 wellfield was put into operation in 2005, and until 2010, water was continuously pumped at that location. Due to the high concentration of nitrates, this location was also put into a state of rest in 2011, after which it is only occasionally included in the water supply system. Although after 2011, the number of groundwater samplings is slightly lower, this analysis was also carried out for the wellfield's passive operation period. The time series diagram of the RAPS data of monthly precipitation amount and monthly average nitrate concentration shows two sub-periods: the first until 2010 and the second after that year (Figure 11). There is a downward trend in the first sub-period and an upward trend in the second, and the RAPS data curves are well aligned with each other The Varaždin 2 wellfield was put into operation in 2005, and until 2010, water was continuously pumped at that location. Due to the high concentration of nitrates, this location was also put into a state of rest in 2011, after which it is only occasionally included in the water supply system. Although after 2011, the number of groundwater samplings is slightly lower, this analysis was also carried out for the wellfield's passive operation period. The time series diagram of the RAPS data of monthly precipitation amount and monthly average nitrate concentration shows two sub-periods: the first until 2010 and the second after that year ( Figure 11). There is a downward trend in the first sub-period and an upward trend in the second, and the RAPS data curves are well aligned with each other throughout the observed period. The value of the correlation coefficient in the first subperiod is 0.80, and in the second, 0.68 indicates a strong or good correlation ( Figure 12). Because the points on the dispersion diagram do not discern two data sets, the correlation coefficient for the entire period was calculated and amounted to 0.74. From the obtained results, it can be concluded that the amount of precipitation affects the concentration of nitrates in the groundwater at this location as well. the water supply system. Although after 2011, the number of groundwater samplings is slightly lower, this analysis was also carried out for the wellfield's passive operation period. The time series diagram of the RAPS data of monthly precipitation amount and monthly average nitrate concentration shows two sub-periods: the first until 2010 and the second after that year ( Figure 11). There is a downward trend in the first sub-period and an upward trend in the second, and the RAPS data curves are well aligned with each other throughout the observed period. The value of the correlation coefficient in the first subperiod is 0.80, and in the second, 0.68 indicates a strong or good correlation ( Figure 12). Because the points on the dispersion diagram do not discern two data sets, the correlation coefficient for the entire period was calculated and amounted to 0.74. From the obtained results, it can be concluded that the amount of precipitation affects the concentration of nitrates in the groundwater at this location as well.  At the Vinokovščak wellfield, the analysis was carried out from 1999, when the wellfield was put into operation, until 2017. The time series of RAPS data throughout the observed period is well aligned, and two sub-periods are discerned: the first until 2008 and the second after that year. The first sub-period is characterized by a downward trend, and the second by an upward trend ( Figure 13). The value of the correlation coefficient in the first sub-period is 0.92; in the second, this value is slightly lower and is 0.84, indicating a strong correlation (Figure 14). Because two data sets cannot be discerned from the points on the dispersion diagram, the correlation coefficient for the entire period was calculated and amounted to 0.86. Thus, the results of this analysis show that at the location of the Vinokovščak wellfield, the amount of precipitation significantly affects the nitrate concen- At the Vinokovščak wellfield, the analysis was carried out from 1999, when the wellfield was put into operation, until 2017. The time series of RAPS data throughout the observed period is well aligned, and two sub-periods are discerned: the first until 2008 and the second after that year. The first sub-period is characterized by a downward trend, and the second by an upward trend ( Figure 13). The value of the correlation coefficient in the first sub-period is 0.92; in the second, this value is slightly lower and is 0.84, indicating a strong correlation ( Figure 14). Because two data sets cannot be discerned from the points on the dispersion diagram, the correlation coefficient for the entire period was calculated and amounted to 0.86. Thus, the results of this analysis show that at the location of the Vinokovščak wellfield, the amount of precipitation significantly affects the nitrate concentration in the groundwater. served period is well aligned, and two sub-periods are discerned: the first until 2008 and the second after that year. The first sub-period is characterized by a downward trend, and the second by an upward trend ( Figure 13). The value of the correlation coefficient in the first sub-period is 0.92; in the second, this value is slightly lower and is 0.84, indicating a strong correlation (Figure 14). Because two data sets cannot be discerned from the points on the dispersion diagram, the correlation coefficient for the entire period was calculated and amounted to 0.86. Thus, the results of this analysis show that at the location of the Vinokovščak wellfield, the amount of precipitation significantly affects the nitrate concentration in the groundwater.

Conclusions
This paper used the RAPS method to establish the connection between the nitrate concentration in the groundwater at the observed wellfields and the amount of precipitation. The analysis results show that the amount of precipitation affects the concentration of nitrates in almost all observed wellfields. However, only at the Bartolovec 2 wellfield was no connection observed between the observed data sets. Based on the regression analysis, a positive relationship between the two data groups was established, which means that the nitrate concentration in the groundwater increases with precipitation. The values of the correlation coefficients in almost all observed locations show a good or strong correlation, so in the upper aquifer, the values are in the interval of 0.81-0.95. A greater connection between the two sets of data was recorded in the upper aquifer layer compared to the lower one, which is understandable considering that this layer is under significantly greater influence of anthropogenic sources of contamination on the surface. Considering the long-standing problem in the water supply of the observed area, primarily the high concentration of nitrates, high-quality drinking water will be one of the main prerequisites for the region's further development. To preserve or improve the quality of groundwater at the observed wellfields, it will be important to adhere to the principles that harmonize agricultural activity with ecological and economic principles, i.e., the way of land use, which can best achieve the goals of economic and environmental sustainability, i.e., the service of so-called good agricultural practices.

Conclusions
This paper used the RAPS method to establish the connection between the nitrate concentration in the groundwater at the observed wellfields and the amount of precipitation. The analysis results show that the amount of precipitation affects the concentration of nitrates in almost all observed wellfields. However, only at the Bartolovec 2 wellfield was no connection observed between the observed data sets. Based on the regression analysis, a positive relationship between the two data groups was established, which means that the nitrate concentration in the groundwater increases with precipitation. The values of the correlation coefficients in almost all observed locations show a good or strong correlation, so in the upper aquifer, the values are in the interval of 0.81-0.95. A greater connection between the two sets of data was recorded in the upper aquifer layer compared to the lower one, which is understandable considering that this layer is under significantly greater influence of anthropogenic sources of contamination on the surface. Considering the long-standing problem in the water supply of the observed area, primarily the high concentration of nitrates, high-quality drinking water will be one of the main prerequisites for the region's further development. To preserve or improve the quality of groundwater at the observed wellfields, it will be important to adhere to the principles that harmonize agricultural activity with ecological and economic principles, i.e., the way of land use, which can best achieve the goals of economic and environmental sustainability, i.e., the service of so-called good agricultural practices.

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

Conflicts of Interest:
The authors declare no conflict of interest.