Salles Et Al-2018-Environmental Earth Sciences
Salles Et Al-2018-Environmental Earth Sciences
Salles Et Al-2018-Environmental Earth Sciences
net/publication/327203787
CITATIONS READS
0 29
6 authors, including:
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Lucas De Queiroz Salles on 24 August 2018.
ORIGINAL ARTICLE
Abstract
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 carbon‑
ate rocks outcrop, with high density of photolineaments and 2.62 sinkholes/km2. 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.
Keywords Karst · Carbonate rocks · Chapada Diamantina · Collapse and subsidence · Brazil
13
Vol.:(0123456789)
593 Page 2 of 11 Environmental Earth Sciences (2018) 77:593
in situations of direct anthropogenic interference (Lei et al. photolineament density and sinkhole density, and summa‑
2001; Galvão et al. 2015). rized in the Hazard Index. This knowledge is paramount
Terrain collapse and subsidence in karst areas are usu‑ for the protection of the physical environment and territo‑
ally influenced by various natural and/or human factors: (1) rial development, which should be applied as a subsidy for
degree of karstification; (2) increase of water input in the elaborating a land use and occupation master plan.
soil (rain water, floods, pipe leakages, irrigation, among oth‑
ers); (3) water level reduction (due to paleoenvironmental
changes, groundwater extraction, and tectonics); (4) water Site description
damming; (5) epikarst erosion and carving; (6) underground
excavations; (7) vegetation removal; (8) distance from faults Iraquara is located in the state of Bahia, Brazil, 340 km west
and fractures; and (9) soil thawing (Hu et al. 2001; He et al. from Salvador, the state capital (Fig. 1a). According to the
2003; Zhou 2006; Galve et al. 2009a, b; Keqiang et al. 2009; last national census, this municipality presents a popula‑
Cooper and Gutiérrez 2013; Gutiérrez et al. 2014; Galvão tion of 25,279 inhabitants (IBGE 2017), distributed over
et al. 2015). Catastrophic collapses directly over cavities an area of 991 km2. Demographic growth was estimated at
occur in less than 1% of cases. In most reports, remobiliza‑ approximately 10% between 2010 and 2016 (IBGE 2016).
tion of unconsolidated cover occurs towards the interior of Urban occupation is mainly observed in the central portion
caves, bulging the surface (Ford and Williams 2007). of the municipality (Fig. 1a). Studies conducted based on
Several cases of collapse and subsidence events in karst satellite images showed that an occupation boom occurred
areas have been reported around the world, for example in the region in 2006 (Fig. 1b).
in Tangshan, China (Hu et al. 2001), in the Ebro Valley, The main economic activities are limited to tomato, ole‑
Zaragoza, Spain (Galve et al. 2009a, b), and in Sango, Ten‑ aginous plants and umbu crops, and to cave tourism. The
nessee, USA (Siska et al. 2016). In Brazil, the municipali‑ region is located within one of the most expressive land‑
ties of Lapão and Itaetê (Bahia), Mairinque and Cajamar scape sites of Brazil, the Chapada Diamantina, where nature
(São Paulo), Sete Lagoas and Lagoa Santa (Minas Gerais), tourism represents an important economic activity. For
Almirante Tamandaré and Bocaiuva do Sul (Paraná), and Auler and Ferrant (1996), the gallery networks that exist in
Teresina (Piauí) are examples of areas where geotechnical the region compose one of the most relevant speleological
processes led to countless environmental and material dam‑ sites in the country, with 99 cavities recorded so far (CANIE
ages (Silva 1998; Pereira 1998; Santos 2008). However, the 2017), or 0.12 caves per km2. This is possibly the location
distribution and probability of occurrence of new collapse with the highest density of underground galleries per unit
and subsidence phenomena represent an important gap in of area in Brazil.
knowledge for most of these municipalities. Geologically, the area is located within the central pan‑
Geological and geomorphological characterization and orama of the state of Bahia. Its evolution model follows
spatialization of information about known events are impor‑ an aulacogen structural pattern (Souza et al. 1993) called
tant actions towards the mitigation of risks associated with Western Chapada Diamantina Thrust and Fold Belt. In this
collapse and subsidence events in karts terrains, to analyze region, Mesoproterozoic siliciclastic lithostratigraphic unit
areas with highest probability of new occurrences (Gutiérrez outcrops are represented, from base to top, by the Tomba‑
et al. 2008). Susceptibility maps represent the probability dor, Caboclo, and Morro do Chapéu formations, composing
of sinkholes either occurring or not in specific areas (Galve the Chapada Diamantina Group, and a thick Neoproterozoic
et al. 2011). The purpose of delimiting zones in these types carbonate–pelite sequence that composes the Una Group
of maps is to express the various levels of nucleation of new (São Francisco Supergroup), where the karst terrains found
sinkholes, which represents an important tool for the man‑ in Iraquara are developed. This sequence lies discordantly
agement of land use and occupation and of water resources upon meta-sediments of the Chapada Diamantina Group
(Hu et al. 2001; Galvão et al. 2015; Siska et al. 2016). and is composed, from base to top, of the Bebedouro and
Founded in 1962, the municipality of Iraquara, which Salitre formations. Covering this set of rocks are outcrops
was the object of the study, is located in a region with high of Cenozoic detrital covers, which are related to successive
density of underground galleries (Auler and Ferrant 1996; erosion cycles that sculpted the relief (Misi 1979; Bonfim
IBGE 2016). The variety of geoforms, associated with disor‑ et al. 1985; Pedreira et al. 1987; Souza et al. 1993, 2002;
derly growth and lack of urban planning, has promoted high Pedreira 2004).
geotechnical risks to the municipality. The objective of the The Salitre Formation is represented by the Nova Amé‑
paper was to delineate zones with varying degrees of hazard rica and Jussara units. The Nova América unit is composed
susceptibility, which result from the collapse and subsidence of calc-siltite, dololutite, dolarenite, and algal mudstone,
of sinkholes, in the municipality of Iraquara, Bahia, Brazil. deposited in an intertidal to subtidal environment with
To do so, a new method was proposed associating lithology, periodical exposure. On the other hand, the Jussara unit is
13
Environmental Earth Sciences (2018) 77:593 Page 3 of 11 593
Fig. 1 a Location of the study area, in geographic coordinates; b urban occupation growth of the municipality of Iraquara, between 1986 and
2016, defined using images from the Landsat satellite; and c digital elevation model defined from an ASTER image, with resolution of 30 m
composed of oolitic calcarenites, which are often oncolitic, by smooth-profile depressions associated with slow subsid‑
deposited in a subtidal environment, and, locally, in an inter‑ ence caused by the removal and infiltration of detrital covers
tidal to supratidal environment (Pedreira 2004). (Fig. 2).
Regarding geomorphology, Iraquara can be divided
into two domains: (1) mountainous region, where meta-
sediments of the Chapada Diamantina Group outcrop, with Materials and methods
topography ranging between 800 and 1200 m, and slope
varying between 8 and 30%; and (2) karst plateau, where Three datasets were used to define the hazard susceptibility
morphological features composed of smoothly sloped and model associated with collapses in Iraquara. These included
wavy ramps predominate, with slopes varying between 0 geological data, a karst phenomenon map, and visual analy‑
and 12% (Fig. 1c). sis in the field. Lithology, lineament density, and karst fea‑
Karst depressions correspond to a set of sinkholes or ture density maps were elaborated from these data. Maps
uvalas that are concentrated along stretches with NW–SE were later integrated using Hazard Index map algebra. These
preferential orientation and, subordinately, N–S orientation, methods are described as follows.
occupying 1.64% of the total area of the municipality. In this
region, karst depressions can be classified into two groups Elaboration of basemaps
(Cruz Junior 1998): (1) collapse sinkholes, which corre‑
spond to enclosed depressions, with steep profiles; and (2) Geological data comprised: (1) lithological map (1:100,000),
suffosion sinkholes (Ford and Williams 1989), represented elaborated from the IBGE/SEI (2003) database; and (2)
13
593 Page 4 of 11 Environmental Earth Sciences (2018) 77:593
Fig. 2 a Collapse sinkhole possibly associated with overexploitation of groundwater and use in tomato crops (See detail in a1); since the col‑
lapse an installation and exploration of a tubular well 200 m from the site occured; b suffosion sinkhole (Ford and Williams 2007)
lineament density map (Fig. 3a), produced using the “Line high, for values above the mean, plus standard deviation; and
Density” tool of ArcMap 10.1 software. Lineaments were (3) medium, for values between intervals 1 and 2 (Fig. 3).
delimited by aligned crests, valley bottoms, and drainages,
according to the interpretation of shadowed relief images Hazard model associated with terrain collapse
which were produced based on azimuths 0°, 90°, 180°, and and subsidence in karst areas
270°.
A karst phenomenon map was developed to build a broad The purpose of this model was to divide the study area into
cartographic inventory. The map was built on the extraction sub-regions with various collapse or subsidence hazard sus‑
of geoforms from aerial photography analysis at a scale of ceptibilities. Three factors were used to create the hazard
1:60,000. Aerial photographs were then scanned and vec‑ model based on Galve et al. (2008), Hu et al. (2001), Galvão
tored in the software ArcMap 10.1. The sinkhole density et al. (2015), and Siska et al. (2016): (1) lithological; (2)
map consisted of an interpolation of central points (unit of fracturing density; and (3) sinkhole density. Hazard factors
sinkhole) using the “Density” tool of the ArcMap 10.1 soft‑ and their attributes are summarized in Fig. 4, according to
ware (Fig. 3b). their degree of influence on hazard susceptibility. The Haz‑
A class distribution analysis was conducted to define ard Index (HI) comprises the sum of the previously men‑
patterns regarding sinkhole and fracturing indexes. This tioned factors (Fig. 4).
analysis consisted of extracting the values of the previously The lithological hazard factor was based on the rocks
mentioned indexes from an equidistant sampling mesh of mapped in the field and their relationships with the num‑
150 m for lineage density (Fig. 3a1), and 85 m for the sink‑ ber of occurrences of sinkholes per square kilometer. Thus,
hole density map (Fig. 3b1). These values were taken from three categories were defined regarding this hazard factor:
the extraction of points from the lowest cells of raster images (1) areas where the outcrop carbonate rocks from the Salitre
(pixel) using the “Conversions Tools” of the ArcMap 10.1 Formation received a hazard score of “10”; (2) carbonate
software. Indexes were divided into three classes: (1) low, rocks covered by unconsolidated sediments that received a
for values below the mean, minus standard deviation (σ); (2) score of “5”; and (3) siliciclastic rocks from the Chapada
13
Environmental Earth Sciences (2018) 77:593 Page 5 of 11 593
Fig. 3 a Photolineament density map extracted from a raster image; elaborated from aerial photographs at a scale of 1:60,000; and (b1)
(a1) class distribution histogram regarding lineament density and class distribution histogram regarding sinkhole density and relation‑
relationship with the degree of Hazard Index; b sinkhole density map, ship with the degree of Hazard Index
Diamantina Group and Bebedouro Formation that received degree of terrain surface instability and also favors the coa‑
a hazard score of “0”. lescence of these superficial karst features, forming uvalas in
The same principle was used to quantify the hazard factor areas where the number of sinkholes per unit of area is high.
of lineament density. Therefore, three classes were deter‑ The scores calculated based on the HI (Fig. 4) resulted in
mined regarding the fracturing density hazard factor (Fig. 4). dimensionless values. Thus, a descriptive statistical analy‑
The hazard factor associated with sinkhole density is sum‑ sis was conducted to evaluate the empiric distribution of
marized in three categories: high, medium, and low, attribut‑ data and establish hazard intervals, based on all possible
ing scores of 1, 2, and 4, respectively. This factor reflects the combinations of the degree of hazard for terrain collapse
13
593 Page 6 of 11 Environmental Earth Sciences (2018) 77:593
and subsidence, calculated through the HI equation (Fig. 4). The carbonate rocks correspond to limestone and dolo‑
Quartiles (the separatrix measurement that divides the data mite of the Salitre Formation are located in the central por‑
sampled into four different groups—each one presenting tion of Iraquara, occupying approximately 36.6% of the area
25% of data) of the possible HI scores were divided on an (IBGE/SEI 2013). Petrographically, the carbonate rocks of
Excel spreadsheet. The relationships between the calculated the region present an extremely low primary porosity. How‑
quartiles and the HI associated with collapse or subsidence ever, when present, porosity is associated with the process
of the karst relief were thus determined. of calcite dolomitization. The effective permeability of the
aquifer is limited to secondary porosity (faults and fractures)
and the formation of tertiary porosity (cavities) (White 1999,
2002).
Results
Recent unconsolidated sediments cover these units. This
cover is indiscriminately distributed over the municipal‑
Results consisted of a geologic map, factors associated with
ity, occupying approximately 19.6% of its area (IBGE/SEI
terrain collapse and subsidence, and the Hazard Index.
2013).
Fig. 5 Geologic map of the municipality of Iraquara. The silici‑ ate rocks of the Salitre Formation outcrop in the central portion of
clastic rocks of the Bebedouro Formation and Chapada Diamantina the area; Cenozoic detrital covers occur indiscriminately (IBGE/SEI
Group occur in the W and E portions of the municipality; the carbon‑ 2013)
13
Environmental Earth Sciences (2018) 77:593 Page 7 of 11 593
Fig. 6 a Cryptoalgal mudstone of the Nova América Formation, mudstone detailing the contact of the OM-rick mudstone overlapped
deposited in a supratidal environment; b petrographic slide, under by carbonate mud; d detail of micrite dolomitization, under plane
plane polarized light and ×10 objective, showing levels of rock disso‑ polarized light and ×20 objective
lution (in black) enriched in organic matter (OM); c dolomitized algal
over the siliciclastic rocks of the Chapada Diamantina Group The geologic hazard factor associated with sinkhole den‑
and Bebedouro Formation. sity (Fig. 7c) showed two regions that presented the highest
Photolineaments possibly consist of foliations, fractures, concentrations of sinkholes in Iraquara. The first was located
faults and interactions between stratification planes and the in the central–southern portion of the municipality, while the
surface. These relief features can be partially masked by the second was in the central–northern portion, where part of
presence of sedimentary covers (Cheema and Islan 1994). the town is located. In most karst areas, the evolution of the
The lineament density map (Fig. 7b) showed that they were system is directly related to superficial processes, and the
heterogeneously distributed over the area of the municipal‑ characterization of exokarsts is one of the main resources to
ity. However, the largest unit of high photolineament density understand underground karsts (Ford and Williams 2007).
was observed in the southern portion of the study area. Thus, a larger network of surface sinkholes suggests more
The unit of low lineament density occupied 17.9% of the evolved subterraneous karsts and, therefore, higher suscep‑
area susceptible to surface bulging. This region presented tibility to collapse and subsidence in these areas (Zhou et al.
0.35 sinkhole unit per k m2. In turn, the moderate photo‑ 2003; Gao et al. 2005; Kemmerly 2006).
lineament density unit occupied approximately 63.8% of the
area subject to terrain collapse and subsidence, with 0.67 Hazard Index associated with terrain collapse
sinkhole unit per km2. Finally, the unit of high lineament and subsidence
density occupied approximately 18.3% of the area subject
to collapse and subsidence, with 2.4 sinkhole units per km2 The Hazard Index (Fig. 8) was defined through the sum
mapped (Fig. 7b). This analysis demonstrated a direct cor‑ of the geologic hazard factors (Fig. 7a), lineament density
relation in the region between photolineament density and (Fig. 7b), and sinkhole density (Fig. 7c). The area that
the presence of sinkholes. was classified as having high hazard susceptibility (in red)
13
593 Page 8 of 11 Environmental Earth Sciences (2018) 77:593
Discussion
13
Environmental Earth Sciences (2018) 77:593 Page 9 of 11 593
Fig. 8 Distribution of classes of the risk model, sinkholes, and the moderate risk associated with collapse and subsidence of the area.
urbanized area of the municipality. The town area, located in the cen‑ However, the settlements of Iraquara are located in areas of low to no
tral portion of the map, is established on units that presented high to risk associated with terrain collapse and subsidence
the intersection between fault planes, fractures and rock bed‑ drainage network (Gao et al. 2005; Kemmerly 2006). This
ding planes and the terrain surface—interpreted as local dis‑ would consequently lead to a higher percentage of empty
continuities of the rocky massif. The presence of sinkholes spaces in the rocky massif. This fact presents a close cause
was favored in regions that presented high photolineament and effect relationship with photolineament densities and
densities possibly because of the association with higher evolution of kart systems, increasing the degree of risk asso‑
permeability and secondary porosity of rocks. Thus, authi‑ ciated with terrain collapse and subsidence in karst areas.
genic acidic water percolation is favored in these regions, The data listed and discussed in the study allow the
accelerating the process of void opening and rock karstifi‑ understanding that the nucleation of surface collapses in the
cation. A similar fact was observed by Hu et al. (2001) in municipality of Iraquara may also be associated with: (1) the
the municipality of Tangshan (China), where these authors occurrence of subsurface discontinuities; (2) carbonate rock
observed high probabilities for the nucleation of new sink‑ dissolution; (3) degree of karst evolution and underground
holes near active fault zones. voids; and (4) the presence of sedimentary cover over car‑
In addition, high rock permeability may facilitate, over bonate rocks; locally, the process may also be associated
time, the lixiviation and transport of unconsolidated mate‑ with (5) groundwater exploitation. The risk associated with
rial over rocks towards the interior of karst ducts, forming groundwater removal is not static. It can vary depending
depressions (suffosion sinkholes) on the relief. This process on the level of exploitation of aquifers, and the number and
occurs gradually and slowly, but can progress faster due to proximity of wells. Therefore, poor water resource manage‑
human interference. However, even when the process occurs ment can aggravate the risk associated with the formation of
slowly and gradually, it can cause countless material losses, collapse or subsidence sinkholes, since it can either reduce
as observed in the reports by Pereira (1998), Silva (1998), the support of walls inside rocky massifs or increase the
Galve et al. (2008), and Santos (2008). gradient and ability to transport sediments towards the inte‑
A close relationship was observed between index factors rior of voids.
associated with sinkhole density (Fig. 7c) and photolinea‑ This integrated and low-cost analysis produced satisfac‑
ment density (Fig. 7b). Areas that presented high sinkhole tory results for the comprehension and spatialization of ter‑
density could be related to regions with higher development rain collapse and subsidence in the municipality of Iraquara.
of tertiary porosity, forming a more evolved underground These analyses can be extrapolated to other karst regions
13
593 Page 10 of 11 Environmental Earth Sciences (2018) 77:593
where there is low groundwater exploitation. The disadvan‑ Bonfim LFC, Rocha AJD, Pedreira AJ, Morais Filho JC, Guimarães
tages of this approach consist of not considering the relation‑ JT, Tesch NA (1985) Projeto Bacia de Irecê. Relatório final.
CPRM. Três volumes. Contrato CPRM/SME/CPM
ship between sinkhole nucleation and groundwater removal CANIE—Cadastro Nacional de Informações Espeleológicas (2017)
and soil thickness over carbonate rocks. http://www.icmbio.gov.br/cecav. Accessed 09 Feb 2017
Cheema TJ, Islan MR (1994) Comparison of cave passage ways
with fracture traces and joints in the black hills region of South
Dakota. Natl Speleol Soc Bull 56:96–103
Conclusion Cooper AH, Gutiérrez F (2013) Dealing with gypsum karst prob‑
lems: hazards, environmental issues and planning. Treat is eon
A Risk Index and the zoning of areas that present potential Geomorphology—Karst Geomorphol 6:451–462
Cruz FW (1998) Aspectos geomorfológicos e geoespeleologia do
for terrain collapse and subsidence were proposed based on carste da região de Iquara, centro norte da Chapada Diamantina.
the analysis of geological, lineage density, and sinkhole den‑ Instituto de Geociências, Universidade de São Paulo, Estado da
sity data in the municipality of Iraquara. The results obtained Bahia. São Paulo
were effective regarding the comprehension and spatializa‑ Elhatip H (1995) The influence of karst features on environmental
studies in Turkey. Environ Geol 31(1/2):27–33
tion of risk areas. This information can support territorial Ford DC, Williams PW (2007) Karst hydrogeology and geomorphol‑
planning and preventive urban planning, indicating high- ogy, 2nd edn. Wiley, Geography
priority zones to implement corrective measures, delimiting Ford DC, Williams PW (1989) Karst hydrogeology and geomorphol‑
land use, occupation, and preservation in Iraquara. ogy, 1st edn. Wiley, Geography
Galvão P, Halihan T, Hirata R (2015) Evaluating karst geotechnical
The susceptibility associated with terrain collapse and risk in the urbanized area os Sete Lagoas, Minas Gerais, Brazil.
subsidence phenomenon is understood in the study as the Hydrogeol J 22:1499–1513
result of the combination of various factors (geological, Galve JP, Bonachea J, Remondo J (2008) Development and valida‑
structural, karst characteristics, groundwater circulation and tion of sinkhole susceptibility models in mantled karst settings.
A case study from the Ebro Valley evaporite Karst (NE Sapin).
use), which makes it difficult to securely anticipate where Eng Geol 99:185–197
and when terrain collapse or subsidence will occur. How‑ Galve JP, Gutiérrez F, Lucha P, Guerrero J, Bonachea J, Remondo J,
ever, when the relationships among these factors are under‑ Cendrero A (2009a) Probabilistic sinkhole modeling for hazard
stood, it is possible to delimit zones with varying levels of assessment. Erath Surface Process Landforms 34:437–452
Galve JP, Gutiérrez F, Remondo J, Bonachea J, Lucha P, Cendrero A
susceptibility to surface collapsing. (2009b) Evaluating and comparing methods of sinkhole suscep‑
Satellite image analysis showed that the demographic tibility mapping in the Ebro Valley evaporate karst (NE Spain).
boom in the municipality of Iraquara occurred recently, Geomorphology 11:160–172
possibly beginning in 2006. An increase in groundwater Galve JP, Remondo J, Gutiérrez F (2011) Improving sinkhole hazard
models incorporating magnitude—frequency relationships and
resource exploitation in the region may locally accelerate nearest neighbor analysis. Geomorphology 134:157–170
the processes of terrain collapse or subsidence. Therefore, Gao Y, Alexander EC, Barnes RJ (2005) Karst database implemen‑
it is paramount to elaborate a management plan for ground‑ tation in Minnesota: analysis of sinkhole distribution. Environ
water use. Geol 47:1083–1098
Goldsheider N, Drew D (2007) Methods in karst hydrology. Taylor
Based on the results obtained, some actions can be taken & Francis Group, London
to avoid losses and damages related to surface bulging: (1) Gutiérrez F, Guerrero J, Lucha P (2008) Quantitative sinkhole hazard
avoid new constructions in areas that present high risk; (2) assessment: a case study from the Ebro Valley evaporate alluvial
control the number of tubular wells dug in areas that present karst (NE Spain). Nat Hazards 45:211–233
Gutiérrez F, Parise M, De Waele J, Jourde H (2014) A review on
moderate to high risk; (3) control large-scale irrigation to natural and human-induced geohazards and impacts in karst.
reduce the input of extra water in the system; and (4) pro‑ Earth Sci Rev 138:61–88
mote educational workshops to raise public awareness of the He K, Liu C, Wang S (2003) Karst collapse related to over-pumping
risks associated with terrain collapse and subsidence and and a criterion for its stability. Environ Geol 43:720–724
He H, Jia Y, CHEN W, Wang R, Luo H (2004) Evaluation of karst
how to act in case of catastrophic events. collapse risks induced by over-pumping and karst groundwater
resource protection in Zaozhuang region, China. Environ Earth
Acknowledgements This word was funded by the National Council for Sci 71:3443–3454
Scientific and Technological Development—CNPq. Special thanks go Hu RL, Yeung MR, Lee CF, Wang SJ, Xiang JX (2001) Regional
to Professor Cícero Paixão, for providing the petrographic slides and risk assessment of karst collapse in Tangshan, China. Environ
to the members of ther Núcleo de Estudos Hidrogeológicos e do Meio Geol 40:1377–1389
Ambiente—NEHMA. Huggett RJ (2007) Fundamentals of geomorphology. Routledge Fun‑
damentals of Physical Geography, 2nd edn, p. 483
IBGE/SEI (2003) Cartas Plani-Altimétricas do Estado da Bahia.
Escala 1:100.000
References IBGE/SEI (2013) Instituto Brasileiro de Geografia e Estatística
[Brasilian Institute of Geography an Statistics] and Sistema
Auler AS, Farrant AR (1996) A brief introduction to karst and caves Eletrônico de Informações (SEI) [Electronic Information
in Brazil. Univer Britl Speloca Soc n 20:187–200
13
Environmental Earth Sciences (2018) 77:593 Page 11 of 11 593
System] Geological Map. Available in http://www.portalseib Sallun Filho W (2012) Subsidencia e colapso em terrenos cársticos.
ahia.saeb.ba.gov.br/. Accessed Mar 2017 Desastres naturais: conhecer para prevenir. Instituto Geológico
IBGE-Instituto Brasileiro de Geografia e Estatística (2016) http://www. de São Paulo
ibge.gov.br/home. Accessed 08 Feb 2017 Santos AR (2008) Terrenos Cárstico: Áreas de risco geológico para a
IBGE (2017) Instituto Brasileiro de Geografia e Estatística [Brasilian Engenharia e para o meio ambiente. Available in http://notivias.
Institute of Geography anStatistics] Basic municipal information. ambientebrasil.com.br. Accessed 13 Dec 2016
Available in http://www.cidades.ibge.gov.br. Accessed Nov 2017 Silva AB (1998) Abatimento de solo na cidade de Sete Lagoas, Minas
Kemmerly PR (2006) Modeling doline populations with logistic growth Gerais, Brasil. Revista Brasileira de Água Subterrânea 12:57–66
functions. Earth Surf Process Landf 32:587–601 Siska PP, Goovaers P, Hung IK (2016) Evaluating susceptibility of
Keqiang H, Guangming Y, Yaoru L (2009) Paleo-karst collapse pillars Karst dolines (sinkholes) for collapse in Sango, Tennessee, USA.
in northern China and their damage to the geological environmen‑ Prog Phys Geogr 40:579–597
tal. Environ Geol 58:1029–1040 Souza SL, Brito PCR, Silva RWS (1993) Estratigrafia, Sedimentologia
Lei M, Jiang X, Li I (2001) New advances of karst collapse research in e Recursos Minerais da Formação Salitre na Bacia de Irecê, Bahia
China. Geotechinal na Environmental Applications of karst Geol‑ (Série Arquivos Abertos 2). CBPM, Salvador
ogy and Hydrology, Balkema, Lisse, pp 51–145 Souza SL, Toledo LAA, Brito PCR, Fróes RJB, Silva RWS (2002)
Misi A (1979) Grupo Bambuí no Estado da Bahia. Geologia e Recursos Análise faciológica e metalogenética da bacia de Irecê, Bahia.
Minerais do Estado da Bahia. Salvador SME/CPRM V1:119–154 Salvador, CBPM. Relatório Final, Companhia Baiana de Pesquisa
Pedreira AJ (2004) O Supergrupo Espinhaço na Chapada Diamantina Mineral
Centro-Oriental, Bahia: Sedimentologia, Estratigrafia e Tectônica. White WB (1999) Conceptual models for karstic aquifers. In: Palmer
Tese de Doutorado, Inst. de Geociências, Universidade de São AN, Palmer MV, Sasowsky ID (eds), Karst modeling: Karst
Paulo Waters Institute Special Publication
Pedreira AJ, Rocha AJD, Costa IVG, Moraes Filho JC (1987) Projeto White WB (2002) Karst hydrology: recent developments and open
Bacia de Irecê II. Relatório final. CPRM, Salvador questions. Eng Geol 65:85–105
de Pereira RGFA (1998) Caracterização geomorfológica e geoespele‑ Zhou W (2006) Drainage and flooding in karst terranes. Environ Geol.
ológica do Carste da Bacia do Rio Una, borda Leste da Chapada https://doi.org/10.1007/s00254-006-0365-3
Diamantina (Município de Itaetê. Estado da Bahia). Dissertação Zhou W, Beck BF, Adams A (2003) Application of matrix analysis in
de Mestrado, Inst. de Geociências, Universidade de São Paulo delineating sinkhole risk areas along highway (I-70 near Freder‑
Salles LQ (2017) Hidrogeologia e risco geológico em rochas carbonáti‑ ick, Maryland). Environ Geol 44:834–842
cas proterozoicas: porção central da Chapada Diamantina, Bahia,
Brasil. Masters dissertation. Universidade Federal da Bahia
13