Next Article in Journal
Rainwater Harvesting for Sustainable Agriculture in High Water-Poor Areas in the West Bank, Palestine
Next Article in Special Issue
Can Soil Hydraulic Parameters be Estimated from the Stable Isotope Composition of Pore Water from a Single Soil Profile?
Previous Article in Journal
Designing Human Rights for Duty Bearers: Making the Human Rights to Water and Sanitation Part of Everyday Practice at the Local Government Level
Previous Article in Special Issue
Application of a Self-Organizing Map of Isotopic and Chemical Data for the Identification of Groundwater Recharge Sources in Nasunogahara Alluvial Fan, Japan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Stable Water Isotopes to Improve Conceptual Model of Alluvial Aquifer in the Varaždin Area

1
Croatian Geological Survey, 10 000 Zagreb, Croatia
2
Meteorological and Hydrological Service, 10 000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Water 2020, 12(2), 379; https://doi.org/10.3390/w12020379
Submission received: 9 December 2019 / Revised: 27 January 2020 / Accepted: 28 January 2020 / Published: 30 January 2020
(This article belongs to the Special Issue Use of Water Stable Isotopes in Hydrological Process)

Abstract

:
To understand groundwater flow and geochemical processes within an aquifer, it is necessary to set up a conceptual model of the aquifer. To accomplish this, different methods are used, and one of them is an isotopic technique. The study area is located in the Varaždin area (NW Croatia). The aquifer represents the main source of potable water for the town of Varaždin and the surrounding settlements. The conceptual model of the alluvial aquifer has to be set up prior to creating a groundwater flow and transport model. Measurements of ratios δ18Ο and δ2H in ground- and surface waters and precipitation samples were carried out. The relationship between ratios δ18Ο, δ2H, and d-excess for local precipitation in the study area showed that precipitation originates from the Atlantic air masses, although during the colder periods of the year, influence of the Mediterranean air masses was not negligible. The monitored period was warmer and wetter than average. Evaporation was observed at all monitored surface waters, but the largest rate was at the location of a gravel pit in Šijanec. The isotopic composition of the precipitation and groundwater showed a good correlation due to the isotopic homogenization of groundwater along the flow path.

1. Introduction

To understand groundwater flow and geochemical processes within an aquifer, it is necessary to set up a conceptual model of it. Different methods are used, and one of them is an isotopic technique, which is often used successfully to help elucidate hydrological studies [1]. Knowledge about the isotopic ratios of oxygen (δ18Ο) and hydrogen (δ2H) in atmospheric precipitation and groundwater is important for hydrological, hydrogeological, climatological, and meteorological applications [1,2,3,4,5,6,7,8] because it can provide information on the mean recharge elevation of the aquifer, the mean residence time, water–rock interactions, etc.
The study area is located in the Varaždin area (NW Croatia). The aquifer represents the main source of potable water for the town of Varaždin and the surrounding settlements. The favorable climate, topography, and available groundwater have insured intensive agricultural practices involving the application of large amounts of synthetic fertilizers and manure that have subsequently led to high nitrate concentrations in the Varaždin aquifer. High concentrations of nitrate have caused the shutting down of the Varaždin pumping site. To determine the behavior of nitrates, a groundwater flow and transport model will be used. A conceptual model of the alluvial aquifer, therefore, has to be set up, and to improve this model, measurements of δ18Ο and δ2H in the ground- and surface waters and precipitation samples were carried out.
The research is still ongoing and, in this paper, only the results of two-year measurements are elaborated in detail. The main goal of this paper is to provide an overview of the spatial and temporal variability in δ18Ο and δ2H values in precipitation and in the ground and surface waters, and this information will be used to determine recharge areas.

2. Geological, Hydrogeological, and Climatic Setting of the Study Area

2.1. Geological and Hydrogeological Setting of the Study Area

The study area of the Varaždin aquifer system is located in NW Croatia, upstream of the town of Varaždin in the valley of the Drava River and it covers an area of approximately 200 km2 (Figure 1). The Varaždin alluvial aquifer is composed of sediments of the Quaternary age, deposited during the Pleistocene and Holocene eras as a result of accumulation processes of the Drava River [9]. The alluvial deposits consist primarily of gravel and sand with occasional lenses and interbeds of silt and clay (Figure 2). The thickness of the aquifer varies from less than 5 m in the NW part to about 65 m in the SE part of the study area (Figure 2, cross-section A–A’). The Varaždin aquifer is an unconfined aquifer, and the groundwater is in direct contact with the surface water: the Drava River and the Plitvica stream, the derivation channel, the accumulation lake Varaždin (Varaždin Lake), and the gravel pit Šijanec (Figure 2). The direction of groundwater flow is generally NW–SE and is parallel to the direction of the Drava River flow. The covering layer of the aquifer is not continuously developed (Figure 2), which makes the aquifer vulnerable to contamination from the surface. Favorable hydrogeological conditions enabled the development of two pumping sites—Varaždin and Vinokovšćak (Figure 2). High concentrations of nitrate have caused the shutting down of the pumping site Varaždin, but Vinokovšćak is still in operation.

2.2. Climatic Setting of the Study Area

The Varaždin meteorological station (46.3° N, 16.137° E, 167 m a.s.l.) and the surroundings have a characteristic precipitation regime with more precipitation during the summer [11]. Because of these features, the local climate is categorized in the Cfb group of the Köppen–Geiger classification system, which is known as “warm-temperate climate” or “marine west coast climate.” The study area is the former [12].
The coldest month is January, with average temperatures of 0.0 °C, while July and August are the warmest months with average temperatures of 20.9 and 20.1 °C, respectively (Table 1). Mean annual temperature is 10.6 °C. The inter-annual variability is the largest in February and January (standard deviation (sd) in Table 1), meaning that average differences from the mean are largest at the end of the winter.
According to the data from the last climate normal period (1981–2010), the lowest average precipitation amounts were during the cold part of the year, and in January the mean precipitation was 38.7 mm (Table 2). From June to September, the precipitation amounts were on average larger than 80 mm, with the highest amount in September (98.3 mm) and the second-highest in June (96.1 mm). The average annual precipitation was 832 mm. Inter-annual variability was largest in January (coefficient of variation (cv) = 0.81).
Based on the long-term data records for 1951–2018, the existence of seasonal trends was tested by a Mann–Kendall test at the 0.05 significance level, and Sen’s slope was calculated to determine the trend value [13]. There is significant warming in all seasons, with Sen’s slope ranging from 1.8 °C/100 years in autumn to 3.7 °C/100 years in summer (Table 3).
There are no statistically significant changes in seasonal monthly precipitation in the 1951–2019 period (Table 4), even though there is a slight indication of drying in summer and in spring, while autumn and winter show a slight tendency of becoming wetter.

3. Materials and Methods

3.1. Water Sampling

Monthly composite precipitation was sampled in the Miko family courtyard in the village Hrašćica (46.3° N, 16.292° E; 177 m a.s.l.) in the period from June 2017 until June 2019. In the field, sample was poured into a 1 L plastic bottle with a tight-fitting cap.
In addition, ground- and surface water sampling was conducted on a monthly basis in the period from June 2017 to June 2019 for chemical and isotopic analyses. Ten observation wells (nine of them are located in the recharge area of the Varaždin pumping site and one in the recharge area of the Vinokovšćak pumping site) and four surface waters (Drava River, the Plitvica stream, Varaždin Lake, which is an accumulation, and a gravel pit in Šijanec) were sampled. The water depth was measured prior to pumping. The observation wells were sampled after stabilization of the parameters EC (electrical conductivity), T (water temperature), pH, and O2 (dissolved oxygen content). Depths of the observation wells are given in Table 5. Samples were poured into a 50 mL plastic bottle with a tight-fitting cap. All samples were measured in the laboratory immediately upon returning from the field.
Meteorological parameters (precipitation and air temperature) used here are from the main Varaždin meteorological station, located in the vicinity of the installed rain gauge.

3.2. Stable Isotope Analyses

The δ18Ο and δ2H were determined using Picarro L2130i (Santa, Clara, USA) in the Hydrochemical Laboratory of the Croatian Geological Survey. The instrument uses CRDS (Cavity Ring-Down Spectroscopy) technology [14]. All measurements were checked with Picarro’s standards (Depleted −29.6 ± 0.2 δ18Ο; −235 ± 1.8 δ2H; Mid −20.6 ± 0.2 δ18Ο; −159 ± 1.3 δ2H; Zero 0.3 ± 0.2 δ18Ο; 1.8 ± 0.9 δ2H), which were checked periodically against the International Atomic Energy Agency (IAEA) standards: Vienna Standard Mean Ocean Water 2 (VSMOW2) and Standard Light Antarctic Precipitation 2 (SLAP2). Measurement precision was ± 0.3 o/oo for δ18Ο and ± 1 o/oo for δ2H.
It is generally known that all isotopic results are expressed as per the international measurement standard, VSMOW2 [15,16].
A global relationship between δ2H and δ18Ο has been observed [16] and called the Global Meteoric Water Line (δ2H = 8·δ18Ο + 10).
The deuterium excess (d-excess [17]) was calculated for each sample as follows:
d-excess (o/oo) = δ2H – 8 δ18Ο.
In 2003, researchers verified this value by using global maps derived by interpolation from more than 340 stations [18] and, in 2005, this was reconfirmed using IAEA-GNIP datasets at 410 stations [19]. Higher values of d-excess are caused by intense evaporation of seawater in conditions of moisture deficit [20].
The local meteoric water line (LMWL) has been calculated using three types of linear regression analysis, two of which are recommended by the IAEA [21]: ordinary least squares regression (OLSR) and reduced major axis (RMA) regression. The new one takes into account the amount of precipitation, using precipitation weighted least squares regression (PWLSR) [22].

4. Results and Discussion

4.1. Precipitation and Temperature

The climatological conditions from June 2017 to June 2019, the period of the isotope sampling, were compared to a climate normal 1981–2010 (Figure 3).
The summer of 2017 was drier and warmer from the average of 1981–2010. Autumn 2017 was mostly wetter than normal, with September 2017 being the wettest month of the period, and it was also colder than average in September. This period, with larger precipitation than average, continued until May 2018, and it was also warmer, except in February and March 2018. From June 2018 to March 2019, precipitation was mostly around or smaller than average, and it was warmer. May 2019 was the second wettest month during this period, and it was also colder than average. June 2019 was again warmer than average.

4.2. Stable Isotopes in Precipitation

The mean stable isotope δ18Ο and δ2H values and the associated d-excess are shown in Table 6.
The stable isotope δ18Ο values varied from −14.91 to −4.5 o/oo, and δ2H varied from −108.1 to −25.1 o/oo (Figure 4, Table S1). The lowest δ18Ο and δ2H values were observed in winter and highest were observed in summer. The d-excess varied from 4.1 to 12.9 o/oo. The d-excess shows the influence of the Atlantic air masses. Nevertheless, the influence of the Mediterranean air masses in the study area was observed during the autumn and winter months (Figure 4). The Mediterranean air masses (precipitation) are characterized by a higher d-excess than the Atlantic air masses [20]. This was observed in References [23,24] in the continental part. Atypical climatological conditions during the observed period had influenced variations of monthly isotopic composition. A sudden change in the air temperature and/or precipitation amount during the season influenced the variation of the monthly isotopic composition of the rain. For example, May 2019 was colder and wetter than average (even than May 2018), and δ18Ο and δ2H values were automatically more negative. In addition, the lowest δ18Ο and δ2H values were measured in the coldest month, which was February 2018 (Figure 3 and Figure 4). It was observed that the isotopic composition of the precipitation in the study area reflects climatological conditions well.
The measured stable isotope δ18Ο and δ2H values were weighted by the amount of precipitation at the Varaždin meteorological station for the observed period. However, there were no large differences between the measured stable isotope δ18Ο and δ2H values and weighted by the amount of precipitation. Because of this, they are not discussed here.
The calculated LMWL for the period from June 2017 to June 2019 is:
OLSR δ2H = (7.54 ± 0.12) δ18Ο + (5.00 ± 1.00), n = 23
RMA δ2H = (7.56 ± 0.11) δ18Ο + (5.17 ± 1.04), n = 23
PWLSR δ2H = (7.55 ± 0.13) δ18Ο + (4.85 ± 1.13), n = 23.
It was observed that all three methods yielded a very similar slope value and axis intercept (b value) of the LMWL, which was supported by very similar measured and weighted values.
In addition, calculated data were compared with data published in Reference [25]. There is a difference between these two slopes values for 0.09 o/oo, and it can be concluded that the OLSR values calculated from the measured data and the published meteoric water line of the study area are not different in terms of slope. However, there is a large difference between these two lines in the axis intercept values for 2.6 o/oo, and the published LMWL is slightly below the measured one (Figure 5). There are several reasons for that: shorter monitored period in our research; very untypical and extreme climatological conditions during our monitored period; and different measurement techniques (our samples were measured using CRDS technology and published were measured using Isotope Ratio Mass Spectrometry (IRMS) technology).

4.3. Stable Isotopes in Ground and Surface Waters

The minimum, maximum, and average isotopic composition of the surface- and groundwaters are given in Table 7 and Table 8 together with their average d-excess.
The measured δ18Ο values in the groundwater varied from −11.47 to −8.26 o/oo, and the δ2H values varied from −81.7 to −58.5 o/oo (Table 7). The measured δ18Ο values in the surface water varied from −12.12 to −3.36 o/oo, and the δ2H values varied from −81.7 to −38.8 o/oo (Table 8).
The correlation between the δ18Ο and δ2H measured values of the groundwater is shown in Figure 5 and indicates that this relationship has a slope of 7.14. Using a Student’s t-test according to Reference [26], a good relationship between groundwater and precipitation was observed. Generally, an isotope relationship between δ18Ο and δ2H with a slope of about 8 is normally observed for precipitation [16]. Since the relationship between the isotopic composition of precipitation and groundwater is good, it can be concluded that groundwater is recharged by precipitation. Values that are slightly more negative were measured in the SPV-11 well and the private well, while at observation wells PDS-5, PDS-6, PDS-7, and P-1529, values are almost identical (Table 7). The highest δ18Ο and δ2H values were measured at observation wells P-4039 and P-2500 (Table 7). The calculated average d-excess values varied from 9.6 to 10.7 o/oo, indicating the influence of recharge by precipitation with signatures of the Atlantic air masses and good homogenization of groundwater along the flow path. This was observed at SPV-11, PDS-6, PDS-7, the private well, P-1530, and P-1529. However, depending on hydrodynamic conditions (low/high water levels), the vicinity of the river or lake, and the depth of the observation well, it was observed that wells, especially the shallower ones and/or those closer to the river and lake, showed high variation in d-excess values. These values were higher than 11 o/oo, indicating recharge by surface waters and faster recharge by precipitation. This was observed at P-1556, PDS-5, P-2500, and P-4039.
The measured δ18Ο and δ2H values of the surface waters distributed around the LMWL shown in Figure 6 indicate a relationship between δ18Ο and δ2H for surface waters with slopes of 5.73 at Varaždin Lake, 6.32 at Drava River, and 4.77 at the gravel pit in Šijanec, indicating an influence of evaporation. A slope from 4 to 6 is attributed to waters with a significant rate of evaporation relative to the input [16]. It was observed that the evaporation process was strongest at the location gravel pit in Šijanec (Figure 7). Nevertheless, the gravel pit was used for fish farming. Because of this activity (resulting in an extra nutrient load due to fish feeding), a low water level, a high load of nutrients, high temperatures, and algae bloom occurred every summer, which had a significant influence on the isotopic and chemical features of this water. The winter-measured values of the δ2Hand δ18O of the Drava River and the summer-measured values of Varaždin Lake are above the LMWL (Figure 6). For the Drava River, this can be explained by the fact that the larger part of the recharge area of the Drava River is situated far upstream of the study area and is under the influence of a different climate, and the influence of the recharge area in the study area is small. Varaždin Lake is recharged by the Drava River, especially in the late winter and springtime when isotopic values are more negative in the river. Since the lake has a high volume, turnover in the lake takes some time. In addition, the Plitvica stream, like the Drava River, has its recharge area in a mountain area where the climate is different, and because of that, more negative values are measured in the wintertime. During the late spring/summer period, the discharge of the stream is low. In the watercourse of the stream, small connected ponds are formed where evapotranspiration is present and, because of that, some values are below the LMWL.
To connect the measured values, a simplified statistical correlation method was used, and results are shown in the correlation matrix in Table 9.
No statistical connection was observed between either the groundwater, the waters of Drava River, Varaždin Lake, or the Plitvica stream with water from the gravel pit in Šijanec. This is partly because there was no observation immediately downstream from the gravel pit. Moreover, it has a very small volume, and its influence is thus limited to its immediate surroundings. Furthermore, the gravel pit is not connected to the river, stream, or lake; consequently, the isotopic composition differs (Figure 1). In addition, a very weak correlation was observed between the waters from the P-2500 and P-4039 wells, which are in the vicinity of the Plitvica stream. This weak correlation is attributed to the drainage roll of the Plitvica stream in this part of the aquifer. A higher correlation was observed between both Varaždin Lake and the Drava River waters and the groundwater of the observation wells. A high correlation between observation wells on the right side of the intake/drain channels indicates a homogenization of the groundwater source (a mixture of the precipitation, river, and lake waters). The left side is different because of the influence of Varaždin Lake. Namely, the Drava River flows from the lake and has the same isotopic composition as the lake water (Figure 1). The river recharges the aquifer as a consequence of groundwater abstraction at the Vinokovšćak pumping site (Figure 1), and the influence of the local precipitation is minor.

5. Conclusions

Although δ18Ο and δ2H values of ground- and surface waters and precipitation were measured for only 24 months, the following conclusions are proposed:
(a)
Meteorological conditions during the observed period were quite different from the average. The summer of 2017 was drier and warmer, while the autumn was wet and cold. This period had more precipitation than average, and this continued until May 2018. It was also warmer, except in February and March 2018. From June 2018 to March 2019, precipitation was mostly around or smaller than average, and it was warmer. May 2019 was the second wettest month during this period and was also colder than average. June 2019 was again warmer than average.
(b)
The isotopic composition of the local precipitation in the study area varied according to the variable meteorological conditions and relationship between δ18Ο and δ2H, and the d-excess showed that precipitation originated from the Atlantic air masses, although during the colder parts of the year, the influence of the Mediterranean air masses was not negligible.
(c)
Evaporation was observed at all surface waters, but the largest rate was at the location of the gravel pit in Šijanec. Moreover, it has a very small volume; therefore, its influence is limited to the immediate surroundings of the aquifer.
(d)
The isotopic composition of precipitation and groundwater showed a good correlation due to the isotopic homogenization of groundwater along the flow path.
(e)
The isotopic composition of the groundwater source on the right side of the intake/drain channels indicates a homogenization of the groundwater source (a mixture of precipitation, river, and lake waters), whereas the left side is different because of the higher influence of Varaždin Lake.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/12/2/379/s1, Table S1: Determined isotopic composition of the precipitation by month for monitored period.

Author Contributions

T.M. contributed to conceptualization, isotopic data interpretation, and writing. I.K. contributed to the hydrogeology and geology description of the study area and graphical editing. M.P.T. contributed to statistical analysis and visualization of climatological data. O.L. contributed by polishing the article. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Croatian Scientific Foundation (HRZZ), grant number HRZZ-IP-2016-06-5356.

Acknowledgments

The authors would like to thank the Miko family for taking care of the rain gauge, as well as the Croatian Meteorological and Hydrological Service for providing climatological data for the Varaždin meteorological station.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Dotsika, E.; Lykoudis, S.; Poutoukis, D. Spatial distribution of the isotopic composition of precipitation and spring water in Greece. Global Planet. Change 2010, 71, 141–149. [Google Scholar] [CrossRef]
  2. Rozanski, K.; Araguas-Araguas, L.; Gonfiantini, R. Isotopic patterns in modern global precipitation. Continental Isotopic Indicators of Climate. Am. Geophys. Union Monogr. 1993, 78, 1–36. [Google Scholar]
  3. Vreča, P.; Krajcar Bronić, I.; Horvatinčić, N.; Barešić, J. Isotopic characteristics of precipitation in Slovenia and Croatia: Comparison of continental and maritime stations. J. Hydrol. 2006, 330, 457–469. [Google Scholar] [CrossRef]
  4. Kanduč, T.; Mori, N.; Kocman, D.; Stibilj, V.; Grassa, F. Hydrogeochemistry of Alpine springs from North Slovenia: Insights from stable isotopes. Chem. Geol. 2012, 300/301, 40–54. [Google Scholar] [CrossRef]
  5. 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. of the Tot. Environ. 2013, 458–460, 508–516. [Google Scholar] [CrossRef] [PubMed]
  6. Lukač Reberski, J.; Marković, T.; Nakić, Z. Definition of the river Gacka springs subcatchment areas on the basis of hydrogeological parameters. Geologia Croatica. 2013, 66, 39–53. [Google Scholar] [CrossRef]
  7. Gross, E.; Andrews, S.; Bergamaschi, B.; Downing, B.; Holleman, R.; Burdick, S.; Durand, D. The Use of Stable Isotope-Based Water Age to Evaluate a Hydrodynamic Model. Water 2019, 11, 2207. [Google Scholar] [CrossRef] [Green Version]
  8. Cervi, F.; Dadomo, A.; Martinelli, G. The Analysis of Short-Term Dataset of Water Stable Isotopes Provides Information on Hydrological Processes Occurring in Large Catchments from the Northern Italian Apennines. Water. 2019, 11, 1360. [Google Scholar] [CrossRef] [Green Version]
  9. Prelogović, E.; Velić, I. Kvartarna tektonska aktivnost u zapadnom dijelu Dravske potoline. Geol. vjesnik 1988, 41, 237–253. [Google Scholar]
  10. Larva, O. Ranjivost vodonosnika na priljevnom području varaždinskih crpilišta (Aquifer vulnerability at catchment area of Varaždin pumping sites—in Croatian). Ph.D. Thesis, University of Zagreb, Zagreb, Croatia, December 2008. [Google Scholar]
  11. 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]
  12. Nimac, I.; Perčec Tadić, M. New 1981–2010 Climatological Normals for Croatia and Comparison to Previous 1961–1990 and 1971–2000 Normals. Proceedings from GeoMLA Conference; University of Belgrade—Faculty of Civil Engineering: Belgrade, Serbia, 2016; pp. 79–85. [Google Scholar]
  13. Wilks, D.S. Statistical Methods in the Atmospheric Sciences, 2nd ed.; Academic Press: London, UK, 2006; p. 704. [Google Scholar]
  14. Available online: http://www.picarro.com/technology/cavity_ring_down_spectroscopy (accessed on 3 June 2019).
  15. Clark, I.D.; Fritz, P. Environmental Isotopes in Hydrogeology; CRC Press: Boca Raton, FL, USA, 2013; ISBN 978-1-4822-4291-1. [Google Scholar]
  16. Craig, H. Isotopic Variations in Meteoric Waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef] [PubMed]
  17. Dansgaard, W. Stable isotopes in precipitation. Tellus 1964, 16, 436–468. [Google Scholar] [CrossRef]
  18. Bowen, G.J.; Revenaugh, J. Interpolating the isotopic composition of modern meteoric precipitation. Water Resour. Res. 2003, 39, 10. [Google Scholar] [CrossRef]
  19. Gourcy, L.L.; Groening, M.; Aggarwal, P.K. Stable oxygen and hydrogen isotopes in precipitation. Isotopes in the Water Cycle: Past. Present and Future of Developing Science. Dordrecht, The Netherlands, 2005 Aggarwal, P.K.; Gat, J.R. and Froehlich, K.F.O., Eds.; Springer, 39–51. Available online: https://www.springer.com/gp/book/9781402030109 (accessed on 30 January 2020). [CrossRef]
  20. Gat, J.R.; Carmi, I. Evolution of the isotopic composition of the atmospheric water in the Mediterranean Sea area. J. Geophys. Res. 1970, 75, 3039–3048. [Google Scholar]
  21. IAEA. Statistical treatment of data on environmental isotopes in precipitation. Technical Report Series 331. International Atomic Energy Agency, Vienna. 1992. Available online: https://www.iaea.org/publications/1435/statistical-treatment-of-data-on-environmental-isotopes-in-precipitation (accessed on 30 January 2020).
  22. Hughes, C.E.; Crawford, J. 2012: A new precipitation weighted method for determining the meteoric water line for hydrological applications demonstrated using Australian and global GNIP data. J. Hydrol. 2012, 464–465, 344–351. [Google Scholar] [CrossRef]
  23. Vreča, P.; Krajcar Bronić, I.; Leis, A. Demšar. M. Isotopic composition of precipitation at the station Ljubljana (Reaktor). Slovenia—period 2007–2010. Geologija 2014, 57, 217–230. [Google Scholar]
  24. Bottyán, E.; Czuppon, G.; Weidinger, T.; Haszpra, L.; Kármán, K. Moisture source diagnostics and isotope characteristics for precipitation in east Hungary: implications for their relationship. Hydrol. Sci. J. 2017, 62, 2049–2060. [Google Scholar]
  25. Hunjak, T.; Lutz, H.O.; Roller-Lutz, Z. Stable isotope composition of the meteoric precipitation in Croatia. Isot. Environ. Health Stud. 2013, 49, 336–345. [Google Scholar] [CrossRef] [PubMed]
  26. Boschetti, T.; Cifuentes, J.; Iacumin, P.; Selmo, E. Local Meteoric Water Line of Northern Chile (18 S–30 S): An Application of Error-in-Variables Regression to the Oxygen and Hydrogen Stable Isotope Ratio of Precipitation. Water 2019, 11, 791. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geographical position of the study area with locations of the sampling points; the transects A–A’ and B–B’ correspond to the representative hydrogeological cross-sections shown in Figure 2. The general groundwater flow direction is defined by the water heads and stable isotopes in the study area.
Figure 1. Geographical position of the study area with locations of the sampling points; the transects A–A’ and B–B’ correspond to the representative hydrogeological cross-sections shown in Figure 2. The general groundwater flow direction is defined by the water heads and stable isotopes in the study area.
Water 12 00379 g001
Figure 2. Schematic hydrogeological cross-sections across the study area (modified according to Reference [10]).
Figure 2. Schematic hydrogeological cross-sections across the study area (modified according to Reference [10]).
Water 12 00379 g002
Figure 3. Monthly anomaly of precipitation (a) and air temperature (b) in the period from June 2017 to June 2019 based on 1981–2010 climate normal.
Figure 3. Monthly anomaly of precipitation (a) and air temperature (b) in the period from June 2017 to June 2019 based on 1981–2010 climate normal.
Water 12 00379 g003
Figure 4. Monthly variations of (a) δ18Ο (o/oo) values, (b) δ2H (o/oo) values, and (c) d-excess (o/oo) in precipitation at the rain gauge in Hrašćica.
Figure 4. Monthly variations of (a) δ18Ο (o/oo) values, (b) δ2H (o/oo) values, and (c) d-excess (o/oo) in precipitation at the rain gauge in Hrašćica.
Water 12 00379 g004
Figure 5. Distribution of groundwater δ18Ο (o/oo) and δ2H (o/oo) values around the local meteoric water line (LMWL).
Figure 5. Distribution of groundwater δ18Ο (o/oo) and δ2H (o/oo) values around the local meteoric water line (LMWL).
Water 12 00379 g005
Figure 6. Distribution of surface waters δ18Ο (o/oo) and δ2H (o/oo) values around the LMWL.
Figure 6. Distribution of surface waters δ18Ο (o/oo) and δ2H (o/oo) values around the LMWL.
Water 12 00379 g006
Figure 7. Distribution of δ18Ο (o/oo) values in surface waters over the monitored time.
Figure 7. Distribution of δ18Ο (o/oo) values in surface waters over the monitored time.
Water 12 00379 g007
Table 1. Monthly and annual average air temperature in Varaždin. Statistics in rows refer to maximum, mean, and minimum average monthly and annual temperature in the 1981–2010 period, and the standard deviation (sd).
Table 1. Monthly and annual average air temperature in Varaždin. Statistics in rows refer to maximum, mean, and minimum average monthly and annual temperature in the 1981–2010 period, and the standard deviation (sd).
Values010203040506070809101112Annual
max (°C)5.86.19.914.018.723.822.824.518.013.89.54.712.1
mean (°C)0.01.66.110.916.019.120.920.115.610.65.31.110.6
sd (°C)2.53.22.11.41.41.41.21.51.31.52.31.80.8
min (°C)−6−4.50.48.012.216.718.518.012.58.10.8−2.79.0
Table 2. The monthly and annual precipitation sum in Varaždin. Statistics in rows refer to maximum, mean, and minimum monthly and annual precipitation in the 1981–2010 period, the standard deviation (sd), and the coefficient of variation (cv).
Table 2. The monthly and annual precipitation sum in Varaždin. Statistics in rows refer to maximum, mean, and minimum monthly and annual precipitation in the 1981–2010 period, the standard deviation (sd), and the coefficient of variation (cv).
Values010203040506070809101112Annual
max (mm)145.4124.6100.6121.3144.2199.9183.7211.7186.1202181.5169.91200.3
mean (mm)38.740.854.961.067.996.181.987.298.378.068.059.1832.0
sd (mm)31.426.726.634.131.140.341.053.946.348.042.035.8131.4
cv0.810.650.480.560.460.420.50.620.470.620.620.610.16
min (mm)3.30.32.14.923.832.115.34.825.62.419.617.1559.7
Table 3. Trend analysis of long-term temperature data for the period 1951–2019. A Mann–Kendall p-value of <0.05 indicates a significant trend. Sen’s slope is a value of that trend. A Sen’s slope p-value of <0.5 indicates a significant value of a trend.
Table 3. Trend analysis of long-term temperature data for the period 1951–2019. A Mann–Kendall p-value of <0.05 indicates a significant trend. Sen’s slope is a value of that trend. A Sen’s slope p-value of <0.5 indicates a significant value of a trend.
SeasonMann–Kendall p-ValueSen's Slope (°C/100y)Sen's Slope p-Value
Autumn1.26 x 10-21.83 x 10-21.26 x 10-2
Spring1.49 x 10-63.30 x 10-21.49 x 10-6
Summer2.89 x 10-93.67 x 10-22.89 x 10-9
Winter1.75 x 10-33.10 x 10-21.75 x 10-3
Table 4. Trend analysis of long-term precipitation data for the period 1951–2019. Parameters are the same as in Table 3.
Table 4. Trend analysis of long-term precipitation data for the period 1951–2019. Parameters are the same as in Table 3.
SeasonMann–Kendall p-ValueSen's SlopeSen's Slope p-Value
Autumn1.64 x 10-18.00 x 10-11.64 x 10-1
Spring4.88 x 10-1−2.66 x 10-14.88 x 10-1
Summer8.46 x 10-2−9.72 x 10-18.46 x 10-2
Winter9.46 x 10-11.36 x 10-29.46 x 10-1
Table 5. Depths of the observation wells.
Table 5. Depths of the observation wells.
Observation WellElevation (m a.s.l.)Depth (m)
Private well Hrašćica176.0015.0
PDS-5178.3631.0
PDS-6184.0725.0
PDS-7175.7142.5
P-1529187.328.0
P-1530183.727.5
P-1556193.0315.6
P-2500167.815.20
P-4039167.768.0
SPV-11177.6940.0
Table 6. Yearly averaged isotopic composition of precipitation at the rain gauge in Hraščica (177 m a.s.l.).
Table 6. Yearly averaged isotopic composition of precipitation at the rain gauge in Hraščica (177 m a.s.l.).
Yearδ18Ο (o/oo)δ2H (o/oo)d-excess (o/oo)
2017 (VI-XII)−7.78−52.98.7
2018 (I-VI)−9.96−70.49.3
2018 (VII-XII)−8.09−57.17.6
2019 (I-VI)−9.68−67.310.1
Table 7. Minimum, maximum, and averaged isotopic composition of the groundwater by sampled well.
Table 7. Minimum, maximum, and averaged isotopic composition of the groundwater by sampled well.
Observation Wellδ18Ο (o/oo)δ2H (o/oo)d-excess (o/oo)
minmaxaveragesdminmaxaveragesdaverage
Private well Hraščica−11.47−9.06−10.000.65−81.7−65.4−70.44.79.6
P-1529−11.20−8.87−9.760.64−77.6−61.2−68.24.49.9
P-1556−10.85−8.62−9.620.64−76.0−58.8−66.34.810.2
SPV-11−11.79−9.13−10.320.63−81.2−62.3−71.94.110.7
P-2500−10.51−8.42−9.310.56−71.3−58.6−64.13.510.4
P-4039−10.33−8.26−9.460.54−71.9−58.5−65.43.910.2
PDS-5−10.97−8.76−9.740.61−77.2−61.0−67.44.310.5
PDS-6−10.93−8.95−9.790.55−77.7−61.1−68.24.610.1
PDS-7−10.97−8.91−9.790.58−77.5−60.9−68.14.510.2
Table 8. Minimum, maximum, and averaged isotopic composition of surface water.
Table 8. Minimum, maximum, and averaged isotopic composition of surface water.
Observation Pointδ18Ο (o/oo)δ2H (o/oo)d-excess (o/oo)
minmaxaveragesdminmaxaveragesdaverage
Plitvica stream−10.15−7.86−9.110.52−69.2−56.9−63.32.89.6
Drava River−11.64−8.99−10.210.73−80.1−61.1−70.74.811.0
Varaždin Lake−12.12−8.10−10.330.87−81.7−59.1−72.45.210.2
Gravel pit in Šijanec−9.47−3.36−6.671.31−65.7−38.8−51.16.52.3
Table 9. Correlation matrix between ground and surface waters for δ18Ο (o/oo) values.
Table 9. Correlation matrix between ground and surface waters for δ18Ο (o/oo) values.
MatrixPrivate WellP-1529P-1556SPV-11P-2500P-4039PDS-5PDS-6PDS-7Plitvica StreamDrava RiverVaraždin LakeGarvel Pit Šijanec
Private well1.00
P-15290.931.00
P-15560.960.911.00
SPV-110.760.760.801.00
P-25000.910.940.890.731.00
P-40390.750.780.780.660.881.00
PDS-50.930.940.930.780.920.781.00
PDS-60.950.980.920.730.930.770.931.00
PDS-70.930.970.950.700.930.820.960.961.00
Plitvica Stream0.040.170.170.000.070.030.180.170.221.00
Drava River0.670.590.580.360.480.450.480.620.53−0.041.00
Varaždin Lake0.560.480.570.850.550.530.510.470.54−0.310.511.00
Gravel pit Šijanec-0.63−0.62−0.67−0.67−0.63−0.48−0.66−0.56−0.62−0.29−0.20-0.291.00

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/w12020379

AMA Style

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(2):379. https://doi.org/10.3390/w12020379

Chicago/Turabian Style

Marković, Tamara, Igor Karlović, Melita Perčec Tadić, and Ozren Larva. 2020. "Application of Stable Water Isotopes to Improve Conceptual Model of Alluvial Aquifer in the Varaždin Area" Water 12, no. 2: 379. https://doi.org/10.3390/w12020379

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