Tracer methods are useful for investigating groundwater travel times and recharge rates and analysing impacts on groundwater quality. The most frequently used tracers are stable isotopes and tritium. Stable isotopes of oxygen (δ¹⁸O) and hydrogen (δ²H) are mainly used as indicators of the recharge condition. Tritium (³H) is used to estimate an approximate mean groundwater age. This paper presents the results of an analysis of stable isotope data and tritium activity in Croatian groundwater samples that were collected between 1997 and 2014 at approximately 100 sites. The composition of the stable isotopes of groundwater in Croatia originates from recent precipitation and is described using two regional groundwater lines. One of them is applied to groundwater accumulated in the aquifers in the Pannonian part of Croatia and the other is for groundwater accumulated in the Dinaric karst of Croatia. The isotope content shows that the studied groundwater is mainly modern water. A mix of sub-modern and modern water is mostly accumulated in semi-confined porous aquifers in northern Croatia, deep carbonate aquifers, and (sub)thermal springs.

The Varaždin aquifer represents the main source of water for public supply, agricultural, and industrial purposes in the Varaždin County in NW Croatia. In the last decades, this area has experienced contamination of groundwater with nitrates. This study describes the conceptualization of the Varaždin aquifer for the purpose of developing numerical model of groundwater flow and nitrate transport. Within the study, three important elements are defined: aquifer geometry, recharge from precipitation, and other boundary conditions. 3D aquifer model revealed that Varaždin aquifer consist of three layers: upper aquifer, semipermeable interlayer, and lower aquifer. The Wetspass-M model was used for the assessment of spatial and temporal distribution of water balance components for the period 2008–2017. Results of the model indicate that the average annual precipitation is distributed as 34% groundwater recharge, 21% surface runoff, and 45% actual evapotranspiration. The maps of equipotential lines show the behavior of the aquifer system and define boundary conditions, i.e., recharge and discharge areas of the aquifer: an inflow boundary from Drava River and accumulation lake Varaždin on the northwest and north, no flow boundary on the west and south, and an outflow boundary on the east.

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 δ¹⁸Ο and δ²H in ground- and surface waters and precipitation samples were carried out. The relationship between ratios δ¹⁸Ο, δ²H, 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.

This paper presents the first study that assesses the vulnerability and risk of contamination of groundwater in the Varaždin aquifer system. The alluvial aquifer system is mostly unconfined with favorable hydrogeological features. Out of three wellfields, two still operate at full capacity, while the Varaždin wellfield, once the major source of drinking water, has been abandoned due to high concentrations of nitrates in the groundwater. Index-based methods are employed to assess groundwater vulnerability—two DRASTIC-based methods, standard and P-DRASTIC; two SINTACS-based methods, with normal and severe weighting strings; and the GOD method. Hazard is evaluated according to recommendations from the EU COST 620 action, while the risk intensity of the resource is calculated using the results of vulnerability and hazard assessments. The results reveal that for all vulnerability models, the resulting maps have a similar distribution pattern of vulnerability classes, with the high vulnerability class prevailing. However, notwithstanding the generally high groundwater vulnerability, a moderate resource contamination risk prevails as a consequence of a dominantly low hazard index. The validation of the groundwater vulnerability models demonstrates a weak correlation between the vulnerability indices and mean nitrate concentrations in groundwater. Conversely, a significantly higher correlation coefficient (0.58) is obtained when the groundwater vulnerability index is replaced by the resource risk intensity index, indicating that the results of resource risk intensity assessments are superior to groundwater vulnerability results in predicting the level of groundwater contamination.

This study investigates groundwater levels (GWLs) in the alluvial aquifer of the Sava River valley, located in the north-western part of Croatia. It provides the first quantitative assessment of groundwater levels using machine learning in this part of Europe. Groundwater levels from 1998 to 2017 were predicted using support vector regression (SVR). The input variables were initially monthly data on two basic elements that influence groundwater dynamics (precipitation and the Sava River levels). Later, GWLs from the previous month (GWL-1) were added as an additional predictor. Results demonstrated that the SVR model effectively predicts groundwater levels. Introducing GWL-1 reduced RMSE and MAE values by more than 47% and 46%, respectively, while increasing the R² value by over 36%. The improvement was more pronounced farther from the Sava River, since GWLs near the river are more directly influenced by river stage fluctuations, diminishing the impact of GWL-1. Compared to the existing regional numerical model (NM), the SVR model outperformed the NM with improvements of approximately 12% to 76% across performance indicators. Our findings suggest that the SVR model provides a reliable method for predicting groundwater levels at specific observation wells, making it a valuable tool for applications such as forecasting groundwater availability for farmers during dry periods and flood risk assessment during periods of heavy rainfall.