Mapping the susceptibility to sinkholes in coastal areas, based on stratigraphy, geomorphology and geophysics

Sinkholes and land subsidence are among the main coastal geologic hazards. Their occurrence poses a serious threat to the man-made environment, due to the increasing density of population, pipelines and other infrastructures along the coasts, and to the catastrophic nature of the phenomena, which generally occur without any premonitory signs. To assess the potential danger from sinkholes along the coast, it is important to identify and monitor the main factors contributing to the process. This article reports a methodology based on sequential stratigraphic, hydrogeological and geophysical investigations to draw up a susceptibility map of sinkholes in coastal areas. The town of Casalabate situated in the Apulia region (southern Italy), affected by a long history of sinkhole phenomena, is here presented as an example. The approach proposed is based on sequential stratigraphical, geomorphological and geophysical surveys to identify the mechanisms of sinkhole formation and to provide a zonation of the areas in which further sinkhole phenomena may likely occur. Interpretation of the ground penetration radar and electrical tomography profiles has enabled us to identify the potentially most unstable sectors, significantly improving the assessment of the sinkhole susceptibility in the area. The proposed methodology is suitable to be exported in other coastal areas where limestone bedrock is not directly exposed at the surface, but covered by a variable thickness of recent deposits.     

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     

Development and morphometry of sinkholes in coastal plains of Apulia, southern Italy. Preliminary sinkhole susceptibility assessment

Evolution of coastlines in karst areas may be strongly controlled by dissolution processes which favour the development of surface and subsurface landforms. The generation of caves in these environments is commonly favoured by the mixing between fresh and brackish waters. The sinkholes resulting from the upward propagation of the caves may interfere with the anthropogenic environment and cause damage to human elements (property and activities). To highlight the often underestimated importance of karst phenomena in coastal areas, we have analyzed a coastal stretch of Apulia, in southern Italy. The study area, covering an extension of about 6 km², is situated in the Ionian coast, and presents several interesting karst landforms that are generally connected to caves. Tens of sinkholes were mapped through field surveys, multi-year aerial-photographs (dating back to the 1940s) and archival research. We have performed a morphometric analysis of the sinkholes. The analysis describes the main parameters of the sinkholes (area, length, width, and depth), and the control exerted by the main discontinuity systems in the area. The detrimental effects derived from interaction between human environment and these karst landforms is also under consideration. A sinkhole susceptibility map, which may provide useful information for planners, developers and the insurance industry has eventually been produced through application of a decision tree model.     

Origin and distribution of different types of sinkholes in the plain areas of Southern Italy

Sinkholes constitute a significant risk in many karst areas, and may even threat human safety. Collapse sinkholes that occur catastrophically without showing premonitory signs may result in severe economic losses and casualties. In the last years, research on sinkholes and related detrimental effects has significantly increased in Italy, in the aftermath of remarkable events.

Aimed at analysing the peculiar conditions which lead to sinkhole occurrence in Southern Italy, a set of cases in the plain areas of Campania, Apulia and Calabria is discussed. The considered regions show a wide variety of environmental conditions. In Campania, the plains are underlain by alluvial deposits with intercalations of volcaniclastic sediments. Sinkholes are generally located along the Tyrrhenian margin of the carbonate massifs or within intramontane Apennine basins. In Apulia, a flat and elongated peninsula, most of the cases occur on calcarenites overlying limestone bedrock along the coast. In Calabria, one of the most seismogenic Italian regions, the surveyed cases seem to be attributable mainly to earthquake-induced liquefaction.

The article provides a first glance on the variety of sinkholes in the plain areas of Southern Italy, to highlight the possibility of further subsidence events in the considered regions, as well as in other comparable areas of the country.     

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

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

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

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,造成这一现象的主要原因是气爆、地下洞穴（溶洞、土洞）垮塌。      

Spatial distribution, morphometry and activity of La Puebla de Alfindén sinkhole field in the Ebro river valley (NE Spain): applied aspects for hazard zonation

A highly active collapse sinkhole field in the evaporitic mantled karst of the Ebro river valley is studied (NE Spain). The subsidence is controlled by a NW-SE trending joint system and accelerated by the discharge of waste water from a nearby industrial state. The morphometry, spatial distribution and temporal evolution of the sinkholes have been analysed. The volume of the sinkholes yields a minimum estimate of average lowering of the surface by collapse subsidence of 46 cm. The clustering of the sinkholes and the tendency to form elongated uvalas and linear belts, in a NW–SE direction have a predictive utility and allow the establishment of criteria for a hazard zonation. With the precipitation record supplied by a pluviograph and periodic cartographic and photographic surveys the influence of heavy rainfall events on the triggering of collapses has been studied.     

Identification, prediction, and mitigation of sinkhole hazards in evaporite karst areas

Sinkholes usually have a higher probability of occurrence and a greater genetic diversity in evaporite terrains than in carbonate karst areas. This is because evaporites have a higher solubility and, commonly, a lower mechanical strength. Subsidence damage resulting from evaporite dissolution generates substantial losses throughout the world, but the causes are only well understood in a few areas. To deal with these hazards, a phased approach is needed for sinkhole identification, investigation, prediction, and mitigation. Identification techniques include field surveys and geomorphological mapping combined with accounts from local people and historical sources. Detailed sinkhole maps can be constructed from sequential historical maps, recent topographical maps, and digital elevation models (DEMs) complemented with building-damage surveying, remote sensing, and high-resolution geodetic surveys. On a more detailed level, information from exposed paleosubsidence features (paleokarst), speleological explorations, geophysical investigations, trenching, dating techniques, and boreholes may help in investigating dissolution and subsidence features. Information on the hydrogeological pathways including caves, springs, and swallow holes are particularly important especially when corroborated by tracer tests. These diverse data sources make a valuable database—the karst inventory. From this dataset, sinkhole susceptibility zonations (relative probability) may be produced based on the spatial distribution of the features and good knowledge of the local geology. Sinkhole distribution can be investigated by spatial distribution analysis techniques including studies of preferential elongation, alignment, and nearest neighbor analysis. More objective susceptibility models may be obtained by analyzing the statistical relationships between the known sinkholes and the conditioning factors. Chronological information on sinkhole formation is required to estimate the probability of occurrence of sinkholes (number of sinkholes/km² year). Such spatial and temporal predictions, frequently derived from limited records and based on the assumption that past sinkhole activity may be extrapolated to the future, are non-corroborated hypotheses. Validation methods allow us to assess the predictive capability of the susceptibility maps and to transform them into probability maps. Avoiding the most hazardous areas by preventive planning is the safest strategy for development in sinkhole-prone areas. Corrective measures could be applied to reduce the dissolution activity and subsidence processes. A more practical solution for safe development is to reduce the vulnerability of the structures by using subsidence-proof designs.     

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.     

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.     

Detection of sinkholes using 2D electrical resistivity imaging in the Cheria Basin (north–east of Algeria)

Sinkhole collapse is one of the main limitations on the development of karst areas, especially where bedrock is covered by unconsolidated material. Studies of sinkhole formation have shown that sinkholes are likely to develop in cutter (enlarged joint) zones as a result of subterranean erosion by flowing groundwater. Electrical resistivity imaging or tomography (RESTOM) is well suited to mapping sinkholes because of the ability of the technique for detecting resistive features and discriminating subtle resistivity variations. Two-dimensional electrical resistivity tomography surveys were conducted at two sinkhole sites near Cheria city where limestone is covered by about 10 m of clayey soils. A Wenner transect was conducted between the two sinkholes. The electrode spacing was 2 m. The length of transect is about 80 m. The survey results suggest that RESTOM is an ideal geophysical tool to aid in the detection and monitoring of sinkholes and other subsurface cavities.     

Models and mechanisms of drilling-induced sinkhole in China

The present paper distinguishes the four general types of collapsed sinkholes induced by drilling in regions where non-indurated sediments (clay, silt, and sand) overlie unknown caves or open voids in underlying lithified karst rocks. These sinkhole models are classified into balance arch, hourglass, vibration, and subjacent drilling-induced erosion forms. A mechanical model was built for each type to assess the collapse probability. Drilling operations using boreholes facilitate rapid, turbulent, and erosive flow. These activities resulted in the 26 sinkholes that have been attributed to constructions in China in the last several years. Awareness of potential collapse conditions should allow construction projects to avoid high-risk settings that result in economic losses, environmental concerns, and life-threatening accidents caused by such rapid collapse of the land surface. Awareness and close monitoring of high-risk conditions during borehole monitoring should minimize the hazard.     

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.     

Groundwater fluxes into a submerged sinkhole area,Central Italy,using radon and water chemistry

The groundwater contribution into Green Lake and Black Lake (Vescovo Lakes Group), two cover collapse sinkholes in Pontina Plain (Central Italy), was estimated using water chemistry and a ²²²Rn budget. These data can constrain the interactions between sinkholes and deep seated fluid circulation, with a special focus on the possibility of the bedrock karst aquifer feeding the lake. The Rn budget accounted for all quantifiable surface and subsurface input and output fluxes including the flux across the sediment–water interface. The total value of groundwater discharge into Green Lake and Black Lake (∼540 ± 160 L s⁻¹) obtained from the Rn budget is lower than, but comparable with historical data on the springs group discharge estimated in the same period of the year (800 ± 90 L s⁻¹). Besides being an indirect test for the reliability of the Rn-budget “tool”, it confirms that both Green and Black Lake are effectively springs and not simply “water filled” sinkholes. New data on the water chemistry and the groundwater fluxes into the sinkhole area of Vescovo Lakes allows the assessment of the mechanism responsible for sinkhole formation in Pontina Plain and suggests the necessity of monitoring the changes of physical and chemical parameters of groundwater below the plain in order to mitigate the associated risk.     

The application of GPR and electrical resistivity tomography as useful tools in detection of sinkholes in the Cheria Basin (northeast of Algeria)

Sinkhole collapse is one of the main limitations in the development of karst areas, especially where bedrock is covered by unconsolidated material. Studies of sinkhole formation have shown that sinkholes are likely to develop in cutter (enlarged joint) zones as a result of subterranean erosion by flowing groundwater. Ground-penetrating radar (GPR) and electrical resistivity imaging or tomography (RESTOM) are well suited to mapping sinkholes because of the ability of these two techniques for detecting voids and discriminating subtle resistivity variations. Nine GPR profiles and two-dimensional electrical resistivity tomography have been applied, with relative success, to locate paleo-collapses and cavities, and to detect and characterize karst at two sinkhole sites near Cheria City where limestone is covered by about 10 m of clayey soils. The survey results suggest that GPR and RESTOM are ideal geophysical tools to aid in the detection and monitoring of sinkholes and other subsurface cavities.     

Conceptual site models for sinkhole formation and remediation

Extensive research has been done on investigation, monitoring, risk evaluation and land management to prevent formation of sinkholes and subsidences in karst terranes. Little emphasis is, however, given to the various processes and methodologies with respect to their remediation. Sinkholes are a surface symptom of complicated erosion and deformation processes that occur on the surface and in the subsurface. The specific method or combination of methods used to mitigate a sinkhole or subsidence depends on the complexity of the conceptual site model composed of influencing factors and triggering mechanisms, depth and lateral extent of instability, impact on existing infrastructure, and existing and reasonably anticipated land uses. Many sinkhole and subsidence rehabilitation methods including dynamic compaction, construction of inverted filters, compaction grouting and construction of water plugs are presented; however, each sinkhole must be treated uniquely when a decision is made on the most appropriate rehabilitation method. A case study was presented to illustrate how the conceptual site model and remedial alternative analysis approach were used to determine the preferred method for sinkhole remediation.     

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.     

覆盖型岩溶平原区岩溶塌陷脆弱性和开发岩溶地下水安全性评价——以广西黎塘镇为例

黎塘研究区塌陷脆弱性较高,为科学地开发该区的岩溶地下水,以满足工农业生产、生活需要,本文针对黎塘镇塌陷的主控触发因素为岩溶地下水动态变化这一特点,采用现场实时监测岩溶地下水动态变化特征和室内塌陷模拟试验相结合的方法,将监测结果与当地塌陷的临界条件进行对比来确定开发岩溶地下水的安全性.结果表明,在当前开采条件下,塌陷脆弱性分区与稳定分区并非完全重合:部分高脆弱和较高脆弱区,在目前开采条件下仍处于稳定状态,部分低脆弱区则处于危险状态,并为研究区十多年来塌陷的发育、分布特征所支持印证.据此认为,在塌陷脆弱性较高的地区未必就需要绝对禁止对岩溶地下水的开发,地下水的开发应以不增加当地塌陷的频率、加剧当地岩溶塌陷的危害为前提.