Analysis of the Erosion Law of Karst Groundwater Using Hydrogeochemical Theory in Liulin Spring Area , North China

The comprehensive geological, hydrogeological and hydrogeochemical model of the Liulin karstic spring area in the eastern limb of the Ordos syncline was established by a combination of chemical thermodynamics, chemical kinetics and hydrogeology. The study area was divided into four zones based on the saturation indices of calcite, dolomite and gypsum, which were computed by the groundwater-chemical simulation software PHREEQC (a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations), with consideration of the geological and hydrogeological conditions and hydro-geochemical reactions. The weight and volume modulus of carbonate rocks and sulphate rocks in each zone were calculated by the method of correlation analysis to evaluate the dissolution law of karst groundwater. The results showed that in the zone I (saturation index of calcite βc ≤ 1) the dissolution of calcite was the major geochemical reaction, the weight modulus of calcite was higher than that of dolomite and gypsum, and the pore space generated by the dissolution of calcite was one order of magnitude larger than that of dolomite and gypsum. In zone II (saturation index of calcite βc > 1 to saturation index of dolomite βd ≤ 1) the corrosion moduli were all smaller than that in zone I, the solubility of dolomite and gypsum increased, and calcite reached saturation. The space occupied by the calcite sediment was less than that dissolved by dolomite and gypsum. In zone III (saturation index of dolomite βd > 1 to saturation index of gypsum βg ≤ 1), calcite and dolomite had reached saturation, accompanied by dedolomitization, and the amount of gypsum dissolution increased obviously. The conclusions indicate that the hydrogeochemical environment plays an important role in mineral dissolution.


Introduction
Carbonate rocks are widely distributed in China, and its acreage accounts for about a third of the land area.In northern China, the location of karst spring is not only an important resource for local water supply but is also a famous tourist site [1].Hence, the government sector and local residents have been attaching importance to water quality problems for many years.Water quality is closely associated with hydrochemical compositions [2], and carbonate rock dissolution plays an important role in the variation of such hydrochemical compositions.Hence, the study of the law of carbonate rock dissolution is of significant importance for obtaining a more comprehensive understanding of the groundwater environment.Scholars at home and abroad usually relate to qualitative and quantitative analysis of karstification and karst dissolution when they analyze the hydrogeology, engineering geology, and mineral geology problems in karst spring area [3].It is necessary for analysis of karst hydrogeology characteristics to make a quantitative analysis and calculation of karst dissolution.For a long time, problems about magnitude and distribution of carbonate rocks karst dissolution-especially in the field of karst hydrogeochemistry and groundwater mixing-have been discussed by many scholars [4][5][6][7][8].
Currently, the methods for determining carbonate rock dissolution quantities mainly include the following.Firstly, to measure the karstic dissolution quantities repeatedly on the surface of the karst with micro erosion measurement instruments in the research field and then to calculate the actual amount of dissolution quantities [9,10].Mottershead [9] had measured the eroded quantities of weathered bedrock for seven years and computed the average decrease rate (0.625 mm/a).Secondly, to put the standard test pieces of carbonate rocks into the air, ground, and soil layer in the study filed or the indoor specific experimental facility [11], then to measure the weight of the test pieces after a period of time and calculate the dissolved quantities during this process based on the quantity difference of test pieces before and after the carbonate rocks dissolution.Plan [12] investigated the influences of altitude, subsoil vs. sub-aerial exposure, vegetation, karst morphology, soil humidity, sample lithology, and sample surface morphology to carbonate dissolution by measuring the mass difference of carbonate tablets exposed to natural dissolution for one year.Thirdly, to calculate the total dissolved quantities (X, viz.X = T.Q) of the catchment basin in the karst area by measuring the outlet flow (Q) of the spring or underground river mouth and the concentrations (T) of mineral that carried by water [13][14][15].Utilizing the formula of Corbel-which was one of the methods-Li et al. [16] computed the dissolved quantity of karst and analyzed the relationship of several factors such as temperature, annual rainfall, lithological characteristics and geomorphologic conditions based on a set of data of hydrochemistry points and the corresponding runoff depth of these points in Hongshui river basin.
In the 1990s, Cao et al. [17] proposed a whole set of theories and methods to determine the dissolved quantity of karst.However, so far these methods had not been applied widely in practical problems.The dissolution law of underground karst water in the Liulin spring area are calculated and analyzed by Cao's theory and method in this paper.It will provide some basis for the formation and evolution of karst groundwater characteristics.
The remainder of the paper is organized as follows.Section 2 describes the information about Liulin spring area, four zones which divided by saturation indices are discussed in Section 3. Section 4 presents the main results and discussions of dissolution law in three zones.The summary and conclusions are given in Section 5.

Physical Geography
Liulin spring is one of the famous karst springs in Shanxi Province, China, and expose dispersedly in Sanchuan river valley in the village of Xuejiawan-Zhaidong, 3 km from Liulin county (showed in Figure 1).The study area is west of the middle of Lvliang Mountain, there are high and middle mountains where bedrock is exposed, and low and middle mountains with loess hilly regions in the east and west of the research area respectively.The climate in the area is continental semi-arid with the characteristics of a windy spring, hot summer, intensive precipitation in fall and drought and cold winter.The average value of temperature and precipitation are 9.2 • C and 506 mm respectively.The rivers in the area belong to the Yellow River system.Beichuan, Dongchuan and Nanchuan Rivers collect as the Sanchuan River at the Jiaokou County, then flows into the Yellow River.The basin area of the Sanchuan River is 4161 km 2 , and its average annual runoff is 194 million m 3 /a.

Geology and Hydrogeology
Regular Precambrian metamorphic rock outcroppings can be found in the northern, eastern, and southern parts of the spring area, while Cambrian-Ordovician carbonate strata outcroppings are found mostly in the south-eastern and central parts of the spring area.Carboniferous-Triassic clastic rocks are mainly outcropped in the western and central regions of the spring area.Tertiary and Quaternary rickle overlies various ages of the bedrock.From a geotectonic point of view, the spring area is situated in a uniclinal structure which tilts to the west with a 2-to 8-degree dip angle, within which secondary folds and fractures develop.The geologic structure dominates the directions of the groundwater flow.
The northern and eastern boundaries of spring area are composed of drainage divides in metamorphic rock mountains, whereas the southern and south-eastern boundaries consist of subsurface divides in exposed areas of carbonate rocks.In the western part of the study area, carbonate rocks incline toward the syncline core of the Ordos Basin and are covered by thickbedded clastic rocks for hundreds of meters.There is no discharge path for groundwater.Hence the

Geology and Hydrogeology
Regular Precambrian metamorphic rock outcroppings can be found in the northern, eastern, and southern parts of the spring area, while Cambrian-Ordovician carbonate strata outcroppings are found mostly in the south-eastern and central parts of the spring area.Carboniferous-Triassic clastic rocks are mainly outcropped in the western and central regions of the spring area.Tertiary and Quaternary rickle overlies various ages of the bedrock.From a geotectonic point of view, the spring area is situated in a uniclinal structure which tilts to the west with a 2-to 8-degree dip angle, within which secondary folds and fractures develop.The geologic structure dominates the directions of the groundwater flow.
The northern and eastern boundaries of spring area are composed of drainage divides in metamorphic rock mountains, whereas the southern and south-eastern boundaries consist of subsurface divides in exposed areas of carbonate rocks.In the western part of the study area, carbonate rocks incline toward the syncline core of the Ordos Basin and are covered by thick-bedded clastic rocks for hundreds of meters.There is no discharge path for groundwater.Hence the groundwater is considered stagnant in this region.This boundary is defined as the lines where the roof depth of carbonate rocks reaches 1000 m.Cambrian and Ordovician marine carbonate rocks constitute the main aquifer for the water supply, especially the middle series of the Ordovician, which consists of limestone, dolomite, and some of gypsum.The upper and the lower confining beds are composed of Carboniferous-Triassic clastic rocks and Precambrian metamorphic rocks, respectively.The karst groundwater system is mainly recharged in exposed areas of carbonate rocks, as defined by the recharge area.Leakage of river water from the Sanchuan River and its tributaries in exposed areas of carbonate rocks serve as additional recharge sources.In the center of the spring area, carbonate rocks are covered by thick-bedded clastic rocks which are composed of an interbed of sandstone and mudstone, preventing the precipitation recharging the aquifer.This serves as the flow-through area, with the waters flowing from the north, east, and south parts of study area towards the Liulin spring.There are two major discharge ways in this groundwater system: spring discharge and well exploitation.Places where the springs emerge in the Sanchuan River valley constitute the dominant discharge regions of the groundwater system.

Sampling and Analysis
Sample collection and component analysis of groundwater is the basis of the following analysis, the precision of which may directly influence reliability and rationality of the results.The groundwater samples were collected along predominant groundwater flow paths.The testing of water chemical components including temperature, pH value, Eh value, HCO 3 − , SO 4 2− , Cl, Na + + K + , Ca 2+ , Mg 2+   and TDS (Total Dissolved Solids) was accomplished in the water environmental monitoring center of Shanxi province.Temperature, pH value and ORP (Oxidation-Reduction Potential) value were measured in situ using potable REX PHB-4 pH meters and a 501 ORP pole.Water samples for cation and anion analysis were filtered through 0.45 µm membranes in the field and collected in 500 mL polyethylene bottles.HCO 3 − was measured using the Gran titration method in the sampling day.Samples for cation analysis were acidified using 1:1 nitric acid to pH < 2 in field.The cations (K + + Na + , Ca 2+ , Mg 2+ ) and anions (Cl − , SO 4 2− ) were determined by atomic absorption spectrometry (PERSEE TAS990F) and ion chromatography (CIC300), respectively.The uncertainties of measurements were all within 5%.

Hydrogeochemical Partition and the Construction of Its Model
The minerals which formed the karst spring aquifer system mainly are calcite, dolomite, gypsum, etc.By studying the nature and characteristics of these minerals' dissolution and precipitation, the corrosion properties of carbonate rocks and sulphate rocks can be understood very well [18,19].The hydrogeochemical reaction proceeded in the system of solid and liquid phases that consisted of different mineral compositions, and the chemical reactions in the hydrogeochemical reaction would change along the direction of groundwater flow [20].
The saturation index (β) can not only determine the direction of dissolution and precipitation of different minerals but is also an important index for dividing the spring area into different hydrogeological zones and different hydrogeochemical environments.In the study area, the features of geology, hydrogeological conditions, hydrogeochemical reactions in the spring area, and the saturation indices of water points in the research area-which are calculated by the water chemical equilibrium simulation software PHREEQC-are comprehensively considered for dividing the four zones shown in Figure 1 (viz.βc ≤ 1, βc > 1 ~βd ≤ 1, βd > 1 ~βg ≤ 1 and βg > 1, where, βc, βd and βg represent the saturation index of calcite, dolomite and gypsum respectively).The hydrogeological-hydrogeochemical model for Liulin spring, as shown in Figure 2, is set up based on these horizontal zones.

βc ≤ 1 (Zone I)
This zone, where the limestone in the layer of Cambrian and Ordovician is exposed or halfexposed and the surface karst developed well, is mainly distributed in the east, north and south of the study area.The groundwater form in zone I is phreatic water with the characteristics of vertical seepage, atmospheric precipitation recharge for groundwater, and the total dissolved solids (TDS) equals 260 mg/L.Calcite, dolomite and gypsum are in the dissolved state.The phreatic water surface and the line where βc equals 1 are the upper bound and lower bound of zone I respectively.Since the upper bound (phreatic water surface) is in contact with the air in the soil directly, this zone is called the CO2 open system.The karst water chemistry reactions, including the reaction of CO2 dissolved in phreatic water and the reaction between H + generated by CO2 dissolution and the calcite and dolomite in carbonate rocks, nevertheless, have little effect on the process of gypsum dissolution.The main equations of chemical reaction models in zone I are as follows: CO2 + H2O = H2CO3

βc ≤ 1 (Zone I)
This zone, where the limestone in the layer of Cambrian and Ordovician is exposed or half-exposed and the surface karst developed well, is mainly distributed in the east, north and south of the study area.The groundwater form in zone I is phreatic water with the characteristics of vertical seepage, atmospheric precipitation recharge for groundwater, and the total dissolved solids (TDS) equals 260 mg/L.Calcite, dolomite and gypsum are in the dissolved state.The phreatic water surface and the line where βc equals 1 are the upper bound and lower bound of zone I respectively.Since the upper bound (phreatic water surface) is in contact with the air in the soil directly, this zone is called the CO 2 open system.The karst water chemistry reactions, including the reaction of CO 2 dissolved in phreatic water and the reaction between H + generated by CO 2 dissolution and the calcite and dolomite in carbonate rocks, nevertheless, have little effect on the process of gypsum dissolution.The main equations of chemical reaction models in zone I are as follows: Water 2018, 10, 674 This zone is distributed in the Sanchuan River and the valley of its main tributaries.The main aquifer is covered by quaternary system strata.The velocity of karst water circulation is relatively quick because of the shallow buried layer.In this zone, the direction of groundwater flow transformed from vertical to horizontal, the value of TDS is 334 mg/L and calcite is in the saturated state, while dolomite and gypsum are in the dissolved state.The line where βc equals 1 and the line where βd equals 1 are respectively the upper and lower bounds of zone II.The form of groundwater transforms from phreatic water to confined water which CO 2 cannot directly pass into, so the hydrogeochemical environment of this zone is a relatively closed system.The karst water chemistry reaction mainly consumes CO 2 that dissolved in zone I, including the precipitation of calcite and dissolution of dolomite and gypsum.The main equations of chemical reaction models of zone II are as follows: The zone III, whose aquifer is covered by the Carboniferous, Permian and Triassic strata, is mainly distributed in the midwest and northwest of the Liulin spring area.It should be noted that the aquifer in zone III has the following characteristics: TDS equals 846 mg/L, there is little variation in the range of water level, as well as slower groundwater circulation, flow velocity and permeability.Mineral ions in carbonate rocks have all reached the supersaturated state then precipitated from the water, leaving sulphate minerals to continue dissolving.The upper and lower bound of zone III are respectively the line where βd equals 1 and the line where βg equals 1.In addition to the precipitation of carbonate rocks and dissolution of sulphate rocks, the karst chemistry reaction is also accompanied by the reactions of de-sulphidation and de-dolomization.The main equations of chemical reaction models in zone III are as follows: This region, located in the west of the Liulin spring area, not only has a relative good sealing performance but also is a reduction zone.The permeability of the aquifer in this area becomes worse and the groundwater circulation becomes very slow; zone IV is a stagnant area, in which the mineral ions of carbonate rocks and sulphate rocks are in the supersaturated state.The main equations of the chemical reaction models in zone IV are as follows: Water 2018, 10, 674 7 of 15

Calculation and Discussion of Corrosion Modulus
The corrosion modulus which consists of the weight modulus (M) and volume modulus (N) in each karst zones mentioned above are calculated by the method of correlation analysis with combining the simplified model of geological-hydrological geochemistry and the theory of the chemical thermodynamics and chemical kinetics.The weight modulus is defined as the weight of minerals dissolved by a unit volume of water, while the volume modulus is defined as the void volume dissolved by a unit volume of groundwater flowing through a certain karst area, or is filled by certain ions' precipitation.In the process of calculating for M and N, the dissolution content of some ions in certain zones is computed by deducting the average ionic content in the rainfall from the total ionic content in groundwater.The calculated results are shown in Table 1.It is necessary to mention that three zones (I, II and III) are analyzed in this section, since the mineral dissolution of carbonate rocks and sulphate rocks mainly occurs in these three zones, and all of the ions in zone IV are in the supersaturated state without dissolution, which means that corrosion modulus cannot be calculated.According to the antecedent model, three minerals (calcite, dolomite, gypsum) can be totally dissolved in zone I. Based on the material equivalent principle of chemical combination and decomposition, the value of Ca 2+ computed by water quality analysis is supposed to be the total quantities of Ca 2+ (T Ca2+ ) dissolved from calcite, dolomite and gypsum.The amount of Ca 2+ and SO 4 2− dissolved from gypsum is equivalent, so that this part of Ca 2+ can be expressed as (T Ca2+/SO42− ).Based on the previous analysis, the content of Ca 2+ dissolved from carbonate rocks, which has correlativity with the content of HCO 3 − , can be computed by deducting T Ca2+/SO42− from T Ca2+ .
Water 2018, 10, 674 8 of 15 Assuming that the contents of gypsum, calcite and dolomite dissolved into groundwater are m 1 , m 2 and m 3 respectively, the dissolution equations of three kinds of minerals are presented as follows: Combining Equations ( 1)-( 3) and balancing the coefficient, the following equation is obtained: Equation ( 4) reveals that the product dissolved by sulphate has no effect on the dissolution of carbonate rocks, so the content of ion dissolved from sulphate and carbonate rocks could be calculated independently.In this zone, little ion is accumulated in the runoff flow and the human and other natural factors have a smaller influence on the aquifer.Therefore, the hypothesis that the dissolved quantities of Ca 2+ and SO 4 2− are identical is reasonable, namely: where, y 1 denotes the average content of SO 4 2− (mmol/L), X 1 represents the average content of T Ca2+/SO42-(mmol/L), and A 1 and B 1 are a constant that usually equal 1 and 0, respectively.Based on material equivalent principle of decomposition, a linear relationship exists between the quantities of HCO 3 − , Ca 2+ and Mg 2+ , so the correlation equations are established as follows: y 3 = A 3 X 2 + B 3 (7) where, y 2 and y 3 are the contents of Ca 2+ dissolved from calcite and Mg 2+ dissolved from dolomite (mmol/L), X 2 denotes the average content of HCO 3 − (mmol/L), A 2 and A 3 are the regression coefficients, B 2 and B 3 are constant.The values of dissolved minerals are utilized to draw the relationship plots of Ca 2+ with SO 4 2− , Ca 2+ (Ca 2+ − SO 4 2− ) with HCO 3 − and Mg with HCO 3 − (shown in Figure 3).The following correlation equations can be observed in Figure 3. Combining Equations ( 1)-( 3) and balancing the coefficient, the following equation is obtained: Equation ( 4) reveals that the product dissolved by sulphate has no effect on the dissolution of carbonate rocks, so the content of ion dissolved from sulphate and carbonate rocks could be calculated independently.In this zone, little ion is accumulated in the runoff flow and the human and other natural factors have a smaller influence on the aquifer.Therefore, the hypothesis that the dissolved quantities of Ca 2+ and SO4 2− are identical is reasonable, namely: y1 = A1X1 + B1 (5) where, y1 denotes the average content of SO4 2− (mmol/L), X1 represents the average content of TCa2+/SO42-(mmol/L), and A1 and B1 are a constant that usually equal 1 and 0, respectively.Based on material equivalent principle of decomposition, a linear relationship exists between the quantities of HCO3 − , Ca 2+ and Mg 2+ , so the correlation equations are established as follows: where, y2 and y3 are the contents of Ca 2+ dissolved from calcite and Mg 2+ dissolved from dolomite (mmol/L), X2 denotes the average content of HCO3 − (mmol/L), A2 and A3 are the regression coefficients, B2 and B3 are constant.The values of dissolved minerals are utilized to draw the relationship plots of Ca 2+ with SO4 2− , Ca 2+ (Ca 2+ − SO4 2− ) with HCO3 − and Mg with HCO3 − (shown in Figure 3).The following correlation equations can be observed in Figure 3. y2 = 0.4012X2 − 0.1921 y3 = 0.0603X2 + 0.3328 To substitute A1X1 for m1, A2X2 for m2, A3X2 for m3, Equation ( 4) changes to: To substitute A 1 X 1 for m 1 , A 2 X 2 for m 2 , A 3 X 2 for m 3 , Equation (4) changes to: where, ∑C ≤ B 2 Ca 2+ + B 3 Mg 2+ , ∑C can be seen as the constant in Equations ( 6) and (7), and is also the product of carbonate rock dissolution and has nothing to do with the dissolution of sulphate rocks.
The following comprehensive chemical equation is established in the form of Equation ( 4): where, B 2 Ca 2+ represents the total content of Ca 2+ dissolved from calcite and dolomite, B 3 represents the content of Ca 2+ dissolved from the dolomite, and (B 2 − B 3 ) represents the content of Ca 2+ dissolved from calcite.The Equations ( 8) and ( 9) are combined to obtain the following equations: The Equations ( 10)-( 12) are employed to compute the weight modulus of gypsum, dolomite and calcite (M g , M d , M c ): M g = 0.5A 1 X 1 CaSO 4 •2H 2 O = 86.09X 1 = 86.09× 0.15 = 12.914(g/m 3 ) M d = 0.25(A 3 X 2 + B 3 )CaMg(CO 3 ) 2 = 46.1(A 3 X 2 + B 3 ) = 46.1 × (0.0603 × 3.65 + 0.3328) = 25.489(g/m 3 ) The total weight modulus (M 1 ) is the sum of each mineral modulus.Namely: The volume modulus (Ng, N d , N c ) of gypsum, dolomite and calcite are computed as follows: Then the total volume modulus (N 1 ) is computed as: In the exposed area of carbonate rocks where βc ≤ 1, since the atmospheric precipitation has directly contact with aquifer the velocity of water exchange becomes rapid, which intensify the extent of water dissolution of the karst .Calcite, dolomite and gypsum are in the dissolved state, where the dissolved quantities of gypsum in the unit volume of water is minimum; M g equals 12.914 g/m 3 and N g equals 5.615 × 10 −6 , and the dissolved quantities of dolomite are such that M d equals 25.489 g/m 3  and N d equals 9.103 × 10 −6 .The dissolved quantities of calcite in the unit volume of water is maximal with M c equal to 47.016 g/m 3 and N c equals 1.741 × 10 −5 .The total weight modulus and volume modulus of the three minerals are 85.418 g/m 3 and 3.213 × 10 −5 , respectively.Evidently, calcite dissolution, whose dissolved quantity is larger than that of dolomite and gypsum, is more likely given priority in this zone.The mineral dissolution lead to the new space increasing, and the analysis which is mentioned above indicates that the volume modulus of new space created by calcite dissolution is one order of magnitude larger than that of dolomite and gypsum dissolution.y 5 = A 5 X 5 + B 5 (17) where, y 5 and X 5 denote the content of SO 4 2− and Ca 2+ , which are equivalent (mmol/L).A 5 and B 5 are constants of 1 and 0, respectively.
According to the values of dissolved minerals, the relationship plots of Ca 2+ with SO 4 2− and Mg with HCO 3 − are drawn (Figure 4).The following correlation equations can be observed from Figure 4.
where, y5 and X5 denote the content of SO4 2− and Ca 2+ , which are equivalent (mmol/L).A5 and B5 are constants of 1 and 0, respectively.
According to the values of dissolved minerals, the relationship plots of Ca 2+ with SO4 2− and Mg with HCO3 − are drawn (Figure 4).The following correlation equations can be observed from Figure 4. y4 = 0.7429X4 − 2.0786 The following equation can be obtained by putting the Equations ( 16) and ( 17) into Equation (15): where ∑C ≤ B4Ca 2+ + B4Mg 2+ is the product of carbonate rock dissolution and has nothing to do with the dissolution of sulphate: The following equation can be obtained by putting the Equations ( 16) and ( 17) into Equation (15): Equations ( 18) and ( 19) are combined to obtain the following equations: The total volume modulus (N 2 ) can be calculated as: In zone II βc >1 > βd ≤ 1, the dissolved quantity of gypsum in the unit volume of water is still minimal, M g equals 21.523 g/m 3 and N g equals 9.358 × 10 −6 .The dissolved quantity of dolomite is increased, which can be seen from M d (36.030 g/m 3 ) and N d (1.278 × 10 −5 ).The dissolution of calcite begins to be subside since it has already reached the supersaturated state and the values of M c and N c equal −51.625 g/m 3 and −1.912 × 10 −5 respectively.The total weight modulus and volume modulus of the three minerals are 5.928 g/m 3 and 3.105 × 10 −6 .Evidently, dolomite dissolution, whose dissolved quantity is larger than that of gypsum, is more likely given priority in this zone while calcite starts to precipitate.Generally speaking, the new space created by the dissolution of dolomite and gypsum is larger than the space filled by the precipitation of calcite, which can be seen from the value of N 2 (3.105 × 10 −6 > 0).

βd >1 ~βg ≤ 1 (Zone III)
In this zone, the carbonate rocks cannot be dissolved by groundwater any more, and the chemical reactions mainly consist of gypsum dissolution, calcite precipitation and the reactions of desulphurization and dedolomitization, which can be presented as the following equations: The solubility of gypsum is high, and dedolomitization needs to consider some conditions.In the Equation ( 23), the exchange quantity of Ca 2+ (in Equation ( 23)) is far lower than the quantity of Ca 2+ (in Equation ( 22)), that is to say, m 6 > m 7 .Combining the Equation ( 22) and (23) then obtains the following equation: According to Equation (24), the relationship plot is drawn (Figure 5), and the following correlation equations can be observed from Figure 5. quantities of precipitation of carbonate rocks minerals was greater than that of sulphate rocks, but a tiny amount of new pore space would still be produced by the interaction between the groundwater and karst aquifer.

Conclusions
This paper considers calcite, dolomite and gypsum as representative minerals of karst rocks in Liulin spring and combines the geology, hydrogeological conditions and hydrogeochemical reactions within the different zones which were divided by saturation indices of calcite, dolomite and gypsum to develop the hydrogeological-hydrogeochemical model of Liulin spring.The two prime issues can be summarized: (1) In zone I where βc ≤ 1, three kinds of minerals (viz.calcite, dolomite and gypsum) are in the dissolved state.The corrosion modulus of calcite is maximum and the second is that of dolomite, and the value of corrosion modulus of gypsum is minimum; in zone II where βc ~1 > βd ≤ 1, the value of corrosion modulus of calcite is negative, which indicates that the reaction of calcite transforms dissolution to precipitation.The corrosion modulus of dolomite and gypsum increased, and the corrosion reaction was dominated by carbonate rocks dissolution; in the third zone where βd > 1 ~βg ≤ 1, calcite and dolomite achieved a supersaturated state and the corrosion modulus of gypsum increased remarkably.(2) From zone I to zone III, the corrosion modulus of calcite decreased gradually along the direction of groundwater flow, however, the corrosion modulus of dolomite increased at first then decreased, and the corrosion modulus of gypsum was augmented with the increase in aquifer depth.The dissolution law of karst groundwater in the Liulin spring area shows that the hydrogeochemical environment plays an important role in mineral corrosion.However, there is less work in the hydrodynamic analysis of the spring system, and more efforts in this regard should be made in the future.

4. 2 . 3 −
βc >1 ~βd ≤ 1 (Zone II)When it comes to zone II, the dissolved quantities of Ca 2+ and CO3  2− are in the supersaturated state, so that the dissolution of calcite does not occur in this zone but the dissolution of dolomite and gypsum still exist.Based on the same principle as zone I, the comprehensive chemical equations (in zone II) about the dissolution of dolomite and gypsum are as follows: m 4 CaMg(CO 3 ) 2 = m 4 Ca 2+ + m 4 Mg 2+ + 2m 4 (15), 2m 4 CaCO 3 (i.e., m 4 CaCO 3 + m 4 Ca 2+ + m 4 CO 3 2− ) represents the precipitation of calcite when m 5 Ca 2+ ≥ m 4 CO 3 2− and the extra quantity of (m 5 − m 4 ) Ca 2+ saves in the karst water, while (m 4 + m 5 ) CaCO 3 (i.e., m 4 CaCO 3 + m 5 Ca 2+ + m 5 CO 3 2− ) represents the quantity of calcite precipitation when m 5 Ca 2+ > m 4 CO 3 2− and the extra quantity of (m 4 − m 5 ) CO 3 2− saves in the karst water.After deducting the content of Ca 2+ dissolved from gypsum from total content in zone II, the rest is the sum content of Ca 2+ which consists of the part inputted from zone I and the other part dissolved from dolomite in zone II.It is difficult to determine the quantity of CaCO 3 which formed by Ca 2+ and CO 3 2− , so the values of A 4 and B 4 are computed by the relationship between Mg 2+ and HCO based on the equation as follows: y 4 = A 4 X 4 + B 4 (16) where, y 4 and X 4 are the content of Mg 2+ and HCO 3 − (mmol/L).A 4 and B 4 are the coefficient of regression and constant, respectively.Both the content of Ca 2+ and SO 4 2− dissolved from gypsum and accumulated from zone I are still few.Assuming that the dissolved quantities of Ca 2+ and SO 4 2− are identical, the correlation equation of Ca 2+ and SO 4 2− is presented as follows:

Table 1 .
Mineral ion concentration and the dissolved quantity of each water sample point.