Management and mitigation of sinkholes on karst lands: an overview of practical applications

When karst lands are developed for residential or commercial purposes, sinkholes are to be properly managed or mitigated. This paper discusses various types of sinkholes and their possible formation mechanisms, presents some guidelines for sinkhole management, and summarizes several engineering measures for sinkhole remediation. Because a sinkhole is often not an isolated feature but a component of an integrated groundwater and surface water system, knowledge of the sinkhole geology and hydrology is essential in selection of the proper management tools and remediation options. Adverse effects including flooding, additional sinkholes, and groundwater contamination can result from poorly managed or improperly repaired sinkholes.     

Formation mechanism of large sinkhole collapses in Laibin,Guangxi, China

On June 3, 2010, a series of karst sinkholes occurred at Jili village surrounded by Gui-Bei highway, Wu-Ping highway and Nan-Liu High-Speed Railway in Laibin, Guangxi, China. The straight-line distances from an large sinkhole pit, 85 m in diameter and 38 m in depth, to the above mainlines are 200, 600 and 500 m, respectively. Several investigation methods including geophysical technology, borehole and well drilling, groundwater elevation survey and hydrochemistry analysis of groundwater were used to study the formation mechanisms of these sinkholes. Based on the results, the spatial distribution of the Jili underground river was confirmed with a strike of SN along the middle Carboniferous limestone bedrock and the Quaternary deposits controlled the sinkhole formation. In addition, both historical sinkhole events and analysis of the groundwater–air pressure monitoring data installed in the underlying karst conduit system indicate that sinkholes in this area are more likely induced by extreme weather conditions within typical karst geological settings. Extreme weather conditions in the study area before the sinkhole collapses consisted of a year-long drought followed by continuous precipitation with a daily maximum precipitation of 442 mm between May 31 and June 1, 2010. Typical geological conditions include the Jili underground river overlain by the Quaternary overburden with thick clayey rubble. Especially in the recharge zone of the underground river, a stabilized shallow water table was formed in response to the extreme rainstorm because of the presence of the thick clayey rubble. When the underground conduit was flooded through the cave entrance on the surface, air blasting may have caused the cave roof collapse followed by formation of soil cavities and surface collapses. Borehole monitoring results of the groundwater–air pressure monitoring show that the potential karst sinkhole can pose threats to Shanbei village, the High-Speed Railway and the Wu-Ping highway. Local government needs to be aware of any early indicators of this geohazard so that devastating sinkholes can be prevented in the future. The results also suggest that groundwater–air pressure monitoring data collected both the Quaternary deposits and the bedrock karst system provide useful indicators for potential sinkhole collapses in similar karst areas where sinkholes usually occur during rainy season or karst groundwater level is always under the rockhead.     

Sinkhole distribution in a rapidly developing urban environment: Hillsborough County, Tampa Bay area, Florida

Sinkhole formation in Florida is a common event. The Florida karst plain is significantly altered by human development and sinkholes cause considerable property damage throughout much of the state. We present in this paper a morphometric analysis of karst depressions in the Tampa Bay area, and the relation with the known distribution of sinkholes. We selected the Tampa Bay area because it is particularly susceptible to the evolution of karst depressions in relation with development of the built-up environment. Karst depressions were mapped from the 1:24,000 USGS topographic maps and a morphometric analysis was performed by using parameters such as shape, circularity index, perimeter, area, length, width, and orientation. Maps showing the distribution of depression density, and the sectors with greatest areas of karst depression were produced using a GIS. These results were compared with data compiled from the database of sinkhole occurrences in Florida maintained by the Florida Geological Survey. Our analysis demonstrates that the distribution of new sinkhole occurrences differs from the distribution of existing sinkholes, indicating that there are processes acting today that are influencing karst landscape formation that are different from those acting in the past.     

Analysis of a database for anthropogenic sinkhole triggering and zonation in the Naples hinterland (Southern Italy)

Anthropogenic sinkholes are generally due to the collapse of man-made underground caves and represent a major threat, especially in urban contexts, where they could cause damage to people, buildings and lifelines. The hinterland of Naples (Campania, Southern Italy) is one of the most urbanized areas of Italy and is historically affected by frequent sinkhole phenomena. This study analyzes a database of both man-made caves and sinkholes collected by the authors over several years. The aim is to comprehend the predisposing and triggering factors of sinkholes in order to attempt a zonation of occurrence. The predisposing factor resulted to be the presence of a preexisting network of caves within the Campanian Ignimbrite tuff and their geometric features. Generally, the shallower the depth of the chamber roof and the lower the thickness of the vault, the higher is the frequency of sinkholes. Furthermore, an intrinsic fragility is represented by the access shafts of vertical wells, usually filled in and abandoned. Meanwhile, the main triggering mechanism identified consists in saturation of the subsoil, due to water leaks coming from buried sewage and water pipelines. The macrozonation of sinkhole occurrence shows that the highest class is achieved where the tuff is shallow and both man-made caves and historical sinkholes are present.

     

Roadway construction in karst areas: management of stormwater runoff and sinkhole risk assessment

In karst-rich regions, it is inevitable that roadways cross karst landscapes. Road building across such terranes faces environmental and engineering challenges because of impacts on water quality from stormwater runoff and concerns of sinkhole collapse under or near roadways. When highway runoff drains rapidly into subsurface conduit networks through open sinkholes and/or sinking streams, the impact of the runoff on the karst aquifer can be qualitatively evaluated by mixing cell models. Formulation of a comprehensive stormwater runoff management plan prior to roadway construction can minimize the associated adverse impacts. The commonly used best management practices help manage the stormwater runoff effectively in some sites. Site-specific management plans are preferable for other sites because of concerns of flooding and land stability. Proactive measures should be taken to identify areas of the greatest sinkhole collapse risk along the proposed route and the associated groundwater drainage patterns.     

Geophysical prediction and following development sinkholes in two Dead Sea areas, Israel and Jordan

Geophysical methods—seismic refraction (SRFR), electrical resistivity tomography (ERT), and microgravity—were applied to the Dead Sea (DS) sinkhole problem in the Ein Gedi area at the earlier stage of the sinkhole development (1998–2002). They allowed determining the sinkhole formation mechanism and localizing the sinkhole hazardous zones. The SRFR method permitted to delineate the underground edge of a salt layer at the depth of 50 m. The salt edge was shaped like the sinkhole line on the surface. It was concluded that the sinkhole development is linked to the salt edge. Geoelectrical quasi-3D mapping based on the ERT technique detected large resistivity anomalies with 250–300 m² diameter and 25–35 m deep. The Ein Gedi area has been also mapped by the use of Microgravity method. The residual Bouguer gravity anomaly map shows negative anomalies arranged along the edge of the salt layer. Those gravity anomalies overall are very similar in plan to the resistivity distribution in this area. The results of forward modeling indicate that both high resistivity and residual gravity anomalies are associated with a subsurface decompaction of the soil mass and deep cavity at the sinkhole site. Following monitoring of the sinkhole development carried out by the Geological Survey of Israel confirmed our suggestions. The drilling of numerous boreholes verified the location of the salt edge. Geographical Information System (GIS) database testifies that during 2003–2009 new sinkholes are continuing to develop along the salt edge within a narrow 50–100 m wide strip oriented approximately in north–south direction (slightly parallel to the shoreline). No promotion in west–east direction (perpendicularly to the DS shoreline) was observed in Israel. Collapse of sinkholes and their clustering have been occurred within the area of high resistivity anomaly and negative residual gravity anomaly. Similar studies carried out at the Ghor Al-Haditha area (Jordan) have shown that sinkholes there are also arranged along the winding line conforming to the salt edge. In this area sinkholes are slowly moved to the Dead Sea direction. Results of geophysical studies in numerous DS sites indicate similar sinkhole development. It allowed generating of the sinkhole formation model based on ancient (10,000–11,000-year old) salt belt girding the Dead Sea along its shores     

Relations between sinkhole density and anthropogenic contaminants in selected carbonate aquifers in the eastern United States

The relation between sinkhole density and water quality was investigated in seven selected carbonate aquifers in the eastern United States. Sinkhole density for these aquifers was grouped into high (>25 sinkholes/100 km²), medium (1–25 sinkholes/100 km²), or low (<1 sinkhole/100 km²) categories using a geographical information system that included four independent databases covering parts of Alabama, Florida, Missouri, Pennsylvania, and Tennessee. Field measurements and concentrations of major ions, nitrate, and selected pesticides in samples from 451 wells and 70 springs were included in the water-quality database. Data were collected as a part of the US Geological Survey (USGS) National Water-Quality Assessment (NAWQA) Program. Areas with high and medium sinkhole density had the greatest well depths and depths to water, the lowest concentrations of total dissolved solids and bicarbonate, the highest concentrations of dissolved oxygen, and the lowest partial pressure of CO₂ compared to areas with low sinkhole density. These chemical indicators are consistent conceptually with a conduit-flow-dominated system in areas with a high density of sinkholes and a diffuse-flow-dominated system in areas with a low density of sinkholes. Higher cave density and spring discharge in Pennsylvania also support the concept that the high sinkhole density areas are dominated by conduit-flow systems. Concentrations of nitrate-N were significantly higher (p < 0.05) in areas with high and medium sinkhole density than in low sinkhole-density areas; when accounting for the variations in land use near the sampling sites, the high sinkhole-density area still had higher concentrations of nitrate-N than the low sinkhole-density area. Detection frequencies of atrazine, simazine, metolachlor, prometon, and the atrazine degradate deethylatrazine indicated a pattern similar to nitrate; highest pesticide detections were associated with high sinkhole-density areas. These patterns generally persisted when analyzing the detection frequency by land-use groups, particularly for agricultural land-use areas where pesticide use would be expected to be higher and more uniform areally compared to urban and forested areas. Although areas with agricultural land use and a high sinkhole density were most vulnerable (median nitrate-N concentration was 3.7 mg/L, 11% of samples exceeded 10 mg/L, and had the highest frequencies of pesticide detection), areas with agricultural land use and low sinkhole density still were vulnerable to contamination (median nitrate-N concentration was 1.5 mg/L, 8% of samples exceeded 10 mg/L, and had some of the highest frequencies of detections of pesticides). This may be due in part to incomplete or missing data regarding karst features (such as buried sinkholes, low-permeability material in bottom of sinkholes) that do not show up at the scales used for regional mapping and to inconsistent methods among states in karst feature delineation.     

Logistic regression model for sinkhole susceptibility due to damaged sewer pipes

The occurrence of anthropogenic sinkholes in urban areas can lead to severe socioeconomic losses. A damaged underground sewer pipe is regarded as one of the primary causes of such a phenomenon. This study adopted the best subsets regression method to produce a logistic regression model that evaluates the susceptibility for sinkholes induced by damaged sewer pipes. The model was developed by analyzing the sewer pipe network as well as cases of sinkholes in Seoul, South Korea. Among numerous sewer pipe characteristics tested as explanatory variables, the length, age, elevation, burial depth, size, slope, and materials of the sewer pipe were found to influence the occurrence of sinkhole. The proposed model reasonably estimated the sinkhole susceptibility in the area studied, with an area value under the receiver-operating characteristics curve of 0.753. The proposed methodology will serve as a useful tool that can help local governments to choose a cavity inspection regime, and to prevent sinkholes induced by damaged sewer pipes.     

Drainage and flooding in karst terranes

Flooding in karst terranes is a commonly occurring geo-hazard. It causes damage to property, businesses, and roadways. It can lead to the formation of cover-collapse sinkholes and groundwater contamination. Generally, three types of flooding or their combinations are related to karst: recharge-related sinkhole flooding, flow-related flooding, and discharge-related flooding. Understanding of the type of flooding is essential for solving the flooding problem. Areas prone to karst flooding should be recognized, and restrictions and laws on land use should be implemented. Runoff and erosion control plans should address the unique characteristics of karst features. Digging out clogged sinkholes, creating retention basins, or installing Class V Injection Wells are possible solutions to improve drainage of storm water. Solutions to flooding problems in karst areas should also be coordinated with the water quality control to prevent groundwater contamination.     

Catastrophic subsidence: An environmental hazard,shelby county,Alabama

Induced sinkholes (catastrophic subsidence) are those caused or accelerated by human activities These sinkholes commonly result from a water level decline due to pumpage Construction activities in a cone of depression greatly increases the likelihood of sinkhole occurrence Almost all occur where cavities develop in unconsolidated deposits overlying solution openings in carbonate rocks. Triggering mechanisms resulting from water level declines are (1) loss of buoyant support of the water, (2) increased gradient and water velocity, (3) water-level fluctuations, and (4) induced recharge Construction activities triggering sinkhole development include ditching, removing overburden, drilling, movement of heavy equipment, blasting and the diversion and impoundment of drainage Triggering mechanisms include piping, saturation, and loading Induced sinkholes resulting from human water development/management activities are most predictable in a youthful karst area impacted by groundwater withdrawals Shape, depth, and timing of catastrophic subsidence can be predicted in general terms Remote sensing techniques are used in prediction of locations of catastrophic subsidence. This provides a basis for design and relocation of structures such as a gas pipeline, dam, or building Utilization of techniques and a case history of the relocation of a pipeline are described     

Karst database implementation in Minnesota: analysis of sinkhole distribution

This paper presents the overall sinkhole distributions and conducts hypothesis tests of sinkhole distributions and sinkhole formation using data stored in the Karst Feature Database (KFD) of Minnesota. Nearest neighbor analysis (NNA) was extended to include different orders of NNA, different scales of concentrated zones of sinkholes, and directions to the nearest sinkholes. The statistical results, along with the sinkhole density distribution, indicate that sinkholes tend to form in highly concentrated zones instead of scattered individuals. The pattern changes from clustered to random to regular as the scale of the analysis decreases from 10–100 km² to 5–30 km² to 2–10 km². Hypotheses that may explain this phenomenon are: (1) areas in the highly concentrated zones of sinkholes have similar geologic and topographical settings that favor sinkhole formation; (2) existing sinkholes change the hydraulic gradient in the surrounding area and increase the solution and erosional processes that eventually form more new sinkholes.     

The hydrologic effects of intense groundwater pumpage in east-central Hillsborough County,Florida, U.S.A.

Well problems, water shortages, local flooding, and induced sinkholes have been periodic problems for residents in east-central Hillsborough County, Florida. This agricultural area has experienced dramatic short-term water-level declines in the Floridan aquifer from seasonal groundwater withdrawals. The sudden declines in the potentiometric surface have been caused from intense irrigation pumpage, primarily for frost and freeze protection and fruit setting. Citrus and strawberry crops are protected from occasional freezes by the application of warm groundwater to maintain minimum soil temperatures of 32°F(0°C). Local residents with inadequately constructed wells lose their source of water when the potentiometric surface is lowered to depths where their wells do not function. Some residents have lost their water supply for a week or more, and many have incurred damage to their pumps. The drawdown of the potentiometric surface in some areas has induced sinkholes causing property damage for some residents and concern for others. In addition, the high application rates for frost and freeze irrigation have created run-off problems resulting in local flooding to some residents. Fortunately, there has been no damage to resident homes from the flooding or sinkholes. This report summarizes the area's hydrogeology, and the consequences of heavy freeze irrigation. A finite-difference, numerical model is used to quantify the regional impacts to the potentiometric surface of the Floridan aquifer from estimated irrigation pumpage. Management recommendations for alleviating the problems are also discussed in the report.     

Karst unit mapping using geographic information system technology,Mower County,Minnesota, USA

Mower County is in southeastern Minnesota in an area underlain by sedimentary bedrock. These rocks are karsted Middle Devonian and Middle Ordovician limestone and dolomite. The karst features of the county were inventoried. These features include sinkholes, disappearing streams, caves, dry valleys, and springs. Previous karst mapping efforts for other counties in Minnesota have produced sinkhole probability maps and a springshed map. In Mower County, we have developed a new type of karst map using a geographic information system (GIS) to produce a karst unit map. Karst units are discrete three-dimensional bodies in which solution of the bedrock has resulted in the integration of surface water and groundwater. The field mapping and hydrologic investigations were done with conventional methods. The karst unit delineation was done using GIS technology, which allowed us to examine the county's karst using two-dimensional and three-dimensional views. Many different overlays of the karst elements were combined to better understand the landscape dynamics. Ultimately, it was the overlay of the karst features, hydrologic information, and depth-to-bedrock mapping on a shaded relief landscape morphology map which allowed us to best delineate the individual karst units.     

地下空洞引起土质地层地陷的形成机制与预测研究

地下空洞引起的土质地层地陷是地质灾害的一种常见形式。由于地陷产生的原因非常复杂,目前相关研究还处于起步阶段。通过地陷成因和实例分析,首先确认了软土和粗粒土地层的地陷破坏模式。然后对软土地层圆柱形地陷,采用土体塑性极限平衡理论导出了计算地陷安全系数的公式,并据此分析了影响地陷的相关因素。对粗粒土地层地陷的圆形塌陷漏斗,基于散体流动理论提出了粗粒土流失判据,推导了地陷安全系数计算公式;通过分析粗粒土地层地陷形成规律,提出了基于土层厚度和地陷漏斗最大沉降的地陷预测与评估方法。最后,采用数值模拟方法分析了软土和粗粒土地层地陷时的地层变形规律,对所提地陷模型的合理性进行了说明。本文研究成果可为认识土质地层地陷的形成机制和地陷防治提供参考。     

Sinkhole formation and subsidence along the Dead Sea coast,Israel

More than 4,000 sinkholes have formed since the 1980s within a 60-km-long and 1-km-wide strip along the western coast of the Dead Sea (DS) in Israel. Their formation rate accelerated in recent years to >400 sinkholes per year. They cluster mostly in specific sites up to 1,000 m long and 200 m wide, which align parallel to the general direction of the fault systems associated with the DS Rift. The abrupt appearance of the sinkholes reflects changes to the groundwater regime around the shrinking DS. The eastward retreat of the shoreline and the lake-level drop (1 m/year in recent years) cause an eastward and downward migration of the fresh/saline groundwater interface. Consequently, a subsurface salt layer, which was previously enveloped by saline groundwater, is gradually being invaded and submerged by relatively fresh groundwater, and cavities form due to the rapid dissolution of the salt. Collapse of the overlying sediments into these cavities results in sinkholes at the surface. An association between sinkhole sites and land subsidence is revealed by interferometric synthetic aperture radar (InSAR) measurements. On a broad scale (hundreds of meters), subsidence occurs due to compaction of fine-grained sediments as groundwater levels decline along the retreating DS shoreline. At smaller scales (tens of meters), subsidence appears above subsurface cavities in association with the sinkholes, serving in many cases as sinkhole precursors, a few weeks to more than a year before their actual appearance at the surface. This paper overviews the processes of sinkhole formation and their relation to land subsidence.     

The influence of urbanization of sinkhole development in central Pennsylvania

The karsted limestone valleys of central Pennsylvania contain two populations of sinkholes. Solution sinkholes occur in the Champlainian limestone units along the margins of the valleys. Solution sinkholes are permanent parts of the landscape and, although a nuisance to construction, do not present other problems. The second population is the suffosional or soil-piping sinkholes These occur on all carbonate rock units including the Beekmantown and Gatesburg dolomites that comprise the two principal carbonate aquifers in the valley. Suffosional sinkholes are the principal land-use hazard. Suffosional sinkholes are transient phonomena. They occur naturally but are exacerbated by runoff modifications that accompany urbanization Suffosional sinkholes are typically 1.5–2.5 m in diameter depending on soil thickness and soil type. The vertical transport of soil to form the void space and soil arch that are the precursors to sinkhole collapse is through solutionally widened fractures and cross-joints and less often through large vertical openings in the bedrock. The limited solution development on the dolomite bedrock combined with soil thickness, seldom greater than 2 m, limits the size of the sinkholes. All aspects of suffosional sinkhole development are shallow processes: transport, piping, void and arch formation, and subsequent collapse take place usually less than 10 m below the land surface Factors exacerbating sinkhole development include pavement, street, and roof runoff which accelerates soil transport Such seemingly minor activities as replacing high grass and brush with mowed grass is observed to accelerate sinkhole development. Dewatering of the aquifer is not a major factor in this region     

广西来宾市良江镇吉利村岩溶塌陷成因机制分析

2010年6月,一场极端暴雨事件过后,广西来宾市良江镇吉利村发生大规模岩溶塌陷,一个月内陆续产生了17处塌陷坑,其中3个塌陷坑逐渐扩大发展,最终形成一长约170 m,深约38 m的大塌陷坑,被称为“广西第一塌”。此次塌陷事件影响面积约1 500 m2,导致70户298人搬离家园。为了弄清吉利村岩溶塌陷的成因机制,在塌陷发生地区开展了地面调查、三维激光扫描、地球物理勘探、水文地质工程地质钻探、岩溶塌陷动力条件监测、降雨量监测等多项工作。分析结果表明:吉利岩溶塌陷为自然塌陷,是多种不良因素共同作用的结果,这些因素包括特殊的“水-土-岩”组合和极端暴雨,其中极端暴雨引起的地下水气压力的急剧变化是其主要动力因素;特殊的地质条件,决定了该区岩溶塌陷发生的特殊性,分析其形成过程,认为主要受气爆效应、水击效应和渗流潜蚀效应的综合作用。      

不同补给条件下裂隙-管道介质间水流交换的示踪试验研究

为探究岩溶含水系统中裂隙-管道介质间不同补给方式下水流交换情况,在实验室内借助高精度摄像仪和有色示踪剂对其进行可视化观测。设计多组试验对单独补给裂隙、共同补给、裂隙和落水洞分开同时补给三种补给方式下水流交换情况进行探讨,从而定性地分析水流交换情况。试验结果表明:落水洞对裂隙水扩散的影响随裂隙距离落水洞的距离增大而减小;相较于共同补给情况,单独补给裂隙情况下落水洞对裂隙水扩散影响较大;在分开同时补给情况下裂隙与落水洞补给水量相同时,仅呈现裂隙对落水洞的单向补给。     

Subsidence hazard avoidance based on geomorphological mapping in the Ebro River valley mantled evaporite karst terrain (NE Spain)

In the valley of the Ebro River to the southeast of the city of Zaragoza (NE Spain), the dissolution of evaporite sediments (gypsum, halite and Na-sulphates) which underlie alluvial deposits gives rise to numerous sinkholes. These sinkholes are a potential hazard to human safety, particularly where they develop in a catastrophic way. Even slow-developing sinkholes are problematic, as they damage urban and agricultural infrastructure, necessitating costly repairs and vigilant maintenance. To assist in developing avoidance strategies for these hazards, the factors controlling sinkhole occurrence have been assessed using geomorphological maps produced from aerial photographs for 1956 and 1981. Important controls on sinkhole development are found to include underlying geological structure (manifest in preferred orientations of sinkholes on the azimuths N130-150E and N30-40E), and the presence of glauberite in the groundwater flow path, which apparently promotes accelerated gypsum dissolution. Perhaps surprisingly, alluvium thickness does not appear to significantly correlate with the density of sinkholes on the floodplain in this area. The maps for 1956 and 1981 reveal that both human activity and natural processes can serve to obscure the true density of sinkhole development. For instance, a large number of sinkholes which were conspicuous in 1956 have since been back-filled by farmers. In the most fluvially active zone of the Ebro valley (the meander belt), the relatively low density of sinkholes compared with adjoining zones suggests that subsidence is being masked by morpho-sedimentary dynamic processes (aggradation and erosion). Careful geomorphological mapping for different time periods yields a much more accurate impression of the frequency of sinkhole development than would be gained from surveying currently visible sinkholes in the area of interest.     

湖南宁乡大成桥充水矿山疏干区岩溶系统水气压力监测及突变特征

湖南煤炭坝煤矿是我国著名的大水矿山,开采标高已降至-400 m,每小时排水量11 000 m3。长期的疏干排水,已形成面积达219.19 km2的地下水降落漏斗。以往的调查资料表明,矿山疏干诱发的岩溶塌陷主要发生在疏干开始、地下水位下降到基岩面附近、矿坑发生突水突泥灾害等3个阶段,在疏干区应该无诱发因素导致岩溶塌陷的发生。但近年来煤炭坝煤矿疏干区大成桥的岩溶塌陷特征表明,尽管相关部门采取了一系列处置措施,但岩溶塌陷仍然频繁发生。为了进一步揭示该地区岩溶塌陷发育规律,使处置决策更加科学化,作者所在团队选择大成桥为典型区,通过岩溶系统地下水气压力监测,分析诱发岩溶塌陷的动力来源及其变化规律。结果表明,疏干区岩溶管道裂隙系统中水气压力突变现象极为明显,年最大变幅达54.72 m,瞬时水气压力变化速率最大为70.6 cm/min,造成这一现象的主要原因是气爆、地下洞穴（溶洞、土洞）垮塌。