 Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system  Zoran Nakić1, Stanko Ružičić2, Kristijan Posavec1, Marta Mileusnić2, Jelena Parlov1, Andrea Bačani1 and Goran Durn2 1 Department of Geology and Geological Engineering, Faculty of Mining, Geology and Petroleum Engineering, University of Zagreb, Pierottijeva 6, 10000 Zagreb, Croatia; (zoran.nakic@rgn.hr) 2 Department of Mineralogy, Petrology and Mineral Resources, Faculty of Mining, Geology and Petroleum Engineering, University of Zagreb, Pierottijeva 6, 10000 Zagreb, Croatia doi: 10.4154/gc.2013.05 Geologia Croatica 66/1 55–76 16 Figs. 2 Tabs. Zagreb 2013 Geologia CroaticaGeologia Croatica and fast growth of the City of Zagreb have considerably af- fected the quality of groundwater in this area (NAKIĆ et al., 2001, BAČANI et al., 2002), with increasingly progressive groundwater pollution being observed over the last twenty- five years. High concentrations of pollutants including atra- zine, potentially toxic metals and chlorinated hydrocarbons confirm the impact of pollution sources on groundwater quality (NAKIĆ et al., 2010). In order to identify a coherent approach to mitigation of the impacts of different pressures, (point and non-point source 1. INTRODUCTION Groundwater is the only source of potable water for the City of Zagreb. The Zagreb aquifer system, where groundwater is naturally stored, is currently under threat due to decreas- ing water quantity and increasing concentrations of pollu- tants. Groundwater levels in the aquifer have particularly decreased over the last 40 years reaching such minimal levels as to cause problems to the City’s water supply during pe- riods of drought (POSAVEC, 2006). Industrial developments Ab sTRA CT This paper presents a conceptual model of the Zagreb aquifer system. The conceptual model can be applied to ground- water status risk assessment and pollution risk assessment at the local scale, thus satisfying both environmental and preventative/limitation objectives of the Water Framework Directive (WFD) and Groundwater Directive (GWD). Its main purpose is to apply a risk assessment procedure, according to the WFD requirements, and to serve as a founda- tion for setting up a numerical model of flow in both the saturated and unsaturated zones in order to identify pressure and impact effects on groundwater quality. The model is divided into two parts, taking into account the WFD require- ment to assess a risk for a wide range of source-pathway-receptor relationships. The Global Conceptual Model (GCM) provides insight into the processes and pressures at the level of the groundwater system. It contains the geological and hydrogeological characterization of the Zagreb aquifer system and the description of the most significant point and diffuse sources and pathways of pollution and processes influencing pollutant behaviour in saturated and unsat- urated zone of the groundwater system. The main pollutants of the Zagreb aquifer system are potentially toxic met- als, nitrates, pesticides, pharmaceuticals and chlorinated aliphatics. A Local Conceptual Model (LCM) supports pa- rameterization of the whole groundwater system through the description of heterogeneities and flow and solute parameters of the system components at two sites representing local conditions in the saturated (Stara Loza) and un- saturated (Kosnica) zones. This concept can be regarded as an effective tool for groundwater management of the groundwater system and its compartments and for communicating the conditions in complex groundwater systems with experts, policy makers and general public in an understandable way. Keywords: Global conceptual model (GCM), Local conceptual model (LCM), risk and status assessment, groundwater quality, Zagreb aquifer system Geologia Croatica 66/1Geologia Croatica 56 of pollution with respect to groundwater quality and over- exploitation of the groundwater), a risk assessment was ap- plied, as defined in CIS document No. 26 (COMMISSIOn of tHe europeAN CommuNItIeS, 2010). Risk assess- ment tries to draw a causal chain linking the origin of a haz- ard or pressure (e.g. an identified or estimated loading of polluting substance) along an environmental pathway to con- sequences for human health or the environment. This is gen- erally known as a source-pathway-receptor (SPR) paradigm (CommISSIoN of tHe europeAN CommuNItIeS, 2010). In the context of this paper, pollution risk is related to achieving the environmental objectives of the Water Frame work Directive (WFD, 2000/60/EC) and not in the per spective of risk to human health. In order to assess the risk of not reaching the goal of good groundwater quality in the Zagreb aquifer system, the Driver, Pressure, State, Im- pact, Response (DPSIR) analytical framework has to be fol- lowed, according to recommendations from the Analysis of Pressures and Impacts – Guidance Document No 3 (COM- mISSIoN of tHe europeAN CommuNItIeS, 2003). According to this approach, Driver is an anthropogenic ac- tivity that can affect the environment; Pressure is the direct effect of a driver, and Impact is an environmental effect of a Pressure. State is the environmental condition resulting from the Pressure and Responses are the measures taken for improving the State. To apply this analytical framework, a conceptual model is needed in the form of simplified representations, or work- ing descriptions, of the groundwater system under investiga- tion, in this case the Zagreb aquifer system. A new Groundwa- ter Directive (GWD, 2006/118/eC) recognized the signi ficance of the conceptual model in the groundwater risk and status assessment, which represents, by definition, the current un- derstanding of the groundwater system based on the know- ledge of its natural characteristics (GrAtH et al., 2007). the aim of building a conceptual model is to gain an adequate understanding of the groundwater system in order to identify pressures on the system and the relevant physical processes within it. According to the CIS guidance No. 17 (COMMIS- SIoN of tHe europeAN CommuNItIeS, 2007a), a conceptual model should describe the relevant geological characteristics, groundwater and surface water flow condi- Figure 1: Geographical location of the Zagreb aquifer system and locations of small scale experiment sites. Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 57 tions, hydrogeochemical (including terrestrial and aquatic ecosystems) processes, anthropogenic activities (relevant land uses) and their interactions. The accuracy and complex- ity of the models increases with the amount of, and confi- dence in, the available environmental information, so they become more effective and reliable descriptions of the sys- tem (BALDerACCHI et al., 2011). The main focus of this paper is to describe a conceptual model that can be applied according to the WFD require- ments in assessing pollution risk caused by point and non- point pollution sources with respect to groundwater quality in the Zagreb aquifer system at different scales, ranging from the scale of the groundwater system down to a single site, representing local conditions in the saturated (Stara Loza site) and unsaturated zones (Kosnica site). The philosophy behind this approach is to start develop- ing adequate groundwater management of the whole Zagreb aquifer system based on available knowledge of the system, and then to gain new information where needed to improve such understanding and the effectiveness of the measures applied for mitigating the risk of not achieving the environ- mental objectives of the WFD and GWD. In this respect, a number of iterations are expected, until the improved con- ceptual model can describe the measured data in a consistent way and with sufficient accuracy. the goal is to obtain a hy- drogeological conceptual model of the Zagreb aquifer sys- tem which will be used not only for describing, but also for quantifying the relevant geological and hydrogeological characteristics, flow and boundary conditions, hydrogeo- chemical processes, including natural groundwater quality and processes influencing contaminant behaviour in the un- saturated and saturated zones. 2. REsEARCH AREA 2.1. Geographical position and climate The study area (extent of the Zagreb aquifer system) is situ- ated in northwest Croatia (Fig. 1) and includes the City of Zagreb and its surroundings, covering an area of 350 km2. The broader region is characterized by highly variable lithology, pedological features and land use. The study area consists of a large alluvial plain bordered to the north by the Mt Medvednica mountain chain. The alluvial plain has two marked geomorphological features: the raised sealed terrace of the Sava River (varying in width down the rivers length), and a Holocene terrace. Numerous meanders of the Sava ri­ ver inundated fluvial cones and numerous bowl­shaped de- pressions abound in the alluvial plain. the climate of the City of Zagreb is classified as a mod- erately continental climate (Cfwbx in the Köppen climate classification system) with four separate seasons. Summers are warm, and winters are cold, without a discernible dry season. the average temperature in winter is −0.5 °C and the average temperature in summer is 22.0°C. the average rain- fall on the basis of ten years observations (1998–2008) is 928 mm year–1 (CROATIAn METEOROLOGICAL AnD HyDroLoGICAL ServICe). Figure 2: Simplified exploration core according to ŠIMUNIĆ & BASCH (1975). Geologia Croatica 66/1Geologia Croatica 58 of a wide variety of soil types. the floodplain soils, formed in materials eroded from sediments present in the catchment basin, are represented by specific mineralogical, chemical and textural characteristics. The investigated area of the Zagreb aquifer system en- compasses predominantly three pedologic units: Fluvisols, Stagnic Podzoluvisols (Pseudogley) and Eutric Cambisols, developed on Holocene deposits (BoGuNovIĆ et al., 1998). fluvisols (AC profile) are found on alluvial plains, river fans and valleys (FAO, 1998). This type of soil is character- ized by weak horizon differentiation, although a distinct Ah- horizon may be present. The texture of these soils is mainly loam, rarely clayey loam, while the structure is crumb to blocky (romIĆ et al., 2005)1. Stagnic podzoluvisols (Ap­BA­Bwr­Cr profile) is a hy- dromorphic soil which belongs to the Stagnosol class. Stag- nosols occur on flat to gently sloping land in cool temperate to subtropical regions with humid to perhumid climate con- ditions (fAo, 2007). In the study area, this type of soil de- veloped from Luvisols, i.e. it is secondary in origin, with a silty-clay-loam texture and a granular and very unstable structure (romIĆ et al., 2005)1. eutric Cambisols (A­Bv­r profile) belong to the cambic class of soils, characterized by a cambic horizon. Cambisols are characterized by slight or moderate weathering of the parent material and by the absence of appreciable quantities of illuviated clay, organic matter, Al and/or Fe compounds (fAo, 2007). In the study area, this type of soil has a humus­ accumulative horizon which is mainly mollic and noncalcar- eous, developed by pedogenc processes from rendzine on marls and Holocene deposits. the texture of this soil is silty­ clay to clay-loam (mainly loamy), developed on marls and Holocene deposits, respectively (romIĆ et al., 2005)1. The Fluvisol type of soil is developed at small scale Ko- snica and Stara Loza experiment sites, and is texturally mainly loamy, in some parts clayey with loam. Structure of this type of soils is mainly granular. In some parts, soil ho- rizons are red to black in colour, due to Fe and Mn oxides enrichments. 3. CONCEPTUAL MODEL DEVELOPMENT The primary purpose in developing the conceptual model of the Zagreb aquifer system was to serve as a foundation for setting up a numerical model of flow and transport in the sa­ turated and unsaturated zones, in order to identify pressure and impact effects on groundwater quantity and quality. How ever, over time, its purpose has been extended and it can now be regarded as a tool for communication with ex- perts, policy makers, or the general public, and can provide the basis for the derivation of measures for groundwater pro- tection. With relation to the WFD objectives, the conceptual A small scale experiment site in the unsaturated zone is located near the village of Velika Kosnica, about eight hun- dred metres from the right bank of the Sava river (45°46’ n; 16°05’ E) (Fig. 1). The location is close (straight line dis- tances) to Zagreb airport (3.5 km), the Zagreb bypass (250 m) and the largest Croatian landfill site Jakuševec (4.5 km). A small scale saturated zone experiment site is located in the area of the abandoned water abstraction site at Stara Loza on the left bank of the Sava river (45° 47’ 37’’ N, 15° 52’ 33’’ E) (Fig. 1). 2.2. Geological setting The geology of the studied area is described in the explana- tory notes for Basic Geological map of Sfry 1:100 000 (BASCH, 1983; ŠIKIĆ et al., 1979; pIKIJA, 1987). The Zagreb aquifer system is composed mainly of Qua- ternary sediments, which have been studied in detail by ŠI­ muNIĆ & BASCH (1975). four horizons of marshy­loess and three horizons of fluviatile deposits were recognised in a 120 m long exploration core (Fig. 2). ŠIKIĆ el al. (1979) and ŠIKIĆ & BASCH (1975) con- sider that detailed division of the Pleistocene deposits of the western part of the Zagreb environs is difficult, apart from division into genetic types, which are pondy-loess, continen- tal loess without carbonates, alluvial deposits of the third Sava river terrace and cave­fillings. Lower pleistocene de- posits are predominantly composed of yellowish-red, yel- lowish-orange and yellowish-brown, clayey silts/silty clays with sporadic lenses and interbeds of gravelly-sands, up to a few decimetres in thickness. Highly variable coefficients of sorting, compared to the average values, possibly indicate that loess, mixed with local materials (loessoid) has contrib- uted to the sediment composition (veLIĆ & DurN, 1993). The Middle Pleistocene unit is relatively uniform in composition. While the lower and middle part is predomi- nantly composed of grey coloured sands, the upper part com- prises grey coloured or red to yellowish-brown mottled silt and clay sized material. During warmer periods, the top of these sediments were exposed to pedogenic processes (veLIĆ & DurN, 1993; veLIĆ & SAftIĆ, 1996). The Upper Pleistocene unit is characterized by frequent lateral changes of gravels, sands, silts and clays. In the upper part, varves are recognized, as alternating, millimetre thick light coloured, fine­grained sands and dark grey to black silts. the Holocene is composed of pale, yellowish­ grey coloured gravels and sands in which limestone cobbles prevail. Both small scale experiment sites are situated in the Hol- ocene sediments. 2.3. Pedological settings The large heterogeneity of the parent material, combined with the climate and geomorphology, has caused the development 1 romIĆ, D., romIĆ, m., DoLANJSKI, D., StrIČevIĆ, I., oNDrAŠeK, G., mAurovIĆ, N., KoNDreS, N., HuSNJAK, S. & HeNGL, t. (2005): održivost agro­ekosustava na području Grada Zagreba s obzirom na onečišćenost teškim metalima. [the sustainability of agro-ecosystems of City of Zagreb with respect to potentially toxic metals contamination -in Croatian] Unpubl. report. Faculty of Agri- culture, University of Zagreb, Zagreb. Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 59 model of the Zagreb aquifer system is intended to assess the risk of failing to achieve good groundwater status. In setting up the conceptual model, separate components have been taken into account. Firstly, subsurface character- istics and dominant processes in the saturated and unsatu- rated zones of the Zagreb aquifer system have been address ed and incorporated in the model. Then, hydrogeological pa- rameters are quantified and variables are introduced that quan tify the exchange of fluxes between different compart- ments, e.g. exchange of water between the shallow aquifer and the Sava River. The current stage of construction focuses on qualifying the significance of point and non­point pollu- tion sources on groundwater deterioration. It is recognized that in evaluating the anthropogenic activities and their in- teractions, anthropogenic and non- anthropogenic effects need to be clearly distinguished. Hence, relevant natural occurring substances in the groundwater system are identified and the importance of the natural background load assessment is highlighted. It is also recognized that, in relation to the ob- jectives set, there is a knowledge gap in understanding the behaviour of the groundwater system at small and medium scales, with regard to its heterogeneity and the determination of groundwater pathways and mixing rates. According to the WFD requirements, risks need to be assessed for a wide range of source-pathway-receptor rela- tionships, at scales from the very local through the medium scale, (for example, impact on individual abstraction bore- holes or terrestrial ecosystems) to whole groundwater bodies (CommISSIoN of tHe europeAN CommuNItIeS, 2010). Hence, further improvement of the conceptual model of the Zagreb aquifer system is made according to the rec- ommendation of CIS guidance no. 15 (COMMISSIOn OF tHe europeAN CommuNItIeS, 2007b). the model is divided into two parts: a) Global Conceptual Model (GCM), which provides an insight into the processes at the level of groundwater system and b) Local Conceptual Model (LCM), which supports parameterization of the groundwater system through small scale experiments conducted at two sites: Ko- snica, which provides an insight into the behaviour of con- taminants in the unsaturated zone and Stara Loza, which pro vides an insight into small scale heterogeneity. The development of both types of conceptual models and the current understanding of the groundwater system and its compartments are presented here. This paper further in- vestigates the sources of contamination and contaminants in the urban water cycle that are directly linked to the deterio- ration of groundwater status in the area of the City of Zagreb. Due to multiple pollution sources and the wide range of po- tential hazards in this highly urbanized region, it is quite dif- ficult to estimate exactly the risks of not meeting the WfD objectives. However, risks can be managed to some extent by using information about the sources of contamination and contaminant behaviour in the subsurface. 3.1. Concept of groundwater status and risk assessment Conceptual models are considered as a vital tool to support the implementation of all aspects of groundwater require- ments of the WFD and GWD (CommISSIoN of tHe eu- ROPEAn COMMUnITIES, 2009). It is clear that both the groundwater status and risk assessment should be based on a well developed conceptual model of a groundwater system. In the status assessment, a clear understanding of the behaviour of pollutants which characterize bodies of ground- water as being at risk of not meeting the WFD objectives is needed, which is also linked to the pressure and impact anal- ysis (risk assessment), in relation to the validation of article 5 of the WFD. If the water body fails to meet its objectives, or is at risk of failing to meet its objectives, then the cause of the failure (i.e. the pressure or combination of pressures) must be investigated (CommISSIoN of tHe europeAN COM MUnITIES, 2003). the potential benefit of the conceptual model of the Za- greb aquifer system is not only in using it to assess the risk of failing to achieve good groundwater status of the system (groundwater bodies within the aquifer system), but also to estimate the risk of pollution at a local scale, thus satisfying the prevent/limit objectives in the WfD and GWD. However, in case of such localized pollution it is necessary to use the conceptual model to assess the overall risk for the ground- water system/body and to take measures to limit the pollu- tion. This approach is followed in estimating the vulnerabil- ity of the Zagreb aquifer system to potentially toxic metal pollution at the Kosnica site. The idea is to estimate the leak- age potential of the potentially toxic metal pollution load through the unsaturated zone at a local scale and to use risk assessment procedures to estimate the concentration range of selected potentially toxic metals in the soil and the un- saturated zone, which may be considered as an acceptable addition of human influence to the natural background con- centration load that doesn’t represent a threat to the relevant receptor, (in this case groundwater itself). By using numeri- cal modelling of transport through the unsaturated zone, the potential of this zone to retard potentially toxic metals may be as sess ed, as a basis for subsequent calculation of back- ground concentrations in the unsaturated zone. This approach is in line with the concept of REIMAnn & GARRETT (2005), who proposed the introduction of ambient back- ground concentration, which occurs under slightly altered conditions, when elevated levels of element concentrations in soil or water result from long-term human impact such as agriculture, industry, and urbanization and are no longer nat- ural. 3.2. Global conceptual model (GCM) A Global Conceptual Model (GCM) of the Zagreb aquifer system provides an overview of the processes which occur at the level of the groundwater system. It contains the geo- logical and hydrogeological characterization of the aquifer system, including the definition of boundary conditions and conclusions of the sensitivity analysis of the factors influ- encing groundwater recharge. To distinguish between natu- ral and man-made concentrations, the natural background loads of selected inorganic parameters in the groundwater are identified, reflecting natural processes unaffected by hu- man activities. The concept of the risk assessment implies Geologia Croatica 66/1Geologia Croatica 60 that the impact of human activity on the groundwater system should be defined, specifically those impacts that threaten the achievement of the objectives of the WFD. To that end, the most significant sources and pathways of pollution and processes influencing pollutant behaviour in saturated and unsaturated zone of the Zagreb aquifer system are identified and specifically outlined here. 3.2.1. Geological model of alluvial aquifer system The aquifer system comprises two Quaternary aquifers (figs. 3 and 4). Hydrogeologically, the Quaternary deposits are divided into three basic units: the overburden of clay and silt; a shallow Holocene aquifer of medium­grain gravel mixed with sands; and deeper aquifers from the Middle and Upper Pleistocene, with frequent lateral and vertical altera- tions of gravel, sand and clay. 3.2.2. Hydrogeological characterization The Zagreb aquifer system (Fig. 1) is an alluvial uncon- fined aquifer with a water table connected to the Sava river (Fig. 3). Sava River represents the main source of recharge. the general groundwater flow direction is from west to east, or south-east. Although recharge also occurs from precipita- Figure 3: Schematic hydrogeological cross-section of the Zagreb aquifer system. Figure 4: 3D geological model of the Zagreb aquifer system. Figure 5: Negative trends of groundwater levels. Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 61 tion, alteration in the Sava river water levels dominantly in- fluence the changes in groundwater levels across the whole aquifer. The quantity of groundwater is diminishing continu- ously (Fig. 5). The main reasons for lowering of groundwa- ter levels are (1) extensive riverbed erosion due to upstream Sava river regulation in the Republic of Slovenia and gravel exploitation from the river; (2) embankment of the Sava river which stopped occasional flooding of the area and potential infiltration to groundwater; (3) excessive abstraction for mu- nicipal and industrial needs; and (4) prolonged periods of drought. Groundwater levels have already reached minimum levels (defined as the upper elevations of the well screens) on some well fields, causing water scarcity during droughts (POSAVEC, 2006). Analysis of the water quantity data of the Zagreb aqui- fer system showed that groundwater levels are declining on average at 1–2 m every 10 years (BAČANI et al., 2010). the total annual abstraction from the Zagreb well fields exceeds the annual renewable groundwater reserves (which from 1997 – 2007 averaged about 107×106 m3 per year), while pumping in the same period averaged about 125×106 m3 per year. Groundwater abstraction rates which exceed renewable groundwater reserves are recovered from the permanent ground- water reserves. The volume of permanent groundwater reser- ves in 1977 was 1.81×109 m3 and by 2007 was 1.68×109 m3, which is a decrease in permanent groundwater reserves of 7% (BAČANI et al., 2010). The Sava River, which is the largest river in Croatia, di- vides the Zagreb aquifer system into two parts. The river, (which is the main source of groundwater recharge within the aquifer system), is in direct hydraulic connection with the shallow aquifer, which has extremely high values of hy- draulic conductivity (up to 3000 m/day). Head contour map analysis (POSAVEC, 2006) showed that during high Sava river water levels, infiltration to ground­ water occurs over all parts of the flow while during medium and low water levels the river drains groundwater on some parts of the flow. Groundwater levels are also strongly affected by the small river dam. An average drop of the Sava river water level amounts to 0,4 m/km, while the drop of the water level down- stream from the dam can be up to 6 m over a rather small distance of tens of metres. This strongly affects the aquifer flow directions and groundwater levels in near vicinity of the river dam (Fig. 6). The thickness of the unsaturated zone of the Zagreb aq- uifer system varies from 8 metres in the nW part to 2 metres Figure 6: River dam location and head contour maps showing the impact of the dam on groundwater flow direction in its near vicinity. Geologia Croatica 66/1Geologia Croatica 62 Figure 7: Isopach map of the unsaturated zone according to RUŽIČIĆ et al. (2011). Figure 8: Boundary conditions of the Zagreb aquifer system. Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 63 in the Se part (fig. 7). An isopach map of the unsaturated zone is derived from differences of groundwater level and terrain elevations. Groundwater level data were collected from 406 wells by the State meteorological and Hydrologi- cal Service and Water Supply and Drainage Company. The zone consists predominantly of gravels and sands, with spo- radic clays. The upper part is dominated by silty to sandy material, with the lower part being mostly gravel. In some parts this material is interbedded with clay layers. Hydraulic conductivity of the unsaturated zone of the Zagreb aquifer system varies from 1.26 m/day in the upper part to 1015 m/day in the lower part. This zone is predomi- nantly composed of quartz, calcite, dolomite and plagioclase in the upper part, while clay minerals represent the dominant mineral phases in the lower part. 3.2.3. Boundary conditions the Zagreb aquifer is an aquifer with no flow boundary to the north, an inflow boundary to the west and south and an outflow boundary to the east (fig. 8). the Dirichlet bound- ary condition specifies the value of the head at the boundary. the Neumann boundary condition specifies the flux normal to the boundary and the Cauchy boundary condition speci- fies the head and the flux normal to the boundary. Determi- nation of inflow/outflow/no flow boundaries was based on analysis of head contour maps for high, medium and low groundwater levels (POSAVEC, 2006). The results showed that there is no significant inflow or outflow through the nort­ hern aquifer boundary, possibly due to the mainly imperme- able deposits of Mt. Medvednica which results in predomi- nantly surface flow. Inflow quantification through the southern boundary is rather difficult due to lack of data and observation wells. Bas ed on head contour map analysis and water balance anal- ysis it can be assumed that some inflow through southern boundary does exist but with rather different fluxes across the region. Boundaries within the aquifer are represented by the Sava River which is the dominant boundary condition (Fig. 8). In order to identify and quantify the relevance of the im- pact of land­use and climate effects on groundwater flow in general, sensitivity analysis of the recharge boundary condi- tion, i.e. infiltration from precipitation, was carried out on existing numerical groundwater flow model of the Zagreb aquifer. for simulation of groundwater flow, a numerical method, i.e. moDfLoW code was identified and used as an appropriate simulation model according to defined goals. A numerical model for groundwater flow was established, ca­ librated and validated for the whole of the Zagreb aquifer. Sensitivity analysis of the recharge boundary condition, i.e. infiltration from precipitation, carried out on an existing Figure 9: Regression models showing zones of higher impact of the Sava River. Geologia Croatica 66/1Geologia Croatica 64 numerical groundwater flow model of the Zagreb aquifer, has shown that recharge from precipitation does not present the relevant boundary condition. the Sava river was defined as the dominant and the most relevant boundary condition. Spatial zonation of areas where the Sava River has a high er impact on groundwater levels was analyzed using re- cession curve models. The analysis of groundwater level time series using recession curve models was performed with the Master Recession Curve Tool (POSAVEC et al., 2006). In total, 278 master recession curves were obtained for the 278 observation wells analysed. Analysis of the spatial distribution of the selected regres- sion model showed that logarithmic regression prevails in parts of the aquifer near the river Sava, while polynomial regression prevails in other parts of the aquifer (Fig. 9). These results are logical and reasonable with respect to the changes in groundwater levels which occur faster in the vi- cinity of the Sava River. In other parts of the aquifer where such strong boundaries do not exist, groundwater level changes occur less rapidly. Principal component analysis (PCA) – following the pro cedure by WInTER et al. (2000) – can be used to define hydrodynamic zones within the regional groundwater sys- tem with identical groundwater table changes or fluctuation patterns in time. This procedure was implemented in the Za- greb case study, as part of the sensitivity analysis, pointing to aquifer recharging conditions and enabling conclusions on prevalent groundwater flow conditions in parts of the Za- greb aquifer system to be drawn. Results of the PCA, implemented in the Zagreb case stu- dy example (NAKIĆ, 2005), enabled characterization of the groundwater fluctuation pattern at the locations of observa- tion wells, revealing diversity of hydrodynamic characteris- tics in different parts of the groundwater system (Fig. 10). On the right bank of the Sava River the groundwater fluctuation pattern is similar for the shallow and deeper lay- ers, and it mainly reflects the impacts of the regional ground- water flow from the west towards the east. At the edge of the aquifer system, the impact of the Sava River is less and the water wave propagation is significantly lagging in time, while in the zone located immediately next to the Sava River, intensive fluctuations in the groundwater tables reflect the rapid and intensive impacts of changes in river water levels on the groundwater tables. Further evidence of the Sava River – groundwater inter- action has been provided by NAKIĆ et al., (2004), who de- lineated hydrogeochemically homogenous areas within the Zagreb aquifer system, using the pHreeQC groundwater modelling tool (pArKHurSt et al, 1980) and Cluster anal- ysis, more specifically the Ward method of hierarchical tree clustering. They revealed the very similar or identical macro chemical composition of groundwater in the area bordering Figure 10: Hydrodynamic zones within the Zagreb aquifer system according to NAKIĆ (2005). Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 65 the river, which supports the evidence of strong impacts of the Sava River water levels on the groundwater levels adja- cent to the river. these results were confirmed by hydrogeo- chemical investigations, conducted by vLAHovIĆ et al. (2009). It was proven that the direct exchange between the Sava River and the aquifer weakens with distance away from the river, while the difference in hydrogeochemical charac- teristics between the Sava River water and groundwater in- creases. To identify the pattern of groundwater – Sava River wa- ter interaction, comprehensive investigations of isotopic composition in precipitation, groundwater and river water were carried out. HorvAtINČIĆ et al., (2011) and PAR- LOV et al. (2012) used measurements of stable isotopes of oxygen and hydrogen in precipitation, surface waters and groundwater in order to determine the relative contributions of the Sava River and precipitation on groundwater recharge along the left (HorvAtINČIĆ et al., 2011) and right (PAR- LOV et al. 2012) banks of the Sava river respectively. The results of these investigations are combined in Fig- ure 11, which indicates the relationship between δ2H and δ18O in: the Sava River and groundwater at the left bank of the Sava River in 2010 (HorvAtINČIĆ et al., 2011); ground- water and Bundek lake at the right bank of the Sava River in 2011 (PARLOV et al., 2012). For illustration purposes, long term records of Zagreb precipitation (1980–2003) and the Local Meteoric Water Line (LMWL) (vreČA et al., 2006) are also included (Fig. 11). It can be seen that all groundwater values lay very close to the correlation line of the Sava River showing its direct influence on the groundwater. Differences between particu- lar piezometers indicate the temporal and spatial variability in groundwater-surface water (Sava River) interactions. Ac- cordingly, mixing of water from the Sava River and ground- water is occurring in the near vicinity of the Sava River and it weakens further away. Analyses of the recharge boundary condition presented here show that recharge from precipitation does not repre- sent an important boundary condition in the Zagreb aquifer system. Due to the identified high impact of the Sava river levels on groundwater levels, the Sava river is defined as the dominant and the most relevant boundary condition in the numerical model. Changes in the Sava River levels are naturally related to climate scenarios, but primarily in the upper parts of the flow where the Sava river is fed by sur- face runoff and by groundwater drainage. Since the course of the Sava River through the Zagreb area is characteristic of the valley river, it doesn’t drain the whole Zagreb aquifer, yet it controls the processes of recharge and drainage. There- fore in the case of the Zagreb aquifer, (i.e. an alluvial aquifer with high hydraulic conductivity and a water table connected to the stream with significant oscillation in water levels due to lack of regulation), the stream i.e. the Sava River repre- sents the dominant boundary condition. 3.2.4. Natural groundwater quality and ambient background values Development of industry and the fast growth of the City of Zagreb during the last three decades have considerably affected the quality of the groundwater in the Zagreb aquifer system. The Zagreb area accommodates an extremely high number of diverse diffuse and point pollution sources caus- ing degradation of groundwater quality and posing a partic- ular threat for its future use. Although human activity may lead to water composi- tions which exceed drinking water standards, water quality limits may be also breached for elements such as F, As and Fe by entirely natural processes, the result of geochemical conditions existing in aquifer or due to the specific geology in the area (EDMUnDS & SHAND, 2008). Groundwater pollution is most easily recognised by the presence of artifi- cial organic compounds, but it is difficult to recognise in- cipient pollution in most groundwaters. In order to be sure that pollution is occurring, natural spatial and temporal characteristics and trends of the chem- ical substances must be ascertained. Under such conditions, it is particularly important to make a clear distinction be- tween geogenic and anthropogenic influences in groundwa- ter quality. This is aided by determination of the geochemi- Figure 11: Relationships between δ2H and δ18O in the waters of the Zagreb area. Geologia Croatica 66/1Geologia Croatica 66 cal background of naturally occurring substances which re flect natural processes unaffected by human activities. the new GWD (2006/118/EC), which introduces the term back- ground level, as: „the concentration of a substance or the value of an indicator in a body of groundwater correspond- ing to no, or only very minor, anthropogenic alterations to undisturbed conditions” particularly recognises this. Although the term geochemical background was intro- duced much earlier in the scientific terminology by HAWKeS & WEBB (1962), who defined it as the normal abundance of an element in barren earth material, human influences on the planet as a whole have been so pervasive for such a long time that it may be pointless to attempt to determine pre- settlement background values, because concentrations from presumed pristine areas in these locations are likely to have also been affected (KeLLy & PAnnO, 2008). An option is to define ambient background values under slightly altered conditions, when elevated levels of element concentrations in soil or water result from long-term human impact, such as agriculture, industry and urbanization, and are no longer natural (REIMAnn & GARRETT, 2005). The geochemical background is a response to variable recharge rates, rainfall composition, water-gas interactions in soil and the unsaturated zone, as well as the processes (min- eral dissolution, redox reactions, ion-exchange reactions and mixing) that have taken place along flow paths. under natural conditions, groundwater quality is determined by the sum of soil­modified atmospheric inputs plus water­rock interactions taking place at the soil-bedrock interface and from longer-term reactions taking place in the saturated zone. In order to fully understand the evolution of groundwater chemistry along a flow path, in depth knowledge of the background and pollu- tion loads in soils and unsaturated zone is needed. Research into soil geochemistry in the area of the Za- greb aquifer system was carried out in different professional and scientific projects at different scales. for the „Geochem- ical Atlas of the Republic of Croatia”, a total of 640 topsoil (0–25 cm) samples were collected in the Central Croatian region using systematic sampling on a 5km grid (HALAmIĆ & MIKO, 2009). numerous studies of metal contamination of soils have been made in the area of Zagreb and Zagreb County (NAmJeSNIK, 1994; MIKO et al., 2001; romIĆ, 2002; romIĆ & romIĆ, 2003; romIĆ et al., 20045; HA­ LAmIĆ & MIKO, 2009; SOLLITTO et al., 2010). A total of 420 topsoil (0–15 cm) samples were collected, using system- atic sampling on a 2km grid, with sampling density increas- ing to 1 km near the urban part of Zagreb (NAmJeSNIK, 1994). SOLLITTO et al. (2010) collected 916 topsoil (0–20 cm) samples in the Zagreb County. Compendious results for potentially toxic metals in soils from these studies, together with maximum permissible concentrations in agricultural soils defined by Croatian government regulation (OG, 2010) are presented in Table 1. It is evident that in some areas of Zagreb concentrations of cadmium, chromium, cooper, mer- cury, nickel, lead and zinc in soil exceed maximum permis- sible concentrations. the vertical distribution of trace elements in soil profiles has been extensively studied, mainly in order to define the mobility of certain contaminants from soil to groundwater. About 70 vertical profiles of soil were studied in the area of Zagreb and Zagreb County. The Department of Field Crops, Forage and Grassland from the Faculty of Agriculture car- ried out agricultural soil trace element contamination studies in the framework of several projects (romIĆ et al., 20036; romIĆ et al., 20045; romIĆ et al., 20051). romIĆ (2002) investigated twelve profiles of the most abundant soil types. profile depths and thicknesses of layers from which the sam- ples were taken vary, depending on soil type and morpho- logical features. The results of the aforementioned investi- gations show that cadmium, iron, manganese and nickel have increased concentrations with depth, while concentrations of lead, zinc, chromium and copper vary. Table 1: Geochemical topsoil data for the City of Zagreb and surrounding area from different studies (data from geochemical atlases and maximal per- missible concentrations from Croatian legislation are for comparison). Element Official GAZETTE (2010) NAMJESNIK (1994) SOLLITO et al. (2010) Geochemical Atlas of the Republic of Croatia (HALAMIĆ & MIKO, 2009) Croatia Zagreb Zagreb county Central Croatia Croatia maximal permissible concentration Min Max Min Max Med Min Max Med Min Max Med * ** *** Cd (mg/kg) 0.5 1.0 2.0 0.01 22.5 0.02 4.94 0.31 0.2 9.4 0.2 0.2 15.5 0.4 Cr (mg/kg) 40 80 120 9 161 11.5 400 51.2 28 524 74 18 524 88.2 Cu (mg/kg) 60 90 120 – – 3.64 1335 23.5 3 248 19 3 429 25.4 Hg (μg/kg) 500 1000 1500 18 1938 – – – 5 4535 50 5 4535 60 Ni (mg/kg) 30 50 75 4 651 0.7 488 29.7 12 427 33 9.2 427 47.5 Pb (mg/kg) 50 100 150 5 546 1 216 19.6 14 217 27 10 699 33 Zn (mg/kg) 60 150 200 12 1250 27.1 479 70.7 28 477 73 23 1432 88 Legend: * – sandy soil; ** – silty-loamy soil; *** – clayey soil. Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 67 Since the WFD (2000/60/EC) came into force, several hydrogeochemical investigations of the Zagreb aquifer sys- tem have been undertaken with the aim of characterising the dominant geochemical processes and sources influencing na­ tural groundwater chemistry. NAKIĆ (2001) and BAČANI et al. (2002) used statistical methods, trend analysis and mul- tivariate statistical analyses, in order to discriminate between natural and man-made concentrations on a qualitative basis. BrKIĆ et al. (2003) studied the intrinsic vulnerability of the Zagreb aquifer system, while NAKIĆ et al. (2004) outlined hydrogeochemically homogenous areas within the Zagreb aquifer system by multivariate statistical analyses and geo- chemical modelling tools, showing that it is possible to diffe- rentiate a shallow Holocene alluvial aquifer from the deep er Middle/Upper Pleistocene lacustrine-marsh aquifers. In order to distinguish between geogenic and anthropo- genic influences in groundwater quality of the Zagreb aqui- fer system, background levels of inorganic substances were calculated on several occasions (Table 2), using model-based objective statistical methods (NAKIĆ et al., 2007; NAKIĆ et al., 2008). Pre-settlement background values of elements were not determined due to the extensive temporal impact of human activity in the Zagreb area. An option used was to define ambient background values under altered conditions, according to the concept set out by REIMAnn & GARRETT (2005). For this purpose, the Visual Basic macro BACK- GROUnD, created by NAKIĆ et al. (2007) was applied as it uses algorithms that incorporate model-based objective statistical methods (the iterative 2-σ technique and the cal- culated distribution function), as defined by mAtSCHuL- LAT et al. (2000). Overall, the ranges of ambient background values for the selected chemical substances reflects local geological and geochemical conditions in the aquifer system as well as the slight impact of anthropogenic sources, e.g. agricultural activities or atmospheric inputs due to emissions from the burning of fossil fuels. With respect to nO3, ambient back- ground includes sources such as soil organic matter from the residue of fertilized and unfertilized crops, products of com- bustion, and evaporation of ammonia (NH3) from synthetic and organic n fertilizer and livestock waste, in addition to naturally derived nO3. At first glance, it seems that ranges of ambient back- ground values for particular chemicals, e.g. for sulphate are wide, however comparison to the range of estimated back- ground values from European aquifers (see e.g. SHAND & EDMUnDS, 2008), suggests that values reflect the specific groundwater chemistry in the Zagreb aquifer system. Similar evidence can be provided for the other elements and compounds under consideration. Variable ambient back- ground ranges of iron and manganese are due to rapidly chang ing concentrations of these elements across redox boun daries in the Zagreb aquifer system. As a consequence, concentrations may vary over four to five orders of magni- tude across the shallow and deeper aquifers. 3.2.5. Sources and pathways of pollution, pollutants and their behaviour Groundwater of the Zagreb aquifer system is facing in- creasing pressure from human activities and there is a threat of contamination by diffuse and point sources. Pathways from sources to receptors are outlined in Figure 12. The main point pollution sources in the area of the Za- greb aquifer system include the Jakuševec city landfill site, marshalling yard, wastewater treatment plant, septic tanks, illegal waste depositories, gravel pits, industrial facilities, Zagreb airport, fuel stations and domestic households. The main diffuse pollution sources are traffic, leaky sewage pipes and agriculture (use of liquid manure, composted materials and agrochemicals such as fertilizers and pesticides). These sources can generally be divided into active (permanent and temporary) and potential sources of pollution. The perma- nent sources of pollution include the location of wastewater discharges, wastewater treatment plant, septic tanks, leaky sewage system, semi-legal disposal of municipal and indus- trial waste, abandoned gravel pits covered with municipal and construction waste. Temporary sources of pollution in- clude storm water discharge locations, tailings disposal in open pit mines, manure and other waste materials dumps from farms, storm drainage system wastewater, and other agricultural activity (BAČANI et al., 2010)2. Potential sources of pollution can occur due to failures, malfunctions or neg- ligence (e.g. active and abandoned gravel pits). Pollutants, causing regional or local pollution which, de- pending on soil properties and hydrogeological conditions, Table 2: Local ambient background values of selected chemical parameters in groundwater of the Zagreb aquifer system. Chemical Parameters Normal range of ambient background values Upper limits of ambient background values Method of calculation Source Nitrate (mg N/l) 0.0 – 4.80 4.8 iterative 2-s technique NAKIĆ et al., 2007 0.0 – 5.9 5.9 calculated distribution function Dissolved oxygen (mg/l O2) 1.4 – 8.4 1.4* iterative 2-s technique NAKIĆ et al., 2008 Sulphate (mg/l SO4) 14.2 – 47.6 47.6 iterative 2-s technique NAKIĆ et al., 2008 Chloride (mg/l Cl) 1.7 – 35.3 35.3 iterative 2-s technique NAKIĆ et al., 2008 Iron (mg/l Fe) 0.0 – 11.4 11.4 iterative 2-s technique NAKIĆ et al., 2008 Manganese (mg/l Mn) 0.0 – 0.5 0.5 iterative 2-s technique NAKIĆ et al., 2008 * For dissolved oxygen a lower limit of ambient background value is applicable to determine the outliers below background fluctuation (low values of dissolved oxygen as indicators of extreme oxygen consumption) Geologia Croatica 66/1Geologia Croatica 68 may reach groundwater, can be divided into natural substan­ ces: (1) potentially toxic metals; (2) radionuclides; (3) road salt; (4) nutrients (mainly nitrogen); and synthetic substances: (1) pesticides; (2) chlorinated aliphatics and petroleum hydro- carbon contaminants; (3) organic wastewater contaminants (BALDerACCHI et al., 2011). potentially toxic metals, ni- trates, pesticides, pharmaceuticals and chlorinated aliphatics are the main pollutants of the Zagreb aquifer system. Pathways of pollutants from the sources of pollution to the environmental receptors (atmosphere, soil, vegetation, Sava River, Quaternary alluvial sediments) ending in the groundwater of the Zagreb aquifer system are shown in Fig. 12. They include: fallout from the atmosphere onto vegeta- tion, soils, sediments and surface waters; runoff from the soil/sediment surfaces and other urban surfaces to the Sava river and other surface waters; flooding of the Sava river onto the Zagreb floodplain; weathering of rocks and sediments; leaching from point sources (e.g. Jakuševec landfill) and pol- luted soils and sediments. Potentially toxic metals Groundwater concentrations of lead exceed the maximum permitted concentrations in some locations though the over- all trend is one of decreasing levels as a consequence of a reduction in the use of leaded petrol. In some wells, the con- centration of manganese is extremely high (BAČANI et al., 2007)3. The occurrence and concentration of Mn in ground- water is controlled by many factors, the main ones being rock geochemistry, water chemistry and microbiological activity. Water chemistry, in particular pH, redox potential (eh), dis- solved oxygen (DO), and dissolved organic carbon (DOC) influence mobilisation of mn and control its speciation and concentration in the water environment. Micro-organisms can play an important role in Mn mobilisation in the envi- ronment and can both enhance or inhibit concentrations in groundwater. The impact of microbiological activity on Mn behaviour in water is frequently evident in the accumulation of Mn oxides on pipe surfaces. There is also the threat of potentially toxic metals leach- ing from point sources and polluted soils and sediments to the Figure 12: The most significant anthropogenic sources (grey boxes), pathways (arrows) and receptors (white and black boxes) of pollution in the area of the Zagreb aquifer system. Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 69 groundwater. NAKIĆ et al. (2007), in their study of ground­ water quality at the Jakuševec landfill indicated such leach- ing of potentially toxic metals into both soil and groundwater. Different studies of soils and sediments in Zagreb aqui fer area indicated elevated concentrations of following trace ele- ments: Cu, Hg, Zn, pb, Cd, Ni, As. HALAmIĆ et al. (2003) showed that the spatial distribution of Cu in topsoils of the Zagreb area was mostly related to the winegrowing areas, i.e. to the anthropogenic input of this metal into the topsoil. They also assumed that anomalous values of As and Cd are mainly of natural origin (mineralization). Greater amounts of Hg in topsoil of Zagreb area were probably caused by fos- sil fuel combustion and airborne pollution from industrial plants (NAmJeSNIK et al., 1992). ALJINovIĆ et al. (1995)4 highlighted higher concentrations of lead in the topsoil of Kosnica and assumed traffic to be the source. elevated con- centrations of Zn, Pb and Cd in topsoils of the Zagreb aqui- fer system are a consequence of the atmospheric deposition of particles from urban sources such are industrial emissions, traffic, waste disposal, heating plants, etc. (romIĆ et al., 2004)5. Higher concentrations of Hg, Cd and Ni as observed in the vertical profile of soils in water protected areas of the Zagreb aquifer system (romIĆ et al., 2005)1 are probably related to airborne deposition and traffic. Anomalous con- centrations of nickel in the vicinity of the highway and the airport result from anthropogenic influences due to fuel com- bustion (romIĆ & romIĆ, 2003). High nickel concentra- tions are probably partly connected to the morphogenetic characteristics of the wider region, primarily basic and ul- trabasic magmatic rocks of the surrounding mountain range. Nitrates Sporadically high concentrations of nitrates are registered in the groundwater beneath the City of Zagreb, especially at the mala mlaka well field. It was confirmed that periodically high nitrate concentrations are a consequence of agricultural production (NAKIĆ et al., 2001; petoŠIĆ et al., 2012). Ni- trogen is present in soils as the nitrate ion (nO3–), ammo- nium ion (NH4+) and as a component of soil organic matter. nitrate leaching from soil depends on the amount, frequency and intensity of precipitation, soil properties, crop type and crop development stage, evaporation, soil tillage practices, and nitrogen fertilization regime (vIDAČeK et al., 1996, 1999; NemetH, 2006; JoSIpovIC et al., 2006; NemČIĆ et al., 2007). romIĆ et al. (2003)6 found that high concen- trations of nitrates result from flushing of previously unsatu- rated layers during periods of high rainfall. Soil solution in arable soil contains 20–40 mg nO3-n/l, rising more than tenfold after fertilization (romIĆ et al, 2003)6. As nitrates are very soluble in water, leaching to groundwater is highest when precipitation exceeds evapotranspiration and for irri- gated land at the beginning of the vegetation growing period. Although intensive agriculture (nitrogen and organic fer- tilizer) is the main potential nitrate pollution source in the Zagreb area, other sources include leakage from septic tanks and the sewage system, atmospheric deposition and other urban waste, as well as accident situations. Pesticides Pesticides present in the groundwater of the Zagreb aquifer system originate mainly from agricultural use and can be dispersed in the environment by diffuse and point sources. Diffuse pesticide pathways into groundwater are mainly leach- ing through the soil and the unsaturated zone. Point sources mainly concern farmyard runoff, septic tank leakage, acci- dental spills, plumes from landfills and poor handling of pes- ticides. GoJmerAC et al. (2006) found high concentrations of atrazine in the groundwater of the Velika Gorica well field. Groundwater samples from observation wells during 2000–2010 show ed that atrazine occasionally exceeded the maximum permissible concentration beneath agricultural areas (BA ČA NI & poSAveC, 2011)7. Soils with high pesticide concentrations occurred at Kosnica (romIĆ et al., 1995)8. this contamination resulted from infiltration through the un- saturated zone which has a high filtration coefficient favour- ing the undisturbed transport of contaminants. The continued presence of these pollutants over time reflects the persistence of atrazine and other pesticides. Agricultural pesticide inputs into the soil and ground- water are dependent on various factors. At the field scale, inputs into soil depend on spraying, crop interception and wash-off, losses through spray drift or volatilization pro- cesses (BALDerACCHI et al., 2011). pesticide pathways to groundwater mostly depend on: (1) sorption processes, which control the amount of pesticide available for leaching into the soil and aquifer, and (2) on transport processes through the soil and the unsaturated zone (JArvIS, 2007; KoHNe et al., 2009b). Such processes are strongly depend- ent on the physical, physicochemical and biological proper- ties of soils. Atrazine for example, can be detected in ground- water several years after the last application (BARAn et al., 2008; MORVAn et al., 2006; TAPPE et al., 2002) because it has slow remobilization from the soil and the unsaturated zone. Chlorinated aliphatics Chlorinated aliphatics are some of the most widespread ground water contaminants in the industrialised world (KÄST- nER, 1991; WILSOn, 1996). Because of spills and acci- dents, these hydrocarbons have become important contami- nants in both soil and groundwater (SCHAerLAeKeNS et al., 1999). The contamination associated with chlorinated aliphatics and petroleum hydrocarbon compounds is typi- cally of a point source such as: landfill with municipal and industrial wastes; leaky sewage collecting leachate from small industrial activities (e.g. Dry-cleaning; metal degreas- ing); petrochemical activities; leakage of underground stor- age tanks. Some of them are primary compounds in indus- trial processes, some are by-products of degradation. Many are recalcitrant, but under favourable conditions (especially anaerobic) they can be transformed and degraded through microbially mediated processes which include reductive de- chlorination, dehydrochlorination and dichloroelimination (ferGuSoN & pIetArI, 2000). Tetrachlorethylene (PCE) concentration in all wells with in the Zagreb aquifer system is relatively high, and in some cases Geologia Croatica 66/1Geologia Croatica 70 exceeding the maximum permitted concentrations (BAČANI et al., 2007)3. Mobility of tetrachlorethylene (PCE), a dense non-aqueous phase liquid (DnAPL), is described as moderate (FETTER, 1988). In the presence of water, PCE will form a separate phase. The half-life degradation rate in groundwater is estimated to be between 1 – 2 years, based on aqueous aer- obic biodegradation (HoWArD et al., 1991), but may be con- siderably longer under certain conditions. Organic wastewater contaminants Recent research, carried out in the Zagreb area, within the framework of the implementation of EU directives, is partly related to the determination of the qualitative status of the groundwater (BAČANI & poSAveC, 2011)7, and it con- firmed that part of the aquifer system is endangered from or­ ganic pollutants. Organic wastewater contaminants include pharmaceuticals, surfactants, flame retardants, plasticizers and sterols. Pharmaceuticals represent the most studied class and are a diverse group of chemicals used as medicinal prod- ucts for both humans and animals (HeBerer, 2002), as well as bioactive food supplements, cosmetics, fragrances, sun- screen agents, insect repellent and antiseptics (ELLIS, 2006). After application at hospitals or households, pharmaceuti- cals are excreted either in their parent form or as metabolites, and enter the aquatic systems via various pathways (FEnT et al., 2006). The main sources of contamination are sewage- system leaching, sewage sludge, septic tanks, re-use of re- claimed domestic wastewater for irrigation and manure, con- taminated landfill and accidental release. most of them are not eliminated during wastewater treatment (JoNeS et al., 2005). Due to the polar structure of most pharmaceutical compounds they are not significantly adsorbed in the subsoil and may leach into the groundwater aquifers from contami- nated surface water (HeBerer et al., 1997). Sediment samples from Gorjak Creek (a small water- course), were analysed by terZIĆ et al. (2011) in order to identify organic contaminants present in the polar fraction of freshwater sediment. A number of different contaminants were successfully identified including various pharmaceuti- cals (chlorthalidone, warfarin, terbinafine, torsemide, zolpi- dem, azithromycin, desmethylazitromycin, erythromycin and dehydrated erythromycin). AHeL & JeLIČIĆ (2000) showed that solid waste in the Jakuševec landfill contains significant amounts of phen- azone and propylphenazone. The attenuation of phenazone concentrations in the unsaturated zone and in groundwater seems to be much faster than the attenuation of propylphen- azone. It can be assumed that the retention of phenazone com pounds by aquifer sediments is low due to their high solubility in water which indicates that abiotic and/or biotic transformations are involved. Radionuclides Although the manufacturing, application and disposal of man -made radionuclides are under strict regulatory control and under normal conditions only minute amounts of radio- activity are released to the environment, nuclear accidents can cause large scale radioactive pollution. HorvAtINČIĆ et al. (2011) reported significant increase of tritium activity in the Sava river in June 2010, as well as in groundwater (wells in the area of petruševec) with a delay for 3–5 months related to the Sava River. This increase was explained by the release of tritiated water from the Krško Nuclear power plant 30 km upstream of Zagreb. Low level radioactive waste in a slag dump near the TETO thermal plant, as well as low ra- dioactive waste from hospitals, also presents potential threats to groundwater. Salts The main anthropogenic source of chloride is naCl used on roads to remove ice at temperatures around 0°C. other an- thropogenic sources of chloride can be dust-binding on roads, storage of salt, snow dumps, leakage from sewage and waste deposits, fertilizers and chemicals used at water and sewage treatment plants (LUnDMARK, 2003). Since chlo- ride is a mobile compound, it will be transported to ground- water. Chloride stored in the unsaturated zone will, during snowmelt, be flushed into groundwater and so increase its chloride concentration (NyStÉN, 1998). 3.3. Local conceptual model (LCM) – supporting GCM parameterization Local Conceptual Model (LCM) supports development of the Global Conceptual Model (GCM) through characteriza- tion of two sites representing local conditions in the satu- rated (Stara Loza) and unsaturated (Kosnica) zones. As pre- viously stated, the main goal is to estimate the risk of pol lution on a local scale through the description of heterogeneities of Figure 13: Sketch of the unsaturated zone profile at Kosnica. Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 71 the groundwater system and through identification of the rel- evant flow and solute parameters in both the unsaturated and saturated zones. Only the design of site experiments is presented here together with an outline of the main results, since both ex- periments are still in progress. 3.3.1. Case studies 3.3.1.1. Experiments in the unsaturated zone The aim of the unsaturated zone experiment is the iden- tification of solute transport parameters, necessary for model- ling the behaviour of potentially toxic metals. The catchment area of the Kosnica well field was chosen for the ex periment because of its accessibility and results of previous investiga- tions (romIĆ et al, 19958, 20045; ALJINovIĆ et al. (1995)4. A solute transport model will be the basis for the risk assess- ment of pollution with potentially toxic metals at a local scale. A prediction transport model of potentially tox ic metals through the unsaturated zone is based on results of field and laboratory experiments. Detail characterization of the unsaturated zone profile (5 m) at Kosnica, together with sedimentological, chemical and mineralogical analyses of soil and sediment samples were the basis for the experimental setup. The lower part of the profile consists of gravels with a sand component, while the upper part is dominated by gravels with silty to sandy mate- rial (Fig. 13). Fluvisols are developed on the top of the un- saturated zone profile. the texture of this type of soil is main ly loamy, in some parts clayey with loam, while the struc ture is mainly granular. Soil pH increases while electrical conduc- tivity generally decreases with depth. The mineral composi- tion of the profile samples includes quartz, calcite, dolomite and plagioclase in the lower part and clay minerals mainly in the upper part. Both the longitudinal dispersion coefficient (DL) and the distribution coefficient (Kd) are the main para­ meters needed for modelling dispersion and sorption processes, respectively. These processes are obtained from labora tory re- sults and will be verified by field and laboratory ex periments. In order to obtain the sorption parameters (distribution coefficient, Kd) of Cd, pb and Zn, a laboratory batch exper- iment was performed on different soil samples taken from a pedological burrow. Multi element solutions of different con- centrations were prepared from Pb(nO3)2, ZnCl2 and CdCl2. Current results from the sorption experiment showed that Pb had a slightly higher affinity for sorption on these soil tex- tures than Cd and Zn (ruŽIČIĆ et al., 2012). to obtain transport parameters and confirm this experi- ment, a field dispersion coefficient determination will be set up. Calcium chloride solution will be used because of its re- la tively low sorption potential and its suitability (SeuNtJeNS et al., 2001; AKHtAr m.S. et al., 2003; DIŞLI, e., 2010). Monitoring of parameters such as water content, electri- cal conductivity, soil water tension and concentrations of chlorides in percolating water is necessary for the calibration of models. therefore, field measurement instruments are in- stalled in a pedological burrow located in the proximity of the profile. Water content is measured using time domain reflectrometry (tDr) equipment; electrical conductivity (EC) using EC-probes (Fig. 14a) and soil water tension us- ing tensiometers (Fig. 14b). The soil column laboratory experiment will also be used for determination of both the longitudinal dispersion coef- ficient (DL) and the distribution coefficient (Kd). the longi- tudinal dispersion coefficient will be determined using break- through curves of chloride concentration. After the dispersion ex periment, soil columns will be irrigated by metal solutions 2 BAČANI, A., poSAveC, K., pArLov, J. (2010.): prva faza izrade programa mjera za zaštitu i sanaciju u zonama zaštite izvorišta. [The first phase of the criterion program for the protection and rehabilitation of the water source protection zones – in Croatian].– Unpubl. report. Faculty of Mining, Geology and Petroleum Engineering, University of Zagreb, Zagreb. 3 BAČANI, A., poSAveC, K., NAKIĆ, Z., perKovIĆ D., mILetIĆ, p., HeINrICH­mILetIĆ, m., pArLov, J. & BAZIJANeC, m. (2007.): elaborat zaštitnih zona vodocrpilišta Grada Zagreba – I. faza [Study of water well protection zones for City of Zagreb – first phase].– Unpubl. report. Faculty of Mining, Geology and Petroleum Engineering, University of Zagreb, Zagreb. 4 ALJINovIĆ, D., mIKo, S., DurN, G., tADeJ, N., KApeLJ, S. (1995): Ispitivanje dosega prodiranja površinskog zagađenja kroz krovinu. [Investigation of surface contamination leaching through unsaturated zone – in Croatian].– Unpubl. report. Faculty of Mining, Geology and Petroleum Engineering, University of Zagreb, Zagreb. 5 romIĆ, D., romIĆ, m., DoLANJSKI, D., StrIČevIĆ, I., oNDrAŠeK, G., mAurovIĆ, N., KoNDreS, N., muStAĆ, I., HuS- NJAK, S., HeNGL, t. (2004): Stanje onečišćenja tala na prostoru Zagrebačke županije [The soil pollution state of the Zagreb county area – in Croatian] Unpubl. report. Faculty of Agriculture, University of Zagreb, Zagreb. 6 romIĆ, D., romIĆ, m., DoLANJSKI, D., StrIČevIĆ, I., oNDrAŠeK, N., HuSNJAK, S., petoŠIĆ, D., BAŠIĆ, f., KISIĆ, I., meSIĆ, m., JurAČAK, J., DurN, G., BrKIĆ, Ž. & BoroŠIĆ, J. (2003): poljoprivredna proizvodnja na vodozaštitnim područjima zagrebačke županije i Grada Zagreba [Agricultural production in water protection areas of Zagreb County and City of Zagreb – in Croa- tian].– Unpubl. report. Faculty of Agriculture, University of Zagreb, Zagreb. 7 BAČANI, A. & poSAveC, K. (2011): Implementacija okvirne direktive o vodama eu u dijelu koji se odnosi na kvalitativno stanje voda na zagrebačkom i samoborskom području. [the implementation of the Water framework Directive of the eu waters in the part relating to the status of water quality in Zagreb and the Samobor area].– Unpubl. report. Faculty of Mining, Geology and Petroleum Engineering, Uni- versity of Zagreb, Zagreb. 8 romIĆ, D., romIĆ, m., StrIČevIĆ, I., petoŠIĆ, D. & KLAČIĆ, Ž. (1995): Istraživanje prodora potencijalnih polutanata kroz solum tla do podzemne vode na području „Kosnice“ [Study of potentially pollutants breakthrough through soil to groundwater at location „Kos- nica“].– Unpubl. report. Faculty of Agriculture, University of Zagreb, Zagreb. Geologia Croatica 66/1Geologia Croatica 72 prepared from Pb(nO3)2, ZnCl2 and CdCl2 for determination of the distribution coefficients (Kd). Transport of potentially toxic metals is modelled using HyDruS 1D (ŠImŮNeK et al., 1998). figure 15 represents development of the solute transport prediction model using parameters obtained by laboratory and field experiments. Mo dels 2 and 3 will be compared to evaluate the validity of soil column experiments in modelling solute transport for the field scale. It should be stated that the experiment at Kosnica can support the parameterization of the soil and unsaturated zone of a large part of the groundwater system. About 88% of the area is covered by fluvisol which is also developed at Kos- nica. thus, field and laboratory experiments of fluvisol at Kosnica will provide the basis for a risk assessment of pol- lution with potentially toxic metals not only at the local scale, but also on a groundwater system scale. 3.3.1.2. Experiments in the saturated zone the Stara Loza site is a well field previously used for Croatia’s capital (Zagreb) public water supply (Fig. 16). It has not operated since 1997. the well field has 5 pumping wells distanced approximately 500 – 1000 m from the Sava river. 15 head observation wells (11 operating) and 7 con- centration observation wells (4 operating) are concentrated in the surrounding area and are still in use for monitoring of ground water levels and quality. Tracer experiments in the saturated zone at Stara Loza were conducted in order to determine longitudinal and trans- versal dispersion parameters as well as the effective (seep- age) velocity. The obtained parameters would then be used as the input for the prediction models of contaminant trans- port in the Zagreb aquifer. The experiment was designed as a natural gradient tracer test, with one injection and fourteen observation wells, setup in two rows 25 and 50 m from the injection well, perpendicular to the ground water flow direc- tion. Due to the expected transverse dispersivities of 0.25 – 1.25 m and associated transverse spread of 7.0 – 15.0 m for the first row and 10.0 – 22.0 m for the second row of obser- vation wells, the distance between wells in the first and second rows was set to 1.5 m and 3.0 m, respectively. the artificial tracer used in the experiment was uranine (Na­fluores cein) which was chosen because of its good sorption behaviour and small mass (LEIBUnDGUT et al., 2009). The mass of the tracer was calculated with respect to water volume and the detection limit of the tracer and the sampling frequency were calculated with respect to expected breakthrough times. for the first row of observation wells 25 m from the injec- tion well, the sampling started 6 hours after tracer injection, and 12 hrs for the second row. The sampling campaign lasted Figure 14: Pedological barrow: a) installed TDR and EC probes; b) installed tensiometers. Figure 15: A scheme representing develop- ment of the solute transport prediction model, using parameters obtained by laboratory and field experiments. Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 73 for 40 days, during which 1598 groundwater samples were taken. Analysis of the tracer experiment results have shown that a major part of the uranine tracer bypassed the first row of wells possibly through small palaeomeanders or palaeo- channels, but appeared in the second row of observation wells. This clearly suggests pronounced local scale heterogeneity of the aquifer and preferential flow paths, but there is a lack of scientific evidence to further support such an interpreta- tion. Therefore further research including geophysics will be conducted. 2D and 3D electrical tomography, which is gen- erally applied in geophysical research of aquifers (ŠumANo­ vAC, 2007), is planned in order to gain a better understanding of lateral changes in lithological composition and hope fully detect preferential flow paths or palaeomeanders i.e. palaeo- channels, which would give further insight into probable tracer migration paths. 4. CONCLUsIONs The main purpose of developing a conceptual model of the Zagreb aquifer system is to apply a risk assessment proce- dure through linking the origin of pollution sources along environmental pathways to consequences for the status of groundwater dependent ecosystems (GDE) and groundwater itself. This will provide a baseline for achieving the environ- mental objectives, according to the WFD requirements. To that end, the conceptual model is intended to assess the risk of failing to achieve good groundwater and GDE status through qualitative and quantitative description of real system be- haviour. At the current state of knowledge and understand- ing of the Zagreb aquifer system, it is evident that a number of iterations in the development of the conceptual model will be needed in order to accurately describe all the relevant pro- cesses within the groundwater system and pressures on it. This is particularly recognized in attempts to describe the heterogeneities of the groundwater system components at different scales and to quantify the significance and effects of point and diffuse pollution sources on the system. Taking into consideration knowledge gaps in the under- standing of groundwater systems, a coherent approach to as- sess a risk for a wide range of source-pathways-receptor re- lationships at different scales is applied by dividing the con ceptual model in two parts: a) Global Conceptual Model (GCM), which gives an insight into the processes and pres- sures at the level of the groundwater system and b) Local Conceptual Model (LCM), which supports parameterization of the whole groundwater system through the description of heterogeneities of the system components at a site scale and through identification of the relevant flow and solute param- eters in both the unsaturated and saturated zones. The rea- soning behind this approach is to use the conceptual model not only to assess the risk of failing to achieve good status of the groundwater and GDE bodies within the Zagreb aq- uifer system, but also to satisfy the preventative/limitation objectives in the WFD and GWD, by estimating the risk of pollution at a local scale. Furthermore, this concept helps to extend the primary purpose of the conceptual model by de- veloping an effective tool for groundwater management of the Zagreb aquifer system, particularly through creation of the basis for delineation of groundwater protection measures Figure 16: Stara Loza site. Geologia Croatica 66/1Geologia Croatica 74 at different scales and through communicating conditions in complex groundwater systems with experts, policy makers and the general public in an understandable way. ACKNOwLEDGMENT the study was supported by the european Community 7th framework Project GEnESIS (226536) on groundwater systems. REFERENCEs AHeL, m. & JeLIČIĆ, I. (2000): occurrence of phenazone analgesics in landfill­leachate polluted groundwater, Issues in the analysis of environmental endocrine disruptors. Extended abstracts of 219th national Meeting of the American Chemical Society, San Francis- co, USA, March 26–30, 2000, 109–111. AKHtAr, m.S., SteeNHuIS, t.S., rICHArDS, B.K. & mcBrIDe, M.B. (2003): Chloride and Lithium Transport in Large Arrays of Undisturbed Silt Loam and Sandy Loam Soil Columns.– Vadose Zone Journal, 2, 715–727. doi:10.2113/2.4.715 BAČANI, A., vLAHovIĆ, t. & NAKIĆ Z. (2002): endangered ground- water quality of the Strmec well field.– In: RIBEIRO, L. (ed.): Pro- ceedings, 3rd International Conference on future groundwater re- sources at risk, Lisbon, June 25–27, 2001, 329–336. BAČANI, A., poSAveC, K. & pArLov, J. (2010): Groundwater quan- tity in the Zagreb aquifer.– In: ZUBER, A., KAnIA, J. & KMIECIK, e. (eds.): XXXvIII IAH Congress Groundwater Quality Sustain- ability, Krakow, September 12–17, 2010, 87–92. BALDerACCHI, m., BeNoIt, p., CAmBIer, p., eKLo, o.m., GAr­ GINI, A., GemItZI, A., GureL, m., KLÖve, B., NAKIĆ, Z., preDA, e., ruŽIČIĆ, S., WACHNIeW, p. & trevISAN, m. (2011): Groundwater pollution and quality monitoring approaches at European-level.– Critical reviews in environmental science and technology. doi: 10.1080/10643389.2011.604259 BArAN, N., LepILLer, m., & mouvet, C. (2008): Agricultural dif- fuse pollution in a chalk aquifer (trois fontaines, france): Influ- ence of pesticide properties and hydrodynamic constraints.– Journal of Hydrology, 358, 56 – 69. doi:10.1016.2008.05.031 BASCH, o. (1983): osnovna geološka karta SfrJ 1:100 000. tumač za list Ivanić Grad L33–81 [Basic Geological Map of SFRY 1:100 000, Geology of the Ivanić Grad sheet – in Croatian]– Geološki zavod, Zagreb, Savezni geolološki zavod, Beograd, 66 p. BoGuNovIĆ, m., vIDAČeK, Ž., HuSNJAK, S., SrAKA, m. & pe­ toŠIĆ, D. (1998): Zemljišni informacijski sustav Zagrebačke žu­ panije [Land Information System of Zagreb Country – in Croatian].– Agriculturae conspectus scientificus, 63/4: 219–231. BrKIĆ, Ž., KApeLJ, S., LArvA, o., mArKovIĆ, t. & vLAHovIĆ, t. (2003): Hydrogeology and hydrogeochemistry in the alluvial aq- uifer of the Zagreb area.– rmZ­m&G vol. 50, No.1, 75–78. CommISSIoN of tHe europeAN CommuNItIeS (2003): Anal- ysis of Pressures and Impacts. Guidance document no 3. ISBn 92– 894–5123–8, European Communities, Luxembourg. CommISSIoN of tHe europeAN CommuNItIeS (2007a): Gui­ dance on preventing or limiting direct and indirect inputs in the con- text of the Groundwater Directive 2006/118/EC. Guidance document No 17. technical report – 2007 – 012. ISBN 978­92 ­ 79­06277­3, European Communities, Luxembourg. CommISSIoN of tHe europeAN CommuNItIeS (2007b): Gui­ dance groundwater monitoring. Guidance document no 15. Tech- nical report – 2007 – 002. ISBN 92­79­04558­X, european Com- munities, Luxembourg. CommISSIoN of tHe europeAN CommuNItIeS (2009): Guid- ance on groundwater status and trend assessment. Guidance docu- ment No 18. technical report – 2009 – 026. ISBN 978­92­79­ 11374­1, european Communities, Luxembourg. CommISSIoN of tHe europeAN CommuNItIeS (2010): Guid- ance on risk assessment and the use of conceptual models for Groundwater. Guidance document no 26. Technical Report – 2010 – 042. ISBN 13 978­92­79­16699­0, european Communities, Lux- embourg. DIŞLI, e. (2010): Batch and Column experiments to Support potential- ly toxic metals (Cu, Zn and Mn) Transport Modelling in Alluvial Sediments Between the mogan Lake and the eymir Lake Gölbaşi, Ankara.– Ground Water monitoring & remediation 30/3, 125–139. doi: 10.1111/j.1745­6592.2010.01302.x eDmuNDS, W.m. & SHAND, p. (2008): Natural Groundwater Quality – Summary and Significance for Water resources management.– In: EDMUnDS, W.m. & SHAND, P. (eds.): natural Groundwater Quality, Blackwell Publishing Ltd., 441–462. eLLIS, J. B. (2006): pharmaceutical and personal care products (ppCps) in urban receiving waters.– Environmental Pollution, 144/1, 184 – 189. doi:10.1016/j.envpol.2005.12.018 eLZAHABI, m & yoNG, r.N. (2001): pH influence on sorption char- acteristics of potentially toxic metal in vadose zone.– Engineering Geology, 60, 61–68. doi:10.1016/S0013­7952(00)00089­2 FAO (1998): FAO–UnESCO soil map of the world: revised legend.– World Soil Resources Rep 60. FAO/UnESCO/ISRIC, Rome. fAo (2007): fAo – uNeSCo World reference Base for Soil resourc- es, First update. FAO/UnESCO/ISRIC, Rome. feNt, K., WeStoN, A., & CAmINADA, D. (2006): ecotoxicology of human pharmaceuticals.– Aquatic toxicology, 76/2, 122–159. doi:10.1016/j.aquatox.2005.09.009 ferGuSoN, J.f. & pIetArI, J.m.H. (2000): Anaerobic transforma- tions and bioremediation of chlorinated solvents.– Environ. Pollut., 107, 2, 209–215. doi:10.1016/S0269­7491(99)00139­6 FETTER, C.W. (1988): Applied hydrogeology. 2nd edition.– Merrill Pub, Columbus, 592 p. GoJmerAC, t., oStoJIĆ, Z., pAuKovIĆ, D., pLeADIN, J. & Žu­ rIĆ, m. (2006): evaluation of Surface Water pollution with Atra- zine, an Endocrine Disruptor Chemical, in Agricultural Areas of turopolje, Croatia, Bull. environ.– Contam. toxicol. 76/3, 490–496. doi: 10.1007/s00128­006­0947­5 GrAtH, J., WArD, r., SCHeIDLeDer, A. & QuevAuvILLer, p. (2007): report on eu guidance on groundwater monitoring devel- oped under the common implementation strategy of the water frame work directive.– Journal of environmental monitoring, 9, 1162–1175. doi: 10.1039/B710665N HALAmIĆ, J., GALovIĆ, L. & ŠpArICA, m. (2003): potentially tox- ic metal (As, Cd, Cu, Hg, pb and Zn) Distribution in topsoil De- veloped on Alluvial Sediments of the Drava and Sava Rivers in nW Croatia.– Geol. Croat., 56/2, 215–232. HALAmIĆ, J. & mIKo, S. (eds.) (2009): Geochemical Atlas of the re- public of Croatia.– Croatian Geological Survey, Zagreb, 87 p. HAWKeS, H.e. & WeBB, J.S. (1962): Geochemistry in mineral explo- ration.– Harper and row, New york. HeBerer, th. & StAN, H.J. (1997): Determination of clofibric aci­ dand n-(phenylsulfonyl)-sarcosine in sewage, river and drinking water.– Int. J. environ. Anal. Chem., 67, 113–124. doi:10.1080/ 03067319708031398 HeBerer, t. (2002): occurrence, fate, and removal of pharmaceutical residues in the aquatic environment: a review of recent research data.– toxicology Letters, 131, 5–17. doi:10.1016/S0378­4274(02) 00041-3 HorvAtINČIĆ, N., BAreŠIĆ, J., KrAJCAr BroNIĆ, I., oBeLIĆ, B., KÁrmÁN, K. & fÓrIZS, I. (2011): Study of the bank filtered groundwater system of the Sava River at Zagreb (Croatia) using Nakić et al.: Conceptual model for groundwater status and risk assessment – case study of the Zagreb aquifer system Geologia Croatica 75 isotope analyses.– Central European Geology, Vol. 54/1–2, 121– 127. doi: 10.1556/CeuGeol.54.2011.1–2.12 HoWArD, p.H., BoetHLING, r.S., JArvIS, W.f., meyLAN, W.m. & mICHALeNKo, e.m. (1991): Handbook of environmental Degradation Rates. PRInTUP, H.t. (ed.).– Lewis publishers, Chel- sea, MI. JArvIS, N. J. (2007): A review of non­equilibrium water flow and sol- ute transport in soil macropores: principles, controlling factors and consequences for water quality.– european Journal of Soil Science, 58/3, 523–546. doi: 10.1111/j.1365–2389.2007.00915.x JoNeS, o.A.H., vouLvouLIS, N., & LeSter, J.N. (2005): Human pharmaceuticals in wastewater treatment processes.– Critical Re- views in environmental Science and technology, 35/4, 401–427. doi:10.1080/10643380590956966 JoSIpovIC, m., KovAČevIĆ, v., ŠoStArIĆ, J., pLAvŠIĆ, H. & LIovIC, J. (2006): Influences of irrigation and fertilization on soy- bean properties and nitrogen leaching.– Cereal Research Commu- nications, 34/1, 513–516. doi: 10.1556/CRC.34.2006.1.12 KÄSTnER, M. (1991): Reductive dechlorination of tri- and tetrachlo- roethylenes depends on transition from aerobic and anaerobic con- ditions.– Appl. environ. microbiol., 57, 2039–2045. KeLLy, W.r. & pANNo, S.v. (2008): Some Considerations in Apply- ing Background Concentrations to Ground Water Studies.– Ground Water, 47, no. 6: 790–792. doi:10.1111/j.1745­6584.2008.00467.x KoHNe, J. m., KoHNe, S. & ŠImŮNeK, J. (2009b): A review of mod- el applications for structured soils: b) pesticide transport.– Journal of Contaminant Hydrology, 104, 36–60. doi:10.1016/j.jcon- hyd.2008.10.003 LIN, y.p., teNG, t.p., & CHANG, t.K. (2002): multivariate analysis of soil potentially toxic metal pollution and landscape pattern in Changhua County in Taiwan.– Landsc. Urban Plan., 62/1, 19–35. doi:10.1016/S0169-2046(02)00094-4 LUnDMARK, A. (2003): Predicting environmental impact of de-icing salt – a modelling approach.– Report 2003–10–23, Dept. of Land and Water resources engineering, KtH. mAtSCHuLLAt, J., otteNSteIN, r. & reImANN, C. (2000): Geochemical background – can we calculate it?– Environmental Geology, 39/9, 990–1000. doi: 10.1007/s002549900084 mIKo, S., HALAmIĆ, J., peH, Z. & GALovIĆ, L. (2001): Geochem- ical Baseline Mapping of Soils Developed on Diverse Bedrock from Two Regions in Croatia.– Geol. Croat., 54/1, 53–118. morvAN, X., mouvet, C., BArAN, N., & GutIerreZ, A. (2006): Pesticides in the groundwater of a spring draining a sandy aquifer: temporal variability of concentrations and fluxes.– Journal of Con- taminant Hydrology, 87, 176–190. doi:10.1016/j.jconhyd.2006. 05.003 NAKIĆ, Z., BAČANI, A. & vLAHovIĆ, t. (2001): Antropogeni utje- caj na podzemne vode crpilišta Strmec (Zagreb, Hrvatska) [Anthro- pogenic influence on groundwater of the Strmec well field (Zagreb, Croatia) – in Croatian with an English Abstract].– Rudarsko-geo- loško­naftni zbornik, 13, 25–38. NAKIĆ, Z. & HorvAt, S. (2002): Analiza hidrokemijskih podataka primjenom multitivarijatnih statističkih metoda na priljevnom po­ dručju crpilišta petruševec [Petruševec well field catchment area hydrochemical data analysis by multivariate statistical methods – in Croatian].– Hrvatske vode, 40, 271–290. NAKIĆ, Z., poSAveC, K. & pArLov, J. (2004): Hydrogeochemical model of the heterogeneous aquifer system: the case of Zagreb, Croatia.– In: Groundwater flow understanding from local to region- al scales, Congress XXXIII – 70 ALHSuD, Zacatecas City. NAKIĆ, Z. (2005): određivanje hidrodinamičkih zona [Determination of hydrodynamic zones – in Croatian, with an English Abstract].– Građevinar, vol. 57, No. 6, 413–421. NAKIĆ, Z., prCe, m. & poSAveC, K. (2007): utjecaj odlagališta ot- pada Jakuševec­prudinec na kakvoću podzemne vode [Impact of the Jakuševec-Prudinec waste disposal site on groundwater quality – in Croatian, with an english Abstract].– rudarsko­geološko­ naftni zbornik,19, 35–45. NAKIĆ, Z., poSAveC, K. & BAČANI, A. (2007): A visual Basic Spread­ sheet Macro for Geochemical Background Analysis.– Ground Wa- ter 45, 642–647. doi: 10.1111/j.1745­6584.2007.00325.x NAKIĆ, Z., poSAveC, K., KoNtreC, p., BAČANI, A. & pArLov, J. (2008): Geochemical background values of selected chemical pa­ ra meters in groundwater of Zagreb aquifer system.– In: TAnIGU- CHI, m., yoSHIoKA, r., SKINNer, A., AureLI (eds.): A. XXXvI IAH Congress, Integrating Groundwater Science and Hu- man Well-being, Toyama. NAKIĆ, Z., poSAveC, K. & pArLov, J. (2010): model­based objec- tive methods for the estimation of groundwater geochemical back- ground.– AQuA mundi, 1, 65–72. doi: 10.4409/Am­010­10­0010 NAmJeSNIK, K., pALINKAŠ, L., mIKo, S., DurN, G., poLIĆ, D. & KvrGNJAŠ, L. (1992): Lead, zinc, nickel and mercury in soil along road sides and of some rural and urban parts of Zagreb, Cro- atia.– rudarsko­metalurški zbornik, 39, 93–112. NAmJeSNIK, K. (1994): Distribucija teških metala u tlima grada Za- greba i njegove okolice. [Potentially toxic metals distribution in the soils of Zagreb and surrounding – in Croatian].– Master thesis. Uni- versity of Zagreb Faculty of Science and Faculty of Mining, Geol- ogy and Petroleum Engineering, Zagreb, 63 p. NemČIĆ JureC J., meSIĆ m., BAŠIĆ f., KISIĆ I. & ZGoreLeC, Ž. (2007): Nitrate concentration in drinking water from wells at three different locations in northwest Croatia.– Cereal Research Com- munication, 35/2, 533–536. doi: 10.1556/CrC.35.2007.2.92 NemetH, t. (2006). Nitrogen in the soil­plant system nitrogen balanc- es.– Cereal Research Communications, 34/1, 61–64. NyStÉN, t. (1998): transport processes of road salt in quaternary for- mations.– In: NyStÉN, t., & SUOKKO, T. (eds.): De-icing and dust binding – risk to aquifers.– Proceedings of an International sym posium, NHp report No. 43, Nordic Hydrological programme, Helsinki, 31–41. oG (2010): pravilnik o zaštiti poljoprivrednog zemljišta od onečišćenja [Regulation on protection of agricultural land from pollution – in Croatian].– official Gazette nr. 32/10. pArKHurSt, D.K., tHoreNStoN, D.C. & pLummer, N.L. (1980): pHreeQe – a computer program for geochemical calculations. U.S. Geological Survey Water Resource Investigations, 80–96. pArLov, J., NAKIĆ, Z., poSAveC, K. & BAČANI, A. (2012): origin and dynamics of aquifer recharge in Zagreb area.– In: MORELL, M. (ed.): International Conference on water, climate and environ- ment, Balwois 2012, Ohrid, Republic of Macedonia. petoŠIĆ, D., HuSNJAK, S., petoŠIĆ, B. & muStAĆ, I. (2012): Im­ pact of Agriculture Upon Groundwater Pollution in the City of Za- greb Water-Conservation Area.– In: MORELL, M. (ed.): Interna- tional Conference on water, climate and environment Balwois 2012, Ohrid, Republic of Macedonia. pIKIJA, m. (1987): osnovna geološka karta SfrJ 1:100 000. tumač za list Sisak L33–93 [Basic Geological Map of SFRY 1:100 000, Geo- logy of the Sisak sheet – in Croatian].– Geološki zavod Zagreb (1986), Savezni geolološki zavod, Beograd, 64 p. poSAveC, K. (2006): Identification and prediction of minimum ground- water levels of Zagreb alluvial aquifer using recession curve mod- els.– Unpubl. Dissertation, University of Zagreb, Faculty of Mining, Geology and Petroleum Engineering, 89 p. reImANN, C. & GArret, r.G. (2005): Geochemical background – concept and reality. Science of the total environment, 350, 12–27. doi:10.1016/j.scitotenv.2005.01.047 Geologia Croatica 66/1Geologia Croatica 76 romIĆ, m. (2002): Sadržaj, oblici i preraspodjela imisije teških metala u poljoprivrednim tlima šireg područja Zagreba [Content, forms and redistribution of potentially toxic metals imision in the agricultural soils of broader area of Zagreb – in Croatian].–Doctoral thesis. Fac- ulty of Agriculture, University of Zagreb, Zagreb, 219 p. romIĆ, m. & romIĆ, D. (2003): potentially toxic metals distribution in agricultural topsoils in urban area.–Environmental Geology, 43, 7, 795–805. doi:10.1007/s00254­002­0694­9 ruŽIČIĆ, S., mILeuSNIĆ, m. & DurN, G. (2011): unsaturated zone of Zagreb aquifer.– In: BANJAC, n. (ed.): Geology in Digital Age, proceedings of the 17th meeting of the Association of european Geological Societies. Belgrade, 2011, 203–208. ruŽIČIĆ, S., mILeuSNIĆ, m. & poSAveC, K. (2012): Building Con- ceptual and Mathematical Model for Water Flow and Solute Trans- port in the Unsaturated zone at Kosnica Site, Mining-Geological- Petroleum Engineering Bulletin, Vol. 25, 21–31. SCHAerLAeKeNS, J., mALLANtS, D., ŠImŮNeK, J., vAN Ge- NuCHteN m. th. & feyeN, J. (1999): Numerical simulation of transport and sequential biodegradation of chlorinated aliphatic hy- drocarbons using CHAIN_2D.­Hydrological processes, 13, 2847– 2859. doi: 10.1002/(SICI)1099–1085(19991215) SeuNtJeNS, p., mALLANtS, D., torIDe, N., CorNeLIS, C. & GEU ZEnS, P. (2001): Grid lysimeter study of steady state chloride transport in two Spodosol types using TDR and wick samplers.– Journal of Contaminant Hydrology, 51, 13–39. doi:10.1016/S0169­ 7722(01)00120­6 SHAND, p. & eDmuNDS, W.m. (2008): the Baseline Inorganic Chem- istry of European Groundwaters.– In: EDMUnDS, W.m. & SHAND, P. (eds.): natural Groundwater Quality, Blackwell Publishing Ltd., 22–59. SoLLItto, D., romIĆ, m., CAStrIGNANo, A., romIĆ, D. & BA­ KIĆ, H. (2010): Assessing potentially toxic metal contamination in soils of the Zagreb region (northwest Croatia) using multivariate geostatistics.– Catena 80, 182–194. doi:10.1016/j.catena.2009.11.005 ŠIKIĆ, K. & BASCH, o. (1975): Geološka zbivanja od paleozoika do kvartara u zapadnom djelu Zagrebačke regije. [Geological events from Paleozoic to Quaternary Period in west part of Zagreb region – in Croatian].-2nd scientific congress section for application of ge- ology, geophysics and geochemistry. Scientific Advisory Board for petroleum JAZu, Zagreb, A/5, 69–86. ŠIKIĆ, K., BASCH, D. & ŠImuNIĆ, A. (1979): osnovna geološka kar- ta SfrJ 1: 100.000. tumač za list Zagreb, l 38–80 [Basic Geologi- cal Map of SFRY 1:100 000, Geology of the Zagreb sheet – in Cro- atian].– Institut za geološka istraživanja, Zagreb (1972), Savezni geo lo ški zavod, Beograd, 81 p. ŠImŮNeK, J., ŠeJNA, m. & vAN GeNuCHteN, m. tH. (1998): the HyDruS­1D software package for simulating the one­dimension- al movement of water, heat, and multiple solutes in variably-satu- rated media, version 2.0. IGWmC­tpS­70, International Ground Water Modeling Center, Colorado School of Mines, Golden, Colo- rado, 162 p. ŠImuNIĆ, A. & BASCH, o. (1975): Stratigrafija kvartarnih sedimena- ta Zagrebačkog posavlja. [Quartary sediments stratigraphy of the Zagreb Posavlje region – in Croatian].– Geol. vjesnik, 28, 153–164. tAppe, W., GroeNeWeG, J., & JANtSCH, B. (2002): Diffuse atra- zine pollution in German aquifers.– Biodegradation, 13/1, 3–10. doi: 10.1023/A:1016325527709 terZIĆ, S. & AHeL, m. (2011): Nontarget analysis of polar contami- nants in freshwater sediments influenced by pharmaceutical indus- try using ultra-high-pressure liquid chromatography-quadrupole time­of­flight mass spectrometry.– env. pollution, 159, 2, 557–566. doi: 10.1016/j.envpol.2010.10.009 veLIĆ J. & DurN, G. (1993): Alternating Lacustrine­marsh Sedimen- tation and Subaerial exposure phases During Quaternary: prečko, Zagreb, Croatia.– Geol. Croat., 46/1, 71–90. veLIĆ J. & SAftIĆ, B. (1996): Dubinskogeološki odnosi područja smet lišta „Jakuševec“ – čimbenik sanacije [Deep geological rela- tionship of city landfill „Jakuševac” - recovery factor – in Croatian] Waste managment, Iv. International symposium, Zagreb, 197–205. vLAHovIĆ, t., BAČANI, A. & poSAveC, K. (2009): Hydrogeoche­ mical stratification of the unconfined Samobor aquifer (Zagreb, Croatia).– environ. Geol., 57/8, 1707–1722. doi: 10.1007/s00254­ 008-1452-4 vIDAČeK, Ž., BoGuNovIĆ, m., SrAKA, m. & HuSNJAK, S. (1996): Water discharges and nitrates from some soils of the Sava river valley. 6. Gumpensteiner Lysimetertagung „Lysimeter in Di- enste des Grundwassershutzes”.– BAL, Gumpenstein, 111–115. vIDAČeK, Ž., SrAKA, m., ČoGA, L. & mIHALIĆ, A. (1999): Ni- trates, potentially toxic metals and herbicides in soil and water of Karasica­vučica catchment area.– Agriculturae Conspectus Scien- tificus, 64/2, 143–150. vreČA, p., KrAJCAr BroNIĆ, I., HorvAtINČIĆ, N. & BAreŠIĆ J. (2006): Isotopic characteristics of precipitation in Slovenia and Croatia: Comparison of continental and maritime stations.– J. Hy- drology, 330, 457–469. doi: 10.1016/j.jhydrol.2006.04.005 WeDepoHL, K. H. (1995): the composition of the continental crust.– Ge- ochimica et Cosmochmica Acta, 59, 7, 1217–1232. doi:10.1016/0016­ 7037(95)00038­2 WILSoN, J. (1996): Chlorinated solvents biodegradation, epA, Ground- water Water Curr., 16, 1–4. Manuscript received September 18, 2012 Revised manuscript accepted November 19, 2012 Available online February 28, 2013