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.     

A genetic classification of sinkholes illustrated from evaporite paleokarst exposures in Spain

This contribution analyses the processes involved in the generation of sinkholes from the study of paleokarst features exposed in four Spanish Tertiary basins. Bedrock strata are subhorizontal evaporites, and in three of the basins they include halite and glauberite in the subsurface. Our studies suggest that formation of dolines in these areas results from a wider range of subsidence processes than those included in the most recently published sinkhole classifications; a new genetic classification of sinkholes applicable to both carbonate and evaporite karst areas is thus proposed. With the exception of solution dolines, it defines the main sinkhole types by use of two terms that refer to the material affected by downward gravitational movements (cover, bedrock or caprock) and the main type of process involved (collapse, suffosion or sagging). Sinkholes that result from the combination of several subsidence processes and affect more than one type of material are described by combinations of the different terms with the dominant material or process followed by the secondary one (e.g. bedrock sagging and collapse sinkhole). The mechanism of collapse includes any brittle gravitational deformation of cover and bedrock material, such as upward stoping of cavities by roof failure, development of well-defined failure planes and rock brecciation. Suffosion is the downward migration of cover deposits through dissolutional conduits accompanied with ductile settling. Sagging is the ductile flexure of sediments caused by differential corrosional lowering of the rockhead or interstratal karstification of the soluble bedrock. The paleokarsts we analysed suggest that the sagging mechanism (not included in previous genetic classifications) plays an important role in the generation of sinkholes in evaporites. Moreover, collapse processes are more significant in extent and rate in areas underlain by evaporites than in carbonate karst, primarily due to the greater solubility of the evaporites and the lower mechanical strength and ductile rheology of gypsum and salt rocks.     

Reliability of dipole-dipole electrical resistivity tomography for defining depth to bedrock in covered karst terranes

 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. Because of the irregular distribution of pinnacles and cutters on the bedrock surface, uncertainties arise when "hit-or-miss" borehole drilling is used to locate potential collapse sites. A high-resolution geophysical technique capable of depicting the details of the bedrock surface is essential for guiding the drilling program. Dipole-dipole electrical resistivity tomography (ERT) was used to map the bedrock surface at a site in southern Indiana where limestone is covered by about 9 m of clayey soils. Forty-nine transects were conducted over an area of approximately 42,037 m². The electrode spacing was 3 m. The length of the transects varied from 81 to 249 m. The tomographs were interpreted with the aid of soil borings. The repeatability of ERT was evaluated by comparing the rock surface elevations interpreted from pairs of transects where they crossed each other. The average difference was 2.4 m, with a maximum of 10 m. The discrepancy between interpreted bedrock-surface elevations for a transect intersection may be caused by variations in the subsurface geology normal to the transect. Averaging the elevation data interpreted from different transects improved the ERT results. A bedrock surface map was generated using only the averaged elevation data at the transect junctions. The accuracy of the map was further evaluated using data from four exploratory boreholes. The average difference between interpreted and actual bedrock surface-elevations was less than 0.4 m. The map shows two large troughs in the limestone surface: one coinciding with an existing sinkhole basin, while the other is in alignment with a small topographic valley. Because sinkholes were observed at the same elevation interval in similar valleys in the vicinity, the delineated trough may have implications for future land use at the site. Received: 4 January 1999 · Accepted: 8 March 1999     

Incipient vertical traction carpets within collapsed sinkhole fills

Small vertically oriented traction carpets are reported from the collapsed sandy fills of 100 m deep Devonian limestone sinkholes underlying the Lower Cretaceous Athabasca oil sands deposit in north‐eastern Alberta, Western Canada. Dissolution of 100 m of underlying halite salt beds caused cataclysmic collapse of the sinkhole floors and water saturated sinkhole sand fills to descend very rapidly. Turbulent currents flushed upper sinkhole fills of friable sandstone blocks and disaggregated sand and quartz pebble for tens of metres. Laminar deposits with inverse grading accumulated as many as six to eight curvilinear entrained pebble streaks, 10 to 30 cm long, vertically impinged against the sides of descending collapse blocks. These deposits were initiated as vertically oriented early stage traction carpets that interlocked fine sand grains and inversely graded overlying pebbles entrained below the dilute overlying turbulent flows. Vortexes that flushed these sinkhole fills and induced these depositional processes may have lasted only seconds before the very rapid descents abruptly halted. Some of the fabrics were suspended vertically in‐place and preserved from unlocking and obliteration. These small fabrics provide insight into the instability and ephemeral character of the transition from strong gravity‐driven grain falls to very early stages of traction carpet formation. These short‐lived deposits of very thin sand layers resulted from sufficient incipient frictional freezing that grain interlocking overcame, however briefly, the strong gravity drives of the vertical falls that would have otherwise dispersed grains and obliterated any organized fabric patterns. Tenuous frictionally locked grains were also suspended at the centres of hyperbolic grain fall flows that briefly developed between turbulent flow eddies, some of which were fortuitously preserved. Some of these suspended grain locking zones passed downward onto the relatively more stable surfaces of the rapidly descending block surfaces. The morphogenesis of these early stage traction carpets differ from more fully developed deposits elsewhere because of their short‐lived transport, dynamic instability and vertical orientation.     

Subsidence and uplift of Sidoarjo (East Java) due to the eruption of the Lusi mud volcano (2006–present)

Global positioning system (GPS) and satellite-based InSAR (Interferometric Synthetic Aperture Radar) measurements of the subsidence and uplift of a populated area of Sidoarjo, East Java are due to the eruption of the Lusi mud volcano (2006–present). These data are the first direct quantitative measurements of deformation due to the growth of a mud volcano edifice. The GPS data were recorded over periods of a few hours to several months and show that between June 2006 and September 2007, the earth’s surface has been subsiding at rates of 0.1–4 cm/day. Maximum rates of subsidence occurred in an area 300–400 m to the northwest of the main mud volcano vent. Horizontal displacements were 0.03–0.9 cm/day and were also towards this area. In general uplifts of up to 0.09 cm/day were recorded in areas outside of the edifice. Changes in elevation measured using satellite imagery (InSAR technique) provide regional datasets of subsidence and uplift. They confirm that during the first year a roughly circular area was undergoing sag-like subsidence centered to the northwest of the main vent and that uplift was occurring 3–4 months after the initiation of the eruption due to the movement Watukosek fault system. Subsidence occurred due to the weight of mud and man-made dams and the collapse of the overburden due to removal of mud from the subsurface. Assuming constant rates of subsidence of 4 cm/day, then in the centre of the edifice there would be up to 44 m of subsidence in 3 years, and up to 16 m in 10 years. The mud volcano is now in a self-organizing state with new fluid conduits forming as a result of the collapse. An erratum to this article can be found at     

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     

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.     

Evaporite karst and sinkholes: a synthesis on the case of Camaiore (Italy)

In 1995 a sinkhole suddenly formed at Camaiore (Tuscany, Italy), causing destruction or heavy damages to several houses and resulting in the evacuation of many people. To understand the causes, for the formation and evolution of the collapse, surface and underground geologic features were investigated and reconstructed on the basis of geologic and geognostic surveys. The sinkhole area is underlain by thick alluvial deposits that cover a bedrock consisting of the Calcare cavernoso formation. This formation results from hydration and dissolution of Triassic evaporites and has a characteristic spongy and vacuolate texture. The bedrock contains karst cavities, generally filled by breccia and/or alluvial materials. Thus, the sinkhole disaster could be ascribed to deep collapse of a cave in the bedrock, and might be considered a distant effect of ancient karst phenomena in evaporites.     

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.     

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.     

佛山地铁塌陷InSAR时序监测及机理分析

2018年2月7日,位于广东省佛山市禅城区的地铁2号线在盾构施工中发生塌陷事故,造成11人死亡、1人失踪、8人受伤,直接经济损失超过5000万元。为深入分析此次事故成因,本文基于自2017年3月~2019年1月期间的56景Sentinel-1A数据,利用SBAS-InSAR技术获取了研究区的时空形变信息。结果发现塌陷区及其附近区域在监测期间存在持续的地面沉降,形变速率达到30 mm·a-1以上。通过对事发地的实地调查和形变特征分析,并结合当地地质资料推测了塌陷形成的机理:供水管道下方的软土存在不均匀沉降,使水管产生裂缝导致管道内水外渗,进而致使还未达到胶装凝固点的管片产生裂缝,最终引起隧道和地面坍塌。研究结果可以为今后盾构施工中塌陷的监测和预警工作提供理论依据。     

Monitoring severe aquifer-system compaction and land subsidence in Taiwan using multiple sensors: Yunlin,the southern Choushui River Alluvial Fan

During 1992–2007, excessive pumping of groundwater caused large-scale aquifer-system compaction and land subsidence in the Choshui River Alluvial Fan, especially in the area of Yunlin county. The subsidence impedes surface-water runoff and endangers the operation of Taiwan High Speed Rail. Leveling, Global Positioning System (GPS), multi-level compaction monitoring well, and Differential Interferometric Synthetic Aperture Radar (DInSAR) are used to study the extent of subsidence in Yunlin and its mechanism. These sensors complement each other in spatial and temporal resolutions. A leveling network totaling 434 km in length was deployed to derive subsidence at every 1.5 km along the routes, and the result is accurate to few mm and shows a basin-like subsidence pattern centering at Tuku Township. Four multi-level compaction monitoring wells, co-located with GPS pillars, detect compactions at different depths, showing that the aquifer-system compaction (the cause of subsidence) occurs mostly below depths >200 m, where reduction of groundwater pumping is most needed. The vertical displacements from GPS and leveling agree to within 1 cm, and are larger than the cumulative compaction detected by the compaction-monitoring wells, suggesting that compaction also occurs below 300 m (the depth of the wells). The vertical displacements derived using DInSAR and 8 ENVISAT SAR images agree with the leveling result to 1–2 cm.     

岩溶矿床疏干区地下水位恢复对岩溶塌陷作用机制的研究——以湖南宁乡大成桥为例

截止2016年底,大成桥疏干区在地下水位恢复过程中共发生了17起岩溶塌陷,文章结合地下水动态监测、水文地质钻探成果以及现场调查资料等,研究地下水恢复特征及其对岩溶塌陷的作用机制,为研究区塌陷的预警防治提供科学依据。结果表明:研究区岩溶塌陷集中在河流两岸阶地和麻枣冲沟中;岩溶塌陷发生顺序为先补给区后径流区,与地下水恢复动态过程密切相关;区内地下水恢复过程可分为4个阶段,不同阶段岩溶塌陷机理可分为渗流潜蚀-气爆效应,负压吸蚀效应,软化增荷-挤压-吸蚀效应和软化-崩解效应等4种力学效应组合,不同力学效应作用下的致塌模式有所不同。其中地下水位在基岩面附近4~5 m范围波动更容易引发岩溶塌陷。目前补给区地下水位标高63.95~71.8 m,高于基岩面5~6 m;径流区为42.8~57.7 m,位于基岩面以下4~17 m,因此须密切观测研究区地下水位变化。      

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.     

A seismological study of shallow weak micro-earthquakes in the urban area of Hamburg city,Germany, and its possible relation to salt dissolution

In the night from 8/9 April 2009, shortly after midnight on Maundy Thursday before Easter, several people in Gross Flottbek, Hamburg, felt unusual strong ground shocks so that some of them left their houses in fear of earthquake shaking. Police and fire brigade received phone calls of worried residents. A few days later, Internet pages were published where people reported their observations. On 21 April 2009 at about 8 p.m. local time, a second ground-shaking event was felt. Damage to buildings or infrastructure did not occur to our knowledge. The Institute of Geophysics, University of Hamburg, installed from 22 April to 17 May 2009 three temporal seismic stations in the epicentral area. Seismological data from two nearby stations at the Deutsches Elektronen-Synchrotron at 1 km and the Geophysical Institute at 7 km distance were collected and integrated to the temporal network. The events occurred above the roof of the shallow Othmarschen Langenfelde salt diapir, in an area known for active sinkhole formation and previous historic ground-shaking events. The analysis of the seismological data shows that three shallow micro-earthquakes occurred from 8 to 21 April at a depth of about 100 m, the largest one with a moment magnitude of about M _W 0.6. Depth location of such shallow events is difficult with standard methods and is here constrained by waveform modeling of surface waves. Earthquakes occurring in soft sediments within the uppermost 100 m are a rare phenomena and cannot be explained by standard models. Rupture process in soft sediments differs from those on faults in more competent rock. We discuss the rupture and source mechanism of the events in the context of previous historic shocks and existing sinkhole and deformation data. Although the event was weak, the rupture duration of 0.3 s was unusual long. Three possible models for the generation of repeated ground-shaking events in Gross Flottbek are developed and discussed, implying quit different hazards for subsidence, ground motion, and sinkhole formation. Our favored model postulates that roof failure occurs in an existing soil cavity beneath the epicenter at a depth of about 100 m. Other models bearing a smaller geo-hazard cannot be disproved with the data available, but future geophysical experiments may be planned to resolve this question.     

Results of geodetic and geotechnical monitoring of subsidence for Taiwan High Speed Rail operation

Taiwan High Speed Rail (THSR), which began operations in January 2007, passes through an area in Yunlin County where the largest cumulative subsidence measured during 1992–2006 exceeds 100 cm. Leveling benchmarks, GPS pillars and multi-level monitoring wells were deployed in this area to collect detailed subsidence data from October 2003 to 2006. Leveling is carried out on both ground benchmarks and survey bolts attached to THSR columns. Minimum constraint solutions of leveling networks produce estimated heights accurate to a few mm. Special attention is paid to code smoothing, ionospheric, tropospheric and ocean tidal loading (OTL) effects, so that height estimates from GPS are optimal. Leveling and GPS-derived height changes are consistent to 1 cm, and show that from Stations 210 to 240K of TSHR, the subsidence is bowl shaped. Measurements of sediment compaction in specific depth intervals at three monitoring wells indicate that most of the subsidence is caused by sediment compaction at depths from 50 to 300 m. The major compaction occurs in the interval 220–300 m and is attributed to ground water withdrawal. Large angular deflections as determined from subsidence measurements are detected at some columns, but are below the upper bound (1/1,000) of tolerance specified in the safety code. With the current subsidence and sediment compaction, no significantly reduced loading capacity of the columns is expected to occur. For a safe THSR operation, subsidence and sediment-compaction monitoring should be continued, and current ground water withdrawal in Yunlin must be reduced or stopped.     

Using differential SAR interferometry to map land subsidence: a case study in the Pingtung Plain of SW Taiwan

Synthetic aperture radar (SAR) interferometry (InSAR) is a geodetic tool widely applied in the studies of earth-surface deformation. This technique has the benefits of high spatial resolution and centimetre-scale accuracy. Differential SAR interferometry (DInSAR) is used to measure ground deformation with repeat-pass SAR images. This study applied DInSAR and persistent scatterers InSAR (PSInSAR) for detecting land subsidence in the Pingtung Plain, southern Taiwan, between 1995 and 2000. In recent years, serious land subsidence occurred along coastal regions of Taiwan as a consequence of over-pumping of underground water. Results of this study revealed that the critical subsidence region is located on the coast near the estuary of Linpien River. It is also found that subsidence was significantly higher during the dry season than the wet season. The maximum annual subsidence rate of the dry season is up to −11.51 cm/year in critical subsidence region and the vertical land movement rate is much slower during the wet season. The average subsidence rates in wet and dry seasons are −0.31 and −3.37 cm/year, respectively. As a result, the subsidence rate in dry seasons is about 3 cm larger than in wet seasons.     

Analysis of sinkhole occurrences over abandoned mines using fuzzy reasoning: a case study

In Korea most of the old mine workings were worked with room and pillar method or sublevel caving method and today they possess great possibility of surface subsidence especially for shallow depth mines. In most of the cases, mine roadways, rooms and pillars are irregular in shape and information about the local geology is uncertain. For these reasons, it is difficult to standardize the estimating method of subsidence especially sinkhole type over abandoned mine area. This paper describes the application procedure for the fuzzy reasoning techniques to analyze the possibility of sinkhole occurrences over abandoned mines. This technique is implemented in software which can simplify the analysis procedure and present the possibility of sinkhole subsidence without having precise information about local geological/mining conditions. This technique has been applied to forecast sinkhole possibilities at Bonghwa site where a massive sinkhole has already been occurred.     

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.     

Sinkhole hazard assessment in Minnesota using a decision tree model

An understanding of what influences sinkhole formation and the ability to accurately predict sinkhole hazards is critical to environmental management efforts in the karst lands of southeastern Minnesota. Based on the distribution of distances to the nearest sinkhole, sinkhole density, bedrock geology and depth to bedrock in southeastern Minnesota and northwestern Iowa, a decision tree model has been developed to construct maps of sinkhole probability in Minnesota. The decision tree model was converted as cartographic models and implemented in ArcGIS to create a preliminary sinkhole probability map in Goodhue, Wabasha, Olmsted, Fillmore, and Mower Counties. This model quantifies bedrock geology, depth to bedrock, sinkhole density, and neighborhood effects in southeastern Minnesota but excludes potential controlling factors such as structural control, topographic settings, human activities and land-use. The sinkhole probability map needs to be verified and updated as more sinkholes are mapped and more information about sinkhole formation is obtained.