Next Article in Journal
Body Composition Analysis in Perimenopausal Women Considering the Influence of Vitamin D, Menstruation, Sociodemographic Factors, and Stimulants Used
Previous Article in Journal
Friendly Residential Environments and Subjective Well-Being in Older People with and without Help Needs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogeochemical Response of Cave Drips to Precipitation during Rainfall in a Karst Desertification Region: A Case Study of Shijiangjun Cave, South China

1
School of Karst Science, Guizhou Normal University, Guiyang 550001, China
2
State Engineering Technology Institute for Karst Desertification Control, Guiyang 550001, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Environ. Res. Public Health 2022, 19(23), 15830; https://doi.org/10.3390/ijerph192315830
Submission received: 18 July 2022 / Revised: 3 October 2022 / Accepted: 24 November 2022 / Published: 28 November 2022
(This article belongs to the Section Environmental Science and Engineering)

Abstract

:
Exploring the hydrogeochemistry of cave drip water and its response to precipitation events in karst rocky desertification regions is of great significance to the paleoenvironment reconstruction of the karst desertification process using speleothem. We selected three perennial drip sites in the Shijiangjun Cave, located in Guizhou Province, Southwest China, and carried out high-frequency monitoring and sampling during two rainfalls from 22 to 25 May 2016. The major hydrogeochemical parameters of drip water and their relationships with karst desertification were analyzed. The results show that the hydrogeochemistry of the drip water in the Shijiangjun Cave, characterized by HCO3-Ca·Mg, was dominated by the dissolution of calcareous dolomite. The three drip sites were classified into the delayed response type (W1) and the rapid response type (W2 and W3) based on the response speed of the drip water indicators to precipitation, which were highly influenced by the piston effect and precipitation dilution, respectively. Furthermore, the response sensitivity of the drip water indicators to precipitation was constrained by the desertification degree in the rainy season, specifically, the faster response appeared in the higher desertification degree area. It is essential to select appropriate drip sites and establish an applicable indicator system for the evolutional history reconstruction of karst desertification using speleothems.

1. Introduction

There are large outcropping carbonate rocks and complex three-dimensional geological structures of duality (aboveground and underground) in karst regions, of which the structure soil in the context is distributed in an island-liked pattern [1], resulting in a disproportionate spatial pattern of soil and water resources. The karst eco-environment ultimately has become one of the most typical vulnerable environments all around the world. The vulnerability of the karst eco-environment, as demonstrated by karst rocky desertification, has restricted local social and economic development in the karst region of South China [2,3]. In recent years, many innovative methods have been utilized to monitor environmental changes before and after karst rocky desertification control, including remotely sensed imagery interpretation, the dynamic variations of vegetation and soil physicochemical properties, and so on [3,4,5,6,7,8]. Although the above surface monitoring methods provide promising approaches to the establishment of a karst rocky desertification monitoring system, however, the research on the recent conditions and variations of karst rocky desertification by way of underground medium or indicators has been rarely reported until now. Therefore, building an integrated environmental monitoring system based on the characteristics of dual structure in karst regions is extremely important for karst rocky desertification control evaluation and would be a breakthrough in improving and developing the current monitoring system.
Speleothems in karst caves (stalagmites) are formed due to the interactions of a multiple-surface environment, such as precipitation, an overlying ecosystem, karst seepage water, and other systems [9,10], and they record abundant information on the paleoclimate and paleoenvironment [11,12,13]. Coplen et al. previously analyzed the stalagmite δ18O and δ13C in the Devils Hole in Nevada and reconstructed the paleoclimate and paleovegetation changes in the past 500,000 years in this region [14]. A consensus has been reached that the stalagmite δ13C obviously responds to the CO2 δ13C in the soil overlying the carbonate rock layer, and subsequent studies have focused on using the stalagmite δ13C to trace the changes in the paleoenvironment and the recent eco-environment [15,16,17]. Liu et al. compared the stalagmite δ13C in Guizhou Province, Southwest China, and found that the stalagmite δ13C in the Shijiangjun Cave was close to that of the surrounding bedrock and far higher than that in the Zhijin Cave [18], the evolutional history of karst rocky desertification in the past 800 years was further discussed in terms of the stalagmites δ13C and δ18O in these two caves, which proved the influence of karst rocky desertification degree on stalagmites δ13C. Cave drip water is an important medium of speleothems formation, which exhibits different geochemical characteristics in various karst rocky desertification environments. For instance, drip water that developed in severe karst rocky desertification conditions was commonly characterized by lower Ca2+ and Mg2+ concentrations than that in a mild desertification environment, as well as variable amplitudes in their ionic concentrations [19]. Although previous studies have presented valuable results, there are still obvious differences in the seepage paths and hydrodynamic properties in the formation process of drip water, hydrochemical compositions of drip water, as it described via the indicators, reflects the constraint by various kinds of hydrogeochemical processes [20], which will certainly influence the response degree of drip water indicators to the overlying environment of karst rocky desertification, and lead to various interpretations in revealing environmental changes by drip water. Therefore, it is critical to systematically sort out and deepen the existing work.
Generally, the signal transmission of groundwater from the earth’s surface to the cave is highly influenced by vegetation type, soil depth, bedrock thickness, cave pattern, drip water location, and hydrological and meteorological conditions [21,22]. However, due to the unique binary karst geological structure, the response of groundwater to atmospheric precipitation is relatively more sensitive than that in non-karst regions [3]. The investigation into hydrogeochemical characteristics of cave drip water under the condition of precipitation is not only conducive to the identifications of hydrological processes (recharge, storage, and migration) of seepage water in the upper vadose zone [23,24], but is also significant to the environment signal interpretation on the basis of speleothems for the reason that the formation process and influencing factors of cave speleothems could be better understood [25]. The drip sites of the Shihua Cave in north Beijing, China, were classified into two groups (rapid and delayed response type) according to the response characteristics of the drip water to precipitation [26]. The response of the drip water to precipitation in the Grotta di Ernesto Cave in northeast Italy was characterized by three patterns: rapid seasonal drip water, slow seasonal drip water and slow drip water [27]. Research on the Liangfeng Cave in Southwest China revealed that information on extreme precipitation events could be recorded by the hydrochemistry and stable isotope variations of drip water, and the intensity of precipitation played a key role in constraining the Ca2+ and HCO3 concentrations [28]. Owing to the complexity of hydrological processes from the soil layer to the surface vadose zone, it is difficult to fully understand the historical information of each drip site. However, the frequency and contemporary hydrogeochemistry of drip water under precipitation conditions provide some clues to clarify the hydrological processes of karst groundwater, such as the water-rock interaction degree, the water retention time, the piston effect, and so on [29]. Therefore, the Shijiangjun Cave under the conditions of moderate karst rocky desertification was selected as the study area. Three drip sites were monitored and sampled at high frequency during the rainy season, aiming to interpret the particular response characteristics of drip water hydrogeochemistry to precipitation in typical karst rocky desertification areas and to explore which indicators of drip water under different hydrologic conditions could potentially reflect the surface environment. This research is expected to provide theoretical references to the paleoenvironment reconstruction of karst rocky desertification using speleothem and karst rocky desertification control at the current stage.

2. Materials and Methods

2.1. Study Area

The Shijiangjun Cave (26°17′03′′ N, 106°03′54′′ E) is located in the central part of the Guizhou Karst Plateau, Southwest China (Figure 1a), which is the watershed between the Yangtze River Basin and the Pearl River Basin. It developed in the thick Middle-Upper Triassic calcareous dolomite strata (T2p-l) with a small amount of evaporative minerals (mainly halite), sandstones and shales in the interlayers [30]. Tectonically, it is located at the southwestern edge of the Sichuan-Guizhou Meridional Structural System, and it developed in the west wing of the Maolipo Anticline, the formation dips are generally from 10° to 25° [30]. The area around the cave system is characterized by a typical subtropical monsoon climate, with an average annual precipitation of 1290~1380 mm and annual temperature of 14~16 °C, and the rainy season is usually from May to October. The geomorphic assembled forms are mainly fengcong (peak cluster)-valley and fenglin (peak forest)-depression, and poljes are quite common for this karst plateau. The regional base level of erosion is close to the bottom of the depression, and the groundwater type is mainly carbonate fissure-cavity water. Shrub grass is the main vegetation type in the surrounding region, with most of it normally distributed in the low-lying areas of mountains. Moderate and intense karst rocky desertification developed in the peak body, and the karst rocky desertification in the overlying environment of the Shijiangjun Cave was graded as moderate due to the steep peak body and the high rock exposure rate of more than 50%.
The Shijiangjun Cave, which is more than 500 m in length was formed in the Shanpen Period (a physiographic stage from the Early Miocene to the Pleistocene), during which period the earth’s crust was stable and most cave channels were horizontally developed [30]. The only cave entrance appears at the bottom of the peak body on the edge of the depression (Figure 1d), and the cave halls are distributed in a beaded-like form. The cave width is from 3.1 to 48 m and the cave height is from 1.8 to 20.3 m. The rock thickness over the cave ceiling is from 79.7 to 118.2 m. Many vertical fissures have developed in the cave roof, including tufted stalactites, sodastraws, and stone curtains due to the chemical precipitation effect of drip water. There are abundant drip sites and well-developed speleothems, such as stalagmites, sodastraws, stalactites, stalagnates, stone curtains, stone flags, and so on. The drip water forms several discontinuous water flows in the cave passage and flows underground. The average annual CO2 concentration in the cave is 409 ppm, the annual temperature is 15.16 °C, and the annual humidity is 93%. There are small amounts of shrubs and grasses growing in the peak body, and only a few small trees growing in deep soil layers and large fissures, however, there are many bare rocks and thinner soil layers overlying the cave. Most soils were washed into the cave space through the cavernous fissures and cracks in the past, forming a 2~4 m sedimentary layer at the bottom of the cave tunnel, and most huge stalagmites developed on the humid sediments (Figure 1e).
Three perennial drip sites, numbered sequentially as W1, W2, and W3 and located in the middle section of the cave tunnel (Figure 1b,c), were selected as the monitoring objects in this study, the drop heights of which are 7.1, 20.1, and 13.1 m, respectively. The cave roofs at the three drip sites consist of calcium plate and new carbonate deposits, and several stalagmites have developed below each drip site.

2.2. Field Monitoring and Sampling

Monitoring and sampling were conducted at the three drip sites every 4 h from 23:00, 22 May to 20:00, 25 May 2016 after a long dry period and before the following rainstorm, and only 54 samples were obtained. Real-time temperature and rainfall data were collected using the WM5200 weather station (Vaisala Company, Helsinki, Finland) in the open space near the cave. The drip speed of the drip water was measured by a stopwatch in each sampling group, and the average drip volume, defined as the drip rate in this study (unit: mL·min−1), was calculated by measuring the water quantity with a 100 mL glass cylinder in 5 min. The water temperature, pH, and electrical conductivity (EC) of the water samples were determined in situ by a portable water quality analyzer (HQ40d, HACH Company, Loveland, CO, USA) with the measurement precisions of 0.1 °C, 0.01 and 1 µS·cm−1, respectively. The pH probe was calibrated with buffer solutions with pH value of 4.005, 7.000, and 10.012 before the measurement. The HCO3 concentration was determined by a basicity meter (Merck Company, Darmstadt, Germany) with the test precision of 0.1 mmol·L−1. Water samples were collected in a 50 mL fully rinsed polyethylene bottle, meanwhile, 1:1 HNO3 was added to acidify the sample to pH < 2 in order to ensure the activity of the cations. Then, the polyethylene bottle was sealed with parafilm (Bemis Company, Neenah, WI, USA) for the subsequent determination of the cations. The remaining water samples were put into a 250 mL sampling bottle. Then, all samples were placed in a portable refrigerator in the field and refrigerated at 4 °C. Finally, the water samples were brought back to the laboratory for cryopreservation, and the cations and anions were tested within 48 h. All the above-mentioned sampling bottles were soaked in dilute acid for 24 h before sampling, and then washed with deionized water more than three times.

2.3. Laboratory Measurements

The samples of drip water were brought to the Geochemical and Isotope Laboratory of Southwest University, China, to detect the cation concentrations of K+, Na+, Ca2+, Mg2+, and Sr2+ by using an Optima-2100DV inductively coupled plasma emission spectrometer (ICP-OES, Perkin-Elmer Company, Waltham, MA, USA). The detection precision was 0.001 mg·L−1 with a standard deviation of less than 2%. NO3 and SO42− were analyzed by a UV2450 ultraviolet spectrophotometer produced by the Shimadzu Company (Kyoto, Japan), and the Cl was tested according to the method of AgNO3 titration. All samples were filtered by glass fiber membranes (Whatman Company, Maidstone, UK, 0.45 μm) before the ion determination. In addition, the hydrogeochemical simulation program of WATSPEC 15 (Wigley, Adelaide, Australia) was used to calculate the calcite saturation index (SIC) and dolomite saturation index (SID) of the water samples based on the data sets of the water temperature, pH, and main ions. The data statistics and analyses were completed in IBM SPSS Statistics 20, and cartographic drawing was performed using Origin 2018 software.

3. Results

3.1. Physico-Chemical Characterization of Cave Drip Water

The comparative analysis of the drip water hydrological and hydrochemical parameters during the precipitation process (Table 1) at the three drip sites in the Shijiangjun Cave shows that the drip rate in unit time at each drip site displays a relationship of “W2 < W1 < W3” (1.79 ± 0.25, 2.01 ± 0.12, 2.22 ± 0.31 mL/min, respectively). The drip rate and the variable amplitude had the largest values at W3, and the hydrological conditions of W2 were comparatively more stable during the monitoring period. All the drip water samples presented obvious weak alkalinity, while both W2 and W3 exhibited slightly stronger alkalinity than that of W1 (8.62 ± 0.14). The EC shows a relationship of “W1 > W2 > W3”, and the range of variation in W1 was much larger than that of W2 and W3 (364.56 ± 21.66, 272.61 ± 1.54, 231.39 ± 0.78 μS/cm, respectively). Most cation and anion concentrations had a consistent variation with the EC except for those of K+, Sr2+ and HCO3. The variable amplitudes of most ion concentrations at W1 were higher than those of W2 and W3, showing a similar spatial pattern of EC. Generally, Ca2+ and Mg2+ dominated the cations in the samples collected from the three drip sites, while HCO3 dominated the anions, showing the typical hydrochemical characteristics of karst regions.
The dataset collected at each drip site was processed in the WATSPEC 15 program (Table 1). The results show that both the SIC and SID for the drip water were greater than zero, showing a supersaturated status of calcite and dolomite in the drip water. For the three drip sites, the mean values of SID (1.68 ± 0.23, 1.55 ± 0.12 and 1.49 ± 0.11, respectively) were significantly larger than those of SIC (0.89 ± 0.11, 0.72 ± 0.06 and 0.63 ± 0.66, respectively), which indicates that the calcite precipitation effect in drip water was weaker than that of dolomite and the dolomite had a stronger precipitation tendency than that of calcite. This also reflects that the drip water hydrochemistry mainly inherited the characteristics of dolomite [31], which is consistent with the geological background of calcareous dolomite in the Shijiangjun Cave. Furthermore, all the carbonate saturation indices, including the SIC and SID at the three drip sites follow a rank order of “W1 > W2 > W3”, indicating that there was a relatively better dissolution ability of the drip water of W2 and W3 than that of W1; in other words, W1 was more prone to calcite and dolomite precipitation than W2 and W3.

3.2. Proportions of Cations and Anions

The drip water hydrogeochemical characteristics can also be seen in the following Piper diagram (Figure 2). The proportions of each cation to the total cations were significantly different at the various drip sites. At W1 and W2, the proportions of Ca2+ (57.85~65.21% and 52.22~53.9%) were obviously greater than those of Mg2+ (30.6~34% and 43.19~45.93%). On the contrary, a relatively higher proportion went to Mg2+ at W3 (49.7~52.99%). It can also be seen that, the proportions of K+ for W2 and W3 (0.43~0.73% and 0.56~0.93%) were both lower than those of Na+ (1.19~2.31% and 0.81~1.6%), while the average K+ proportion of W1 (2.44%) was more than twice that of Na+ (1.1%), and the proportion of Sr2+ to total cations was the lowest at all drip sites (0.03%, 0.02%, and 0.02%). In addition, the proportions of anion concentration to the total anions at all three drip sites are basically ranked as “HCO3 > SO42− > Cl > NO3”. Specifically, the anion proportions were 87.33%, 8.96%, 2.37%, and 1.34% for W1, 88.92%, 8.13%, 1.54%, and 1.42% for W2, 94.21%, 3.21%, 1.58%, and 0.99% for W3.
The hydrochemical types of drip water could be classified into the type of HCO3-Ca·Mg based on the analysis results of the Piper diagram, which also shows the particular hydrochemical composition in the dolomite region. The concentrations of Na+ and K+ were relatively lower, which indicates that the silicate weathering contributed limitedly to the cations. The above fact suggests that carbonate weathering plays a dominant role in the hydrochemical composition of drip water, especially the dolomite dissolution. The small amount of Na+ and Cl in the drip water is derived from the evaporite minerals dissolution in the dolomite strata [32]. Comparatively, the rainwater had a limited contribution to the substances in the drip water. According to the previously investigated data (Table 2), the max ion concentrations of the rainwater were relatively lower, for the cations (K+, Na+, Ca2+, Mg2+, and Sr2+), the max values were 0.37, 0.2, 2.8, 0.63 and 0.003 mg·L−1, respectively; and for the anions (HCO3, SO42−, NO3, and Cl), the max values were 12.2, 1.48, 1.01, and 0.38 mg·L−1, respectively. Rainwater mainly affects drip water formation by changing the seepage path and hydrodynamic conditions [33].

4. Discussion

4.1. Precipitation Response and Hydrological Process

According to the data statistics from a small meteorological station, there were two precipitation events in this monitoring period. The first precipitation was from 23:07, 22 May to 14:20, 23 May 2016, and the second was from 0:11, 24 May to 15:42, 24 May 2016. The rainfall amounts were 24.7 and 8.3 mm with the peak values appearing at 7:00, 23 May and 6:00, 24 May, respectively.
Based on the analysis of the relationship between the drip rate and the rainfall, it was found that the drip rate at W1 showed a downward trend during the first precipitation event, and a transient increase at 20:00 on 25 May (Figure 3). It began to rise 4 h after the rainfall peak in the second precipitation event, and a peak value appeared after about 10 h, then declined, and a valley value appeared until 8:00 on 25 May. The drip rate at W2 showed a declining tendency during the previous precipitation, gradually increased in the second precipitation, and then declined after maintaining a stable state for 16 h, and reached its peak around about 20:00 on 24 May. The drip rate at W3 also showed a downward trend after the initial precipitation and increased obviously at 20:00 on 23 May. The first peak value appeared after the second precipitation, and the second one appears at 20:00 on 25 May. Thus, it can be preliminarily concluded that the W1 had the largest hysteresis in responding to the precipitation, while the W2 and W3 were sensitive to respond to precipitation and showed better seepage path connectivity. Both W2 and W3 were classified into the drip site group of the rapid response type.
By comparing the variation trend of the EC and pH in different precipitation periods (Figure 3), we further confirmed the above understanding. The peak value of the EC for W1 appeared 8 h after the beginning of the first precipitation, and the EC started to gradually increase after 2 days of a stable period. It is inferred that there is an overlying thick cave roof above the site of W1. As a result of the influence of the “piston effect”, the “old water”, which was stored in the epikarst before the precipitation [34,35], pioneered to enter the seepage path and formed drip water with a high value of EC. Subsequently, the “old water” flowed into the underground cavity, and the EC of drip water became to stabilize, and started responding to the precipitation in the exoteric environment until the increase in the EC in the later period. The EC of both W2 and W3 exhibited intense fluctuations, declining at least three times during the precipitation period, but the average EC of both sites was lower than that of W1. The variation in the pH values shows a rising tendency after the precipitation, which was relatively lagging behind that of the EC. According to the results of a previous study, the karst groundwater hydrogeochemistry in the process of precipitation is sensitive to the dilution of rainwater [36], and the EC decreased along with the increasing drip rate at W2 and W3, demonstrating that the drip waters at both W2 and W3 in the Shijiangjun Cave were also significantly affected by the dilution effect. In addition, this also indicates that the retention time of the seepage water in strata was relatively shorter, and the water rock interaction was insufficient, suggesting better fissure connectivity above the two drip sites. Therefore, the latter two drip sites responded rapidly to the precipitation events, further confirming the above judgment.
The annual and high-resolution monitoring of drip rate, EC, and other parameters of drip water helps to clearly explain the seasonal variations of drip water and aquifer structure of the karst cave [37]. Therefore, we compiled a data graph set on the physical and main ion parameters of the drip water from July 2014 to June 2015 in the Shijiangjun Cave (Figure 4). The data atlas includes the average monthly rainfall, the drip rate (in mL/min), EC, pH, calcium, magnesium, and bicarbonate of the monthly water samples of the three drip sites. The field sampling of monthly water sample was completed at the end of each month. It is showed that all the drip rates and EC of the three drip sites had positive responses to the precipitation in the monitoring of this hydrological year. Overall, the drip rates and EC of the three drips were at a high level in the rainy season, while relatively lower in the dry season (from November 2014 to April 2015). In comparison, there was a relatively stable trend of the seasonal variation of drip rates and EC for W1, while those of W2 and W3 fluctuated greatly, especially for W3. According to the relationship between the average drip rate and the Coefficient of Variation of the drip rate, the three drip sites can be clearly classified into two groups during both the rainfall monitoring (Figure 5a) and the hydrological year (Figure 5b). The first type includes W2 and W3 with large variable amplitudes of drip rate and obvious response to rainfall, and the other one only includes W1, which had relatively stable variation of the drip rate. Furthermore, the drip rates and EC of W2 and W3 obviously decreased with the rainfall in August 2014, while there was not significant change for W1. The drip rates and EC of the three drip sites showed obvious upward trend with the increasing rainfall after a long dry period in May 2015. Specifically, the variable amplitude of EC for W1 was much larger than those of W2 and W3, which again confirmed the above viewpoint that W1 was significantly affected by the piston effect. Although the EC of W2 and W3 also increased, the overall level (about 250 μS/cm) was much lower than that of W1 (350 μS/cm), which was related to the abovementioned rainwater dilution and insufficient water rock interaction, indicating that the karst tunnels in aquifers above W2 and W3 had stronger connectivity than that of W1.

4.2. Response of Drip Water Geochemistry to Precipitation

The results of the tracer experiment using NaCl, related to the hydrodynamic process of cave drip water responding to precipitation [38], suggest that the drip water in the Shijiangjun Cave responded to the atmospheric precipitation in a timely manner and indicate that the rainwater could quickly pass through the strata in the cave roof, so that the retention time of underground water was shortened and the water rock interaction intensity was weakened owing to the vertical fissures in the cave roof. Then, a substantial number of the main chemical compositions of the drip water were directly derived from the eluviation of the overlying soil layer. As reported in previous research, Ca2+ and Mg2+ in the drip water of the Shijiangjun Cave were restricted by the development degree of the overlying soil [20]. Specifically, the Ca2+ and Mg2+ concentrations in the drip water increased with the thickness of the soil layer; in other words, the thicker soil layer contributed more Ca2+ and Mg2+ to the drip water. As the dominant cations of the drip water, the Ca2+ and Mg2+ variations were essentially in accordance with that of the EC from the monitoring data at the three drip sites during the rainfall process in this study. The Ca2+ and Mg2+ concentrations at W1 gradually stabilized after a brief fluctuation, then presented a slight upward trend, indicating that the surrounding rocks were fully dissolved by the groundwater with the extension of the water retention time, so that the response of the drip water to the precipitation was relatively more delayed. The Ca2+ and Mg2+ concentrations at W2 and W3 fluctuated greatly, with a rapid increase after each precipitation (Figure 3), indicating that the drip water was mainly “new water” formed in the shallow soil layer after each precipitation event and that the Ca and Mg in the soil layer were quite vulnerable to be leached and lost by rainwater [39].
The seasonal variations of Ca2+ at the three drip sites were relatively stable in the monthly monitoring of a hydrological year, there were not obvious peak and valley values of calcium ions during the observation (Figure 4), which is different from those of the major ions at the perennial drip sites of other caves in South China. For example, the high-frequency study of drip water during rainfall at six drip sites in the Liangfeng Cave, Guilin, showed that Ca2+ had a large variation (40~120 mg/L) in the dry season from October, 2015 to January, 2016 [28], it was pointed out that the rainfall in the dry season was the major factor affecting the hydrochemistry of drip water. However, the Ca2+ concentrations of the drip water from the Shijiangjun Cave had remained at a relatively lower level throughout the hydrological year, which may be affected by the prior calcite precipitation. In contrast, the seasonal fluctuations of Mg2+ concentrations were significantly larger than those of Ca2+, and the variation trends were close to those of pH in rainy season, which indicates that the high temperature and groundwater full of CO2 promoted the dissolution of calcium dolomites in rainy season.
Generally, K+ and Na+ in natural surface water are mainly from atmospheric precipitation, silicate weathering, and evaporite dissolution. However, those in the cave drip water are largely contributed by the input of soil layer [40]. Accordingly, because the bedrock in the Shijiangjun Cave is mainly composed of calcareous dolomite and the surrounding environment is hardly influenced by human activities, the K+ in the drip water is mainly derived from the eluviation of the overlying soil, while the atmospheric precipitation contributes less to the K+ concentration (see Section 3.2). In addition to the dissolution of halite, soil leaching is also one of the sources of Na+. Since the K+ and Na+ concentrations at the three drip sites were relatively lower and influenced by the dilution effect of rainwater to various degrees, their responses of to the variations in rainfall are not obvious. Typically, clay minerals in the soil layer have strong adsorption properties to the element of Sr. It was pointed out that the Sr2+ in the drip water was relatively active, and its concentration variation was closely related to the precipitation outside in previous studies [41], which could also be found in the fluctuations of the Sr2+ concentration at W2 and W3 with the precipitation events and is consistent with the results of our study. Furthermore, the spatial pattern of the hydrogeochemical indicators shows that most of the cation concentrations at W1 were higher than those at W2 and W3 (Figure 6), indicating that the response sensitivity of the drip water indicators to precipitation was influenced by the differences in the material sources for the various drip waters.
As one of the main participants in karstification, bicarbonate occupies a large proportion of the total anions and TDS at the three drip sites. However, the dynamic changes in the HCO3 were not consistent with that of the EC during the precipitation, especially at the drip site of W1. The HCO3 of W1 showed a better steady behavior during the first precipitation, as suspected from the previous studies [34,38], and there was a thick soil layer above the drip site of W1. It started to decline gradually in the second precipitation during the first several hours, then increased until 4:00, 25 May and decreased again 12 h later. The rainwater, flushing through the soil layer in the karst region, carries a great amount of CO2 and flows into the karst fissures or karst tunnels, forming drip water with supersaturated carbonate. Actually, this resulted in a decreased trend of Ca2+ with the increasing Mg/Ca, indicating that the calcite saturation and precipitation had occurred before the seepage water drops (prior calcite precipitation) (Figure 7), which could also be verified by the carbonate suture structures distributed in the cave roof, especially for W2 and W3. Meanwhile, the fissures developing in the cave roof strengthened the vertical connectivity of the strata, resulting in carbonate-saturated water prone to CO2 degassing, which led to the effect of prior calcite precipitation effect. Furthermore, both the SIC and SID of W1 were generally large compared to those of W2 and W3 (Table 1), and in Figure 7, the data points of Ca2+ vs. Mg/Ca for W1 tended to be closer to the weak end of prior calcite saturation than those of W2 and W3, which can be summarized as follows: the W1 could be characterized by different hydrogeochemical processes compared to W2 and W3.
The geochemical behavior of chlorine, which is mainly from sea salt deposition, is relatively conservative in the natural process of the water cycle, therefore, Cl is usually used as a potential indicator and tracer to explore the response of drip water to atmospheric precipitation [42]. Although Cl is not the main component of anions in the drip water in the Shijiangjun Cave, it had a significant response to both of the two precipitation events due to the characteristics of a single source, demonstrating the dilution effect of rainwater on the hydrogeochemistry of drip water (Figure 3). Particularly, there were two peak values of Cl for W1 during the monitoring period, indicating the response of drip water to precipitation influenced by the piston effect, while the following lower values after each peak shows the impact of rainwater dilution. It is worth mentioning that the response degrees of different drip sites are subject to the thickness of the cave roof and multiplicate flow paths. The variation trends of SO42− and NO3 in drip water were not consistent throughout the precipitation process in this study. The SO42− concentration is characterized by an intense fluctuation with a slightly increasing trend and shows a better response to atmospheric precipitation. Previous studies have pointed out that the sulfur element in drip water is not directly derived from the atmosphere, but contributed by the sulfur compounds (the coal seam or coal soot) in the cave roof [43,44]. The overlying soil environment is another source of sulfate in cave drip water. Specifically, the sulfur fixed in the soil layer is oxidized and forms sulfate before the precipitation [45]. The leaching effect of the rain increases the SO42− concentration and maintains a high level until the preferential flow disappears. The changes in nitrate in drip water in the Shijiangjun Cave were obviously different from that of sulfate. The NO3 concentration of W1 presents an obvious downward trend with fewer fluctuations, while those of W2 and W3 reached a peak after the first precipitation, and gradually increased after the second one. Generally, the NO3 concentration followed the rank order of “W1 > W2 > W3” at the different drip sites from an overall perspective (Figure 6) due to the better soil development at W1. The nitrate in the drip water was mainly derived from the overlying soil nitrification. Therefore, the changes in the NO3 in the drip water mainly reflect the different leaching degrees of the soil under the condition of precipitation.

4.3. Influence of Karst Rocky Desertification on the Response of Drip Water Indicators to Precipitation

The classification standard list of the degree of karst rocky desertification involves multiple parameters, such as the vegetation coverage, the rate of outcrop rock, soil depth, slope, and so on [2]. As a result, the karst rocky desertification grade is actually used as a comprehensive evaluation of the karst rocky desertification degree at the regional scale. Furthermore, the karst region itself is characterized by discontinuous vegetation and soil layer, and the heterogeneity of the eco-environment displayed more obviously in moderate and intense karst rocky desertification regions. In particular, the strong aggregation effect of the bare rocks not only leads to an inhomogeneous intra-regional distribution of soil and vegetation, but also results in the differences in the soil physicochemical properties and the plant species diversity [46]. Accordingly, the numerous drip water developed in the same cave experienced different overlying environments in the formation processes, and the surface karst rocky desertification degree above each drip site is different. As mentioned above, the overlying bedrocks for W2 and W3 are relatively thinner, with developed fissures and well-connected groundwater channels; therefore, the drip waters carry more environmental information about karst rocky desertification driven by rainwater. For instance, the cations of K+ and Sr2+ in the drip water can better reflect the leaching degree of the overlying soil, as well as the mineral nutrients in the soil under the background of karst rocky desertification to some extent. However, for W1, the elevation of the drip site is relatively lower than those of W2 and W3, the overlying soil layer and strata are thicker, with a milder karst rocky desertification degree, and the flow paths are of poor connectivity. Therefore, this reflects the combined effect of the karst rocky desertification environment and bedrocks on the dynamic variation in the indicators of cave drip water [34,38], which results in a slightly lagging response to the precipitation event, as shown by the indicators. Especially, the Ca2+ and Mg2+ are also contributed by the detained water with full water rock interaction in the karst strata, from which the “piston effect” resulted in the less sensitive response of those hydrogeochemical indicators to karst rocky desertification and seriously interfered with the accuracy of the indicators reflecting the external eco-environment.
In order to further clarify the effect of the karst rocky desertification environment on the response of drip water indicators to atmospheric precipitation, the response time of the drip water in the rainy season for the Shijiangjun Cave was investigated and compared with those of other drip waters developed in different karst rocky desertification environments (Table 3). It is necessary to state that all the selected karst caves are located in the subtropical environment in South China to ensure a scientific and effective comparison, and the grade of karst rocky desertification overlying each cave is summarized according to the original description in the reference literature, as well as the response time data collection. The response time of drip water to the external atmospheric precipitation decreases with increases in the karst rocky desertification grade. Specifically, the drip water developed in the moderate-intense grade of karst rocky desertification generally responds to precipitation events rapidly within 1~2 days. Nevertheless, the drip water developed in the non-desertification and mild desertification regions requires at least a month to respond to precipitation events. The above-mentioned difference is mainly considered to be due to the overlying soil [38]. In general, the perennial type of drip water is principally supplied by atmospheric precipitation during the rainy season [24]. Comparatively, the soil layer developed in the moderate-intense grade of karst rocky desertification is thinner, with a poorer capacity for water conservation, such as the Xiaoyan Cave, Guangxi, South China. The soil structure is vulnerable and easily destroyed by the splash erosion of raindrops. As a result, the rainwater can quickly flow into the karst fissures and tunnels through the overlying soil layer. Furthermore, the carbonate rocks located in regions with intense degrees of karst rocky desertification are also characterized by a high level of dissolution and developed fissures [34], which leads to a relatively rapid response of the drip water to precipitation events.

4.4. Implications for the Paleoenvironment Reconstruction of Karst Rocky Desertification Using Speleothem

During the formation of drip water, the rainwater carries a large amount of calcium, magnesium minerals, and organic matter when it flows through the soil surface and into the underground space (i.e., pores, fissures, or caves). As a result, the formation of cave drip water is also a process of nutrient leakage and loss in karst regions. The soil and water leakage causes the loss of surface water and soil resources, leading to the damage of vegetation living conditions, a more fragile karst ecosystem, and the aggravation of karst rocky desertification [49]. Furthermore, the karst rocky desertification regions have unstable soil structure and an intense infiltration capacity of rainwater, in which the seepage water fills the karst fissures and dissolves the bedrocks, accelerating the development of fissures and pipelines, and transporting soil particles and nutrients to the underground space. Ultimately, it is inevitable that ecological degradation forms a vicious cycle in karst rocky desertification regions.
The surrounding environment of the Shijiangjun Cave is characterized by the steep slope of the peak body, high exposure rate of overlying rocks, a lower vegetation coverage rate and thinner soil layers. The fissures and karrens are extremely developed, which results in the strong connectivity of the flow paths and the main flow types of preferential flow and matrix flow [38] so the hydrogeochemical parameters of drip water rapidly respond to the external atmospheric precipitation. The retention time of moisture is relatively short and the water rock interaction is relatively weak in the strata due to the vertically developed fissure structures. Some geochemical indicators, such as SO42− and NO3, inherit the signals of the overlying environment directly. The results from a previous investigation indicated that the SO42− in drip water could be used as a tracer to indicate the surface biomass [43,45,50], which provides possibility and a reference for the reconstruction of the paleoenvironment using speleothems in the karst rocky desertification areas. However, the contents of cations (Ca2+, Mg2+ and Sr2+) in drip water coming from bedrock dissolution are inevitably and seriously affected by many hydrogeochemical factors, such as dilution effect, water sources, and prior calcite precipitation. The information carried by the cations in drip water may help researchers to understand the evolution process of karst rocky desertification. Therefore, they can explain the surface environment status and evolution process of karst rocky desertification on the premise of tracer mass sources so as to avoid the uncertainty generated from the interpretation of the overlying environment through the analysis of the hydrogeochemical indicators.
Some researchers have paid attention to the scientific values of organic indicators in their research on paleoclimate and paleoenvironment by stalagmites [51,52,53,54]. The organic matter in cave drip water is mainly derived from the overlying soil and vegetation, which maintain stable properties during the migration of seepage water. The signals carried by organic matter can effectively conserve the environmental information of overlying vegetation and soil. This preferably avoids the uncertainty in the response to climatic or environmental changes by using the stable isotopic composition of inorganic carbon and has the potential advantages of recording the changes in the climate and surface environment, as well as the important reference values in reconstructing the formation and evolutional history of karst rocky desertification using speleothems. Furthermore, the selection of appropriate drip sites is also of great importance for the speleothem−based paleoenvironment reconstruction of karst rocky desertification, in addition to the establishment of an appropriate indicator system. Overall, contemporary drip water that responds rapidly to precipitation is not suitable for the interpretation of paleoenvironmental information in karst rocky desertification regions reflected by the indicators in speleothems.

5. Conclusions

The hydrogeochemistry of drip water in the Shijiangjun Cave is controlled by the chemical dissolution of calcareous dolomite and shows the weak alkalinity of typical karst water. Both calcium and magnesium ions are dominant in the cations, and bicarbonate is dominant in the anions. All hydrochemical types of the selected three drip waters are classified as the HCO3-Ca·Mg type. According to the dynamic variations in the hydrogeochemistry of drip water over time, it is concluded that the response of W1 to precipitation was more delayed than those of W2 and W3; thus the three drip sites are grouped into two categories, including the delayed (W1) and rapid response types (W2 and W3), which resulted from the influence of different hydrogeochemical processes during the investigation. Specifically, the chemical compositions of W1 are constrained by the “piston effect”, and those of W2 and W3 are significantly affected by the rainwater dilution effect owing to the well-connected seepage flow paths. The relatively lower Ca2+ concentrations in the drip water of the Shijiangjun Cave during the hydrological year could be related to the effect of prior calcite precipitation, which needs to be verified by conducting a continuous and scientific observation in our following monitoring.In addition to the bedrock dissolution, the soil layer is one of the essential element sources for drip water, and the response sensitivity of the hydrogeochemical indicators for drip water to the external atmospheric precipitation is restricted by the karst rocky desertification degree during the rainy season. The influence of the overlying karst rocky desertification environment on the response of drip water to precipitation in the rainy season could be summarized as follows: the lower the degree of the karst rocky desertification, the more the response of the drip water is delayed; the higher degree of the karst rocky desertification, the faster the drip water responds to precipitation. This requires selecting more drip sites and conducting long-term and high-frequency observation for verification in the future. For the Shijiangjun Cave, several kinds of parameters, such as sulfate and nitrate in the drip waters of W2 and W3, have probably inherited the signals of the overlying soil environment and have the potential of indicating the degree of surface karst rocky desertification depending on the drip water hydrogeochemistry. However, a more accurate conclusion must be achieved on the basis of the stable isotope determination of nitrogen or sulfur in further exploration. Ultimately, it is worth noting that the selections of the drip sites and indicators in revealing the evolutional history of karst rocky desertification using speleothems should be paid more attention.

Author Contributions

Conceptualization, X.L. and K.X.; methodology, X.L.; software, X.L. and Y.L.; formal analysis, X.L. and Y.L.; investigation, X.L. and Y.L.; data curation, X.L. and Y.L.; writing—original draft preparation, X.L.; writing—review and editing, K.X., Y.L. and X.L.; visualization, X.L. and Y.L.; project administration, K.X.; funding acquisition, K.X. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the China-UNESCO Program of World Heritage Nomination and Conservation (No. Tianhe 20200728), the Key Project of Science and Technology Program of Guizhou Province (No. 5411 2017 Qiankehe Pingtai Rencai), the Science and Technology Program of Guizhou Province (No. Yiban 185 2021 Qiankehe Jichu-ZK) and the Ph.D. Research Startup Foundation of Guizhou Normal University in 2019 (No. 2 2019 GZNUD).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We appreciate the anonymous reviewers for their invaluable comments and suggestions on this manuscript, and we are particularly grateful to Wang Yubang from Qiyanqiao Town, Anshun City, for his warm help in the cave monitoring, we would also like to thank Zhang Yu and Zhan Zhaojun from School of Geographical Sciences, Southwest University (China), for their enthusiastic help in ion determinations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cao, J.H.; Yang, H.; Zhang, C.L.; Wu, X.; Bai, B.; Huang, F. Characteristics of structure and material cycling of the karst critical zone in Southwest China. Geol. Surv. China 2018, 5, 1–12, (In Chinese with an English Abstract). [Google Scholar]
  2. Xiong, K.N.; Li, P.; Zhou, Z.F.; An, Y.L.; Lyu, T.; Lan, A.J. Typical Study on RS and GIS of Karst Rocky Desertification; Geological Publish House: Beijing, China, 2002; pp. 17–28. (In Chinese) [Google Scholar]
  3. Huang, Q.H.; Cai, Y.L. Spatial pattern of karst rock desertification in the Middle of Guizhou Province, Southwestern China. Environ. Geol. 2007, 52, 1325–1330. [Google Scholar] [CrossRef]
  4. Wang, S.J.; Liu, X.M.; Zhang, D.F. Karst rocky desertification in southwestern China: Geomorphology, landuse, impact and rehabilitation. Land Degrad. Develop. 2004, 15, 115–121. [Google Scholar] [CrossRef]
  5. Tang, Y.Q.; Li, J.; Zhang, X.H.; Yang, P.; Wang, J.X.; Zhou, N.Q. Fractal characteristics and stability of soil aggregates in karst rocky desertification areas. Nat. Hazards 2013, 65, 563–579. [Google Scholar] [CrossRef]
  6. Zhou, Z.D.; Gao, T.; Zhu, Q.; Yan, T.T.; Li, D.C.; Xue, J.H.; Wu, Y.B. Increases in bacterial community network complexity induced by biocharbased fertilizer amendments to karst calcareous soil. Geoderma 2019, 337, 691–700. [Google Scholar] [CrossRef]
  7. Zhang, X.M.; Yue, Y.M.; Tong, X.W.; Wang, K.L.; Qi, X.K.; Deng, C.X.; Brandt, M. Eco-engineering controls vegetation trends in southwest china karst. Sci. Total Environ. 2021, 770, 145160. [Google Scholar] [CrossRef]
  8. Qian, C.H.; Qiang, H.Q.; Qin, C.Y.; Wang, Z.; Li, M.Y. Spatiotemporal evolution analysis and future scenario prediction of rocky desertification in a subtropical karst region. Remote Sens. 2022, 14, 292. [Google Scholar] [CrossRef]
  9. Tadros, C.V.; Treble, P.C.; Baker, A.; Fairchild, I.; Hankin, S.; Roach, R.; Markowska, M.; McDonald, J. ENSO–cave drip water hydrochemical relationship: A 7-year dataset from south-eastern Australia. Hydrol. Earth Syst. Sci. 2016, 20, 4625–4640. [Google Scholar] [CrossRef] [Green Version]
  10. Cheng, H.; Zhang, H.W.; Zhao, J.Y.; Li, H.Y.; Ning, Y.F.; Kathayat, G. Chinese stalagmite paleoclimate researches: A review and perspective. Sci. China Earth Sci. 2019, 10, 1489–1513. [Google Scholar] [CrossRef]
  11. Dreybrodt, W. The kinetics of calcite precipitation from thin films of calcareous solutions and the growth of speleothems: Revisited. Chem. Geol. 1981, 32, 237–245. [Google Scholar] [CrossRef]
  12. Baker, A.; Genty, D.; Dreybrodt, W.; Barnes, W.L.; Mockler, N.J.; Grapes, J. Testing theoretically predicted stalagmite growth rate with recent annually laminated samples: Implications for past stalagmite deposition. Geochim. Cosmochim. Acta 1998, 62, 393–404. [Google Scholar] [CrossRef]
  13. Fairchild, I.J.; Treble, P.C. Trace elements in speleothems as recorders of environmental change. Quat. Sci. Rev. 2009, 28, 449–468. [Google Scholar] [CrossRef]
  14. Coplen, T.B.; Winograd, L.J.; Landwehr, J.M.; Riggs, A.C. 500,000-years stable carbon isotopic record from Devils Hole, Nevada. Nature 1994, 263, 361–365. [Google Scholar] [CrossRef]
  15. Andrews, J.E.; Riding, R.; Dennis, P.R. The stable-isotope record of environmental and climatic signals in modern terrestrial microbial carbonates from Europe. Palaeogeogr. Palaeoclimatol. Palaeoecol. 1997, 129, 171–189. [Google Scholar] [CrossRef]
  16. Faimon, J.; Ličbinská, M.; Zajíček, P. Relationship between carbon dioxide in Balcarka Cave and adjacent soils in the Moravian Karst region of the Czech Republic. Int. J. Speleol. 2012, 41, 17–28. [Google Scholar] [CrossRef]
  17. Li, Y.D.; Yang, Y.; Jiang, X.Y.; Zhao, J.Y.; Sun, Z.; Shi, X.; Tian, N.; Yang, Y.Y.; Li, J.C.; Duan, J.W. The transport mechanism of carbon isotopes based on 10 years of cave monitoring: Implications for paleoclimate reconstruction. J. Hydrol. 2021, 592, 125841. [Google Scholar] [CrossRef]
  18. Liu, Z.Q. Applications of Geochemical Proxies in Speleothem to the Study on Evolution and Impact factor of Karst Desertification in Central Western Guizhou during Modern and Contemporary Period. Ph. D. Thesis, Southwest University, Chongqing, China, 2008. [Google Scholar]
  19. Duan, Y.F.; Liu, Z.Q.; He, Q.F.; Lv, X.F.; Jiang, J.J. Geochemistry of drip water in cave of rocky areas: A case study of Shijiangjun Cave in Guizhou Province. J. Southwest China Norm. Univ. (Nat. Sci. Ed.) 2015, 40, 95–101, (In Chinese with an English Abstract). [Google Scholar]
  20. Zhou, Y.C.; Wang, S.J.; Xie, X.N.; Luo, W.J.; Li, T.Y. Significance and dynamics of drip water responding to rainfall in four caves of Guizhou, China. Sci. Bull. 2005, 50, 155–162. [Google Scholar] [CrossRef]
  21. Domínguez-Villar, D.; Fairchild, I.J.; Baker, A.; Carrasco, R.M.; Pedraza, J. Reconstruction of cave air temperature based on surface atmosphere temperature and vegetation changes: Implications for speleothem palaeoclimate records. Earth Planet. Sci. Lett. 2013, S369–370, 158–168. [Google Scholar] [CrossRef]
  22. Rau, G.C.; Cuthbert, M.O.; Andersen, M.S.; Baker, A.; Rutlidge, H.; Markowska, M.; Roshan, H.; Marjo, C.E.; Graham, P.W.; Acworth, R.I. Controls on cave drip water temperature and implications for speleothem-based paleoclimate reconstructions. Quat. Sci. Rev. 2015, 127, 19–36. [Google Scholar] [CrossRef]
  23. Sheffer, N.A.; Cohen, M.; Morin, E.; Grodek, T.; Gimburg, A.; Magal, E.; Gvirtzman, H.; Nied, M.; Isele, D.; Frumkin, A. Integrated cave drip monitoring for epikarst recharge estimation in a dry Mediterranean area, Sif Cave, Israel. Hydrol. Process. 2011, 25, 2837–2845. [Google Scholar] [CrossRef]
  24. Guo, X.J.; Jiang, G.H.; Gong, X.P.; Yin, J.J.; Wu, X. Recharge processes on typical karst slopes implied by isotopic and hydrochemical indexes in Xiaoyan Cave, Guilin, China. J. Hydrol. 2015, 530, 612–622. [Google Scholar] [CrossRef]
  25. Guo, X.J.; Gong, X.P.; Yuan, D.X.; Jiang, G.H.; Cao, J.H.; Lin, Y.S.; Andrew, K.F.; Chen, C.J. Response of drip water temperature to climate variability: A case study in Xiaoyan Cave, southwest China. Hydrol. Sci. J. 2019, 64, 873–884. [Google Scholar] [CrossRef]
  26. Ban, F.M.; Pan, G.X.; Zhu, J.; Cai, B.; Tan, M. Temporal and spatial variations in the discharge and dissolved organic carbon of drip waters in Beijing Shihua Cave, China. Hydrol. Process. 2008, 22, 3749–3758. [Google Scholar] [CrossRef]
  27. Miorandi, R.; Borsato, A.; Frisia, S.; Fairchild, I.J.; Richter, D.K. Epikarst hydrology and implications for stalagmite capture of climate changes at Grotta di Ernesto (NE Italy): Results from long-term monitoring. Hydrol. Process. 2010, 24, 3101–3114. [Google Scholar] [CrossRef]
  28. Wu, X.; Pan, M.C.; Zhu, X.Y.; Cao, J.H.; Zhang, M.L. Effect of extreme precipitation events on the hydrochemistry index and stable isotope compositions of drip water in a subtropical cave, Guangxi, SW China. Carbonates Evaporites 2017, 33, 123–131. [Google Scholar] [CrossRef]
  29. Kelley, J.B.; Rowe, H.; Springer, G.S.; Guo, Y.L. Multi-year cave dripwater frequency and hydrochemical monitoring of three caves in Eastern North America. J. Cave Karst Stud. 2019, 81, 188–202. [Google Scholar] [CrossRef]
  30. Yang, M.D. Characteristics of karst landscape in Guizhou Plateau. J. Guizhou Norm. Univ. 1990, 2, 1–3. (In Chinese) [Google Scholar]
  31. Ford, D.C.; Williams, P.W. Karst Hydrogeology and Geomorphology; John Willy & Sons: Chichester, UK, 2007; pp. 45–52. [Google Scholar]
  32. Li, S.Y.; Lu, X.X.; Bush, R.T. Chemical weathering and CO2 consumption in the Lower Mekong River. Sci. Total Environ. 2014, 472, 162–177. [Google Scholar] [CrossRef]
  33. Kogovsek, J.; Petric, M. Solute transport processes in a karst vadose zone characterized by long-term tracer tests (the cave system of Postojnska Jama, Slovenia). J. Hydrol. 2014, 519, 1205–1213. [Google Scholar] [CrossRef]
  34. Luo, W.J.; Wang, S.J.; Xie, X.N. A comparative study on the stable isotopes from precipitation to speleothem in four caves of Guizhou, China. Chemie der Erde 2013, 73, 205–215. [Google Scholar] [CrossRef]
  35. Huang, Y.M.; Fairchild, I.J.; Borsato, A.; Frisia, S.; Cassidy, N.J.; Dermott, F.M.; Hawkesworth, C.J. Seasonal variations in Sr, Mg and P in modern speleothems (Grotta di Ernesto, Italy). Chem. Geol. 2001, 175, 429–448. [Google Scholar] [CrossRef]
  36. Dreybrodt, W.; Buhmann, D.; Michaelis, J.; Usdowski, E. Geochemically controlled calcite precipitation by CO2 outgassing: Field measurements of precipitation rates in comparison to theoretical predictions. Chem. Geol. 1992, 97, 285–294. [Google Scholar] [CrossRef]
  37. Poulain, A.; Rochez, G.; Roy, J.P.V.; Dewaide, L.; Hallet, V.; Sadelaer, G.D. A compact field fluorometer and its application to dye tracing in karst environments. Hydrogeol. J. 2017, 25, 1517–1524. [Google Scholar] [CrossRef] [Green Version]
  38. Luo, W.J.; Wang, S.J.; Xie, X.N.; Zhou, Y.C.; Li, T.Y. Temporal and spatial variations in hydrogeochemistry of cave percolation water and their implications for four caves in Guizhou, China. Chin. J. Geochem. 2013, 32, 119–129. [Google Scholar] [CrossRef]
  39. Li, J.Y.; Li, T.Y.; Wang, J.L.; Xiang, X.J.; Chen, Y.X.; Xuan, L. Characteristics and environmental significance of Ca, Mg, Sr in the soil infiltrating water overlying the Furong Cave, Chongqing, China. Sci. China Earth Sci. 2013, 56, 2126–2134. [Google Scholar] [CrossRef]
  40. Liu, W.; Zhou, C.; Burnet, J.E.; Brancelj, A. The effect of hydrological and hydrochemical parameters on the micro-distribution of aquatic fauna in drip water in the Velika Pasica Cave, central Slovenia. Ecohydrology 2017, 10, e1835. [Google Scholar] [CrossRef]
  41. Fairchild, I.J.; Borsato, A.; Tooth, A.F.; Frisia, S.; Hawkesworth, C.J.; Huang, Y.M.; McDermott, F.; Spiro, B. Controls on trace element (Sr-Mg) compositions of carbonate cave waters: Implications for speleothem climatic records. Chem. Geol. 2000, 166, 255–269. [Google Scholar] [CrossRef]
  42. Grosbois, C.; Negrel, P.; Fouillac, C.; Grimaud, D. Dissolved load of the Loire River: Chemical and isotopic characterization. Chem. Geol. 2000, 170, 179–201. [Google Scholar] [CrossRef]
  43. Tooth, A.F.; Fairchild, I.J. Soil and karst aquifer hydrological controls on the geochemical evolution of speleothem-forming drip waters, Crag Cave, southwest Ireland. J. Hydrol. 2003, 273, 51–68. [Google Scholar] [CrossRef]
  44. Spötl, C.; Fairchild, I.J.; Tooth, A.F. Cave air control on dripwater geochemistry, Obir Caves (Austria): Implications for speleothem deposition in dynamically ventilated caves. Geochim. Cosmochim. Acta 2005, 69, 2451–2468. [Google Scholar] [CrossRef]
  45. Ban, F.M.; Pan, G.X.; Cai, B.G.; Zhu, J.; Tan, M. Temporal-spatial variation of SO42− concentration of the drip water and its significance in the Shihua cave, Beijing. Carsologica Sin. 2009, 28, 243–248, (In Chinese with an English Abstract). [Google Scholar] [CrossRef]
  46. Sheng, M.Y.; Xiong, K.N.; Wang, L.J.; Li, X.N.; Li, R.; Tian, X.J. Response of soil physical and chemical properties to rocky desertification succession in South China Karst. Carbonates Evaporites 2018, 33, 15–28. [Google Scholar] [CrossRef]
  47. Wang, X.X.; Wu, Y.H.; Shen, L.C. Influences of air CO2 on hydrochemistry of drip water and implications for paleoclimate study in a stream-developed cave, SW China. Acta. Geochim. 2016, 35, 172–183. [Google Scholar] [CrossRef]
  48. Wu, X.; Zhu, X.Y.; Pan, M.C.; Zhang, M.L. Seasonal variability of oxygen and hydrogen stable isotopes in precipitation and cave drip water at Guilin, southwest China. Environ. Earth. Sci. 2014, 2, 3183–3191. [Google Scholar] [CrossRef]
  49. Gao, R.X.; Dai, Q.H.; Gan, Y.X.; Yan, Y.J.; Peng, X.D. The mechanisms of nutrient output through water flow from sloping farmland with slight rocky desertification in a karst region. Environ. Res. Lett. 2020, 15, 094085. [Google Scholar] [CrossRef]
  50. Onac, B.P.; Wynn, J.G.; Sumrall, J.B. Tracing the sources of cave sulfates: A unique case from Cerna Valley, Romania. Chem. Geol. 2011, 288, 105–114. [Google Scholar] [CrossRef]
  51. Baker, A.; Genty, D. Fluorescence wavelength and intensity variations of cave waters. J. Hydrol. 1999, 217, 19–34. [Google Scholar] [CrossRef]
  52. Tatár, E.; Mihucz, V.G.; Zámbó, L.; Gasparics, T.; Záray, G. Seasonal changes of fulvic acid, Ca and Mg concentrations of water samples collected above and in the Beke Cave of the Aggtelek karst system (Hungary). Appl. Geochem. 2004, 19, 1727–1733. [Google Scholar] [CrossRef]
  53. Li, X.L.; Hu, C.Y.; Liao, J.; Bao, L.L.; Mao, Q.X. An improved method for fluorescence analysis of dissolved organic matter in cave drip water. Front. Earth Sci. 2014, 8, 595–598. [Google Scholar] [CrossRef]
  54. Rutlidge, H.; Andersen, M.S.; Baker, A.; Chinu, K.J.; Cuthbert, M.O.; Jex, C.N.; Marjo, C.E.; Markowska, M.; Rau, G.C. Organic characterisation of cave drip water by LC-OCD and fluorescence analysis. Geochim. Cosmochim. Acta 2015, 166, 15–28. [Google Scholar] [CrossRef]
Figure 1. The selected drip sites in the Shijiangjun Cave (a): location of the Shijiangjun Cave in China; (b): longitudinal profile of cave sampling section; (c): planar graph of the sampling section in cave passage; (d): topographic map of the surrounding area, the dark-yellow region represents the planar graph of the Shijiangjun Cave, and the red dotted lines represent the cave sampling section; (e): speleothems and sediments in the Shijiangjun Cave.
Figure 1. The selected drip sites in the Shijiangjun Cave (a): location of the Shijiangjun Cave in China; (b): longitudinal profile of cave sampling section; (c): planar graph of the sampling section in cave passage; (d): topographic map of the surrounding area, the dark-yellow region represents the planar graph of the Shijiangjun Cave, and the red dotted lines represent the cave sampling section; (e): speleothems and sediments in the Shijiangjun Cave.
Ijerph 19 15830 g001
Figure 2. Piper diagram of drip water hydrochemistry.
Figure 2. Piper diagram of drip water hydrochemistry.
Ijerph 19 15830 g002
Figure 3. Variations of hydrogeochemical characteristics and the responses of drip rate to precipitation for three drip sites during precipitation events.
Figure 3. Variations of hydrogeochemical characteristics and the responses of drip rate to precipitation for three drip sites during precipitation events.
Ijerph 19 15830 g003
Figure 4. Time series of major hydrologic and hydrogeochemical characteristics for cave drip water in the Shijiangjun Cave (2014~2015).
Figure 4. Time series of major hydrologic and hydrogeochemical characteristics for cave drip water in the Shijiangjun Cave (2014~2015).
Ijerph 19 15830 g004
Figure 5. Comparisons between drip rates for the three selected drip sites during this rainfall monitoring (a) and the hydrological year from 2014 to 2015 (b).
Figure 5. Comparisons between drip rates for the three selected drip sites during this rainfall monitoring (a) and the hydrological year from 2014 to 2015 (b).
Ijerph 19 15830 g005
Figure 6. Spatial pattern of drip water hydrogeochemical indicators in the Shijiangjun Cave.
Figure 6. Spatial pattern of drip water hydrogeochemical indicators in the Shijiangjun Cave.
Ijerph 19 15830 g006
Figure 7. Correlation between Ca2+ and Mg/Ca at the three selected cave drip water sites.
Figure 7. Correlation between Ca2+ and Mg/Ca at the three selected cave drip water sites.
Ijerph 19 15830 g007
Table 1. Temporal-ordered parameters and descriptive statistics of the drip water hydrogeochemical composition in the Shijiangjun Cave.
Table 1. Temporal-ordered parameters and descriptive statistics of the drip water hydrogeochemical composition in the Shijiangjun Cave.
W1
(n = 18)
Drip Rate
/mL·min−1
pHEC
/μS·cm−1
HCO3
/mg·L−1
SO42−
/mg·L−1
NO3
/mg·L−1
Cl
/mg·L−1
K+
/mg·L−1
Na+
/mg·L−1
Ca2+
/mg·L−1
Mg2+
/mg·L−1
Sr2+
/mg·L−1
SICSID
2.218.41360213.521.563.623.550.950.5940.2220.970.0160.721.34
2.18.41370213.521.283.5551.50.594121.10.0160.731.35
2.058.41384213.521.013.527.092.280.5841.6421.270.0160.741.36
2.038.49354213.521.343.543.550.581.1243.1721.320.0190.831.53
2.028.52442207.420.73.1214.187.320.5839.6420.970.0160.811.52
2.168.78364207.420.143.385.671.220.5938.5320.790.0161.031.96
2.048.78362207.420.33.25.31.10.573920.70.0161.031.96
2.018.78361213.520.53.194.961.040.5639.520.650.0161.051.99
1.998.78351201.320.213.163.550.390.5739.5720.80.0151.031.96
1.898.54348201.320.623.13.550.40.7539.7120.730.0160.821.53
2.048.54353201.320.733.023.550.370.5840.0820.990.0160.821.53
2.248.7435718320.052.943.550.340.5440.1120.930.0160.961.81
2.18.7355170.820.12.853.550.40.741210.0170.911.7
1.858.61354170.820.132.843.550.430.7841.2821.280.0170.831.55
1.838.76351213.520.282.843.550.530.9840.2221.290.0161.041.99
1.958.57360219.620.892.477.83.050.6340.0221.250.0160.881.66
1.848.65355225.721.282.774.961.351.0540.4621.270.0170.971.84
1.918.67381170.821.272.5212.055.870.9441.8422.090.0170.891.67
Min1.838.41348170.820.052.473.550.340.5438.5320.650.020.721.34
Max2.248.78442225.721.563.6214.187.321.1243.1722.090.021.051.99
Average2.018.62364.56202.6620.693.095.51.620.7140.3921.080.020.891.68
S.D0.120.1421.6617.20.50.343.081.970.191.120.340.00080.110.23
W2
(n = 18)
Drip Rate
/mL·min−1
pHEC
/μS·cm1
HCO3/mg·L1SO42/mg·L1NO3/mg·L1Cl
/mg·L1
K+
/mg·L1
Na+
/mg·L1
Ca2+
/mg·L1
Mg2+
/mg·L1
Sr2+
/mg·L1
SICSID
2.038.66274158.615.672.83.550.230.5122.1919.430.0080.61.31
1.958.7275170.815.42.762.130.230.5222.719.70.0080.671.46
1.88.74275170.815.052.751.420.220.5322.9219.890.0080.711.54
1.758.68272152.516.212.724.250.320.6423.7219.960.0090.631.36
1.828.78272170.814.912.73.550.210.5323.4319.820.0080.761.62
1.658.82273170.815.072.744.250.240.5422.4219.720.0080.771.67
1.598.8274176.915.072.83.550.210.532319.80.0080.781.67
1.458.78274176.915.12.852.840.20.5223.2919.90.0080.771.65
1.418.85271170.815.912.722.130.220.5623.0619.820.0080.811.73
1.388.7270176.914.982.74.250.20.5323.6220.040.0080.71.51
1.718.7273158.614.942.642.130.20.5323.3819.890.0080.661.42
2.18.7527418315.042.652.840.220.5723.6119.860.0090.761.62
2.018.7272176.9162.682.840.210.5423.519.80.0080.71.5
2.028.7270176.916.492.713.550.210.5323.2619.730.0080.71.5
2.038.77271170.815.392.692.840.190.6623.6720.130.0090.751.61
2.098.73272170.816.382.72.130.210.5423.6320.110.0080.711.54
1.988.74273176.917.262.732.130.20.5423.5719.850.0090.741.58
1.438.78272176.916.882.712.840.271.0624.8319.90.010.791.67
Min1.388.66270152.514.912.641.420.190.5122.1919.430.010.61.31
Max2.18.8527518317.262.854.250.321.0624.8320.130.010.811.73
Average1.798.74272.61171.4815.652.722.950.220.5823.3219.850.010.721.55
S.D0.250.051.547.80.740.050.850.030.130.580.160.00060.060.12
W3
(n = 18)
Drip Rate
/mL·min−1
pHEC
/μS·cm1
HCO3/mg·L1SO42/mg·L1NO3/mg·L1Cl
/mg·L1
K+
/mg·L1
Na+
/mg·L1
Ca2+
/mg·L1
Mg2+
/mg·L1
Sr2+
/mg·L1
SICSID
2.378.7232170.85.81.913.550.230.3116.8619.280.0060.561.35
2.258.72232164.75.621.930.260.311719.50.0060.561.36
2.188.75232158.65.591.92.840.290.3217.0419.570.0060.581.39
28.71232164.74.991.862.840.260.316.7719.490.0060.551.34
1.728.86232207.46.331.883.550.230.3316.7819.240.0060.771.78
1.648.812311836.431.872.130.250.4919.1319.660.0080.731.65
1.878.82321836.11.812.830.220.351719.40.0070.671.59
1.948.792321835.831.773.550.20.316.5919.270.0060.661.56
2.118.84230164.75.91.812.840.210.5916.8419.150.0070.671.56
2.38.73230176.95.371.893.550.20.3116.6619.140.0060.591.42
2.518.742311836.051.782.840.210.316.7919.380.0060.621.47
2.18.77231164.75.91.792.840.20.3216.5719.270.0060.61.44
2.328.74232176.96.11.793.20.290.41719.70.0070.611.46
2.558.74232176.96.271.83.550.370.5219.1519.80.0080.661.51
2.618.78230176.95.841.862.840.210.316.8419.560.0060.641.52
2.838.76231176.96.671.83.550.220.3116.7619.350.0060.621.48
2.218.75232158.66.211.812.130.220.4917.4419.60.0070.591.4
2.458.762311836.421.861.420.220.3116.3518.810.0060.621.49
Min1.648.7230158.64.991.771.420.20.316.3518.810.010.551.34
Max2.838.86232207.46.671.913.550.370.5919.1519.80.010.771.78
Average2.228.76231.39175.215.971.842.940.240.3617.0919.40.010.631.49
S.D0.310.040.7811.80.410.050.590.040.090.780.240.00060.060.11
Table 2. Previous hydrochemical data compilation of rainwater during rainy seasons in the study area (unit: mg·L−1).
Table 2. Previous hydrochemical data compilation of rainwater during rainy seasons in the study area (unit: mg·L−1).
DatepHK+Na+Ca2+Mg2+Sr2+HCO3SO42−NO3Cl
12 May 20147.220.320.112.270.60.0038.71.40.880.12
17 June 20147.130.250.132.80.540.00210.41.480.450.3
9 August 20147.350.370.22.410.630.00312.21.330.90.38
15 June 20157.360.330.182.60.440.00211.30.891.010.31
14 July 20157.410.280.151.830.50.00311.21.420.640.22
Table 3. Comparison of the drip water response time to precipitation among different grades of karst rocky desertification regions in the rainy season.
Table 3. Comparison of the drip water response time to precipitation among different grades of karst rocky desertification regions in the rainy season.
CaveSampling
Period
LocationVegetation
Type
Soil
Thickness/cm
Karst
Desertification
Grade
Response
Time
References
Shijiangjun2016Anshun,
Guizhou
Shrub, grass5~46Moderate38~42 hThis study
Qixing 2003Duyun, GuizhouShrub50~90Mild27~40 d[20]
Liangfeng 2003Libo, GuizhouPrimordial forest100None36 d[20]
Xiniu 2003Zhenning,
Guizhou
Shrub grass20~60Mild28 d[20]
Xiaoyan2014Guilin, GuangxiShrub grass30~100Moderate-intense<48 h[24,25]
Furong2010Wulong,
Chongqing
Arbor, shrub>30None>30 d[39]
Xueyu2012~2013Fengdu,
Chongqing
Arbor, shrub0~50None30 d[47]
Panlong2012Guilin, GuangxiShrub30~100Mild15 d[48]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lyu, X.; Li, Y.; Xiong, K. Hydrogeochemical Response of Cave Drips to Precipitation during Rainfall in a Karst Desertification Region: A Case Study of Shijiangjun Cave, South China. Int. J. Environ. Res. Public Health 2022, 19, 15830. https://doi.org/10.3390/ijerph192315830

AMA Style

Lyu X, Li Y, Xiong K. Hydrogeochemical Response of Cave Drips to Precipitation during Rainfall in a Karst Desertification Region: A Case Study of Shijiangjun Cave, South China. International Journal of Environmental Research and Public Health. 2022; 19(23):15830. https://doi.org/10.3390/ijerph192315830

Chicago/Turabian Style

Lyu, Xiaoxi, Yuan Li, and Kangning Xiong. 2022. "Hydrogeochemical Response of Cave Drips to Precipitation during Rainfall in a Karst Desertification Region: A Case Study of Shijiangjun Cave, South China" International Journal of Environmental Research and Public Health 19, no. 23: 15830. https://doi.org/10.3390/ijerph192315830

APA Style

Lyu, X., Li, Y., & Xiong, K. (2022). Hydrogeochemical Response of Cave Drips to Precipitation during Rainfall in a Karst Desertification Region: A Case Study of Shijiangjun Cave, South China. International Journal of Environmental Research and Public Health, 19(23), 15830. https://doi.org/10.3390/ijerph192315830

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