Next Article in Journal
Does Crop Rotation Enhance Groundwater Health? A Review of the Winter Wheat Fallow Policy in the North China Plain
Previous Article in Journal
The Analysis on Similarity of Spectrum Analysis of Landslide and Bareland through Hyper-Spectrum Image Bands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Intermittent Karst River: The Case of the Čikola River (Dinaric Karst, Croatia)

1
Faculty of Civil Engineering, Architecture and Geodesy, University of Split, Matice hrvatske 15, HR-21000 Split, Croatia
2
Croatian Geological Survey, Sachsova 2, HR-10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Water 2019, 11(11), 2415; https://doi.org/10.3390/w11112415
Submission received: 22 September 2019 / Revised: 12 November 2019 / Accepted: 14 November 2019 / Published: 17 November 2019
(This article belongs to the Section Hydrology)

Abstract

:
Intermittent and ephemeral streams (IRES) are responsible for transporting about half of the water on Earth’s surface. Their hydrological behavior is different in various landscapes. IRES are found more often in karst terrains than in any other regions, as a consequence of strong and direct interaction between groundwater and surface water. This paper presents a hydrogeological and hydrological analysis of the intermittent Čikola River and Spring catchment, which is located in deeply karstified and developed parts of the Dinaric karst in Croatia. Hydrological calculations determined that the catchment area covers approximately 300 km2 and very probably changes in accordance with rapid variations in groundwater level. The karst spring of the Čikola River is a cave, extracted for a public water supply with four drilled extraction wells. The results of the interrelated hydrological and hydrogeological analysis show interesting phenomena from an intermittent karst spring (cave) and its catchment, flowing downstream through a karst polje with several smaller confluences, then entering a karst canyon (where the river sinks during certain periods), and ending in an estuary before contributing to the larger Krka River. The research presented was based on water balance calculations, climatic and hydrological time series analyses, spring pumping tests, and thorough hydrogeological interpretation.

1. Introduction

In intermittent and ephemeral streams (IRES), surface water is absent during certain times of the year. This type of stream drains about half of the world’s land surface [1]. Phenomena of their drying (timing, frequency, and duration) strongly depend on climate change and anthropogenic actions. They are among the most common and most hydrologically dynamic freshwater ecosystems. Their planetary significance far exceeds the level of our knowledge of them. It should be stressed that streams that periodically cease to flow are not restricted to arid regions. Although IRES make up a high proportion of river networks on Earth, they have been understudied for a long period [2]. Datry et al. [3] emphasize: “After years of obscurity, the science of IRES has bloomed recently and it is being recognized that IRES support a unique and high biodiversity, provide essential ecosystem services and are functionally part of river networks and groundwater systems.”
Intermittent streams do not have continuous flowing water year-round. Intermittence creates unique longitudinal dynamics in flow regimes. Generally, they have flowing water periods during the wet season and are dry during hot summer months. These streams are partially supplied by groundwater rising to the surface as stream baseflow. Intermittent streams dry up periodically as ground water levels decline in drier seasons. Ephemeral streams have less flow than intermittent streams. They are normally dry for most of the year and have flowing water for brief periods in response to intensive rainfall. In the literature, there are many different, more or less similar, definitions and names of IRES, for example [1]: (1) Ephemeral; (2) episodic; (3) temporary; (4) intermittent; (5) seasonal; (6) dryland; (7) interrupted; (8) nonperennial; (9) near perennial; (10) torrential; (11) influence, etc.
Flow cessation and drying can be caused by one or more of the following processes [4]: (1) Transmission loss (seepage through porous beds or through small swallow-holes at the karst stream’s bottom), (2) evapotranspiration, (3) downward shifts in water tables caused by natural and/or anthropogenic factors, (4) hillslope runoff recession, (5) freeze-up (cessation of melting or soil freezing), (6) cessation of karst spring discharge, (7) drops in lake level below outlet elevation, and/or (8) cessation of discharge from reservoirs. Flow intermittence is an extreme form of flow variation, which has strong and quick influence on different processes in its catchment, especially ecological.
Specific hydrological, geomorphological, and ecological processes characterize the hydrological regime of any given IRES. Generally, an open stream’s flow regime is characterized by: (1) Magnitude—the total amount of flow at any given time, (2) frequency—how often flow exceeds or is below a given magnitude, (3) duration—how long flow exceeds or is below a given magnitude, (4) predictability—regularity of occurrence of different flow events, and (5) rate of change or flashiness—how quickly flow changes from one magnitude to another. Detailed knowledge of IRES hydrological regimes is a prerequisite to ensure their sustainable development. The most important prerequisite for interdisciplinary analyses of IRES is to understand how alterations in flow regimes impair IRES biodiversity, functions, and services.
After interruption of flow in some IRES, surface waters can be isolated in pools and may continue to flow through the hyporheic zone below the streambed. After some time, this flow may also cease [5]. One of the main obstacles is that each case is individual, having very different environmental—as well as climatological and hydrological–hydrogeological—characteristics. Considerable work is needed to explain reactions of various ecosystems and/or species to fast and drastic hydrological changes, especially in timing and duration of flow cessation. These investigations are of special importance in the case of analyses of IRES in karst regions, where such changes occur more often than in other geological media. Because of their unique features, the karstic IRES represent especially endangered and vulnerable systems.
Due to very high infiltration rates, overland and surface flows are rare in comparison with non-karst terrains [6]. Flow regimes in open streamflows in karst depend mostly upon the interaction between the groundwater and the surface water [7,8]. They are hydraulically connected through numerous karst forms, which facilitate and govern the exchange of water between the surface and subsurface [9]. Karst landforms in the analyzed catchment and along the rivercourse play predominant roles in occurrence and duration of flow cessation. Generally, karstic IRES dry up, with the number of zero-discharge days increasing from the spring towards the outlet.
Open streamflow in karst very often disappears underground multiple times and emerges again in different karst springs, usually under a different name. Sinking, losing, and underground streams are frequent karst phenomena [10]. Bonacci et al. [11] described a case of one river in Dinaric karst, which, over a length of 106 km, switches among losing, sinking, and underground stream sections eight times. The occurrence of IRES in karst terrain is more the rule than the exception. A losing streamflow is a surface stream that contributes water to the karst groundwater system in localized areas. A sinking surface streamflow can be defined as a surface river or stream flowing onto or over karst that then disappears completely underground through a swallow-hole (ponor or sinkhole), and which may or may not rise again and flow as a resurgent surface river or stream. Infiltration from sinking streams into the karst groundwater system is the most rapid form of recharge for carbonate aquifers [12]. Underground or subterranean streamflows are subsurface karst passages that have the main characteristics of open rivers or streams.
Karstic IRES cease flowing and dry at certain stream courses, strongly depending on surface and underground karst features and groundwater level oscillations. Rapid and substantial rising and falling of karst groundwater level (more than 100 m in a few hours) is a characteristic which strongly and quickly influences their hydrological characteristics. The lack of techniques to assess the extent and significance of underground conduit networks, caves, and jamas (vertical or subvertical caves), as well as their interconnections with surface karst forms, presents a major problem in analyses, protection, and management of these systems [13]. Nevertheless, groundwater flow through the rock mass discontinuity networks should not be neglected, especially in regional scale studies, as it often represents a significant portion of the storage—even in karstic areas. To ascertain the vulnerability of water and ecological resources of IRES in karst is much more difficult than in non-karst regions. Karst IRES are among the least-known ecosystems on Earth. Long-term ecological studies and population size estimates in these systems are very rare [14].
It should be stressed that in all of Earth’s karst regions, there are many IRES. Some cases from Dinaric karst are described in the literature [11,15,16,17,18,19,20,21,22,23]. IRES are common to all Mediterranean karst areas, whereas the soil is quite thin and the infiltration is relatively high [15,16,17,18,19,20,21,22,23]. For Doglioni et al. [19], IRES are morphological elements of karstic low-relief areas, characterized by relatively large and flat transects. They occasionally drain runoff, in particular after extraordinary or extreme rainfall events. A karst IRES loses water as it flows downstream, and contributes to the karst groundwater systems in localized areas. These losses can be massive in particular river sections, whereas in others, they are small and difficult or even impossible to observe without performing especially precise measurements [11]. The main processes causing flow intermittence in karst areas are: (1) Transmission loss; (2) cessation of spring discharge; and (3) groundwater decline. It is also important to mention other environments of similar monitoring and time series were researched [24,25,26,27].
This paper presents a hydrogeological and hydrological analysis of the intermittent Čikola River and Spring catchment, which are located in deep and developed parts of the Dinaric karst in Croatia. The objective of this paper is to provide a comprehensive explanation of the hydrogeological and hydrological characteristics of one very characteristic intermittent karst river from the Dinaric karst area. The purpose is to improve the general understanding of flow processes within the analyzed intermittent river on the basis of existing data. We hope that this example will be useful for the broader scientific community dealing with complex interdisciplinary problems of karstic IRES.

2. Study Area: Geology, Hydrogeology, Hydrography, and Geomorphology

The study area is located in the central part of the Dinaric karst in Croatia. Figure 1 shows a location map of the study area, indicating the Čikola River and Spring, main tributaries (Vrba and Mahnitaš), two temporary springs (Velika and Mala Kanjevača—too close to Čikola to be presented in maps), the permanent spring Torak, three hydrological stations (Ružić, Drniš, Ključice), meteorological station (Drniš), two rain gauging stations (Muć, Otavice), a water pumping site, and some major ponors. The Čikola is a typical karst intermittent river, while its main tributaries, Vrba and Mahnitaš, are ephemeral open streamflows. Mahnitaš is a typical torrential stream. The Čikola River, the largest tributary of the Krka River, emerges in the edge of Petrovo polje under the Svilaja Mtn., runs through the town of Drniš, and drains into the Krka River near Torak. Thanks to its remarkable beauty, the pristine environment of the Čikola River Canyon has been protected since 1965.
Due to well-developed and deep karstification, catchment area boundaries of the Čikola River are not precisely defined. According to the analyses and calculations made by Bonacci et al. [28], the catchment area of Čikola at Ružić hydrological station covers approximately 300 km2, and very probably changes in accordance with very fast groundwater level variations.
Dinaric karst is morphologically a typical karst terrain, known for its deep compressive tectonics, especially in the neotectonic phase, which cause intensive reverse faulting, fracturing of rock masses, and, consequentially, deep and very irregular karstification. Such natural conditions practically exclude a groundwater hydraulic modeling approach, especially on a regional scale, and simultaneous hydrological and hydrogeological analysis gives the best results for understanding the stream and its complex connections with the groundwater. Although in certain karst terrains hydraulic modeling is acceptable [29], for the Dinaric karst regional scale studies, it would be completely senseless. Calculation of the Dinaric karst aquifers hydraulic properties is acceptable only as a rough approximation [30], in subregional or local scale studies. The investigation of karst aquifers is often limited to the determination of the geological setting, combined with the research of surface water occurrences, as springs, ponors, surface streams, and preferential underground flow paths [31,32]. From a hydrogeological point of view, this middle and north Dalmatian karst terrain was researched on a regional scale [33,34,35]. Several studies have addressed similar problems in karst areas (especially Dinaric) [11,15,16,17,20,21,36,37]. The data from several unpublished technical reports were used as well.
In most of the Dinaric karst of Croatia, the strike of regional geological structures is NW–SE, the so-called Dinaric strike. Small-scale “Basic geological maps” [38,39] cover geological construction of the research area. In the hydrogeological map (Figure 2), lithological and stratigraphic members were grouped according to their hydrogeological features:
Quaternary deposits—on the surface, mostly alluvial, with intergranular porosity, connected with karst poljes. Permeability depends on clast sizes (usually mixtures of terra rossa clays and silts with clasts of base rocks, mostly carbonates). These deposits are of low relevance for deeper karst groundwater studies, but in Petrovo polje, there is influence on Čikola stream and spring.
Water-tight rock masses—mostly covered by Quaternary deposits of karst in Petrovo polje, various rock types from Permian igneous rocks (diabase), Triassic sandstones, and similar clastites to Pliocene marls and limestone-marls (lake sediments).
Low permeability carbonates—mostly dolomites and dolomitic limestones, which are less karstified. Cretaceous, Jurassic, and Triassic rocks.
Intermediate permeability carbonates—dolomitic limestones and marly limestones, which are tectonically fractured and karstified to a certain extent. Cretaceous, Jurassic, and Triassic rocks;
High permeability carbonates—highly fractured and intensively karstified rocks, mostly from the Upper Cretaceous (rudist limestones) but also most Eocene–Oligocene “Promina” rocks (all except marl beds), and other karstified limestones.
The tectonics of the area is extensive, and the rocks are highly karstified [40,41,42,43,44,45,46]. Due to compressive tectonics dominated by the reverse faulting of geological structures, karstification reaches deep horizons in the underground and is very irregularly distributed. On the other hand, the Middle Dalmatian plateau where Krka and Čikola canyons were cut is characterized by relatively mild tectonics. The Čikola canyon is up to 130 m deep and is steeply cut, being present in a karst-plateau with low tectonic disturbance, built mostly in Eocene–Oligocene “Promina” or Upper Cretaceous limestone rock mass. Promina is a series of mostly thick (meter-sized) beds of carbonates, conglomerates, and marly limestones. These rocks are mostly subhorizontal and are some of the most important factors in the creation of the recent relief. Because these marly beds are less present, the Promina complex has mostly karstic hydrogeological behavior and resembles neighboring Upper Cretaceous rudist limestones. Carbonates within Promina are calcitic and karstified to about the same extent as neighboring limestones.
According to such spread of lithological composition, most of the terrain is covered by pervious karstified carbonate rocks. In Petrovo polje, under the Quaternary deposits, there are practically impervious clastic (marly) rocks (Figure 3).
In Figure 4, there is a detailed cross-section of the Čikola spring zone. The main spring is a cave, speleologically explored several tens of meters. In continuation, there is also a depression with several estavelles (suffosion sinkholes), which are, during high groundwater levels, springs of Čikola. When the spring dries up, these estavelles are either dry or act as ponors—for some shorter period of time, groundwater from alluvial Quaternary sediments of Petrovo polje sinks into the karst underground. As shown, there is an extraction site with four bored wells, which were drilled directly into the cave under the lowest groundwater level and are used for public water supply with a maximal pumping rate of 180 L/s. Extraction from the four close-by wells, out of the same cave which acts as estavelle, and numerous smaller estavelles nearby, makes any attempt of hydraulic modeling impossible, even if the results would be presented within an order of magnitude.

3. Study Area Climate

The climate of the study area is predominantly Mediterranean, with influence of continental and mountain climates at the northern upper part of the catchment. In the Čikola River catchment, there are: (1) A meteorological station, Drniš (1957–2018, missing 1991–1997), and (2) two rain gauging stations, Muć (1991–2018) and Otavice (1998–2018).
The annual rainfall measured at the meteorological station Drniš in the available period (Figure 5a) varied within a very large range, between 1581.8 mm (1976) and 630.1 (2011), with an average value of 1027.1 mm. The linear regression line in Figure 5a shows a decreasing trend. Due to missing data for seven years (1991–1997), there is no utility in testing its statistical significance.
Areal rainfall distribution follows a well-developed morphology. Precipitation events are most abundant at Muć station. Over the 21-year period (1998–2018), when all three rain gauging stations were active, average annual rainfall at Muć was 1261.6 mm, at Otavice 991.1 mm, and at Drniš 951.9 mm. July and August are the driest, while November and December are the wettest months. The rainfalls frequently (few times in the year) reach a greater intensity than 50 mm/h.
The time series of mean annual air temperatures measured at Drniš meteorological station (1957–2018, missing 1991–1996) depicted in Figure 5b indicate that air temperatures varied over a broad range, between 12.2 °C (1976 and 1978) and 14.6 °C (2018), with a mean value of 13.3 °C. Linear and second-order polynomial trend lines show existence of an increasing trend, which is very intensive in the recent period of 22 years (1996–2018). If this increasing trend continues, it could have a strong influence on the Čikola hydrological regime, especially in combination with declines in annual rainfall.
Minimum and maximum measured mean daily air temperatures vary between −9.8 °C and 32.9 °C, respectively. Maximum instantaneous temperature can reach 37 °C. The hottest month is July, and the coldest is January [28].
Previous analyses indicate that in the recent period, annual rainfall has declined, while air temperature has tended to increase. If these patterns continue, water availability and the ecological regime in the intermittent Čikola River and its catchment may be significantly endangered.

4. Čikola Karst River Hydrology

Along the Čikola River’s relatively short watercourse (39 km), there are three hydrological stations (Table 1). Detailed and reliable hydrological analysis is inhibited by different durations of their discharge time series, as well as by gaps in their operation. Hydrological station Ružić, located 6.6 km from the river spring, has had the longest period of measurement (58 complete years). The other two downstream stations, Drniš and Ključice, have had 32 and 28 complete years, respectively.
Figure 6a shows three time series of the mean annual discharges, Qmean (m3/s), measured at the three hydrological stations along the Čikola River course: (1) Ružić (1947–2017, missing 1981, 1991–2002), (2) Drniš (1979–2017, missing 1991–1997), (3) Ključice (1987–2017, missing 1992–1994 and 2008). The figure shows similar hydrological behavior of mean annual discharges at the three hydrological stations along the river course. Over a period of 15 years (2003–2017, missing 2008), when all three hydrological stations were active, average mean annual discharges were the following: (1) QRužić-mean-av = 4.4 m3/s, (2) QDrniš-mean-av = 5.0 m3/s, (3) QKljučice-mean-av = 4.9 m3/s. In all cases, t-tests show no statistically significant differences between average mean annual discharges (see third column in Table 2).
The situation is different in the case of maximum annual discharges at the three hydrological stations. Figure 6b shows three time series of the maximum annual discharges, Qmax (m3/s), measured at the three hydrological stations along the Čikola River course. Over a period of 15 years (2003–2017, missing 2008), when all three hydrological stations were active, average annual discharges were the following: (1) QRužić-max-av = 40.2 m3/s, (2) QDrniš-max-av = 49.4 m3/s, (3) QKljučice-max-av = 53.2 m3/s. In two cases, t-tests show no statistically significant differences between average mean annual discharges. The fourth column of Table 2 presents the results of t-tests for average maximum annual discharges between the three pairs of hydrological stations along the Čikola River over a period of 15 years (2003–2017, missing 2008), when all three hydrological stations were active. The results of the t-tests show that Ružić differed significantly (p < 0.01) from both downstream stations Drniš and Ključice in average maximum annual discharges. The difference between average maximum annual discharges between Drniš and Ključice is not statistically significant.
Although intensive precipitation events are frequent in the Čikola catchment, maximum discharges are not adequately high, which can be explained by the influence of karst features and limited outflow capacity of karst springs [47].
Figure 7 shows average monthly discharges, Qmonth-av (m3/s), measured at the three hydrological stations along the Čikola River course during 18 years in the period 1987–2017 (missing 1991–2002 and 2008). Maximum values of average mean monthly discharges occur in February, and minimum values in August. It is interesting to notice that during January and February, the average mean monthly discharges in Ključice are higher than at the two upstream stations, Ružić and Drniš, indicating that when groundwater levels are high, the karstic catchment of the canyon drains groundwater directly into the stream. In March, April, and December, average mean monthly discharges at Ključice are higher than at Ružić, but lower than at Drniš. In every other month from May to November, average mean monthly discharges at Ključice are lower than at the two upstream stations—river water sinks into the karst rock mass, mostly reappearing at the Torak spring in its estuary.
At Ružić hydrological station (1947–2017, missing 1981, 1991–2002), the number of dry days per year, N, varied over a broad range between N = 0 (1947, 1948, 1951, 1957, 1959, 1963, 1966, 1968, 1969, 1972, 1978, 1979, 1980, 1982, 1984) and N = 184 (2003), with an average value of N = 56.3. At Drniš hydrological station (1979–2017, missing 1991–1997), the number of dry days per year, N, varied over a broad range between N = 0 (1983, 1984, 1986, 1989, 1999) and N = 178 (2006), with an average value of N = 75.9. At Ključice hydrological station (1987–2017, missing 1992–1994 and 2008), the number of dry days per year, N, varied over a broad range between N = 93 (2014) and N = 317 (1989), with an average value of N = 211.9 [28]. Figure 8 shows average number of dry days, Nmonth-av (day), during the months monitored at the three hydrological stations along the Čikola River course over 18 years in the period 1987–2017 (missing 1991–2002, 2008). At Ružić station, drying never occurred in January, February, or March. At Drniš station, drying never occurred in March, April, or May. At the most downstream station, Ključice, flow cessation occurred every month.

5. Relationship between Groundwater and Surface Water

The entrance into the Čikola spring cave is situated at ca. 279 m a.s.l. When groundwater starts to appear at the surface, it slowly fills the depression with estavelles (Figure 3), which start to act as springs as well. At the end of this spring zone is a weir (low dam), and when water starts to spill over it significantly, the river starts to flow. This happens at the elevation of 285 m a.s.l. After that, the small increase in groundwater level in the Čikola spring results in significant increase of discharge at the Ružić hydrological station (Figure 9), and after that, at Drniš and Ključice stations as well. After the point when the Čikola spring zone starts to spring over the weir, all water quantities from Vrba and other smaller springs in the area become negligible.
Over the period 3 February 2009–8 November 2010 (645 days), continuous measurement of groundwater level (GWL) was performed at one of the extraction wells drilled in the Čikola spring/cave (Figure 4). During that period, the minimum GWL was 266.91 m a.s.l. and the maximum value was 290.38 m a.s.l. The relationship between the mean daily groundwater level, GWL (m a.s.l.), and the mean daily discharge, Q (m3/s), measured at Ružić (dark blue color), Drniš (red color), and Ključice (green color) hydrological stations during that period is presented in Figure 9. Emergence of surface flow at Ružić hydrological station occurs when GWL at the extraction site exceeds ~274 m a.s.l. At the two downstream hydrological stations, Drniš and Ključice, surface flow occurs when GWL exceeds ~278.5 and ~284 m a.s.l., respectively (Table 1). Groundwater pumping quantity is similar every day (working only a few hours during the night with quantities below 180 L/s), and its influence on groundwater levels proves to be very low, especially with minimal GWL, and therefore can be neglected. Although the Čikola spring starts to discharge groundwater from the cave at some 279 m a.s.l., the real discharge starts when spring-zone water exceeds the top of a small weir—ca. 285 m a.s.l. River flow measured at all three stations below that level at spring is quite low (Figure 9), especially at the Ključice station. Since the other two stations (Ružić and Drniš) are situated outside of the karstic areas, in Petrovo polje Quaternary sediments that overlay impervious rock masses, it can be presumed that water in the Čikola River during such periods (levels at extraction site between 274 and 285 m a.s.l.) is connected with smaller tributaries, such as Vrba, and groundwater from the alluvial sediment of the polje itself. As the measured discharge in the riverbed in such terms rarely exceeds a few cubic meters per second, the data are probably associated with heavy rains in the polje and/or rapid drops in water level during the beginning of the recession period.
Three time series of turbidity S (NTU), groundwater temperature T (°C), and electrical conductivity EC (μS/cm) were measured using a datalogger (fluorimeter) GGUN-FL during one of the two pumping tests (the second pumping test was shorter and is not presented here) that were done on the Čikola spring (extraction wells) (Figure 10). Turbidity ranged from Smin = 4.53 NTU to Smax = 80.97 NTU, electrical conductivity from ECmin = 391 μS/cm to ECmax = 397 μS/cm, and groundwater temperature from Tmin = 11.01 °C to Tmax = 12.49 °C. The changes in EC and groundwater temperature during the pumping test were rather low, and values are typical for karst groundwater in that area. Variations in turbidity were associated with the turbulent start of the pumping test, but after less than a day, they became very low (up to six NTU) and steady.
The pumping test was conducted using maximal pumping quantity in all four wells, 180 L/s, without disruption (except marked period, Figure 10). The test started during very low GWL, ca. 267 m a.s.l., at the end of a long dry period—the riverbed was completely dry. In the first few days, GWL declined with an almost linear trend, but very slightly. On 25 September, rain started in the spring’s catchment, and in less than one day, pumping became negligible, GWL rose significantly (over 20 m during one day), and the river started to flow. Of course, the pumping test was stopped. During both pumping tests, groundwater temperature was equivalent to the average air temperature in the catchment area, ca. 11–12.5 °C, while EC values approximately 395 μS/cm were typical for Dinaric karst groundwater. EC increases with increasing water hardness [32,48]. The karst groundwater conductivity measured at several karst springs in France ranged between 220 and 470 μS/cm, with a mode of 380 μS/cm. During the pumping test, total drawdown was approximately 5 m, followed by a slight increase in the EC value and T. After the first heavy rain in the catchment area, the spring reacted quickly.
One of the most challenging issues when researching karst systems on a regional scale is the positioning of water divides, i.e., delineation of catchment areas. Probably the most powerful tool for this delineation is a groundwater tracing experiment (see connections proven by tracing in Figure 2). The catchment delineation in such karst systems is especially difficult in intermittent systems, when a spring or river dries out. The catchment area is strongly associated with geological construction, and the influence of topography is minor. As mentioned, the hydrologically determined catchment area of the Čikola River at Ružić station is approximately 300 km2, while, for the previously described geological reasons, it cannot exceed 250 km2. The absence of hydrological stations in higher karst areas of Svilaja Mtn. is probably the reason. It is also important to point out that in significant parts of the Dinaric karst, divides are so-called zonal divides, meaning they change position according to hydrological conditions, groundwater levels. Strong karst springs of neighboring catchments preclude the possibility of increasing the Čikola spring catchment area; only common use of hydrological and hydrogeological approaches can lead to increase in knowledge and improvement in the positioning of groundwater divides (they cannot practically be connected with topography). The groundwater divide in the karst plateau between the Krka and Čikola rivers is probably closer to Čikola, but it changes position significantly according to hydrological conditions (unfortunately, there are no observed GWL data in Svilaja Mtn., nor in mentioned karst plateau; no monitoring boreholes at all). All of this area is usually considered to be a catchment of the Torak spring, but during high GWL, it strongly influences Čikola, as is recorded at the Ključice hydrological station (Figure 6).
Krka canyon is deeper (up to 200 m) than the Čikola canyon (up to 130 m). Both canyons were cut in the same geological structure, mostly Eocene Promina carbonate clastites and Cretaceous carbonates, which are both intensively karstified. Krka’s stream is permanent and much stronger, due to stronger contributories and a much bigger catchment area, while Čikola flows at higher altitude and that is the main reason for more emphasized sinking within its own riverbed after entering karst terrain (canyon), except during periods of very high GWL.
In karst terrains, groundwater and surface water constitute a single dynamic system. Determination of the catchment area and boundaries in karst environment presents an extremely difficult and complex task, which very often remains unsolved. A karst basin represents a complex system of water circulation governed by numerous extremely different surfaces and underground karst forms. Locations, dimensions, and spreading of underground karst phenomena are generally unknown. In such heterogeneous and anisotropic systems, water may circulate in unexpected directions, mainly depending on groundwater level. The interdependency between surface water and groundwater in karst occasionally or permanently changes in time and space. These changes are caused by natural processes, as well as anthropogenic interventions, which have been very intensive in the last hundred years or so and are generally uncontrolled. Due to previous briefly specified reasons, the discovery, explanation, and quantification of water circulation in karst media is not an easy task, and in some cases is even impossible.

6. Conclusions

In the case of an intermittent karst river, the specificity of the karst geological setting strongly influences complex interrelationships between the surface stream and groundwater. Detailed hydrogeological and hydrological monitoring and investigations can help elucidate and better explain these processes. The example of the Čikola River provides good evidence of this fact. An exclusively hydrogeological or exclusively hydrological approach would not be adequate for interpretation, as presented. Further improvements can be accomplished with more detailed monitoring of both river and groundwater, which would include continuous measurements of levels and physical–chemical characteristics, especially during storm events and the end of recession. Delineation of catchments and subcatchments can only be achieved with several iterations among hydrogeologists and hydrologists. Due to ongoing climate changes/variations, the monitoring must be more detailed, considering both spatial and temporal distribution, and must cover hydrological, hydrogeological, and ecological processes in the analyzed area. Previous analyses have revealed strong and dangerous increasing trends of air temperature and decreasing trends of precipitation over the last 20 years. This paper shows that only a multi-methodical approach is possible in terrains as such, where hydraulic modeling is not an option. Nevertheless, within further research, a thorough program of hydrochemical, isotopic, and geophysical methods should be included to achieve better understanding of the system.
To better understand and manage IRES of crucial importance is to bring together hydrologists, hydrogeologists, biogeochemists, ecologists, modelers, environmental economists, social researchers, and stakeholders to form a research network for synthesizing the fragmented, recent knowledge on IRES, improving our understanding of IRES and translating this into science-based, sustainable management of river networks. Drying affects flow continuity through time and flow connectivity across longitudinal, lateral, and vertical dimensions of space.
The conclusions of the presented Čikola River and Spring research are comprised of several points: (1) The spring (cave) has limited discharge capacity, and after it reaches its maximal rate, several supplemental springs occur in the vicinity; (2) the reaction of the spring on rainfall in hinterland is rapid and abundant; (3) the river flows in three different environments—through the karst polje, through the canyon, and in estuary; (4) groundwater and surface water are interrelated; and (5) after appearing in polje, the discharges have to increase (Figure 10) to provide enough water for the river to flow through the canyon.
There are multiple areas in which a lack of knowledge hampers the effective management of intermittent streams [49]: (1) Hydrologic processes, (2) the effects of flow stress on ecological processes, (3) restoring degraded IRES, and (4) assessment of the condition of IRES. For both hydrology and ecology, there is a paucity of detailed scientific data on IRES, especially these located in karst regions. Traditional flow-gauging systems have vastly underestimated the number of rivers and streams flowing intermittently [3]. Unfortunately, for more detailed and more reliable conclusions, the analyzed case requires measures from more hydrological stations and a dense network of piezometers drilled in karst aquifer (and Petrovo karst polje as well) along the Čikola River watercourse, as well as precipitation measuring gauges within the catchment areas. Despite this obvious shortage, the example discussed in this paper is very characteristic for practically all IRES in karst terrains and can usefully help other colleagues dealing with the analyses of karstic IRES.
Larned et al. [4] and Boulton [50] suggested six primary objectives for IRES conservation: (1) Preservation or restoration of their aquatic–terrestrial habitat mosaics, (2) preservation or restoration of their natural intermittency, (3) identification of flow requirements for highly valued species and processes, (4) adoption of conservation measures that explicitly acknowledge the importance of the spatial arrangement of ephemeral streams to the entire drainage network, (5) protection of the contributions of ephemeral streams to local and regional aquifer recharge and water quality, and (6) recognition of the cumulative importance of the brief periods of flow to aquatic and terrestrial ecological processes.
Climate change may exert a strong effect on any IRES, and particularly on the Čikola River. Intensification of global warming will cause earlier initiation or later completion of flow cessation. Additionally, water pumping can influence extension of drying duration. A definite conclusion is that water and ecological resources in this vulnerable karst catchment are under stress. These effects are presently difficult to predict.
It is realistic to expect that the occurrence of spring and river drying in the nearby future will be much more frequent and severe. As a result, conflicts between water supplies and the environment will very likely be stronger and with more dangerous consequences. To help mitigate negative consequences, Dragoni et al. [51] insist on a better understanding of complex interrelations between surface water and groundwater, as well as the influence of future anthropogenic actions. We hope that this paper can help in achieving better results in understanding and management of the extraordinarily valuable and complex systems of karstic IRES.

Author Contributions

O.B. provided conceptualization, hydrological investigation, supervision, and validation, and wrote the original draft. J.T. provided data curation, hydrogeological analysis and conceptualization, and data interpretation, and wrote the original draft. T.R.-B. provided supervision and validation, and wrote the original draft. T.F. provided data curation and visualization, and wrote the original draft.

Funding

The hydrological part of the research received no external funding. The hydrogeological part of the research was partially funded by the Croatian Geological Survey’s fundamental project, Basic Hydrogeological Map of Croatia, and partially by Croatian Waters who funded the pumping tests presented.

Acknowledgments

The authors wish to express their gratitude to all authors of the numerous technical reports used for this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Datry, T.; Bonada, N.; Boulton, A. (Eds.) Intermittent Rivers and Ephemeral Streams-Ecology and Management; Academic Press: London, UK, 2017; ISBN 9780128039045. [Google Scholar]
  2. Datry, T.; Singer, G.; Sauquet, E.; Jorda-Capdevilla, D.; Von Schiller, D.; Subbington, R.; Magand, C.; Pařil, P.; Miliša, M.; Acuña, V.; et al. Science and management of intermittent rivers and ephemeral streams (SMIRES). Res. Ideas Outcomes 2017, 3, e21774. [Google Scholar] [CrossRef]
  3. Datry, T.; Bonada, N.; Boulton, A. General Introduction. In Intermittent Rivers and Ephemeral Streams-Ecology and Management; Datry, T., Bonada, N., Boulton, A., Eds.; Academic Press: London, UK, 2017; pp. 1–20. ISBN 9780128039045. [Google Scholar]
  4. Larned, S.T.; Datry, T.; Arscott, D.B.; Tockner, K. Emerging concepts in temporary-river ecology. Freshw. Biol. 2010, 55, 717–738. [Google Scholar] [CrossRef]
  5. Costigan, K.H.; Kennard, M.J.; Leigh, C.; Sauquet, E.; Datry, T.; Boulton, A.J. Flow regimes in intermittent rivers and ephemeral streams. In Intermittent Rivers and Ephemeral Stream-Ecology and Management; Datry, T., Bonada, N., Boulton, A., Eds.; Academic Press: London, UK, 2017; pp. 51–78. [Google Scholar] [CrossRef]
  6. Bonacci, O. Karst Hydrology with Special Reference to the Dinaric Karst; Springer: Berlin, Germany, 1987. [Google Scholar]
  7. Keshavarzi, M.; Baker, A.; Kelly, B.F.J.; Andersen, M.S. River-groundwater connectivity in a karst system, Wellington, New South Wales, Australia. Hydrogeol. J. 2017, 25, 557–574. [Google Scholar] [CrossRef]
  8. Chu, H.; Wei, J.; Wang, R.; Xin, B. Characterizing the interaction of groundwater and surface water in the karst aquifer of Fangshan, Beijing (China). Hydrogeol. J. 2017, 25, 575–588. [Google Scholar] [CrossRef]
  9. Katz, B.G.; DeHan, R.S.; Hirten, J.J.; Catches, J.S. Interactions between ground water and surface water in the Suwannee river basin, Florida. J. Am. Water Resour. Assoc. 1997, 33, 1237–1254. [Google Scholar] [CrossRef]
  10. Bonacci, O. Surface waters and groundwater in karst. In Karst Aquifers-Characterisation and Engineering; Stevanović, Z., Ed.; Springer: Berlin, Germany, 2015; pp. 149–169. [Google Scholar]
  11. Bonacci, O.; Željković, I.; Galić, A. Karst rivers particularity: An example from Dinaric karst (Croatia/Bosnia and Herzegovina). Environ. Earth Sci. 2013, 70, 963–974. [Google Scholar] [CrossRef]
  12. Hess, J.W.; Wells, S.G.; Quinlan, J.F.; White, W.B. Hydrogeology of the South-Central Kentucky karst. In Karst Hydrology Concepts from the Mammoth Cave Area; White, W.B., White, E.L., Eds.; Van Nostrand Reinhold: New York, NY, USA, 1989; pp. 15–63. [Google Scholar]
  13. Worthington, S.R.H. Diagnostic hydrogeologic characteristics of a karst aquifer (Kentucky, USA). Hydrogeol. J. 2009, 17, 1665–1678. [Google Scholar] [CrossRef]
  14. Jovanović Glavaš, O.; Jalžić, B.; Bilandžija, H. Population density, habitat dynamic and aerial survival of relict cave bivalves from genus Congeria in the Dinaric karst. Int. J. Speleol. 2017, 46, 13–22. [Google Scholar] [CrossRef]
  15. Cucchi, F.; Giorgetti, F.; Marinetti, E.; Kranjc, A. Experiences in monitoring Timavo River (Classical Karst). In Tracer Hydrology 97; Kranjc, A., Ed.; Balkema: Rotterdam, The Netherlands, 1997; pp. 213–218. [Google Scholar]
  16. Cucchi, F.; Luca, Z. Underground Timavo river monitoring (Classical karst). Acta Carsologica 2002, 31, 75–84. [Google Scholar] [CrossRef]
  17. Bonacci, O.; Andrić, I. Sinking karst rivers hydrology: Case of the Lika and Gacka (Croatia). Acta Carsologica 2008, 37, 185–196. [Google Scholar] [CrossRef]
  18. Baily-Comte, V.; Martin, J.B.; Jourde, H.; Screaton, E.J.; Pistre, S.; Langston, A. Water exchange and pressure transfer between conduits and matrix and their influence on hydrodynamics of two karst aquifers with sinking streams. J. Hydrol. 2010, 386, 55–66. [Google Scholar] [CrossRef]
  19. Doglioni, A.; Simeone, V.; Giustolisi, O. The activation of ephemeral streams in karst catchments of semi-arid regions. CATENA 2012, 99, 54–65. [Google Scholar] [CrossRef]
  20. Terzić, J.; Stroj, A.; Frangen, T. Hydrogeologic investigation of karst system properties by common use of diverse methods: A case study of Lička Jesenica springs in Dinaric karst of Croatia. Hydrol. Process. 2012, 26, 3302–3311. [Google Scholar] [CrossRef]
  21. Terzić, J.; Marković, T.; Lukač Reberski, J. Hydrogeological properties of a complex Dinaric karst catchment: Miljacka Spring case study. Environ. Earth Sci. 2014, 72, 1129–1142. [Google Scholar] [CrossRef]
  22. Vilenica, M.; Stanković, V.M.; Sartori, M.; Kučinić, M.; Mihaljević, Z. Environmental factors affecting mayfly assemblages in tufa-depositing habitats of the Dinaric Karst. Knowl. Manag. Aquat. Ecosyst. 2017, 418, 14. [Google Scholar] [CrossRef]
  23. Bonacci, O.; Oštrić, M.; Roje-Bonacci, T. Water resources analysis of the Rječina karst spring and river (Dinaric karst). Acta Carsologica 2018, 47, 123–137. [Google Scholar] [CrossRef]
  24. Turconi, L.; Coviello, V.; Arattano, M.; Savio, G.; Tropeano, D. Monitoring mud-flows for investigative and warning purposes: The instrumented catchment of Rio Marderello (North–Western Italy). In Engineering Geology for Society and Territory; Lollino, G., Arattano, M., Rinaldi, M., Giustolisi, O., Marechal, J.C., Grant, G.E., Eds.; Springer: Basel, Switzerland, 2015; pp. 85–90. [Google Scholar]
  25. Lega, M.; Casazza, M.; Turconi, L.; Luino, F.; Tropeano, D.; Savio, G.; Ulgiati, S.; Endreny, T. Environmental Data Acquisition, Elaboration and Integration: Preliminary Application to a Vulnerable Mountain Landscape and Village (Novalesa, NW Italy). Engineering 2018, 4, 635–642. [Google Scholar] [CrossRef]
  26. Cristino, M.C.; Doglioni, A.; Simeone, V. Notes on morphological characters of ephemeral streams in Apulia region, south Italy. In Engineering Geology for Society and Territory; Volume 3, River Basins, Reservoir Sedimentation and Water Resources; Lollino, G., Arrattano, M., Rinaldi, M., Giustolisi, O., Marechal, J.-C., Grant, G.E., Eds.; Springer: New York, NY, USA, 2014. [Google Scholar]
  27. Cristino, M.C.; Doglioni, A.; Simeone, V. Preliminary results about morphological characteristics of ephemeral stream catchments in Apulia region, South Italy. Rend Online Soc. Geol. Ital. 2013, 24, 70–72. [Google Scholar]
  28. Bonacci, O.; Terzić, J.; Roje-Bonacci, T. Hidrološka analiza krške rijeke Čikole. Hrvat. Vode 2018, 26, 281–291. (In Croatian) [Google Scholar]
  29. Fiorillo, F.; Malik, P. Hydraulic Behavior of Karst Aquifers (editorial). Water 2019, 11, 1563. [Google Scholar] [CrossRef]
  30. Terzić, J.; Šumanovac, F.; Buljan, R. An assessment of hydrogeological parameters on the karstic island of Dugi otok, Croatia. J. Hydrol. 2007, 343, 29–42. [Google Scholar] [CrossRef]
  31. Bakalowicz, M. Karst groundwater: A challenge for new resources. Hydrogeol. J. 2005, 13, 148–160. [Google Scholar] [CrossRef]
  32. Ford, D.; Williams, P. Karst Hydrogeology and Geomorphology; Wiley and Sons: Chichester, UK, 2007. [Google Scholar]
  33. Fritz, F. Utjecaj recentnog okršavanja na zahvaćanje voda (Influence of recent karstification on water extraction). Geološki Vjesn. 1991, 44, 281–288. [Google Scholar]
  34. Fritz, F.; Pavičić, A.; Renić, A. Hydrogeological study of the Trogir-Šivenik. Drniš-Knin region, croatian geological survey. Zagreb Unpublished work. 1984. (In Croatian) [Google Scholar]
  35. Kapelj, S.; Kapelj, J. Hydrogeological Exploratory Works for the Sanitary Protection Zones of the Čikola Spring; Technical Report; Geotehnički fakultet: Varaždin, Croatia, 2009. (In Croatian) [Google Scholar]
  36. Fritz, F.; Pavičić, A. Hydrogeologically “hanging” parts of the Zrmanja and Krka rivers. In Proceedings of the “VII jugoslavenski simpozij o hidrogeologiji i inženjerskoj geologiji”, Novi Sad, Serbia, 22–26 June 1982; pp. 115–121. (In Croatian). [Google Scholar]
  37. Fritz, F.; Renić, A.; Pavičić, A. Hydrogeology of the Hinterland of Šibenik and Trogir, Croatia. Geol. Croat. 1993, 46, 291–306. [Google Scholar]
  38. Velić, I.; Vlahović, I. (Eds.) Tumač geološke karte Republike Hrvatske 1:300,000 (Guidebook for the basic geological map of the Republic of Croatia, 1:300,000); Croatian Geological Survey: Zagreb, Croatia, 2009; ISBN 978-953-6907-10-6. [Google Scholar]
  39. Ivanović, A.; Sikirica, V.; Marković, S.; Sakač, K. Basic Geological Map, Scale 1:100,000, Drniš Sheet; Croatian Geological Survey: Zagreb, Croatia; Federal Geological Survey of SFRJ: Belgrade, Serbia, 1977. [Google Scholar]
  40. Herak, M.; Bahun, S.; Magdalenić, A. Pozitivni i negativni utjecaji na razvoj krša u Hrvatskoj (Positive and negative influences on the development of the Karst. Krš Jugosl., akad. znan. i umjetn. Zagreb 1969, 6, 45–78. (In Croatian) [Google Scholar]
  41. Pamić, J.; Gušić, I.; Jelaska, V. Geodynamic evolution of the Central Dinarides. Tectonophysics 1998, 297, 251–268. [Google Scholar] [CrossRef]
  42. Prelogović, E.; Aljinović, B.; Bahun, S. New data on structural relationships in the Northern Dalmatian Dinaride Area. Geol. Croat. 1995, 48, 167–176. [Google Scholar]
  43. Vlahović, I.; Tišljar, J.; Velić, I.; Matičec, D. Evolution of the Adriatic carbonate platform: Palaeogeography, main events and depositional dynamics. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2005, 220, 333–360. [Google Scholar] [CrossRef]
  44. Tari, V. Evolution of the Northern and Western Dinarides: A Tectonostratigraphic Approach; Stephan Mueller Special Publication Series; European Geosciences Union: Munich, Germany, 2002; Volume 1, pp. 223–236. [Google Scholar]
  45. Schmid, S.M.; Bernoulli, D.; Fügenschuh, B.; Matenco, L.; Schefer, S.; Schuster, R.; Tischler, M.; Ustaszewski, K. The Alpine-Carpathian-Dinaridic orogenic system: Correlation and evolution of tectonic units. Swiss J. Geosci. 2008, 101, 139–183. [Google Scholar] [CrossRef]
  46. Španiček, J.; Ćosović, V.; Mrinjek, E.; Vlahović, I. Early Eocene evolution of carbonate depositional environments recorded in the Čikola Canyon (North Dalmatian Foreland Basin, Croatia). Geol. Croat. 2017, 70, 11–25. [Google Scholar] [CrossRef]
  47. Barberá, J.A.; Andreo, B. Hydrogeological processes in a fluviokarstic area inferred from the analyses of natural hydrogeochemical tracers: The case study of eastern Serranía de Ronda (S Spain). J. Hydrol. 2015, 523, 500–514. [Google Scholar] [CrossRef]
  48. Bakalowicz, M.; Mangin, A. L’aquifère karstique. Sa définition, ses caractéristiques et son identification. Mémoires Hydrogéologique Série Société Géologique de France 1980, 11, 71–79. [Google Scholar]
  49. Bond, N.; Cottingham, P. Ecology and Hydrology of Temporary Streams: Implications for Sustainable Water Management; eWater Technical Report; eWater Cooperative Research Centre: Canberra, Australia, 2008. [Google Scholar]
  50. Boulton, A.J. Conservation of ephemeral streams and their ecosystem services: What are we missing? Aquat. Conserv. Mar. Freshw. Ecosyst. 2014, 24, 733–738. [Google Scholar] [CrossRef]
  51. Dragoni, W.; Giontella, C.; Melillo, M.; Cambi, C.; Di Matteo, L.; Valigi, D. Possible response of two water systems in Central Italy to climatic changes. In Advances in Watershed Hydrology; Moramarco, T., Barbetta, S., Brocca, L., Eds.; Water Resources Publications: HIghhlands Ranch, CO, USA, 2015; pp. 397–424. [Google Scholar]
Figure 1. Research area position map with most important features.
Figure 1. Research area position map with most important features.
Water 11 02415 g001
Figure 2. Hydrogeological map of the research area.
Figure 2. Hydrogeological map of the research area.
Water 11 02415 g002
Figure 3. Hydrogeological and geomorphological profile of the Čikola riverbed (a) and the spring’s catchment ((b) modified after [35]).
Figure 3. Hydrogeological and geomorphological profile of the Čikola riverbed (a) and the spring’s catchment ((b) modified after [35]).
Water 11 02415 g003
Figure 4. Schematized hydrogeological profile of the Čikola spring site.
Figure 4. Schematized hydrogeological profile of the Čikola spring site.
Water 11 02415 g004
Figure 5. Time series of annual precipitation (a) and mean annual air temperature (b) measured at the meteorological station Drniš over the period 1957–2018 (1991–1997 missing data due to war activity).
Figure 5. Time series of annual precipitation (a) and mean annual air temperature (b) measured at the meteorological station Drniš over the period 1957–2018 (1991–1997 missing data due to war activity).
Water 11 02415 g005aWater 11 02415 g005b
Figure 6. Time series of the mean annual discharges (a) and maximum annual discharges (b) measured at the three hydrological stations along the Čikola River course: (1) Ružić (1941–2017—missing 1981, 1991–2002), (2) Drniš (1970–2017—missing 1991–1997), (3) Ključice (1987–2017—missing 1992–1994, 2008).
Figure 6. Time series of the mean annual discharges (a) and maximum annual discharges (b) measured at the three hydrological stations along the Čikola River course: (1) Ružić (1941–2017—missing 1981, 1991–2002), (2) Drniš (1970–2017—missing 1991–1997), (3) Ključice (1987–2017—missing 1992–1994, 2008).
Water 11 02415 g006aWater 11 02415 g006b
Figure 7. Average monthly discharges measured at the three hydrological stations along the Čikola River course over the period 1987–2017 (missing 1991–2002, 2008).
Figure 7. Average monthly discharges measured at the three hydrological stations along the Čikola River course over the period 1987–2017 (missing 1991–2002, 2008).
Water 11 02415 g007
Figure 8. Average dry days, N, in the months monitored at the three hydrological stations along the Čikola River course over the period 1987–2017 (missing 1991–2002, 2008).
Figure 8. Average dry days, N, in the months monitored at the three hydrological stations along the Čikola River course over the period 1987–2017 (missing 1991–2002, 2008).
Water 11 02415 g008
Figure 9. Relationship between the mean daily groundwater level, GWL (m a.s.l.), and the mean daily discharges, Q (m3/s), measured at Ružić, Drniš, and Ključice hydrological stations over the period 3 February 2009–8 November 2010.
Figure 9. Relationship between the mean daily groundwater level, GWL (m a.s.l.), and the mean daily discharges, Q (m3/s), measured at Ružić, Drniš, and Ključice hydrological stations over the period 3 February 2009–8 November 2010.
Water 11 02415 g009
Figure 10. Pumping test in 2009: Time series of turbidity S (NTU), groundwater temperature T (°C), and electrical conductivity EC (μS/cm) measured every 2 min during the period 16 September 2010 at 6 h 23 min 49 s–25 September 2010 at 6 h 52 min 24 s (868 data). * Period when datalogger was out of water due to technical problems. Groundwater levels on separate axis (pumping rate constant, 180 L/s), red arrow shows the end of the pumping.
Figure 10. Pumping test in 2009: Time series of turbidity S (NTU), groundwater temperature T (°C), and electrical conductivity EC (μS/cm) measured every 2 min during the period 16 September 2010 at 6 h 23 min 49 s–25 September 2010 at 6 h 52 min 24 s (868 data). * Period when datalogger was out of water due to technical problems. Groundwater levels on separate axis (pumping rate constant, 180 L/s), red arrow shows the end of the pumping.
Water 11 02415 g010
Table 1. Main characteristics of the three analyzed hydrological stations along the Čikola River.
Table 1. Main characteristics of the three analyzed hydrological stations along the Čikola River.
Station NameDistance from the Spring—L (km)Elevation (m a.s.l.)Available PeriodMissing Data
Ružić6.6273.7811947–20171981; 1991–2002
Drniš19.0260.6001979–20171991–1997
Ključice30.581.3101987–20171992–1994; 2008
Table 2. Results of t-tests for average mean and maximum annual discharges between the three pairs of hydrological stations during common periods of their operation.
Table 2. Results of t-tests for average mean and maximum annual discharges between the three pairs of hydrological stations during common periods of their operation.
Pair of StationsCommon Period of WorkMean DischargesMaximum Discharges
Ružić–Drniš26 years (1979–2017, missing 1981 and 1991–2002)p-value = 0.330p-value < 0.01
Ružić–Ključice18 years (1987–2017, missing 1991–2002 and 2008)p-value = 0.599p-value < 0.01
Drniš–Ključice23 years (1987–2017, missing 1991–1997 and 2008)p-value = 0.728p-value = 0.222

Share and Cite

MDPI and ACS Style

Bonacci, O.; Terzić, J.; Roje-Bonacci, T.; Frangen, T. An Intermittent Karst River: The Case of the Čikola River (Dinaric Karst, Croatia). Water 2019, 11, 2415. https://doi.org/10.3390/w11112415

AMA Style

Bonacci O, Terzić J, Roje-Bonacci T, Frangen T. An Intermittent Karst River: The Case of the Čikola River (Dinaric Karst, Croatia). Water. 2019; 11(11):2415. https://doi.org/10.3390/w11112415

Chicago/Turabian Style

Bonacci, Ognjen, Josip Terzić, Tanja Roje-Bonacci, and Tihomir Frangen. 2019. "An Intermittent Karst River: The Case of the Čikola River (Dinaric Karst, Croatia)" Water 11, no. 11: 2415. https://doi.org/10.3390/w11112415

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

Article Metrics

Back to TopTop