Next Article in Journal
Water-Use Strategies and Habitat Adaptation of Four Tree Species in Karstic Climax Forest in Maolan
Next Article in Special Issue
Mobile Sources Mixing Model Implementation for a Better Quantification of Hydrochemical Origins in Allogenic Karst Outlets: Application on the Ouysse Karst System
Previous Article in Journal
Reliability Treatment of Silicon in Oilfield Wastewater by Electrocoagulation
Previous Article in Special Issue
Hydrogeochemical Characteristics of Karst Areas: A Case Study of Dongzhuang Reservoir Area in Jinghe River
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Change Analysis of Karst Landforms, Hydrogeological Conditions and Effects of Tunnel Excavation on Groundwater Environment in Three Topography Grades in China

1
State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China
2
China Railway Eryuan Engineering Group Co., Ltd., Chengdu 610031, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(1), 207; https://doi.org/10.3390/w15010207
Submission received: 24 November 2022 / Revised: 28 December 2022 / Accepted: 28 December 2022 / Published: 3 January 2023
(This article belongs to the Special Issue Hydrogeology and Geochemistry of Karst Aquifers)

Abstract

:
One-third of the Earth in China is formed by Karst topography, which exposes different Karst landforms in three topography grades from southeast to northwest, corresponding to below several hundred meters for the first grade, one to two thousand meters for the second grade, and more than 4000 m for the in Qinghai-Tibet Plateau. Through the hydrochemical and D-18O stable isotopes of 64 water samples collected along two railway lines and the topography fractal characteristics of three typical Karst areas in different topography grades, the changes in Karst development degree, changes in groundwater activities, and the influence of tunnel excavation effects on groundwater environment were analyzed. The results indicated that: (1) the Karst development degree and the influence of Karst tunnel excavation on the groundwater environment are somehow similar in the first and second grades, while there are significant differences between the slopes area from second to third grade and the third grade area. (2) In detail, the relatively weaker Karst development degree and flow seeping in the second grade relatively weaken the influences of tunnel excavation, including the distribution pattern of water resources, the groundwater flow field, and water circulation, while the tunnel elevation has little room to rise. (3) There are many large faults in the north-southward direction in the third topography grade, and the transportation lines in the eastern-western direction will inevitably encounter them. In the intersection area, the tunnel excavation has great effects on the groundwater environment. (4) The lighter hydrogen and oxygen isotopes are enriched in Karst water from the first grade to the third grade, indicating that the recharge source of Karst water presents obvious elevation effect.

1. Introduction

Karst covers approximately one-third of the Earth in China. As the altitude rises, different Karst landforms are exposed [1]. Topography altitudes are divided into three grades from southeast to northwest in China. The areas of the first grade have an altitude of about several hundred meters and are mostly distributed in the southeast area such as Guangxi Province and Yun-Gui Plateau situates. The areas of the second grade have an altitude of one to two thousand meters. The highest-level area is located in Qinghai-Tibet Plateau situates, which are more than 4000 m high. The distribution of the three grades is shown in Figure 1 [2,3]. In southern China, there are four typical types of Karst landform compositions, including peak-cluster with depression (the funnel), peak-forest with depression, isolate peak (residual hillock) with Karst plain, and Karst hill with depression. The first three types show landform diversity from Guangxi Province to Yun-Gui Plateau [4,5]. Shallow dish with Hillock becomes the main Karst landform composition in the central area of Yun-Gui Plateau [1]. There are two different viewpoints about Karst landforms in Qinghai-Tibet. Some researchers consider the palaeo Karst of Qinghai-Tibet as the covered Karst. They believe most of the landforms have already been destroyed with topographical uplift, and current Karst features include some stone pinnacles and rock walls [6,7,8,9]. While others consider that the recognizable Karst features in Qinghai-Tibet Plateau are small caves and corrosion fissures formed by physical weathering under the condition of Quaternary periglacial climate [10,11]. Both opinions acknowledge that there are few typical Karst landforms in this area.
The huge altitude difference between these three grades becomes an insurmountable barrier for transportation [12]. As a result, tunnel excavation has become a common need for efficient transportation in the forms of expressways and railways in the Karst plateau and mountain areas [13,14], and consequently, there are a large number of tunnels in southwest China [12]. The depths of tunnels are generally tens of meters to hundreds of meters, and some in the world even exceed 1000 m [15]. The depth of the diversion tunnel of Jinping Hydropower Station in China is 2525 m, which has been designed and is to be constructed [16]. In addition, the railway from Sichuan to Lasa situates on the Qinghai-Tibet Plateau, and the design depth of some Karst tunnels is greater than 2000 m. By the end of 2018, there are 3896 tunnels with a total length of 4300 km in carbonate area [12]. The most significant hydrogeological impact of the tunnel excavation in an aquifer is the drainage effect [17]. If the tunnel is located below the regional groundwater level, springs or underground rivers may be drained, and the groundwater level will drop significantly, which can lead to hydrogeological and environmental impacts on groundwater-dependent ecosystems [18].
Karst environments are characterized by unique landforms related to dissolution and dominant subsurface drainage, which can affect the groundwater flow characteristics [19]. In another word, the most significant hydrogeological impacts, such as barrier and drain effects in an aquifer, are influenced by the Karst development degree and the characteristics of groundwater flow [19]. The barrier effect causes the water table to rise in the upper reaches and the water table to drop in the down reaches [20,21]. The drain effect can extract groundwater, reduce the groundwater table, and cause some environmental and geotechnical consequences [22,23,24]. In this study, we investigated three typical hydrogeology systems along two railways, as shown in Figure 1b,c. The three hydrogeology systems cross the Karst area and are distributed in three topography grades. The geochemistry, isotope characteristics, and geomorphology analysis methods were used to analyze the relationship between Karst landforms and the characteristics of ground flow, including the characteristics of groundwater, hydraulic connection, and hydrodynamic conditions. Then we discussed the changes in Karst development degree, the activities of groundwater, and the effects of tunnel excavation on the groundwater environment in different hydrogeology grades. The results can suggest a better tunnel layout to avoid or reduce the impact of tunnel excavation.

2. Methods

2.1. Hydrochemistry and Isotope Analyses

In August 2011 and August 2018, a series of samples, consisting of some springs, surface water, underground river water, and two warm springs, were collected in the areas where railway lines pass. The locations of samples are shown in Figure 1b. Three typical tunnels, i.e., Pingguo Tunnel, Dongfeng Tunnel, and Batang tunnel, are located in the three topography grades with different Karst landforms. The samples were collected in the strata, including Silurian, Devonian, Carboniferous, Permian System, and Triassic.
All indoor analyses were rapidly performed in the aqueous chemistry facility operated by the Chengdu University of Technology. Anions and cations were analyzed using DIONEX (ICS-1100) and ICP-OES (ICAP6300), respectively. The stable isotopic composition of hydrogen and oxygen was determined by the Isotopic Laboratory of the China Geological Survey. The standard deviation of the isotopic analyses (MAT 253 stable isotope ratio mass spectrometer) is ±0.2‰ for δ18O and ±2‰ for δD.

2.2. Fractal Characteristics of the Karst Landforms

The morphological parameter statistics and fractal characteristics were used to study the Karst development degrees, including morphological parameter statistics of Karst peak, as well as the fractal dimension of Karst depression and contour. The Karst development progress can be divided into three periods, i.e., juvenile stage, maturity stage, and old stage, corresponding to different Karst morphological features, namely, peak–cluster, peak–forest, and isolate peak, as shown in Figure 2 [25,26]. The shape of the curve of N(δ) versus δ can be used to study the Karst development stage. Based on the different height standards (δ), the number of peaks (N(δ)) exceeding the height standard can be counted, and the curve of N(δ) versus δ can be constructed.
The number of Karst peaks (N(δ)) based on the height standard (δ) for different Karst development stages; N1, N2, and N3 represent the number of Karst peaks for the juvenile stage, maturity stage, and old stage, respectively; δ1, δ2, and δ3 represent the height standards of Karst peaks for Juvenile stage, maturity stage, and old stage, respectively.
Relationship curves between the number of Karst peaks of the height standard (δ) for different Karst development stages are represented.
The morphology of Karst depression and contour is divided into three grades (Figure 3 and Figure 4) [25]. The smoother the depression boundary, the stronger the Karst development. The smoother the contour, the weaker the Karst development. Based on the morphology datum, the perimeter and area of Karst depressions can be measured, and Equation (1) can be used to obtain the fractal dimension of Karst depressions.
D = 2 LogL ( δ ) LogA ( δ )
where D is the fractal dimension, δ is the length standard, L(δ) is the perimeter of the depression based on the different length standards, and A(δ) is the area of depression based on the square area δ2. According to the range of D values from 1 to 2, the Karst development can be divided into three stages. The strongest stage has the smallest value of D, and the weakest stage has the largest value of D.
Based on the different length standards (δ), the numbers of segmentations (N(δ)) can be determined, and the fractal dimension is defined as the gradient of line relationship between lnN(δ) and δ. The higher the gradient, the more complex the topography, and the stronger the Karst development.

3. Description of the Study Site

The study was carried out along the area of three typical railway tunnels within two railway lines (Figure 1b,c). Line ①, Yun-Gui railway, is located from the lowest to the second topography grade, and Line ②, Western Sichuan traffic corridor, is located in the slope zone between the second and highest topography grades. The geological characteristics of these three railway tunnels are summarized as follows:
Pingguo Tunnel in railway line ① passes a zone with a typical landscape of peak-cluster and depression (Figure 5a), and the elevation of this area ranges from ~100 to 800 m. The rivers flow around the Karst mountains, following the direction of strata extension. The Devonian and Carboniferous strata are located in the eastern area, and the Permian stratum is located in the western area and is cut by the Youjiang River. The total thickness of Devonian and Carboniferous strata is about 500~1600 m. Limestone dominates and is interbedded with dolomite, which is the main lithology in this area. Some springs with a large volume of water appear near the banks of the river, and there are two Karst tunnels with several ponors lying above. Faults are relatively limited in the region, and their extension direction is consistent with the strata, thus making them a potential groundwater runoff channel.
Dongfeng Tunnel in railway line ① passes the zone with shallow depressions and dwarf hills (Figure 5b), and the elevation of this area is about 1500~2600 m. This area is mainly covered by the early Paleozoic Cambrian and Devonian strata, whose lithology is dolomite and impure limestone. Groundwater is buried more deeply in the middle area. Some springs with a volume of 1.2 to 10 L/s appear where rivers cut the rocks at the boundary of this system. The groundwater in this area discharges more dispersedly. The northeast striking regional fault forms a lower terrain, and the lithological variations become the confining boundary, and most of the groundwater discharges nearby.
The tunnel lithology information is shown in Figure 6. Batang Tunnel crosses the plateau canyon landscape (Figure 5c), and the altitude difference can reach 2500 m. In this area, some big glacially eroded depressions and sharp peaks appear, and some long and wide dissolution fissures are also found in rocks. However, Karst caves and holes are rare. Paleozoic strata of crystalline limestone and marble are the main strata here. Two big fracture systems extend from north to south in this area, which has strong water-bearing characteristics. A large regional fault forms multiple parallel branches in the area, as well as most the groundwater transport along the fault. When the rivers and valleys cut the fracture system, some springs with big water flow are found. About 70% of the Karst phenomenon, including springs with a volume of 5 to 500 L/s and small Karst caves, are located at the fault or close to it.

4. Results

4.1. Chemical and Isotope Characteristics of Groundwater

The saturation index of calcium carbonate (SI-CaCO3) of 64 samples is shown in Figure 7. The value of SI-CaCO3 ranges from −2.448 to 0.7429. The SI-CaCO3 of the surface river and underground river are lowest in zone ⑤. Both surface and subsurface rivers are unsaturated, and SI-CaCO3 of the underground river is higher than that of the surface river. SI-CaCO3 of groundwater samples decreases from southeast to northwest and can be divided into four zones. The SI-CaCO3 of zone ① is greater than zero, indicating that the Karst groundwater in the south area can be saturated. The SI-CaCO3 of zones ②, ③, and ④ are less than zero, indicating the Karst groundwater samples in this area, especially the groundwater samples from the northwest area, are unsaturated. The lowest SI of this area is about −0.75, which is higher than the SI-CaCO3 of rivers.
The relationship between δ18O and δD is plotted in Figure 8, which also shows the global meteoric water line with δ D = 8 δ 18 O + 10 [27], the Southwestern China meteoric line with δ D = 1.96 δ 18 O + 9.53 [28], and the meteoric water line in the middle-lower region of the Yarlung Zangbu River with δD = 7.54δ18O + 15.92 [29]. The trend line of all samples is also plotted in this figure, which is δ D = 7.704 δ 18 O + 3.95 , and the value of R2 is 0.96.
Most of the samples are located near the globe meteoric water line, and even two warm springs of No. 37 and No. 38, δ18O values do not show positive shifts, and the δ18O well responds to elevation change. The elevation effect of δ18O (Equation (2)) [30] can infer the supply altitude of groundwater. The distribution of δ18O vs. δD is clearly in three zones.
H * H 0 100 × ( δ 18 O * δ 18 O 0 ) = 0.27
where H* and H0 represent the altitude of the analysis sample and standard sample, respectively; δ18O* and δ18O0 come from the analysis sample and standard sample, respectively.
The elevation difference between the supply area and the appearing area is shown in Figure 9. From the figure, the elevation difference increases from the southeast to the northwest. The southeast area locates in the lowest topography grade, and the elevation ranges from 194 m to 419 m, with a difference of smaller than 500 m. The northwest area locates in the highest topography grade, and the largest elevation difference can reach 2213 m. The middle area locates in the second topography grade with an elevation gap of up to 704 m, and about half of the springs’ elevations are more than 500 m.
ΔH is the difference between supply and appearing elevations of groundwater, and H is the appearing elevation of groundwater.

4.2. The Equilibrium of Ca2+

(1)
Pingguo System
Pingguo Tunnel is located in the northwest area of Pingguo Province county, as shown in Figure 5(a-1). This area lies in the lowest topography grade area of China. The insoluble stratum and rivers are considered boundaries of the hydrogeological system through which the tunnel passes. The area of this system is about 141 km2. According to the water-bearing characteristics of groundwater, this area can be divided into two regions, i.e., the southeastern region with strong groundwater bearing characteristics and the northwestern region with weak groundwater bearing characteristics. The typical Karst landscapes in this area include peak forests, isolated peaks, and depressions. Two underground rivers are lying in the south area, with several geological windows connected with them. The water samples in this area include some springs and groundwater from the geological windows, as shown in Figure 5(a-1).
This area is divided into three subsystems of groundwater, as shown in Figure 5(a-1). The saturation index of calcium carbonate has the same trend as that of the concentration of Ca2+. In the sub-hydrogeology system, No. 5 and No. 7 are spring water, and No. 3, No. 4, and No. 6 are underground water collected through the geological windows. The concentration of Ca2+ in the samples is in the order of No. 3 < No. 4 < No. 5 < No. 6 < No. 7, while the SI-CaCO3 is in the order of No. 2 < No. 3 < No. 4 < No. 7 < No. 6. These trends indicate that the flow direction is from the northwestern region to the southeastern region, as shown in Figure 5(a-1). It supports the subsystems division of geology. The layout of the two tunnels is shown in Figure 10a. The north layout crosses the subsystem, while the south layout does not cross the subsystem.
(2)
Dongfeng System
Dongfeng Tunnel is located in the west area of Mile County in Yunnan Province, as shown in Figure 5(b-1). This area lies in the second topography grade area of China. The insoluble rock and topography watershed are considered the boundaries of the hydrogeological system through which the tunnel passes. The area of this system is about 739.5 km2. The groundwater-bearing characteristics are used to divide this area into three parts, and the strata of Carboniferous and Devonian are divided into two zones. The deeply buried groundwater is in the middle zone of this water-bearing stratum, while shallowly buried groundwater is in the surrounding zone. The stratum of Permian has strong groundwater-bearing characteristics and lies in the east area near the strata of Carboniferous and Devonian. The water samples in this area include 14 springs including two warm springs, and most of them appear in the stratum of Permian, and the boundary between the Permian stratum and Quaternary is shown in Figure 5(b-1).
In this hydrogeological system, the saturation index of Calcium carbonate is all below zero, ranging from −0.6168 to −0.0558. In every subsystem, the SI and the concentration of Ca2+ do not have obvious differences, and both are lower than those of the hydrogeology system of the Pingguo System.
(3)
Batang System
Batang Tunnel is located in the western area of Batang County, close to the boundary between Sichuan Province and Tibet, as shown in Figure 5c. This area lies in the slope zone from the second to first topography grade area. The insoluble rock and rivers are considered the boundaries of the groundwater system, through which the tunnel passes. The area of this system is about 1521.2 km2. The groundwater-bearing characteristics can be divided into three grades. The strata of Carboniferous and Devonian have the strongest groundwater-bearing characteristics, followed by the Silurian stratum, and the stratum of Permian has the weakest characteristics. Thirteen water samples were collected in this area, including springs and surface water. Most of the samples were from the strata of Carboniferous and Devonian along the faults, and some samples were collected from the boundary area between the stratum of Permian and the insoluble rock, as shown in Figure 5c.
In this hydrogeological system, the saturation index of calcium carbonate is below zero, ranging from −1.3161 to −0.2733. This hydrogeological system has a significantly lower saturation index than the hydrogeology system of the Dongfeng System. The SI-CaCO3 in subsystem Ⅱ-2 has the trend of No. 57 < No. 56 < No. 53 < No. 52. In contrast, in subsystem Ⅱ-1, the SI-CaCO3 has the trend of No. 62 < No. 61 < No. 60 < No. 59. In these subsystems, the concentration of calcite trended opposite to SI-CaCO3.

4.3. Fractal Characteristics of Typical Karst Areas

In this study, the topography data of SRTMDE.M 90M were used to count the area and circumference of the polygon connecting the Karst peaks and saddle zones around Karst depressions. There are 99 Karst depressions in the Pingguo hydrogeological system, among which 66 are in the southeastern area, and 33 are in the northeastern area. δ and δ2 were used to measure the circumference and area, and the relationship between lnL(δ) and lnA(δ2) is shown in Figure 11a. The gradients of the line of lnL(δ) and lnA(δ2) are the fractal dimensions. The total fractal dimension, the southeastern fractal dimension, and the northwestern fractal dimension are 1.3818, 1.4122, and 1.3561, respectively, and the determination coefficients are 0.8499, 0.8847, and 0.8016, respectively.
There are 63 Karst depressions in the Dongfeng hydrogeological system, among which 36 are in the groundwater deeply buried area, and 27 are in the groundwater shallowly buried area. δ and δ2 were used to measure the circumference and area, and the relationship between lnL(δ) and lnA(δ2) is shown in Figure 11b. The total fractal dimension, the fractal dimension in the deeply buried area, and the fractal dimension in the shallowly buried area are 1.4628, 1.43412, and 1.5159, respectively, and the determination coefficients are 0.903, 0.9185, and 0.902, respectively.
The relationships between ln(δ) and ln(N(δ)) in these three typical study areas are shown in Figure 12. The values of δ are 50 m, 100 m, 300 m, 500 m, 1000 m, 3000 m, 5000 m, and 7000 m. Additionally, N(δ) comes from all elevation lines in the study area. According to the fractal dimension in subsystems, the gradients of the straight line are 1.501 and 1.394 in the south area and north area, respectively, in the Pingguo System. In the Dongfeng System, the gradients are 1.126 and 1.123 in shallowly and deeply buried areas, respectively. In the Batang system, the gradients are 1.063, 1.083, and 1.101 in different strata.

5. Discussions

5.1. Karst Development Characteristics

Typical Karst landforms in the Pingguo System are peak forests, depressions, and isolated hills [4,5]. In the Karst tunnel shown in Figure 5(a-1), there are a few springs with considerable volume, a Karst well, and geological windows. In the area of 141 km2, there are 99 depression polygons and more than 870 Karst peaks, as shown in Figure 5(a-2). More than two-thirds of them are located in the southeast area where the Karst tunnel passes through. The density of peaks and depressions, and the underground river conduits’ appearance, all show strong Karst development. In addition, the results indicated that the topography environment has a more significant influence on the Karst development in this area than the strata [31]. The typical landform in Dongfeng System is shallow plate depressions and hillocks. In the area of 739.5 km2, there are 63 depressions and more than 597 peaks, as shown in Figure 5(b-2). In this area, the density of peaks and depressions is lower than Pingguo System. There are more springs in this area than in Pingguo System. A larger number of springs with small volumes indicate that the groundwater doesn’t intensively excrete. Based on the density of the Karst depressions and peaks and the occurrence of springs in this area, it can be inferred that the Karst development is weaker than that in the Pingguo System. The elevation of the Batang system ranges from 2800 m to 5200 m, the average temperature is about 5 °C, and more than half of the area is higher than the snow line. Currently, there are few Karst depressions in the carbonate area and only a few big valleys distributed in the high elevation area, probably formed by glacier erosion [11,32,33,34]. There are many springs in this area, and the volume of springs ranges from 0.1 L/s to 500 L/s. All springs emerge from the fracture and faults. The real number of the springs is bigger than shown in Table 1 because most areas are challenging to test. Since no typical Karst landforms emerge, it can be inferred that the Karst development is weak in this area, maybe because it is a Karst in a phase of greater youth.
The fractal characteristics can be used to study the Karst development [25,35,36,37], and the relationship between the number of peaks and the elevation is used to evaluate the Karst development stage [25,35]. As shown in Figure 12, the relationship in the Pingguo area is similar to the type (2) curve in Figure 2b, which has a slow decreasing trend. However, the relationship in the Dongfeng area is similar to the type (3) curve, with a gentle incline first, then a steep incline, and finally a gentle incline. As shown in Figure 2a,b, curve (3) shows that there are a large number of Karst peaks, and the mountain and the foundation are integrated. Curve (2) shows that the peak bodies are independent, and the basement is integrated. The result shows that Pingguo system is more evolved than Dongfeng system, and the Karst development of the Pingguo System is stronger than that of the Dongfeng System.
The gradient of the relationship line between lnL(δ) and logA(δ) represents the fractal dimension of Karst depressions [25,35]. The fractal dimension of the Pingguo System is lower than that of the Dongfeng System, as shown in Figure 11. The smaller the fractal dimension, the stronger the Karst development. The fractal dimension in the southeast area is lower than that in the northwest area, which corresponds to the groundwater activity that most springs with big volumes appear in the northwest area (Figure 5(a-1)). In the Dongfeng System, the fractal dimension in the area with shallowly buried groundwater is lower than that in the area with deeply buried groundwater, probably because the springs emerge in the shallowly buried area where the groundwater flow converges. Since it is not a typical Karst depression, the fractal dimension of the Batang system was not obtained.
Considering all the elevation lines with an elevation difference of 50 m shown in Figure 10, the relationship between lnδ and lnN(δ) was determined, as shown in Figure 13. The gradient of trend lines is considered as the fractal dimensions of elevation lines. The gradients of different systems are in the order of Pingguo System > Dongfeng System > Batang system. For the Pingguo System and Dongfeng System, the Karst landform is Karst depressions, peak-cluster, Karst holes, and groundwater rivers, and the contours are influenced by Karst actions. The bigger the gradient of trend lines, the more complex the contours, and the stronger the Karst development. Therefore, the Karst development in the Pingguo System is stronger than that in the Dongfeng System. In the Batang system, the obvious Karst landforms are Karst depressions, while Karst holes do not appear. Thus, the Karst action has a weak impact on topography, and the gradient of the trend line is the smallest. The strong glaciation and wind erosion in the high elevation area will cause complex contours, but their effects on the contours are weaker than the Karst action.

5.2. Difference in the Groundwater Environment in Three Topography Stages

The saturation index of calcium carbonate and concentration of Ca2+ can be used to analyze the seeping conditions [38,39,40]. The SI-CaCO3 of all groundwater samples (Table 1) are higher than that of surface water samples (in Figure 7), suggesting the water–rock reaction during the groundwater flow. In a closed environment, CO2 in groundwater cannot be replenished in time, and saturation is often not reached in this condition. The SI-CaCO3 of springs shows a reducing trend from southeast to northwest, as shown in Figure 5(a-1) and Figure 7. In the Pingguo System, a few springs can achieve the calcite saturation phase. It means that this area has a relatively open condition, which is caused by enough CO2 supply through many Karst depressions, geological windows, and Karst wells. Along the Yun-Gui railway from east to west, the springs appear from the first stage of topography to the second stage. Near the first stage, a small number of springs in the Pingguo System reach calcite saturation, but in the area further from the first stage, most springs cannot reach saturation, as shown in Figure 7. This phenomenon might be caused by the lower air temperature and less rainfall in the second topography stage area [31,41]. Adapting to this environment, Karst landform becomes shallow plate depressions and hillocks [42,43], the Karst development is weaker, and groundwater is often buried relatively deep in the supply area, as shown in Figure 5(b-1). It means that the environment for groundwater flow is relatively closed, and the supply of CO2 is insufficient. When springs appear in the slope area from the second to the third topography stage, all springs do not reach saturation shown in Figure 7, and the SI-CaCO3 of No.62 is reduced to −0.9495. In this area, typical Karst landforms, such as Karst depressions and big Karst holes, are not found. Many springs with considerable capacity appear from the slope area instead of the bottom of the valley. In this area, groundwater flows are controlled by large faults and dense cracks. Especially, along the large faults, the groundwater can flow for a long way in a relatively closed environment.
The distribution of δD versus δ 18O is influenced by the topographical uplift. The δD and δ 18O decrease from southeast to northwest, and concentrate into three zones, as shown in Figure 8, which correspond to the three topography stages. All samples are located near the trend line, which is similar to the southwestern China meteoric line, and the R2 can reach 0.96. This suggests that there is no “18O shift” in these samples, and the elevation effect of δD versus δ 18O can be used to infer the supply elevation of groundwater wells [30]. The elevation difference indicates the hydrodynamics of groundwater. Additionally, the major ions and water molecule isotopes could be used to investigate the hydrodynamic function of Karst aquifer system [44,45]. Overall, the trend of ΔH is in the following order: the third stage > the second stage > the first stage, which shows that the higher topography area may provide stronger hydrodynamics and extend the groundwater flow path. The value of ΔH increases weakly from the first stage to the second stage and is concentrated in these two regions with respect to the third stage area. The ΔH values in the third stage range from ~500 m to 2500 m. This is caused by the length and the width of fractures and faults, and the large faults control the big volume springs with long flow paths.

5.3. Influence of Tunnels on the Groundwater Environment

Tunnels’ influence on the groundwater environment is related to not only the Karst development and the characteristics of groundwater flow, but also the tunnel layout [12,18,46]. In addition, the impacts of tunnel excavation on the groundwater include the distribution pattern of water resources, the groundwater flow field, water circulation, groundwater hydrogeochemistry, and water environment [22,47,48]. All these effects occur when the elevation of the tunnel is lower than the groundwater level.
In the first stage topography area, such as Pingguo System, according to the analysis from Section 5.1 and Section 5.2, it can be inferred that the Karst is strongly developed and groundwater flows rapidly in this area. It means that the tunnel construction crossing through this area can influence the environment easily and strongly, but different tunnel layouts lead to different influences on the groundwater environment. In Figure 5(a-1), there are two types of tunnel layouts with similar tunnel elevations. The northern layout passes through the hydrogeologic unit I, while the southern layout passes through the hydrogeologic unit II. The presence of springs, the Karst landforms, and the SI-CaCO3 all imply that there is a quicker flow, a more groundwater-bearing stratum, and a lower groundwater level. Whether the north or south layout is used, some portion of the tunnel will be below the groundwater level. During the tunnel construction period, tunnel drainage destroys the aquifer structure, and the groundwater will excrete into the tunnel. As a result, the groundwater level will drop and the springs and the groundwater of geological wells maybe disappear, or the volume of springs may decrease in the hydrogeological unit. All these influences are weaker in the hydrogeological unit Ⅱ. In the second-stage topographical area, Karst development, groundwater flow condition, and groundwater-bearing are weaker. All these conditions might cause the vadose zone thinner than the first grade, so tunnels are usually below the groundwater level, and their construction will lower the groundwater level more strongly. If the tunnel does not pass through the discharge area as shown in Figure 5, the spring volume will decrease to some extent but will not dry up. If the tunnel passes through an area covered with the Quaternary aquifer, the Karst collapse will be caused by a drop in the ground level.
In the first grade, the relief of the terrains is relatively small, and the Karst groundwater level is deep and gentle. When the tunnel passes through in the form of shallow burial, the length of the section below groundwater level is limited, and the influence on the Karst water environment is weak, which may lead to the reduction of the flow of some spring points. In the deep-buried tunnel form, the tunnel may be located below the groundwater level for a long distance and reveal the main Karst pipeline, causing the reduction of the groundwater level in a large area and even dry off the springs or underground rivers. The north line of Pingguo Tunnel is parallel to the main runoff direction of s5, s7 springs. Since part of the tunnel is lower than the groundwater level, the flow rate of the springs may be sharply reduced. While the south plan does not cross the runoff area of spring s5, s7, it mainly affects the discharge of the spring s2 by passing through the runoff area of s2. Due to the high Karst development in this area, the groundwater level fluctuates gently. The influence of tunnel drainage on groundwater environment can be effectively reduced or even eliminated by raising the tunnel elevation or avoiding the main runoff channel of groundwater.
In the second grade, the topography fluctuates significantly, and the depth of the groundwater level fluctuates somewhat with the topography. Dongfeng Tunnel is completely below the groundwater level when it passes through, and its distribution direction is basically consistent with the direction of groundwater migration, which will have an impact on spring points within unit II and unit III. If the tunnel direction cannot be changed, the change in tunnel elevation could only change the influence degree of groundwater environment.
In the third grade, the topographic relief is huge and the buried depth of the groundwater level is relatively shallow, but changes with the topographic relief significantly, tunnel construction might be facing a huge head pressure. When the tunnel passes through the slope area from the second to the third stage, as shown in Figure 5c, the typical Karst landforms are not present in the Batang system. The groundwater is bearing in rocks’ fractures, and the springs with giant volumes are generally controlled by the large fracture and can flow a long distance. The intense topographical uplift in this area is caused by the collide between the Indian Ocean plate and Eurasian plate, which caused a large number of north–south fractures [49]. If the railway line goes from the second to the third terrain phase, there is no way to avoid this type of fracture. In this area, the groundwater level is high due to the presence of fracture-bearing strata, and tunnels are often located below the groundwater level. The tunnel construction can cause the groundwater level to decrease. Especially, in the area the big fracture passes through, the influence can extend several kilometers [50]. In contrast, in other areas without big faults, the groundwater level and the spring volume only have a weak decrease. Figure 5c shows different layouts with similar elevations in the Batang system, one near the boundary of the river and one in the watershed area. Both cannot avoid meeting with big faults. However, the south layout with a big depth is located in the groundwater recharge area, where the finite depth of fracture development makes weak seeping in deep depth and the influence of the tunnel on the groundwater environment is weak. Relatively speaking, the north layout with a small depth is located in the outlet area of groundwater, where the big fault area has been a strong bearing area for groundwater, and the range of groundwater level is reduced, especially for the fault extending orientation.
Figure 14 shows the relationship between the tunnels’ elevation and groundwater level in three topography grades, and it can be inferred that: in the first grade, groundwater level is low and gentle, while the topographic elevation difference is within 500 m, the water level difference might be less than 100 m, makes it effective to control the tunnel groundwater environment influence by changing the elevation of tunnel. In the second grade, groundwater level changes with the relief variation of the terrain, while the topographic elevation difference ranging from 1000 to 2000 m, and the water level difference could reach hundreds meters, making it hard to control the influence only by changing the tunnel elevation. In the third grade, the huge topographic relief with more than 2000 m of altitude difference also brings a huge water level gap with more than 1000 m. Under such conditions, the variation of tunnel elevation can only extremely limited reduce the impact on groundwater environment

6. Conclusions

From the chemical, isotopic, and topography data discussed in this article, Karst development, Karst landforms, and groundwater activities are obviously different in the three topography grade areas of China. The influence of Karst tunnel excavation on the groundwater environment is similar to some extent, and some small differences are present as the terrain rises from the southeast to the northwest area, as shown in Figure 15. In addition, the tunnel layout is a notable factor that may cause detailed differences in impacts.
With the topographical decrease from the second to the first topography area, the Karst development becomes stronger and more complex to some extent. The Karst density and morphology, groundwater confluence, and excretion also suggest that the water-bearing of strata and the flow seeping are stronger in the lower topography. However, if the tunnel elevation is lower than the groundwater level, the influence of groundwater on the environment, including the distribution pattern of water resources, the groundwater flow field, and water circulation, is similar, as shown in Figure 14a,b. In detail, if the tunnel passes through the area near the excretion area of groundwater, depending on the development of Karst, there is a possibility that springs will dry up in the lower topography, and the impact of a drop in the groundwater level might be more significant in this area. Rising the tunnel elevation and avoiding some special sub-hydrogeology systems can effectively reduce the influence. In addition, the relatively weaker Karst development and flow seeping make the vadose zone thinner, and there is less room to raise the tunnel elevation in the second topography area.
The topography of the slope area is obviously different from the second grade to the third grade, and there are also obvious differences in the area of the third grade. Typical Karst landforms do not appear in most carbonate rocks area. The groundwater bearing and flow are controlled by the concentrated fracture and giant faults. Based on the huge altitude difference in this area, most tunnels must be below the groundwater level. The strong topography uplift in this area is caused by the collide between the Indian Ocean plate and the Eurasian plate, which leads to a large number of north–south orientation fractures that are not avoidable for transportation routes from east to west. Tunnel excavation encounters giant faults in this area, causing the range of groundwater level to decrease for several miles along the faults extending direction and also diminishing springs controlled by the giant faults. If the tunnel layout is located near the watershed, the deeply buried environment closes the faults, and water hardly seeps into the tunnel, effectively reducing the influence.

Author Contributions

Conceptualization, Q.Z.; Methodology, Q.Z.; Formal analysis, X.L.; Investigation, C.Q. and L.Y.; Writing–original draft, Y.T.; Writing–review & editing, J.Q.; Project administration, M.X. and D.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the Natural Science Foundation of China’s NSFC Grant (No. 41402223).

Data Availability Statement

All data, models, and code generated or used during the study appear in the submitted article.

Acknowledgments

The authors are grateful to the Sichuan Institute of Geological Engineering Investigation Group Co., Ltd.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, S.J.; Zhang, X.B.; Bai, X.Y. An outline of karst geomorphology Zoning in the karst areas of southern China. Mt. Res. 2015, 33, 641–648, (In Chinese with English abstract). [Google Scholar]
  2. Wang, P.X. Deformation of Asia and global cooling: Searching links between climate and tectonics. Quat. Sci. 1998, 3, 213–221, (In Chinese with English abstract). [Google Scholar]
  3. Ying, Y.W.; Tian, W.S.; Ju, J.H. Spatial-Temporal characteristics of temperatures on different terrains of southwest China. Adv. Clim. Change Res. 2010, 6, 429–435, (In Chinese with English abstract). [Google Scholar]
  4. Zeng, Z.X.; Cheng, M.H.; Yao, Q.Y. The karst morphological characters in south China. Carsologica Sinica 1982, 1, 27–32, (In Chinese with English abstract). [Google Scholar]
  5. Ren, M.E. Geomorphic Inheritance and the Development of Tower Karst; The Commercial Press: Beijing, China, 1987; p. 82. [Google Scholar]
  6. Zhang, D. A morp342.hological analysis of Tibetan limestone pinnacles: Are they remnants of tropical karst tower and cones? Geomorphology 1996, 15, 79–91. [Google Scholar]
  7. Li, D.W.; Cui, Z.J.; Liu, G.N. Existence of Palaeokarst on Tibet plateau and its correlation with the karst of the eastern district. Carsologica Sin. 1999, 18, 309–318, (In Chinese with English abstract). [Google Scholar]
  8. Cui, Z.J.; Li, D.W.; Feng, J.L.; Liu, G.N.; Li, H.J. The covered karst, weathering crust and karst(double-level) planation surface. Sci. China 2002, 45, 365–379. [Google Scholar] [CrossRef]
  9. Hao, Y.; Cao, B.; Zhang, P.; Wang, Q.; Li, Z.; Yeh, T.-C.J. Differences in karst processes between northern and southern China. Carbonates Evaporites 2012, 27, 331–342. [Google Scholar] [CrossRef]
  10. Wang, F.B. Some problems of karst on Qinghai-Xizang Plateau. Mt. Res. 1991, 9, 65–72, (In Chinese with English abstract). [Google Scholar]
  11. Zhu, X.W. The nature of Tibet plateau karst and the query concerning “relict fenglin karst”. Carsologica Sin. 1994, 13, 220–227. [Google Scholar]
  12. Lv, Y.X.; Jiang, Y.J.; Hu, W.; Cao, M.; Mao, Y. A review of the effect of tunnel excavation on the hydrology, ecology, and environment in karst areas: Current status, challenge, and perspectives. J. Hydrol. 2020, 586, 124891. [Google Scholar] [CrossRef]
  13. Vincenzi, V.; Gargini, A.; Goldscheider, N. Using tracer tests and hydrological observations to evaluate effects of tunnel drainage on groundwater and surface waters in the Northern Apennines (Italy). Appl. Hydrogeol. 2009, 17, 135–150. [Google Scholar] [CrossRef]
  14. Liu, J.; Shen, L.; Wang, Z.; Duan, S.; Wu, W.; Peng, X.; Wu, C.; Jiang, Y. Response of plants water uptake patterns to tunnels excavation based on stable isotopes in a karst trough valley. J. Hydrol. 2019, 571, 485–493. [Google Scholar] [CrossRef]
  15. Masset, O.; Loew, S. Quantitative hydraulic analysis of pre-drillings and inflows to the Gotthard Base Tunnel (Sedrum Lot, Switzerland). Eng. Geol. 2013, 164, 50–66. [Google Scholar] [CrossRef]
  16. San, Z.G.; Zhou, C.H.; Yong, G. The Karst Hydrogeology Condition and Engineering Effects on Yanlongjiang Hydropower Station; Science Press: Beijing, China, 2018; p. 5. [Google Scholar]
  17. Vázquez-Suñé, E.; Sanchez-Vila, X.; Carrera, J. Introductory review of specific factors influencing urban groundwater, an emerging branch of hydrogeology, with reference to Barcelona, Spain. Hydrogeol. J. 2005, 13, 522–533. [Google Scholar] [CrossRef]
  18. Kim, Y.-Y.; Lee, K.-K.; Sung, I. Urbanization and the groundwater budget, metropolitan Seoul area, Korea. Hydrogeol. J. 2001, 9, 401–412. [Google Scholar] [CrossRef]
  19. Gutiérrez, F.; Parise, M.; De Waele, J.; Jourde, H. A review on natural and human-induced geohazards and impacts in karst. Earth-Science Rev. 2014, 138, 61–88. [Google Scholar] [CrossRef]
  20. Shin, J. Numerical modeling of coupled structural and hydraulic interactions in tunnel linings. Struct. Eng. Mech. 2008, 29, 1–16. [Google Scholar] [CrossRef]
  21. Deveughèle, M.; Zokimila, P. Impact of an impervious shallow gallery on ground water flow. Bull. Eng. Geol. Environ. 2010, 69, 143–152. [Google Scholar] [CrossRef]
  22. Li, H.; Kagami, H. Groundwater level and chemistry changes resulting from tunnel construction near Matsumoto City. Environ. Geol. 1997, 32, 76–84. [Google Scholar]
  23. Vincenzi, V.; Gargini, A.; Goldscheider, N.; Piccinini, L. Differential Hydrogeological Effects of Draining Tunnels Through the Northern Apennines, Italy. Rock Mech. Rock Eng. 2014, 47, 947–965. [Google Scholar] [CrossRef]
  24. Butscher, C.; Huggenberger, P.; Zechner, E. Impact of tunneling on regional groundwater flow and implications for swelling of clay–sulfate rocks. Eng. Geol. 2011, 117, 198–206. [Google Scholar] [CrossRef]
  25. Wang, D.; Xu, M. Analysis on Morphologic Features of the Peak-Cluster Depression in Qiubei, Southeast Yunnan. Log. Sin. 2010, 29, 239–245, (In Chinese with English abstract). [Google Scholar]
  26. Zhu, X.W. New considerations on characteristics and evolution of fenglin karst. Carsologica Sin. 1991, 10, 171–178, (In Chinese with English abstract). [Google Scholar]
  27. Craig, H. Isotopic Variations in Meteoric Waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef] [PubMed]
  28. Zhang, H.P.; Liu, E.K.; Wang, D.S.; Jia, Y.K.; Sun, J.C. Composition of stable isotopes of precipitation and controlling factors in China. Bull Inst. Hydrogeol. Eng. Geol. CAGS 1991, 7, 101–110, (In Chinese with English abstract). [Google Scholar]
  29. Wang, J.; Liu, T.C.; Ying, G. Characteristics of isotope distribution in precipitation in the middle-lower reaches of Yarlung Zangbo River. Geol. Geochem. 2000, 28, 63–76, (In Chinese with English abstract). [Google Scholar]
  30. Yu, J.; Yu, F.; Liu, D. Hydrogen and oxygen isotopic compositions of meteoric waters in the eastern part of China. Chin. J. Geochem. 1987, 6, 367–371. [Google Scholar] [CrossRef]
  31. Zuo, J.H.; Luo, Q.; Guo, X.M. Influence on karst development by the three topography grades in China. Sci. Technol. Inf. 2012, 18, 2010–2013. (In Chinese) [Google Scholar]
  32. Gao, Q.Z.; Cui, Z.J.; Tao, Z.; Liu, G.N.; Hong, Y. The nature, formation age and genetic environment of the Palaeokarst on Qinghai-Xizang Plateau. Acta Geogr. Sin. 2002, 57, 267–274, (In Chinese with English abstract). [Google Scholar]
  33. Franco, C.; Furio, F. Karst landforms in Friuli Venezia Giulia: From Alpine to coastal karst. Landsc. Landf. Italy 2017, 12, 147–156. [Google Scholar]
  34. David, J.A.E.; Nigel, A.; Emrys, P. Glacial geomorphology of neutral hills uplands, southeast Alberta, Canada: The pro-cess-form imprints of dynamic ice streams and surging ice lobes. Geomorphology 2020, 35, 1–40. [Google Scholar]
  35. Gao, W.; Qi, J.H.; Xu, M.; Li, Y.X.; Wang, N.F. Influence of fractal characteristics of a surface karst landform on tunnel construction. Mod. Tunn. Technol. 2016, 53, 35–43, (In Chinese with English abstract). [Google Scholar]
  36. Li, C.S.; Liao, Y.K.; Ding, J.F. Collapse and water inflow characteristics of a tunnel in metamorphic schist in the nappe structure region. Mod. Tunn. Technol. 2017, 54, 24–44, (In Chinese with English abstract). [Google Scholar]
  37. Li, C.S.; Liao, Y.K.; Ding, J.F. Feasibility of research on chemical characteristics of karst groundwater by using fractal theory. J. Eng. Geol. 2018, 26, 464–472, (In Chinese with English abstract). [Google Scholar]
  38. Shen, Z.L. Hydrogeochemistry; Geology Press: Beijing, China, 1993; p. 42. [Google Scholar]
  39. Peyraube, N.; Lastennet, R.; Denis, A. Geochemical evolution of groundwater in unsaturated zone of karstic massif, using the Pco2-Sic relationship. J. Hydrol. 2012, 430–431, 13–24. [Google Scholar] [CrossRef]
  40. Peyraube, N.; Lastennet, R.; Denis, A.; Malaurent, P.; Villanueva, J.D. Interpreting CO2–SIc relationship to estimate CO2 baseline in limestone aquifers. Environ. Earth Sci. 2014, 72, 4207–4215. [Google Scholar] [CrossRef] [Green Version]
  41. Li, D.W.; Cui, Z.J.; Liu, G.N. Feature and origin of covered karst on hunan, Guangxi, Guizhou, Yunnan and Tibet. J. Mt. Sci. 2000, 18, 289–295, (In Chinese with English abstract). [Google Scholar]
  42. Wang, Y.; Wang, Z.X. Patterns of karst geomorphologica combinations in areas with rich groundwater. Carsologica Sin. 2015, 34, 314–324, (In Chinese with English abstract). [Google Scholar]
  43. Wang, Y. Vertical zoning of groundwater runoff system in karst plateau. Carsologica Sin. 2018, 37, 1–8, (In Chinese with English abstract). [Google Scholar]
  44. Hamed, Y.; Ahmadi, R.; Demdoum, A.; Bouri, S.; Gargouri, I.; Hamed, B.D.; Al-Gamal, S.; Laouar, R.; Choura, A. Use of geochemical, isotopic, and age tracer data to develope models of groundwater flow: A case study of Gafsa mining ba-sin-Southern Tunisia. J. Afr. Earth Sci. 2014, 100, 418–436. [Google Scholar] [CrossRef]
  45. Hamed, Y.; Dassi, L.; Ahmadi, R.; Ben Dhia, H. Geochemical and isotopic study of the multilayer aquifer system in the Moulares-Redayef basin, southern Tunisia/Etude géochimique et isotopique du système aquifère multicouche du bassin de Moulares-Redayef, sud tunisien. Hydrol. Sci. J. 2008, 53, 1241–1252. [Google Scholar] [CrossRef]
  46. Yuan, D.X. Discussion on karst environment system. Carsol. Sin. 1988, 3, 179–186, (In Chinese with English abstract). [Google Scholar]
  47. Chae, G.T.; Yun, S.T.; Choi, B.Y.; Jo, Y.; Mayer, B.; Kim, Y.; Lee, J.Y. Hydrochemistry of urban groundwater, Seoul, Korea: The impact of subway tunnels on groundwater quality. J. Contam. Hydrol. 2008, 101, 42–52. [Google Scholar] [CrossRef] [PubMed]
  48. Butscher, C. Steady-state groundwater inflow into a circular tunnel. Tunn. Undergr. Sp. Technol. 2012, 32, 158–167. [Google Scholar] [CrossRef]
  49. Wang, C.Y.; Chen, W.P.; Wang, L.P. Temperature beneath Tibet. Earth Planet. Sci. Lett. 2014, 375, 326–337. [Google Scholar] [CrossRef]
  50. Ford, D.C.; Williams, P.W. Karst Geomorphology and Hydrology; Chapman and Hall: New York, NY, USA, 1989. [Google Scholar]
Figure 1. (a) Distribution map of provinces in China and Study area. (b) Three altitude grades of topography in the study area ranging from 200 m to over 5000 m. (c) Water sample location and tunnel distribution.
Figure 1. (a) Distribution map of provinces in China and Study area. (b) Three altitude grades of topography in the study area ranging from 200 m to over 5000 m. (c) Water sample location and tunnel distribution.
Water 15 00207 g001
Figure 2. (a) Relationship between the number of Karst peaks and Karst development stages. (b) Relationship between the height standards and Karst development stages.
Figure 2. (a) Relationship between the number of Karst peaks and Karst development stages. (b) Relationship between the height standards and Karst development stages.
Water 15 00207 g002
Figure 3. The morphology of Karst depressions for different Karst stages 1, 2, and 3 represent three Karst development stages from strong to weak degree.
Figure 3. The morphology of Karst depressions for different Karst stages 1, 2, and 3 represent three Karst development stages from strong to weak degree.
Water 15 00207 g003
Figure 4. The morphology of contour for different Karst stages 1, 2, and 3 represent three Karst development stages from strong to weak degree.
Figure 4. The morphology of contour for different Karst stages 1, 2, and 3 represent three Karst development stages from strong to weak degree.
Water 15 00207 g004
Figure 5. The typical hydrogeological system for three topography grades and tunnel layout. (a-1). Layout of the Pingguo tunnel and the hydrogeology systems located in the first topography area; (a-2). The polygon distribution jointing the Karst peaks around Karst depressions. (b-1). Layout of the Dongfeng tunnel and the hydrogeology systems in the second topography area; (c). Layout of the Batang tunnel and the hydrogeologic system in the third topography area. (a-2,b-2) are used to obtain the fractal dimension representing the Karst development degree.
Figure 5. The typical hydrogeological system for three topography grades and tunnel layout. (a-1). Layout of the Pingguo tunnel and the hydrogeology systems located in the first topography area; (a-2). The polygon distribution jointing the Karst peaks around Karst depressions. (b-1). Layout of the Dongfeng tunnel and the hydrogeology systems in the second topography area; (c). Layout of the Batang tunnel and the hydrogeologic system in the third topography area. (a-2,b-2) are used to obtain the fractal dimension representing the Karst development degree.
Water 15 00207 g005
Figure 6. Lithology and water-bearing characteristics of study areas traversed by three typical railroad tunnels.
Figure 6. Lithology and water-bearing characteristics of study areas traversed by three typical railroad tunnels.
Water 15 00207 g006
Figure 7. SI- CaCO 3 of calcium carbonate in samples. ① SI ranges from 0 to 0.25; ② SI ranges from −0.25 to 0; ③ SI ranges from −0.5 to−0.25; ④ SI ranges from −0.5 to −0.75.
Figure 7. SI- CaCO 3 of calcium carbonate in samples. ① SI ranges from 0 to 0.25; ② SI ranges from −0.25 to 0; ③ SI ranges from −0.5 to−0.25; ④ SI ranges from −0.5 to −0.75.
Water 15 00207 g007
Figure 8. Relationship between δ18O (V−SMOW, permill) and δ2H (V−SMOW, permill).
Figure 8. Relationship between δ18O (V−SMOW, permill) and δ2H (V−SMOW, permill).
Water 15 00207 g008
Figure 9. Elevation difference between supply area and appearing area of water samples.
Figure 9. Elevation difference between supply area and appearing area of water samples.
Water 15 00207 g009
Figure 10. The altitude distribution and the contours in three typical study areas (the data come from the SRTMDEM 90M given by Geospatial Data Cloud, and these data are used to obtain the fractal Dimensions). (a-1,a-2) are for Pingguo System, which is divided into the strong groundwater bearing area and mild groundwater bearing area. (b-1,b-2) are for Dongfeng System, which is divided into deeply buried groundwater and shallowly buried groundwater. (c-1,c-2) are for Batang system, which is divided into a strong water-bearing stratum (C + D), a moderate water-bearing stratum (S), and the weak water-bearing stratum (P1-2).
Figure 10. The altitude distribution and the contours in three typical study areas (the data come from the SRTMDEM 90M given by Geospatial Data Cloud, and these data are used to obtain the fractal Dimensions). (a-1,a-2) are for Pingguo System, which is divided into the strong groundwater bearing area and mild groundwater bearing area. (b-1,b-2) are for Dongfeng System, which is divided into deeply buried groundwater and shallowly buried groundwater. (c-1,c-2) are for Batang system, which is divided into a strong water-bearing stratum (C + D), a moderate water-bearing stratum (S), and the weak water-bearing stratum (P1-2).
Water 15 00207 g010
Figure 11. Fractal relationship between area and circumference of polygon connecting the Karst peaks and saddle zones around Karst depressions (shown in Figure 5(a-2,b-2)). (a) Pingguo System; (b) Dongfeng System.
Figure 11. Fractal relationship between area and circumference of polygon connecting the Karst peaks and saddle zones around Karst depressions (shown in Figure 5(a-2,b-2)). (a) Pingguo System; (b) Dongfeng System.
Water 15 00207 g011
Figure 12. (a) Relationship between the height standards and Karst development stages (From Figure 2b). (b) Karst development phrase of the Pingguo System and Dongfeng System.
Figure 12. (a) Relationship between the height standards and Karst development stages (From Figure 2b). (b) Karst development phrase of the Pingguo System and Dongfeng System.
Water 15 00207 g012
Figure 13. Relationship between lnδ and ln(N(δ)). (a) Dongfeng System; (b) Batang system; (c) Pingguo System.
Figure 13. Relationship between lnδ and ln(N(δ)). (a) Dongfeng System; (b) Batang system; (c) Pingguo System.
Water 15 00207 g013
Figure 14. Relationship between tunnel, terrain and groundwater level in three grades. (a) Third grade; (b) second grade; (c) first grade.
Figure 14. Relationship between tunnel, terrain and groundwater level in three grades. (a) Third grade; (b) second grade; (c) first grade.
Water 15 00207 g014
Figure 15. Influence pattern of tunnel excavation on groundwater environment of three topography stages in China.
Figure 15. Influence pattern of tunnel excavation on groundwater environment of three topography stages in China.
Water 15 00207 g015
Table 1. Chemical and isotopic compositions of water samples in the study area where railways pass.
Table 1. Chemical and isotopic compositions of water samples in the study area where railways pass.
Sample
No.
Sample
Type a–e
Sampling Site fTDS
(mg/L)
pHT
(°C)
Na+K+Mg2+Ca2+ClSO42−HCO3Charge
Balance (%)
δ18O
(‰)
δ2H
(‰)
Exposed
Altitude (m)
Supply
Altitude (m)
Volume
(L/s)
Hydrochemical
Type
(mg/L)
1S1Pingguo Tunnel-----------−7.60−50.6112531-
2S26066.523.116.21.454.8614017.02128.722980.73−7.60−50.6112531-HCO3·SO4-Ca
3Grw33627.121.212.01.227.7871.311.3428.822290.94−7.74−64.5118477-HCO3-Ca
4Grw43887.221.012.30.853.8981.89.9817.62620.28−7.81−49.9243437-HCO3-Ca
5S53947.222.310.21.042.4388.211.3419.212620.71−7.67−49.2154506-HCO3-Ca
6Grw64127.120.55.920.391.4697.811.813.452810.62−7.32−53.6317553-HCO3-Ca
7Gr74167.222.36.350.744.8696.29.9319.212792.35−7.53−49.1154506-HCO3-Ca.Mg
8S8Xiaozai
Tunnel
3746.922.03.210.2214.5984.315.61.922549.04−9.95−76.912081492-HCO3-Ca
9S93067.122.03.350.319.7372.611.328.81798.69−9.68−72.612001388-HCO3-Ca
10S103567.320.425.32.734.8671.48.511.922419.81−10.3−78.712651630-HCO3-Ca.Na
11R111347.321.35.240.723.9528.98.9317.368.77.49−9.52−69.612641327-HCO3-Ca
12Gr12 Changqinpo
Tunnel
4637.120.532.12.117.380.27.0944.22901.51−8.57−7011001154637HCO3-Ca.Na
13S136527.022.155.24.3915.5697.89.9388.43811.92−10.6−68.81302183111.3HCO3.SO4-Ca.Na
14S149846.720.788.17.6139.412818.4244.04592.06−10.2−67.214771703-HCO3.SO4-Ca.Na.Mg
15W15 6997.322.459.24.2139.977.814.2196.03081.86−10.1−66.914761648-HCO3.Cl.SO4-Ca.Mg
16S16Xingfu
Tunnel
-----------−10.3−78.813901789-
17S17-----------−9.29−74.215121453-
18W184516.322.725.11.3223.877.070.91481051.12−10.2−67.115121746-HCO3.Cl.SO4-Ca.Na Mg
19Gr19 1086.524.213.21.127.308.027.5536.534.22.23---HCO3.SO4-Ca.Na
20R20 77.66.724.511.10.912.438.027.0911.536.63.03---HCO3.SO4-Ca.Na
21Gr21 1026.527.39.00.722.4316.07.3615.3751.31.65−11.3−70.518712112-HCO3.SO4-Ca.Na
22S224427.125.031.32.231.9584.29.9344.22681.53−10.4−70.216431806-HCO3-Ca.Na
23S233647.122.337.23.028.2748.98.2938.42201.96−10.0−70.014591695-HCO3-Ca.Na
24S24-----------−10.4−73.014791826-
25S25-----------−9.78−71.415901618-
26Gr26Xinshao
Tunnel
5747.123.14.120.2923.810512.01.924270.24−10.3−78.7125217802400HCO3-Ca.Mg
27S275177.122.48.100.6516.110117.013.53610.52−9.50−72.7124516391050HCO3-Ca
28S285217.120.31.210.0726.892.217.015.43690.30−9.86−75.51245170528.3HCO3-Ca.Mg
29S293977.019.72.100.0319.572.114.232.72560.33--22.0HCO3-Ca
30S30Dongfeng
Tunnel
4507.120.221.31.724.8688.214.91.973170.95−10.5−81.013971992330HCO3-Ca.Mg
31S314417.120.98.110.663.4099.414.72.023120.36−10.4−81.313971947240HCO3-Ca.Na
32S32-----------−10.5−81.213981963-
33S334327.020.45.310.344.8697.814.61.893080.27−10.2−69.51467187232.2HCO3-Ca.Na
34S344937.021.314.01.233.4010711.342.293540.68−9.48−67.21460158712.8HCO3-Ca
35S35----------0.00−10.2−72.914241861-
36S363947.020.51.110.035.3591.414.65.762760.06−10.6−84.414152020-HCO3-Ca
37Ws373927.142.012.01.093.4084.214.33.842730.77−10.9−86.414122116-HCO3-Ca
38Ws384807.147.010.30.828.2710015.14.043420.36−10.6−74.1143520271440HCO3-Ca
39S395237.122.45.210.449.2411621.317.33540.36--1.60HCO3-Ca
40S404767.122.76.210.706.3210614.91.923410.44−10.8−75.5160920792560HCO3-Ca
41S414357.019.54.230.237.3096.213.41.923120.54−9.44−70.41561157139.6HCO3-Ca
42S424087.220.913.01.1311.276.217.03.842860.88−10.8−78.915222106268HCO3-Ca
43S433717.221.312.11.012.4380.29.939.612560.81−10.8−82.7159021080.30HCO3-Ca
44S44Shiling
Tunnel
4547.121.011.40.974.3898.68.5113.53170.66−10.6−80.9173020021300HCO3-Ca
45S454557.022.39.230.794.8610015.05.763200.37−10.0−81.216911795100HCO3-Ca
46S46Laoshishan
Tunnel
4317.122.72.200.125.351008.5126.92880.22- 34.19HCO3-Ca
47S474037.220.211.20.957.3084.214.236.52490.74−10.8−79170817426.00HCO3-Ca
48S485477.720.59.140.7312.211411.338.43610.52---HCO3-Ca
49S493527.319.12.010.0715.664.18.517.682540.34−12.2−89.21852228014.6HCO3-Ca
50S502397.319.07.320.336.3245.77.915.761660.98−11.6−83.117882042278HCO3-Ca
51S515137.319.011.30.929.241105.6778.82980.60---HCO3-Ca
52S52Batang tunnel341 7.811.06.622.3112.253.41.1521.0244−5.15−17.7−13026794892-HCO3-Ca.Mg
53S53334 7.410.00.530.8824.640.70.9115.4251−3.82−17.0−12827864312-HCO3-Ca.Mg
54S543137.711.33.530.8719.048.90.9032.92071.08−17.4−13129034492-HCO3-Ca.Mg
55S55461 7.111.118.21.168.4976.81.4329.6325−5.49−17.6−13039364558-HCO3-Ca.Mg
56S56201 8.012.50.040.1111.031.10.687.59151−3.52−17.1−13133634369-HCO3-Ca.Mg
57S57227 8.012.20.110.2012.528.50.847.19178−10.88−18.1−13529844726-HCO3-Ca.Mg
58R58110 8.013.02.120.124.0123.60.519.1970.98.12−18.0−13330644704-HCO3-Ca.Mg
59S59338 7.515.02.850.723.8177.61.0538.32140.02---HCO3-Ca.Mg
60S60288 7.512.01.060.599.5257.81.0430.51870.20---HCO3.SO4-Ca
61S61283 7.514.32.170.667.1159.51.0332.31800.35---HCO3.SO4-Ca
62S62183 7.512.92.120.192.7941.21.008.911271.89---HCO3-Ca.Mg
63S63263 7.514.50.010.567.8756.10.9935.51620.73---HCO3.SO4-Ca
64S64375 7.513.26.201.1212.079.83.5483.31893.44---HCO3.SO4-Ca
Notes: a. S is the sample from springs; b. Grw is the groundwater sample obtained through the geological window; c. Gr is the groundwater sample from ground river; d. R is the water sample from surface river; e. Ws is the spring sample from warm spring; f. Water samples are classified based on nearby railway tunnels.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tang, Y.; Zhang, Q.; Qi, J.; Xu, M.; Li, X.; Qu, C.; Yi, L.; Wang, D. Change Analysis of Karst Landforms, Hydrogeological Conditions and Effects of Tunnel Excavation on Groundwater Environment in Three Topography Grades in China. Water 2023, 15, 207. https://doi.org/10.3390/w15010207

AMA Style

Tang Y, Zhang Q, Qi J, Xu M, Li X, Qu C, Yi L, Wang D. Change Analysis of Karst Landforms, Hydrogeological Conditions and Effects of Tunnel Excavation on Groundwater Environment in Three Topography Grades in China. Water. 2023; 15(1):207. https://doi.org/10.3390/w15010207

Chicago/Turabian Style

Tang, Yige, Qiang Zhang, Jihong Qi, Mo Xu, Xiao Li, Chenhao Qu, Lei Yi, and Dong Wang. 2023. "Change Analysis of Karst Landforms, Hydrogeological Conditions and Effects of Tunnel Excavation on Groundwater Environment in Three Topography Grades in China" Water 15, no. 1: 207. https://doi.org/10.3390/w15010207

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