Aires Etal 2018 BarragemMariana Landusechange
Aires Etal 2018 BarragemMariana Landusechange
Aires Etal 2018 BarragemMariana Landusechange
discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/320947381
CITATIONS READS
0 87
5 authors, including:
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Changes in land use and land cover as a result of the failure of a mining tailings dam in Mariana, MG,
Brazil View project
All content following this page was uploaded by Uilson Ricardo Venâncio Aires on 09 November 2017.
Changes in land use and land cover as a result of the failure of a mining MARK
tailings dam in Mariana, MG, Brazil
⁎
Uilson Ricardo Venâncio Airesa, Bismarck Soares Matos Santosb, Clívia Dias Coelhoa, ,
Demetrius David da Silvaa, Maria Lúcia Calijurib
a
Universidade Federal de Viçosa, Departamento de Engenharia Agrícola, Campus Universitário, Viçosa, MG, 36570-000, Brazil
b
Universidade Federal de Viçosa, Departamento de Engenharia Civil, Campus Universitário, Viçosa, MG, 36570-000, Brazil
A R T I C L E I N F O A B S T R A C T
Keywords: The failure of the mining tailings dam in Mariana, Minas Gerais, Brazil, was the largest disaster of its type in
Environmental disaster Brazil, generating serious environmental and socioeconomic problems. The aim of this study was to identify
Geographic information system potential sites for the disposal of the tailings that accumulated along the watercourses affected by the dam failure
Multicriteria analysis in the municipalities of Mariana and Barra Longa, as well as to quantify the area affected by the tailings. To
quantify the area affected by the tailings dam failure, the Land Change Modeler (LCM) tool was used to verify the
changes in land use and land cover that occurred in the periods before and after the environmental impact. For
the identification of the suitable sites for tailings disposal, multicriteria analysis was performed using the Idrisi
Selva® software and considering three scenarios to identify potentially affected areas through the application of
the weighted linear combination (WLC) and ordered weighted average (OWA) techniques. In scenario 1, the
WLC technique was applied, considering average risk and high compensation. In scenarios 2 and 3, the OWA
technique was applied, with high risk and high compensation for scenario 2 and low risk and high compensation
for scenario 3. The highest suitability value for tailings disposal was identified in scenario 2, which obtained the
maximum value of 213 on a scale of 0–255, while for the scenarios 1 and 3, the suitability values were 158 and
124, respectively. The largest site for the tailings dam identified in the analysis for the municipally of Mariana
has a size of 0.25 km2 and obtained a suitability value of 207. In the municipally of Barra Longa, the largest site
identified has a size of 0.037 km2, with maximum suitability value of 209. In relation to the changes in land use
and cover, it was verified that the vegetated areas, composed mainly of riparian forest and pasture, presented
greater losses as a consequence of the dam failure, with reductions of 11.99 and 4.81 km2, respectively, de-
monstrating that the dam rupture disaster has caused serious damage to the environment and the population
living in the region.
⁎
Corresponding author.
E-mail address: [email protected] (C.D. Coelho).
http://dx.doi.org/10.1016/j.landusepol.2017.10.026
Received 24 April 2017; Received in revised form 1 September 2017; Accepted 13 October 2017
0264-8377/ © 2017 Published by Elsevier Ltd.
U.R.V. Aires et al. Land Use Policy 70 (2018) 63–70
Fig. 1. Location of the area affected by the tailings dam failure used in this study.
water and the death of several aquatic species (Davies, 2002). A major life, increasing biochemical oxygen demand, and contributing to diffuse
disaster of this type also occurred in 2003 in the Republic of Macedonia pollution and contamination by pathogens (Bilotta and Brazier, 2008;
(Europe), releasing approximately 100,000 m3 of damaging tailings Stutter et al., 2017).
with elevated amounts of heavy metals (Ag, As, Cd, Cu, Mo, Pb, Sb and Aiming toward the environmental recovery of the affected areas, it
Zn) through the Kamenica River valley, which affected the water is crucial to remove the mining tailings that have been accumulating
quality of Lake Kalimanci (Vrhovnik et al., 2013; Rico et al., 2008). along the river banks and other affected areas and ensure their ade-
There are historical records of environmental disasters involving quate disposal so that this material is not carried to the watercourses,
dams or tailings pile failure in Brazil; however, none of them had the mainly through surface runoff in the rainy season. In this context,
same impacts as the Fundão dam failure on November 5th, 2015, which multicriteria analysis is a tool that can help in the decision-making
belongs to the Germano mining complex in the municipality of process in identifying potential sites for tailings disposal, since it takes
Mariana, MG and was the largest environmental disaster of this type in into account several criteria that delimit the most appropriate areas.
the country (MMA, 2016; ANA, 2016). Such analyses have been used for several purposes, such as the identi-
The dam contained 50 million m3 of iron tailings, of which 34 fication of sites for the implementation of sanitary landfills, which takes
million m3 were released into the environment. This amount reached into account criteria that are similar to those of tailings disposal
663 km of rivers and streams in the Rio Doce river basin in Minas Gerais (Demesouka et al., 2014; Ersoy and Bulut, 2009; Felicori et al., 2016;
and Espírito Santo states (MMA, 2016). Approximately 80 species of Gbanie et al., 2013).
fish are identified as native to this basin, of which 11 are considered Multicriteria analysis is a mathematical tool used to investigate a
endangered and 12 are considered endemic to the affected region series of alternative answers to a given question based on competing or
(IBAMA, 2015). The endemic species are the most threatened because conflicting priorities, providing a structured method for comparative
their development is dependent on the natural characteristics of the analysis and assisting in the decision-making process (Martinelli et al.,
environment that they inhabit (Burlakova et al., 2011). 2014). In this analysis, the criteria are divided into constraints or fac-
In addition to environmental damage, the failure of the dam also led tors, of which the constraints allow the imposition of absolute limita-
to severe socioeconomic losses. The district of Bento Rodrigues, one of tions to a given geographic space according to Boolean logic, that is,
the most affected areas, situated 5 km downstream from the dam, was values of 1 or 0 represent yes or no, respectively, to identify appropriate
practically buried by water, mud and debris. Nineteen fatalities were areas and eliminate unfit areas. Factors are relative limitations that
identified, with 600 people living at the site (ANA, 2016), and an es- determine the degree of fitness of a region of interest. They are con-
timated 207 of the 251 properties in the area were destroyed (IBAMA, tinuous in nature and represent the variation in aptitude within an
2015). established range (Campos et al., 2014; Felicori et al., 2016).
The failure of the dam also resulted in impacts to water quality, In this context, the aim of this work was to apply multicriteria
which impaired the public water supply in municipalities and for sev- analysis for the definition of potential areas that can be used for the
eral other water users around the affected area. Other uses of water disposal of mining tailings resulting from the rupture of the Fundão
resources, such as electric power generation, industrial use, leisure and dam, which accumulated on the river banks of the Mariana and Barra
fishing, were also harmed (ANA, 2016). For example, at the water Longa municipalities in the state of Minas Gerais, Brazil. In addition, it
quality monitoring station in Linhares, in the state of Espírito Santo, the is intended to verify the extent of the area affected by the tailings
observed turbidity value was on the order of 2500 NTU (CPRM, 2015). through comparison in relation to the land use and cover prior to the
A high concentration of suspended sediments in watercourses sig- environmental impact.
nificantly reduces light penetration, affecting the dynamics of aquatic
64
U.R.V. Aires et al. Land Use Policy 70 (2018) 63–70
2. Methods 2.3. Factors and criteria used to identify areas suitable for tailings disposal
2.1. Study area Elevation data were extracted from the digital elevation model
(DEM), with a spatial resolution of 30 m, generated from satellite in-
The work was carried out in the municipalities of Mariana and Barra formation (ASTER sensor) from the Shuttle Radar Topography Mission
Longa, state of Minas Gerais, Brazil, covering a total area of 1578 km2 (SRTM) (USGS, 2016). The DEM was used to obtain the information
(Fig. 1). The Fundão dam was built in the subdistrict of Bento Ro- related to the slope and hydrography of the study area.
drigues, located 35 km from the center of Mariana city and 124 km To obtain the hydrography of the area, the software MapWindow
from Belo Horizonte city. The second affected municipality was Barra was used, which enabled the generation of the hydrologic conditioned
Longa, in the Zona da Mata region, 172 km from Belo Horizonte. digital elevation model (HCDEM) through the IPH – Hydrotools tool.
This procedure was necessary to validate the hydrography extracted
from the DEM (more detailed) with the official Brazilian hydrography.
2.2. Changes in land use and land cover due to the failure of the tailings The dates related to the vegetation information, the geological
dam faults, the soil maps and the conservation units were provided by the
Spatial Data Infrastructure of Minas Gerais state (IEDE, 2016) at a scale
The land use and land cover in the region affected by the tailings of 1:50,000.
dam failure were obtained from Landsat 8 sensor OLI (Operational Land The road and highway database was extracted from the TrackSourse
Imager) images with a spatial resolution of 30 m, provided by the Project (2016) and compared to Google Earth Pro® images to confirm
United States Geological Survey (USGS, 2016) (https://earthexplorer. the locations of the identified roads and highways.
usgs.gov/). In this way, the definition of the potential areas for tailings disposal
An image from October 2015 was used as a representation of the was conducted based on the combination of information regarding
area prior to the environmental impact (which occurred on November paved roads, hydrography, slope, urban areas, geological faults, soils,
5th, 2015), and an image from November 2015 was used as a re- vegetation and conservation units.
presentation of the area after the environmental impact. Road information was used to determine which areas were most
To improve the resolution of the Landsat images, the multispectral suitable for tailings disposal in order to reduce transport costs.
bands were fused with the Panchromatic band to obtain the final re- Consequently, the areas closer to the roads resulted in the greatest
solution of 15 m. This procedure allows the attainment of more detailed suitability in the analysis results.
information about the study area. In the case of hydrography, the restriction of potential areas to a
For the classification of the land use and land cover, the maximum minimum distance of 200 m (buffer) from watercourses according to
likelihood supervised classification method was applied. Six classes of the Brazilian Standard (NBR 13896–Solid Wastes Landfill – Project,
land use and land cover were identified: Rock Outcropping, Pastures, Installation and Operation Criteria Procedure) was adopted as a cri-
Water Resources, Tailings, Vegetation and Urban Area, which were the terion in order to reduce the risks of these materials again occurring in
most visible in the selected images and of the greatest relevance for this the rivers. In this way, the further away from these areas, the greater
study. the suitability of the results obtained (ABNT, 1997).
The Land Change Modeler (LCM) tool in the Idrisi Selva® 17.0 In terms of slope values, potential sites were restricted to sites with a
software was used to verify the changes in land use and land cover due maximum slope of 30% in accordance with NBR 13896. Consequently,
to the failure of the tailings dam. To verify the accuracy of the classi- the suitability of the sites increased with a decrease in slope (ABNT,
fication, the Kappa index (K) was calculated (Eq. (1)). The samples 1997).
considered in the assembly of the confusion matrix were verified with A buffer of 500 m was established around the urban area for tailings
Google Earth Pro®, disposal, since this distance is also required for the construction of solid
waste landfills (nonhazardous) (ABNT, 1997) to avoid public health
y y
problems through the inhalation of fine particles of these materials that
n. ∑ x ii − ∑ (x i +. x+i )
can be carried by the wind.
i=1 i=1
K= y According to the Brazilian Standard (NBR 8419–Introduction of
n2 − ∑ (x i +. x+i ) Projects for The Sanitary Filling of Urban Solid Wastes), it is re-
i=1 (1)
commended that the construction of sanitary landfills includes a buffer
y of 200 m from geological faults; thus, this distance was also adopted as
where ∑ x ii is the sum of the diagonal of the confusion matrix, n is the a restriction in the definition of potential areas for the disposal of
i=1 mining tailings (ABNT, 1992).
total number of samples representing the soil use and cover observed in
The NBR 8419 recommends the implementation of the sanitary
the field, and xi+ and x+i are the sum of the rows and columns of the
filling of urban solid wastes only in areas with a permeability of less
confusion matrix, respectively. The interpretation of the values ob-
than 10−6 cm s−1 and an unsaturated zone with a thickness greater
tained for K were made based on the classification proposed by Landis
than 3 m. There is no information about the groundwater depth in the
and Koch (1977) according to Table 1.
study area; for this reason, it was decided to use sites with latosol soil,
which are deeper, with previous compaction work to reduce soil per-
Table 1
Interpretation of K values for the classification of the land use and land
meability.
cover. For vegetation (native or planted) and conservation units, it was
Source: Landis and Koch (1977). determined that the mining tailings disposal areas should be located
outside of vegetation areas, aiming toward the preservation of the en-
Values of Kappa Interpretation
vironment.
<0 No agreement Table 2 presents, in a synthetic way, the criteria established for the
0.00–0.19 Poor agreement analysis.
0.20–0.39 Fair agreement In the comparison of criteria at different scales, it is necessary to
0.40–0.59 Moderate agreement
normalize the dataset through pertinence functions, which obey fuzzy
0.60–0.79 Substantial agreement
0.80–1.00 Almost perfect agreement logic, converting the data to a scale that ranges from 0 to 1. In this
paper, linear functions were used, as the authors consider that these
65
U.R.V. Aires et al. Land Use Policy 70 (2018) 63–70
Table 2 combination of results from the analysis, so that the site with the best
Variables considered in the analysis and the respective criteria and constraints. suitability for tailings disposal has the highest score (Blachowski, 2015;
Source: Adapted from ABNT/NBR 13896: 1997–Solid Wastes Landfill: Project,
Romano et al., 2015).
Installation and Operation Criteria Procedure and ABNT/NBR 8419: 1992–Introduction
of Projects for Sanitary Filling of Urban Solid Wastes. The WLC technique does not always provide an adequate analysis
because it assumes average risk, not considering factors with a greater
Theme Function Control Points degree of limitations than others and presenting deficiencies when
applied in more complex areas. On the other hand, the OWA technique,
a b c d
in addition to allowing for the compensation of factors through the
Roads Decreasing 15 m 15 m 15 m 17,750 m application of weights, also controls the level of risk (Eq. (2)) of the
linear analysis and the degree of compensation (Eq. (3)), since it involves two
Hydrography Increasing linear 200 m 17,870 m 17,870 m 17,870 m sets of weights: factor importance criterion weights and order weights.
Slope Decreasing 0% 0% 0% 30%
In this way, risk levels depend on the position of the weights ordered in
linear
Urban Area Decreasing 500 m 500 m 500 m 30,560 m the ranking and the magnitude of their values, with higher values lo-
linear cated in the first positions representing lower risk and, on the other
Geological Faults Increasing linear 60 m 17,800 m 17,800 m 17,800 m hand, higher values in the last positions resulting in greater risks
Soil – Scale (0–255)
(Romano et al., 2015; Calijuri et al., 2007).
Vegetation – Scale (0–255)
Conservation – Scale (0–255) n
1 ⎡ ⎤
Units Risk = . ∑ (n − 1). Oi⎥
n − 1 ⎢ i=1 (2)
⎣ ⎦
66
U.R.V. Aires et al. Land Use Policy 70 (2018) 63–70
Fig. 4. Classification of the land use and the land cover: (a) before and (b) after the tailings dam failure.
67
U.R.V. Aires et al. Land Use Policy 70 (2018) 63–70
basin into stretches according to the degree of environmental impact where the maximum suitability value was 213. Scenarios 1 and 3 ob-
(Renova, 2017). tained suitability values of 158 and 124, respectively. The low value of
The impacts caused by the rupture of the tailings dam indicate the suitability identified in scenario 3 is related to the weights assigned to
great need for the increased surveillance and monitoring of dams; each factor, indicating that the low risk condition is not suitable for the
Morgenstern et al. (2016) note several failures in the construction analysis of the study area; this is because the best solution is not always
mechanisms of the Fundão dam as it increased in volume, resulting in associated with lower risk (Eastman, 2001).
its failure. It is estimated that there are approximately 126 tailings For the definition of the suitable areas for the disposal of the tail-
dams in Brazil that may present structural problems, evidencing the ings, the criterion employed by Calijuri et al. (2002) in the identifica-
need for effective management strategies in large-scale mining opera- tion of areas for the disposal of urban solid waste was used, in which
tions (Garcia et al., 2017). areas with suitability greater than 200 fulfilled the needs of the project.
Thus, we selected areas with a suitability value higher than 200, which
3.2. Multicriteria analysis used to identify the suitable areas for tailings are found in scenario 2.
disposal In Fig. 6, it is possible to verify the locations of the suitable areas for
the tailings disposal in the municipalities of Mariana and Barra Longa,
Based on the criteria established in Table 2, the areas suitable for with an emphasis on the larger areas.
tailings disposal were obtained considering a scale of suitability ranging The site (Area 1) identified for the disposal of the tailings dam in the
from 0 to 255, where the closer to the maximum value, the more sui- analysis is located approximately 1.0 km of the old Bento Rodrigues
table is the site of interest. district and 2.5 km from the site where the failure of the tailings dam
The best result was found using the information for scenario 2, occurred. The area selected in this municipality has a size of 0.25 km2,
68
U.R.V. Aires et al. Land Use Policy 70 (2018) 63–70
Fig. 6. Suitable areas for tailings disposal obtained in the municipalities of Mariana and Barra Longa using the analysis performed with the scenario 2 conditions.
and the maximum suitability reached at this location was 207. The conditions of the material, such as the water content, and the possibility
selected area is also located in one of the places most affected by the of carrying out compacting operations, considering the lack of areas for
tailings dam failure. residue disposal, particularly because of the difficulty in obtaining en-
The site (Area 2) selected for tailings disposal in the municipality of vironmental licenses (Nery et al., 2014). An alternative to this type of
Barra Longa is also the one that presented the largest size, at 0.037 km2. storage can be dry stacking, which is a method used to replace the
This site is located approximately 3 km from the urban area of Barra construction of dams. This methodology aims to use a filter to remove
Longa and obtained a maximum suitability of 209. the moisture from the tailings. In a study by Gomes et al. (2016), the
The removal of the remaining material in the river channels and in use of this methodology demonstrated the viability of the dehydration
the most affected areas is essential to avoid the continuous con- of iron ore tailings, in which a final moisture content of 15% was ob-
tamination of the watercourses. For comparison, the accumulation of tained. The results are applicable to dry stacking, in which the final
the remaining tailings along the river banks of up to 40 cm thick was moisture content of these materials is recommended to be between 20%
observed (Gomes et al., 2017). In the preliminary analysis of the che- and 25% (Bowker and Chambers, 2015).
mical and cytotoxic profile, highly mobilized fractions of Al, As, Ba, Fe,
Mn, Pb, and Sr were observed in the tailings samples, which had the 4. Conclusions
potential for cytotoxicity and DNA damage (Segura et al., 2016).
The definition of suitable areas for the disposal of waste can help in Based on the objectives proposed in this study and the results ob-
the current proposal for the management plan for the waste in the Rio tained, the following can be concluded:
Doce river basin. The management plan for the tailings proposed by
Renova (2017) includes the use of techniques to remove the tailings • It was possible to successfully identify potential areas for the dis-
from affected areas, such as the excavation, dredging and manual re- posal of mining tailings using multicriteria analysis.
moval of the tailings in the most sensitive areas, including the affected • The analysis of changes in land use and occupation due to the
region in the Bento Rodrigues district. In this context, the suitable areas rupture of the tailings dam provides information that can help in the
identified in the multicriteria analysis can provide essential information quantification of the impact caused and the adoption of measures to
for this project. recover the affected areas.
The appropriate disposal of the tailings that accumulated in the • The sites that presented the greatest potential for the disposal of the
affected areas was also one of the main measures for the recovery of the mining tailings were identified near the urban area of Barra Longa
areas affected by the Aznalcóllar dam rupture in Spain in 1998, with the city and in the Bento Rodrigues district in the municipality of
disposal of the tailings in suitable places that allowed the treatment of Mariana.
this waste (Olíes et al., 2012).
We recommended that the tailings disposal be conducted using the
stacking procedure, where the volume to be stored depends on the
69
U.R.V. Aires et al. Land Use Policy 70 (2018) 63–70
70