Next Article in Journal
Development and Verification of the Available Number of Water Intake Days in Ungauged Local Water Source Using the SWAT Model and Flow Recession Curves
Previous Article in Journal
Biogeochemistry of Mediterranean Wetlands: A Review about the Effects of Water-Level Fluctuations on Phosphorus Cycling and Greenhouse Gas Emissions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geochemical Characteristics of Alluvial Aquifer in the Varaždin Region

1
Croatian Geological Survey, 10 000 Zagreb, Croatia
2
Croatian Waters, Central Water Management Laboratory, 10 000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Water 2021, 13(11), 1508; https://doi.org/10.3390/w13111508
Submission received: 4 May 2021 / Revised: 19 May 2021 / Accepted: 24 May 2021 / Published: 27 May 2021
(This article belongs to the Section Hydrology)

Abstract

:
The variation in the major groundwater chemistry can be controlled by dissolution and precipitation of minerals, oxidation-reduction reactions, sorption and exchange reactions, and transformation of organic matter, but it can also occur as a result of anthropogenic influence. The alluvial aquifer represents the main source of potable water for public water supply of the town Varaždin and the surrounding settlements. Sampling campaigns were carried out from June 2017 until June 2019 to collect groundwater samples from nine observation wells. Major cations and anions, dissolved organic carbon and nutrients were analyzed in the Hydrochemical Laboratory of Croatian Geological survey. The sampled waters belong to the CaMg-HCO3 hydrochemical type, except the water from observation well P-4039 that belongs to NaCa-HCO3 hydrochemical type. It was identified that groundwater chemistry is mainly controlled by hydrogeological environment (natural mechanism), but anthropogenic influence is not negligible. The results of this research have significant implications on sustainable coexistence between agricultural production and water supply.

1. Introduction

Natural waters acquire their chemical characteristics both by dissolution and by chemical reactions with solids, liquids and gases, with which they come into contact during the various phases of the hydrological cycle [1]. The chemical composition of groundwater is often used to investigate groundwater residence time, origin, flow direction and anthropogenic or natural contamination [2,3,4,5,6,7,8,9,10]. The variation in the major cations and anions of groundwater can be controlled by dissolution and precipitation of minerals, oxidation-reduction reactions, sorption and exchange reactions, and transformation of organic matter. In addition, major cations and anions can be added to the aquifer systems as a result of anthropogenic influence; for example calcium, magnesium, sodium, chloride and potassium are present in sludge, waste water, manure [11,12,13].
In recent decades, high nitrate concentrations emerged as a globally growing problem for drinking and agricultural purposes [14,15,16]. The adverse health effects of high nitrate levels in drinking water have been well documented, including gastric cancer, non-Hodgkin’s lymphoma, and methemoglobinemia [17,18,19]. In some parts of Croatia, uncontrolled and extensive agricultural production is causing the pollution of groundwater with nitrate. An example of an area where high nitrate concentrations are present is the alluvial aquifer in the northwestern part of Croatia, in the Varaždin region. To determine the impact of the hydrogeological environment and humans on groundwater chemical features of the alluvial aquifer, geochemical investigations were performed. The aquifer represents the main source of potable water for public water supply of the town of Varaždin and the surrounding settlements. The most important activity in the region is agricultural production (plantation of wheat, maize, cabbage, poultry and dairy farming). The present research describes the geochemical reactions that influence the chemical composition of groundwater, and evaluates the key controlling processes within the alluvial aquifer. The main aim of the paper is to assess the natural and human influence on the groundwater chemistry in the study area, in order to assure sustainable coexistence between agricultural production and water supply, and prevent further groundwater quality deterioration.

2. Description of the Study Area

The study area is located in the northwestern part of Croatia, in the Varaždin region. It belongs to the Black Sea catchment area. The aquifer, which is situated in the Drava river lowland, is characterized by intergranular porosity. The topography is characterized by wide flatlands surrounded by hills. The Drava River presents the aquifer boundary in the northwest and north, and in the southeastern part of the study area, the Plitvica stream flows (Figure 1).
The study area has a characteristic precipitation regime with more precipitation during the summer [20]. Consequently, the local climate is categorized in the Cfb group according to the Köppen–Geiger classification system, which is known as “warm-temperate climate” or “marine west coast climate.” The study area is the former [21]. The mean annual temperature is 10.6 °C, with January as the coldest month (average temperatures of 0.0 °C), and July and August as the warmest months (average temperatures of 20.9 and 20.1 °C, respectively) [9]. According to the data from the last climate normal period (1981–2010), the average annual precipitation was 832 mm, with a lower average precipitation amount during the cold part of the year (minimum in January with 38.7 mm), and a higher average precipitation amount during the warm part of the year (maximum in September with 98.3 mm) [9]. The favorable climate, topography and available groundwater have enabled intensive agricultural practices, including the application of synthetic fertilizers and manure that has subsequently led to high nitrate concentrations in the Varaždin aquifer.
The aquifer is composed of gravel and sand with variable portions of silt [22,23,24]. It was formed during the Pleistocene and Holocene as a result of accumulation processes of the Drava River [25]. At the utmost northwestern part of the study area, the aquifer thickness is less than 5 m, and it gradually increases in downstream direction, reaching its maximum of roughly 50 m in the eastern part of the study area (Figure 1). It is noticed that particle size changes going from the northwestern part downstream, i.e., the size of gravel and sand particles gets gradually smaller as result of the decrease in energy of the Drava River. Deposits of gravel and sand show stratification in some places, which is characterized by a sudden change in the size of pebbles, or an increased amount of sandy component [26]. Gneiss and quartz pebbles prevail, but there are also pebbles of basic and neutral eruptive rocks; limestone, dolomite, etc. [27]. The main constituents of the sand are quartz, feldspars and carbonate minerals, and it contains significant amounts of heavy minerals such as garnet, epidote, amphibole, rutile, kyanite, etc. [28,29]. Along with gravel and sand in the study area, there are also oxbow deposits that were deposited in the old Drava riverbed, where the still water environment of sedimentation remained for a long time. Various fine sediments were deposited, such as silt, clay and organic matter, forming distinctive facies of oxbow [28,29]. Remains of oxbows have been observed in the area of Strmec, Petrijanec, Otok Virje, Svibovec and Sračinec.
In the southeastern part of the study area, near Varaždin town, a tiny aquitard composed of clay and silt appears, dividing the aquifer into two hydrogeological units. The aquitard has regional significance, especially downstream outside the study area, but not so much in the study area, due to its small thickness. The covering layer of the aquifer is not continuously developed throughout the entire study area. In the central part and near the Drava River it rarely exceeds 50 cm, while often it completely disappears. Such conditions are favorable if they are considered from the aspect of aquifer recharge, but at the same time, tiny covering layers makes the aquifer quite vulnerable. The aquifer of the study area is unconfined, and is recharged by precipitation infiltration through unsaturated zones and by surface water percolation [9,30]. The general groundwater flow direction is NW–SE and is parallel to the Drava River (Figure 1). It is noteworthy that the groundwater flow net has been significantly changed since the building of a hydroelectric power plant in 1970s. Namely, prior to this intervention, the groundwater had flown towards the Drava River, which had represented the discharge zone, and now it is the recharge zone. After the construction of the Varaždin accumulation lake, the pressure head layout changed, leading to percolation of the lake water to the aquifer. At the same time, in the vicinity of the derivation channel, the groundwater level is lowered because the channel is deeply cut into the aquifer (Figure 1). Another discharge zone is the Plitvica stream, which drains the aquifer most of the time and recharges it only in high water level conditions.

3. Materials and Methods

Groundwater sampling campaigns were carried out from June 2017 until June 2019. Samples were collected from alluvial aquifer by pumping 9 observation wells (8 piezometric wells and 1 private well). The depths of observation wells and their screens (if known) are given in Table 1.
In situ parameters such as temperature, pH, dissolved oxygen (DO) and electrical conductivity (EC) were measured in the field using a WTW multi-probe. Alkalinity was also measured in the field by titration with 1.6 N H2SO4, using phenolphthalein and bromocresol green-methyl red as indicators, and then converted to the equivalent HCO3 concentrations. Samples for analysis of cations and anions were filtered through 0.45 µm cellulose membrane filters into the HPDE 500 mL bottles prior measuring on Ion Chromatographer Dionex ICS 6000, while low concentrations of NH4+, NO2 and PO43−-P were analyzed using spectrophotometer HACH DR 9000. Samples for measurement of dissolved organic carbon (DOC) were, in the field, collected into 100 mL dark glass bottles and analyzed using HACH QBD1200 analyzer. Samples were kept in the portable refrigerator during transport to the laboratory and analyzed in the evening of the same day. The ion balance errors for the analyses were checked by the relative deviation from charge balance (Δmeq = 100 × (Σmeq+ − Σmeq−)/(Σmeq+ + Σmeq−) < ±5%) [31,32]. Concentrations of dissolved metals in water were measured using inductively coupled plasma-mass spectrometry on Agilent 8900 ICP-MS Triple Quad with solution of 30 µg L−1 Ge, Y, In and Tb as internal standards according to HRN EN ISO 17294-2:2016 norm [33]. All measurements were performed in quintuplets. Quality control of the ICP-MS method was performed by the analysis of the elements of interest in certified reference material Anas-38 (Inorganic Venture) at the beginning and after analyzing each series of samples. Calibration lines for each element and internal standards were made using Agilent multi-element calibration standard solutions and internal standard mix solution. Before the analysis, the samples were filtered through a 0.45 µm filter on the field, and acidified with ultra-pure 6 N HNO3 acid. The PHREEQC software was used to determine saturation indices and CO2 pressure [34]. The determination of the redox state within the aquifer was performed using McMahon and Chapelle’s methodology [35]. The correlation diagrams and calculation of correlation coefficients were determined using MS Excel tool.

4. Results

The average, minimum and maximum values of the analyzed physicochemical parameters and metal concentrations in the groundwater samples are presented in Table 2a,b The calculated redox conditions of groundwater are given in Table 3.
In situ parameters are presented in Table 2a. The EC values ranged from 490 to 1091 μS/cm, with the highest values measured in water samples from P-4039 and the lowest values in SPV-11. The highest EC values in water from P-4039 are a consequence of a high concentration of dissolved solids, especially sodium and chloride ions (Table 2a). On the other hand, the water from SPV-11 has the lowest values because of the influence of the Drava River on the alluvial aquifer (dilution effect). The groundwater temperatures ranged from 9.4 to 16 °C. The pH values of groundwater ranged from 6.86 to 7.76, meaning that the waters are mildly acid to alkaline. The DO ranged from 0.2 to 7.1 mg/L. Calculated redox category follows the DO values. Low DO values in the aquifer are accompanied by mixed (oxic-anoxic) conditions and if there is a deficiency in DO, anoxic conditions prevail (Table 3).
The order of dominance ions among cations in waters of observation wells PDS-5, PDS-6, PDS-7, P-1529, P-1556, P-2500, SPV-11 and the private well is Ca2+ > Mg2+ > Na+ > K+, while in water of observation well P-4039 is Na+ > Ca2+ > Mg2+ > K+ (Table 2a). In a case of anions, the order is HCO3 > NO3 > SO42− > Cl in waters of wells PDS-5, PDS-6, PDS-7, P-1529, P-2500 and the private well, but for wells P-1556, SPV-11 and P-4039, it differs. Nitrate concentration in wells PDS-5, PDS-6, PDS-7, P-1529, P-2500 and the private well exceed the maximum contaminant level (MCL) of 50 mg/L NO3 [36] for most of the observed time due to agricultural practices and waste water from surrounding settlements. The order of dominance ions among anions for wells P-1556 and SPV-11 is HCO3 > SO42− > NO3 > Cl, and for well P-4039 is HCO3 > Cl > SO42− > NO3. SO42− concentrations did not exceed MCL of 250 mg/L in the analyzed samples from all wells, but Cl concentrations in water from P-4039 occasionally exceed MCL value of 250 mg/L because of the seasonal de-icing of roads, and sewage waters from semipermeable septic tanks (Table 2a).
Concentrations of nitrite and ammonia in all samples did not exceed MCL values of 0.5 mg/L, ranging from below detection limit to 0.24 mg/L for NO2, and from below detection limit to 0.16 mg/L for NH4+ (Table 2a). Orthophosphate concentrations occasionally exceed MCL value of 0.3 mg/L in waters from wells P-1529, P-1556, P-2500, PDS-5 and PDS-7 (Figure 2).
Table 2b and Figure 3 show that water samples from wells P-1556, P-4039 and PDS-7 contain high concentrations of some heavy metals, occasionally exceeding MCL values. High heavy metal concentrations are attributed to weathering of oxbow sediment, which contains heavy metals, combined with anthropogenic influence. In the rest of the wells, concentrations of heavy metals are very low. Dissolved iron concentrations ranged from 1.1 to 2265 µg/L, manganese concentrations from 0.08 to 237 µg/L, and zinc concentrations from 3.5 to 5665 μg/L (Table 2b and Figure 3). In addition, the waters of well P-4039 show high concentrations of Pb and Ni (Table 2b).

5. Discussion

According to the major ionic composition, sampled waters belong to the CaMg–HCO3 hydrochemical type, except the water from observation well P-4039, which belongs to the NaCa-HCO3 hydrochemical type (Figure 4). Such hydrochemical type of water is a consequence of the dissolution and weathering of carbonate (limestone, dolomite) and silicate minerals (micas, feldspar, etc.) that build aquifer sediments. Since weathering rates of limestone and dolomite are up to 80 and 12 times faster than silicate weathering rates [37], carbonate dissolution mainly dominates major ionic composition and presents the first geochemical process. The influence of silicate weathering, which is the second geochemical process, was analyzed by the bivariate mixing diagram of Na+-normalized Mg2+ versus Na+-normalized Ca2+ (Figure 5). There is no pronounced silicate weathering in the studied waters, but it was observed that the catchment areas of wells P-4039, P-2500 and occasionally P-1529 and P-1556 indicate the influence of silicate minerals weathering.
In order to determine which carbonate mineral predominantly weathers, molar ratio Mg2+/Ca2+ was used. Overall, it is observed that calcite dissolution is dominant over dolomite dissolution (Figure 6). However, in the catchment areas of wells SPV-11, P-4039, P-1529 and PDS-5, the dolomite dissolution is more pronounced than the calcite dissolution because the Mg2+/Ca2+ ratios are over 0.33 value (Figure 6).
The third geochemical process, cation exchange, was observed, taking into account the bivariate diagram of (Ca2+ + Mg2+) − (HCO3 + SO42−) versus Na+ − Cl (Figure 7). Concentrations of bivalent cations (Ca2+ and Mg2+) that may have been involved in exchange reactions were corrected by subtracting equivalent concentrations of associated anions (HCO3 and SO42−) that would be derived from other processes (e.g., carbonate or silicate weathering, where calcite and anorthite (Ca-feldspar) produce similar molar concentrations of Ca2+ and HCO3, but no SiO2 [38,39]). Similarly, Na+ that may be derived from the aquifer matrix can be accounted for by assuming that Na+ contributions of meteoric origin would be balanced by equivalent concentrations of Cl [40]. For active cation exchange taking place in the aquifer, the slope of this bivariate plot should be −1 [41]. Since the cation exchange process is well-pronounced in the catchment of well P-4039 and in the catchments of other wells is masked, two subfigures are given: Figure 7a showing all wells, and in Figure 7b, well P-4039 is left out. It is observed that the cation exchange process is not pronounced in the catchment of the rest of the wells (Figure 7b).
Mineral equilibrium calculations for groundwater are useful in predicting the presence of reactive minerals in the groundwater system and estimating mineral reactivity. By using the saturation index approach, it is possible to predict the reactive mineralogy of the subsurface from groundwater data, without collecting the samples of the solid phase and analyzing the mineralogy [42]. This approach was used, and the saturation indices (SI) of calcite and dolomite and partial pressure of CO2 were calculated. If the groundwater is saturated (SI > 0) with respect to the calcite and/or dolomite minerals, precipitation of calcite and dolomite minerals is possible. On the other hand, if the groundwater is undersaturated (SI < 0) with respect to minerals, dissolution would continue. Most of the time, the sampled groundwater is saturated with respect to calcite and undersaturated with respect to dolomite (Figure 8). Occasionally, especially during summer periods when water levels are decreasing, groundwater is saturated with respect to dolomite. Currently, partial pressure of CO2 is very low and enables precipitation. On the other hand, when water levels increase due to the rainy season, partial pressure of CO2 increases due to the flushing of the surface and unsaturated zone. The SI of both minerals decreases and becomes negative for dolomite and lower or negative for calcite.
Higher partial pressure of CO2 is connected with an increase in DOC in groundwater (Figure 9). In catchments of the observation wells P-1529, PDS-7, PDS-5, and PDS-6, it is observed that higher DOC concentrations in water are accompanied by higher partial pressure of CO2, as a consequence of the flushing of the organic matter from the soil and unsaturated zone into the aquifer. However, in the catchments of the observation wells SPV-11, P-1556, and P-2500, changes of DOC concentrations do not significantly affect the partial pressure of CO2, which is mainly controlled by the dissolution of carbonate minerals. In the aquifers with carbonate matrix, it is observed that increasing chemical weathering of carbonate minerals is related to increasing CO2 in groundwater [43].
The fourth process that influences the geochemical evolution of groundwater is anthropogenic influence, which is recognized through agricultural and urban activities. In Figure 10, Na+/Cl ratios show that values are mostly scattered around halite line and the most of samples are shifted to Cl side, indicating the influence of the waste water and manure. When cation exchange process is dominant, the values shift to the Na+ side. The source of the halite in the study area is also not natural, but anthropogenic. Halite is used during the winter period for de-icing of the roads.
Another indicator of anthropogenic influence is high nitrate concentration (Figure 11). Usually, the sources of nitrate are fertilizers, organic and mineral [44,45,46]. Although agricultural production is dominant in the research area, there has been a decrease in agricultural surfaces and the application of fertilizers in the past 10–15 years, followed by an increase in the urban area by 12% [47]. The construction of the sewerage network did not follow the urbanization of the study area, and nitrate pollution may also occur due to the discharge of waste water into the ground. In the catchments of observation wells that are situated in or close to urban area, a positive correlation between nitrate and chloride and a good connection between nitrate and phosphates was observed (Table 2a). It is generally known that source of phosphates is waste water [48]. In addition, bromides were observed in waters of those observation wells, especially during the wet period of the year, and Br/Cl ratio confirms the influence of the waste water (Table 2a). The highest concentrations of nitrate are observed in the middle of the study area, where the intensive agricultural production and urban areas exist. The observation wells that are close to the Drava river have low nitrate concentrations, because the river recharges the alluvial aquifer [9,30] and causes a dilution effect. Moreover, sessional periodic variations of nitrate and chloride were observed (Figure 11). During the winter season and early spring, high chloride concentrations are measured in groundwater samples due to the de-icing of roads. Conversely, high nitrate concentrations are measured during intensive agricultural production and irrigation, which occurs in the late spring–summer season.
In addition to these two anthropogenic indicators, sulfate is also a relevant indicator because it can be released into the groundwater as part of domestic waste water [49]. However, the highest concentrations in groundwater are usually from natural sources such as gypsum, anhydrite, oxidation of sulfide minerals, etc. [49]. In the research area, the origin of sulfate is natural, because concentrations in all catchments are similar.

6. Conclusions

The alluvial aquifer in the Varaždin region is an important groundwater source for human consumption and the dependent ecosystem. Therefore, it is vital to ensure the sustainable use of this valuable water resource. The conducted research, based on the chemical analyses of groundwater samples from nine observation wells, identified four main processes that influence the groundwater chemistry:
(a)
The dissolution and precipitation of carbonate minerals represents the main mechanism controlling the groundwater chemistry. Although the aquifer is composed of carbonate and silicate minerals, carbonate dissolution is dominant against silicate weathering, due to the great difference in their weathering rates. Most of the time, sampled groundwater is saturated with respect to calcite, which enables the precipitation of calcite, and undersaturated with respect to dolomite.
(b)
The cation exchange process is well documented in the catchment area of well P-4039, while other observation wells do not show the signs of this process.
(c)
The transformation of organic matter is observed in the catchment area of the observation wells P-1529, PDS-7, PDS-5, and PDS-6. High DOC concentrations in water are followed by high partial pressure of CO2, which is a consequence of flushing organic matter from the soil and unsaturated zone into the aquifer.
(d)
An anthropogenic influence is recognized through high nitrate concentrations in groundwater. The application of synthetic fertilizers and manure in agricultural production is considered the main source of nitrate contamination. However, changes in land use and recent urbanization caused a more significant impact of waste water on nitrate content in the Varaždin aquifer.

Author Contributions

I.K.: investigation, chemical analyses, conceptualization, writing—review and editing; T.M.: funding acquisition, investigation, chemical analyses, conceptualization, data interpretation, writing—original draft preparation; M.Š.M.: chemical analyses; K.M.: chemical analyses and Methods writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Croatian Scientific Foundation (HRZZ) under grant number HRZZ-IP-2016-06-5365 and by Young Researchers Career Development Project—Training of New PhDs—HRZZ and ESF.

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.

References

  1. Stumm, W.; Morgan, J.J. Aquatic Chemistry: Chemical Equilibria and Rates in Natural Waters, 3rd ed.; John Wiley & Sons, INC: New Jersey, NJ, USA, 1995; p. 625. [Google Scholar]
  2. Malcom, R.; Soulsby, C. Hydrogeochemistry of groundwater in coastal wetlands: Implications for coastal conservation in Scotland. Sci. Total. Environ. 2001, 265, 269–280. [Google Scholar] [CrossRef]
  3. Robins, N.S. Groundwater quality in Scotland: Major ion chemistry of the key groundwater bodies. Sci. Total. Environ. 2002, 294, 41–56. [Google Scholar] [CrossRef] [Green Version]
  4. Edmunds, W.M.; Shand, P.; Hart, P.; Ward, R.S. The natural (baseline) quality of groundwater: A UK pilot study. Sci. Total. Environ. 2003, 310, 25–35. [Google Scholar] [CrossRef] [Green Version]
  5. Rosen, M.; Jones, S. Controls on the chemical composition of groundwater from alluvial aquifers in the Wanaka and Wakatipu basins, Central Otago, New Zealand. Hydrogeol. J. 2007, 6, 264–281. [Google Scholar]
  6. Marković, T.; Miko, S.; Kapelj, S.; Buljan, R.; Larva, O.; Peh, Z. Behaviour of metals and nutrients in soils and groundwater of a karst polje. J. Geochem. Explor. 2006, 88, 124–129. [Google Scholar] [CrossRef]
  7. Marković, T.; Brkić, Ž.; Larva, O. Using hydrochemical data and modelling to enhance the knowledge of groundwater flow and quality in an alluvial aquifer of Zagreb, Croatia. Sci. Total. Environ. 2013, 458–460, 508–516. [Google Scholar] [CrossRef] [PubMed]
  8. Marković, T.; Borović, S.; Larva, O. Geochemical characteristics of thermal waters of Hrvatsko zagorje. Geol. Croat. 2015, 68, 67–77. [Google Scholar] [CrossRef]
  9. Marković, T.; Karlović, I.; Perčec Tadić, M.; Larva, O. Application of Stable Water Isotopes to Improve Conceptual Model of Alluvial Aquifer in the Varaždin Area. Water 2020, 12, 379. [Google Scholar] [CrossRef] [Green Version]
  10. Brkić, Ž.; Briški, M.; Marković, T. Use of hydrochemistry and isotopes for improving the knowledge of groundwater flow in a semiconfined aquifer system of the Eastern Slavonia (Croatia). Catena 2016, 142, 153–165. [Google Scholar] [CrossRef]
  11. Pandey, S.K.; Singh, A.K.; Hasnain, S.I. Hydrochemical characteristics of meltwater draining from Pindari glacier, Kumon Himalaya. J. Geol. Soc. India 2001, 57, 519–527. [Google Scholar]
  12. McGinleya, P.M.; Masarika, K.C.; Gotkowitzb, M.B.; Mechenicha, D.J. Impact of anthropogenic geochemical change and aquifer geology on groundwater phosphorus concentrations. Appl. Geochem. 2016, 72, 1–9. [Google Scholar] [CrossRef]
  13. Sajil Kumar, P.J. Hydrogeochemical and multivariate statistical appraisal of pollution sources in the groundwater of the lower Bhavani River basin in Tamil Nadu. GeologyEcolog. Landsc. 2020, 4, 40–51. [Google Scholar] [CrossRef] [Green Version]
  14. Hallburg, G.R.; Keeney, D.R. Nitrate. In Regional Groundwater Quality; Alley, W.M., Ed.; Van Norstrand Reinhold: New York, NY, USA, 1993; pp. 297–322. [Google Scholar]
  15. Spalding, R.F.; Exner, M.E. Occurence of nitrate in groundwater: A review. J. Environ. Qual. 1993, 22, 392–402. [Google Scholar] [CrossRef]
  16. Nolan, B.T. Relating nitrogen sources and aquifer susceptibility of nitrate in shallow ground waters of the United States. Ground Water 2001, 39, 290–299. [Google Scholar] [CrossRef]
  17. Walton, B. Survey of literature relating infant methemoglobinemia due to nitrate contaminated water. Am. J. Public Health 1951, 41, 986–996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Winneberger, J.H.T. Nitrogen, Public Health and the Environment; Ann Arbor Science Publishers Inc.: Ann Arbor, MI, USA, 1982. [Google Scholar]
  19. World Health Organization. Health Hazards from Nitrates in Drinking Water, Report on a WHO Meeting in Copenhagen, Denmark, 5–9 March 1984; Regional Office for Europe: Copenhagen, Denmark, 1985; Volume 3, pp. 49–66. [Google Scholar]
  20. Zaninović, K.; Gajić-Čapka, M.; Perčec Tadić, M. Klimatski Atlas Hrvatske. Climate Atlas of Croatia: 1961–1990, 1971–2000; Državni hidrometeorološki zavod: Zagreb, Croatia, 2008; p. 200. [Google Scholar]
  21. Nimac, I.; Perčec Tadić, M. New 1981–2010 Climatological Normals for Croatia and Comparison to Previous 1961–1990 and 1971–2000 Normals. In Proceedings of the GeoMLA Conference, Belgrade, Serbia, 21–24 June 2016; pp. 79–85. [Google Scholar]
  22. Babić, Ž.; Čakarun, I.; Sokač, A.; Mraz, V. On geological features of quaternary sediments of Drava basin on Croatian territory. Geol. Vjesn. 1978, 30, 43–61. [Google Scholar]
  23. Urumović, K. On quaternary aquifer complex in Varaždin area. Geol. Vjesn. 1971, 24, 183–191. [Google Scholar]
  24. Urumović, K.; Hlevnjak, B.; Prelogović, E.; Mayer, D. Hydrogeological conditions of Varaždin aquifer. Geol. Vjesn. 1990, 43, 149–158. [Google Scholar]
  25. Prelogović, E.; Velić, I. Quaternary tectonic activity in western part of Drava basin. Geol. Vjesn. 1988, 41, 237–253. [Google Scholar]
  26. Gazarek, M.; Crnički, J.; Premur, V.; Kreč, D. Granulometrijski sastav šljunka i pijesaka i teški minerali u pijescima dravskog bazena od Ormoža do Đurđevca. Rudarsko Geološko Naftni Zbornik 1990, 2, 67–73. [Google Scholar]
  27. Crnički, J. Mineralne Sirovine Varaždinske Regije i Njihovo Privredno Značenje. Varaždinski Zbornik 1181–1981, Varaždin, Croatia, 1–3 October 1981; Mohorovičić, A., Ed.; Jugoslavenska akademija znanosti i umjetnosti: Zagreb, Hrvatska, 1983. [Google Scholar]
  28. Mioč, P.; Marković, S. Osnovna Geološka Karta SFRJ 1:100.000. In List Čakovec L33–57 (Basic Geological Map of SFRY, Scale 1:100.000. Sheet Čakovec L33–57—in Croatian); Geological Survey, Zagreb and Geological, Geotechnical and Geophysical Survey: Ljubljana, Slovenia, 1998. [Google Scholar]
  29. Mioč, P.; Marković, S. Osnovna Geološka Karta SFRJ 1:100.000. In Tumač za list Čakovec L33–57 (Basic Geological Map of SFRY, Scale 1:100.000. The Map Book for Sheet Čakovec L33–57—in Croatian); Geological Survey, Zagreb and Geological, Geotechnical and Geophysical Survey: Ljubljana, Slovenia, 1998. [Google Scholar]
  30. Karlović, I.; Marković, T.; Vujnović, T.; Larva, O. Development of a Hydrogeological Conceptual Model of the Varaždin Alluvial Aquifer. Hydrology 2021, 8, 19. [Google Scholar] [CrossRef]
  31. Mandel, S.; Shiftan, Z.L. Groundwater Resources: Investigation and Development; Academic Press: New York, NY, USA, 1981; p. 269. [Google Scholar]
  32. Domenico, P.A.; Schwartz, F.W. Physical and Chemical Hydrogeology; John Willey and Sons: New York, NY, USA, 1990; p. 824. [Google Scholar]
  33. HRN EN ISO 17294-2:2016. Determination of Selected Elements Including Uranium Isotopes. Available online: https://repozitorij.hzn.hr/norm/HRN+EN+ISO+17294-2%3A2016 (accessed on 21 April 2021).
  34. Parkhurst, D.L.; Appelo, C.A.J. Description of Input and Examples for PHREEQC Version 3—A Computer Program for Speciation, Batch-Reaction, One-Dimensional Transport, and Inverse Geochemical Calculations; U.S. Geological Survey Techniques and Methods, Book 6; USA; 2013; Chapter A43; p. 497. Available online: https://pubs.usgs.gov/tm/06/a43/ (accessed on 21 April 2021).
  35. McMahon, P.B.; Chapelle, F.H. Redox processes and water quality of selected principal aquifer systems. Ground Water 2008, 46, 259–271. [Google Scholar] [CrossRef]
  36. OG 125/17. Rulebook on Compliance Parameters, Methods of Analysis, Monitoring and Water Safety Plans for Human Consumption. Available online: https://narodne-novine.nn.hr/clanci/sluzbeni/2017_12_125_2848.html (accessed on 21 April 2021).
  37. Meybeck, M. Global chemical weathering from surficial rocks estimated from river dissolved loads. Am. J. Sci. 1987, 287, 401–428. [Google Scholar] [CrossRef]
  38. Faure, G. Principles and Applications of Geochemistry, 2nd ed.; Prentice Hall: New Jersey, NJ, USA, 1998; p. 600. [Google Scholar]
  39. Mukherjee, A.; Bhattacharya, P.; Shi, F.; Fryar, A.E.; Mukherjee, A.B.; Xie, Z.M.; Jacks, G.; Bundschuh, J. Chemical evolution in the high arsenic groundwater of the Huhhot basin (Inner Mongolia, PR China) and its difference from the western Bengal basin (India). Appl. Geochem. 2009, 24, 1835–1851. [Google Scholar] [CrossRef]
  40. McLean, W.; Jankowski, J.; Lavitt, N. Groundwater quality and sustainability in an alluvial aquifer, Australia. In Groundwater, Past Achievements and Future Challenges, Proceedings of the XXX IAH congress, Cape Town, South Africa, 26 November–1 December 2000; Balkema: Rotterdam, The Nederland, 2000; pp. 567–573. [Google Scholar]
  41. Jankowski, J.; Acworth, R.I.; Shekarforoush, S. Reverse ion exchange in deeply weathered porphyritic dacite fractured aquifer system, Yass, New South Wales, Australia. In Proceedings of the 9th International Symposium on Water-Rock Interaction, Taupo, New Zealand, 30 March–3 April 1998; pp. 243–246. [Google Scholar]
  42. Deutsch, W.J. Groundwater Geochemistry: Fundamentals and Application to Contamination, 1st ed.; CRC Press: Boca Raton, FL, USA, 1997; p. 232. [Google Scholar]
  43. Macpherson, G.L.; Roberts, J.A.; Blair, J.M.; Townsend, M.A.; Fowle, D.A.; Beisner, K.R. Increasing shallow groundwater CO2 and limestone weathering, Konza Prairie, USA. Geochim. Cosmochim. Acta 2008, 72, 5581–5599. [Google Scholar] [CrossRef]
  44. Keeney, D. Sources of Nitrate to Ground Water. Crit. Rev. Environ. Contr. 1986, 32, 1039–1050. [Google Scholar] [CrossRef]
  45. Kaushal, S.S.; Groffman, P.M.; Brand, L.E.; Elliott, E.M.; Shields, C.A.; Kendall, C. Tracking Nonpoint Source Nitrogen Pollution in Human-Impacted Watersheds. Environ. Sci. Technol. 2011, 45, 8225–8232. [Google Scholar] [CrossRef]
  46. Jung, H.; Koh, D.-C.; Kim, Y.S.; Jeen, S.-W.; Lee, J. Stable Isotopes of Water and Nitrate for the Identification of Groundwater Flowpaths: A Review. Water 2020, 12, 138. [Google Scholar] [CrossRef] [Green Version]
  47. Jogun, K.; Pavlek, K.; Belić, T.; Buhin, S.; Malešić, N. Land cover changes in northern Croatia from 1981 to 2011. Hrvat. Geogr. Glas. 2017, 79, 33–59. [Google Scholar] [CrossRef] [Green Version]
  48. Kacaroglu, F.; Gunay, G. Impacts of Human Activities on Groundwater Quality of an Alluvial Aquifer: A Case Study of the Eskisehir plain, Turkey. Hydrogeol. J. 1997, 5, 60–70. [Google Scholar] [CrossRef]
  49. Fernando, W.A.M.; Ilankoon, I.; Syed, T.H.; Yellishetty, M. Challenges and opportunities in the removal of sulphate ions in contaminated mine water: A review. Miner. Eng. 2018, 117, 74–90. [Google Scholar] [CrossRef]
Figure 1. Geographical position of the study area with locations of the groundwater sampling points and photographs of pumping the observation wells at characteristic locations (P-1556, PDS-5, SPV-11, P-4039). The general groundwater flow direction is defined by the head contours for medium water levels measured on 14 October 2013.
Figure 1. Geographical position of the study area with locations of the groundwater sampling points and photographs of pumping the observation wells at characteristic locations (P-1556, PDS-5, SPV-11, P-4039). The general groundwater flow direction is defined by the head contours for medium water levels measured on 14 October 2013.
Water 13 01508 g001
Figure 2. Distribution of ammonia, nitrite and orthophosphate in sampled waters by wells.
Figure 2. Distribution of ammonia, nitrite and orthophosphate in sampled waters by wells.
Water 13 01508 g002
Figure 3. Distribution of dissolved iron and zinc concentrations in sampled waters by well.
Figure 3. Distribution of dissolved iron and zinc concentrations in sampled waters by well.
Water 13 01508 g003
Figure 4. Piper diagram of sampled groundwater.
Figure 4. Piper diagram of sampled groundwater.
Water 13 01508 g004
Figure 5. Bivariate mixing diagram of Na+-normalized Mg2+ versus Na+-normalized Ca2+.
Figure 5. Bivariate mixing diagram of Na+-normalized Mg2+ versus Na+-normalized Ca2+.
Water 13 01508 g005
Figure 6. Molar ratio Mg2+/Ca2+ versus Ca2+.
Figure 6. Molar ratio Mg2+/Ca2+ versus Ca2+.
Water 13 01508 g006
Figure 7. Bivariate diagram of (Ca2+ + Mg2+) - (HCO3 + SO42−) versus Na+ − Cl: (a) all wells; (b) all wells except well P-4039.
Figure 7. Bivariate diagram of (Ca2+ + Mg2+) - (HCO3 + SO42−) versus Na+ − Cl: (a) all wells; (b) all wells except well P-4039.
Water 13 01508 g007
Figure 8. Monthly variation of calcite and dolomite saturation indices and partial pressure of CO2 in groundwater.
Figure 8. Monthly variation of calcite and dolomite saturation indices and partial pressure of CO2 in groundwater.
Water 13 01508 g008
Figure 9. Relationship between concentrations of DOC and logpCO2 in sampled waters.
Figure 9. Relationship between concentrations of DOC and logpCO2 in sampled waters.
Water 13 01508 g009
Figure 10. Na+ vs. Cl in the sampled groundwater.
Figure 10. Na+ vs. Cl in the sampled groundwater.
Water 13 01508 g010
Figure 11. Distribution of sulfate, nitrate and chloride in sampled groundwater.
Figure 11. Distribution of sulfate, nitrate and chloride in sampled groundwater.
Water 13 01508 g011
Table 1. Observation wells information.
Table 1. Observation wells information.
Observation WellElevation
(m a.s.l.)
Depth of the Well
(m)
Depth Interval of the Screen (m)
Private well Hrašćica176.0015.05–15
PDS-5178.3631.013.7–19.7
PDS-6184.0725.011.7–17.7
PDS-7175.7142.529.3–32.3
P-1529187.328.0n.a. 1
P-1556193.0315.6n.a. 1
P-2500167.815.20n.a. 1
P-4039167.768.0n.a. 1
SPV-11177.6940.024.5–35.8
1 information about screen interval not available.
Table 2. (a) The average, minimum and maximum values of the analyzed physicochemical parameters in the groundwater samples. (b) The average, minimum and maximum values of metal concentrations in the groundwater samples.
(a)
(a)
EC
(µS/cm)
T (°C)pHO2 (mg/L)HCO3 (mg/L)PO43−-P (mg/L)NH4+ (mg/L)NO2 (mg/L)Cl (mg/L)SO42− (mg/L)NO3 (mg/L)TN (mg/L)Br (mg/L)Ca2+ (mg/L)Mg2+ (mg/L)Na+ (mg/L)K+ (mg/L)DOC (mg/L)SiO2 (mg/L)
Private well min67310.86.913.1342<0.01<0.01<0.019.321.041.810.6<0.1010118.97.23.90.3011.4
max74713.87.408.54140.220.050.0243.134.091.920.94.013323.819.75.30.8314.4
average71313.17.226.73820.040.020.0119.026.658.414.42.511120.79.24.50.4512.5
P-1529min75510.36.861.5388<0.01<0.01<0.0114.121.040.111.2<0.1010819.913.55.10.1910.1
max81414.57.369.25120.540.130.2437.331.877.117.54.015328.317.512.54.618.4
average78712.67.156.14370.100.050.0423.125.855.714.02.911924.615.16.11.312.7
P-1556min6589.406.930.6381<0.01<0.01<0.015.79.45.32.2<0.1010517.25.24.20.439.2
max87716.07.457.85120.760.070.0432.934.835.66.7<0.1015422.521.310.15.415.1
average74313.17.154.94420.100.030.0214.723.817.44.4<0.1012319.69.95.81.211.8
SPV-11min49012.07.031.4249<0.01<0.01<0.015.324.48.02.2<0.1073.915.52.40.940.1310.3
max49614.47.567.72780.140.070.0224.233.218.95.00<0.1082.620.613.13.80.6016.1
average49412.37.432.22670.050.030.019.127.311.22.7<0.1075.816.83.91.60.5912.7
P-2500min69610.87.010.9325<0.01<0.01<0.0118.715.027.511.3<0.1099.419.011.51.00.5210.5
max80217.97.559.24050.640.160.0465.646.813731.35.813822.748.12.91.215.6
average73713.77.296.53610.070.040.0134.631.366.415.93.611220.622.31.70.7112.4
P-4039min76611.07.010.2238<0.01<0.01<0.0175.86.6<0.10<1.0<0.1082.610.747.11.50.699.6
max109114.77.714.24100.350.140.0627940.221.92.44.214726.51573.72.422.0
average97513.27.401.53240.080.040.0217028.25.11.42.810721.981.52.81.313.1
PDS-5min66111.86.913.6322<0.01<0.01<0.018.813.042.512.1<0.1010018.75.311.40.2510.7
max69413.97.459.94560.820.070.0222.467.521047.73.014729.915.988.80.5219.7
average68312.67.288.33870.120.030.0114.727.383.020.01.511322.76.231.30.3513.4
PDS-6min70811.46.927.3260<0.01<0.01<0.019.011.038.511.8<0.1086.317.97.01.50.3010.8
max74413.07.459.14560.220.040.0224.030.012327.94.015021.79.73.02.529.9
average73012.47.178.73810.050.030.0114.722.265.415.42.512019.68.02.20.6114.4
PDS-7min73011.26.996.6280<0.01<0.01<0.0111.419.052.512.6<0.1011219.54.30.400.2510.9
max77713.07.7610.94351.10.040.0468.035.918041.04.015125.415.08.92.623.5
average75812.37.328.93740.140.020.0220.427.796.723.23.012321.46.71.90.6415.1
(b)
(b)
As
(µg/L)
Cd
(µg/L)
Cr (µg/L)Cu (µg/L)Fe (µg/L)Li
(µg/L)
Mn (µg/L)Mo (µg/L)Ni (µg/L)Pb (µg/L)Sr
(µg/L)
Zn (µg/L)
Private well min0.11<0.010.280.491.10.680.080.51<0.020.0613513.2
max0.440.060.577.319.83.91.70.991.10.4621188.4
average0.180.020.431.406.62.10.540.770.510.2018140.2
P-1529min0.07<0.010.250.552.11.20.060.360.220.1117128.0
max0.360.100.896.121.44.59.80.641.11.726778.0
average0.160.050.452.37.32.72.70.510.520.5323545.9
P-1556min0.34<0.01<0.020.302081.20.290.330.050.10242104
max1.10.050.263.54833.62370.713.20.603712134
average0.810.030.080.963492.71580.521.10.27303557
SPV-11min0.10<0.010.060.071.7<0.01<0.050.530.120.052433.5
max0.390.060.461.341.81.10.980.955.30.4132713.5
average0.190.020.180.579.30.690.420.780.690.172966.2
P-2500min0.09<0.010.030.151.80.140.270.120.110.051267.6
max0.300.080.422.51842.18.80.344.60.7621276.3
average0.150.020.131.019.41.31.60.210.610.1917218.6
P-4039min0.120.01<0.020.131350.9511.10.280.060.11146427
max0.850.090.444.122652.930.50.9446.74.62515600
average0.270.030.111.17441.819.80.522.61.22052873
PDS-5min0.11<0.010.370.111.50.100.290.19<0.020.06143205
max0.390.050.636.734.02.91.60.341.70.60223543
average0.170.020.510.698.21.60.680.260.320.19188330
PDS-6min0.08<0.010.270.152.00.890.250.16<0.020.41192265
max0.350.110.663.29.64.24.50.710.841.7309627
average0.140.030.470.884.92.61.10.250.340.76260466
PDS-7min0.05<0.010.150.0210.20.312.20.290.040.01116169
max0.330.090.821.659.32.716.50.581.31.52685665
average0.110.040.490.7234.11.67.80.440.440.482272169
Table 3. The calculated redox conditions of groundwater.
Table 3. The calculated redox conditions of groundwater.
Observation WellGeneral Redox CategoryRedox Process
Private well HrašćicaOxicO2
P-1529OxicO2
P-1556Mixed (oxic-anoxic)O2-Fe(III)/SO4 or O2-Mn(IV)
SPV-11OxicO2
P-2500OxicO2
P-4039Anoxic or Mixed (oxic-anoxic)NO3-Fe(III)/SO4 or O2-Fe(III)/SO4
PDS-5OxicO2
PDS-6OxicO2
PDS-7OxicO2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Karlović, I.; Marković, T.; Šparica Miko, M.; Maldini, K. Geochemical Characteristics of Alluvial Aquifer in the Varaždin Region. Water 2021, 13, 1508. https://doi.org/10.3390/w13111508

AMA Style

Karlović I, Marković T, Šparica Miko M, Maldini K. Geochemical Characteristics of Alluvial Aquifer in the Varaždin Region. Water. 2021; 13(11):1508. https://doi.org/10.3390/w13111508

Chicago/Turabian Style

Karlović, Igor, Tamara Marković, Martina Šparica Miko, and Krešimir Maldini. 2021. "Geochemical Characteristics of Alluvial Aquifer in the Varaždin Region" Water 13, no. 11: 1508. https://doi.org/10.3390/w13111508

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop