• Included in CSCD
  • Chinese Core Journals
  • Included in WJCI Report
  • Included in Scopus, CA, DOAJ, EBSCO, JST
  • The Key Magazine of China Technology

2022 Vol. 41, No. 2

Display Method:
Low Impact Development construction of peak cluster depression based on regulation of epikarst zone
JIANG Guanghui, LIU Fan, WANG Qigang, GUO Fang
2022, 41(2): 165-173. doi: 10.11932/karst2021y32
Abstract:
Most of the bare karst areas environments with the widest distribution in southern China have peak cluster depression landforms. With high incidence of many environmental problems such as rocky desertification, flood, water shortage, and water pollution, the karst environment there is proved to be very fragile. For the bare karst environment, it has a very developed subterranean river drainage networks, causing most of the water resources to enter the aquifer, and then flow and discharge rapidly through the underground karst conduits. For a long time, as a karst living environment, the lack of water and soil resources have become the main cause of environmental degradation.The introduction of the concept of Low Impact Development (LID) is necessary and timely in the karst peak-cluster depression area, based the rich research and practice in southern China. It is believed to be effective for promoting the utilization of water resources, reducing human disturbance, and covering the other environmental problems, such as rocky desertification and risk of underground river pollution. The LID methods based on the regulation of the epikarst zone and sinkholes are mainly suitable for guiding the development of karst land in dolines and depressions. It is realized by solving the questions of how to plan the land use and to design the rural infrastructure according to the special hydrogeological conditions and epikarst hydrological process. The concept and technology of LID in view of urban facilities will give guidance for the construction of a holistic environmental remediation project in a rural area, including water resources storage, pollutant treatment, and land restoration for rocky desertification. It is very helpful to promote rural revitalization in the remote and fragile ecological environment area of karst rocky mountains. Karst dolines naturally have topographic and hydrogeological conditions that promote rapid rainfall infiltration. Controlled by monsoon climate, the karst environment is characterized by distinct rainy season and dry season. Due to the lack of soil cover, storm water forms various types of runoff on the rocky slope, which flows into the underground river eventually through epikarst zone and sinkholes. The epikarst zone with a high degree of karstification has high hydraulic conductivity and regulation and storage function with a large number of small caves. It is the main water storage space on the slope and forms the epikarst spring. The bottom of the depression is flat with continuous soil cover. The flood is discharged by the sinkhole, but sometimes the groundwater inversely enters the karst depression through the sinkhole, resulting in waterlogging. Karst depressions are widely used as living space and agricultural planting. According to its terrains, development practice and disturbance intensity, the land of karst doline can be divided into three kinds of sensitive areas, and they are, (1) sloping farmland with risk of soil erosion, (2) epikarst springs for drinking water source with risk of pollution, and (3) the sinkholes with risk of contaminant input. The LID model would try to achieve five management goals in the three sensitive areas.(1) The goal of LID in the sloping farmland is to minimize water consumption by crop irrigation and to protect precious and thin soil from erosion, which technically can be achieved by requiring a filtering system such as building terraced fields, interception basin and straw mulching, so as to regulate soil moisture and slope flood runoff. (2) The management goals of the sensitive areas of epikarst springs are the utilization of rainfall-runoff resources and the ensurence of drinking water safety, which technically can be achieved by implementing rainwater harvesting projects and setting up water source protection area. (3) Contaminant interception and absorption facilities, such as buffer zone and sedimentation and filtration ditch, should be built in the sensitive area of sinkhole for achieving the goal of stopping contaminants entering groundwater. The LID model for peak cluster, or “LID based on epikarst”, which is characterized by rainfall -runoff management in the epikarst zone, has been explored and summarized in the karst rocky desertification region of South China. The model requires reserving the water resources in the woodland area at the top of the peak in the process of development planning and utilization of karst depression, for enhancing the recharge by infiltration through epikarst zone. In the intensive utilization of land at the bottom of karst dolines, the buffer zone would be reserved for controlling flood discharge by sinkholes, and a low elevation green belt is set with filled karst grooves for promoting natural infiltration and purification. In the transition zone between “peak” and “doline”, the scope and boundary conditions of water confluence at epikarst spring should be clearly defined, and the boundary of the water source protection zone should be determined. At the same time, it is necessary to increase the utilization of rainwater to supplement the water shortage of ecosystems and living. Although LID model based on epikarst is raised from the rocky desertification control in southern China, it could be possibly extend to the karst areas in other regions such as Southeast Asia, where peak-cluster is also well developed.
Niangziguan spring group in Shanxi Province and its water source
TANG Chunlei, LIANG Yongping, JIN Hua, SHEN Haoyong, ZHAO Chunhong, WANG Zhiheng, XIE Hao, ZHAO Yi, WANG Shina
2022, 41(2): 174-182. doi: 10.11932/karst20220201
Abstract:
As a main natural discharge point of karst groundwater, Niangziguan spring group consists of 11 springs and is located in Niangziguan town, Yangquan City, Shanxi Province. It is distributed on the floodplain and terrace of the Taohe river, the Wenhe river and the Mianhe river, with an exposed elevation of 360 m to 392 m. Among these springs, Chengjia spring is the highest and Jinqu spring is the lowest. The annual average flow of spring group is 9.48 m3·s−1 and the flow differences of springs are large, among which the difference of Wulong spring is the largest. Niangziguan spring group is the largest karst spring in northern China, with a catchment area of 7,436 km2. Previous studies believe that karst water in the spring area flows from north, west and south to Niangziguan area, and the chert mass or strip dolomite of lower Ordovicianin Niangziguan area relatively rises to the aquifuge, and is eroded and exposed by the Taohe River, which makes the karst groundwater overflow to the surface and form a spring group. The main aquifer is made up of gypsum-bearing carbonate rock of middle Ordovician. However, because Niangziguan springs are not from the same source, their hydrochemistry and isotope are different. Therefore, it is necessary to study the source of Niangziguan spring group. This study is conducted with the methods of hydrochemistry, isotope and hydrogeological section, and is concluded that there are two aquifers and three subsystems in Niangziguan spring area. The aquifer consists of middle Ordovician limestone aquifer and middle upper Cambrian dolomite one. Subsystem distribution is comprised of Ordovician karst water system in the west, and Ordovician karst water system and middle upper Cambrian karst water system both in the East. Chengxi spring and Chengjia spring are exposed on the marl of lower Majiagou formation of middle Ordovician, and their aquifers consist of fissure and fissure water of Middle Ordovician limestone. Because the marl of lower Majiagou formation is uplifted in this region, the surface spring-the eastern Ordovician karst water system-is formed. Mainly supplied by the water-bearing formation of middle upper Cambrian, Podi spring, Wulong spring, Shibanmo spring, Shuiliandong spring and Weizeguan spring are pressure-bearing ascending springs and belong to the middle upper Cambrian karst water system in the east.
Influence process of thermal structure variations of a karst water reservoir on dissolved inorganic carbon and its stable carbon isotope
LI Dashuai, WU Shaojing, LI Jianhong, ZHANG Tao
2022, 41(2): 183-196. doi: 10.11932/karst2022y02
Abstract:
Dalongdong Reservoir, a typical groundwater recharge reservoir in karst area of southwest China, is located at 23°30′01″−23°40′08″N and 108°30′02″−108°36′04″E in Shanglin county, Nanning City, Guangxi. This large reservoir was built mainly for irrigation and power generation by blocking caves, sinkholes and some fissures in a karst valley in 1958. It covers a rainwater collecting area of 310 km2, totaling storage capacity of 151 million m3. Approximately 95% of the recharge water comes from two karst subterranean streams. The principal strata in the reservoir catchment are carbonate rocks of Carboniferous (C1) and Devonian (D2d3) era. Influenced by the subtropical monsoon climate, the average temperature in the reservoir area is 21℃ and the average rainfall is 1,837.3 mm. The rainy season in this area is from April to September. Samples were collected respectively by the plankton method at the upper, middle and downstream monitoring sites in January, April, July, October and December, 2018. Based on DIC concentration, δ13CDIC value, water temperature, DO, SIC, partial pressure of carbon dioxide, carbon dioxide diffusion flux and other indicators, the thermal structure change on DIC concentrations and isotope composition δ13CDIC in each water layer of reservoir and its influencing factors are systematically studied. The gradient differences of DIC concentrations and isotopic composition δ13CDIC in different water layers and their causes are also discussed.In order to provide scientific and technological support for accurate assessment of carbon budget of land and water, the carbon cycle process of karst reservoirs was deeply revealed. The results show that Dalongdong reservoir presents periodic thermal structure changes in the mixing stage−stratification stage−in a hydrological year. The thermal stratification begins to appear in April, gradually presents complete thermal stratification in July, and gradually disappears after October.Besides, thermal stratification is an important driving force of DIC concentration and δ13CDIC value. After the formation of thermal stratification, DIC concentrations in the upper, middle and lower reaches show the same characteristics as those of the δ13CDIC values at the same layer, and DIC concentrations increase with depth, while δ13CDIC values become negative. DIC concentrations and δ13CDIC values are respectively 3.22 mmol·L−1 and −9.15‰ in surface water, 3.43 mmol·L−1 and −9.70‰ in thermocline, and 4.32 mmol·L−1 and −11.89‰ in bottom water. After thermal stratification, the vertical, horizontal and seasonal variations of DIC concentrations and δ13CDIC values gradually disappear, and finally show their homogenization.Finally, carbon dioxide degassing and photosynthesis dominate the water-air interface in the surface water layer. Degradation of organic matter is the main process in thermocline, and carbonate precipitation can be mostly found in the bottom water layer. The change of these three processes with the reservoir thermal structure may drive the variation of DIC concentrations and isotope gradient G (DIC) and G13CDIC), which shows the characteristics of bottom water layer < surface water layer < thermocline layer.
Metallogenic regularity of strontium-rich drinking natural mineral water in Devonian carbonate area in southern Hunan: Taking strontium-rich drinking mineral water in Xinxu town, Xintian county as an example
ZHOU Xin, WANG Can, ZHENG Pengfei, YAO Tengfei, WU Zhengqing, MI Maosheng, QIN Zuohui, LI Yang
2022, 41(2): 197-209. doi: 10.11932/karst20220202
Abstract:
The study area is located in Xinxu, Xintian county, the middle of southern Hunan, in which Devonian carbonate rocks are widely distributed in the south of N26°30'. 44 samples have been taken from some boreholes and springs in different lithologic distribution areas in different seasons for many years. The sampling sites are located in the marl distribution area of Shetianqiao formation of Devonian interplatform basin (D3s). Testing statistics of these samples show that the strontium content of groundwater in ZK13-10, ZK1601 and ZK1602 boreholes is 1.66 mg·L−1, 6.76 mg·L−1 and 0.41 mg·L−1 respectively, which can prove the existence of strontium-rich drinking natural mineral water.Based on the analysis of lithology, structure, hydrogeological conditions, hydrogeochemical origin and influencing factors of strontium enrichment, the metallogenic regularity of strontium-rich mineral spring water in Devonian carbonate distribution area is preliminarily revealed. According to the statistics of rock and mineral analysis and test results of 80 groups of rock cores at different depths of 5 boreholes, the average content of strontium in marl of Shetianqiao formation of Devonian interplatform basin (D3s) is 524×10−6, peaking at 800×0−6, much higher than that in other rocks.The high abundance of strontium is the material basis for the formation of strontium-rich mineral water in the study area. Primarily, the content variety of strontium in groundwater is closely related to groundwater circulation conditions. Three spots of strontium-rich groundwater in the study area mainly come from strontium-rich water in the zone of deep subsurface flow. High circulation depth, relatively closed groundwater circulation and weak hydrodynamic conditions are the main constraints for the formation of strontium-rich mineral spring. Because the groundwater recharge runoff is restricted by topographic and geomorphic conditions, seepage conditions, fault structure location and water control properties, the farther the distance of groundwater recharge runoff is, the higher the circulation depth will be, and in turn the more conducive will be to the enrichment of strontium and the formation of deep strontium-rich water. The water temperature is another important factor affecting the enrichment of strontium in groundwater because it affects the dissolution and ion replacement of minerals containing strontium, and it also affects the CO2 content and corrosiveness in groundwater, thus affecting the enrichment of strontium in groundwater. Meanwhile, the circulation depth and geothermal field will control the groundwater temperature.The metallogenic regularity of Xinxu strontium-rich drinking mineral water shows that the marl distribution area of Shetianqiao formation of Devonian interplatform basin (D3s) is a favorable area for people to look for the strontium-rich mineral spring in southern Hunan. In particular, the hanging wall of water-blocking faults in gentle karst hills as well as the location close to magmatic rocks and regional thermal conductivity faults are the target areas.
Evaluation of rocky desertification comprehensive control based on orchard intercropped by herbage: A case study of the Xingyi City, Guizhou Province
CHENG Hanting, LIU Guodao, DONG Rongshu, WANG Xiaomin, LI Qinfen, ZHOU Xiaohui, ZHANG Xianbo
2022, 41(2): 210-219. doi: 10.11932/karst20220203
Abstract:
An evaluation of rocky desertification comprehensive control is very important to optimize and improve strategy for rocky desertification control modes. This paper aims to (Ⅰ) analyze the relationship among agricultural structure, industrial resources situation, economy benefit and ecology benefit; (Ⅱ) evaluate the comprehensive benefit of karst rocky desertification control by orchard-herbage intercropping, and further find its own advantages and disadvantages; (Ⅲ) provide ecological economic and sustainable development reference for rocky desertification control in the southwest karst area of China.A comprehensive prevention and control of rocky desertification area has been implemented in combination with tropical fruit tree intercropped by legume herb in Tianfang village, Nanpanjiang town, Xingyi City of Guizhou Province,where the rocky desertification is serious. To evaluate the comprehensive benefit of the demonstration area, we investigated the agricultural structure, industrial resources situation, and the benefits of economy and ecology before (2014) and after (2017 and 2020) karst rocky desertification control, and then adopted correlation analysis, Redundancy Analysis (RDA) and Structure Equation Modeling (SEM).By converting cultivated land to orchard (mango and Macadamia intercropped with herbage or green manure) the agricultural industrial structure ,industrial resources development,economic benefit and ecological benefit of the demonstration area showed a good trend after 6 years of rocky desertification comprehensive control. Correlation analysis showed that the agricultural industrial structure and industrial resources development situation were significantly correlated with economic benefit and ecological benefit, except for farming/forest-grass area ratio and commodity output value of agricultural products. RDA analysis showed that there are five variables which are the important contributing factors to the changes of economic benefit and ecological benefit such as, proportion of orchard land, migrant worker population, economic contribution rate of migrant workers, commodity output value of agricultural products farming/forest-grass area ratio SEM further showed that the rocky desertification control are able to couple the coordinated development of ecological environment and industrial resources, and at the same time achieve ecological restoration and economic benefit.Conclusions are drawn as follows, (Ⅰ) The rocky desertification control in orchards intercropped by fruit tree has increased the number of labor force and migrant workers,and improved the economic benefits by increasing commodity output value of agricultural products and migrant workers income. (Ⅱ) The ecological environment has been improved by increasing forest and grass coverage,soil erosion and water loss and surface runoff loss in the comprehensive control area of rocky desertification has been reduced. (Ⅲ) The rocky desertification control by orchard covered with fruit tree and intercropped by crops controlling rocky desertification will improve both economic benefit and ecological benefit. Overall, rocky desertification control by fruit tree intercropping orchard is an effective way to achieve sustainable development in the southwest karst area of China The study results provide a basis for the comprehensive control of karst rocky desertification.
Comparison of soil calcium and magnesium fractions transport in classic karst and non-karst region, Guilin
LIANG Jianhong, CUI Xudong, WEN Laiyan, LIU Ding, YI Chenxu, HUANG Kezun, WANG Jun
2022, 41(2): 220-227. doi: 10.11932/karst20220204
Abstract:
The exposed karst area in the world is about 17.9 million km2, mainly distributed in China, Vietnam, France, Slovenia, the United States and other countries. Covering a total area of 780,000 km2, the karst are in south China is most widely distributed in the world, mainly in the most area of 8 provinces/cities such as Yunnan, Guizhou, Guangxi, Hunan, Hubei, Chongqing, Sichuan and Guangdong. The southern karst area represented by Guilin is a typical karst area in China, where soil is formed by the physical, chemical and biological effects of carbonate rocks in the geological background. Due to the richness of calcium (Ca) and magnesium (Mg) in carbonate rocks, alkaline soil rich in Ca and Mg is distributed in karst areas. But non-karst areas are primarily clastic rock, where soil mainly exists acidic. Besides, Ca in soil is an essential mineral nutrient for plant growth and metabolism. Soil Ca and its morphological changes affect the physicochemical properties of soil, and also affect the uptake of calcium and other nutrients by plants. Magnesium is also an essential mineral element for plant growth, ranking the fourth after nitrogen, phosphorus and potassium. Therefore, plants need to take in Mg, depending on the content and supply mechanism of Mg in soil.With a modified Tessier sequential extraction procedure, soil profiles from typical karst and non-karst areas in Guilin are analyzed to determine different fractions of Ca and Mg,the ion-exchange fraction (including water-soluble state), carbonate-bound fraction, humic acid-bound fraction, Fe-Mn oxides bound fraction, strong organic matter-bound fraction (including some sulfides) and residual fraction. The migration characteristics of Ca and Mg elements in soil from these two types of areas are explored. The vegetation type of sampling sites is mainly forest. During the sampling, a soil profile with 1 m wide, 2 m long and 1.5 m deep was dug. The profile is stratified from the surface to 140 cm, and divided into 7 horizons at intervals of 20 cm. Then soil samples are collected from the bottom layer by layer. A total of 7 soil samples in the limestone soil profile are collected from the depth of 105 cm at intervals of 15 cm. The results showed that in the limestone soil profile of the karst area, the pH value increases with the depth of the profile, while the contents of soil organic matter, total nitrogen (N) and total phosphorus (P) decrease. Except the parent material, the P content of soil is related to other factors in the process of soil formation, among which bioaccumulation and soil formation are important. The weak migration of P, together with the limited availability of limestone soil makes the P content difficulty to leach down from the profile, hence the P content of the surface layer is higher than that of the lower layer. Calcium mainly exists in the exchange fraction and its content is similar in each soil profile, while Mg mainly exists in the residual fraction. There is an accumulation when the percentage of humic acid-bound fraction decreases with the depth of the soil profile. The maximum value of the ion-exchange fraction of Ca occurs in the middle profile at the depth of 45-60 cm. In the karst area, Mg is mainly present in the residual fraction (97.45%) in the limestone soil profile. But Mg content is low in the ion-exchange fraction including the water-soluble fraction (0.61%), the carbonate-bound fraction (0.11%) and the humic acid-bound fraction , accounting for 1.80% of the total amount in the latter three fractions. Compared with the typical karst area, the total Ca content is lower in acidic soil of the non-karst area. The soil profile from 0 cm to 60 cm is dominated by strong organic bound fraction, and the deeper profile (>60 cm) by residual fraction. Calcium in the topsoil in non-karst areas is mainly present in the ion-exchange fraction (including the water-soluble fraction), the Fe-Mn oxide-bound fraction and the strong organic-bound fraction (including some sulfide). The dominance of the ion-exchange fraction (including the water-soluble fraction) indicates the relatively high mobility of Ca. In non-karst areas. The effective content Mg is low, mainly in humic acid-bound fraction and residual fraction. Similar to karst areas. Mg in non-karst areas also accumulates in the soil to some extent. In conclusion, the study explores the morphological migration of Ca and Mg in soil profiles in karst and non-karst areas, and provides scientific support for soil improvement and risk assessment of Ca and Mg leaching in soil.
Spatial variation analysis of soil organic matter in karst area
ZHANG Chunlai, LU Laimou, YANG Hui, HUANG Fen
2022, 41(2): 228-239. doi: 10.11932/karst20220205
Abstract:
As an important parameter for soil quality assessment, soil organic matter (SOM) can provide nutrients for crops, strengthen soil fertility, and improve soil physicochemical properties. Besides, SOM also plays an important role in restricting the activity of heavy metal elements and pesticide residues, and in regulating CO2 in soil. Therefore, the study on spatial distribution of SOM, its influence factors and prediction by GIS and geostatistical analysis is important in agricultural practice, environmental management and the measurement of soil carbon storage.The content of SOM is easily affected by the factors such as geological background, climate change, agricultural activities and land use change. Hence the complexity and non-stationarity of the spatial variability of SOM will make it difficult for the quantitative prediction. The high heterogeneity of SOM and insufficient sampling data in the karst area also leads to low accuracy of spatial prediction. Guangxi is one of the areas with the most widely distributed karst areas in China. Typical karst landforms are developed in the northern part of Mashan county, with a range of land use and parent materials. Different patterns of land use include paddy field, dry land, shrub land, forest land, orchard and grassland, accounting for 10.55%, 29.39%, 29.49%, 12.86%, 0.70% and 0.66% respectively. Paddy fields, consisting of paddy soil, are located in karst valleys; dry lands, consisting of red soil, brown lime soil, and alluvial soil, are distributed at the bottom of slopes and karst depressions on both sides of the valley; shrubland and forest land, mainly consisting of brown lime soil and red soil, spread over mountains. A total of 441 SOM data, with 8 high-density topsoil (0-20 cm) samples per km2 on average, were obtained in the geochemical survey of land quality in the northern karst area of Mashan county, Guangxi in 2017. The data facilitates the spatial prediction of highly heterogeneous SOM in karst areas.In this paper, qualitative variables such as land uses and soil types were converted into quantitative variables through the assignment method of dummy variables, and terrain factors were used as auxiliary variables to meet the requirement of geostatistical regression analysis. With the purpose of exploring the applicability of geostatistical SOM mapping in karst areas, this study respectively established five method models—ordinary Kriging (OK), regression Kriging (RK), Geographically Weighted Regression Kriging (GWRK), Ordinary Kriging of Median Centralization (MC_OK) and Mean Modified Ordinary Kriging (MM_OK)—to compare the prediction results, and select the most suitable prediction model for the karst area.. Results show that the SOM content in the study area ranges from 0.81% to 5.03%, with an average of 2.36%, and the coefficient of variation is 37.3%. The spatial distribution is moderate spatial variation. Results of geostatistical analysis by GS+9.0 show that the fitted theoretical variogram models are all exponential ones. Both nugget value and abutment value are less than 25%, which indicates that there is a certain spatial autocorrelation, and the SOM content is mainly affected by structural factors; hence correlation interpolation analysis can be carried out.In addition, the spatial variation of SOM in the karst area is jointly affected by land use, soil types and topographical factors. The areas with high content of SOM are located respectively in the karst areas and paddy fields where lime soil is distributed in the northwest, west and east of the study area, and the areas with low content are located in the alluvial land along the Hongshui river in the north of the study area. The SOM content shows significantly positive correlation with paddy fields because the reducing environment of paddy fields is conducive to the preservation and accumulation of SOM. But SOM content shows significantly negative correlation with dry land due to the oxidative environment of this land type where frequent disturbance and rapid degradation of SOM will occur.Finally, RK, GWRK, MM_OK and MC_OK can be used for SOM prediction and mapping in karst areas because of their good interpretability. The results of internal verification show that the root-mean-square of GWRK is the smallest, and the average standard error is close to root-mean-square, indicating the better fitness of GWRK model. The results of external verification demonstrate that the mean absolute error and root mean squares error of GWRK are the smallest, suggesting the highest precision and the best stability of GWRK model. Meanwhile, MC_OK shows its best accuracy. Combined with auxiliary variable factors such as land use, soil types, and topography, GWRK model can effectively eliminate the influence of spatial variation factors and overcome the spatial non-stationarity of SOM content in karst areas, thereby improving the stability and precision of the SOM prediction model, while MC_OK model can improve the prediction accuracy.
Laboratory simulation study on soil surface loss and underground leakage in the epikarst fissure zone
YAN Yiquan, LIU Qi, DENG Dapeng, WANG Han
2022, 41(2): 240-248. doi: 10.11932/karst20220206
Abstract:
The study area, Guanling-Zhenfeng Huajiang demonstration area of karst plateau canyon, covers 51.62 km2 with the karst area of 45.39 km2. It is a typical area of karst plateau canyon on the Guizhou plateau with an elevation difference between 370 and 1,355 meters. In this study, a simulation device was developed by researchers to investigate soil surface loss and underground leakage in the epikarst fissure zone. By simulating indoor rainfall and controlling the rainfall intensity, slope, fracture width and fracture occurrence, a 90-minute experiment was carried out to probe into the influence of soil surface loss and underground leakage. During the experiment, the data such as run-up time of underground leakage, water and soil surface loss and underground leakage was obtained every 4 minutes. EXCEL 2016 was used for statistical analysis of the data, and Origin 2017 for mapping.Results show that soil surface loss is mainly affected by rainfall intensity and slope. The lose increases with the intensifying of rainfall and the rising of slope, reaching the highest at 30° slope gradient. Besides, soil underground leakage is mainly affected by fissure width, occurrence and slope, but not obviously affected by rainfall intensity. The leakage is positively correlated with the fissure width, and most likely to occur when the fissure direction is 30° to the slope direction. The slope gradient is negatively correlated with the occurrence of soil underground leakage. During the rainfall, the underground leakage rate changes greatly, experiencing a decrease after an increase and finally reducing to 0. Surface loss is the main form of soil loss on the karst slope without sinkholes and funnels, and soil underground loss contributes less than 5% of the soil loss. Large karst pipelines such as sinkholes and funnels are the main channels of soil underground leakage. These experimental results can provide theoretical reference to the control of soil and water loss as well as the ecological restoration of karst sloping land in southwest China.
Change and its prediction of landscape patterns in the watershed of typical karst lake wetland
WEN Guojing, WANG Yan, LIU Yungen, HOU Lei
2022, 41(2): 249-258. doi: 10.11932/karst2021y41
Abstract:
As a research focus of landscape pattern change, the change and its prediction of landscape patterns in the watershed of karst lake is explored with the purpose of providing a scientific basis for the ecological safety in karst watershed in this study. The data of remote sensing image in the watershed of Puzhehei karst lake was collected in 1990, 1995, 2000, 2005, 2010, and 2015 respectively. With the remote sensing and geographic information system, those images were classified and interpreted based on the actual situation of this area. In addition, landscape patterns and the dynamic characteristics were analyzed and the pattern of future wetland landscape was simulated and predicted with CA-Markov model. Results show that landscape patterns of the watershed of Puzhehei karst wetland changed significantly from 1990 to 2015, during which construction land and agricultural land continued to increase, but woodland, wetland and unused land continued to decrease. The number of blocks (NP) rose from 861 to 889. Generally speaking, agricultural land is increasing year by year, but the growth rate is slow and the wetland landscape area is decreasing. Besides, from 1990 to 2015 the degree of landscape fragmentation in the Puzhehei watershed was generally augmented,and the number of landscape patches, patch density and number index of maximum patches also increased. The number of patches reached a maximum of 936, and the patch density index increased from 2.596 in 1990 to 2.822 in 2000, and then decreased to 2.689 in 2015. The number of maximum patches increased from 58.345 in 1990 to 62.036 in 2015, which means that the impact of human activities on the landscape pattern in the Puzhehei watershed made the landscape structure more complex and the fragmented. Finally, the dominant patch area of landscape types gradually increased from 1990 to 2015. The diversity index decreased from 1.064 to 0.966, and the edge density increased from 47.552 to 48.063. But the landscape shape index decreased from 30.175 to 29.625. Due to the continuous reduction of forest land and wetland in 2000, the landscape shape index increased and both the edge density and landscape shape index peaked at 49.987 and 31.012 respectively, which made the landscape structure complex. In the next 20 years, the landscape type change in Puzhehei watershed will generally show an increase in construction land, agricultural land and wetland, but a decrease in forest land, garden land and unused land. This trend can be reflected from the previous data: construction land increased from 2.79% to 2.97%, agricultural land from 60.12% to 60.74%, wetlands from 6.67% to 7.02%, but of forest land dropped from 26.70% to 26.40%. Accordingly, the land use pattern of Puzhehei watershed will change significantly. By 2030, the existing construction land will be expanded mainly around the scenic area of sinkholes as well as new in the villages and new towns around the National Wetland Park. Constrcution land will spread to the surrounding areas along with the main traffic roads. It is predicted that wetland, construction land and agricultural land will experience the largest change. The research findings can provide suggestions and reference for the ecological protection of the watershed Puzhehei karst lake.
Study and prospect of karst collapse columns and their water inrush in the coalfield of North China
XU Guangquan, ZHANG Haitao, ZHOU Jisheng, LI Xu, WANG Minhua, LIU Mancai
2022, 41(2): 259-275. doi: 10.11932/karst20220207
Abstract:
Karst collapse columns (KCCs) are formed in the process of geological evolution in North China coalfield. Their roots are developed in the Ordovician or Cambrian carbonate rocks, and pass upward through Carboniferous, Permian, Triassic strata, and even enter the Jurassic, Cretaceous and Quaternary loose rock strata. Over the past 80 years, more than 10,000 mines have been exposed in 39 mining areas in the coalfields of North China, resulting in more than 20 major water inrush accidents. Therefore, the study on KCCs is of important theoretical and practical significance for the systematic understanding of karst hydrogeological conditions of mines in North China and of prevention and control of karst water disaster as well.This study systematically summarized the spatial distribution characteristics of KCCs in the coalfields of North China, mainly focusing the areas along the Fenhe river, the western foothills of Taihang mountain, the foothills of eastern and southeastern Taihang mountain, the foothills of southern Yanshan mountain, and the west side of the Tanlu fault zone. The plane shapes of KCCs are mostly circular and elliptical, and their cross-sections are mainly cone-shaped. KCCs are usually filled and cemented by broken debris of wall rock and secondary minerals, so most of the collapsed columns in the coalfields do not conduct (contain) water, and only few cases of water inrush occur there. KCCs usually form in such basic conditions as soluble rocks and special stratigraphic structures, multi-stage and multi-source fluid erosiveness, geological tectonic evolution, and paleokarst groundwater flow. Five theories on the genesis mechanism have been put forward, namely “karst gravity collapse”, “gypsum dissolution collapse”, “vacuum absorption collapse”, “hydrothermal origin” and “groundwater internal circulation”.According to the current research on the formation and evolution of KCCs in the Huainan coalfield in the southern margin of North China, this study shows that consisting of Ordovician or Cambrian carbonate rocks, KCCs move through fluid migration channels of faults and fractures generated by multi-stage tectonic movements in this area, and are accommodated with the paleokarst formed by dissolution in the carbonate rocks. After a long-term process of dissolution, transportation, and collapse, KCCs are finally formed with the joint effect of the self-gravity of the rock mass, the in-situ stress, and the vacuum negative pressure in the cave. The formation of KCCs is related to the multi-stage tectonic movements in Mesozoic, and the Yanshan movement might be the key stage of its formation and development.In this study, the water inrush mechanism of KCCs in the mining process is reviewed from the aspects of mining-induced water inrush mechanism, seepage transformation mechanism, and water inrush model. Four modes of water inrush of KCCs in working face or roadway, as well as the mechanical mechanism and seepage transformation mechanism of water inrush are analyzed, and the reasons for the sudden and hysteresis of water inrush are also explained. The technical route for KCC detection and prevention are determined by geophysical prospecting, geochemical prospecting, roadway prospecting, drilling, and hydrogeological conditions.At present, the prevention and control of KCC water disaster is mainly carried out from two aspects,exploration and treatment. In terms of exploration, several methods can be used. The technology of three-dimensional seismic exploration can be adopted to preliminarily determine the impacted area of collapse columns. The comprehensive geophysical method can be used to detect the water-bearing area of KCCs. The exact location, shape and hydrogeological characteristics of KCCs can be determined by methods of drilling, the test of pumping (draining) water, geochemical exploration and others. In terms of treatment, two methods can be mainly adopted,water-proof coal (rock) pillar retention and grouting reinforcement/water blocking to cut off the hydraulic connection between the KCCs water channel and the water inrush source. This technology has been successfully applied in the treatment of KCCs in the mining areas subject to karst water disaster such as Hebei and Anhui.As to the research gap and deficiency in the prevention of water disaster caused by KCCs, the following suggestions are put forward,(1) Under the guidance of earth system science, the study on evolution of KCCs can be conducted combined with the formation mechanism of KCCs. Meanwhile, a multidisciplinary method involving petrology, structural geology, hydrology, geomorphology, karst hydrogeology, GIS, RS, and others can also be applied based on the time axis of evolution of regional geological structure.(2) For the precise detection, data and interpretation of small hidden collapse columns, efforts should be made to perfect the basic theory of full-space physics, to develop new equipment, new technologies and new methods, to further improve the precision, resolution, anti-interference competence of geophysical exploration, and to increase the detection depth.(3) Research on the KCCs model of “dynamic monitoring-prediction model-water inrush mechanism” should be conducted. Firstly, a dynamic monitoring system to obtain parameters such as stress, deformation, acoustic emission, rupture (wave velocity), water pressure (position), water chemistry, and water temperature during the water inrush process of KCCs can be established, and an indicator monitoring instrument of stable multi-parameter can be developed. Secondly, with the help of big data and data processing technology of cloud platform, a multi-parameter integrated water inrush prediction model of KCCs can be established based on the dynamic monitoring system characterized by information transmission and the existing water-rock coupling numerical model.(4) The prevention and control technology system of KCCs water disaster should be established and improved. This system will mainly include the technical subsystem such as basic theoretical research on KCCs formation, exploration of KCC basic characteristics, prevention, management and emergency planing of KCC water disaster. In short, the establishment of the system is of important theoretical and practical significance for the research of KCCs and water disaster control.
Risk assessment of tunnel flood based on the weighting of index reliability measurement and set pair cloud
JIANG Yingli, CUI Jie, WANG Jingmei, ZHANG Yanlong
2022, 41(2): 276-286. doi: 10.11932/karst20220208
Abstract:
Tunnel flood assessment is a nonlinear and complex system with uncertainty. Scientific and reasonable flood grading evaluation and corresponding prevention and treatment measures have become the primary problem to be solved in the construction and operation of karst tunnels. Firstly, aimed at the uncertainty, fuzziness and randomness of evaluation index parameters, the correlation between indexes by the multi-evidence correlation coefficient of Jousselme distance are analyzed in this study. The dynamic weighting theory of index reliability measurement based on Jousselme distance is also put forward. With the weighting method, the accurate weight of each index can be obtained according to the reliability classification of indexes. Thus, the dynamic weighting of the whole system by the measured values of different cases and indexes is realized, and the risk of deviation of evaluation results caused by the error or error of measured values of indexes in practice is reduced; thus, the robustness of the evaluation model will be enhanced.Secondly, based on the theory and idea of uncertain artificial intelligence, and with the combination of the eigenvalues in cloud theory and the set pair theory, the cloud theory to optimize the set pair connection degree is used to obtain risk assessment indexes of karst tunnel flood. The index cloud connection degree and index reliability measurement are weighted by the dynamic weighting method to obtain the degree of system comprehensive cloud connection, and the risk value is obtained by the expected weighted average of the grading evaluation interval. At the same time, the corresponding grade cloud map is generated to determine the flood grade of karst tunnel. Accordingly, the flood state of karst tunnel is determined, and the visualization of flood grade determination is realized.Thirdly, the karst tunnel flood occurs in a complex system composed of both the tunnel and the external environment; therefore, its risk assessment is affected by many factors, and the evaluation index system should be established in a comprehensive and concise way. From the perspective that atmospheric precipitation is the main source of flood in karst tunnel, five correlation indexes-annual precipitation, infiltration coefficient, catchment area, permeability coefficient and unit water inflow-are constructed as evaluation indexes of set pair cloud of karst tunnel flood in this study. On this basis, the grading standard of each evaluation index and its corresponding cloud map of flood grade are established. In order to verify the accuracy and effectiveness of the set pair cloud model and its evaluation index system constructed based on the weighting of the index reliability measurement, the data of six typical karst tunnel samples for the model test is collected. Consistency between evaluation results and those of other methods proves the reliability and effectiveness of the model proposed in this study. Finally, this model has been applied to the flood accident of Pishuangao karst tunnel of Beijing-Zhuhai expressway, the results of which are also consistent with the actual situation of the project. At the same time, the treatment measures corresponding to Grade I have achieved effective flood control. The result shows that the set pair cloud model constructed in this study based on the weighting of the index reliability measurement takes into account the correlation between the evaluation indexes, the uncertainty of the evaluation system and the ambiguity of the evaluation index grade. The model can improve the accuracy of flood risk assessment in karst tunnels, and the evaluation process of this model is highly maneuverable. In fact, this model can provide a new method for rapid and effective analysis of flood in karst tunnels, and can also provide reference and guidance for the prediction and prevention of flood in karst tunnels in China.
Karst leakage and its sealing at Wanzi reservoir in Luoping county on the plateau of eastern Yunnan
DUAN Qiaowen, YU Fuyou, ZHANG Tianbai, HE Wei, DUAN Chunlin
2022, 41(2): 287-297. doi: 10.11932/karst20220209
Abstract:
Built on the Dagan river-a three-tier tributary of the Nanpan river from the Zhujiang river system, Wanzi reservoir is located at Wanzi village on the edge of Luoping basin 13 km from the southwest of Luoping county of Yunan Province. It is a typically karst-damaged reservoir in eastern Yunnan, with its lateral leakage toward lower adjacent valley. The landform of this reservoir belongs to erosion-corrosion alpin and medium-height mountains of eastern Yunan structure. The control drainage area of Wanzi reservoir is 104 km2 with the average runoff of 8,843×104 m3. This medium-sized reservoir covers a designed irrigation area of 2,133×104 m2 with its designed height of 23.5 m and normal capacity of 1,046×104 m3. It has experienced four times of leakage sealing since its construction in 1959, but the leakage reaches more than 0.5 m3·s−1 at medium water level with water supply of only 141.7×104 m3.In order to provide the theoretical foundation for further successful sealing and even the expanding of reservoir capacity, we investigate the main causes of the failure of leakage sealing. Meanwhile, we also profoundly analyze the factors of leakage according to the previous studies as well as the data about Wanzi reservoir such as stratum lithology, geological structure, geological conditions of hydrogeological engineering, characteristics of karst development, conditions of flow concentration, the failure of sealing, etc. Results show that because the study area is located on the strong karst development zone controlled by stratum lithology and structure, the vacuum suction erosion and the groundwater erosion caused by frequent and substantial water-level fluctuation result in corrosion fissures and the scouring of holes. The additional stress, composed of the water storage pressure and the negative pressure caused by the vacuum due to the coverage and enclosure of reservoir bottom, is larger than the shear strength, or larger than the strength bonding karst cave fillings and cave walls. Consequently, there exist failures of leakage sealing, sinkholes and the repeated karst collapse in the same leakage area, all of which are the main causes of leakage.In spite of several relatively separate leakage belts inside the reservoir area, all the leakage belts outside converge into the karst pipelines on both sides of the reservoir. These karst pipelines characterized by bedded development are narrow due to the holdup of clastic rocks. The leakage sealing spots scattering over a large area bringing huge engineering work and the failure of leakage protection. However, the water-tight curtain properly laid out outside the reservoir area can not only eradicate the leakage, but it can also effectively raise the groundwater level to improve storage conditions. In the way of laying out water-tight curtain, the reservoir capacity can be expanded more than 3 times as large as the previous one, equal to 3 other medium-sized reservoirs.
Study on physical simulation of spatial distribution identification of karst underground pipeline based on time-lapse 3D resistivity inversion
LIU Wei, ZHOU Qiyou, PAN Xiaodong, HE Changxiang
2022, 41(2): 298-307. doi: 10.11932/karst2022y04
Abstract:
The spatial distribution of karst underground river pipelines is complex and changeable; therefore, it is of great significance for us to carry out the earth science work in karst areas to develop an effective detection technology for the identification of pipelines. Based on the characteristics that karst underground river pipelines change greatly with seasons and time, the spatial distribution identification of underground river pipelines is studied by the inversion technology of time-lapse 3D resistivity. The results of pilot physical model experiments in limestone medium show that, the spatial distribution of simulated three-dimensional underground river pipeline is well highlighted according to the time-delay inversion imaging of resistivity data respectively collected from water-filled pipelines in rainy season and air-filled pipelines in dry season, The area with large resistivity change is in good agreement with the actual spatial position of pipeline in transverse and longitudinal directions. The closer the resistivity is to the center of the pipeline, the more significant the resistivity change is. In terms of the electrical sounding curve, the resistivity presents a low-high-low feature, and the resistivity contrast between pipeline space and other depth intervals is significant. In terms of the electrical sounding gradient curve, both the gradient value of pipeline space near the overlying clay layer and the value underlying limestone layer are large. The gradient value is positive near the upper contact surface and negative near the lower contact surface. The inversion effect with water-filled pipeline is inferior to that of time-lapse inversion. When the simulated pipeline is filled with water, the pipeline space presents relatively low resistivity. The resistivity value at the bottom of the pipeline is greater than that of the upper part, and the actual spatial position of the pipeline can be effectively reflected. The resistivity presents a high-low-high curve on the electrical sounding curve. However, on the electrical sounding gradient curve, the gradient value of the pipeline space near the overlying clay layer is negative, and it is positive near the underlying limestone layer. It is difficult to effectively identify the karst pipeline through the inversion images only generated by air-filled data. The resistivity rises monotonously and the resistivity value in the pipeline space changes more quickly than that in other depth intervals on the electrical sounding curve, while on the electrical sounding gradient curve, the resistivity gradient value in the pipeline space is slightly larger than that in other depth intervals, but the anomaly in the pipeline space is very weak compared with the time-lapse inversion results and water filling inversion results. It is suggested that when the spatial distribution of underground river pipelines in the field is studied, the 3D resistivity data should be collected in wet season and dry season, respectively, and the time-lapse 3D resistivity inversion image, electric sounding curve and electric sounding gradient curve can be comprehensively analyzed for the research. The subordinate choice is to carry out 3D detection of karst pipeline in wet season and analyze the data comprehensively, so as to obtain a practical detection technology that can effectively locate the underground river pipeline whether it is filled with water or free of water, which can be popularized and applied in practice.
Application of transient electromagnetism and cross-hole CT imaging to detect karst distribution and morphological characteristics﹕A case study of Jinan,Shandong Province
DAI Fangyuan, GAO Yang, SU Qingwei, HU Tao, GENG Fuqiang, DONG Yanan
2022, 41(2): 308-317, 328. doi: 10.11932/karst20220210
Abstract:
The study area is located in the southern mountainous area of Jinan City, which is an exposed-covered karst mountainous area with low hilly landform surrounded by mountains on the east, west and south sides. The overall terrain is high in the south and low in the north, and is generally a north inclined monoclinic structure with exposed strata of Cambrian,Ordovician and Quaternary. Precipitation in this area is mainly concentrated in July, August and September. The catchment area of nearby rivers is limited and cut off for a long time.The development of surface karst in the study area is ordinary, but the underground shallow karst is relatively developed. There are two types of groundwater,pore water and karst water. The karst water aquifer corresponds to the limestone of Cambrian-Ordovician Sanshanzi formation in which karst is relatively developed. The depth of study area is 50 m with dolomite as its the main lithology. There are two kinds of karst development in this area,small karst caves with the pore diameter greater than 20 cm, fast footage and low sampling rate and dissolved pores with the pore diameter less than 20 cm, normal footage, and high sampling rate. Caves and pores are basically filled with cohesive soil.This study is aimed to calculate the dispersion coefficient and eigenvalue of apparent resistivity and apparent absorption coefficient in karst development area, and to explore a method to analyze the development degree and morphological characteristics of karst cave. The research methods include Transient Electromagnetic Method (TEM), cross-hole electromagnetic tomography and drilling. In TEM, the central loop device, multi-turn small loop emission and point measurement are adopted. The distance between fixed sending points and measuring points of cross-hole electromagnetic tomography is 1.0 m, and data is collected in the four frequency bands of 6, 8, 10 and 12.In this study, the transient electromagnetic profile is compared with the karst development section exposed by drilling, and 120-ohm meter is determined as the boundary value for interpreting the karst development area, so as to delineate the karst distribution area. Then the drilling is carried out for verification, and its result is in conformity with that of TEM. Imaging detection of the cross-hole electromagnetic tomography in the verification hole and nearby boreholes is conducted to further verify and characterize the karst development between these two types of holes. By comparing TEM with drilling, the distribution range, the average value and the characteristic value of apparent resistivity in karst development area are statistically analyzed to delineate the distribution range of karst. Then the imaging of cross-hole electromagnetic tomography is used to precisely divide the distribution and development form of karst. Results show that apparent resistivity of small karst caves in TEM is relatively scattered, but it is larger and concentrated in dissolved pores. Under the same conditions of karst cave, both average resistivity and maximum apparent resistivity of non-filled karst caves are larger than those of the caves filled with cohesive soil.The results of these two geophysical prospecting methods show the same degree of dispersion. The degree of dispersion of small karst caves data is greater than that of dissolved pores, but the performance of data average value and characteristic value is different,with transient electromagnetic method, the average and characteristic values of apparent resistivity in small karst caves are less than those of dissolved pores; however, with cross-hole electromagnetic tomography, they are greater than those in dissolved pores.
Research progress and prospect of karst geomorphology in China based on digital elevation model
BI Benteng, YANG Chen, LI Jingwen, JIANG Jianwu, ZHOU Lixin
2022, 41(2): 318-328. doi: 10.11932/karst20220211
Abstract:
Digital elevation model (DEM) contains rich morphological information because the process of digital terrain modeling and extracting a series of terrain parameters based on DEM data with various spatial resolutions, known as Digital Terrain Analysis (DTA), has achieved a number of advances in topographic analysis, morphological modeling, distribution patterns and evolution patterns of karst landforms. This paper summarizes the new progress in spatial analysis of DEM based on the key technical methods of applying DEM to the previous Chinese studies on karst landform in recent years. The future research of DEM in karst geomorphology is also predicted in order to expand its application and provide technical ideas for the study of karst geomorphology from topology to the process and mechanism of topographic evolution.In this study, the applicability of the data model and the method of determining the optimal scale of analysis are discussed. Both DEM and DTA are obviously scale-dependent, and it is indisputable that the slope scales, watershed scales, and regional scales of karst topography require different resolution of DEM. Higher spatial resolution does not necessarily guarantee better automatic karst detection. Therefore, the main features of regional topography must be considered. When local terrain parameters are calculated and extracted, determining an optimal analysis window is necessary to reflect the completeness of the landform and guarantee the representativeness of calculation results.For the identification and classification of karst landforms, four effective methods for automatic extraction of karst landform units with DEM data and the advantages and disadvantages of these methods are analyzed. These methods are respectively based on cell statistics of composite factors, the extraction of terrain feature points, contour trees, and positive and negative topography. The first method is to divide the study area into a series of grids, calculate the topographic attributes of each cell grid, and construct an identification model including the establishment of peak forests, peak clusters and other karst geomorphic units. The second one is to construct a spatial relationship model of terrain feature points of different geomorphic units, and then automatically extract these units. The third one is to extract geomorphic contour lines by building a contour line recognition model for classification. The fourth one can reduce and simplify the complex and diverse terrains into positive and negative terrains to effectively highlight the morphological differences, and then divide them according to the characteristics of geomorphic units such as peak forests, peak clusters, etc.As for the morphology and pattern analysis of karst landform, the typical index parameters used to characterize the spatial pattern of karst landform are sorted out. In addition, the main applications of DEM in the ecological environment of karst area have been outlined. DEM can be used for topographic statistics, spatial pattern, spatial estimation and 3D visualization of meteorology, hydrology, soil and habitat, and then for the analysis of ecological environment changes in karst areas.However, the present classification system of karst landform is mainly focused on abstract expression and qualitative description, which can hardly quantify the detailed information of spatial structure of karst landform. This kind of system also cannot meet the demand of digitalization and intellectualization. In addition, with good recognition effect, the current methods for morphological and quantitative expression of karst landform have been widely applied to typical karst landforms such as peak forests, peaks and depressions. But for atypical karst landforms, DTA has rarely been tried.In conclusion, the prospective direction of DEM-based karst geomorphology research has been proposed in this study. First of all, a more effective quantitative classification system needs to be constructed for karst geomorphology and its topographic elements. A digital indicator system of karst geomorphology should be established to describe its spatial and geological properties in terms of spatial distribution relationships, change patterns, etc. Secondly, more comprehensive geomorphological information should be extracted from regional, watershed and slope scales to reveal the relationship between each indicator and karst development in the study on mechanism and evolution of karst development. Moreover, it is necessary to integrate remote sensing with other multi-source data to enrich the geological attributes of DEM and to delve into more information in soil, vegetation, geological structure and other fields.