Next Article in Journal
Exploring Effective Built Environment Factors for Evaluating Pedestrian Volume in High-Density Areas: A New Finding for the Central Business District in Melbourne, Australia
Next Article in Special Issue
Land Cover Mapping from Colorized CORONA Archived Greyscale Satellite Data and Feature Extraction Classification
Previous Article in Journal / Special Issue
Tracking Deformation Processes at the Legnica Glogow Copper District (Poland) by Satellite InSAR—I: Room and Pillar Mine District
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Land
Volume 10
Issue 6
10.3390/land10060654
land-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tracking Deformation Processes at the Legnica Glogow Copper District (Poland) by Satellite InSAR—II: Żelazny Most Tailings Dam

by
Paolo Mazzanti
1,2,
Benedetta Antonielli
1,*,
Alessandra Sciortino
3,
Stefano Scancella
2 and
Francesca Bozzano
1,2
1
Department of Earth Sciences & CERI Research Center on Geological Risks, Sapienza University of Rome, P. le A. Moro 5, 00185 Rome, Italy
2
NHAZCA S.r.l., Via Vittorio Bachelet, 12, 00185 Rome, Italy
3
Department of Civil, Construction-Architectural and Environmental Engineering (DICEAA), University of L’Aquila, Via Giovanni Gronchi, 18, 67100 L’Aquila, Italy
*
Author to whom correspondence should be addressed.
Land 2021, 10(6), 654; https://doi.org/10.3390/land10060654
Submission received: 18 May 2021 / Revised: 3 June 2021 / Accepted: 8 June 2021 / Published: 18 June 2021
(This article belongs to the Special Issue Geomatics for Resource Monitoring and Management)
Browse Figures
Versions Notes

Abstract

:
The failures of tailings dams have a major negative impact on the economy, surrounding properties, and people’s lives, and therefore the monitoring of these facilities is crucial to mitigate the risk of failure, but this can be challenging due to their size and inaccessibility. In this work, the deformation processes at Żelazny Most tailings dam (Poland) were analyzed using satellite Ad-vanced Differential SAR Interferometry (A-DInSAR) from October 2014 to April 2019, showing that the dam is affected by both settlements (with a maximum rate of 30 mm/yr), and horizontal sliding in radial direction with respect to the ponds. The load of the tailings is pushing the dam forward along the glacio-tectonic shear planes located at depth, in the Pliocene clays, causing horizontal displacements at a rate up to 30 mm/yr, which could lead to a passive failure of the dam. The measured displacements have been compared with the ones observed by in situ data from the 90s to 2013, available in the literature. The outcomes indicate that intense localized deformations occur in the eastern and northern sectors of the dam, while the western sector is deforming evenly. Moreover, although the horizontal deformation had a slowdown from 2010 until 2013, it continued in 2014 to 2019 with recovered intensity. The upper and the recent embankments are affected by major settlements, possibly due to a lower consolidation degree of the most recent tailings and a larger thickness of compressible materials.

1. Introduction

Tailings dams are earth dams without an impermeable core, designed for permanent containment of mineral waste. One of the most used construction methods is the upstream method, by which the dam banks grow during the tailing’s storage activity at a rate of about 1–1.5 m/year. This method is the most cost-effective [1,2], as it consists in the construction of a foundation base very resistant to load and of the so-called starter dam, then the tailings are discharged into the basin until it is saturated and the placing of the new material occurs within the existing impoundment, that is, the new embankments are founded on tailings [3,4]. These latter ones consist of solid or liquid waste represented by a mixture of fine particles, in some cases toxic or radioactive, deriving from mechanical separation processes of the mineral of economic interest from the gangue.
As its height and therefore the load increases, this kind of dam is more susceptible to geotechnical failure, or it can evolve in conditions of instability, as the new embankments are founded on tailings, which are in most cases highly erodible due to their fine particle size and can contain unwanted minerals and process chemicals [5].
In the past 20 years, the number of tailings dam failures has doubled [6] with dramatic consequences due to the toxicity of the stored tailings and the large run-out distances (up to hundreds of km), causing losses and major damage to health and environment. An upward trend in high-consequence failure events since 1990 has been also highlighted by Bowker and Chambers [7].
Azam and Li [8] summarized the major tailings dam failures from 1910 to 2010, reporting that the annual failure frequency for major tailings dams is about 0.12%, which is more than two orders of magnitude greater than the annual failure frequency for large water retention dams [9]. Major failures occur with a frequency that ranges from two to five times per year, whereas minor failures occur around 35 times per year [10]. The most common result of failure is bursting of the dam which surrounds the facility and release of a huge mass of liquefied tailings mixed with water in the form of a flood wave.
Stefaniak and Wróżyńska [10] and Bowker and Chambers [7] described some examples of the over 100 documented catastrophic failures of tailings dams in the world over the last 50 and 100 years, respectively, and Rico et al. [4] conducted a study on 147 cases of tailings dam failures worldwide, 26 of which had occurred in Europe. The most important causes of failure described on literature are the following: issues related to foundations, slope instability, overfilled tailings pond, mine subsidence or deformations of the dam body, unusual heavy rainfall, water pressure excess in tailings, piping/seepage, seismic liquefaction caused by earthquakes, structural issues, management, or operational errors. In many cases, no single precise cause for the failures has been identified, and various overlapping factors could combine [10]. Lyu et al. [2] also summarized the influence of the failure height, the building type, the geographical location, and time distribution of the tailings dam failures, and interestingly highlighted that upstream dams have a high probability (58%) of damage, and that the majority (85%) of tailings dam failures have occurred in dams of less than 45 m high.
On 5 November 2015, Samarco’s tailings dam breached, destroying two villages, killing 19 people, and seriously polluting 650 km of the Rio Doce River in the Brazilian state of Minas Gerais. In the same state, and with the same mine owner, in January 2019, the Corrego do Feijão dam at Brumadinho collapsed, killing about 248 people. The 86 m height tailings dam, constructed with the upstream method, catastrophically failed, releasing 12 million m3 of mineral waste. Catastrophic failures at Samarco and Brumadinho in 2015 and 2019, respectively, brought mine tailings disasters onto the global level in unprecedented ways, and several studies have developed these themes following these events (e.g., [6,11,12,13,14,15,16]).
Another very serious failure event took place in 1985, at two tailings dams nearby the town of Trento in Northern Italy. The collapse triggered a disastrous mudflow [17] that destroyed the villages of Stava and Tesero, killing 268 people. In this case, human factors contributed most to the failure; the dam’s superstructure which did not meet high safety standards [18] and the failure of the two dam bodies can be explained through excess amount of water in the tailings pond and liquefaction by infiltration phenomenon [19,20].
The assessment of tailings dam stability involves the study of the hydrogeological settings, tailing material characterization, deposition methods and of the possible trigger factors leading to dam failure (e.g., earthquakes, extreme rainfall, man-made factors [21]). However, the monitoring of deformations of the dam body is also of primary importance, even if it is challenging due to the size and inaccessibility of these facilities, which are among the world’s largest engineered structures.
In-situ geotechnical instruments, levelling, total stations, or GNSS and visual inspections represent the conventional labor-intensive and time-consuming monitoring systems generally used at tailings dams. However, traditional techniques are conducted on a point-by-point basis with limited spatial coverage [22]. Satellite remote sensing technologies such as Synthetic Aperture Radar (SAR) Interferometry (InSAR), can overcome this limitation, allowing for the analysis of wide areas, and ability to obtain information from the past, using data archives. This consolidation technique can be an effective tool for the engineers in charge of structural health monitoring of dams (e.g., [23,24,25,26,27]) and tailings dams [14,15,16,22,28]. Specifically, Advanced Differential SAR Interferometry (A-DInSAR) can provide information about the ongoing and past deformations with millimeter accuracy and can be used for the prediction of catastrophic events (e.g., [12,15]), by the detection of precursory deformation at the dam structure before the failures occur.
This paper focuses on the deformation phenomena that affect the tailings dam of Żelazny Most (Figure 1) using satellite A-DInSAR. The Żelazny Most Tailings Storage Facility (TSF) is the largest sump reservoir in Europe built using the upstream method (Figure 1d, [1,29]) and is one of the largest in the world. It stores the waste derived from copper extraction activities of the adjacent Legnica Glogow Copper District (LGCD, Poland) the most important copper mines of Lower Silesia, managed by the KGHM Polska Miedź S.A. For more details on the mining basin and on the superficial deformation phenomena that affect it, see Antonielli et al., submitted to this journal. More than 99% of the original ore extracted in the copper mines of Lubin, Rudna, and Polkowice-Sieroszowice (Figure 1c) represents waste: everyday an average of 80,000 tons of waste derived from the mechanical separation processes are produced [29] and then discharged in the basin, by means of flotation.

2. Study Area

The morphological depression in the southeastern portion of the Dalkowskie Hills indicated the site as suitable for the location of this tailings dam. The construction of the Żelazny Most TSF started in 1974, and the discharge of tailings started in 1977. The ring dams surrounding the disposal area are 14.3 km long, and the area of the basin is 16,380,000 m2, where more than 600 million m3 of tailings have been deposited [30]. The dam’s height has varied over time from 40 to 71.5 m. In this type of construction, the coarse fraction of the tailings (i.e., sand), that is separated from the fine fraction (i.e., silt and clay), is used to build the shell in a series of lifts [29]. The Żelazny Most TSF is formed of a selected sandy material, deposited on the so-called beaches of the storage facility. Only a part of this material is deposited, and most of it is used to build the embankment [10].
The construction technique causes the stability of the dam to vary over time as an embankment is added to the dam, posing a relevant threat for the surrounding areas. TSF is adjacent to urbanized and populated areas: Tarnówek village is 600 m to the west of the site boundary, Rudna village is 1 km to the east, and Grodowiec and Rynarcice villages are 1.5 km to the north and to the south, respectively, of the dam boundary. The structure is continuously monitored and adapted to avoid rupture. KGHM Polska Miedź has been preparing future tailings storage methodology to ensure deposition until 2042 (1 billion m3; [30]). The development of the Żelazny Most TSF has been planned by constructing the Southern Extension (609 ha) adjoining the main dam basin to the south.
The study site is in the Lower Silesia region (south-western part of Poland) which is a part of the Fore-Sudetic Monocline, where metamorphic rocks and volcanogenic-sedimentary beds outcrops. In particular, the Fore-Sudetic Monocline is composed of a thick epi-Hercynian sequence of sediments and rocks of Carboniferous, Rotliegende, and Zechstein age [4]. Copper ore of the LGCD mine is located at the base of Zechstein formations, and the thickness of ore layers varies from 0.1 to 15 m [4].
Tectonic recognition of faults of a large stratigraphic throw in the area of the KGHM copper ore mine reveals that the main faults strike NW–SE [31,32]. The dominant fault system in the mining area coincides with the fault zone of Rudna and Lubin, which represent two parallel structures striking NW-SE. Rudna fault is considered the most representative among the faults which pose a seismic hazard and is located about 5 km south from Żelazny Most tailings dam.
The Żelazny Most tailings dam was built on a Quaternary and Tertiary clay substrate [1,33]. The storage facility is located in a morphological depression in the southeastern portion of the moraine hills of Dałków [30,34], surrounded by frontal moraine-dammed hills of the Riss glaciation. More detail on Quaternary stratigraphy and paleogeography of the study area can be found in [35,36,37,38]
In particular, the dam is located on clayey and sandy lacustrine deposits Pleistocene aged, which rest on a thick layer of lacustrine clays with medium-high plasticity, interbedded with thin peaty or sandy layers.
During the Pleistocene, Central Europe was affected by at least three glaciations with relative progradation of glacial masses towards the south of Poland, where the study area is located [1,29,39]. The thickness of the ice masses was at least 1000 m and the weight caused deformation phenomena (the so-called glacio-tectonics) that could have affected the subsoil up to hundred meters in depth [1,29]. The thrusting (horizontal or inclined) caused shear surfaces in the Pliocene clays (Figure 2), which have a low strength, referred to as the residual strength [30]. The shear stress caused by the passage of glaciers during the cooling events generated several sub-planar shear surfaces that formed in the high plasticity clays. The orientation of the major shear surfaces within Pliocene clays is expected to follow the general direction of the Quaternary ice movement, i.e., NW-SE direction to the southeast.
The Quaternary deposits consist of sands and gravels, deposited as outwash from the Pleistocene icesheets. The glaci-tectonics thrusting pushed these large lenses of gravels over the Neogene clays. The resulting complex geometry of the Quaternary gravels outlines the drainage boundaries to the masses of Neogene clay and controls the drainage of the waters below the dam [30].

3. Methods

InSAR (Interferometric Synthetic Aperture Radar [40,41] can be considered one of the most powerful satellite remote sensing techniques used for the measurement of displacements of the Earth’s surface. It is based on the comparison between RADAR images acquired in different times over the same area. The Advanced Differential SAR Interferometry (A-DInSAR) approach consists of the combination of many SAR images, allowing for the retrieval of the temporal evolution of displacements of objects on the ground during the investigated time period [42,43]. The use of this Multi-Image technique allows for the identification of the deformation process, past or in progress, and to estimate its evolution over time and space, so it can be profitably exploited to investigate the spatial and temporal evolution of ground deformation processes with millimetric accuracy.
In this study, 411 Sentinel-1 SAR scenes, in ascending and descending orbital geometries covering the study area, have been analyzed using a Multi-Image A-DInSAR technique. Sentinel-1 (European Space Agency, ESA) is a two-satellite constellation carrying onboard a C-band SAR sensor (5.40 GHz frequency and about 5.55 cm wavelength), operating since the second half of 2014. The acquired stacks cover a period of about 4 years and 7 months (from October 2014 to April 2019). We used Interferometric Wide beam mode in Single Polarization (IW) SAR images, characterized by a spatial resolution of approximately 5 × 20 m, a wavelength of 5.6 cm (C-band), and 250 km wide images (Table 1).
A-DInSAR analysis were performed to acquire information regarding past displacements that affected the ground surface over the entire investigated area since late 2014. The multi-temporal interferometric analyses have been carried out using Persistent Scatterers Interferometry (PSI) technique [41,42,43,44,45] and with proprietary procedures implemented in SARPROZ software [46]. The accuracy of the PSI results depends on many parameters, such as the number of SAR images, the spatial distribution of PSs, climatic and atmospheric conditions, distance from the reference point, and the quality of PSs within the area of interest. Specifically, to reduce the decorrelation effects, all images were related to a single master image, which was selected by considering the perpendicular and temporal baselines. After all images were co-registered to the master image, the maps of the radar parameter were generated: the reflectivity map (a multi-temporal amplitude value for each pixel) and the amplitude stability index map (the coefficient of variation of the amplitude). These maps were used as quality estimators for the selection of persistent scatterer candidates (PSCs) in the PSI workflow. Then, the standard “star” graph (Figure 3; [43]) has been used to connect all slave images to a single master image and thus to generate interferograms. This graph represents all the SAR images according to their temporal baseline (X axis) and their normal baseline (Y axis), both for ascending and descending orbital geometries, and allows to generate a network by connecting all the images (slave images) to one (master image; Figure 3).
A network of PSCs was created to estimate preliminary height and displacement parameters. This step is needed to detect and remove the atmospheric phase screen (APS), starting from the residual phase components. After APS estimation and removal, the final estimation of the velocity and residual height for the Persistent Scatterers (PS), was performed on a much larger set of points, as the final step of the PSI procedure. For each PS the velocity values along the Line of Sight (LOS), the displacement time series, and the height, were computed (relative to a reference point identified in a stable region). The resulting PSs have been selected by applying a high temporal coherence threshold, thus selecting only pixels characterized by a temporal coherence >0.6, in order to obtain reliable time series. It is worth mentioning that displacements are measured along the satellite LOS, the path between the satellite sensor and the target (e.g., look angle in Table 1). The results of the procedure are represented by velocity of displacement maps (mm/yr) in the LOS direction for both ascending and descending datasets.
The Up-Down and East-West velocities of displacement have been retrieved by combining the ascending and descending orbital geometries [47,48]. An application for vector decomposition owned by NHAZCA S.r.l. was used for this purpose. After a discretization of the LOS maps in a grid of hexagonal cells, the vertical and horizontal velocity rates can be achieved by the combination of at least one PS from the ascending and one PS from the descending. Through the vector decomposition of the ascending and descending LOS vectors, the velocities are projected along the vertical and horizontal directions (up-down and east-west). The horizontal component of the motion along the north-south direction cannot be properly measured, because the SAR satellites fly in the near-polar orbit, approximately parallel to the north-south direction and perpendicular to the LOS [48,49]. Thus, in the non-polar regions, the N–S surface displacements contribute very little to the A-DInSAR LOS measurements.

4. Results

The A-DInSAR velocity maps of the Żelazny Most TSF, for the period October 2014–April 2019, are shown in Figure 4. The panels a and b show the LOS velocity, in ascending and descending orbital geometry respectively, while c and d the east-west and up-down velocity maps, respectively. For the LOS velocity maps, the used convention for the color bar is as follows: green PSs represent the velocity range between −2.5 and 2.5 mm/yr, which has been considered no deformation; PSs moving away from the satellite are in yellow, orange, and red, while PSs moving toward the satellite are from cyan to blue. For the east-west and up-down velocity maps, the green PSs represent no displacement; yellow to red PSs are moving westwards or downward, respectively, while cyan to blue PSs are moving eastward or upwards, respectively.
As regards the horizontal displacements (east-west direction in Figure 4c), it is worth noting that, along the eastern embankment of the dam, PSs are moving eastwards with velocity up to 30 mm/yr, while along the western embankment, PSs are moving westwards, with similar maximum rate of displacement (up to −30 mm/yr). Therefore, in these sectors of the dam, the embankment is deforming outwards, i.e., the velocity vectors are radial with respect to the center of the basin. Along the northern sector of the dam, instead, the east-west deformation rates are very small, with values ranging between 5 and −5 mm/yr. The horizontal rates in the north-south direction, as mentioned in the previous paragraph, cannot be properly measured by InSAR [49,50].
As regards the vertical displacements (Figure 4d), PSs moving downwards (subsidence) are in yellow to red colors, while there are no areas with PSs moving upwards (uplift). The velocity values measured along the vertical direction show that several parts of the dam ring are affected by subsidence or settlement (hereinafter considered as synonyms). The southern part of the dam shows lower subsidence (from −5 to −15 mm/yr), while in all other sectors of the dam, subsidence is more intense, especially in the innermost and higher rings. In the eastern sector, the subsidence rate is up to −30 mm/yr in the innermost part of the tailings dam ring, with peak values of −50 mm/yr in the central part of the eastern sector. Subsidence decreases up to values of −5 mm/yr towards the base of the dam structure, i.e., in the so-called starter dam. In the western sector of the tailings dam a similar pattern of deformation can be observed: in the innermost rings the subsidence is up to −30 mm/yr, while in the starter dam the settlements are in the order of −6 mm/yr. In the northern sector of the dam, the maximum velocity is about −35 mm/yr in the inner portion of the dam and decreases in the starter dam.
In order to analyze the trend of horizontal and vertical displacement over time, the time series of displacement has been computed for each measurement point. Some interesting examples are shown in Figure 5: for each sector of the dam, four time series selected at different highs of the dam body are represented (location of the selected points is also given in the map panels). For each dam sector, Points P1 are located near the beach and the ponds, in a higher position; points P2 and P3 are selected in the middle highs, and P4 are located on the starter dam. The time series of vertical displacements (Figure 5a) highlights that the points P1 show steep trends, as they are affected by intense deformation. The deformation pattern increases quite steadily over time (from October 2014 to April 2019), with a linear trend. The displacements become more and more moderate, going down in altitude at points P2, P3, and P4 (less steep time series).
As regards the time series of east-west displacements (Figure 5b), the eastward displacements are shown in blue, and westward displacements are in red. Also in this case, the trend of the time series shows larger cumulative displacements values for points located near the beach and the ponds (P1ew and P2ew) compared to the points at the foot of the dam (P3ew and P4ew). However, such difference in the horizontal deformation patterns is very weak in the western sector of the dam, where the displacement towards the west (or towards the outside), seems to take place as an en bloc movement.
The displacement pattern in space has been investigated through two cross sections parallel to the dam body, along the eastern and western embankments (Sections 1 and 2, located in Figure 6a,b). The cross sections run at a middling altitude of the dam body. Each curve of the section graphs (Figure 6c–f) represents the A-DInSAR displacement profile of the dam at that height, for each analyzed year. In particular, the curves of the graphs show the cumulative deformation that affected those sectors of the dam up to September of each year, and therefore the used time span is from October 2014 to September 2018. Greater distances between the graph curves represent higher velocities of displacement.
Section 1, which concerns the eastern dam, shows that the eastern perimeter of the dam does not deform uniformly, but the central zone of the section shows more intense displacement, both in the horizontal (Figure 6c) and vertical directions of movement (Figure 6d), with maximum cumulative displacement of 80 mm eastwards and 100 mm downward, in about four years. As discussed above, higher deformation values are expected at higher altitudes of the dam, and lower values are expected at the base of the dam.
Section 2, along the western dam, instead shows a uniform pattern of deformation both in the horizontal (Figure 6e) and vertical directions (Figure 6f), with maximum horizontal cumulative displacements similar with respect to the eastern sector (about 80 mm), but in the opposite direction (i.e., westwards), and the cumulative subsidence is around 60 mm in about four years.

5. Discussion

The Żelazny Most dam structure is affected by a complex deformation pattern that involves both the vertical and horizontal directions. As regards the vertical displacements, the deformation corresponds to subsidence (or settlements) occurring on the dam ring with different intensity depending on the position and altitude along the dam body: higher in the upper part and lower at the base of the dam ring. Such differential movement also concerns horizontal displacements for the eastern sector. This pattern of deformation is represented on a schematic cross section (cross-section XVIE of Jamiolkowski [1] located in Figure 6a,b) in Figure 7, where both settlements and horizontal displacement are represented for two time intervals: from 2008 to 2013, with values deduced from Jamiolkowski [1] and from 2014 to April 2019, with data measured by A-DInSAR technique. The original values of cumulated displacements provided by Jamiolkowski [1] covered different time frames for each benchmark represented on Figure 7 a,b, as each benchmark was monitored since different dates: for example, the benchmark BM-208 located at the starter dam (location in Figure 6a,b), was monitored from 1981 until 2013, while BM-208.7, on top of the dam, was added later and was monitored from 2007 until 2013. To make these displacement values comparable to each other and roughly comparable also with those retrieved by the A-DInSAR four-year analysis, the average annual velocity for each benchmark has been calculated, and the related displacement value has been obtained for a time interval of four years.
Figure 7c and d show the superficial deformation achieved by A-DInSAR technique from October 2014 to April 2019 at the points P1, P2, P3, and P4, located in Figure 5.
In both time intervals 2008 to January 2013, and October 2014–April 2019, the stronger settlement affects the highest and the most recent part of the storage facility, although the intensities of the deformation are different, and are more moderate in the second time interval, probably due to stabilization measures implemented and slowdown in dam construction operations.
In general, it is possible to conclude that on the starter dam the deformation is moderate, while it increases towards the higher, inner, and therefore most recent parts of the dam. Such spatial distribution of subsidence rates could be linked to the age of construction and therefore to a different consolidation degree of the tailings, but it could also be linked to internal zoning of the dam materials and to differences in their thicknesses. Stefanek et al. [30] described the deposition system adopted for the tailings dam construction. These authors state that the coarse tailings remain on the beach, while fine-grained (dusty) tailings flow into the pond, forming an impermeable layer underneath the water table, useful to strongly reduce seepage beneath the pond. A scheme of tailings zoning in the site’s embankment is given in Figure 8, as well as the list of the selected geotechnical parameters of the soil for each zone of the dam structure [30]. The range of subsidence A-DInSAR velocities, corresponding to each part of the dam section (eastern and the western sectors of the dam which show a similar pattern of vertical deformation), are also reported above the profile in Figure 8.
According to Stefanek et al. [30], the material type that forms the starter dam is sand mechanically condensed (2 and 2’ in Figure 8). In this basal part of the dam ring, the subsidence rate can range from 5 to 10 mm/yr. Zone I consists of sandy tailings, not mechanically condensed (3 and 3’ in Figure 8) and Zone II, which is the most inner part of the dam structure (4 and 4’ in Figure 8) is made of uncondensed dusty tailings. The parts of the embankment affected by more intense subsidence (20–30 mm/yr) are founded on greater thicknesses of sandy tailings not mechanically condensed or on uncondensed dusty tailings (note the thickness of the Zone II in Figure 8). Therefore, the more pronounced subsidence rates in the highest sectors of the dam could be linked both to the most recent age of construction and to the larger thickness of compressible materials in the inner part of the dam body.
As regards the horizontal deformation rates, the eastern sector of the dam moves towards the east and the western one moves towards the west, with vector intensity of more than 30 mm/yr in both cases. Northern and Southern sectors of the dam do not show evident horizontal deformations in the east-west direction. It could be possible that these sectors undergo horizontal deformations in the north-south direction along which the classical AD-InSAR technique cannot provide measurements [48,49].
By the analysis of the horizontal displacements that occur along the eastern and western sectors, the dam structure appears to deform towards the outside, with velocity vectors approximately radial to the center of the basin, and about the same intensity in both sectors of the basin.
The velocities recorded with A-DInSAR along the dam were compared with those acquired by Jamiolkowski [1], who analyzed inclinometer measurements located along the dam. Along the eastern side of the dam, 26 deep inclinometers were placed and aligned along two cross sections, the XVIE and XVIIIE (located in Figure 6a,b). These measurements highlight the occurrence of horizontal displacements localized at precise depths, where horizontal deformation takes place along sliding planes [1]. According to the author, the movement in the dam foundation is indeed concentrated along sub-horizontal shear planes, formed through glacio-tectonic phenomena in the Pliocene clays. Such planes can be described as ‘horizontal thrusts’ where the shear strength is close to residual, thus controlling the stability of the dam [1,30,51]. Therefore, the inclinometers measured the horizontal displacements of the top layers of Pliocene high plasticity clays, which are induced by tailings load. The inclinometer installed at the benchmark-208 (BM-208 located at the cross-section XVIE in Figure 6a), shows two active shear planes, at elevations 75 m and 45 m asl [1].
In Table 2 are listed the horizontal velocities of displacement at BM208, from 1999 to 2013, using data from Jamiolkowski [1], while from 2014 to 2018, the horizontal velocities of displacement have been retrieved by the A-DInSAR analysis, averaging the velocity for each analyzed year. Such east-west velocities are well comparable with the horizontal rates measured by Jamiolkowski [1] for previous years, i.e., 2011–2013. According to Jamiolkowski [1], towards the end of 1995, when the dam crest had reached the height of 40 m, the rate of horizontal surface displacement increased, so that during 2001–2009 displacement rates were between 40 and 50 mm/year. The horizontal displacements were mitigated by some stabilization measures implemented during 2007–2009 at the dam structure, which resulted in a moderate reduction of movement measured at BM208 (i.e., at the starter dam) of about one-third of the previously observed rate of displacement (Table 2).
The time series of displacement of BM208, from 1982 to 2014, achieved by benchmark time series [1], has been extended to 2018 (Figure 9), using the A-DInSAR measurements. The time series since the second half of the 90s until 2011 shows a steady growth of the eastward movement, with a constant velocity. From 2011 to 2014, and therefore until the end of the period analyzed by Jamiolkosky [1] (orange line in Figure 9), the curve reaches a plateau, and the horizontal displacement slows down. From 2015 to 2019, the A-DInSAR analysis shows instead that the cumulative displacement values of the eastern dam begin to slightly grow again (blue line in Figure 9) at the BM208.
It is worth mentioning that the time series of the displacement of BM208 represents the behavior of the base of the dam structure (starter dam), while at higher altitudes of the dam, the intensity of horizontal displacements is considerably greater, as stated above (Figure 7b,d). As discussed in previous works, the geodetic and the inclinometer data show that the starter dam and the foundation soil above the shallower shear plane in clays (45 m asl) move as a single, rigid body [1]. This behavior changes at higher elevations within the dam, where there is a more complex deformation mechanism.
The longitudinal section along the eastern dam (Section 1 in Figure 6a,b) shows how cumulative displacement is distributed along the eastern perimeter and allows for the analysis of the spatial deformation pattern of the east side of the dam. Each curve of the graph shows the cumulative deformation that affected the dam up to September of each year. This timing was chosen in order to compare and combine this graph with the same longitudinal section of Jamiolkowski [1], which represents a plan view of horizontal cumulated displacements from 2004 to 2013. The latter has been extended with the measurements obtained by A-DInSAR for the years 2014–2018, along the longitudinal section at medium height of the dam (Figure 10). Despite some uncertainties about the exact location of the section of Jamiolkowski, the curves show a very similar trend: greater displacements occur in the central part of the eastern dam, and the rates of eastward displacement (that can be deduced by the distance between the curves) are similar. Two remarks can be made about the above mentioned discussion: the first one is that in the eastern sector of the TSF, the dam does not behave like a rigid body and shows areas of greater weakness in the central part; the second remark is that, although the horizontal deformation seems to have had a slowdown from 2010 to 2013, it continued in 2014 to 2018 with recovered intensity, up to almost 370 mm of displacement. It is possible that stabilization and prevention works performed to reduce the lateral sliding of the dam have not been effective.
In the western sector of the dam (which was not described by Jamiolkowski [1]), the deformation rates are very similar to that of the eastern dam, apart from the direction of the horizontal movement, which is to the west, and the spatial pattern of the displacement along the longitudinal section (Figure 6e), which shows a displacement in bloc of the dam body. Therefore, in this sector there are no localized areas of weakness, as is evident for the eastern sector, and the horizontal movement is uniformly distributed along the perimeter of this part of the dam.
According to the above-mentioned hypothesis that the whole embankment structure is moving outwards, the horizontal displacement of the northern and southern sectors of the dam should be, respectively, northward and southward. As aforementioned, the typical limitation of the InSAR technique in measuring displacements in the North-South direction could therefore lead to an underestimation of north and south displacement of the dam. In support of this hypothesis, the north dam started exhibiting horizontal displacements in 2006/2007, and the rates increased with time, according to Jamiolkowski, [1]. These displacements were localized along two sections of the dam: section VIIN and section XVIN (localized in Figure 6a,b), with a total cumulated displacement since 1987 to 2013 of 283 mm and 408 mm, respectively. Even in the northern sector the horizontal deformation did not occur uniformly but affected mainly these two specific areas. As regards the southern sector of the dam basin, on the other hand, we do not have any precise information on the possible occurrence of southward displacements at the embankment.
The displacement of the dam towards the outside of the dam basin can be produced by two combined factors: (1) the load of the tailings pushing on the containment banks, (2) the sliding along sub-horizontal shear planes at depth, in the high plasticity, over-consolidated Pliocene clays, where the shear strength is close to the residual (details on shear strength parameters can be found in Jamiolkosky [1], and reference therein). The high-plasticity clays exhibit a strong spatial variability, discontinuous pattern, and variable lateral extension [1]. The occurrence of some localized areas of weakness, with greater horizontal deformations, observed in the eastern and in the norther parts of the dam embankment (sections XVIE, XVIN and VIIN in Figure 6a,b), could be linked to peculiar settings of the Pliocene clays at depth; otherwise, it could be linked to structural issues of the dam structure.
Horizontal displacements potentially could cause the dam failure, as the development of progressive failure in a weak soil layer in the dam foundation is indeed considered an important mechanism of failure for tailings dam build up with the upstream method [21]. A failure simulation at a Żelazny Most dam crest elevation of 205 m asl. (or dam height of 91 m), has been developed by [21]. This height is expected to be reached in 2023 if the present filling rate is continued. The failure mode would result in an active rupture mechanism that continues in depth into the Pliocene clays down to the deepest measured shear zones. The authors state that this failure mode is pushing the dam forward along the pre-sheared zones in the clays at the dam foundation, ending with a passive failure mechanism in front of the dam.
A further destabilizing factor for the Żelazny Most TSF stability is represented by the effects of the mining-induced seismicity. The intense seismicity documented at LGCD (e.g., [52,53]) counts thousands of small tremors every year; more than 2000 events with magnitudes between 0.9 and 3.5 are recorded yearly in this district [53]. The earthquakes with magnitude >4 could have an influence on the susceptibility of tailings to failure and can represent a triggering factor for such a system, affected by active and continuous deformation processes. Several mining-induced tremors occurred in the analyzed time span near the dam facility: two M 4.0 events in September and October 2015 struck the mine district at 2.8 km and only 400 m from the dam, respectively, and in 2016, four other major events (ML 4.1, 4.3, and 4.5) occurred within a radius of 9 km from the dam. The closest epicenter was that of the ML 4.3 tremor of 17 October, which occurred about 3.7 km away from the dam. However, these important seismic-induced earthquakes that affected the area near the Żelazny Most TSF caused no observable changes in the dam’s deformation pattern analyzed by A-DInSAR. As described in Jamiolkowski [1], various Norwegian Geotechnical Institute technical reports investigated the effects of mining-induced seismicity on the stability of the East dam and of the tailings, and according to these authors, the impact on the stability of the ring dam has not been considered as relevant.

6. Conclusions

Evaluation of tailings dam stability is a complex geotechnical study, involving the hydrogeological settings assessment, the analysis of the tailing material composition, the deposition technology and the possible triggering factors leading to dam failure, like man-made factors but also including earthquakes or heavy rainfalls.
The monitoring of these delicate and important infrastructures is of strategic importance to prevent damage or failures and to plan tailored risk mitigation measures. Satellite remote sensing technologies, such as satellite SAR Interferometry, represent a consolidate and effective tool for monitoring of structures and infrastructures, including the tailings dams, allowing for analysis of the whole basin, and obtainment of temporal information, using data archives.
The deformation processes at Żelazny Most tailings dam have been analyzed by A-DInSAR technique from October 2014 to April 2019, showing that the storage facility is affected by both vertical deformation (settlements or compaction phenomena) and horizontal sliding towards the outside of the dam. The displacements have been compared with the deformation pattern observed by in situ measurements for the previous years by Jamiolkowski [1], especially for the eastern part of the dam.
Subsidence or settlements in the tailings dam constructed by the upstream method, where the new embankments are founded on tailings, is generally expected as its height increases. The settlement rates of the 2014–2019 period is lower than those described by Jamiolkosky [1] for the previous years, but they have a constant growth trend during the analyzed period. The more affected zones are the upper and the recent parts of the TSF due to both a lower consolidation degree of the recent tailings, and a larger thickness of compressible materials in this inner part of the dam structure.
Horizontal displacements at the Żelazny Most TSF are of primary importance in assessing the stability of the dam. The direction of such horizontal displacements is radial from the center of the dam basin, and in the eastern and northern sectors, it is more severe in some localized areas, while the western sector is moving westward evenly. The outward displacements can be produced by two combined factors: the load of the tailings is pushing the dam forward along horizontal slip zones at depth in the glacio-tectonic surfaces located in the over consolidated Pliocene clay, where the shear strength is close to the residual. This mechanism could end with a passive failure in front of the dam, with a major negative impact on the economy, surrounding properties, and people’s lives aggravated by possible large run-out distances. An accurate monitoring of the storage facility is therefore of crucial importance, especially in the dam sectors where intense deformations have been identified. A proper monitoring plan to mitigate the risk of a dam failure could include the use of Terrestrial SAR interferometry ([54,55,56]), which is considered a proper technology for early warning purposes. Due to its high sampling frequency of data collection (in the order of few minutes or even seconds) and its high accuracy of displacement measurement (up to decimal millimeter order), this technique can be effective for high resolution, local-scale real-time monitoring.
A subsequent development and improvement of this study should consist of the accurate estimate of the local displacement field, including possible N-S movements, overcoming the intrinsic limitations of the 1-D LOS measurements of the satellite SAR interferometry [48,49]. A straightforward approach is the combination of multiple InSAR LOS measurements from different viewing geometries and different angles of incidence [50,57], but it is also possible to resolve the full 3-D surface displacement vectors merging LOS measurements together with other geodetic observations, mainly GPS data [58,59,60].

Author Contributions

Conceptualization P.M., B.A., A.S., F.B.; A-DInSAR processing, S.S.; methodology and data curation, S.S., B.A., A.S., P.M.; formal analysis and investigation, B.A., A.S., S.S.; writing—original draft preparation, B.A.; writing—review and editing, B.A., A.S., P.M., F.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Sapienza Università di Roma, funding program: “Progetti di Ricerca Medi SAPIENZA 2018”, grant number RM118164367B3FF7, title of the project: “Analisi dei rischi geologici interferenti con le grandi dighe attraverso l’uso integrato di dati telerilevati e realtà virtuale”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jamiolkowski, M. Soil mechanics and the observational method: Challenges at the Zelazny Most copper tailings disposal facility. Géotechnique 2014, 64, 590–618. [Google Scholar] [CrossRef] [Green Version]
  2. Lyu, Z.; Chai, J.; Xu, Z.; Qin, Y.; Cao, J. A Comprehensive Review on Reasons for Tailings Dam Failures Based on Case History. Adv. Civ. Eng. 2019, 2019, 1–18. [Google Scholar] [CrossRef]
  3. Kossoff, D.; Dubbin, W.; Alfredsson, M.; Edwards, S.; Macklin, M.; Hudson-Edwards, K. Mine tailings dams: Characteristics, failure, environmental impacts, and remediation. Appl. Geochem. 2014, 51, 229–245. [Google Scholar] [CrossRef] [Green Version]
  4. Rico, M.; Benito, G.; Salgueiro, A.R.; Díez-Herrero, A.; Pereira, H. Reported tailings dam failures. J. Hazard. Mater. 2008, 152, 846–852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Hancock, G. A method for assessing the long-term integrity of tailings dams. Sci. Total Environ. 2021, 779, 146083. [Google Scholar] [CrossRef] [PubMed]
  6. Armstrong, M.; Petter, R.; Petter, C. Why have so many tailings dams failed in recent years? Resour. Policy 2019, 63, 101412. [Google Scholar] [CrossRef]
  7. Bowker, L.N.; Chambers, D.M. The Risk Public Liability & Economics of Tailings Facility Failures. U.S. Department of Agriculture. Tonto National Forest: Payson, AZ, USA, 2015; 56p, Available online: https://www.resolutionmineeis.us/documents/bowker-chambers-2015 (accessed on 1 May 2021).
  8. Azam, S.; Li, Q. Tailings dam failures: A review of the last one hundred years. Geotech. News 2010, 28, 50–53. [Google Scholar]
  9. Tallini, M.; Spadi, M.; Cosentino, D.; Nocentini, M.; Cavuoto, G.; Di Fiore, V. High-resolution seismic reflection exploration for evaluating the seismic hazard in a Plio-Quaternary intermontane basin (L’Aquila downtown, central Italy). Quat. Int. 2019, 532, 34–47. [Google Scholar] [CrossRef]
  10. Stefaniak, K.; Wró, M. On possibilities of using global monitoring in effective prevention of tailings storage facilities failures. Environ. Sci. Pollut. Res. 2018, 25, 5280–5297. [Google Scholar] [CrossRef]
  11. Owen, J.; Kemp, D.; Lèbre, É.; Svobodova, K.; Murillo, G.P. Catastrophic tailings dam failures and disaster risk disclosure. Int. J. Disaster Risk Reduct. 2020, 42, 101361. [Google Scholar] [CrossRef]
  12. Lumbroso, D.; McElroy, C.; Goff, C.; Collell, M.R.; Petkovsek, G.; Wetton, M. The potential to reduce the risks posed by tailings dams using satellite-based information. Int. J. Disaster Risk Reduct. 2019, 38, 101209. [Google Scholar] [CrossRef]
  13. Lumbroso, D.; Davison, M.; Body, R.; Petkovšek, G. Modelling the Brumadinho tailings dam failure, the subsequent loss of life and how it could have been reduced. Nat. Hazards Earth Syst. Sci. 2021, 21, 21–37. [Google Scholar] [CrossRef]
  14. Du, Z.; Ge, L.; Ng, A.H.-M.; Zhu, Q.; Horgan, F.G.; Zhang, Q. Risk assessment for tailings dams in Brumadinho of Brazil using InSAR time series approach. Sci. Total Environ. 2020, 717, 137125. [Google Scholar] [CrossRef] [PubMed]
  15. Carlà, T.; Intrieri, E.; Raspini, F.; Bardi, F.; Farina, P.; Ferretti, A.; Colombo, D.; Novali, F.; Casagli, N. Author Correction: Perspectives on the prediction of catastrophic slope failures from satellite InSAR (Scientific Reports, (2019), 9, 1, (14137), 10.1038/s41598-019-50792-y). Sci. Rep. 2019, 9, 1–9. [Google Scholar] [CrossRef] [Green Version]
  16. Holden, D.; Donegan, S.; Pon, A. Brumadinho Dam InSAR study: Analysis of TerraSAR-X, COSMO-SkyMed and Sentinel-1 images preceding the collapse. In Proceedings of the 2020 International Symposium on Slope Stability in Open Pit Mining and Civil Engineering, Australian Centre for Geomechanics, Crawley, WA, Australia, 12–14 May 2020; pp. 293–306. [Google Scholar]
  17. Pirulli, M.; Barbero, M.; Marchelli, M.; Scavia, C. The failure of the Stava Valley tailings dams (Northern Italy): Numerical analysis of the flow dynamics and rheological properties. Geoenviron. Disasters 2017, 4, 279. [Google Scholar] [CrossRef] [Green Version]
  18. .Chandler, D.W.; Grantham, D.W. Minimum audible movement angle in the horizontal plane as a function of stimulus frequency and bandwidth, source azimuth, and velocity. J. Acoust. Soc. Am. 1992, 91, 1624–1636. [Google Scholar] [CrossRef]
  19. Luino, F.; De Graff, J.V. The Stava mudflow of 19 July 1985 (Northern Italy): A disaster that effective regulation might have prevented. Nat. Hazards Earth Syst. Sci. 2012, 12, 1029–1044. [Google Scholar] [CrossRef] [Green Version]
  20. Sammarco, O. A Tragic Disaster Caused by the Failure of Tailings Dams Leads to the Formation of the Stava 1985 Foundation. Mine Water Environ. 2004, 23, 91–95. [Google Scholar] [CrossRef]
  21. Kalsnes, B.; Jostad, H.P.; Nadim, F.; Hauge, A.; Dutra, A.; Muxfeldt, A.; Mikos, M.; Tiwari, B.; Yin, Y.; Sassa, K. Tailings Dam Stability. In Proceedings of the Advancing Culture of Living with Landslides; Mikos, M., Tiwari, B., Yin, Y., Sassa, K., Eds.; Springer Science and Business Media LLC: Cham, Switzerland, 2017; pp. 1173–1180. [Google Scholar]
  22. Mura, J.C.; Gama, F.F.; Paradella, W.R.; Negrão, P.; Carneiro, S.; De Oliveira, C.G.; Brandão, W.S. Monitoring the Vulnerability of the Dam and Dikes in Germano Iron Mining Area after the Collapse of the Tailings Dam of Fundão (Mariana-MG, Brazil) Using DInSAR Techniques with TerraSAR-X Data. Remote Sens. 2018, 10, 1507. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, T.; Perissin, D.; Rocca, F.; Liao, M. Three Gorges Dam stability monitoring with time-series InSAR image analysis. Sci. China Earth Sci. 2011, 54, 720–732. [Google Scholar] [CrossRef]
  24. Tomás, R.; Cano, M.; García-Barba, J.; Vicente, F.; Herrera, G.; Lopez-Sanchez, J.M.; Mallorquí, J. Monitoring an earthfill dam using differential SAR interferometry: La Pedrera dam, Alicante, Spain. Eng. Geol. 2013, 157, 21–32. [Google Scholar] [CrossRef]
  25. Di Martire, D.; Iglesias, R.; Monells, D.; Centolanza, G.; Sica, S.; Ramondini, M.; Pagano, L.; Mallorquí, J.J.; Calcaterra, D. Comparison between Differential SAR interferometry and ground measurements data in the displacement monitoring of the earth-dam of Conza della Campania (Italy). Remote Sens. Environ. 2014, 148, 58–69. [Google Scholar] [CrossRef]
  26. Milillo, P.; Burgmann, R.; Lundgren, P.; Salzer, J.; Perissin, D.; Fielding, E.; Biondi, F.; Milillo, G. Space geodetic monitoring of engineered structures: The ongoing destabilization of the Mosul dam, Iraq. Sci. Rep. 2016, 6, 37408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Antonielli, B.; Caporossi, P.; Mazzanti, P.; Moretto, S.; Rocca, A. InSAR & Photomonitoringtm for Dams and Reservoir Slopes Health & Safety Monitoring. In Proceedings of the Twenty-Sixth Congress on Large Dams, Commission Internationale des Grands Barrages, Vienna, Austria, 4–6 July 2018; pp. 2565–2577. [Google Scholar] [CrossRef]
  28. Rizzo, P.; Milazzo, A.; Collection, S. Lecture Notes in Civil, Engineering European Workshop on Structural Health Monitoring; Springer International: Palermo, Italy, 6–9 July 2020; Volume 1, ISBN 9783030645939. [Google Scholar]
  29. Jamiolkowski, M.; Carrier, W.D.; Chandler, R.J.; Hoeg, K.; Swierczynski, W.; Wolski, W. 1st Za Chieh-Moh distinguished lecture: The geotechnical problems of the second world largest copper tailings pond at Zelazny Most, Poland. Geotech. Eng. 2010, 41, 1–16. [Google Scholar]
  30. Stefanek, P.; Engels, J.; Wrzosek, K.; Sobiesak, P.; Zalewski, M. Surface tailings disposal at the Żelazny Most TSF, today and into the future. In Proceedings of the 20th International Seminar on Paste and Thickened Tailings, Beijing, China, 15–18 June 2017; pp. 213–225. [Google Scholar]
  31. Butra, J.; Kudełko, J. Rockburst hazard evaluation and prevention methods in Polish copper mines. Cuprum 2011, 61, 5–20. [Google Scholar]
  32. Burtan, Z. The influence of regional geological settings on the seismic hazard level in copper mines in the Legnica-Głogów Copper Belt Area (Poland). E3S Web Conf. 2017, 24, 1004. [Google Scholar] [CrossRef] [Green Version]
  33. Korzec, A.; Świdziński, W. Dynamic response of Zelazny Most tailings dam to mining induced extreme seismic event occurred in 2016. MATEC Web. Conf. 2019, 262, 01001. [Google Scholar] [CrossRef]
  34. Strzelecki, T.; Maniecki, Ł.; Bartlewska-Urban, M. Building a Numerical Model of the Filtration Flow in the Żelazny Most Flotation Tailings Disposal Lake. Stud. Geotech. Mech. 2015, 37, 93–104. [Google Scholar] [CrossRef] [Green Version]
  35. Marks, L.; Dzierżek, J.; Janiszewski, R.; Kaczorowski, J.; Lindner, L.; Majecka, A.; Makos, M.; Szymanek, M.; Tołoczko-Pasek, A.; Woronko, B. Quaternary stratigraphy and palaeogeography of Poland. Acta Geol. Pol. 2016, 66, 410–434. [Google Scholar] [CrossRef]
  36. Lindner, L.; Gozhik, P.; Marciniak, B.; Marks, L.; Yelovicheva, Y. Main climatic changes in the Quaternary of Poland, Belarus and Ukraine. Geol. Q. 2004, 48, 97–114. [Google Scholar]
  37. Badura, J.; Przybylski, B.; Zuchiewicz, W. Cainozoic evolution of lower Silesia, SW Poland: A new interpretation in the light of sub-Cainozoic and sub-Quaternary topography. Acta Geodyn. Geomater. 2004, 1, 7–29. [Google Scholar]
  38. Wiśniewski, A.; Adamiec, G.; Badura, J.; Bluszcz, A.; Kowalska, A.; Kufel-Diakowska, B.; Mikołajczyk, A.; Murczkiewicz, M.; Musil, R.; Przybylski, B.; et al. Occupation dynamics north of the Carpathians and Sudetes during the Weichselian (MIS5d-3): The Lower Silesia (SW Poland) case study. Quat. Int. 2013, 294, 20–40. [Google Scholar] [CrossRef]
  39. Gogolewska, A.B.; Strzeszynska, J. Factors Influencing Rock Burst Hazard in Deep Copper Ore Mine, SW Poland. IOP Conf. Ser. Earth Environ. Sci. 2019, 362, 012023. [Google Scholar] [CrossRef] [Green Version]
  40. Massonnet, D.; Feigl, K.L. Radar interferometry and its application to changes in the Earth’s surface. Rev. Geophys. 1998, 36, 441–500. [Google Scholar] [CrossRef] [Green Version]
  41. Hanssen, R.F. Satellite radar interferometry for deformation monitoring: A priori assessment of feasibility and accuracy. Int. J. Appl. Earth Obs. Geoinf. 2005, 6, 253–260. [Google Scholar] [CrossRef]
  42. Ferretti, A.; Prati, C.; Rocca, F. Permanent scatterers in SAR interferometry. IEEE Trans. Geosci. Remote. Sens. 2001, 39, 8–20. [Google Scholar] [CrossRef]
  43. Kampes, B.M. Radar Interferometry: Persistent Scatterer Technique; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  44. Ferretti, A.; Prati, C.; Rocca, F. Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry. IEEE Trans. Geosci. Remote. Sens. 2000, 38, 2202–2212. [Google Scholar] [CrossRef] [Green Version]
  45. Crosetto, M.; Monserrat, O.; Devanthéry, N.; Cuevas-González, M.; Barra, A.; Crippa, B. Persistent Scatterer Interferometry Using Sentinel-1 Data. ISPRS Int. Arch. Photogramm. Remote. Sens. Spat. Inf. Sci. 2016, XLI-B7, 835–839. [Google Scholar] [CrossRef] [Green Version]
  46. Perissin, D.; Wang, Z.; Wang, T. The SARPROZ InSAR tool for urban subsidence/manmade structure stability monitoring in China. In Proceedings of the ISRSE, Sidney, NSW, Australia, 10–15 April 2011. [Google Scholar]
  47. Samieie-Esfahany, S.; Hanssen, R.F.; Van Thienen-visser, K.; Muntendam-bos, A.; Samiei-Esfahany, S.; Hanssen, R.F.; Thienen-visser, K.; Van Muntendam-bos, A. On the effect of horizontal deformation on InSAR subsidence estimates. In Proceedings of the Proceedings Fringe 2009 Work, Frascati, Italy, 30 November–4 December 2009; pp. 1–7. [Google Scholar]
  48. Wright, T.J.; Parsons, B.E.; Lu, Z. Toward mapping surface deformation in three dimensions using InSAR. Geophys. Res. Lett. 2004, 31, 1–5. [Google Scholar] [CrossRef] [Green Version]
  49. Hu, J.; Li, Z.-W.; Ding, X.; Zhu, J.; Zhang, L.; Sun, Q. Resolving three-dimensional surface displacements from InSAR measurements: A review. Earth Sci. Rev. 2014, 133, 1–17. [Google Scholar] [CrossRef]
  50. Fuhrmann, T.; Garthwaite, M.C. Resolving Three-Dimensional Surface Motion with InSAR: Constraints from Multi-Geometry Data Fusion. Remote Sens. 2019, 11, 241. [Google Scholar] [CrossRef] [Green Version]
  51. Skau, K.S.; Andresen, L.; Jostad, H.P.; Fornes, P.; Grimstad, G.; Page, A. Stability and deformations of Żelazny Most—One of the world’s largest deponies for copper tailings. In Proceedings of the Fjellsprengningsdagen, Bergmekanikkdagen, Geoteknikkdagen, Oslo, Norway, 21–22 November 2013; Tekna. ISBN 978-82-8208-037-8. [Google Scholar]
  52. Lizurek, G.; Lasocki, S. Clustering of mining-induced seismic events in equivalent dimension spaces. J. Seism. 2014, 18, 543–563. [Google Scholar] [CrossRef]
  53. Rudziński, Ł.; Cesca, S.; Lizurek, G. Complex Rupture Process of the 19 March 2013, Rudna Mine (Poland) Induced Seismic Event and Collapse in the Light of Local and Regional Moment Tensor Inversion. Seism. Res. Lett. 2016, 87, 274–284. [Google Scholar] [CrossRef] [Green Version]
  54. Luzi, G. Ground based SAR interferometry: A novel tool for Geoscience. In Geoscience and Remote Sensing, New Achievements; Imperatore, P., Riccio, D., Eds.; InTech: Rijeka, Croatia, 2010; pp. 1–26. [Google Scholar]
  55. Monserrat, O.; Crosetto, M.; Luzi, G. A review of ground-based SAR interferometry for deformation measurement. ISPRS J. Photogramm. Remote Sens. 2014, 93, 40–48. [Google Scholar] [CrossRef] [Green Version]
  56. Mazzanti, P. Toward transportation asset management: What is the role of geotechnical monitoring? J. Civ. Struct. Health Monit. 2017, 7, 645–656. [Google Scholar] [CrossRef]
  57. Raucoules, D.; De Michele, M.; Malet, J.-P.; Ulrich, P. Time-variable 3D ground displacements from high-resolution synthetic aperture radar (SAR). Application to La Valette landslide (South French Alps). Remote Sens. Environ. 2013, 139, 198–204. [Google Scholar] [CrossRef] [Green Version]
  58. Catalao, J.; Nico, G.; Hanssen, R.; Catita, C. Merging GPS and Atmospherically Corrected InSAR Data to Map 3-D Terrain Displacement Velocity. IEEE Trans. Geosci. Remote Sens. 2011, 49, 2354–2360. [Google Scholar] [CrossRef] [Green Version]
  59. Fuhrmann, T.; Cuenca, M.C.; Knöpfler, A.; van Leijen, F.; Mayer, M.; Westerhaus, M.; Hanssen, R.; Heck, B. Estimation of small surface displacements in the Upper Rhine Graben area from a combined analysis of PS-InSAR, levelling and GNSS data. Geophys. J. Int. 2015, 203, 614–631. [Google Scholar] [CrossRef] [Green Version]
  60. Samsonov, S.; Tiampo, K. Analytical Optimization of a DInSAR and GPS Dataset for Derivation of Three-Dimensional Surface Motion. IEEE Geosci. Remote Sens. Lett. 2006, 3, 107–111. [Google Scholar] [CrossRef]
Land 10 00654 g001 550
Figure 1. Żelazny Most Tailings dam and adjacent Legnica Glogow Copper District (LGCD); (a) location of the study area; (b) aerial photo of Żelany Most; (c) location of the tailings dam and borders of LGCD mines in yellow (from http://dm.pgi.gov.pl/, accessed on 2 May 2021); (d) schematic cross-sections of a ring dam (modified form [1], reprinted with permission from [1]. 2014, ICE Publishing).
Figure 1. Żelazny Most Tailings dam and adjacent Legnica Glogow Copper District (LGCD); (a) location of the study area; (b) aerial photo of Żelany Most; (c) location of the tailings dam and borders of LGCD mines in yellow (from http://dm.pgi.gov.pl/, accessed on 2 May 2021); (d) schematic cross-sections of a ring dam (modified form [1], reprinted with permission from [1]. 2014, ICE Publishing).
Land 10 00654 g001
Land 10 00654 g002 550
Figure 2. Section XVIE of the Żelazny Most Tailings dam (eastern part of the dam; modified form [1], reprinted with permission from [1]. 2014, ICE Publishing). Glacio-tectonic shear surfaces are represented in red.
Figure 2. Section XVIE of the Żelazny Most Tailings dam (eastern part of the dam; modified form [1], reprinted with permission from [1]. 2014, ICE Publishing). Glacio-tectonic shear surfaces are represented in red.
Land 10 00654 g002
Land 10 00654 g003 550
Figure 3. STAR graph generated by referring all images to a single master according to the temporal baseline (X axis) and normal baseline (Y axis) for ascending and descending orbital geometries. Every dot represents an image, whereas every line represents an interferogram.
Figure 3. STAR graph generated by referring all images to a single master according to the temporal baseline (X axis) and normal baseline (Y axis) for ascending and descending orbital geometries. Every dot represents an image, whereas every line represents an interferogram.
Land 10 00654 g003
Land 10 00654 g004 550
Figure 4. A-DInSAR velocity of displacement maps (mm/yr) of the Żelazny Most tailings dam for the period October 2014–April 2019; (a,b) velocity in LOS ascending and descending orbital geometry respectively; green PSs are considered as stable; yellow-red PSs are moving away from the satellite; cyan-blue PSs are moving toward the satellite; (c) velocity in east-west direction; cyan-blue PSs are moving eastwards; yellow-red PSs are moving westwards; (d) velocity in up-down direction; yellow-red PSs are moving downward.
Figure 4. A-DInSAR velocity of displacement maps (mm/yr) of the Żelazny Most tailings dam for the period October 2014–April 2019; (a,b) velocity in LOS ascending and descending orbital geometry respectively; green PSs are considered as stable; yellow-red PSs are moving away from the satellite; cyan-blue PSs are moving toward the satellite; (c) velocity in east-west direction; cyan-blue PSs are moving eastwards; yellow-red PSs are moving westwards; (d) velocity in up-down direction; yellow-red PSs are moving downward.
Land 10 00654 g004
Land 10 00654 g005 550
Figure 5. Time series of displacement for the Żelazny Most TSF; (a) time series of up-down displacements for the northern, southern, eastern and western sectors of the dam; downward displacements are in red; (b) time series of east-west displacements for the eastern and western sectors of the dam; eastward displacements are in blue, westward displacement are in red. The location of the PSs selected for processing the time series is also shown in the map panels.
Figure 5. Time series of displacement for the Żelazny Most TSF; (a) time series of up-down displacements for the northern, southern, eastern and western sectors of the dam; downward displacements are in red; (b) time series of east-west displacements for the eastern and western sectors of the dam; eastward displacements are in blue, westward displacement are in red. The location of the PSs selected for processing the time series is also shown in the map panels.
Land 10 00654 g005
Land 10 00654 g006 550
Figure 6. Longitudinal sections of cumulative displacement along the eastern dam (Section 1) and the western dam (Section 2); sections and dam sectors location are overlaid on (a) the east-west velocity map and (b) the up-down velocity map; (c) east-west displacement along the Section 1; (d) up-down displacement along the Section 1; (e) east-west displacement along the Section 2; (f) up-down displacement along the Section 2.
Figure 6. Longitudinal sections of cumulative displacement along the eastern dam (Section 1) and the western dam (Section 2); sections and dam sectors location are overlaid on (a) the east-west velocity map and (b) the up-down velocity map; (c) east-west displacement along the Section 1; (d) up-down displacement along the Section 1; (e) east-west displacement along the Section 2; (f) up-down displacement along the Section 2.
Land 10 00654 g006
Land 10 00654 g007 550
Figure 7. Schematic cross-section XVIE (location in Figure 6a,b) on the eastern dam; vectors represent the cumulative settlement (mm) at different heights of the dam body; (a,b) settlements and east-west displacements, respectively, from 2008 to January 2013 at 7 benchmarks (deduced from [1]); (c,d) settlements and east-west displacements respectively, from October 2014 to April 2019, measured by A-DInSAR, at points P1, P2, P3, and P4 (located in Figure 5).
Figure 7. Schematic cross-section XVIE (location in Figure 6a,b) on the eastern dam; vectors represent the cumulative settlement (mm) at different heights of the dam body; (a,b) settlements and east-west displacements, respectively, from 2008 to January 2013 at 7 benchmarks (deduced from [1]); (c,d) settlements and east-west displacements respectively, from October 2014 to April 2019, measured by A-DInSAR, at points P1, P2, P3, and P4 (located in Figure 5).
Land 10 00654 g007
Land 10 00654 g008 550
Figure 8. Scheme of tailings zoning in the site’s embankment and selected geotechnical parameters of the soil form [30]; 2’, 3’ and 4’ are for fully saturated material; γ: unit weight; φ’: internal friction angle; c’: effective cohesion; AD-InSAR subsidence velocities measured in this work are also reported above the profile.
Figure 8. Scheme of tailings zoning in the site’s embankment and selected geotechnical parameters of the soil form [30]; 2’, 3’ and 4’ are for fully saturated material; γ: unit weight; φ’: internal friction angle; c’: effective cohesion; AD-InSAR subsidence velocities measured in this work are also reported above the profile.
Land 10 00654 g008
Land 10 00654 g009 550
Figure 9. Time series of horizontal surface displacements of benchmarks BM-208; orange line: benchmark time series from [1]; blue line: A-DInSAR time series.
Figure 9. Time series of horizontal surface displacements of benchmarks BM-208; orange line: benchmark time series from [1]; blue line: A-DInSAR time series.
Land 10 00654 g009
Land 10 00654 g010 550
Figure 10. Longitudinal section of cumulative horizontal displacement along the eastern dam (Section 1, location in Figure 6a); curves from 2004 to 2013 are from [1]; curves from 2014 to 2018 are from AD-InSAR measurements.
Figure 10. Longitudinal section of cumulative horizontal displacement along the eastern dam (Section 1, location in Figure 6a); curves from 2004 to 2013 are from [1]; curves from 2014 to 2018 are from AD-InSAR measurements.
Land 10 00654 g010
Table
Table 1. Description of the analyzed dataset.
Table 1. Description of the analyzed dataset.
Sentinel-1AscendingDescending
Nr of scenes206205
Start date15 November 201407 October 2014
End date05 May 201914 April 2019
Look angle (°)39.0538.96
Table
Table 2. Horizontal surface displacement rates of benchmark BM-208 located in the easter sector of the dam; the velocity is an averaged value for each analyzed year.
Table 2. Horizontal surface displacement rates of benchmark BM-208 located in the easter sector of the dam; the velocity is an averaged value for each analyzed year.
Horizontal Velocity of
BM-208 (mm/yr)		YearData Source
28.31999–2000BM-208 (Jamiolkowski, 2014)
19.52000–2001BM-208 (Jamiolkowski, 2014)
41.52001–2002BM-208 (Jamiolkowski, 2014)
37.12002–2003BM-208 (Jamiolkowski, 2014)
37.82003–2004BM-208 (Jamiolkowski, 2014)
382004–2005BM-208 (Jamiolkowski, 2014)
58.92005–2006BM-208 (Jamiolkowski, 2014)
21.82006–2007BM-208 (Jamiolkowski, 2014)
41.42007–2008BM-208 (Jamiolkowski, 2014)
42.32008–2009BM-208 (Jamiolkowski, 2014)
48.62009–2010BM-208 (Jamiolkowski, 2014)
15.32010–2011BM-208 (Jamiolkowski, 2014)
132011–2012BM-208 (Jamiolkowski, 2014)
92012–2013BM-208 (Jamiolkowski, 2014)
-2013–2014No data
12.682014–2015A-DInSAR (this work)
10.62015–2016A-DInSAR (this work)
12.712016–2017A-DInSAR (this work)
10.232017–2018A-DInSAR (this work)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mazzanti, P.; Antonielli, B.; Sciortino, A.; Scancella, S.; Bozzano, F. Tracking Deformation Processes at the Legnica Glogow Copper District (Poland) by Satellite InSAR—II: Żelazny Most Tailings Dam. Land 2021, 10, 654. https://doi.org/10.3390/land10060654

AMA Style

Mazzanti P, Antonielli B, Sciortino A, Scancella S, Bozzano F. Tracking Deformation Processes at the Legnica Glogow Copper District (Poland) by Satellite InSAR—II: Żelazny Most Tailings Dam. Land. 2021; 10(6):654. https://doi.org/10.3390/land10060654

Chicago/Turabian Style

Mazzanti, Paolo, Benedetta Antonielli, Alessandra Sciortino, Stefano Scancella, and Francesca Bozzano. 2021. "Tracking Deformation Processes at the Legnica Glogow Copper District (Poland) by Satellite InSAR—II: Żelazny Most Tailings Dam" Land 10, no. 6: 654. https://doi.org/10.3390/land10060654

APA Style

Mazzanti, P., Antonielli, B., Sciortino, A., Scancella, S., & Bozzano, F. (2021). Tracking Deformation Processes at the Legnica Glogow Copper District (Poland) by Satellite InSAR—II: Żelazny Most Tailings Dam. Land, 10(6), 654. https://doi.org/10.3390/land10060654

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop