ACTA CARSOLOGICA XXVI/2 32 373-385 LJUBLJANA 1997 SPRING CHEMOGRAPH ANALYSIS - THE INFLUENCE OF THAW EFFECT AND DISPERSED POLLUTION IMPULSES (CRACOW-CZESTOCHOWA UPLAND, POLAND) ANALIZA KEMOGRAMOV IZVIROV - VPLIV EFEKTA TALJENJA SNEGA IN IMPULZOV RAZPRŠENEGA ONESNAŽEVANJA (KRAKOWSKO-CZESTOCHOWSKO VIŠAVJE, POLJSKA) ANDRZEJ TYC' Izvleček UDK 556.33/.34(438) Andrzej Tyc: Analiza kemogramov izvirov - vpliv efekta taljenja snega in impulzov razpršenega onesnaževanja (Krakowsko-Czestochowsko višavje, Poljska) Napajanje, tok in vskladiščenje podzemne vode so v karbonatnih vodonosnikih zelo spremenljivi. Skladno z različnimi stopnjami zakraselosti je napajanje v razponu od koncentriranega do razpršenega. Članek predstavlja rezultate proučevanja onesnaženja glede na različne načine napajanja v karbonatnih vodonosnikih, s pomočjo monitoringa temperature vode in koncentracij glavnih ionov. Na podlagi analize kemogramov kraških izvirov in ocene tradicionalnih modelov napajanja, avtor predlaga nov model. Taljenja snega se odraža v kemizmu voda SO/" in NO^", celokupna trdota, specifična električna prevodnost). Ključne besede: hidrologija krasa, kraški izvir, kemograf, Poljska. Abstract UDC 556.33/.34(438) Andrzej iyc: Spring cheniograph analysis - the influence of thaw effect and dispersed pollution impulses (Cracow-Czestochowa Upland, Poland) Groundwater recharge flow and subsurface water storage can be extremely variable in carbonate aquifers. Owing to different degrees of karstification, recharge can range from concentrated to dispersed. The paper outlines the results of an investigation of pollution in respect of different recharge situations in carbonate aquifers including of monitoring of water temperature and ion concentration. According to the chemograph analysis, the author suggests modifications to the traditional models. The thaw effect is demonstrated by water chemistry (Ca^+, SO/", NO^", total hardness, specific conductivity). Key words: karst hydrology, karst spring, chemograph, Poland. ' Department of Geomorphology, University of Silesia, ul. B^dzinska 60, 41-200 SOSNOWIEC, POLAND INTRODUCTION For several decades, karst hydrology has used different indirect methods to investigate groundwater systems in karst aquifers. These have included analysis of karst spring hydrographs and chemographs (Mangin 1975; Bakalowicz 1979; Bonacci 1987; Ford & Williams 1989; Bonacci 1993). The analysis of a spring runoff hydrograph during individual flushes has successfully been applied in hydrology, in respect that it enables us to classify karstic aquifers according to the degree of karstification they exhibit and to the state of development of the groundwater circulation system. With the general growth of hydrochemical investigations in karst areas, this kind of analysis has advanced to a stage where detailed spring chemographs have been drawn up for short spring flushes (Jakucs 1959; Bonacci 1987). Certainly, this method of indirect hydro-logical investigation has gained many advocates in recent times and its popularity is, no doubt, due to the increasing availability of automatic measuring stations in the field and the computer in the laboratory. Spring chemographs may be used in the cases of both a single flush and in the analysis of seasonal variability of spring hydrochemical conditions (see also, e.g. White & Stellmack 1968; Shuster & White 1971; Tyc 1992; Reeder & Day 1993; Tyc 1994; Lackey & Krothe 1996). Research into changes of karstic water chemistry where influenced by anthropogenic pollution, which have been carried out in several karst areas, has provided much new data for the interpretation of spring chemographs. On the one hand, pollution represents a valuable chemical marker, whereas, on the other, it makes the evaluation of karst aquifer activity even more complicated than it usually is. In cold-temperate karst areas, where there is a regular annual snow cover, thaw flushes are an important factor in influencing the hydrological and hydrochemical regime of a karst aquifer. With respect to the author's own observations carried out in the central part of the Cracow-Czestochowa Upland, and by analysing the results of the research from other areas (e.g. Reeder & Day, 1993; Lakey & Krothe, 1996), it became apparent that the generally accepted interpretation of karst spring chemographs should be modified for those areas where thaw flushes occur and where the groundwater environment is considerably polluted. Classical chemograph analyses Long-term research in different karst areas all over the World shows that there is a positive relationship between spring regimes, the types of circulation in the aquifer and the temporal variability of some physical features and chemical composition of spring water. The theoretical basis for the application of spring hydrographs analysis gave L. Jakucs (1959) (who investigated springs in Kolmos, Hungary); W.B. White & J.A. Stellmack (1968), W.B. White (1969), and E.T. Shuster & W.B. White (1971) (the Nittany Valley, Pennsylvania); R. Gospodaric & P. Habič (1976) (the Dynaric karst); S.R.H. Worthington (1991) (the Rocky Mountains of the U.S.A.) and, particularly, A. Mangin (1975), M. Bakalowicz (1979) and M. Bakalowicz & A. Mangin (1980) (the Aliou, Baget and Fontestorbes systems in the Pyrhenees, France). Recently, the latter investigation has been enlarged by A. Pulido-Bosch in the Torcal de Antequera System of Spain (e.g. Padilla et al. 1995). These authors all assume that karst aquifers are heterogeneous and represent all possible combinations and organisational stages of runoff which can range between the diffuse flow and conduit flow systems. The hydrological and hydrochemical regimes of springs (spring yields, temperature and specific conductivity and the temporal variability of the chemical properties of spring water) are in a direct relationship with the structure and condition of water circulation in the karst aquifer. This regime is modified by the hmiting conditions of the aquifer. The variability of karst spring discharge relates to the variabihty of precipitation, the intensity of thawing and the volume of aquifer recharge by allochthonous surface streams. Therefore, the hydrological and hydrochemical regime of springs is mainly dependent on the recharge conditions. According to L. Jakucs (1959), it is important to recognise the relationship between the nature of the karst aquifer recharge and the spring response; springs of small variability are associated with autochthonous recharge (percolation and infiltration) whereas there is another type represented by springs of great variability, where, usually, allochthonous recharge by surface streams is a strong influence. According to S.R.H. Worthington (1991), the location of a spring at a particular altitude in relation to neighbouring springs which drain the same karst system, may also be an important factor in influencing variability of discharge. From these considerations, it follows that any classification of springs in a karst area must contain the foUowing: (1) full-flow springs, which drain the whole aquifer (2) underflow springs, which mainly drain basic underground runoff (3) overflow springs, which mainly drain flush runoff, and (4) underflow-overflow springs, which represent a transitional type. Springs which drain a whole aquifer are very rare in carbonate karst areas; they normally occur only in those areas where there are no overflow springs and the single spring is then situated at the lowest topographical point in the whole catchment. S.R.H. Worthington (1991) also assumed that each overflow spring must have a complementary underflow spring; however, in as much as these are seldom recorded in any survey of a karst area, they are often difficult to find. In such cases, the basal runoff goes beyond the aquifer, without surface discharge (therefore it is not possible to measure it). Thus, it is the overflow springs which tend to be the most frequently investigated (and, also, supposedly, the commonest) in carbonate karst terrains. The best-known and commonly apphed analysis of karst spring chemograph Rain [mm] Q [m'/s] C [mS/m] CaCOs [mg/l] Turbidity [mg/l] l)i.spl5o>ö>a)oj Q [l/s] 19,0 1S.0 17,0 16,0 15,0 14,0 13,0 12,0 11,0 lO.O 9 9 9 5) o) 5 s 8 8 $80 cj) O) O) a> g g ro uo ? 8 i i i O) ö> O) ' 0> 0> Fig. 4: Variability of specific conductivity (C rnS/rn) and nitrates (NO^^' mg/l) in water of Klucze II (a), Klucze I (B) and Ryczowek II (C) springs during thaw season of 1991. Q-mean daily discharge in Ryczowek II. In the karst areas which are not influenced by human activity, an increase in the yield is accompanied by a decrease of TH values and mineralisation. In the case of the Cracow-Czestochowa Upland, the spring chemographs are more complicated. Maximum discharge takes place during thaws and, at these times, there is an obvious increase in the values for specific conductivity, total hardness and ionic content which come mainly from pollutants. Detailed measurements during the thaw period were made at the Ryczowek II (a spring which has a fairly consistent flow), Klucze I (one with variable flows) and Klucze II (one which has consistent all-year-round flows) in the thaw period of 1991 (Fig. 4). This chemograph showed that the accepted view concerning changes of mineralisation and the amount of dissolved and salts during thaws is not wholly acceptable. The traditional view holds that, only in the case of Summer precipitation recharge does meteoric water infiltrate into the karst aquifer quickly. However, this does not adequately relate to the minerahsation and chemical composition of spring water which comes from thaw or thaw-rain recharge, which is an important element of spring discharges in such situations. This may be explained by considering in detail the relationship between the runoff hydrograph and the zones of the carbonate aquifer which are drained during the individual phases of a flush. During thaws, when a spring flush is preceded by the long contact of the water (i.e. that which comes from melting snow covers) and the ground, the increase of minerlaisa-tion and total hardness is associated in the spring recharge with water from the the epikarst zone and the upper part of the vadose zone (this is the converse of previously accepted interpretation of spring chemograph). In the early phases of thawing (and also during winter thaws), water which migrates through the epikarst zone has a specific chemical properties, which, later, influences the "thaw effect" in water which discharges from the springs. This is well illustrated by the analyses of the soil water and water from weathered zone below soil cover at Ksi^ža Gora, in Ryczowek (456.2 m a.s.l.), which forms the recharge zone of the Ryczowek II Spring. The water sampled directly after the melting of an unusually thin snow cover here in February 1990 shows physical and chemical properties which are quite different from those exhibited by water which percolates through the soils and remains in the epikarst zone in a late spring (Tyc 1994, 1996) (Table 1). In case of the water sampled in Fubruary 1990, very large levels of Ca-+ (5.5 mv/1), total hardness (5.9 mv/1) and mineralisation (C = 69.1 mS/m) were recorded; these are all considerably larger than the maximal values noted during the whole study period at the Ryczowek Springs, whereas the content of hydrocarbonates in the water studied was small (0.4 mv/1) and the contents of other anions very large (SO/-, 2-4 mv/1; NO^', 1.1 mv/1). These are much larger than the values obtained from the Upper Jurassic Aquifer groundwater in other periods, whether before or after thaws. It is concluded, therefore, that content and total hardness of the epikarst water in these periods are more closely associa- C [mS/m] TH [mg/11 Ca^^ fmg/l] [mg/11 HCO3 fmg/ll SO/ fmg/n CI [mg/ll NO,-rmg/11 THAWING SEASON Snow 5,3 8,0 8,0 0,0 17,1 11,0 2,5 3,2 Epikarstic zone (dispersed flow) 69,1 118,0 110,0 4,8 24,4 194,8 14,0 69,3 Spring water (fed by dispersed flmv) 50,8 108,0 100,0 4,8 210,5 60,5 16,8 31,5 LATE SUMMER Rain water 6,0 8,4 6,8 1,0 3,1 10,1 3,5 5,7 Epikarstic zone (dispersed flaw) 31,0 54,0 48,0 3,6 48,2 75,4 3,5 23,9 Spring water (fed by- dispersed flmi<) 47,4 92,0 88,8 1,92 208,6 32,2 13,0 18,3 Table 1: Chemical comosition of water in recharge and discharge zones of Ryczowek II spring. Water samples represent thawing and late summer seasons of 1990 (after A. Tyc, 1996). ted with sulphates and nitrates than with the CaCO, forthcoming from rock dissolution. This, in turn, shows that the spring water of a thaw flush maximum contains mainly and Mg^^ salts (probably sulphates and nitrates) which are associated with processes taking place in the soils and in the epikarst zone, and which are not the results of solution processes. One result of the research described is the realisation that there is an artificial variabihty of the total hardness in the hydrological cycle. The natural changes in total hardness, which are influenced by environmental factors associated with solution of carbonates, overlap with the changes influenced by the cycles of pollution supplied to the aquifer. In a cold-temperate climatic zone, spring thaws play an important role. Fig. 5 shows the changes of total hardness values (TH) against a background of pH for the Ryczowek Spring.This shows that the pH-TH relationship takes two forms: (i) the pH and the TH both increase or decrease concomitantly, (ii) a pH decrease is accompanied by a TH increase. The first form is typical of groundwater in carbonate rocks. The reaction of the water reflects the level of saturation of dissolved carbonates; an increase of CaCO^ reflects a more alkaline environment and vice versa. The second type of relationship shows that there are cases in the hydrological cycle when a decrease in the water reaction is demonstrably associated with an increase in the total hardness. This is probably caused by an increase in the amount of pollution in the groundwater and the content of calcium salts which come from outside the local H^O—»CO^^CaCO^ system. This phenomenon Fig. 5: Variability of pH and total hardness (TH) in water in Ryczowek II spring in the period of 1989-1990. occurs during the runoff of the water from thaw recharge (and also during the Winters). This is also an element of the so-called "thaw effect" where human influence on the karst system is considerable. This research suggests that the natural hydrochemical regime of the groundwater in the karst of the Cracow-Czestochowa Upland is influenced by human impact during periods of thaw and thaw-rain recharges in the epikarst zone and in associated soil covers (the "thaw effect"). This is illustrated in the chemograph of the thaw period for 1991 and the third week of December 1988, when a sudden thaw took place (thereby increasing the values of some chemical components in most of the springs studied (see Fig. 3). CONCLUSIONS The analysis of spring chemographs in the area of the Cracow-Czestochowa Upland shows that previously held views on changes of the dissolved calcium and magnesium salts during thaw discharges should be modified. This view seems also to be true in the case where spring recharge by summer rainwater takes place following meteoric water infiltrates relatively quickly into the karst aquifer and when the conduit system of water circulation predominates. However, it does not apply in the case where mineralisation and chemical composition of that spring water which comes from thaw or thaw-rain recharge takes place, i.e., normally in the Spring seasons of cold-temperate climates. It is emphasised here that, hitherto, classical elaborations of this problem concerned mainly the areas where rapid thaws have not occurred. Moreover, the karst in those areas was not significantly affected by human impact, a factor which appreciably affects the interpretation of karst aquifer behaviour. The monitoring of karst springs using automatic recording devices, which is being carried out in many parts of the World, should provide much new information about this problem. The springs of the Cracow-Czestochowa Upland arc currently being monitored in this way. REFERENCES BAKALOWICZ M., 1979: Contribution de la geochemie des eaux ä la connaissance de I'aquifere karstique et de la karstification. These Doct. Univ. Pierre et Marie Curie - Paris 6, 269 p., Paris. BONACCI O., 1987: Karst Hydrology. With special reference to the Dinaric Karst. Springer, 184 p., Berhn-Heidelberg-New York. BONACCI O., 1993: Karst springs hydrographs as indicators of karst aquifers. Hydrological Sciences Journal, 38, 1,2, pp. 51-62, Amsterdam. FORD, D., WILLIAMS, P, 1989: Karst Geomorphology and Hydrology. Unwin Hyman, 601 p., London. GOSPODARIČ R., HABIČ P, (editors) 1976: Underground water tracing. Investigations in Slovenia 1972-1975. Inst. Karst Research, 312 p., Ljubljana. JAKUCS L., 1959: Nue Methoden der Höhlenforschung in Ungarn und ihre Ergebnisse. Die Höhle 10, 4, pp. 88-98, Wien. LAKEY B., KROTHE N.C., 1996: Stable isotopic variation of storm discharge from a perrenial spring, Indiana. Water Resources Research, vol. 32, no. 3, pp. 721-731. MANGIN A., 1975: Contribution ä I'etude des aquiferes karstiques ä partir de I'analyse des courbes de decrues et tarissement. Ann. Speleol., 29, 3: pp. 283-332; 29, 4: pp. 495-601; 30, 1: pp. 21-124, Moulis. PADILLA A., PULIDO-BOSCH A., 1995: Study of hydrographs of karstic aquifers by means of correlation and cross-spectral analysis. Journal of Hydrology 168, pp. 73-79, Amsterdam, REEDER Ph.R, DAY M.J., 1993: Seasonahty of chloride and nitrate contamination in the Southwestern Wisconsin karst. ApHed Karst Geology, Beck (ed.), Balkhema, pp. 53-61, Rotterdam. SHUSTER E.T., WHITE W.B., 1971: Seasonal fluctuations in the chemistry of limestone springs: a possible mean for haracterizing carbonate aquifers. J. HydroL, vol. 14, pp. 93-128, Amsterdam. TYC A., 1992: Action actuelle des processus karstiques dans le karst polygonal du Pan de Guajaibon (Sierra del Rosario, NW de Cuba). [W:] R. Maire, J.N. Salomon: Karst et evolutions chmatiques. Presses Universitaires de Bordeaux, pp. 135-144, Bordeaux. TYC A., 1994: Anthropogenic Impact on Karst Processes in the Olkusz-Zawiercie Karst Area (Silesian-Cracow Upland, South Poland). Wplyw antropopresji na procesy krasowe na przykladzie obszaru Olkusz-Zawier-cie (Wyzyna Slasko-Krakowska). PhD thesis, the University of A. Mickie-wicz, Poznan. TYC A., 1996: The nature of epikarst and its role in dispersed pollution of carbonate aquifers. International Conference on Karst-Fractured Aquifers - Vulnerability and Sustainability. Katowice-Ustroii, June 10-13, 1996, pp. 270-281, Katowice. WHITE W.B., 1969: Conceptual model for carbonate aquifers. Ground Water, 7, 3, pp. 15-21. WHITE W.B., STELLMACK J.A., 1968: Seasonal fluctuations in the chemistry of karst groundwater. Actes 4e Congr. Int. Speleol., Yougoslavie, Ljubljana, 3, pp. 261-267, Ljubljana. WORTHINGTON S.R.H., 1991: Karst Hydrogeolog>' of the Canadian Rocky Mountain. PhD Thesis, Department of Geography, McMaster University, Hamilton. ANALIZA KEMOGRAMOV IZVIROV V PRIMERU VPLIVA EFEKTA TALJENJA SNEGA IN IMPULZOV RAZPRŠENEGA ONESNAŽENJA (NA PRIMERU IZVIROV V KRAKOWSKO-CZESTOCHOWSKEM VIŠAVJU) Povzetek Napajanje, tok in vskladiščenje podzemne vode so lahko v karbonatnih vodonosnikih ekstremno spremenljivi. Skladno z različnimi stopnjami zakraselo-sti je lahko napajanje v razponu od koncentriranega do razpršenega. Članek predstavlja rezultate proučevanja onesnaženja glede na različne načine napajanja v karbonatnih vodonosnikih. V ta namen je bil izpeljan monitoring temperature vode in koncentracij glavnih ionov v nekaterih izvirih višavja Krakow-Czestochowa v južni Poljski. Analize kemogramov kraških izvirov, ki jih običajno srečujemo v mednarodni literaturi, so bile pregledane in ponovno ovrednotene. Tradicionalni modeli napajanja zaradi taljenja snega na kraških območjih zmernega podnebja in človeškega vpliva v zmerno-hladni klimi so bih ocenjeni in predlagane so bile nekatere spremembe. V predlaganem novem modelu se efekt taljenja snega in impulzi onesnaženja odražajo na kemizmu voda - značilna povečanja koncentracij Ca^+ , SO/' in NO,;, celokupne trdote in specifične električne prevodnosti vode v tleh, epikraške vode in vode na izvirih - pri maksimalnih vrednostih napajanja zaradi taljenja snega. To spremlja še relativno nizka vsebnost bikar-bonatov. Koncentracija kalcija in celokupna trdota v staljeni vodi sta bolj povezani s sulfati in nitrati iz onesnaženja kot pa z efektom raztapljanja v sistemu CaCO, - H,O - CO,.