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.     

Development and validation of sinkhole susceptibility models in mantled karst settings. A case study from the Ebro valley evaporite karst (NE Spain)

A preliminary sinkhole susceptibility analysis has been carried out in a stretch 50 km² in area of the Ebro valley alluvial evaporite karst (NE Spain). A spatial database consisting of a sinkhole layer and 27 thematic layers related to causal factors was constructed and implemented in a GIS. Three types of sinkholes were differentiated on the basis of their markedly different morphometry and geomorphic distribution: large subsidence depressions (24), large collapse sinkholes (23), and small cover-collapse sinkholes (447). The susceptibility models were produced analysing the statistical relationships between the mapped sinkholes and a set of conditioning factors using the Favourability Functions approach. The statistical analyses indicate that the best models are obtained with 6 conditioning factors out of the 27 available ones and that different factors and processes are involved in the generation of each type of sinkhole. The validation of two models by means of a random-split strategy shows that reasonably good predictions on the spatial distribution of future dolines may be produced with this approach; around 75% of the sinkholes of the validation sample occur on the 10% of the pixels with the highest susceptibility and about 45% of the area can be considered as safe.     

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.     

Natural and human-induced sinkholes in gypsum terrain and associated environmental problems in NE Spain

The central Ebro Basin comprises thick evaporite materials whose high solubility produces typically karstic landforms. The sinkhole morphology developed in the overlying alluvium has been studied using gravimetry and ground-penetrating radar (GPR) on stream terraces, as well as analyzing the evolution of sinkhole morphologies observed in aerial photographs taken in 1928, 1957, and 1985. The sinkhole morphologies give some idea of possible subsurface processes as well as an indication of the final mechanisms involve in sinkhole development. On stream terraces and cover pediments the most commonly encountered dolines are bowl-shaped in their morphology with both diffuse and scarped edges. In contrast, dolines developed in the gypsiferous silt infilled valleys have a funnel and well-shaped morphology. The diffuse-edged bowl-shaped dolines are developed through the progressive subsidence of the alluvial cover, due to washing down of alluvial particles through small voids and cracks into deeper subsurface caves, resulting in a decrease alluvial density. Future compaction of the alluvial cover will produce surface subsidences. This type of dolines are associated with negative gravity anomalies. In contrast, the scarped-edge dolines are formed by the sudden collapse of a cavity roof. The cavities and cracks formed in the gypsum karst may migrate to the surface through the alluvial deposits by piping, and they may subsequently collapse. In this instance, the cavities can be detected by both gravity and GPR anomalies where the voids are not deeper than 4–5 m from the surface. These processes forming sinkholes can be enhanced by man-induced changes in the groundwater hydrologic regime by both inflows, due to irrigation, ditch losses, or pipe leakages, and by outflows from pumping activities.     

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.     

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     

Paleosubsidence and active subsidence due to evaporite dissolution in the Zaragoza area (Huerva River valley, NE Spain): processes, spatial distribution and protection measures for transport routes

The lowest 17-km long reach of the Huerva River valley, down to its confluence with the Ebro River in Zaragoza city, flows across salt-bearing evaporites of the Ebro Tertiary Basin (NE Spain). Upstream, the horizontally lying Miocene evaporites are interfingered with non-soluble distal alluvial fan facies (shales and sandstones). The proportion of soluble facies in the Huerva River valley increases in a downstream direction towards the basin depocenter. On the basis of the type and magnitude of the paleosubsidence features, the valley has been divided into four reaches. Along reach I, undeformed terrace deposits less than 4 m thick rest on insoluble detrital bedrock. In reaches II and III, dissolution at the alluvium–bedrock boundary has generated local thickening, deformation and paleocollapse structures, which only affect the alluvial mantle. In reach IV, terrace deposits thicken to over 60 m resulting from a large-scale synsedimentary subsidence. In this sector, subsidence locally affects to both the alluvium and the underlying bedrock. This indicates that dissolution acts at the rockhead beneath the alluvial cover (alluvial karst) and within the evaporitic substratum (interstratal karst). The development of an intraevaporitic karst in reach IV is attributed to gypsum and salt dissolution. Irregular terrace gravel bodies (gravel pockets) embedded in a fine-grained matrix associated with paleocollapse structures have been interpreted as liquefaction–fluidization structures resulting from ground acceleration and suction induced by catastrophic collapses. Subsidence is currently active in the region affecting areas with a thin alluvial cover in reaches III and IV. The low subsidence activity in most of Zaragoza city is explained by the presence of thickened (around 50 m) and indurated alluvial deposits. In the surrounding area, numerous buildings in Cadrete and Santa Fe villages have been severely damaged by subsidence. Natural and human-induced subsidence favours the development of slope movements in the gypsum scarp overlooking Cadrete village. Several transport routes including the Imperial Canal (irrigation canal) and the recently completed Madrid–Barcelona high-speed railway are affected by human-induced sinkholes. The paleocollapse structures exposed in the trenches of this railway and a ring road under construction point to hazardous locations underlain by cavities and collapse structures where special protection measures should be applied. Rigid structures are recommended beneath the high-speed railway with sufficient strength to span the larger sinkholes with no deformation. Electronic monitoring devices linked to a warning system can detect subtle subsidence-induced deformations in carriageways or railways. This research demonstrates that the study of the paleokarst helps to understand the processes involved in the present-day subsidence phenomena and their general spatial distribution.     

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.     

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.     

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.     

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.     

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.     

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.     

Morphometric analysis of sinkholes in a karst coastal area of southern Apulia (Italy)

Salento, the southern portion of Apulia region (SE Italy), is a narrow and elongated peninsula in carbonate rocks, with prevailing low coastlines, locally interrupted by high rock cliffs. The long stretches of low coasts are marked by typical karst landforms consisting of collapse sinkholes. As observed in many other karst coastal settings worldwide, development of sinkholes may be particularly severe along the coasts, where both natural and anthropogenic processes contribute to accelerate the dissolution of carbonate rocks and subsidence processes, even influencing the coastline evolution. Following a previous study, where the main features of sinkholes at Torre Castiglione (Taranto province) were investigated and described, and a preliminary susceptibility map produced, at the light of updated data and elaborations in the present paper we perform a detailed morphometric analysis on the sample of identified sinkholes. The main morphometric parameters generally used for sinkhole characterization have been considered in this study: shape of the sinkhole, azimuth and length of maximum and minimum axes, depth, elongation ratio, and distance from the shorelines. Each of them is described, both as individual parameter and in conjunction with the others, in the attempt to identify the main factors controlling development of sinkholes in the area, and their evolution as well. With regard to this latter aspect, beside simple morphometry of the sample of sinkholes at Torre Castiglione, we also focused our attention on the likely relationships existing between distribution and shape of the sinkholes and the tectonic discontinuities. To investigate the matter, a three-stage analysis has been carried out in this study by means of: field measurements of the fractures bounding the sinkholes, field measurement of the long axes azimuth of the elongated sinkholes, comparison of the previously described sets with the strikes of the main regional geological structures. The obtained results show, in addition to the coincidence of the main regional discontinuity systems with the major axis of elongated sinkholes, a clear control exerted by development and evolution of the sinkholes on the formation of coastal inlets and bays. Eventually, the approach here presented may be applied in other karst coastal sinkhole-prone areas, to gain new knowledge on the genesis and evolution of coastal sinkholes, and to properly evaluate the hazard they pose to the anthropogenic environment.     

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.     

Geomorphological hazards related to deep dissolution phenomena in the Western Italian Alps: Distribution, assessment and interaction with human activities

Deep dissolution affects great part of soluble rocks (e.g. gypsum and anhydrite) of the Western Italian Alps. The related superficial phenomena (sinkholes, gravity-induced processes and a local worsening of geomechanical rock properties) are not limited to typical karsts landscape and cause slope instability also affecting populated sites and infrastructures. The paper aims to describe general characteristic of dissolution phenomena, to interpret their conditioning factors and evolutionary stages and to assess possible hazards due to their superficial effects.The search for evidences of deep dissolution leads to the selection of representative sites in the central part of the Western Italian Alps (Piemonte and Valle d'Aosta Region). Detailed geological and geomorphological studies have been used to classify the selected sites by type, size and variable state of activity. Very different evolutionary stages of dissolution phenomena have been interpreted by comparison of case-studies: some are early “embryonic”; others are more evolved, up to typical sinkholes, or even remodelled by other phenomena. Some cases show an extreme complexity in the interactions between corrosion phenomena and other geomorphic processes: slope deformations, from one side, and karst, fluvial and glacial phenomena, to the other. A wide range of movement rates on slope instabilities induced by deep dissolution have been estimated by topographic and geomorphic data. Geochemical data on removed rocks by dissolution indicate 0.4 mm/year values for local subsidence. Historical and technical data indicate low frequency of major dissolution-induced collapses, but highlight widespread damages to tunnels, roads and buildings, especially along slopes.     

Evaluation of susceptibility for terrain collapse and subsidence in karst areas,municipality of Iraquara,Chapada Diamantina (BA), Brazil

The morphological evolution of the karstic systems is associated with a set of physical and chemical processes, triggered by the dissolution of the rocks, related to percolation of groundwater and surface water, which consequently open underground voids and carve out peculiar forms of relief. Due to environmental and geotechnical aspects, this system is naturally more fragile and vulnerable than other natural systems and, therefore, has increasingly received the attention of the scientific community over the past decades. The objective of the study was to delimit zones with varying degrees of susceptibility for collapses and subsidence of sinkholes in the municipality of Iraquara, Chapada Diamantina (BA), Brazil, and to understand their geological and morphological determinant factors. Geological data, karst phenomenon map, and visual analysis in the field were used to categorize zones with different types of susceptibilities to the nucleation of new sinkholes based on a Hazard Index. This index was defined from the sum of geological hazard factors, lineament density, and sinkhole density. The areas that presented the highest susceptibility for terrain collapse and subsidence corresponded to regions where carbonate rocks outcrop, with high density of photolineaments and 2.62 sinkholes/km². Processes associated with terrain collapse and subsidence in karst areas consisted of a combination of various factors, hindering precise predictions. However, zones of different types of susceptibilities to terrain collapse and subsidence can be delimited when the relationships between these processes and their factors are understood. The Hazard Index proposed does not provide quantitative values for the probability of hazard susceptibility, but rather indicates areas that are more susceptible to terrain subsidence and collapse.     

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     

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.     

Sinkholes related to discontinuous pumping: susceptibility mapping based on geophysical studies. The case of Crestatx (Majorca, Spain)

The Crestatx aquifer is the main source of water supply to the Bay of Alcudia, one of the largest resorts on the island of Majorca (Spain). This water has been used since the 1970s using several pumping wells, which draw an annual volume of 1.5 hm³. The seasonal exploitation of this karstic aquifer causes a substantial cone of depression with great variations in the piezometric level (up to 120 m) and dynamic water levels of down to 87 m below sea level. At the end of the 1990s, several sinkholes and subsidence depressions started being detected in a highly karstified area. Twenty subsidence and sinkhole morphologies have been inventoried in an area measuring 70,000 m², with diameters up to 23 m and depth more than 15 m. The intense and continuous rainfall during recent years (2008–2010) has considerably accelerated the process, increasing the dimensions of existing sinkholes and the appearance of new morphologies. By means of electrical tomography techniques, a ground study was carried out. Numerous cavities have been identified in the upper 30 m. Using a map of the surface morphologies and the geophysical profiles, we can determine that the propagation and orientation of the sinkholes lie along three main directions: N30°E, N130°E and N60°E. The first is the most relevant, which is parallel to the main tectonic structures in the area. The interpretation of the electrical profiles has enabled us to identify the potentially most unstable areas, which is an effective tool to assess risk in the area, as there are roads and a housing development nearby. The high, but discontinuous, exploitation of the aquifer is considered the main trigger for these sinkholes and subsidence depressions, as it causes large variations of pressure and accelerates the dissolution process in the underlying rock.