03 Kuhta.indd 41 �AB STRA CT Spring hydrographs and thermographs are a direct refl ection of all processes that occur within aquifer systems. There- fore, they contain signifi cant information about the hydrogeological characteristics of such systems. This article anal- yses hydrographs and thermographs from four springs located in the foothills of the Mt. Biokovo massif in southern Croatia. These springs are recharged by carbonate aquifers. Monitoring of yields and groundwater temperatures, as well as analysing hydrograph recessions, daily discharge and rainfall time series and water temperature dynamics, facilitated the identifi cation of the main properties of the aquifer system located in the hinterland of the individual springs. Signifi cant differences in the recharge mechanisms of individual springs were determined to be a conse- quence of varying geological conditions, degree of karstifi cation and conduit network characteristics. The results suggest that the Cretaceous and Palaeogene basinal carbonate deposits (Kotišina limestones and breccias), the hy- drogeology of which has not yet been studied, have the characteristics of permeable karstic rocks. Keywords: karst hydrogeology, spring hydrograph, time series analysis, Biokovo Hydrodynamic characteristics of Mt. Biokovo foothill springs in Croatia � Mladen Kuhta, Andrej Stroj and Željka Brkić Department of Hydrogeology and Engineering Geology, Croatian Geological Survey, Sachsova 2, 10 000 Zagreb, Croatia; (mkuhta@hgi-cgs.hr, astroj@hgi-cgs.hr, zbrkic@hgi-cgs.hr) doi: 104154/gc.2012.03 Geologia Croatica 65/1 41–51 6 Figs. 3 Tabs. Zagreb 2012 Geologia CroaticaGeologia Croatica 1. INTRODUCTION A spring hydrograph is a direct refl ection of all the physical processes that infl uence groundwater fl ow within an aquifer. Hydrographs of karst springs are characterised by pro- nounced, sudden changes of discharge, and a briefl y delayed response to rainfall (BONACCI 1987, 1993; GOLDS- CHEIDER et al., 2007), which differentiates them from hydro graphs of springs that drain other aquifers, i.e. fi ssured and intergranular aquifers (Fig. 1). Rapid spring responses to rainfall occur mainly as a result of sudden infi ltration to conduit systems related to the function of the subcutaneous (epikarst) zone (MANGIN, 1975; WILLIAMS, 1983, 2008; PRONK et al., 2009). When rainfall ceases, a sudden incre- ase in fl ow is initially followed by a nearly equal sudden de- crease, after which the decrease gradually slows, refl ecting an alteration in the system dynamics. Hydrographs, there- fore, refl ect the characteristics of two main types of fl ow in an aquifer: fast-fl ow, which is determined by conduit sys- tems and rainfall distribution properties, and base-fl ow, which is predominantly controlled by slow drainage of low hydraulic conductivity volumes in the aquifer (ATKINSON, 1977; PADILLA et al., 1994). Thus, the shape of a hydro- graph is a signifi cant indicator of an aquifer’s characteristics and serves as the most important time series in analyses of the hydrogeological characteristics of karstic aquifers. The Makarska Coast area extends along approximately 53 km of the Central Dalmatian coastline along the foothills and slopes of the Mt. Biokovo massif that stretches from Brela in the north to Gradac in the south (Fig. 2). This area is most ly supplied with water from the Regional Water Sup- ply System, which intakes water from the Cetina River near Zadvarje. This system is important during dry periods, when tourism raises the demand for drinking water. In addition, the water supply system of the Makarska Coast includes six smaller water intakes built on the springs along the base of Mt. Biokovo. These intakes include the springs of Baška Voda, Makarski Vrutak, Orašje, Grebice, Podgorski Vrutak and Izbitac in Drašnice. Although their total yield during the summer is far below the actual consumption rate, the cap- tured springs signifi cantly contribute to the stability of the Geologia Croatica 65/1Geologia Croatica 42 public water supply system, especially because the amount of water in the Regional System from the Cetina River is technically limited. During the tourist off-season, and par- ticularly during favourable hydrological periods, the local water intakes alone satisfy the majority of water demands. These springs are especially signifi cant due to their location high above the urban areas, which facilitates gravitational fl ow. The majority of the springs on the southern slopes of Mt. Biokovo have not been investigated, and data on their yield and hydrogeological characteristics were based on the empirical insight of the Makarska Water Utility Company. Thus, during investigations to identify sanitary protection zones for springs used in the public water supply system, special attention was paid to the establishment of a system to continuously monitor yields at the technically adequate water intakes, including the springs of Makarski Vrutak, Orašje, Podgorski Vrutak and Izbitac in Drašnice. These data enabled hydrograph analyses to be made. 2. GEOLOGICAL AND HYDROGEOLOGICAL SETTINGS The Dinaric karst is characterized by a succession of more than 8000 m of predominantly carbonate sediments (TIŠ- LJAR et al., 2002; VELIĆ et al., 2002), that was deposited and exposed to intense tectonic disturbances in several phas es from the Triassic period until the present. The main defor- mation episode began in the Late Cretaceous, when synsedi- mentary tectonics became stronger, reaching a maximum in the Oligocene/Miocene. This led to tangential movements and uplift of the Dinarides, including the Mt. Biokovo area. Due to the strong NE–SW-oriented regional tectonic stress, the main resulting structures, including folds, faults, thrusts and imbricate structures, are of the NW–SE strike (i.e., the Dinaric strike). Later orogenic movements shifted the re- gion al tectonic stress to the N–S and caused wrench tectonic deformations. The main part of the Mt. Biokovo massif consists of shal- low-water carbonates, mostly limestones with subordinate dolomites, deposited in the period from the Lower Juras sic to the Upper Cretaceous (Fig. 3). The basic structural form defi nes a large and prominent Biokovo anticline (MARINČIĆ et al., 1977) in the southwest. This anticline is in tectonic contact with highly deformed limestones that contain pelagic elements and carbonate breccias of the Upper Cretaceous- Palaeogene, (i.e., the Kotišina deposits) and/or a narrow belt of tectonically reduced fl ysch deposits. JELASKA et al. (2003) proposed a name for this basinal domain: the Kotišina- Tilovica trough. Palaeogeographically, this trough represents a NW embayment of the Budva-Cukali basin (CHOROW- ICZ, 1975). Basinal carbonate deposits, were tectonically emplaced into the Dinaridic segment, which had been previ- ously over thrust on them (KORBAR, 2009; KORBAR et al., 2010). In the SW, these carbonate units are in reverse contact or are overthrust over intensely deformed, predominantly clastic fl ysch deposits of the Middle and Upper Eocene, which comprise the foothills of Mt. Biokovo. In several places, gently inclined fl ysch slopes are covered by tens of metres of deluvial and proluvial Quaternary deposits. Recent debris is frequently found below the steep carbonate slop es. The carbonate deposits of the Mt. Biokovo massif are intensely tectonically fragmented and karstifi ed, and in hy- drogeological terms they form a permeable environment in which a deep karst aquifer has developed (Fig. 3). The high karst plateau is dominated by numerous deep dolines, which often contain karst pits. Several hundred pits have been in- Figure 1: Typical hydrograph charac- teristics of a spring draining a) a fi s- sured aquifer and b) a karstic (con- duit) aquifer (Q – discharge, P – rainfall, t – time, Dt – time diff erence between rainfall and hydrograph peak; modi- fi ed from BONACCI, 1993). Kuhta et al.: Hydrodynamic characteristics of Mt. Biokovo foothill springs in Croatia Geologia Croatica 43 Figure 2: Regional catch ments and locations of the investi gat- ed springs. Figure 3: Hydrogeological map of the study area (geological background ac- cording to MARINČIĆ et al., 1977, partly modifi ed). Geologia Croatica 65/1Geologia Croatica 44 vestigated through speleological activities on Mt. Biokovo. Ten of these pits were deeper than 250 m, and four exceeded a depth of 500 m. The deepest pit of Mt. Biokovo is the Am- fora pit, which has been investigated to a depth of 788 m (BAKŠIĆ & JALŽIĆ, 2001; BAKŠIĆ & LACKOVIĆ, 2002). In hydrogeological terms, deep pits are a signifi cant indica- tor of karstifi cation of the vadose zone and a deep ground- water table. The lowest absolute elevation of approximately 790 m a.s.l. is reached in the Vilimova Jama pit. The fl ysch deposits that comprise the SW foothills of the carbonate massif of Mt. Biokovo are lithologically heteroge- neous. However, due to the dominance of marls, this sediment series, when considered in its entirety, is impermeable to wa- ter. In the study area, the fl ysch outcrops appear to extend to a distance of 3 km from the coastline (Fig. 3), and in places reach an altitude of 500 m. Considering their structural posi- tion and their nearly uninterrupted stretch along the coastline and far beyond the study area, the fl ysch deposits form a com- plete hydrogeological barrier to groundwater movement from the carbonate hinterland toward the sea. Thus, they have a sig- nifi cant impact on the regional hydrogeology. According to previous regional hydrogeological studies (BOJANIĆ, 1980)1, groundwater from the Mt. Biokovo mas- sif and a portion of its spacious karstic hinterland drains pri- marily through strong submarine springs: Dupci NW of Brela, and Klokun near Drašnice (Fig. 2). They develop in locations where fl ysch and carbonate rocks come into con- tact below sea level. The springs function as the primary groundwater drains for the greater regional area, which ex- plains the discharge values during rainfall periods of up to several dozen m3/s. The Dupci is considered to be one of the strongest submarine springs of the Adriatic. In addition to these submarine springs, several perma- nent springs with lower discharge values, some of which are captured for use in the public water supply system, are lo- cated along the area of tectonic contact between permeable carbonate rocks and impermeable fl ysch. This research fo- cused on the hydrogeological characteristics of four captured springs located at the foot of steep slopes in the central and SE parts of the Mt. Biokovo massif between Makarska and Drašnice. The recharge of these springs is unrelated to in- fl ows from the regional hinterland. They are recharged from the carbonate deposits of the Mt. Biokovo massif, and they have their own recharge areas in the SW, marginal parts of the regional catchments of the Dupci (Makarski Vrutak and Orašje) and Klokun (Podgorski Vrutak and Izbitac) subma- rine springs. This fact was confi rmed by tracing from the Jama Male Marije pit in the Lađene area of Mt. Biokovo, which proved the groundwater connection with the Orašje and Grebice springs (Fig. 3). The Makarski Vrutak spring is captured by a 10-m long water intake gallery located at 160 m a.s.l. The gallery was built in partly consolidated breccia debris, although it can be assumed that there are fl ysch deposits (E2,3) in the bedrock. The nearest outcrops of carbonate rocks from the Upper Cretaceous (K2 3) and Palaeogene (Pg) ages form the structural unit of the Kotišina-Tilovica trough, located ap- proximately 100 m north of the water intake. The primary source is probably located in the tectonic contact zone be- tween the carbonate aquifer of Mt. Biokovo and the fl ysch, which is covered by recent debris, as indicated by a tempo- rary spring in the gully above the water intake. The Orašje spring is located at the foot of Mt. Bioko- vo’s steep slopes, at an altitude of approximately 465 m a.s.l. Groundwater is extracted by two small intake galleries 35 m apart and then transferred by pipelines to a nearby collection chamber. Although the intake is in an area composed of fl y- sch deposits located approximately 160 m away from its fault contact with the carbonate rocks of Mt. Biokovo, and the groundwater infl ow route is unknown, (draining through the debris or the more permeable sequences within the fl ysch), it is indisputable that the springs derive recharge from a karstic carbonate aquifer in the Mt. Biokovo massif, as con- fi rmed by the aforementioned tracing study. The Podgorski Vrutak spring is captured by a 44-m long water intake gallery at 384 m a.s.l. Although the bed- rock is not visible in the water intake area due to the Qua- ternary cover, it can be assumed that the abstraction site is located where the fl ysch and Palaeogene carbonate deposits come into contact. A large portion of the hinterland of the abstraction site, which reaches the steep slopes of Jurassic rocks at an approximate altitude of 700 m a.s.l., is also cov- ered with thick Quaternary deposits, and so the geological situation is not fully understood. Based on available data for the area upstream of the abstraction site, there appear to be no fl ysch deposits or interruptions in the continuity of the karstic aquifer. The Izbitac spring is captured by a 23-m long water intake gallery located at 20 m a.s.l. The gallery was built along the tectonic contact between the Eocene fl ysch depos- its (E2,3) and the Palaeogene carbonate deposits (Pg). Be- tween the Palaeogene carbonate deposits and the carbonate hinterland, (composed of the Jurassic and Cretaceous car- bonate deposits), where the majority of the spring recharge area is located, there is an upper layer of fl ysch deposits. Al- though tectonically reduced and pinching out to the NW, these layers seem to signifi cantly infl uence the hydrody- namic characteristics of the Izbitac spring. 3. SPRING MEASUREMENTS AND DATA ANALYSIS Measurements of water level and temperature were re- corded for the Makarski Vrutak, Orašje, Podgorski Vrutak and Izbitac springs (using Hobo Water Level Loggers, On- set Computers, USA). In addition to these springs, in the SW foothills of central Mt. Biokovo, the Grebice spring 1 BOJANIĆ, L. (1980): Hidrogeološka studija područja Aržano-Brela-Metković [Hydrogeological Study of Aržano-Brela-Metković area]. Unpubl. report, Croatian Geological Survey, Zagreb, 75 p. Kuhta et al.: Hydrodynamic characteristics of Mt. Biokovo foothill springs in Croatia Geologia Croatica 45 (Fig. 3) is captured for use in the public water supply sys- tem; however, due to the technical characteristics of the water intake structure, it was not possible to establish con- tinuous water level monitoring at this location. 3.1. Spring measurements Measurements were taken by installing square or standard tri- angular weirs in the water collection chamber outlets. Up- stream of the weirs, water level and temperature data loggers were installed and programmed to take hourly measurements. Subsequent processing with appropriate mathematical equa- tions, converted the measured level data into spring discharge fi gures. The monitoring period lasted from 25 November 2007 until 17 December 2008. Direct volumetric measurements were also occasionally made during this period. Despite the millimetre accuracy of the measurement de- vices, the data logger results had certain defi ciencies due to the inability to fulfi l all necessary parameters, including the width and characteristics of the supply canal and the ratio between the weir height and the upstream water depth, at the weir installation sites. These defi ciencies were mostly over- come by correcting the parameters in the applied mathemat- ical formulae for square (KINDSVATER & CARTER, 1959), and triangle weirs (RANTZ et al., 1982). For some springs (especially the Orašje spring), it was not possible to fully calm the water surface upstream of the weir, which caused disturbance (noise) in the results. 3.2. Analysis of hydrograph recessions Methods for analysing the hydrograph recession of individ- ual fl ood pulses were initially used primarily for discharge forecasts in the case of long dry periods, and to determine available groundwater supplies in aquifers (BOUSSINESQ, 1877, 1904; MAILLET, 1905). Afterwards, karst hydroge- ology specifi c methods for recession analysis were devel- oped (DROGUE, 1972; MANGIN, 1975; ATKINSON, 1977; MILANOVIĆ, 1981; BONACCI, 1987, 1993). These methods are used to determine the characteristics of hydrau- lically contrasting aquifer volumes, (highly permeable con- duit networks and low permeability rock masses with fi ssure porosity), as well as to characterise and classify karst aqui- fers based on quantitative values. In this study, the longest available hydrograph recessions for the monitored springs were analysed. Hydrograph reces- sion in base-fl ow conditions is described with Maillet’s ex- ponential equation (MAILLET, 1905), which is probably the most frequent recession equation used in the literature (DE- WANDEL et al., 2003): Qt=Q0e–αt (1) where Qt is the yield (L/s) in time t, Q0 is the yield (L/s) at the beginning of a recession segment (in time t = t0), t is the time (day) from the beginning of recession and a is the base- fl ow recession coeffi cient (1/day), which depends on an aquifer’s characteristics. The total hydrograph recession can be described by ap- plying an exponential equation with one recession coeffi cient only when fast-fl ow spring components are absent. To de- scribe the recession of typical karst springs with marked fast- fl ow and base-fl ow components, several different recession coeffi cients are sometimes used (BONACCI, 1993). The re- cession coeffi cient is gradually decreased in several steps. The value of the highest recession coeffi cient at the begin- ning of the recession diagram is attributed to the drainage characteristics of the aquifer system’s conduit network dur- ing fast-fl ow conditions, whereas the lower recession coef- fi cient, which ends the recession diagram, is attributed to the drainage characteristics of low permeability volumes within a karst aquifer. However, this interpretation is not completely correct due to signifi cant variations in the temporal and spa- tial distribution of infi ltrated rainfall, which impact the hy- grograph’s shape during fast-fl ow conditions (GRASSO & JEANNIN, 1994). This phenomenon explains why there are signifi cant differences in the characteristics of hydrograph recession during fast-fl ow conditions for individual fl ood pulses and why these hydrographs cannot be used to unam- biguously interpret the conduit network characterises of an aquifer. Thus, only the recession components for base-fl ow conditions were isolated and processed in this study. By integrating Maillet’s equation, the dynamic volume V, which is the groundwater volume at the beginning of the recession period (in time t0) that is stored in the low perme- ability medium of an aquifer above spring level, can be de- termined: (2) where the constant c depends on the discharge and the reces- sion coeffi cient units. 3.3. Analysis of the spring hydrograph using auto- correlation and cross-correlation functions Time series analyses in karst hydrogeology are used to de- termine karst system characteristics, to investigate the rela- tionships between spring discharge and rainfall and between discharge, to continuously measure physico-chemical pa- rameters and to assess the interdependence of discharge in several springs (MANGIN, 1984; PADILLA & PULIDO- BOSCH, 1995; LAROCQUE et al. 1998; PANAGOPULOUS & LAMBRAKIS, 2006; KOVAČIČ, 2010). Univariate meth- ods are used to analyse the characteristics and structure of in- dividual time series, while bivariate methods are used to anal- yse the connection between the input and output time series (BOX & JENKINS, 1974). In this study, mean daily discharge series were analysed with the autocorrelation function, while the connection between discharge and rainfall was analysed with the cross-correlation function. The autocorrelation function defi nes the linear depend- ence of successive data values within a time series, depend- ing on their time lag. A time series is compared with itself, Geologia Croatica 65/1Geologia Croatica 46 with a discrete increase of time lag, and the autocorrelation coeffi cient r(k) is calculated for individual time lags k with the following equation: (3) (4) where xi = (x1,…xi,…xN) is a time series of N data points for which m autocorrelation coeffi cients rk = (r0,…rk,…rm) are calculated. The value of m should not exceed N/2 (BOX & JENKINS, 1974) or N/3 (MANGIN, 1984). A graphical pre- sentation of rk values for individual time lags refl ects the character, karstifi cation intensity and retention capacities of an aquifer. MANGIN (1984) terms the time necessary for rk to fall below 0.2 the „memory effect“ which expresses the duration of the system’s response to an input impulse. A lon- ger memory effect indicates a less developed karst conduit network and a higher retention capability of the karst system, while a short-term effect refl ects the system’s low retention capability, with a dense, connected conduit network. The cross-correlation function defi nes the linear depend- ence of output yt and input xt time series, depending on their time lags. The output yt and input xt time series are mutually compared, with a discrete increase of the time lag, and the cross-correlation coeffi cient rxy(k) is calculated for individual time lags k with the following equation: , (5) where: , (6) and σx and σy are the standard deviations of the input and output series, respectively. If the input series is random, then the cross-correlation function rxy(k) presents the system re- sponse to the input impulse. Because rainfall can mostly be considered as a random signal, the cross-correlation function of the rainfall and discharge presents a unit hydrograph. In this study, daily rainfall data from a gauge in Makarska were used. These data represent rainfall at the foot of Mt. Bi- okovo. Although events in the higher parts of the massif are more relevant for aquifer recharge, no daily rainfall gauges have been installed in the higher parts of the massif. 4. RESULTS AND DISCUSSION Spring yield and water temperature data were used to deter- mine the hydrodynamic characteristics of the Makarski Vru- tak, Orašje, Podgorski Vrutak and Izbitac springs. As late summer and early autumn during the monitoring period were dry, exceptionally low water levels and representative values of minimum yield were observed. Conversely, a wet period at the end of 2008 resulted in exceptionally high water lev- els, so relevant maximum yields were also recorded. The spring hydrographs and thermographs (Fig. 4) were created using mean daily values of yield and water temperature dur- ing the monitoring period. The basic statistical data on dis- charge for the investigated springs are shown in Table 1. Due to technical conditions at the water intake of the Makarski Vrutak spring, the installed weir could not reliably measure discharges greater than 200 L/s, which was probably refl ected in the slightly higher occurrence of errors in the maximum yield calculations and the shape of the hydrograph peaks. The more pronounced yield fl uctuations in the hy- drograph for the Orašje spring were caused by an unsettled water surface in the collection basin. The hydrographs (Fig. 4) exhibit the different hydrody- namic characteristics of the investigated springs. Consider- ing the minimum to maximum yield ratios (Table 1), the Makarski Vrutak spring has a very high ratio, typical of karstic aquifers. The ratios are signifi cantly lower for the Podgorski Vrutak and the Orašje springs and very low for the Izbitac spring. The hydrograph of the Makarski Vrutak spring shows abrupt fl ow variations, with a pronounced transition from the fast-fl ow to the base-fl ow regime (Figs. 4a and 5a). Re- cession analysis was conducted for the longest recession pe- riod in the monitoring period, which lasted from the middle of June to the end of October 2008, a total of 142 days. Af- ter the fi rst 40 days of the recession period, the fl ow reces- sion assumed an exponential form with a constant recession coeffi cient (Fig. 5a). The transition to slow exponential re- cession (i.e., to base-fl ow recession) occurred at a low dis- charge rate (12 L/s), which was considerably lower than the spring’s median and mean discharges (Table 1). This result indicates that the discharge of the Makarski Vrutak spring is in the fast-fl ow regime for most of the year. The pronounced discharge fl uctuations in the fast-fl ow regime result from both the rapid infi ltration of precipitation in the epikarst and a well developed system of drainage conduits, which effec- tively convey infi ltrated water to the spring. A fast through- put, with an accompanying low residence time in the system, results in the spring’s high vulnerability to contamination. A relatively low recession coeffi cient during the base-fl ow re- gime refl ects the high retention capability of low permeabil- ity volumes in the aquifer, although their contribution to the total spring discharge is minor. Consequently, the spring’s dynamic volume is lower than those of the Podgorski Vrutak and Izbitac springs (Table 2), despite its higher mean and median discharges (Table 1). The Orašje spring hydrograph is characterised by a clear distinction between the fast-fl ow and base-fl ow regimes and a less dominant fast-fl ow component in the total spring dis- charge, compared to the Makarski Vrutak spring (Figs. 4b and 5b). The analysed recession period lasted for 140 days, of which the fi rst 23 were in the fast-fl ow regime. The base- fl ow recession coeffi cient was the lowest of all investigated springs, but was similar to that of the Izbitac spring. Despite a low recession coeffi cient, the spring system’s dynamic vol- Kuhta et al.: Hydrodynamic characteristics of Mt. Biokovo foothill springs in Croatia Geologia Croatica 47 ume was also the smallest (Table 2) as a consequence of its lowest mean and median discharges (Table 1). The difference between the fast-fl ow and base-fl ow com- ponents of the Podgorski Vrutak spring hydrograph was much less pronounced than those of the Makarski Vrutak and Orašje springs (Figs. 4c and 5c). The base-fl ow was a dom- inant component of the spring’s total discharge. The hy- drograph’s less pronounced changes in the fast-fl ow regime Figure 4: Daily rainfall in Makarska and temperature and discharge values of a) the Makarski Vrutak spring, b) the Orašje spring, c) the Podgorski Vrutak spring and d) the Izbitac spring. Table 1: Basic statistics for daily discharge of the springs during the monitoring period (25 November 2007 – 17 December 2008). Parameter Location Makarski Vrutak Orašje Podgorski Vrutak Izbitac Qmin (L/s) 3.7 1.3 3.7 4.8 Qmax (L/s) 448.2 44.9 73.7 20.9 Qmean (L/s) 112.6 9.9 28.9 12.6 Qstdev (L/s) 129.2 9.0 17.2 4.9 Qmed (L/s) 54.4 6.0 28.9 12.9 Qmin : Qmean : Qmax 0.03 : 1.00 : 3.98 0.13 : 1.00 : 4.54 0.13 : 1.00 : 2.55 0.38 : 1.00 : 1.66 Geologia Croatica 65/1Geologia Croatica 48 could be a result of the lower maximum draining capacity of the system’s conduit network or the signifi cant sediment layer that covers the karst and retards infi ltration. However, the highest base-fl ow recession coeffi cient suggests the ex- istence of a dense conduit network. Because of the fast-fl ow to base-fl ow transition at a relatively high discharge (ap- proximately the magnitude of the mean and median annual discharges, Table 1, Fig. 5c), the system has a proportionally high dynamic volume (Table 2). The characteristics of the Izbitac spring hydrograph are different from those of the other springs. The fast-fl ow dis- charge component is completely absent, and the hydrograph recession can be described with the same exponential rela- tionship from the beginning of the recession period (Figs. 4d and 5d). This, in combination with a low recession coeffi - cient, results in the system’s highest dynamic volume (Table 2). The lack of a fast-fl ow component suggests the absence of a developed karstic conduit network in the spring catch- ment. Given the lithological composition of the spring’s sur- rounding area (Fig. 3), the appearance of fl ysch rocks in the nearby hinterland could also retard the discharge dynamics. Occasional, slight increases in the chloride contents recorded for this spring indicate bathyphreatic recharge conditions, possibly through non-karstifi ed carbonate rocks below the flysch deposits. A good retention capability and longer groundwater residence times make the Izbitac spring the most valuable for purposes of water supply. The time series analysis utilises daily discharge values of the springs and daily rainfall values measured in Makar- ska. The uniform decline of autocorrelation function values for the Izbitac spring discharge series confi rms the presence of only one base-fl ow component, while the presence of both fast-fl ow and base-fl ow components is clearly marked by a transition from the steep to the gentle slope in the shape of Figure 5: Hydrograph recession characteristics for a) the Makarski Vrutak spring, a) the Orašje spring, c) the Podgorski Vrutak spring and d) the Izbitac spring. Table 2: Base-fl ow recession parameters (α – recession coeffi cient; Q0 – base-fl ow discharge at the beginning of recession; V – dynamic aquifer volume). Parameter Location Makarski Vrutak Orašje Podgorski Vrutak Izbitac α (day-1) 0.0111 0.0083 0.0173 0.0085 Q0 (L/s) 16 7 28 20 V (m3) 12.5·104 7.3·104 14.0·104 20.3·104 Kuhta et al.: Hydrodynamic characteristics of Mt. Biokovo foothill springs in Croatia Geologia Croatica 49 the other springs’ functions (Fig. 6a). Consequently, the Iz- bitac spring hydrograph shows the longest memory effect (Fig. 6a, Table 3). The Makarski Vrutak spring’s autocorre- lation function has typically, a two component shape, with a steep slope for the fi rst 25 days and an almost horizontal slope for the next 55 days. Because of fast-fl ow dominance in the Makarski Vrutak discharge dynamics, the transition to an almost horizontal function shape occurs at low autocor- relation values below 0.2 (Fig 6a), resulting in a short mem- ory effect. The Orašje spring’s autocorrelation function shows similar characteristics to that of the Makarski Vrutak spring, with an even lower memory effect (Fig 6a, Table 3). Low memory effect and zero autocorrelation values for the time lags of 24 to 33 days on these springs are largely caused by irregularities in the Orašje discharge measurement data, which is especially pronounced in the low fl ow (i.e., base- flow) regime (Fig. 4). The Podgorski Vrutak function presents an intermediate case between the Izbitac and Makar- ski Vrutak functions, with two distinct fl ow components and a less dominant and more inert fast-fl ow component, com- pared to the Makarski Vrutak function (Fig 6a). The cross-correlation functions show the most rapid re- sponse and the strongest dependence on rainfall for the Oraš- je spring’s discharge series (Fig. 6b, Table 3). The Makarski Vrutak spring’s discharge series had a slightly lower value and a longer time lag of maximum cross-correlation coeffi - cient (Table 3). The general shapes of the functions for both springs are similar, with a rapid fall of values for the fi rst 15 or 16 days, refl ecting the spring’s fast-fl ow response. The low function values then slowly and irregularly approach zero until a time lag of approximately 45 days is reached. The low values following the initial rapid fall of the func- tions indicate the dominance of the fast-fl ow response, whe- reas the approach to zero marks the total duration of the sys- tem response. The cross-correlation function of the Pod gorski Vrutak spring discharge and rainfall shows a longer duration of the fast-fl ow regime (23 days), a lower maximum corre- lation coeffi cient and a signifi cantly longer time lag (Fig. 6b, Table 3). A long time lag of the system response implies a delay of rainfall infi ltration and slower transport through the system, which could be caused by the presence of a sediment karst cover or a less developed conduit network. Again, the function shows a more regular shape and higher coeffi cients after the end of the fast-fl ow regime, which is representative of the Podgorski Vrutak system with its less pronounced dif- ference between the fast-fl ow and base-fl ow responses. The Izbitac spring cross-correlation function shows no fast-fl ow response and a generally low dependence on daily rainfall, with maximum correlation values for time lags of more than 50 days (Fig. 6b). Groundwater temperature data indicate the varying ther- modynamic characteristics of the springs’ aquifer systems (Fig. 4). The Makarski Vrutak spring had the most prono- unced temperature fl uctuations in relation to discharge dy- namics. The total range of the measured temperatures fell within the range 13.6 to 14.1 °C. Flood pulses were followed by temperature drops (Fig. 4a). During November, water waves fi rst occurred after a long recession period, and tem- perature drops were preceded by small increases of short duration. Delayed temperature drops in relation to the start of yield increases during water waves were heterogeneous and lasted, on average, 1 to 2 days. The delays were longer Figure 6: a) Autocorrelation functions of the springs’ daily discharge series and b) cross-correlation functions of the springs’ daily discharge series and daily rainfall series. Table 3: Time series analysis parameters: system memory eff ect, maximum discharge/rainfall cross-correlation coeffi cient and time lag for maximum cross-correlation coeffi cient. Parameter Location Makarski Vrutak Orašje Podgorski Vrutak Izbitac Memory eff ect 23 17 60 70 Crosscorrelation coeffi cient 0.35 0.4 0.27 0.13 Time lag 1.5 1 10 50 Geologia Croatica 65/1Geologia Croatica 50 following longer hydrograph recession periods and shorter during wet periods. Along with sudden temperature changes related to changes in discharge, seasonal fl uctuations are also evident in the spring, with a maximum temperature recorded at the time of low yields in September and the beginning of November. Such thermodynamic characteristics indicate the fast transport of water through the system by well-developed conduits. Temperature increases that precede drops are prob- ably related to the mobilisation of the water contained in low permeability fi ssures due to rapid pressure increases caused by rainfall infi ltration. An alternative explanation of anteced- ent temperature increases is the activation of deep conduit fl ow as the result of an increase in the hydraulic gradient within the system. The temperature of the Orašje spring occasionally chang ed, mostly by dropping, in relation to yield changes, although these changes were small and near the resolution limit of the data loggers (Fig. 4b). As with the Makarski Vru- tak spring, seasonal fl uctuations were evident, with maxi- mum temperatures observed during the summer recession, in the period from July to the end of September. The total range of measured temperatures in this spring was between 12.7 to 13.0 °C, with an average of 12.85 °C. The Podgorski Vrutak and Izbitac springs are character- ised by different, although relatively homogenous water tem- peratures during the entire year (Figs. 4c and 4d). The tem- perature fl uctuations in the Podgorski Vrutak spring ranged from 12.1 to 12.3 °C, while those in the Izbitac spring ranged from 14.9 to 15.0 °C. Such small changes are within the res- olution limits of the instruments and reveal the thermal sta- bility of the recharge system. The stability of the Podgorski Vrutak system could be a consequence of better heat ex- change between rocks and water in conduits with smaller radii, while for the Izbitac spring, the stability is probably the result of the absence of karst conduits in the nearby hin- terland. The minimum water temperature values in the Pod- gorski Vrutak spring were obtained during two short tem- perature drops in January and March. These drops were not temporally coordinated with yield increases; rather, they oc- curred during shorter winter recession periods. The mean annual groundwater temperatures in the ob- served springs ranged from 12.2 to 15.0°C, whereas the mean annual temperatures in the Mt. Biokovo high karst plateau were 6–7 °C. Such temperature differences indicate that the majority of the investigated springs’ local catchments are located in the marginal SW parts of the massif, near the steep slopes along contact between the fl ysch and carbonate de- posits. Water from the major part of the plateau area pre- dominantly drains toward the main discharge areas of sub- marine springs Dupci and Klokun. The Cretaceous and Palaeogene deep basinal deposits (limestones and breccias) on the NW slopes of the Mt. Bi- okovo massif, which according to recent geological investi- gations are considered to be a separate tectono-stratigraphic unit (Kotišina-Tilovica trough), comprise signifi cant hinter- land sections of the Makarski Vrutak and Podgorski Vrutak springs. These deposits haven’t been previously hydrogeo- logically investigated in detail. According to the results of this study, it is signifi cant that in terms of the hydrodynamic characteristics of the springs, they have characteristics of permeable karstic rocks, similar to the platform carbonates of the remaining parts of Mt. Biokovo (High karst, Dinaridic segment). 5. CONCLUSIONS Hydrograph analyses for the Makarski Vrutak, Orašje, Pod- gorski Vrutak and Izbitac springs in the Makarska Coast area indicate the different hydrogeological characteristics of the aquifers in the catchment areas of the springs. The Makarski Vrutak and Orašje springs are typical karst springs. The Pod- gorski Vrutak recession analysis indicates a dense conduit network in the hinterland, but also its lower maximum drain- ing capacity or the signifi cant sediment layer that covers the karst and retards infi ltration. The Kotišina deposits in the hinterland of Makarski Vrutak and Podgorski Vrutak springs have not been previously hydrogeologically analysed, and are assessed to have characteristics of permeable karstic rocks, similar to the platform carbonates of the remaining parts of Mt. Biokovo. The analysis results for Izbitac indi- cate that it is not a typical karst spring. The geological set- tings of the spring’s hinterland (fl ysch deposits) could be the cause of the spring’s dynamics. The results of temperature data analyses additionally confi rm the determined differ- ences in hydrogeological conditions in the recharge areas of the investigated springs. All investigated springs, with the exception of the Iz- bitac spring, are extremely sensitive to surface pollution. However, due to unfavourable topographic conditions that have prevented urbanisation and the fact that the catchment areas of these springs are located within the „Biokovo“ Na- ture Park, there are virtually no anthropogenic impacts on the majority of the sanitary protection zones. Thus, the risk of groundwater pollution can be considered to be low. The specifi c hydrodynamic characteristics of the Izbitac spring indicate a lower sensitivity to pollution in the context of groundwater fl ow velocity and residence time. However, the atypical hydrodynamic behaviour is probably a result of the delaying effect of fl ysch deposits in the spring’s imme- diate hinterland, while the aquifer is karstic and is suscepti- ble to pollution. ACKNOWLEDGEMENT The paper represents part of a hydrogeological investigation within the project entitled the „Basic Hydrogeological Map of the Republic of Croatia” (project no. 181-1811096-3165) which is fi nanced by the Ministry of Science, Education and Sports. The authors wish to express their gratitude to the Ma- karska Water Utility Company for cooperation during fi eld explorations, Nico GOLDSCHEIDER (Karlsruhe Institute of Technology) and an anonymous reviewer for their con- structive comments which substantially improved the sub- mitted manuscript and David DREW (Trinity College Dub- lin) for his linguistic help and useful comments. Kuhta et al.: Hydrodynamic characteristics of Mt. Biokovo foothill springs in Croatia Geologia Croatica 51 REFERENCES ATKINSON, T.C. (1977): Diffuse fl ow and conduit fl ow in limestone terrain in Mendip Hills, Somerset (Great Britain).– J. Hydrol., 35, 93–100. doi:10.1016/0022-1694(77)90079-8 BAKŠIĆ, D. & JALŽIĆ, B. (2001): Jama Amfora [The Amfora pit – in Croatian].– Speleo’zin, 14, 7–9. BAKŠIĆ, D. & LACKOVIĆ, D. (2002): Jama Amfora, –614 m – naj- dublja jama Biokova [The Amfora pit, –614 m – the deepest pit of Mt. Biokovo – in Croatian].– Velebiten, 36, 16–22. BONACCI, O. (1987): Karst Hydrology.– Springer-Verlag, Berlin, 173 p. BONACCI, O. (1993): Karst springs hydrographs as indicators of karst aquifers.– Hydrological Science Journal, 38/1–2, 51–62. BOUSSINESQ, J. (1877): Essai sur la théorie des eaux courantes do mouvement nonpermanent des eaux souterraines.– Acad. Sci. Inst. Fr., 23, 252–260. BOUSSINESQ, J. (1904): Recherches théoriques sur l’écoulement des nappes d’eau infi ltrées dans le sol et sur le débit des sources.– J. Math. Pure Appl., 10, 5–78. BOX, G.E.P. & JENKINS, G.M. (1974): Time Series Analysis: Forecast- ing and Control.– Holden Day, San Francisco, 575 p. CHOROWICZ, J. (1975): Le devenir de la zone de Budva vers le Nord- Ouest de la Yougoslavie.– Bulletin de la Societe Geologique de France, 7/17, 699–709. DEWANDEL, B., LACHASSAGNE, P., BAKALOWICZ, M., WENG, PH. & AL-MALKI, A. (2003): Evaluation of aquifer thickness by analysing recession hydrographs. Application to the Oman ophiolite hard-rock aquifer.– J. Hydrol., 274, 248–269. doi: 10.1016/S0022- 1694(02)00418-3 DROGUE, C. (1972): Analyse statistique des hydrogrammes de decrues des sources karstiques.– J. Hydrol., 15, 49–68. doi:10.1016/0022- 1694(72)90075-3 GAJIĆ-ČAPKA, M., PERČEC TADIĆ, M. & PATARČIĆ, M. (2003): Digitalna godišnja oborinska karta Hrvatske [Digital annual pre- cipitation map of Croatia – in Croatian].– Hrvatski meteorološki časopis, 38, 21–33. GOLDSCHEIDER, N., DREW, D. & WORTHINGTON, S. (2007): In- troduction.– In: GOLDSCHEIDER, N. & DREW, D. (eds.): Meth- ods in Karst Hydrogeology: IAH: International Contributions to Hydrogeology, 26, Taylor and Francis, 1–8. GRANDIĆ, S. (1974): Neke naftnogeološke karatkteristike naslaga Vanj skih Dinarida [Some regional petroleum geological charac- teristics of the External Dinarides deposits – in Croatian].– Nafta, 25/3, 111–120. GRASSO, D.A. & JEANNIN, P.Y. (1994): Etude critique des méthodes d’analyse de la réponse globale des systèmes karstiques. Applica- tion au site de Bure (JU, Suisse).– Bulletin d’Hydrogéologie de l’Université de Neuchâtel, 13, 87–113. JELASKA, V., BENČEK, Đ., CVETKO TEŠOVIĆ, B., ČOSOVIĆ, V., GUŠIĆ, I., IŠTUK, Ž. & MATIČEC, D. (2003): Platform Dynam- ics During the Late Cretaceous and Early Paleogene – External Di- narides, Dalmatia.– In: VLAHOVIĆ, I. & TIŠLJAR, J. (eds.): Eval- uation of Depositional Environments from the Paleozoic to the Quaternary in the Karst Dinarides and Pannomian Basin, 22nd IAS Meeting of Sedimentology, Opatija – September 17–19, 2003, Field Trip Guidebook, 101–107, Zagreb. KINDSVATER, C.E. & CARTER, R.W. (1959): Discharge characteris- tics of rectangular thin-plate weirs.– Transections. American Soci- ety of Civil Engineers, 24/3001. KORBAR, T. (2009): Orogenic evaluation of the External Dinarides in the NE Adriatic region: a model constrained by tectonostratigraphy of Upper Cretaceous to Paleogene carbonates.– Earth-Science Re- views, 96/4, 296–312. doi: 10.1016/j.earscirev.2009.07.004 KORBAR, T., FUČEK, L., GRGASOVIĆ, T., PALENIK, D., NOVAK, V., KOCH, G., BLAGUS, Z., HUĐBER, K. & HODAK, T. (2010): Kompleksna geološka struktura zapadnog Biokova na trasi tunela Sv. Ilija [Complex geological structure of western Biokovo Mt. on the route of Sv. Ilija tunnel – in Croatian].– In: HORVAT, M. (ed.): 4th Croatian Geological Congress, Šibenik – October 14–15, 2010, Abstracts book CD, 392–393, Zagreb. KOVAČIČ, G. (2010): Hydrogeological study of the Malenščica karst spring (SW Slovenia) by means of a time series analysis.– Acta Car- sologica, 39/2, 201–215. LAROCQUE, M., MANGIN, A., RAZACK, M. & BANTON, O. (1998): Contribution of correlation and spectral analyses to the re- gional study of a large karst aquifer (Charente, France).– J. Hydrol., 205, 217–231. doi:10.1016/S0022-1694(97)00155-8 MAILLET, E. (1905): Essais d’hydraulique souterraine et fl uviale.– Li- brairie Sci. Hermann Paris, 218 p. MANGIN, A. (1975): Contribution á l`étude hydrodynamique des aqui- fères karstiques.– PhD These, Université de Dijon, 124 p. MANGIN, A. (1984): Pour une meilleure connaissance des systemes hydrologiques a partir des analyses correlatoire et spectrale.– J. Hy- drol., 67, 25–43. doi:10.1016/0022-1694(84)90230-0 MARINČIĆ, S., MAGAŠ, N., MAJCEN, Ž. & BENČEK, Đ. (1977): Osnovna geološka karta SFRJ 1:100000, list Ploče K 33–35 [Basic geological map of SFRY 1:100000, Omiš sheet – in Croatian].– Inst. geol. istraž. Zagreb (1967–1971). Savezni geol. zavod, Beograd. MILANOVIĆ, P. (1981): Karst Hydrogeology.– Water resources publi- cation, Colorado, 434 p. PADILLA, A., PULIDO-BOSCH, A. & MANGIN, A. (1994): Relative importance of basefl ow and quickfl ow from hydrographs of karst spring.– Ground Water, 32/2, 267–277. PADILLA, A. & PULIDO-BOSCH, A. (1995): Study of hydrographs of karstic aquifers by means of correlation and cross-spectral analy- sis.– J. Hydrol., 168, 73–89. doi: 10.1016/0022-1694(94)02648-U PANAGOPOULOS, G. & LAMBRAKIS, N. (2006): The contribution of time series analysis to the study of the hydrodynamic character- istics of the karst systems: Application on two typical karst aquifers of Greece (Trifi lia, Almyros Crete).– J. Hydrol., 329, 368–376. doi: 10.1016/j.jhydrol.2006.02.023 RANTZ, S.E. et al. (1982): Measurement and computation of stream- fl ow.– Volume 1: Measurement of stage and discharge. U.S. Geo- logical Survey Water Supply Paper no. 2175. SCHMID, S.M., BERNOULLI, D., FÜGENSCHUH, B., MATENCO, L., SCHEFER, S., SCHUSTER, R., TISCHLER, M. & USTASZE- WSKI, K. (2008): The Alpine–Carpathian–Dinaridic orogenic sys- tem: correlation and evolution of tectonic units.– Swiss J. Geosci., 101, 139–183. doi: 10.1007/s00015-008-1247-3 TIŠLJAR, J., VLAHOVIĆ. I., VELIĆ, I. & SOKAČ, B. (2002): Carbon- ate platform megafacies of the Jurassic and Cretaceous deposits of the Karst Dinarides.– Geol. Croat., 55/2, 139–170. VELIĆ, I., VLAHOVIĆ, I. & MATIČEC, D. (2002): Depositional se- quences and palaeogeography of the Adriatic Carbonate Platform.– Mem. Soc. Geol. Ital., 57, 141–151. WILLIAMS, P.W. (1983): The role of the subcutaneous zone in karst hydrology.– J. Hydrol., 61, 45–67. doi:10.1016/0022-1694(83) 90234-2 Manuscript received July 19, 2011 Revised manuscript accepted December 28, 2011 Available online February 25, 2012 Geologia Croatica 65/1Geologia Croatica 52