General Hydrochemistry
According to the detection results, Seven conventional ions( Ca2+, Mg2+, Na+, K+, Cl-, SO42-, HCO3− ) concentration,pH and TDS contour maps were drawn respectively(Fig. 2 ). The runoff-recharging directions of Ordovician limestone water can be divided into two types: one is from the Zhuozishan Mountain to Gongdeer Mountain, and then to south-north runoff, and the other is the runoff from the Zhuozishan Mountain to the south. Obviously, the distribution characteristics of hydrochemical indicators are related with the direction of the groundwater runoff. The concentration of Na++ K+, HCO3-, Cl-, SO42-, TDS, and pH generally increase from the recharge area to the discharge area, whereas Ca2+ and Mg2+ showed the opposite distribution characteristic. Especially in the central region, this distribution is more conspicuous. The contour distribution of various indicators are relatively dense in the northern and southern parts, indicating that the hydrochemical indicators varied greatly.
Hydrochemical Types
Piper trilinear diagram is favourable for determining connections of different dissolved components and types of groundwater according to its hydrochemical characters(Tiwari et al.2017, Shakya et al.2019). As indicated by the Piper diagram(Fig. 3). The hydrochemical types of Ordovician limestone water in the Zhuozishan Coalfield mainly included SO4▪Cl-Ca▪Na, HCO3▪SO4▪Cl-Ca▪Na, HCO3▪Cl-Ca▪Na, and HCO3-Ca▪Na type. According to the direction of groundwater runoff, the changes in hydrochemical types can be divided into two types, as shown in the figure: ①from the Zhuozishan Mountain to the Gongdeer Mountain and then to the north-south directional runoff, with the hydrochemical type changing from SO4▪Cl-Ca▪Na type in the upstream region to HCO3▪SO4▪Cl-Ca▪Na, HCO3▪Cl-Ca▪Na, or HCO3-Ca▪Na; ②the runoff from the Zhuozishan Mountain to the south, and the hydrochemical type is always the SO4▪Cl-Ca▪Na type. According to TDS contour map, the TDS concentration of runoff from the Zhuozishan Mountain to the Gongdeer Mountain and then to the north is 1,200–1,400 mg/L, and that of the runoff to the south ranges from 600–1,000 mg/L. The TDS concentration of the runoff from the Zhuozishan Mountain to the south is 1,200–2,800 mg/L, with a gradual increase. Obviously, the proportion of HCO3 (mg equivalent) increases, whereas the TDS concentration changes only slightly in the runoff process from the Zhuozishan Mountain to the Gongdeer Mountain and then to the south or north. In the runoff process from the Zhuozishan Mountain to the south, the hydrochemical type remains unchanged but the TDS concentration gradually increases, indicating that hydrochemical reactions occur in the groundwater during the runoff process.
Hydrochemical Process
1. Rock weathering action
The Gibbs diagram (Gibbs 1970) was an important means primitively to study the mechanisms surface water evolution, and now it is widely used in groundwater researches (Li et al. 2016a).Generally, the main control factors can be divided into three types: evaporation and concentration dominance, rock weathering dominance, and precipitation dominance(Pehlivan et al. 2020). The analytical results show that groundwater hydrochemical process includes two types:① the runoff from Zhuozishan Mountain and then to north-south directional runoff, and groundwater control action is mainly rock weathering action ② the runoff from Zhuozishan Mountain to the south, and the control action of groundwater transforms from the rock dominance to the evaporation dominance from the upstream to downstream. However, the Ordovician limestone groundwater is deeply buried, which will not be subjected to evaporation and concentration. Therefore, it can be speculated that the reverse cation exchange effect may occur in the runoff direction from Zhuozishan Mountain to the south, resulting in the increase of some ions concentration.
2. Cation Exchange
Cation exchange is also a vital process affecting groundwater chemical feature (Zhang et al.2020; Shi et al. 2017). Typically, two chloro-alkaline indices (CAI-1 and CAI-2) proposed by Schoeller (1965) were important way to research the occurrence of cation exchange, where the units of all ions are meq/L. When CAI-1 and CAI-2 are negative, it indicates that an ion exchange effect occurs, in which Na+ and (or) K+ in water-bearing media are replaced by Ca2+ and (or) Mg2+ in groundwater. Whereas CAI-1 and CAI-2 are positive, the reverse cation exchange effect occurs (Qian et al. 2016). The calculation equations of CAI-1 and CAI-2 can be estimated as follows (Li et al. 2014):
In this study, north and central parts of Zhuozishan Coalfield groundwater samples are drawn in the lower left area, whereas south part samples are mostly in the upper right area of Fig. 5a. This suggests that an ion exchange effect(Eq. (3)) occurs in the runoff from Zhuozishan Mountain to Gongdeer Mountain and then to south-north directional and the reverse cation exchange(Eq. (4)) occurs the runoff from Zhuozishan Mountain to the south. In other words, groundwater cation exchanges of these two parts occur reverse.
To verify the hypothesis of cation exchange reactions, the relationship between (Na+ + K+− Cl−) and [(Ca2+ + Mg2+) – (HCO3−- SO42−)] was examined(Fig. 5b). Na+and Cl− come from halite dissolution, so (Na+ + K+ − Cl−) values shows the amount of Na+ gains or lose to that caused by halite dissolution (Li et al. 2018b). And, [(Ca2+ + Mg2+) – (HCO3−- SO42−)] values shows the amount of Ca2+ and Mg2+ gains or lose to that caused by gypsum, calcite and dolomite dissolution (Farid et al. 2013). Figure 5a illustrates that all groundwater samples are distributed along or closed to the line with − 1 slope indicating that the Ordovician limestone groundwater is all affected by cation exchange. The samples from coal mines located in the northern and central parts of Zhuozishan Coalfield are mainly distributed in Ⅳ quadrant of the coordinate axis, whereas those in the south part are mainly distributed in Ⅱ quadrant of the coordinate axis. Besides, the farther it is from Zhuozishan Mountain (groundwater recharging area), the farther the water sampling point is from the origin of coordinates, which indicates that the farther it is from the recharging area, the stronger the cation exchange effect.
3. Main ion sources
Na++K+ vs Cl− can be used to study the sources of Na+ and K+ in groundwater. When the samples plots near the 1:1 equiline, it shows that Na+ and K+ mainly originate from halite dissolution (Li et al. 2016b). Figure 6a indicates that groundwater samples from the south and north of Zhuozishan Coalfield plot far from the 1:1 line, indicating that Na+ and K+ originate from halite dissolution, barely be affected by cation exchange.
Previous researches have shown that if the sulphate rocks dissolution are the main processes in a groundwater system, the (Ca2++Mg2+)/ SO42− ratio will be or closed to 1(Li et al. 2018c). Figure 6b shows that most of the water samples are distributed near 2:1 above the 1:1 line, indicating that Ca2+ and Mg2+ in the Ordovician limestone groundwater scarcely comes from sulfate minerals, and they mainly comes from carbonate or silicate dissolution
Besides, if the carbonate dissolution is the main processes in a groundwater system, the (Ca2++Mg2+)/ HCO3− ratio will be or closed to 1. Figure 6c indicates that only few groundwater samples plot near the 1:1 line. The south and north of the coalfield water samples plot above the 1:1 aquiline, Whereas samples from the middle part of the coalfield are mainly distributed below the 1:1 line. Combined with Fig. 4b and 4c, Ca2+ and Mg2+ are mainly from the silicate dissolution, with a minority from the other kinds of rocks dissolution.
Mechanisms of Hydeogeochemistry
According to previous sections, the Ordovician groundwater runoff in Zhuozishan Coalfield can be divided into two directions: one is from Zhuozishan Mountain to Gongdeer Mountain, and then to the north-south directional runoff; the other is from Zhuozishan Mountain to the south runoff. And according to the hydrochemical analysis above, we divided Zhuozishan Coalfield into three characteristic regions including “northern”, “central” and “southern” parts. Hydrochemical mechanism of each region is summarized, respectively as follows:
1. Northern part
In the northern part, groundwater flows westward after being recharged from Zhuozishan Mountain. It passes through the Kabuqi syncline to the eastern margin of Gongdeer Mountain, which is blocked by water-blocking fracture, and changes to north-south directional runoff in the syncline. Due to the sudden change in the direction of groundwater runoff, the retention area of groundwater is formed in the Kabuqi syncline. Besides, due to the constant recharge of groundwater, the direction of groundwater runoff in the retention area is disordered (Fig. 7). Thus, groundwater flows to the downstream area for a long time, resulting in serried and non-uniform contours of hydrochemical indices (Fig. 2), and the hydrochemical types and the hydrochemical process are more complicated in turn.
2. Central part
In the central part, groundwater comes from the runoff in the north or Zhuozishan Mountain recharge area, and then it flows downstream to the southern part of Zhuozishan Coalfield along the fault zone. Therefore, the groundwater in the central part has not been blocked by geological structure, and the runoff direction remains unchanged. Thus the contours of hydrochemical indices (Fig. 2) regularly change, and the hydrochemical type is dominated by HCO3▪Cl-Ca▪Na(Fig. 3). The hydrochemical process is controlled by rock weathering(Fig. 3) action and cation exchange ( Fig. 6 )
3. Southern part
In the southern part, groundwater comes from the runoff in the north or Zhuozishan Mountain recharge area. However, due to the gradual deepening of strata in the south and the blocking of Zhengyiguan strike-slip fault, a retention zone has been formed (Fig. 7). Then, it runoff to the southwest Therefore, affected by the above two factors, the contours of hydrochemical indices is dense and regular(. Figure 2).The hydrochemical type is always dominated by SO4▪Cl-Ca▪Na type (Fig. 3), with a stronger reverse cation exchange effect (Fig. 6).