Karst and Caves in Carbonate Rocks, Salt and Gypsum – poster 2013 ICS Proceedings

FIELD MEASUREMENTS OF GYPSUM DENUDATION RATE IN KULOGORSKAYA CAVE SYSTEM

Nikolay Franz1, Sergey Sorokin2,3, Alexandra Alexeeva4, Irina Inshina4, Olga Novysh3, Anton Kazak3 1Arkhangels Speleological Association “Labirinth”, Arkhangels, , [email protected] 2Tver State University, Zhelyabova 33, Tver, Russia, [email protected] 3Arkhangels Speleological Association “Labirinth”, Arkhangels, Russia 4Tver State University, Tver, Russia

Results of measurements of karst denudation rate in natural conditions in Kulogorskaya cave system are reported. Taking into consideration relatively high rate of denudation in gypsum caves we have used simplified MEM measurement method. Observation of hydrological situation suggests that most of the dissolution in caves of the area occurs during a few days of spring flood. Here we report results of measurement at 62 points located at 6 sites for the period from 2003 to 2012. Data analysis gives us average retreat rate of 0.0448 mm/day, while point is underwater. Confidence interval for this value is 0.0379 to 0.0516 mm/day.

1. Introduction annual rainfall 560 mm. The coldest month – January (average long-term temperature -13.0 °C). The warmest Kulogorskaya cave system includes 7 caves, three of which month – July (average long-term temperature +15.4 °C). constitute the largest cave system in Russia in a gypsum, with a total length of surveyed passages more than 17.5 km. Karstifiable rocks of massif are represented by gypsum and dolomite of Permian age. Primary karst rock is gypsum, Cave system is studied by speleologists from Arkhangels, St. forming homogeneous strata ranging from 0.2 to 7 m, Petersburg, Moscow, Tver and other cities under the guidance interleaved with dolomites. Average total width of of Arkhangels Speleological Association “Labirinth”. sedimentary cover above caves is 20 m. This report presents the preliminary results of the Caverns in the massif are well represented in all measurements of dissolution of gypsum in the natural hydrodynamic zones: vertical transfer zone, zone of conditions, which are held from 2003. seasonal fluctuations of water table levels and phreatic zone. Human-passable channels are located only at the zone of water table. In karst dissolution infiltrational, condensational and floodwaters are involved. Local erosion bases for caves in Kulogorskaya cave system is River and Pinega-Kuloy channel.

3. Conditions of karst denudation The most important role in the process of dissolution of gypsum walls in the water-table zone play floodwater penetrating the karst massif through numerous small ponors laid along the Kulogorsky ledge in the area of contact between the karst rocks and the floodplain alluvium. The penetration of floodwater into massif mainly occurs in spring (April–May) when Pinega waters rise by 4–5 meters above the winter low water (in the caves in this period water rises up to 2–3 meters, completely flooding almost all horizontal passages). During this period water flows in Figure 1. Geographical position of Kulogorsky Cave Region. horizontal passages in water-table zone, which are typically dry in other seasons. Absorbed water temperature typically ranges around 1 °C. 2. Geography and geology Autumn floods in the river and caves are less affluent: the Geographically Kulogorskaya cave system is located in the rise of water in the river – up to 2–3 m, in caves – up to 1 north of the East European Plain, in the watershed of Pinega m. The inflow of water into the massif in this period is only (tributary of ) and Kuloi (Fig. 1). through channels of phreatic zone (flooding of caves is caused by raising of the water-table level). Reverse runoff The climate here is quite severe, with low air temperatures after autumn flood occurs on the same path in the opposite and high humidity, rainfall exceeds evaporation. The multi- direction (the horizontal filtering through the thickness of year average annual temperature is +0.2 °C. The average the alluvium in the direction of erosion basis).

185 Karst and Caves in Carbonate Rocks, Salt and Gypsum – poster 2013 ICS Proceedings

Figure 2. Measurement method.

The main feature of the hydrogeology of Kulogorsky massif Rate of Karst Denudation in Russian) were fitted by us only is absence of significant inverse flow of flood waters after in the passages of this zone. passing through the peak of the spring flooding. In other words, floodwaters which penetrate the caves from the floodplain, never pour back by the same horizontal 4. Methods channels. Our current understanding is that at mean water period the water from caves is unloaded by slow horizontal To study the rate of dissolution of gypsum in the natural filtering through the thickness of alluvium in the direction conditions since 2003, we have equipped 15 measurement of the local base level of erosion. sites in five caves. Almost all sites were equipped on the surfaces of the ceilings and walls of passages laid in the To determine the saturation of cave waters with gypsum, massif of white-and fine-grained gypsum with rare clusters we measure its electrical conductivity, which is directly of large gypsum crystals. proportional to the salinity of water – both according to the data of other researchers (Krawczyk and Ford 2006), and Since the expected rate of dissolution of gypsum is higher confirmed by our observations. than that of limestone, we used a simplified, compared to the classic version of the MEM method (Smith, 1977). In 2004 we conducted experiment in K-2 cave on the dissolution of gypsum powder in the water from Each SKD site includes from 1 to 4 rows of control points. underground lake (initial conductivity 2.1 mS) at natural Rows of control points may be oriented horizontally or cave temperature +2.2 °C, which gives us maximal vertically. Each row, in turn, includes from 4 to 7 observation conductivity of 2.5 mS. points located with an interval of 3 cm along a straight line between the two fixed reference points. But these values are not reachable in the natural underground water exchange, which is convincingly shown by the results of long-term observations of hydrochemical karst springs in the river Pinega: conductivity of discharged water is never more than 2.2 mS (88% of maximum mineralization). Since caves are inaccessible to humans during flood, we have used the tanks with automatic closure to determine the salinity of flood waters. These tanks have been installed by us in 2005 at a different distance from ponors. Analysis of the data showed, that conductivity of aggressive flood water, which penetrates the massif is 0.5–0.6 mS (20–24% of maximum). After passing for many hundreds of meters underground transit, these waters are saturated to 1.8–1.9 mS (72–76%). Then, in a few months of slow water level drop, conductivity of underground water reaches 2.2–2.3 mS (88–92%). Thus, during the transit time (in a few days at the peak of the spring flood) in conditions of heavy water exchange, floodwaters are saturated with the products of dissolution of sulphate rocks up to 50% approximately. Over the following months, during slow level decline, water gather additionally about 15% to reach the maximum observed mineralization of 2.2–2.3 mS, or 88–92% of the maximum possible in these climatic conditions. Since, as mentioned above, main floodwater flow occurs in the zone of water-table, we can assume that maximal denudation of cave walls happens in this zone. That is why Figure 3. Location of measurement sites. all karst denudation rate measurement sites (“SKD” for

186 Karst and Caves in Carbonate Rocks, Salt and Gypsum – poster 2013 ICS Proceedings

Figure 4. Hydrograph of spring flood at Pinega river. Data from the state gauging station in Kulogory village in the immediate vicinity of the cave system (Federal State Unitary Enterprise “Centre of Russian water works inventory and state water cadaster”, 2012).

Reference points are stainless steel screws fixed to the wall exceptionally high rate where identified. They are shown with plastic plugs to a depth of 50–60 mm. During the with “high rate” comment in tables 1 and 4, and they are measurement special aluminum rod of square section with not used for any calculations. test holes every 3 cm is imposed on the heads of the screws Using data from the remaining points we calculate mean (Fig. 2). retreat rate for every measurement row and site. Vernier caliper probe is inserted to each hole and distance Average retreat rate for all observation points is 0.0448 between wall surface and outer surface of the rod is measured. mm/day. The confidence interval for this value from 0.0379 On every SKD site we made a special “passport” with the to 0.0516 mm/day. tables of the measured values on each set of control points. The column “conditions” of Tables 1 and 2 provides an We also made an exact measurement of height above sea assessment of denudation conditions for each of the sites. level of each control point row. In the sites marked “Corrosion” there is no signs of water flow. “Erosion” – conditions are suitable for water flow and there are signs of such flow. 5. Results and discussion It is interesting to note that the maximum retreat rate on This report presents the results of observations on six oldest surfaces of walls and the roof was measured at the SKD-8 SKD sites equipped in 2003–2004 and affected by nine site, equipped in one of the passages of Vodnaya (Water) spring floods (2004–2012). Locations of sites are shown in Cave, where the nature of water deposits suggests flood- Fig. 3. Totally 62 points were analyzed. flow velocities above 1 m/sec. and here, respectively, there Obviously, in the existing hydrogeological environment is a very intense water exchange. hypsometric position of the observation point plays a From these data it is clear that the denudation rate of cave significant role. Thus, the higher (in absolute elevation) is channels depends strongly on the local conditions, which is the reference point, the less time it is exposed to dissolving characteristic for gypsum karst (Klimchouk et al. 1996). On effects of flood waters, and vice versa. This is illustrated by the example of SKD-2 site one can see that the dissolution a spring flood hydrograph of Pinega River in 2012 (Fig. 4), rate can vary significantly even at points located in close which is fairly typical. proximity to each other. However, the obtained confidence Therefore, in the analysis of the data, we took into account limits give hope that average rate really characterizes the the absolute height of each of the control points. Taking into process of denudation within all massif. consideration available data on the dynamics of the spring Another trend that can be observed in collected data, is floods of the river Pinega for 2004–2012 years, it was possible higher rate of dissolution of horizontal surfaces (ceilings) to determine the estimated time (in days) of the interaction than of vertical (walls). However, the existing data does not between gypsum surface and aggressive floodwater for each reliably establish that fact. We hope to clarify this point, control point for the entire observation period. after completing the measurements on the other sites. The resulting values of the total retreat (in mm) at each test The resulting average speeds of karst denudation in real point and the data on the duration of flooding of each of them hydrogeological conditions of Kulogorskaya cave system allowed us to determine the retreat rate (Tables 1 and 2). can later be used to build a variety of mathematical models We have used Smirnov-Grabbs criterion on the values of and reconstructions of paleogeography environments and the retreat rate to identify outliers. Four points with stages of development of the caves of massif.

187 Karst and Caves in Carbonate Rocks, Salt and Gypsum – poster 2013 ICS Proceedings

We continue to equip new measurement sites and are Klimchuk A, Cicchi F, Calaforra JM, Aksem S, Finocchiaro F, working on improvement of our techniques. Long-term Forti P, 1996. Dissolution of Gypsum from Field Observations. monitoring of natural karst denudation processes will Gypsum Karst of the World. International Journal of undoubtedly allow us to deepen our knowledge of the laws Speleology, Theme issue, 1996, Vol 25, Issue 0, 37–48. of development of not only Kulogorskaya cave system, but Krawczyk WE, Ford DC, 2007. Correlating specific conductivity also of the entire northern sulfate karst. with total hardness in gypsum karst waters. Earth Surf. Process. Landforms, 32: 612–620. doi: 10.1002/esp.1409. Smith DI, 1977. The Micro Erosion Meter: Its application to the References Weathering of Rock Surfaces. Conservation of Rock Art. Proceedings of the International Workshop on the Conservation Federal State Unitary Enterprise “Centre of Russian water works of Rock Art, Perth, September 1977 (1978), 44–53. inventory and state water cadaster”. http://www.waterinfo.ru/ (in Russian).

Table 1. Measurement data and statistics, part 1.

Depth Depth Total Retreat Avg. rate Avg. rate Height, Days Pecularities in 2004 in 2012 retreat, rate, for row, for location,

Site

Row m a.s.l. underwater

Point

Sufrace mm mm mm mm/day mm/day mm/day

Conditions 1 16.69 25 22.0 23.7 1.7 0.0680 2 16.66 25 27.7 29.4 1.7 0.0680 A 4 16.60 25 35.0 36.4 1.4 0.0560 0.0690 Wall 5 high rate 16.57 25 33.6 36.8 3.2 0.1280 6 16.54 25 29.2 31.3 2.1 0.0840 1 16.36 36 26.5 28.2 1.7 0.0472 2 16.33 36 31.2 32.1 0.9 0.0250 B 0.0299 Wall 3 16.30 36 32.6 33.1 0.5 0.0139 6 crystall 16.21 36 29.1 30.3 1.2 0.0333 1 16.02 46 35.0 38.4 3.4 0.0739 0.0475 2 15.99 46 34.2 36.0 1.8 0.0391

SKD-2

Erosion C 3 15.96 46 28.4 31.2 2.8 0.0609 0.0591 Wall 4 15.93 46 27.8 30.5 2.7 0.0587 5 crystall 15.90 46 25.0 27.9 2.9 0.0630 1 15.35 79 25.0 27.3 2.3 0.0291 2 15.32 79 23.0 25.1 2.1 0.0266 3 15.29 79 25.6 28.5 2.9 0.0367 D 0.0352

Wall 4 15.26 79 23.3 26.6 3.3 0.0418 5 15.23 79 26.2 29.3 3.1 0.0392 6 15.20 79 28.0 31.0 3.0 0.0380 1 15.57 51 18.1 19.9 1.8 0.0353 2 15.57 51 16.2 16.8 0.6 0.0118 A 0.0196

Wall 3 15.57 51 16.0 16.9 0.9 0.0176 4 15.57 51 13.8 14.5 0.7 0.0137 1 16.27 33 18.0 19.2 1.2 0.0364 0.0330

SKD-3 2 16.27 33 17.3 18.3 1.0 0.0303

Corrosion B 3 16.27 33 17.2 18.3 1.1 0.0333 0.0436 Wall 4 16.27 33 18.0 18.8 0.8 0.0242 5 16.27 33 17.7 20.8 3.1 0.0939 1 15.18 80 16.4 17.3 0.9 0.0113 2 15.18 80 12.0 13.1 1.1 0.0138 A 3 15.18 80 17.8 19.5 1.7 0.0213 0.0148 0.0148

Wall

SKD-4 Erosion 4 15.18 80 18.5 19.7 1.2 0.0150 5 15.18 80 20.0 21.0 1.0 0.0125 1 15.34 72 11.2 13.0 1.8 0.0250 2 15.34 72 14.8 16.8 2.0 0.0278 A 0.0378

Wall 3 15.34 72 16.0 19.2 3.2 0.0444 5 15.34 72 25.4 29.3 3.9 0.0542 0.0372 1 15.70 50 14.4 16.5 2.1 0.0420

SKD-5 Erosion 2 15.70 50 18.5 20.1 1.6 0.0320 B 0.0365

Wall 3 15.70 50 16.3 18.6 2.3 0.0460 4 15.70 50 20.1 21.4 1.3 0.0260 1 15.94 44 21.5 22.5 1.0 0.0227 2 15.94 44 21.0 23.0 2.0 0.0455 3 15.94 44 20.0 21.4 1.4 0.0318 A 0.0352 0.0352

Wall 4 15.94 44 20.8 22.3 1.5 0.0341

SKD-6 Erosion 5 15.94 44 21.5 23.0 1.5 0.0341 6 15.94 44 20.8 22.7 1.9 0.0432

188 Karst and Caves in Carbonate Rocks, Salt and Gypsum – poster 2013 ICS Proceedings

Table 2. Measurement data and statistics, part 2.

Depth Depth Total Retreat Avg. rate Avg. rate Height, Days Pecularities in 2004 in 2012 retreat, rate, for row, for location,

Site

Row m a.s.l. underwater

Point

Sufrace mm mm mm mm/day mm/day mm/day

Conditions 1 high rate 15.90 45 14.2 21.9 7.7 0.1711 2 15.90 45 15.2 20.3 5.1 0.1133 A 3 high rate 15.90 45 15.2 23.5 8.3 0.1844 0.0844

Celling 4 15.90 45 14.2 18.6 4.4 0.0978 5 15.90 45 15.0 16.9 1.9 0.0422 1 15.60 52 14.3 18.2 3.9 0.0750 2 15.60 52 15.5 19.8 4.3 0.0827 0.074 B 3 15.60 52 21.0 24.1 3.1 0.0596 0.0723

Wall SKD-8 4 15.60 52 22.6 26.3 3.7 0.0712 5 15.60 52 19.6 23.4 3.8 0.0731 1 15.30 73 20.1 25.5 5.4 0.0740

Erosion, high-speed currents 2 15.30 73 16.1 21.4 5.3 0.0726 C 0.0667

Wall 3 high rate 15.30 73 13.2 23.2 10.0 0.1370 4 15.30 73 16.0 19.9 3.9 0.0534

189