Abstract
The hydrochemical characteristics and variation rules of water in karst areas are important for the basic research of karst field. Because only a few researches focused on the hydrochemical characteristics of humid subtropical dolomite karst basins in the past, the Shanmuhe Basin and its subbasin Huangzhouhe Basin in Shibing Dolomite Karst World Natural Heritage Site in China were selected as the study area. A total of 215 water samples were collected and analyzed. It was found that the hydrochemical type of subtropical dolomite karst basin is HCO3-Ca·Mg type, which is different from that in limestone areas and shale districts, which is HCO3-Ca type for both; EC, Ca2+, Mg2+, HCO3-, Mg2+/Ca2+, and Mg2+/HCO3- have similar spatial variations, which are lower in allogenic water in the upstream in the northwest of the basin and become higher when the water enters the dolomite karst area in the mid and downstream of the basin, and SIC and SID of the allogenic water in the upstream are below zero, while those become saturated after entering the dolomite area in the mid and downstream; the monthly variation of hydrochemistry shows that the seasonal variation of ion concentration is not obvious, and there is a dilution effect in rainy seasons, but the variation of ion concentration is small and relatively stable; higher HCO3- concentration and Mg2+/Ca2+ in the pure dolomite basin are higher than those in the impure dolomite basin; the main influencing factors of the hydrochemical features are geological background especially the lithology, and the water chemistry is mainly controlled by rock weathering and water-rock interaction; rainfall import is low and the low concentrations of K+, Na+, NO3-, Cl-, and SO42- in water indicate that human activities have little effect on the hydrochemistry in the study area.
1. Introduction
The influence of water dissolution of carbonate rock will be directly reflected by the hydrochemical characteristics of water. The chemical weathering of rocks and soils occurs under the leaching of rainwater, and leaching ions enter rivers. Therefore, the hydrochemical composition of rivers has become the main way to study the chemical weathering process of rocks recently and has been widely used in the world. By studying the hydrochemical characteristics of rivers, we can get the exogenous circulation process of chemical weathering, climate, and average chemical elements of upper crust in the continent-river-ocean system [1]. The hydrochemical types, compositions, temporal and spatial variations, and ion sources of surface rivers have been studied in depth [2–5]. As early as the 1970s, some scholars studied the hydrochemical characteristics of karst basins dominated by limestone and gypsum [6–8]. In China, the geochemistry of surface water and groundwater in karst basins was also studied. In surface river basins, the characteristics and sources of hydrochemistry and isotopes were mainly studied [1, 9]; and the regional hydrochemical differences between surface water and groundwater and their causes were discussed from the regional scale [10, 11]. From the perspective of the impact of land use and human activities on hydrochemistry, water quality, and stable isotopes of dissolved inorganic carbon in karst groundwater, more extensive and in-depth studies were conducted [12–17]. Based on preliminary identification of karst distribution and hydrogeochemical characteristics of typical karst areas in Thailand, Zhang et al. [18] analyzed the main characteristics and dynamic conditions of karst development in Thailand and compared them with those in southwest China. The main karst strata in Thailand are Ordovician, Permian, and Triassic limestone. However, most of the researches are focusing on limestone karst areas, and studies of dolomite karst area are insufficient. Researches on the hydrochemistry of dolomite karst in tropical-subtropical climatic conditions are more insufficient. Cowell et al. [19] studied the hydrochemistry of dolomite karst area under climate in temperate zone and pointed out that the hydrochemical type in this area is HCO3-Ca·Mg type; calcite saturation index (SIC) and dolomite saturation index (SID) of surface water are both above 0, which represents mature landforms. Zhang [20] studied the dolomite in Luocheng, Guangxi, China, and pointed out that the hydrochemical type in this area was HCO3-Ca·Mg type; when dolomite dissolves, calcite leaches preferentially. Dolomite is as commonly distributed as limestone in the world [21]. The dolomite and dolomite-limestone interbedding lithology account for 25.7% of the total distribution area of carbonate rocks in China. Therefore, more attention should be given on the research on dolomite karst areas.
Under the background of serious large-scale rocky desertification in karst areas of Guizhou, China, Shibing has become a good place for the study of dolomite karst in China and even in the whole world because of the preservation of the natural dolomite karst process in a relatively pure dolomite geological background with primary forest coverage. Yunnan, Guangxi, Hunan, and Hubei also have large areas of dolomite, but most of them are dolomite with limestone interbedded or limestone dolomite, while Shibing retains pure dolomite with a total thickness of 1100 m. Therefore, Shibing, Guizhou, in China was finally selected as the most typical example of the subtropical karst area developed on dolomite in the world and was successfully inscribed on the World Heritage List in 2014 as an irreplaceable component of the South China Karst. Shibing Karst, which has the best known example of subtropical fengcong karst in dolomite, provides a spectacular fengcong landscape, which is also exceptional because it developed in relatively insoluble dolomite rocks (UNESCO World Heritage Committee, 2014). It is one of the ideal and natural places to study dolomite karst. In this paper, the Shanmuhe (Shanmu River) Basin and its subbasin Huangzhouhe (Huangzhou River) Basin in Shibing, Guizhou Province, were selected as the study areas, and the hydrogeochemical characteristics of dolomite rock weathering were revealed through manual survey and intensive monitoring of hydrochemistry in the basin. On the basis of dividing the underground basin of 1 : 50000 and 1 : 200000 hydrogeological maps in the study area, combined with 1 : 50000 and 1 : 10000 topographic maps, through field hydrogeological survey and hydrochemical analysis, Shanmuhe and Huangzhouhe basins were defined. The main springs and surface water from upstream to downstream of Shanmuhe Basin were sampled jointly by hydrochemical method, the differences analyzed, and the controlling factors of weathering and karstification in dolomite basin discussed as a function of lithology, precipitation, and human activities.
2. Study Area
The study area, including Shanmuhe Basin and its subbasin Huangzhouhe Basin (Figure 1) in Shibing Dolomite Karst World Natural Heritage Site, is located in Shibing County, Guizhou Province, China, which is the hilly slope zone of the transition from Yunnan-Guizhou plateau to western Hunan hills. The total area of Shanmuhe Basin is 242.00 km2, and the area of dolomite outcrop is 192.34 km2, accounting for 80% of the total area of the basin and making it a typical dolomite basin. The area of Huangzhouhe Basin, the subbasin of Shanmuhe Basin, is about 49.19 km2, which is not affected by allogenic water. Except for a small amount of Ordovician limestone and argillaceous limestone in the lower part of the Cambrian Qingxudong Formation, the vast majority of the basin (45.27 km2, 92%) is Cambrian Shilengshui Formation dolomite, trans-Cambrian-Ordovician Loushanguan Formation dolomite, and a small amount of fine-grained dolomite. It has become an ideal basin for the study of subtropical dolomite karst process nationwide and even globally, with a high degree of typicality and representativeness.
3. Methodology
3.1. Water Sample Collection
The purpose is to study the hydrochemical characteristics of subtropical dolomite karst area. According to the geological background and hydrogeological characteristics of the basin and considering the influence of lithology on water, groundwater and surface water were sampled from 39 points (Figure 1 and Table 1) from upstream to downstream of the Shanmuhe Basin from May 2014 to July 2016, including 14 points of groundwater (Nos. 26-39, 1 allogenic groundwater sample, 12 karst groundwater samples) and 25 points of surface water (Nos. 1-25: Nos. 1-22 are karst water samples, and Nos. 22-24 are allogenic water). Two samples of atmospheric precipitation were obtained during the period. Most of these water are sampled monthly or seasonally to understand the basic hydrochemical characteristics, influencing factors, and temporal and spatial variations of the dolomite basin. On the basis of surface water sample collection and analysis, the key point is to select No. 4 outlet of Huangzhouhe Basin for intensive collection and monitoring. A total of 215 water samples were collected.
Sampling method: at the sampling site, water temperature (), pH, conductivity (EC), and dissolved oxygen (DO) were measured by French PONSEL ODEON portable digital water quality measuring recorder, with the accuracy of 0.1°C, 0.01 pH unit, 1 μS/cm, and 0.01 mg/L, respectively. The calibration methods include the following. (1) The temperature sensor is calibrated annually, in a 2-step process: the sensor is placed in a container containing a water bath and ice, and then, the sensor is placed in a medium (air or water in a temperature-controlled bath) at a known temperature. This temperature may be measured using a certified thermometer. (2) Oxygen calibration: using a clean sensor, occasionally check the 0% Sat value by immersing the sensor in a sulphite solution in water (sulphite by weight). If there is a zero error, perform a complete calibration of the sensor. The 2-point calibration is performed using a sulphite solution (for the offset); then, after rinsing and drying, the sensor’s gain is determined by exposing the sensor to humid, oxygen-saturated air. (3) The calibration of the conductivity sensor is a 2-step process: step 1 (offset): the sensor is exposed to air; step 2 (gain): the sensor is placed in a buffer solution of known conductivity. (4) The calibration of the pH sensor is a 2-step process: step 1 (offset): the sensor is placed in a pH 7.01 buffer solution; step 2 (gain): the sensor is placed in a pH 4.01 (or 9.01 or 10.01) buffer solution. The concentration of HCO3- and Ca2+ in water was measured by alkalinity meter and hardness meter produced by Merck Company, Germany, respectively. The accuracy is 0.1 mmol/L and 2 mg/L.
Cationic water samples were collected by 50 mL high-density polyethylene plastic bottles cleaned before, and anionic water samples were collected by 100 mL high-density polyethylene plastic bottles. Among them, cationic water was acidified with 1 : 1 HNO3 until . The cap was tightened and sealed with adhesive tape. The above samples were sealed and stored in the incubator at low temperature about 4°C and then brought back to the laboratory for analysis and determination.
3.2. Water Sample Analysis
The analysis of water samples mainly includes the main physical and chemical indexes (Table 2).
Anions were determined by ICS 90 type ion chromatograph produced by Dionex Company in the United States, while cations were determined by ICP-OES with a detection accuracy of 0.01 mg/L (relative standard ).
Calcite saturation index (SIc), dolomite saturation index (SId), and partial pressure of carbon dioxide () were calculated on Phreeqc and Watspec software by using , pH, and main ion concentrations of water samples.
The above tests were completed in the Water Chemical Analysis Laboratory of Southwest University and the State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences.
4. Results
4.1. Hydrochemical Characteristics
4.1.1. Hydrochemical Types
In order to study the hydrochemical types of surface water and groundwater in the Shanmuhe Dolomite Basin, Shibing, Guizhou, the corresponding Piper triangle diagram was made according to the ion content. The Piper diagram can better illustrate the composition and distribution characteristics of ions in water and can also distinguish the source of major and minor ions in water and the control of rock weathering.
From the Piper diagram (Figure 2) of the study area, HCO3- is the main anion in the karst surface water and groundwater. Its average values are 236.6 mg/L and 336.1 mg/L, accounting for 91.7% and 93.7% of the total anions, respectively. The concentration and proportion of HCO3- in groundwater are higher than those in surface water, which is related to its longer water-rock interaction time. Ca2+ and Mg+ are the main cations in karst surface water and groundwater of the study area. The average concentration of Ca2+ in surface water and groundwater is 45 mg/L and 65 mg/L, respectively, accounting for 63.71% and 65.24% of the total cations, and the average concentration of Mg+ in the surface water and groundwater is 25 mg/L and 32 mg/L, respectively, accounting for 32.50% and 34.84% of the total cations, respectively. The remaining ions only account for less than 3%. Therefore, the chemical types of karst surface water and groundwater in the study area are HCO3-Ca·Mg. The concentrations of Ca2+, Mg2+, and HCO3- are closely related to karstification, reflecting the influence of dolomite dissolution on hydrochemical characteristics in the study area.
The ion content of surface water and groundwater in the allogenic water is very low, with HCO3- as the main anion, the average value of which is 36.6 mg/L and 18.3 mg/L, accounting for 71.96% and 69.48% of the total anion, while the average concentration of Ca2+ as the main cation is 13 mg/L and 6 mg/L, accounting for 74.29% and 64.28% of the total cation, respectively, and the Mg2+ content is very low. The hydrochemical type of the allogenic water in the study area is HCO3-Ca, but the proportion of Na+, K+, NO3-, Cl-, and SO42- is higher than that of karst surface water and groundwater. Only low concentration of Ca2+, NO3-, and SO42- was detected in atmospheric precipitation.
Gibbs’ semilogarithmic coordinate diagrams can be used to analyze the main sources of hydrochemistry of rivers [23]. As shown in Figure 3, the longitudinal coordinate of Gibbs diagram is TDS (logarithmic representation, mg/L), and the horizontal ordinate represents the mass concentration ratio of cationic Na+/(Na++Ca2+) and anionic Cl-/(Cl-+HCO3-), respectively. Almost all the ionic component values of surface water in the world fall into the imaginary frame in the figure. Points with low TDS and high Na+/(Na++Ca2+), Cl-/(Cl-+HCO3-) ratios (close to 1) are distributed in the lower right corner, reflecting regional precipitation control; points with medium TDS and Na+/(Na++Ca2+), Cl-/(Cl-+HCO3-) ratios (<0.5) are located on the left side of the diagram, reflecting the influence of rock weathering; and higher TDS and Na+/(Na++Ca2+), Cl-/(Cl-+HCO3-) ratios (close to 1) reflect the influence of evaporation and concentration.
The diagram indicates that the distribution of TDS in Shibing Shanmuhe Dolomite Basin is as follows: karst surface water 164.4-282.4 mg/L with average 210.0 mg/L, karst groundwater 234.8-341.8 mg/L with average 290.7 mg/L, and allogenic water 21.5-88.8 mg/L with average 42 mg/L. The Na+/(Na++Ca2+) value is 0-0.034 for karst surface water with average 0.009, 0.003-0.057 for karst groundwater with average 0.013, and 0.013-0.212 for allogenic water with average 0.118. The ratio of Cl-/(Cl-+HCO3-) is 0-0.183 for karst surface water with average 0.017, 0.004-0.032 for karst groundwater with average 0.016, and 0.008-0.366 for allogenic water with average 0.087. All water points fall within the range of Na+/(Na++Ca2+) or Cl-/(Cl-+HCO3-) ratios less than 0.5 and are distributed on the left middle side of the figure, reflecting that surface water and groundwater in Shibing dolomite karst basin belong to the “rock weathering type.” Ion composition in water mainly comes from the weathering process of rocks, and weathering of rocks has a significant impact on the hydrochemistry in this area. The ratios of Na+/(Na++Ca2+) or Cl-/(Cl-+HCO3-) in some samples of allogenic water are higher, but the ratio is less than 0.5, reflecting that some water points of allogenic water may also be affected by atmospheric precipitation.
4.1.2. Overall Hydrochemical Characteristics
The composition of water system in Shibing dolomite karst basin is divided into atmospheric precipitation, karst surface water, karst groundwater, and allogenic water.
(1) Atmospheric Precipitation. Two precipitation samples were obtained in the basin. The pH of the precipitation in the basin is 5.0 and 5.3, respectively, with an average value of 5.2. The conductivity of the two rainfalls was 22 μS/cm and 25 μS/cm, respectively, with an average of 23 μS/cm. The concentration of ions in precipitation is very low, even some ions are not detected, only Ca2+ and SO42- with the mean values of 0.03 mmol/L and 0.05 mmol/L, respectively. The concentration of Mg2+, K+, and Cl- is very low, and Na+ and HCO3- cannot be detected.
(2) Karst Surface Water. The karst surface water in the basin was taken from the main trunk and tributaries of the Shanmuhe from the upstream to the downstream of the basin. The variation range of water temperature is 4.1-24.0°C. The standard deviation (STD) and coefficient of variation (CV) are large, indicating that the change of water temperature is obvious, mainly affected by the change of surface atmospheric temperature. The pH ranges from 7.92 to 9.44, with an average of 8.67, belonging to alkaline water. The conductivity ranges from 320 μs/cm to 595 μs/cm, with an average of 374 μs/cm. The variation range is large, but the CV is only 9.15%. This indicates that the conductivity in surface water is relatively stable. The concentrations of Ca2+ and Mg2+ in surface water are 0.70~1.55 mmol/L and 0.63~1.49 mmol/L, respectively. The average values are 1.12 mmol/L and 1.02 mmol/L, respectively. The concentration of the two ions has little change, and the ion concentration is at the same level. The variation range of Mg2+/Ca2+ is 0.64~1.30 with the average of 0.92, which is close to 1, and the variation is not significant, indicating that dolomite dissolution is dominant in the basin. Meanwhile, the mean value of (Ca2++Mg2+)/HCO3- is 0.56, which is equivalent to 2HCO3- = Ca2++Mg2+, indicating that carbonate controls the strong dissolution of dolomite. The contents of K+ and Na+ ions are very low, and the mean sum of them is only 0.03 mmol/L. The anion is mainly HCO3- with a value of 2.90-5.20 meq/L, a mean of 3.87 meq/L, and a CV of 11.86%. This indicates that the temporal and spatial distribution of HCO3- concentration of surface karst water in the basin is relatively stable under the same climate and geological background. SO42-, Cl-, and NO3- ions are all very low, and the sum of them accounts for less than 10% of the anions.
(3) Karst Groundwater. The temperature of karst groundwater in the basin ranges from 11.7 to 23.3°C, with an average of 16.7°C, and the variation range is less than that of surface water. The pH ranges from 7.15 to 8.82, with an average of 7.74. The conductivity ranges from 219 μs/cm to 612 μs/cm, with an average of 496.8 μs/cm. Compared with the karst surface water, the groundwater has lower pH and higher conductivity. On the one hand, lower pH promotes dissolution; on the other hand, groundwater has a relatively long time of rock-water interaction, which leads to higher ion concentration in groundwater. Ca2+, Mg2+, and HCO3- are the main ions in the groundwater, and the average concentrations are 1.57 mmol/L, 1.28 mmol/L, and 5.25 meq/L, respectively. The concentration of Ca2+ is slightly higher than that of Mg2+. The range of Mg2+/Ca2+ varies from 0.49 to 1.0, with an average of 0.81, which may be the result of ratio of limestone dissolution in groundwater source. The average value of (Ca2++Mg2+)/HCO3- is 0.55, which is equivalent to 2HCO3- = Ca2++Mg2+, indicating that carbonate controls the strong dissolution of dolomite in the basin. Similar to karst surface water, K+, Na+, SO42-, Cl-, and NO3- ions in karst groundwater are very low. This is mainly because, on the one hand, Shibing Shanmuhe Basin is a world natural heritage site with good ecological environment and low ion composition of precipitation, which can neglect the input of surface water and groundwater in the basin; on the other hand, the lithology of the study area is dominated by dolomite.
(4) Allogenic Water. The nonkarst area in the basin is mainly located in the northwest part. Its strata and lithology mainly include the clay shale of Cambrian Palang Formation, the carbonaceous shale and sandstone of Bianmachong Formation, the siliceous rock and carbonaceous shale of Sinian Laobao Formation, and the dolomite, argillaceous dolomite, shale, and claystone of Doushantuo Formation. The pH of the allogenic water ranges from 6.68 to 8.21, with an average of 7.54, which is lower than that of the karst water. The conductivity is only 42.9-181.2 μS/cm, which is lower than that of the karst water. Ca2+ and HCO3- are the main ions, and their average values are 0.26 mmol/L and 0.49 mmol/L, respectively. Then, the concentration of Mg2+ and SO42- is 0.08 mmol/L, the average concentration of NO3-+Cl- is 0.09 mmol/L, and Na++K+ is 0.06 mmol/L. The average values of Mg2+/Ca2+ and (Ca2++Mg2+)/HCO3- are 0.34 and 0.77, respectively, which are not consistent with the results of carbonate rock dissolution.
4.2. Variation Rules of Hydrochemistry
4.2.1. Spatial Variation of Hydrochemistry
Dolomite is the main lithology in Shibing Shanmuhe Basin, but there are many different strata in the basin. In the northwest part of the basin, it is influenced by noncarbonate rocks such as sandstone and shale of Cambrian Bianmachong Formation and Palang Formation as well as impure carbonate rocks of Qingxudong Formation and also impure dolomite of Shilengshui Formation. The distribution of impure dolomite in Shilengshui Formation and pure dolomite in Loushanguan Formation is dominant in the whole basin, and the Huangzhouhe subbasin is the typical pure dolomite basin. Because the basin is dominated by surface rivers, there are few big springs and many small springs, but the water quantity of springs is small. Therefore, the authors only discuss the spatial variation characteristics of surface water.
Samples of river water were collected in the main stream and tributaries of the Shibing Shanmuhe Basin. Considering that the main stream, in the lower reaches of the Shibing Shanmuhe Basin, is affected by the river water from the upper reaches, while the tributary river water is not affected by the upper main stream, its hydrochemistry represents the characteristics of the subbasin; the distribution of lithology is also taken into account. Therefore, the hydrochemical data of the tributaries and tributary water sources of the Shanmuhe are presented below for analysis. No. 3, which represents the main stream of Shanmuhe, No. 4, which represents Huangzhouhe Basin, No. 6, which represents Jiangjiatian tributary influenced by pure dolomite, No. 9, which represents the west of Liangcha River, No. 10, which represents the east of Liangcha River, Nos. 15-20 (mean value) and No. 21 of Dolomite Basin in Shandong River, and Nos. 22-24 of allogenic water in the upper and northwest of the basin were selected for analysis (the distribution of sample points is shown in Figure 1).
Figure 4 shows that EC, Ca2+, Mg2+, HCO3-, Mg2+/Ca2+, and Mg2+/HCO3- have similar spatial variation, mainly due to different lithology controls. The average EC concentration of the allogenic water (No. 22-24 sample points) northeast of the upper reaches of the basin was , , and 181 μS/cm, respectively; the average Ca2+ concentration is , , and 0.66 mmol/L; the average Mg2+ concentration is , , and 0.14 mmol/L; the average HCO3- concentration is , , and 1.40 mmol/L; the concentration of these ions is significantly lower than that of the surface river sampling points in the karst area. The concentration of No. 24 sample is slightly higher than that of No. 22 and No. 23, mainly because No. 24 is affected not only by clay shale of the Palang Formation but also by impure carbonate rocks of the Qingxudong Formation. The flow of No. 22 and No. 23 samples reaches No. 10 after passing through the dolomite area, and No. 24 flows into the Shanmuhe at No. 9. The concentrations of EC, Ca2+, Mg2+, and HCO3- are significantly increased, mainly because the SIC and SID of the allogenic water are both lower than 0 and have strong dissolution capacity, which can dissolve the hosting carbonate rocks. SIC and SID of surface rivers in all karst areas are >0.8, indicating that carbonate deposits will occur when they reach supersaturation. In addition, the water at the No. 15-21 sample points in the northeast of the upper reaches of the basin will also affect the ion concentration at the No. 10 sample point. Due to the influence of the upstream allogenic water, the concentrations of EC, Ca2+, Mg2+, and HCO3- at No. 9 and No. 10 sample points are lower than those at other karst surface rivers. Generally, the upstream ion concentration is lower than that in the middle and lower reaches. The small basin where No. 4, No. 6, and No. 21 sample points are located is pure dolomite, and the average HCO3- concentration is , , and , respectively; the average Mg2+/Ca2+ concentration is , , and , respectively. No. 3, 9, and 10 sample points belong to impure dolomite basins. The average HCO3- concentration is , , , and , and the average Mg2+/Ca2+ concentration is , , , and , respectively. These results indicate that the concentrations of HCO3- and Mg2+/Ca2+ in pure dolomite basin are generally higher than those in dolomite and limestone basin. It shows that although the dissolution rate of dolomite is slower than that of limestone, its higher solubility makes it produce more HCO3-.
4.2.2. Monthly Variation of Hydrochemistry
In order to illustrate the hydrochemical monthly variation characteristics of typical dolomite karst basins in subtropical zone, two representative samples (No. 4 and No. 29) were selected, of which No. 4 was the overall outlet of Huangzhouhe Basin, representing the surface water of dolomite karst basin, and No. 29 YYD was a karst spring representing the groundwater of dolomite karst basin. In Figure 5, the main hydrochemical analysis of water samples taken monthly at these two sampling points is reported. According to the analysis above, the contents of K+, Na+, SO42-, Cl-, and NO3- ions in the basin are very low. Therefore, the monthly variation characteristics of Ca2+, Mg2+, and HCO3- and related indicators in the basin are discussed here.
From Figure 5, it can be seen that the pH of the outlet water of Huangzhouhe varies from 8.40 to 8.80, which is alkaline, low in rainy season and high in dry season. The conductivity varies between 397 μs/cm and 424 μs/cm, and it is smaller in rainy season than in dry season, but the variation coefficient is 0.03, which indicates that the dilution effect is not obvious in rainy season. Correspondingly, the changes of Ca2+, Mg2+, and HCO3- are also small, and their concentrations are between 42-54 mg/L, 20-31 mg/L, and 238-347 mg/L. The variation coefficients are 0.07, 0.17, and 0.11, respectively. The variation rule is that the concentration in rainy season is smaller than that in dry season. The saturation indexes SIC and SID in river water have the same trend of change, and the change is not significant, and they reach the state of supersaturation every month. The results show that although the dissolution rate of dolomite is relatively slow, it may still reach dissolution equilibrium at the basin scale. Shibing belongs to the subtropical monsoon climate zone. Rainfall and heat occur in the same period. Generally, high temperature and strong microbial activity in the rainy season promote dissolution and increase ion production, but at the same time, the increase of precipitation will dilute various ions. However, there is little difference between ion concentrations in rainy season and in dry season, and the ion concentration is relatively stable. That is to say, the large variation of flow rate does not cause the large variation of ion concentration, which indicates that the variation of flow rate is much larger than that of ion concentration. On the one hand, water-rock interaction may be the main factor affecting ion concentration; on the other hand, it indicates that there may be chemical steady-state behavior at the outlet of Huangzhouhe Basin.
From Figure 6, it can be seen that the pH of YYD spring water at point No. 29 ranges from 7.40 to 7.89, which is lower than that of the surface rivers. EC, Ca2+, Mg2+, and HCO3- of the spring water are lower in rainy season than in dry season, and the maximum values all occur in dry season, but their variation coefficients are less than 0.1, which is relatively stable in a hydrological year. Saturation indexes SIC and SID have the same trend of change, and the change is not significant, and each month is greater than 0, basically reaching saturation state. It is also evident that the groundwater of the dolomite reaches the dissolution equilibrium.
4.3. Factors Influencing the Formation of Hydrochemistry in Dolomite Basin
4.3.1. Geological Background and Karstification Control
(1) The Geochemical Effect of Carbonate Minerals on Main Ions. The dissolution reaction equation of carbonate dissolution limestone and dolomite can be expressed as follows:
For limestone:
For dolomite:
Limestone and dolomite dissolve to form Ca2+, Mg2+, and HCO3-. In theory, (Ca2++Mg2+) : HCO3- molar ratio is 1 : 2. The samples of karst surface water (Ca2++Mg2+)/HCO3- in the basin are located near the 2 : 1 line of carbonate mineral dissolution reaction. Some of the samples are slightly higher and some are slightly lower (Figure 7). The concentration of HCO3- varies from 2.90 to 5.20 meq/L. Most of the samples of the karst groundwater (Ca2++Mg2+)/HCO3- fall slightly below the 2 : 1 line, and their HCO3- concentration ranges from 1.90 to 6.20 meq/L. Karst surface water and groundwater have similar trends. Most of the water sample points are below the 2 : 1 line of carbonate rock, indicating that the surface water and groundwater in the basin are rich in Ca2+ and Mg2+, and carbonate rock (dolomite, limestone) is not the only source of Ca2+ and Mg2+ compared with HCO3-. The contribution of atmospheric precipitation to Ca2+ is 0.03 mmol/L, 0.003 mmol/L to Mg2+, and 0 to HCO3-. Even after precipitation correction, its contribution to surface water and groundwater is very weak, (Ca2++Mg2+) is still relatively enriched. The hydrochemical analysis in the basin shows that the concentration of SO42- is very low, and there is no sulfur-bearing stratum in the area, and sulfuric acid has little effect on the dissolution of carbonate rocks. The sample points of (Ca2++Mg2+)/HCO3- in the allogenic water tend to the cationic side, which indicates that the contribution of allogenic water to Ca2+ and Mg2+ is more than that to HCO3-.
Because of the existence of limestone and impure dolomite in the basin, Ca2+ comes from both dolomite mineral dissolution and limestone mineral dissolution, while Mg2+ only comes from dolomite mineral dissolution. To clarify the relative contribution of limestone and dolomite dissolution, the relationship between Mg2+/Ca2+ in basin water is shown in Figure 8. When only limestone dissolves, assuming that limestone has small amount of Mg to replace Ca, Mg2+/Ca2+ is much lower than 0.10; when the same amount of limestone and dolomite dissolves, ; the composition of dolomite increases against limestone, and then, Mg2+/Ca2+ approaches 1 [24]. Mg2+/Ca2+ in the karst surface water and groundwater is higher than 0.50. Most of the groundwater has , which is close to 1. A small part is between 0.50 and 0.75, with an average value of . The trend of Mg2+/Ca2+ in karst surface water samples is similar to that of the groundwater. Most of them are located on the line of , and their variation ranges from 0.64 to 1.30, with an average value of . The results show that Ca2+, Mg2+, and HCO3- in karst surface water and groundwater mainly come from dolomite mineral dissolution, and the contribution of limestone mineral dissolution is small. The Ca/Mg ratios are almost all around 1 in the samples of Shilengshui Formation and Loushanguan Formation in the basin, which fully indicates the control of dolomite on hydrochemistry. Mg2+/Ca2+ in karst surface water is higher than that in groundwater. The possible reason is that after the groundwater is exposed to become surface water, CO2 in the water degasses, resulting in CaCO3 deposition, which makes the concentration of Ca2+ in the water decrease relatively.
(2) The Geochemical Effect of Other Lithology on Main Ions. Lithology determines the hydrochemical composition of the basin and is the material basis for the formation of hydrochemistry. Figure 2 shows that the surface water and groundwater in the basin are HCO3-Ca·Mg, and the allogenic water is HCO3-Ca type. The karst water is mainly controlled by dolomite mineral dissolution, while the allogenic water is mainly controlled by shale. Comparing the hydrochemistry of the allogenic surface water (Tangtou River 22, Chayuan River 23) and outcropped spring water (Hetaoping 36) flowing through the northern shale area of the basin with that of surface water (Liangcha River East No. 10) flowing into the dolomite karst area (Figure 9), the main ion composition and saturation index are obviously different. Ca2+ is the main cation in the water flowing through the shale area, and Mg2+ content is very low. Therefore, the Mg2+/Ca2+ ratio in the allogenic water is the lowest (Figure 8), which is mainly related to the fact that the Mg content in rock is lower than that of Ca. Tests of noncarbonate rock samples from the Cambrian Bianmachong Formation in Hetaoping show that CaO and MgO account for 4.63% and 1.01%, respectively, while dolomite and calcite account for 4.62% and 5.75%, respectively. The content of HCO3- in anions is very low with an average of 0.49 mmol/L, and the concentration and proportion of Ca2+, Mg2+, and HCO3- in anions are obviously increased after they flow through dolomite area. Significant differences in major ions also lead to significant differences in conductivity. In terms of SIC and SID, both SIC and SID in water flowing through shale area are negative and have strong dissolution ability, while after flowing through dolomite area, SIC and SID reach 1.04 and 2.05, respectively, reaching supersaturation state.
(a)
(b)
(c)
(d)
4.4. Influence of Atmospheric Precipitation
Surface water and groundwater in Shibing Shanmuhe Basin mainly come from atmospheric precipitation, but due to the relatively low ion composition in precipitation, its main influence is dilution. However, the low pH of atmospheric precipitation in the study area can lead to the decrease of the pH of karst water, thus enhancing the carbonate karst dissolution and increasing the conductivity of the basin. According to the Gibbs diagram above, the hydrochemistry of the basin belongs to rock weathering type. Therefore, the contribution of atmospheric precipitation to the chemical ion composition of basin water is very restricted.
4.5. Influence of Human Activities
The concentrations of K+, Na+, NO3-, Cl-, and SO42- in the water in the study area are very low. Because the study area is a world natural heritage site, there are no serious human pollutions such as industrial activities, but rather only a small area of traditional agricultural activities whose impact is very small. The typical dolomite karst ecosystem has been well preserved in the study area, and a large area of primary forest vegetation and diversified ecosystem has been preserved and developed. Primitive central subtropical evergreen broad-leaved forest ecosystem has been formed in the middle and lower parts of the mountain. Mountainous evergreen hard-leaved forest and coniferous forest have been formed at the top of the ridge and mountain. They belong to the top community type of recessive karst forest vegetation, including woodland coniferous forest shrubs and sparse woodlands that account for more than 87% of the total area of the study area. Therefore, human activities have little effect on the hydrochemistry of the basin.
5. Discussions
(1)The pH of surface water and groundwater in typical dolomite basin in humid subtropical zone is neutral to alkaline, controlled mainly by geological background, and the hydrochemistry of surface water and groundwater is quite different. However, Ca2+, Mg2+, and HCO3- are the main ions in surface water and groundwater in the dolomite karst basin. The hydrochemical type in typical dolomite basin in humid subtropical zone is HCO3-Ca·Mg type, which is the same as the previous study by Zhang [20]. It is different from the hydrochemical types in limestone and shale districts, in which the proportion of Cl-, NO3-, and SO42- is higher, and the hydrochemical type is HCO3-Ca. The hydrochemical type in limestone area is also HCO3-Ca type(2)Controlled by lithology, EC, Ca2+, Mg2+, HCO3-, Mg2+/Ca2+, and Mg2+/HCO3- have similar spatial variations, which are lower in allogenic water in the upstream and northwest of the basin and become higher when the water enters the dolomite karst area in the mid and downstream of the basin. HCO3- concentration and Mg2+/Ca2+ in the pure dolomite basin are higher than those in the impure dolomite basin. Although the dissolution rate of dolomite is slower than that of limestone, its higher solubility results in more HCO3-. SIC and SID of the allogenic water are below zero, which strongly erodes the carbonate rocks in the basin. As a result, Ca2+, Mg2+, and HCO3- of the allogenic water entering the karst area are significantly increased(3)The monthly variation of hydrochemistry shows that the seasonal variation of ion concentration is not obvious, and there is a dilution effect in rainy season, but the variation of ion concentration is small and relatively stable, indicating that water-rock interaction may be the main factor affecting ion concentration(4)TDS and Na+/(Na++Ca2+), Cl-/(Cl-+HCO3-) indicate that surface water and groundwater in the basin are mainly controlled by rock weathering, and rainfall has little effect. (Ca2++Mg2+)/HCO3- samples locate near the 2 : 1 line of mineral dissolution reaction of carbonate rocks, which indicates that the formation and variation of hydrochemistry in surface water and groundwater in karst basin are mainly controlled by carbonate weathering. Atmospheric precipitation and allogenic water have some influence on the input of Ca2+ and Mg2+ in karst water, which makes the ion balance slightly biased towards cations. Low concentrations of K+, Na+, NO3-, Cl-, and SO42- in water indicate that human activities have little effect on the hydrochemistry in the study area
6. Conclusions
The hydrochemical characteristics and variation rules of water in the humid subtropical karst areas in the Shanmuhe Basin and its subbasin Huangzhouhe Basin in Shibing Dolomite Karst World Natural Heritage Site in China were studied. 215 water samples were collected and analyzed; the geochemical data were elaborated by using graphical methods and calculation of saturation indexes. It was found that the hydrochemical type in typical dolomite basin in humid subtropical zone is HCO3-Ca·Mg type, and the water-rock interaction may be the main factor affecting the hydrochemical characteristics through carbonate weathering. The lithology is the main influencing factors of the hydrochemical features of the study area which is less affected by human activities.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
SX was in charge of designing the study topic and writing the manuscript. JW, CZ, and JL were in charge of checking the research results and paper structures. YD and JH performed the experiments, and HX was in charge of mapping.
Acknowledgments
This work was financially supported by the Project of National Key Research and Development Program of China in the 13th Five-Year Plan (2016YFC0502606-01), the National Natural Science Foundation of China (41673129), the Sci & Tech Project of Guizhou Province, China (Qian-ke-he LH [2015]7775), the Science and Technology Cooperation Project of Guizhou Province (No. LH [2017]7059), and the doctoral research project supported by Guizhou Normal University (GZNUD [2017]12).