Application of a Small Baseline Subset Time Series Method with Atmospheric Correction in Monitoring Results of Mining Activity on Ground Surface and in Detecting Induced Seismic Events

Estimation of Mining-Induced Horizontal Strain Tensor of Land Surface Applying InSAR

Horizontal strains related to mining-induced subsidence may endanger infrastructure and surface users’ safety. While directional horizontal strains should be well determined, appropriate solutions for a complete assessment of the terrain surface deformation field are still required. As a result, the presented study examined a new method for calculating horizontal strain tensor based on the decomposition of satellite radar interferometry (InSAR) observations into vertical and azimuth look direction (ALD) displacements. Based on a geometric integral model, we tested our method on experimental data before applying it to an underground copper ore mine in Poland. In the case study, the displacement field was determined using the Multi-Temporal InSAR method on Sentinel-1 data. The model data relative error did not exceed 0.02 at σ = ±0.003. For the case study, land subsidence of up to −167 mm and ALD displacements ranging from −110 mm to +62 mm was obtained, whereas the extreme values of horizontal strains ranged from −0.52 mm/m to +0.36 mm/m at σ = ±0.050 mm/m. Our results demonstrate the high accuracy of the method in determining the horizontal strain tensor. As a result, the approach can broaden the assessment of the environmental impact of land subsidence worldwide.

Induced Seismic Events—Distribution of Ground Surface Displacements Based on InSAR Methods and Mogi and Yang Models

In this article, we present a possible approach to use satellite radar data for a complete description of the formation process of a subsidence trough resulting from an induced seismic event—a mining tremor. Our main goal was to verify whether SAR data allow for the calculation of the basic indicators for the trough (w—subsidence, T—trough slope, K—curvature, u—horizontal displacements, ε—horizontal deformations). We verified the extent to which the Mogi and Yang models can be fitted to match the actual displacements recorded after an induced seismic tremor. The calculations were performed for the Legnica-Glogow Copper Belt (LGCB) area in southwest Poland. Due to intensive mining operations and specific geological and tectonic conditions, the area shows a high level of induced seismic activity. Our detailed analysis focused on four powerful mining tremors: the first tremor occurred on 29 November 2016 (MW3.4), the second on 7 December 2017 (MW3.3), the next on 26 December 2017 (MW3.6) and the last tremor on 29 January 2019 (MW3.7). For each analyzed event, we determined the displacements based on the Differential Interferometric Synthetic Aperture Radar (DInSAR) method and Sentinel 1 synthetic aperture radar (SAR) data from two paths (22 and 73). Additionally, for the period from November 2014 to October 2020, we calculated the displacements using the Small Baseline Subset method (SBAS) time series method. In all cases, the tremor was followed by the development of long-lasting surface deformations. The obtained results allowed us to conclude that it is possible to calculate indicators that result from a specific induced mining event. Considering the full moment tensor and nature of the tremor source, we demonstrated that the Mogi and Yang models can be employed to describe the influence of an induced tremor on the surface in an area of mining activity. We also confirmed the global character of the influence of the reduced troposphere on SAR data calculations. Our conclusions indicate that accounting for the tropospheric correction does not distort horizontal and vertical displacement values in regions influenced by mining activity/tremors.

Land Surface Subsidence Due to Mining-Induced Tremors in the Upper Silesian Coal Basin (Poland)—Case Study

Seismic phenomena threaten land-based buildings, structures, and infrastructure and can transform land topography. There are two basic types of seismic phenomena, namely, tectonic and anthropogenic, which differ mainly in epicenter depth, surface impact range, and magnitude (energy). This article shows how a land surface was changed by a series of seven rock mass tremors of magnitude ML = 2.3–2.6 in March–May 2017. Their immediate cause was the “momentary” acceleration of void clamping, which was activated by local and short-term seismic phenomena caused by human activity. The induced seismic events resulted from the geological structure of the rock mass, which in the specific region of examination was classified as being highly prone to mining tremors. The authors focused on describing vertical surface displacements in the Upper Silesian Coal Basin in the south of Poland. The surface deformations were identified using DInSAR technology, which allows quasi-continuous monitoring of large areas of land surface. The present research used freely available data from the Copernicus Program and seismic data from the European Plate Observing System.

Application of DInSAR and Spatial Statistics Methods in Analysis of Surface Displacements Caused by Induced Tremors

Induced seismicity is one of the negative phenomena caused by anthropogenic activities that include mining of minerals. This phenomenon manifests itself as sudden and unpredictable shocks of rock mass, which can cause surface deformation and damage to ground infrastructure. Until the advent of satellite radar interferometry that enables analysis of historical events, the characteristics of these unexpected surface deformations were difficult to assess. The main aim of the research was the spatial analysis of the geometry of surface displacements caused by eight induced tremors in the Rudna copper mine (SW Poland) and the dependence of deformation characteristics (vertical displacements, extent) on the induced shock energy. For this purpose, Sentinel-1 satellite imagery, the differential radar satellite interferometry (DInSAR) method and geographic information systems (GIS) based spatial statistics were used. Vertical displacements were mapped on the basis of 37 calculated interferograms. Spatial statistics on the pixel-to-pixel level were performed in the GIS Map Algebra environment. In the result, descriptive and spatial statistics characterizing deformations caused by individual shocks were calculated. The average values of vertical displacements ranged from −44 to −119 mm. Strong, statistical correlation between the extent, maximum vertical displacement, and energy values was determined. In addition, geometries of the formed deformation areas were analyzed and presented graphically. The results obtained in this research constitute development of a knowledge base on surface displacements caused by induced tremors in underground copper mining.

UAV Applications for Determination of Land Deformations Caused by Underground Mining

This article presents a case study that demonstrates the applicability of unmanned aerial vehicle (UAV) photogrammetric data to land surface deformation monitoring in areas affected by underground mining. The results presented include data from two objects located in the Upper Silesian Coal Basin in Poland. The limits of coordinate and displacement accuracy are determined by comparing UAV-derived photogrammetric products to reference data. Vertical displacements are determined based on differences between digital surface models created using UAV imagery from several measurement series. Interpretation problems related to vegetation growth on the terrain surface that significantly affect vertical displacement error are pointed out. Horizontal displacements are determined based on points of observation lines established in the field for monitoring purposes, as well as based on scattered situational details. The use of this type of processing is limited by the need for unambiguous situational details with clear contours. Such details are easy to find in urbanized areas but difficult to find in fields and meadows. In addition, various types of discontinuous deformations are detected and their development over time is presented. The results are compared to forecasted land deformations. As a result of the data processing, it has been estimated that the accuracy of the determination of XY coordinates and the horizontal displacements (RMS) in best case scenario is on the level of 1.5–2 GSD, and about 2–3 GSD for heights and subsidence.

Combined Study of a Significant Mine Collapse Based on Seismological and Geodetic Data—29 January 2019, Rudna Mine, Poland

On 29 January 2019, the collapse of a mine roof resulted in a significant surface deformation and generated a tremor with a magnitude of 4.6 in Rudna Mine, Poland. This study combines the seismological and geodetic monitoring of the event. Data from local and regional seismological networks were used to estimate the mechanism of the source and the ground motion caused by the earthquake. Global Navigation Satellite System data, collected at 10 Hz, and processed as a long-term time-series of daily coordinates solutions and short-term high frequency oscillations, are in good agreement with the seismological outputs, having detected several more tremors. The range and dynamics of the deformed surface area were monitored using satellite radar techniques for slow and fast motion detection. The radar data revealed that a 2-km2 area was affected in the six days after the collapse and that there was an increase in the post-event rate of subsidence.

The Three-Segment Control and Measurement of Reliable Monitoring of the Deformation of the Rock Mass Surface and Engineering Structures on the Międzyodrze Islands in Szczecin, NW Poland

The research area is located in north-western Poland. It is the city of Szczecin with a particular emphasis on the Międzyodrze islands. The area of the EcoGenerator Waste Disposal Plant is part of the research area. The analysis of the geological structure of the subsurface layer of Earth’s crust within Szczecin, was carried out with particular emphasis on the EcoGenerator Waste Disposal Plant. The analysis of height changes of the benchmarks, was based on archival materials measured in two campaigns. A detailed recognition of the geological structure in connection with the analysis of changes in the height of the benchmarks was important. This enabled stable benchmarks to be located in several areas of Szczecin. They formed the basis for reliable monitoring of surface deformations of organic and existing sediments within the EkoGenerator Plant. The application of an appropriate three segment control and measurement system. In the area around the EcoGenerator Plant, vertical movements of the area were observed using the InSAR Small Baseline Subset Method. An InSAR analysis is only used here for very broad identification of the moving area. The radar data came from Sentinel 1 A and 1 B satellites. A total of 129 images from 15.11.2014 to 28.07.2019 were used.The results of the analyses conducted, form the basis for discussion and act as a summary of the considerations in this paper.

An Analysis Applying InSAR of Subsidence Caused by Nearby Mining-Induced Earthquakes

Earthquake occurrence is usually unpredictable apart from sites in the vicinity of volcanoes. It is not easy to measure displacements caused by seismic phenomena using classical geodetic methods, which are based on point survey. Therefore, the surveying of ground movements caused by seismic events should be carried out continuously. Nowadays, remote sensing data and InSAR are often applied to monitor ground displacements in areas affected by seismicity. The effects of severe nearby mining-induced earthquakes have been discussed in the paper. The earthquakes occurred in 2017 and had a magnitude of 4.7 and 4.8. The distance between the epicenters of the mining-induced earthquakes was around 1.6 km. The aim of the investigation has been to analyze the spatio-temporal distribution of ground movements caused by the two tremors using the InSAR technique. Superposition of surface displacement has been studied in time and space. The main scientific aim has been to prove that in the areas where high-energy tremors occur, ground movements overlap. Due to proximity between the epicenters, the mining-induced earthquakes caused the formation of a large subsidence trough with the dimension of approximately 1.2 km × 4.2 km and total subsidence of ca. 116 mm. Two-time phases of subsidence were determined with temporal overlapping. The subsidence analysis has enhanced the cognition of the impact of mining-induced seismicity on the kinematics of surface changes. Moreover, the present work supports the thesis that InSAR is a valuable and adequately accurate technique to monitor ground displacements caused by mining induced earthquakes.