Submit manuscript...
MOJ
eISSN: 2573-2919

Ecology & Environmental Sciences

Research Article Volume 9 Issue 2

Assessment of the influence of DTM quality on dam rupture simulation processes

Marcelo Antonio Nero,1 Vinícius Tavares Parreiras de Morais,1 Marcos Antonio Timbó Elmiro,1 Ricardo Alexandrino Garcia,2 Jorge Pimentel Cintra,3 Nestor Alonso Mancipe-Muñoz4

1Departamento de Cartografia, Universidade Federal de Minas Gerais, Brasil
2Departamento de Geografia, Universidade Federal de Minas Gerais, Brasil
3Universidade de São Paulo, Museu Paulista, Brasil
4Universidad Nacional de Colombía, UNC – Facultad de Ingeniería, Colombía

Correspondence: Marcelo Antonio Nero, Universidade Federal de Minas Gerais, Brasil, Departamento de Cartografia, Av. Antônio Carlos, 6627, Belo Horizonte, Minas Gerais, Brasil, Tel 31270901 +55 (31) 3409.5000

Received: March 25, 2024 | Published: April 5, 2024

Citation: Nero MA, de Morais VTP, Elmiro MAT, et al. Assessment of the influence of DTM quality on dam rupture simulation processes. MOJ Eco Environ Sci. 2024;9(2):61-70. DOI: 10.15406/mojes.2024.09.00308

Download PDF

Abstract

Computer programs applied to disaster simulation are widely used and widespread today, taking as input various data types, from specific to the application area to DTMs (Digital Terrain Models). This need for data input and, in particular, data related to relief is very relevant in prediction procedures for forecasting catastrophes, such as the failure of tailings mining deposit dams. Therefore, it is of fundamental importance to know and quantify the quality of this input data in question, in order to effectively serve this application. To this end, in this research, several tests were carried out, using as a reference for best results those obtained with the introduction of DTM from a LIDAR (Light Detection and Ranging) flight survey, this data being used as one of the primary and input into HEC-RAS (Hydrologic Engineering Center-River Analysis System). Subsequently, these same data had their spatial resolution degraded, that is, the pixel size increased, generating models with worse quality for new inputs and obtaining new simulation results of mining dam failures. The test area of the B1 Dam, located in Brumadinho-MG, Brazil, was used as a laboratory, where one of the biggest mining disasters in the world recently occurred and, for which, there are real data from the area affected by the dam collapse. The results obtained demonstrated that the use of an DTM with a spatial resolution of at least 2.5 meters or better, with DTM class A cartographic quality compatible with the most recent Brazilian standard, would guarantee reliable results.

Keywords: dam, digital terrain model, quality, mining, map accuracy standards, rupture

Abbreviations

DTM, digital terrain model; HEC-RAS, hydrologic engineering center-river analysis system; LIDAR, light detection and ranging; MAS, map accuracy standards

Introduction

The process of simulation and prediction of mineral exploration dam failure is extensively covered by several authors, and even similar themes using HEC-RAS (Hydrologic Engineering Center - River Analysis System) are presented in several references around the world.1–12

In the context of the application of this technology and the present work, it is worth highlighting that the input of primary data, that is, topographic and/or cartographic data, is extremely important, in order to collaborate in this process, as stated and the aforementioned subject covered.11–15 Thus, with the DTM (digital terrain model) as one of the main elements, considered this as a primary data of Cartography, this research sought to carry out several simulations and tests, in order to assess which spatial resolution would be sufficient. In this sense, we sought to ensure the presentation of dam rupture scenarios, which were as close as possible to reality and, in order to guarantee the safety of the population living in the surrounding area and the effects of the rupture on the environment of the affected region. . Considering these aspects, quality control of the DTM for each simulation scenario was carried out by adapting the Brazilian standard, as discussed in detail,13–16 which is based on the control method positional quality of digital maps.17

It is also worth highlighting, in terms of the state of the art of quality control specific to DTMs, several reference works, such as the work carried out,18 where the authors analyze the topic over the last three decades;19 which addresses the perspectives and estimates of needs in terms of quality of DTMs through a study with users;20 where comparison criteria are sought for the various global DTMs, taking into account the needs of users;21 a study where the authors, based on aspects of quality control of classes derived from the DTM, analyze the quality of the slope;22 where the authors study the state of the art of DTM generation using a deep learning process and also analyzing the final quality of the generated DTM;23 where the authors used quasi-spline interpolation based on the Bernstein–Bézier coefficient to generate DTM and evaluated the final accuracy;24 apply the buffer method to verify the elevation accuracy of the DTM;25 evaluate the results generated by orbital DTMs for flood analysis and evaluate the coherence of the simulation performed.

In this article, an inference of quality control focused on altimetric data was carried out, considering the application of predictions of rupture of ore tailings dams from open mines. To this end, below, a brief theory will be presented on quality control in terms of altimetry applied to DTMs, the practical development with the definition of the study area and the details of each of the stages, the analysis of the results obtained and finally the final considerations.

Altimetric quality control theory applied to DTM

As in conventional cartography, where displacement is analyzed in the plane coordinates, precision, which is presented based on research carried out.13–17 Once the control points via sampling have been established, they must be interpolated with the input and reference DTMs. Then, the altimetric difference between the input DTMs and the reference DTM are calculated, according to equation (1).

Δ Z i = Z i Z ci MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacqqHuo arcaWGAbqcfa4aaSbaaSqaaKqzadGaamyAaaWcbeaajugibiabg2da 9iaadQfajuaGdaWgaaWcbaqcLbmacaWGPbaaleqaaKqzGeGaeyOeI0 IaamOwaKqbaoaaBaaaleaajugWaiaadogacaWGPbaaleqaaaaa@4724@   (1)

Where: ∆𝑍𝑖= elevation error at point i of the sampling set of𝑛samples; 𝑍𝑖= Altimetry of point i in the DTM under analysis; 𝑍𝑐𝑖= Altimetry of point i in the reference DTM (greater accuracy).

It is expected that the errors found in the sample set are random and have a normal distribution, so that it is not possible to determine the reason for such a difference. The probability of an error greater than three times the standard deviation is small, although possible. Therefore, if these values ​​occur, the points must be inspected and, if an inconsistency is confirmed, excluded from the analyzes as they are considered gross errors.15–17

To complete the analysis, the samples are subjected to two statistical tests, namely: the Normal Z test, which highlights the occurrence of systematic errors, and the Chi-square test (X²), which presents the precision of the sample set in relation to Brazilian regulations.17

The Z Trend Test is applied to a sample set with a large number of elements and reveals whether there are systematic errors in the set. Through equation (2).

Z H = 1 EP *μ* n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsacaWGAb qcfa4aaSbaaSqaaKqzadGaamisaaWcbeaajugibiabg2da9Kqbaoaa laaakeaajugibiaaigdaaOqaaKqzGeGaamyraiaadcfaaaGaaiOkai abeY7aTjaacQcajuaGdaGcaaGcbaqcLbsacaWGUbaaleqaaaaa@452A@   (2)

Where: Z𝐻= the calculated Z value; 𝐸𝑃= S = σ= sample standard deviation; 𝜇= mean of the sample set; 𝑛= number of samples

From the statistical table we obtain 𝑍(1−𝛼/2) based on the sample size and the significance α adopted in the analysis, a theoretical Z value is provided which, compared to the Z𝐻 value, indicates the occurrence or not of systematic errors in the samples. In the present study, a significance level of 0.1 was adopted.

This way we have the theoretical Z, according to equation 3.

Z (1α/2) = Z 0.95 = 1.6449 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadQfajuaGpaWaaSbaaSqaaKqzadGaaiika8qacaaIXaGa eyOeI0IaeqySdeMaai4laiaaikdapaGaaiykaaWcbeaajugib8qacq GH9aqpcaWGAbqcfa4damaaBaaaleaajugWa8qacaaIWaGaaiOlaiaa iMdacaaI1aaal8aabeaajugib8qacqGH9aqpcaqGGaGaaGymaiaac6 cacaaI2aGaaGinaiaaisdacaaI5aaaaa@4DFB@   (3)

Therefore, while the calculated Z value, ZH, is greater than the theoretical Z,𝑍0.95, there will be a systematic error in the altimetric component of the DTMs.

The Chi-square precision test is performed using equation (4) and reveals whether the analyzed data meets the quality standards established by the standard being applied. Therefore, the value of X H,n1 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadIfajuaGdaqhaaWcbaqcLbmacaWGibGaaiilaiaad6ga cqGHsislcaaIXaaaleaajugWaiaaikdaaaaaaa@3F78@ of the sampling point is compared to a calculated value𝑋𝑛2−1, which depends on the sample set n and the level of trust.

X H,n1 2 =(n1)* S H 2 / σ H 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadIfajuaGdaqhaaWcbaqcLbmacaWGibGaaiilaiaad6ga cqGHsislcaaIXaaaleaajugWaiaaikdaaaqcLbsacqGH9aqpcaGGOa GaamOBaiabgkHiTiaaigdacaGGPaGaaiOkaKqbaoaalyaakeaajugi biaadofajuaGdaqhaaWcbaqcLbmacaWGibaaleaajugWaiaaikdaaa aakeaajugibiabeo8aZLqbaoaaDaaaleaajugWaiaadIeaaSqaaKqz adGaaGOmaaaaaaaaaa@5376@   (4)

Where: X H,n1 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadIfajuaGdaqhaaWcbaqcLbmacaWGibGaaiilaiaad6ga cqGHsislcaaIXaaaleaajugWaiaaikdaaaaaaa@3F78@ = X² calculated for the sample;𝑛= number of samples; S H 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadofajuaGdaqhaaWcbaqcLbmacaWGibaaleaajugWaiaa ikdaaaaaaa@3C28@ = the square of the sample standard deviation; σ H 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiabeo8aZLqbaoaaDaaaleaajugWaiaadIeaaSqaaKqzadGa aGOmaaaaaaa@3D13@ = the square of the standard error of the precision class used.

The calculated value X n1,α 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadIfajuaGdaqhaaWcbaqcLbmacaWGUbGaeyOeI0IaaGym aiaacYcacqaHXoqyaSqaaKqzadGaaGOmaaaaaaa@404A@ ,  is obtained by equation (5) due to the number of samples.

X n1,α 2 =df* (1p+z* p ) 3 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadIfajuaGdaqhaaWcbaqcLbmacaWGUbGaeyOeI0IaaGym aiaacYcacqaHXoqyaSqaaKqzadGaaGOmaaaajugibiabg2da9iaads gacaWGMbGaaiOkaiaacIcacaaIXaGaeyOeI0IaamiCaiabgUcaRiaa dQhacaGGQaqcfa4aaOaaaeaajugibiaadchaaKqbagqaaKqzGeGaai ykaKqbaoaaCaaabeqaaKqzadGaaG4maaaaaaa@50C0@   (5)

Where: df=n-1 = degrees of freedom; p= 2 9 *df MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadchacqGH9aqpjuaGdaWcaaGcbaqcLbsacaaIYaaakeaa jugibiaaiMdaaaGaaiOkaiaadsgacaWGMbaaaa@3E71@ auxiliary variable; 𝑧 = −1.645 = value of the normal curve with area of ​​α = 0.1, on the right.

So, if X H,n1 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadIfajuaGdaqhaaWcbaqcLbmacaWGibGaaiilaiaad6ga cqGHsislcaaIXaaaleaajugWaiaaikdaaaaaaa@3F78@ X n1,α 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcLbsaqaaaaa aaaaWdbiaadIfajuaGdaqhaaWcbaqcLbmacaWGUbGaeyOeI0IaaGym aiaacYcacqaHXoqyaSqaaKqzadGaaGOmaaaaaaa@404A@ , then the sample meets the tested accuracy. The concepts used can be better studied.15–31

The quality parameter adopted was the Cartographic Accuracy Standard for Digital Cartographic Products CAS-DCP (PEC-PCD), which establishes 4 classes for framing geospatial data, considering the scale of the data. To do so, it is necessary that 90% of the sample values present an error below the value stipulated in the CAS-DCP (PEC-PCD) and that the altimetric variation between data under analysis and reference data present values less than or equal to the Maximum Error –ME and the Standard Error –SE(EP).32–34 Considering the global context of experiences and applications of positional quality control standards in Cartography, it is worth consulting several additional sources in other countries, such as South Africa,35–38 Australia,39 Spain,40 United States,41–43 Mexico,44,45 Portugal,46 as well as that organized by bodies that bring together several countries, such as the NATO standard (Organization of North Atlantic Treaty47).

As presented in ET-CQDG, the absolute positional accuracy for DTM, DEM and DSM data is obtained from the Maximum Error (ME) and the Standard Error SE (EP). These values applied in Table 1 establish the Cartographic Accuracy Standard for Digital Cartographic Products, CAS-DCP (PEC-PCD),32–34 according to the reference class and scale.

CAS-DCP (3)

1:1.000

1:2.000

1:5.000

1:10.000

1:25.000

1:50.000

1:100.000

1:250.000

CAS

SE
(m)

CAS
(m)

SE
(m)

CAS
(m)

SE
(m)

CAS
(m)

SE
(m)

CAS
(m)

SE (n)

CAS
(11)

SE
(m)

CAS
(m)

SE
(m)

CAS
(111)

SE
(11)

A

0,27

0,17

0,27

0,17

0,54

0,34

1,35

0,84

2,70

1,67

5,50

3,33

13,70

8,33

27,00

16,67

B

3,50

0,33

0,50

0,33

1,00

0,66

2,50

1,67

5,00

3,33

10,00

6,66

25,00

16,66

50,00

33,33

c

0,60

0,40

0,60

0,40

1,20

0,80

3,00

2,00

6,00

4,00

12,00

8,00

30,00

20,00

60,00

40,00

D

0,75

9,50

0,75

0,50

1,50

1,00

3,75

2,50

7,50

5,00

15,00

10,00

37,50

25,00

75,00

50,09

Table 1 Altimetric cartographic accuracy standard for listed points and DTM, DEM and DSM for the production of digital cartographic products. Source: Original from 32, 33, 34

Material and methods

The development of this job followed structured steps, based on the practical methodology based,27–31 which was adapted for the purposes of this work. Thus, the general methodology research was structured according to Figure 1.

Figure 1 General research development methodology. Source: Authors (2024).

Study area

The study area (Figure 2) corresponds to the valley from the headwaters of Ribeirão Ferro-Carvão to its confluence with the Paraopeba River and, from then on, following this watercourse to the stretch where the river crosses the Serra do Funil, on the edge of the municipalities of Brumadinho, Mário Campos and São Joaquim de Bicas, State of Minas Gerais. This state corresponds to one of the largest in terms of territorial extension in Brazil, being the third economy in the GDP ranking, where the most important commercial activity corresponds to the mining sector.

Figure 2 General location of the study area in relation to the city of Belo Horizonte – MG and Brazil.

In the most upstream portion of the valley is located the Córrego do Feijão Mine, operated by the mining company Vale, one of the largest in the world and leader in the sector in Brazil. Among the component structures, it is worth highlighting the B1 and B6 dams with storage capacity of 11.7x106 m³ and 4x106 m³ respectively, located at the head of Ribeirão Ferro-Carvão. Immediately below these structures, there is the ITM – Ore Treatment Installation followed by the railway stop, wagon loading yard for transporting the ore. Continuing with the description of the area, it is worth highlighting the buildings that made up the administrative part of the mine and, further down, two smaller dams B4 and B4A, for containing sediments (tailings).

Following, the course of Ribeirão Ferro-Carvão passes through the district of Córrego do Feijão, by Pousada Rural Nova Estância, followed by a railway bridge. At the end of the river there is the Parque da Cachoeira neighborhood where the Alberto Flores bridge stands out, close to the confluence with the Paraopeba River. In the section considered, the Paraopeba River is bordered by the urban center of the municipality of Brumadinho where there is a bridge on the MG040 highway. In the final stretch of the study area, the river crosses the Serra do Funil.

It is worth highlighting as a historical report that on January 25, 2019, the B1 tailings dam collapsed, after its mass underwent the liquefaction process, releasing the volume of iron ore tailings mud contained in its reservoir into the valley. downstream, described previously.48 The dam failure mode occurred due to liquefaction of the massif, the breach formation time was 6 seconds and the reservoir volume propagation time was around 5 minutes.48

During its propagation, the wave reached a speed of more than 90 km/h and a depth of more than 40 meters.48,49 As a consequence, most of the operational area of ​​the Córrego do Feijão Mine was destroyed, as were the structures and houses along the Ribeirão Ferro-Carvão valley, resulting in 259 people dead and 11 missing.50

The wave reached the administrative part of the mine in less than 2 minutes. Due to the geomorphological characteristics of the valley downstream and with the incorporation of material along the wave path, the speed decreased considerably, taking around 30 minutes to reach the region of the Nova Estância guesthouse. Finally, at around 3:50 pm the wave reached the Paraopeba River, after more than 3 hours of the collapse of Dam 1.49

Figure 3 presents a greater detail of the study area with the main elements observed, and Figure 4 presents an example of the visual appearance of part of the study area before and after the dam collapse.

Figure 3 Detail of the study area.
Source: 13.

Figure 4 Detail of part of the affected area before and after the rupture.
Source: 13.

In the data selection and processing procedures in the Geographic Information System – GIS, the ArcGIS and QGIS programs were used,51–53 including the tools and applications associated with them.

DTM degradation

Furthermore, it is worth highlighting and presenting details of the DTM images, which were generated from the original file to those degraded with the worsening of spatial resolution to carry out practical tests for introduction into HEC-RAS, generating different simulation scenarios. This can be seen in Figure 5, where there is the original image from DTM (C01, 1 meter pixel spatial resolution) to the others (C05, 5 meters; C10, 10 meters; C15, 15 meters; C20, 20 meters; C25, 25 meters), with bilinear interpolation being applied in this process in a GIS program.51–53 In these figures, you can see the difference in visualization of the railway line layout.

Figure 5 Detail of part of the affected area before and after the rupture.
Source: 13.

Sample size, distribution and selection

When locating the sampling points, locations were considered that did not undergo changes due to the rupture of Dam B1 and that were less likely to change their morphological characteristics over time. In this way, preference was given to demarcating control points in easement areas, access roads and railways, identified in the orthophotos used in this study.

In total, 247 points were established to control the vertical accuracy of the DTMs, selected as shown in Figure 6. These sample points were interpolated to obtain the altimetric value in each input DTM, previously presented in Figure 5, and in the reference DTM, coming from of the LIDAR survey.

Figure 6 Sample of control points in the study area.
Source: 13.

Statistics analysis

The synthesis of the descriptive statistics of the data, the trend and precision tests and the quality standard of the evaluated DTMs (C01, 1 meter, C05, 5 meters; C10, 10 meters; C15, 15 meters; C20, 20 meters; C25, 25 meters) based on the previously presented theory are presented in Table 2. In this same Table 2, all the classical statistics developed are presented, considering the corresponding sample of each of the DTMs with different spatial resolutions/pixel size (C01 , 1 meter, DTM_1m; C05, 5 meters, DTM_5m; C10, 10 meters, DTM_10m; C15, 15 meters, DTM_15m; C20, 20 meters, DTM_20m; C25, 25 meters, DTM_25m). In addition, several other variables are presented in this same table, namely: mean, median, mode, standard error, standard deviation, sample variance, kurtosis, mean plus 3 standard deviations, mean minus 3 standard deviations, asymmetry, interval, minimum, maximum, sum and sample size.

Sample

DTM_1m

DTM_5m

DTM_10m

DTM_15m

DTM_20m

DTM_25m

Average

0.052

0.069

0.099

0.141

-0.121

-0.153

Median

0.123

0.099

0.044

0.094

-0.047

-0.097

Fashion

-0.491

-0.503

-0.765

-1.597

1,917

0.25

Standard Error

0.027

0.028

0.032

0.042

0.048

0.059

Standard deviation

0.419

0.431

0.494

0.651

0.747

0.908

Sample variance

0.176

0.186

0.244

0.424

0.559

0.824

Kurtosis

2,264

1,412

1.025

2,660

0.873

2,030

Average + 3*Standard Deviation

1,309

1,362

1,581

2,095

2,121

2,571

Average - 3*Standard Deviation

-1.205

-1.224

-1.382

-1.813

-2.363

-2.876

Asymmetry

-0.483

-0.288

0.2

0.503

-0.233

0.228

Interval

3,518

3,324

3,724

5,181

4,543

6,800

Minimum

-2.181

-2.051

-1.953

-1.666

-2.626

-3.443

Maximum

1,337

1,273

1,771

3,515

1,917

3,357

Sum

12,542

16,579

23,781

33,835

-29.057

-36,614

Sample Size

240

240

240

240

240

240

Table 2 Results of statistical analyses.
Source: Adapted from 13.

Taking a more visual approach, Figure 7 shows the error distribution graphs for each of the DTMs that were used in the dam failure simulations in HEC-RAS processing. It can be seen that, as expected, the data does not adhere to a normal distribution. However, this is an increasingly common pattern in the engineering area and more specifically in terms of geotechnologies, considering data collection with DRONES and GNSS, in terms of positional accuracy, obtaining the closest and most adherent results. the reality.

Figure 7 Distribution of sample errors from each simulation.
Source: 13.

Next, considering the cartographic quality inferences, the concepts of the item referring to the Theory of altimetric quality control applied to the DTM, already presented previously, were applied. Thus, the results in terms of cartographic quality are presented in Table 3, as well as the final classification according to the standards of the Brazilian standard (see based on the parameters in Table 2,32–34).

Sample

DTM_1m

DTM_5m

DTM_10m

DTM_15m

DTM_20m

DTM_25m

Sample Size

240

240

240

240

240

240

Z Test (trend)

1,619

1,622

1,535

1,308

-1.123

-1.415

Z(1-α/2) or Z0.95

1,645

1,645

1,645

1,645

1,645

1,645

SE (CAS-DCP PEC-PCD) (m)

0.5

0.66

1,000

1,670

1,670

1,670

Accuracy Test (X²z,n1)

167,867

101,911

58,286

36,369

47,877

70,625

X²n-1, α

267,412

267,412

267,412

267,412

267,412

267,412

CAS-DCP

D

B

D

B

B

B

Scale

1:2,000

1:5,000

1:5,000

1:10,000

1:10,000

1:10,000

Maximum Error (m)

0.75

1

1.5

2.5

2.5

2.5

Table 3 Results in terms of cartographic quality control of the DTMs and in relation to the Brazilian standard.
Source: Adapted from 13.

Dum break inferences

Before even starting the dissertation on this item, it is worth highlighting that the main objective of this article was to report the application and proposition of an innovative methodology for quality control of DTMs used in the dam failure simulation process and their influence, in particular with respect to HEC-RAS.54–56 In terms of theory and concepts of the parameters and variables that will be reported and evaluated here, there is already a vast list of documentation and scientific articles, with emphasis on,1–14 which are used in hydrology and widely discussed and enshrined in the language of this scientific area. Therefore, there was no need to address these issues in theoretical terms, moving more towards application.

In the process of evaluating the results of dam failure parameters, the objective was to analyze the impacts in terms of the classification of each DTM of different spatial resolutions (pixel size) on the simulation results. To this end, the analysis of the flooded area forecast (performance indicator F) was verified, as well as the values ​​obtained for depth, speed and arrival time of the waste wave. Where the columns with C01, C05, C10, C20 and C25, respectively correspond to the scenarios associated with the DTMs of the different spatial resolutions of 1 m, 5 m, 10 m, 20 m and 25 m, respectively.

The performance indicator F was obtained by comparing the flooded area mapped in the post-rupture orthophoto and the flooded areas calculated in the simulated scenarios. Thus, the results shown in Table 4.13 and plotted in the graphs in Figure 8,13 considering the Bias values ​​all above 1.0, demonstrate that the scenarios tended to overestimate the flood level; as shown by the penalty rates for overestimates, F4, with values ​​around 0.6; The F² index above 0.8 indicates the critical success achieved by all scenarios in predicting flood spots.

 

CO1

CO5

C10

C15

C20

C25

Vies

1.169

1.173

1.167

1.162

1.155

1.134

F2

0.822

0.819

0.821

0.822

0.824

0.828

F3

0.804

0.802

0.802

0.802

0.801

0.799

F4

0.662

0.657

0.662

0.665

0.670

0.685

Table 4 Results of performance indicators F.
Source: 13.

Figure 8 Performance indicator graphs F.
Source: 13.

In terms of classification of affected areas (A-coherent; B-overestimated; C-underestimated), compared to reality, Table 5 is presented.

 

Col (ha)

CO5 (ha)

C10 (ha)

C15 (ha)

C20 (ha)

C25 (ha)

A

264,49

264,56

263,98

263,56

262,97

260,78

B

51,52

52,55

51,29

50,31

49,14

45,12

C

5,86

5,73

6,24

6,58

7,17

9,04

Total

321,87

322,83

321,51

320,45

319,27

314,95

Table 5 Results of performance indicators F in terms of area classification.
Source: Adapted from 13.

Analyzing Table 6, we have to:

  1. 1) the total of areas coherent with the reality of the disaster that occurred (A-coherent, marked in brown), with the highest value found in C05 being 264.56 ha (very close to C01, 264.49 ha) and the least coherent being the area of scenario C25 of 260.78 ha;
  2. in terms of overestimated area (B-overestimated, marked in blue) we have in C05 the value of 52.55 ha (being the most alarming and very close to C01 of 51.52 ha), with the least alarming being C25, pointing to the area of ​​45.12 ha (see blue marking, with a difference of 7.43 ha);
  3. considering the underestimated areas (C-underestimated, marked in red) there would be the highest value and, therefore, the greatest error in decision-making, which could lead to negligence and the occurrence of fatalities, would be the C25 with 9.04 ha; C05, with a value of 5.73 ha (very close to C01 of 5.86 ha), being presented as the most coherent and, therefore, the lowest risk in decision making), indicating a difference of 3.18 ha, that is, of 63.38% increased risk of making mistakes when making a decision;
  4. in terms of what would be real (marked in black) and, considering the spatial resolution, C01 would be the closest to the real area due to the fact that it has the best spatial resolution, having the total area considering A, B and C, of ​​321.87 ha and although C05 had a very close value (322.83 ha, difference of just 0.3%, that is, practically equal); when comparing with the theoretically less precise scenario (C25), due to its worse spatial resolution, with an area of ​​314.95 ha, there is a difference of 6.92 ha), which would represent a difference of 2.14%, the which depending on the area could be considerable.

On the other hand, in visual terms in the representation, in the form of maps of these same scenarios and, considering the scale of representation, as seen in Figure 9, the values ​​presented are very close, being practically imperceptible.

Figure 9 Performance indicator graphs F.
Source: 13.

Continuing the application's inferences, it is clear that the fluid simulated in the scenarios did not travel the entire length of the calculation mesh, reaching the stability condition just after the representative section ST-10 (fluid transport section-10, of the total out of 14 that were considered). The results of depth, speed and wave arrival time in the sections are represented in Tables 6–8, analytically and visually in the graphs in Figures 10–12.

In terms of maximum depth, see Table 6 and Figure 10.

Observing Table 6 and the graph in Figure 10, it can be seen that in relation to maximum depths, the greatest variation occurred in the initial cross section, ST-01, and in the final ones, from ST-07, already in the middle end of the wave propagation section. When observing Table 7 and the graph in Figure 11, in relation to the speeds in the final sections, ST-09 and ST-10, there is greater sensitivity to variation in the terrain scale. Regarding arrival times, which can be seen in Table 8 and the graph in Figure 12, it can be seen that in the first stretches, more specifically in the first 3 kilometers, there was no great variation, with a significant difference being found from the section ST-06. Finally, analyzing the set, maximum speed was the metric that showed the greatest variation between scenarios with an average of 4.2%, followed by maximum depth with 3.75% and finally, the arrival time of the wave with 2%.

Section

Distance from Dum (km)

Maximum depth (m)

Average

Standarddeviation

Coefficient of variation

C01

C05

C10

C15

C20

C25

ST-01                  

0,01

24,84

24,37

23,74

24,16

23,04

21,60

23,63

1,16

4,93%

ST-02

0,78

25,03      

25,00      

24,73      

24,91      

24,91      

25,42      

25,00      

0,23

0,92%

ST-03

1,27

15,28

15,42

15,07

15,61

15.

16.20

15,49

0,39

2,54%

ST-04

1,89

16,90

16,15

15,69

16,68

17,23

16.

16,48

0,55

3,36%

ST-05

3,03

19,83

19,83

19,33

20,01

19,45

19,21

19,61

0,32

1,65%

ST-06

4,13

16,88

16,76

16,51

16,42

16,76

16,76

16,68

0,18

1,07%

ST-07

5,09

14,04

13,54

13,57

12,62

13,98

12,65

13,40

0,63

4,67%

ST-08

7,07

10,70

10,73

10,93

9,75

10,00

9,85

10,33

0,52

5,01%

ST-09

9,31

5,71

6,08

5,57

6,01

6,20

5,03

5,77

0,43

7,46%

ST-10

10,80

7,85

8,07

8,24

8,66

9,01

9,03

8,48

0,50

5,89%

Table 6 Maximum depth analytical table.
Source: 13.

Section

Distance from Dum (km)

Speed (m/s)

Average

Standarddeviation

Coefficient of variation

C01

C05

C10

C15

C20

C25

ST-01                  

0,01

86,32      

85,63      

84,75      

83,00      

81.

79.

83,36

2,72

3%

ST-02

0,78

70,70

72,14

72.

73,84

71.

76.

72,58      

1,86

3%

ST-03

1.

58,04

58,05

59.

59,32

59,41      

60.

58,91

0,82

1%

ST-04

1,89

61,33

61,44

60,65

61,32

61,06

61.80      

61,27

0,39

1%

ST-05

3,03

44,38

44,89

45,32

49,47

46.

47,22

46,16

1,89

4%

ST-06

4,13

27,36

27,76

27,09

29,37

28,98

28,94

28,25

0,96

3%

ST-07

5,09

26,07

25,64

27,32

25,83

26,01

25,72

26,10

0,62

2%

ST-08

7,07

14,80

14,91

15,09

15,05

15.

15,07

14,95

0,14

1%

ST-09

9,31

10,14

11,13

11,66

11,57

13,61

13,11

11,87

1,28

11%

ST-10

10,80

3,89

4,12

4,71

4,50

4.

3,33

4.

0,51

13%

Table 7 Maximum speed analytical table.
Source: 13.

Section

Distance from Dum (km)

Wave arrival time (hh:mm)

Average

Standarddeviation

Coefficient of variation

C01

C05

C10

C15

C20

C25

ST-01                  

0,01

0:01      

0:01      

0:01      

0:01      

0:01      

0:01      

0:01:00      

0:00:00      

0.00%

ST-02

0,78

0:01

0:01

0:01

0:01

0:01

0:01

0:01:00

0:00:00

0.00%

ST-03

1,27

0:02

0:02

0:02

0:02

0:02

0:02

0:02:00

0:00:00

0,00%

ST-04

1,89

0:03

0:03

0:03

0:03

0:03

0:03

0:03:00

0:00:00

0.00%

ST-05

3,03

0:05

0:05

0:05

0:05

0:05

0:05

0:05:00

0:00:00

0,00%

ST-06

4,13

0:08

0:07

0:07

0:07

0:07

0:07

0:07:10

0:00:17

3,88%

ST-07

5,09

0:11

0:11

0:10

0:10

0:10

0:10

0:10:20

0:00:27

4.%

ST-08

7,07

0:20

0:20

0:20

0:20

0:19

0:19

0:19:40

0:00:27

2,26%

ST-09

9,31

0:39

0:39

0:38

0:37

0:37

0:37

0:37:50

0:00:50

2.%

ST-10

10,80

0:53

0:51

0:50

0:49

0:47

0:47

0:49:30

0:01:50

4.%

Table 8 Analytical table of arrival times.
Source: 13.

Figure 10 Maximum depth graph.
Source: 13.

Figure 11 Maximum speed graph.
Source: 13.

Figure 12 Arrival time graph.
Source: 13.

Discussion

In terms of compatible cartographic quality for the application of the research developed in this scientific and practical work in the field of Mining Engineering and Mineral Extraction, the DTMs tested reached the quality standards established in the CAS-DCP (PEC-PCD), compatible with scales from 1:2,000 to 1:10,000, with maximum errors varying between 0.75 and 2.5 meters, considering the statistical margin of gross error (three times the standard error).

Among the metrics analyzed, the maximum propagation speed of the waste flow in the sections was the variable that presented the greatest sensitivity with a variation of 4.21%, followed by the depth with 3.75% and finally, the arrival time of the wave. which presented a 2.07% variation.

Regarding the arrival time, when analyzing the results of the cross sections more specifically, it is clear that in the first initial section of three sections, up to section ST-05, no change in the simulation result was found in this item. corresponding to the arrival time. It was found that the most sensitive change occurred in the last section, in the Paraopeba River, where the arrival time showed a maximum difference of 1 minute and 50 seconds, in relation to the average time of the different scenarios presented. When compared with the real event, it can be said that the model failed to represent the real chronology of the event. One of the possible reasons for this is due to the fact that although HEC-RAS version 6.0 contemplates the propagation of viscous fluids, it still does not have the ability to vary the volumetric concentration of this fluid in time-space, as the wave it incorporates the disturbed material along its path.

The depth results showed an average standard deviation of 50 centimeters, with the greatest depth difference found in the ST-01 section, corresponding to 1.16 meters around the average. This difference occurred precisely in the C25 simulation scenario, where the terrain has worse spatial resolution.

Regarding the maximum speed results, it is notable that the large difference in values ​​occurred in the final section of the simulation, reaching a variation above 10%. When comparing the values ​​of the representative sections and the observed values, a certain coherence can be seen between what occurred and what was simulated.

The impacted area defined in each scenario was analyzed using the performance indicator F, which demonstrates that the results obtained are not biased and that they achieved critical success greater than 80% according to F². The overestimated areas were responsible for most of the errors, as demonstrated by indicators F3 and F4, which penalized the underestimated and overestimated flood areas, respectively, as also discussed and commented in the item called Dum break inferences. In a way, overestimated areas have a less harmful effect than an area impacted in a real rupture event and not predicted in a simulated event, since it is the simulated flood spots that will support the development of emergency action plans. evacuation and rescue.

Conclusion and recommendation

In general, it is possible to state that the scenarios presented a very similar performance in explaining the rupture of Dam 1, and that speed was the most sensitive variable to the change in terrain, without, however, presenting a very discrepant response in relation to the other metrics. analyzed. Thus, scales more detailed than 1:10,000 are capable of representing an event of the magnitude of Dam B1.

However, the results demonstrate an even greater need to deepen the scale analysis of terrain representation for the simulation of dam rupture events. In this sense, similar research applied to other dams with different characteristics and different downstream valleys can help to elucidate a compatible minimum scale for the development of rupture studies.

The results can also serve as a basis for defining guidelines for mapping valleys downstream of dams, optimizing resources without losing in terms of information quality.

Finally, it is worth highlighting that the diagnosis and prediction of the disaster studied here can serve as a good practical methodological basis, providing support for the definition and classification of risk areas with these characteristics. All of this could be associated with a multivariate analysis, in order to classify the risk areas with the highest and lowest priority for evacuation. Additionally, the brigade team and the public authorities could carry out training, awareness and simulation of the evacuation of the area, training the population residing in these areas, in order to prepare them for a real situation and avoid fatalities. In another context, to be studied, new conduct and planning for the recovery of these waste areas could be adopted, such as the drainage of the dammed product, with the use of waste material (for example in the production of bricks) . It is also worth considering the recovery of the environment in the degraded area with the reforestation and restoration of the region's native forest, which is already adopted by some mining companies in Brazil.

Acknowledgments

Acknowledgments for the Analysis and Modeling from Environmental System Post Graduation Program of Federal University of Minas Gerais State - Brazil.

Funding

None.

Conflicts of interest

There is no conflict of interest.

References

  1. Peker İB, Gülbaz S, Demir V, et al. Integration of HEC-RAS and HEC-HMS with GIS in Flood Modeling and Flood Hazard Mapping. Sustainability. 2024;16(3):1226.
  2. Doğan E, Temiz T, Sönmez O, et al. Dam failure analysis according to different methods in HEC-RAS, Gökçe Dam, Turkey. Acta Scientiarum Technology. 2024;46(1).
  3. Tessema BH, Gebremedhn AY, Getahun YS. Dam breach analysis and flood inundation mapping of Dire Dam, using HEC-HMS and HEC-RAS models. Sustainable Water Resources Management. 2024;10(2):45.
  4. Rahman M, Ali Md Shahjahan. Drivers of tidal flow variability in the Pussur fluvial estuary: A numerical study by HEC-RAS. Heliyon, 2024.
  5. Ennouini W, Fenocchi A, Petaccia G, et al. A complete methodology to assess hydraulic risk in small ungauged catchments based on HEC-RAS 2D Rain-On-Grid simulations. Natural Hazards. 2024:1–29.
  6. Kannapiran UM, Bhaskar AS. Flood inundation mapping of upstream region in the Adyar River basin: Integrating hydrologic engineering center's river analysis system (HEC-RAS) approach with groundwater considerations. Groundwater for Sustainable Development. 2024;24:101085.
  7. Nathanael E, Sejati W. Effective government management of flood discharge in drainage channels using HEC-RAS 6.3. 1 application. APTISI Transactions on Management. 2023;7(3):210–220.
  8. Sathya A, Thampi SG, Chithra NR. Development of a framework for sand auditing of the Chaliyar River basin, Kerala, India using HEC-HMS and HEC-RAS model coupling. International Journal of River Basin Management. 2023;21(1):67–80.
  9. Namara WG, Damisse TA, Tufa FG. Application of HEC-RAS and HEC-GeoRAS model for flood inundation mapping, the case of Awash bello flood plain, upper Awash River Basin, oromiya regional state, Ethiopia. Modeling Earth Systems and Environment. 2022;8(2):1449–1460.
  10. Mendes TA, Sousa MBD, Pereira SADS, et al. Use of the HEC-RAS model based on LiDAR information to assess urban flooding. Sanitary and Environmental Engineering. 2022;27:141–157.
  11. Morais VTP, Nero MA, Elmiro MAT, et al. Systematic review on the use of digital terrain models in dam rupture simulations. Int J Hydro. 2022;6(5):197–200.
  12. Peixoto JS, de Moraes MAE, Garcia K, et al. Performance analysis of the lisflood hydrological model in a flood event in the Madeira river basin. Int J Hydro. 2024;8(2):38‒43.
  13. Morais VTP de. Use of digital terrain models with different spatial resolutions from Lidar data in the simulation of the Dam rupture in Brumadinho-Minas Gerais. 2021.
  14. Lima RP, Elmiro MAT, Nero MA, et al. Assessment of digital terrain models in dam break simulation studies. Bulletin of Geodetic Sciences. 2021;27:e2021005.
  15. Ferreira FR. Proposal for evaluating the altimetric accuracy of LIDAR data (Doctoral dissertation). 2019.
  16. Ferreira FR, Cintra JP. Assessment of the altimetric accuracy of Light Detection and Ranging (LIDAR) data. Multidisciplinary Scientific Journal Knowledge Center. 2021;10:14–41.
  17. Cintra JP, Nero MA. New method for positional cartographic quality control in digital mapping. Journal of surveying engineering. 2015;141(3):04015001.
  18. Mesa-Mingorance JL, Ariza-López FJ. Accuracy assessment of digital elevation models (DEMs): A critical review of practices of the past three decades. Remote Sensing. 2020;12(16):2630.
  19. Ariza-López FJ, Mora EGC, Mingorance JLM, et al. DEMs: An approach to users and uses from the quality perspective. International Journal of Spatial Data Infrastructures Research. 2018;13:131–171.
  20. López-Vázquez C, Ariza-López FJ. Global digital elevation model comparison criteria: an evident need to consider their application. ISPRS International Journal of Geo-Information. 2023;12(8):337.
  21. Alba-Fernández MV, Ariza-López FJ, Jiménez-Gamero MD. A new approach to the quality control of slope and aspect classes derived from digital elevation models. Remote Sensing. 2021;13(11):2069.
  22. Ruiz-Lendínez JJ, Ariza-López FJ, Reinoso-Gordo JF, et al. Deep learning methods applied to digital elevation models: state of the art. Geocarto International. 2023;38(1):2252389.
  23. Ariza‐López FJ, Barrera D, Eddargani S, et al. Spline quasi‐interpolation in the Bernstein basis and its application to digital elevation models. Mathematical Methods in the Applied Sciences. 2023;46(2):1687–1698.
  24. Ariza-López FJ, Reinoso-Gordo JF. Edge methods for evaluating elevation accuracy in digital terrain elevation models. Cartographic Magazine. 2021;(103):33–45.
  25. Rau MI, Julzarika A, Yoshikawa N, et al. Application of topographic elevation data generated by remote sensing approaches to flood inundation analysis model. Paddy and Water Environment. 2024:1–15.
  26. Nero MA, Cintra JP, Ferreira GDF, et al. A computational tool to evaluate the sample size in map positional accuracy. Bulletin of Geodetic Sciences. 2017;23:445–460.
  27. Nero MA, Pereira TÁJ, de Paulo Ramos M. Positional quality control of orthoimages from satellite images with good spatial resolution according to PEC-PCD. Revista Contemporânea. 2022;2(6):1235–1251.
  28. Nero MA. Proposals for quality control of cartographic bases with emphasis on the positional component (Doctoral dissertation, University of São Paulo). 2005.
  29. Nogueira Júnior JB. Quality control of cartographic products: a methodological proposal (Doctoral dissertation, Universidade Estadual Paulista). 2003.
  30. Santos ADPD. Cartographic quality control: methodologies for evaluating positional accuracy in spatial data. (Doctoral dissertation, Federal University of Viçosa). 2015.
  31. Ariza-López FJ, García-Balboa J, Rodríguez-Avi J, et al. Guide for the positional accuracy assessment of geospatial data. Pan American Institute of Geography and History, Occasional Publication. 2021;563.
  32. DSG-geographic service board. Technical specification for geospatial data quality control (ET-CQDG). 2016.
  33. DSG-geographic service directorate. Technical specification standard for acquisition of vector geospatial data for the defense of the Land Force (ET-ADGV Defense F Ter), 2nd Edn. 2016.
  34. BRAZIL. Decree No. 89,817 of June 20, 1984. Establishes the regulatory instructions for the technical standards of National Cartography. Brazil. 1984.
  35. DRDLR. Department Rural Development and Land Reform. Standard for the 1:10 000 Orthophoto Map Series. Republic of South Africa. 2018.
  36. DRDLR - Department Rural Development and Land Reform. Standard for the 1:50 000 Orthophoto Map Series. Republic of South Africa. 2023.
  37. DRDLR - Department Rural Development and Land Reform. Standard for the 1:250 000 Orthophoto Map Series. Republic of South Africa. 2017.
  38. DRDLR - Department Rural Development and Land Reform. Standard for the 1:500 000 Orthophoto Map Series. Republic of South Africa. 2017.
  39. ICSM - Intergovernmental Committee on Survey and Mapping. Australian map and spatial horizontal data accuracy standard. Australia. 2009.
  40. AENOR - Spanish Standardization and Certification Association. UNE 148002:2016: Methodology for evaluating the positional accuracy of geographic information. Madrid. 2016. 46 p.
  41. Federal Geographic Data Committee. Geospatial positioning accuracy standards part 3: National standard for spatial data accuracy. US Geological Survey. 1998.
  42. Abdullah Q. The ASPRS positional accuracy standards. Edition 2: The Geospatial Mapping Industry Guide to Best Practices. Photogrammetric Engineering & Remote Sensing. 2023;89(10):581–588.
  43. ASPRS - American Society for Photogrammetry and Remote Sensing. ASPRS positional accuracy standards for digital geospatial data. Photogramm. Eng Remote Sens. 2015;81(3):A1–A26.
  44. INEGI - National Institute of Statistics and Geography. Technical standard for positional accuracy standards. Official Gazette of the Federation on December 23, 2010. 2011.
  45. INEGI - National Institute of Statistics and Geography. Spatial data calidation model: Technical annex 2 - evaluation of spatial data calibration. 2009.
  46. ​​DGT - General Directorate of the Territory. Standards and technical specifications for vector and image topographic cartography. Portugal. 2020.
  47. NSA - Nato Standardization Agency. Evaluation of land maps, aeronautical charts and digital topographic data. STANAG (Standardization Agreement) No. 2215. 7th Edn. (Text not published), NATO, Brussels. 2010.
  48. Robertson PK, Melo L, Williams DJ, et al. Report of the expert panel on the technical causes of the failure of Feijão Dam I. Commissioned by Vale. 2019. 81 p.
  49. Ragazzi L, Rocha M. Brumadinho: the engineering of a crime. Letramento Publisher. 2021.
  50. CBMMG. Tragedy in Brumadinho completes two years and is marked by tributes to the victims. 2021.
  51. ESRI. ArcGIS sesktop. Cut fill, 2021.
  52. ESRI. ArcGIS desktop. Image classification using the ArcGIS Spatial Analyst extension, 2021.
  53. QGIS. org, QGIS [software], Version 3.18.1, QGIS. org. 2020.
  54. USACE- US Army Corps of Engineers. HEC-RAS river analysis system. 2D modeling user's manual. Version 6.0. Davis, CA: US Army Corps of Engineers. Institute for Water Resources, Hydrologic Engineering Center. 2023.
  55. USACE - US Army Corps of Engineers. HEC-RAS river analysis system. 2D modeling user's manual. Version 6.0. Davis, CA: US Army Corps of Engineers. Institute for Water Resources, Hydrologic Engineering Center. 2024.
  56. Brunner GW. HEC-RAS river analysis system 2D modeling user's manual version 6.4.1. US Army Corps of Engineers—Hydrologic Engineering Center. 2023:1–171.
Creative Commons Attribution License

©2024 Nero, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.