Geophysical Research of Secondary Deformations in the Post Mining Area of the Glaciotectonic Muskau Arch Geopark—Preliminary Results

Abstract

Underground and opencast mining adversely affects the surrounding environment. This process may continue even decades after the end of actual mineral extraction. One of the most significant effects of ceased mining are secondary deformations. Safe, new development of post-mining areas requires reliable information on potential deformation risk zones, which may be difficult to obtain due to a lack of necessary data. This study aimed to investigate and understand the secondary deformation processes in the underground mining area of the former “Babina” lignite mine, located in the unique glaciotectonic environment of the Muskau Arch, in western Poland. A combination of GIS-based historical mapping, geophysical 2D/3D microgravimetry, and Electrical Resistivity Tomography (ERT) measurements allowed the identification of subsidence-prone areas and the determination of potential factors of sinkhole development. The latter are associated with anthropogenic transformation of rock mass and hydrogeological conditions, by shallow underground mining. The results confirmed that multi-level mining of coal deposits in complex and complicated glaciotectonic conditions cause discontinuous deformations, and may be hazardous as long as 50 years after the end of mining operations.

1. Introduction

Post-mining areas, resulting from both underground and open-cast operations, experience negative consequences even decades after the end of mining operations [1,2,3,4,5]. These effects include, among other things: flooding of the subsided areas [6,7]; sinkholes and other forms of discontinuous deformations, as a result of sudden destruction of underground workings, especially shallow ones [8,9]; secondary deformations (uplift), as a result of the restoration of groundwater levels [10,11,12,13]; erosion of mining waste heaps; fires related to coal mining waste heaps; pollution caused by waste heaps [14,15,16,17]; difficulties with the introduction of vegetation and degradation of vegetation conditions due to unfavorable soil characteristics; and subsidence and flooding [18,19,20,21].

These processes pose a threat for new development and use of post-mining areas. As precise information on the location and extent of former excavations is often incomplete or lost, it is difficult to properly plan measures to counteract and protect against the potential, and often sudden, secondary changes of post-mining areas. In addition, as these areas are no longer monitored, e.g., by mine surveying services, the extent and intensity of these processes remain largely unknown.

Secondary subsidence is one of the most threatening processes in post-mining areas. It is deformation of the land surface caused by the readjustment of rocks above voids left by underground mining activity, and is usually sudden. The collapse of a mine working could be the result of factors such as: changes in hydrogeological conditions related to water infiltration into the rock mass; a loss of support for old, shallow excavation workings; weakening of the rock mass properties, as a result of weathering and rheological processes; development of voids, as a result of displacement of backfill material in improperly liquidated shafts; dynamic loading of the rock mass with communication vibrations; or mining tremors [22,23,24,25].

Studies of post-mining areas with a focus on establishing of the extent of former underground mining activity, mapping present-day deformations, modelling and predicting the potential for secondary deformations, and risk assessment of post-mining areas have been conducted, with different measurement, processing, and analytical methods. Land deformation measurements in post-mining areas are carried out with terrestrial (e.g., levelling, geophysical), as well as satellite and/or airborne techniques [8,25,26,27,28,29,30,31,32,33,34].

Geophysical surveys are an important aspect of studies of post-mining areas, as they enable the continuous and non-invasive measurement of specific parameters of the rock mass. These parameters, depending on the chosen technique, are usually determined from the processed data. Ground-based surveys may take one of many forms [35,36], e.g., seismic surveys based on measurement of the propagation of seismic waves produced by human-generated impact or vibration in the ground; gravity surveys using a gravimeter to measure the gravity field to determine variations in rock density; electromagnetic surveys, which involve examination of changes in ground conductivity induced by electromagnetic fields; Ground Penetrating Radar (GPR) measurements; or resistivity methods, based on measurements of the variation of electrical properties of the ground. These methods have been applied for many years to define lithological boundaries and tectonic features, as well as to locate post-mining voids and underground infrastructure in the rock mass, such as drifts or shafts, and the results aid in the assessment of risk of ground subsidence. A review of geophysical methods for the detection and mapping of underground mine workings has been presented by [37]. Noteworthy examples of the application of these methods in studies of post-mining areas include:

• Location of abandoned transport and ventilation shafts in a post-mining lead-zinc area in the USA, with seismic reflection and GPR methods [38];

• Location of underground coal mine workings in India to assess the risk of crown hole subsidence occurrence, with multi-electrode resistivity imaging (RI) and Ground Penetrating Radar (GPR) [39];

• Characterization of derelict coal mine workings in the UK, that included locating an abandoned coal mine shaft [40];

• Detection of post-mining voids using the microgravity technique in karstic environment [41];

• Long-term microgravity and geotechnical monitoring of relict salt mines in the UK [42];

• Estimation of sinkhole occurrence risk in post-mining areas of shallow hard coal exploitation in Poland, with gravimetric and seismic measurements [43];

• Assessment of mine collapse (sinkhole) risk due to voids in underground shafts and the surrounding rock mass [24];

• Detection of abandoned strontianite mine workings in Germany, with a combination of seismic refraction tomography and resistivity tomography [44];

• Identification of voids from underground coal mining which pose a threat of sinkhole occurrence on the surface in Upper Silesia (Poland), with GPR methods [45];

• Mapping of subsidence and sinkholes related to post-mining voids resulting from underground coal mining in South Africa, with High-Resolution Seismic Reflection surveys [46];

• Mapping of abandoned mine underground workings and unconsolidated zones in abandoned coal mining areas in Bulgaria, with application of electrical resistivity tomography (ERT) [47];

• Analysis of the occurrence of voids or loosening zones in the rock mass near an abandoned mining shaft [48];

• Combination of geodetic, GNSS, and geophysical measurements for complex monitoring of mining areas [49].

Geophysical methods have advantages and disadvantages. The main advantage is the possibility of “seeing without drilling”. All of the mentioned geophysical methods provide the distribution of petrophysical parameters below the surface in the form of cross-sections and maps. For studies of post-mining areas, the gravity and seismic methods are the most valuable, because density and/or velocity models can be obtained. Changes in the density of rocks are directly related to the remains of underground works. Resistivity methods are useful for hydrogeology, as electrical parameters depend on the conductivity of the subsurface. Considering the disadvantages of geophysical methods, the problem of resolution should be mentioned here. The more a priori information (for instance, from existing boreholes), the better the obtainable results from the modelled data. However, for so-called “shallow geophysics”, very trustable information can be provided from surface geophysics, without boreholes. As proved in the mentioned publications, a combination of geophysical methods can be used for the development of models describing the underground situation in post-mining areas.

The research presented in this paper aimed at creating a digital inventory of the selected post-mining workings, located in a unique postglacial environment, and mapping the extent of the rock mass disturbed by underground brown coal mining, with the objective of identifying weaker areas susceptible to present-day mine collapse and secondary land deformation. For these purposes, a 3D geospatial database of underground workings, based on available historic mining documentations, was created in GIS, and 2D/3D microgravimetry, followed by electrical resistivity tomography (ERT) geophysical measurements, was performed, and their results were analyzed and interpreted. The study area, methodology, and results, with their discussion, are described in Section 2, Section 3 and Section 4, respectively.

2. Study Area

The study area, which is the former brown coal mine “Friendship of Nations—Babina”, is located in western Poland, close to the border with Germany (Figure 1). The mine was in operation from 1921 until 1973. It is a unique post-mining area, as, initially, the brown coal was mined through a system of shallow underground workings that followed the complex geometry of coal deposits, a result of glaciotectonic processes, and were accessed through a system of inclined shafts. The depth of the underground workings reached up to 100 m below the terrain level. Since the 1960s, the deposits were also mined from the surface in several open pits. The overlying waste rocks were deposited in external and internal waste heaps. After the end of mining, due to economic reasons, the area was subjected to reclamation measures aimed at the restoration of vegetation. These activities included, in the first phase, the levelling of steep slopes in the pits and on the heaps. In the second phase, the pits were filled with water and trees were planted in the remaining area. The second stage was never finished. Details of the reclamation activities are described in [50]. Four study fields were established in the scope of the research project within the border of the former mining field, and the presented study concerns geophysical surveys carried out in test field three (Figure 1). This test field was established in the area of underground mining activity.

In terms of geology and geomorphology, the study area is located in the south-east part of the glaciotectonic structure called the Muskau (Muzakow) Arch. It is a terminal moraine in the shape of a regular arc, open to the north, with an arm span of about 20 km, length of approximately 40 km, and width of around 5–6 km. The central part the arch is divided by the valley of the Nysa Luzycka river, which is also the border between Poland and Germany. The structure was formed as a result of the load exerted by a moving glacial ice sheet during the Marine Oxygen Isotope Stage (MIS) 12, also known as the Elsterian glaciation [51]. The original moraine, that reflected the shape of an ice lobe, was eroded and modelled during the following MIS 6 (Salle) glaciations. In the effects of glacial erosion processes, the height of the moraine was greatly reduced, and deeper glaciotectonic structures, dating back to the MIS 12 period, were exposed in numerous places. The ablation moraines, often perpendicular to the older glacial structures, are a trace of the glacial accumulation from this period.

The Middle and Upper Miocene deposits, deformed during the MIS 12 and MIS 6 glaciations, were tilted and raised from a depth of about 150 m below ground level. The deformation forms of the Paleogene and Neogene deposits are of different sizes, ranging from over a dozen meters to several kilometers. These are mainly continuous compression structures, such as folds and diapirs, as well as discontinuous structures that are dominant in the area, such as overhangs and glaciotectonic scales (Figure 2). These glaciotectonic structures are clearly visible in the present-day relief of the terrain. Distinct geomorphological forms include elongated lowering of the terrain in the zones of disturbed brown coal deposits. Their origin is related to the process of compaction of brown coal caused by weathering. The width of the lowered zones reaches 30 m, with a length of up to 2 km and a depth between 3 and 5 m. They are filled with Holocene alluvial deposits and peats. These subsided forms are divided by deposits of ablation moraines. Together, these represent a stripe-like arrangement of the geological relief and landscape composition [52,53].

The specific study area is test field three, which is located in a glaciotectonic complex called the Nowe Czaple zone, and is dominated by fold-type structures. In the morphology of this area, natural and anthropogenic zones of lowered ground surface are clearly visible. These set the V shaped course of brown coal deposit outcrops, which is most probably a large fold structure consisting of two (northern and southern) arms, with the area between them elevated and made up of quaternary sand and gravel formations. Within the two arms of the V shaped fold structure, multi-level underground mining of brown coal deposits was conducted in the period of 1963–1972. In the present-day topography of the area, above the complex architecture of underground drifts and shafts, elongated subsidence basins, partly covered with water, are clearly visible [53]. The area of test field three is shown in Figure 3.

3. Materials and Methods

The research methodology adopted for the study of the condition of ground disturbed by underground mining in the selected test field consisted of the following stages:

• Development of 3D geodatabase of underground workings;

• Setting out geophysical surveying lines and sites;

• Geophysical 2D microgravimetry surveying of profile lines;

• Geophysical 3D microgravimetry surveying of test sites;

• Geophysical surveying with electrical resistivity tomography;

• Analysis and interpretation of the condition of ground disturbed by underground mining.

3.1. Development of the Geodatabase

The geodatabase was prepared with ESRI ArcGIS Pro software, licensed to the Wroclaw University of Science and Technology, based on archival cartographic documentation comprising of 67 scanned paper maps representing underground and open-pit workings. The scanned maps were obtained from the State Mining Authority archives, and are drawn at a scale of 1:1000. The preserved and available maps do not cover the entire area of the study (all test fields). However, for the analyzed test field (three), the documentation is complete. The scanned and georeferenced maps were used to create vector-type features representing location and 3D geometry of the former underground workings. An example map is shown in Figure 4, and a visualization of the underground workings in this test field is shown in Figure 3. In addition, the location of geological profiles made prior to the start of the mining operation in the area, obtained from the Polish Geological Institute–National Research Institute archives, was included in the geodatabase. Other data used for this study included DEM with a 1 m resolution from the year 2011, as well as Database of Topographical Objects at a scale 1:10,000, both downloaded from the Head Office of Geodesy and Cartography geoportal.

3.2. 2D and 3D Geophysical Microgravimetry Surveys

Microgravimetry is a method used for the search and recognition of shallow underground forms of small size, often of anthropogenic origin. The results of microgravimetric studies are used to detect voids and discontinuous deformations. Conclusions concerning geological structure of the studied area are developed in the form of maps and density models showing subsurface distribution of rock masses. Details of the method can be found in [37,41,42,54], and example applications of the method in studies of post-mining areas are stated in the Introduction section.

Based on the analysis of the location of former underground operations, in relation to present-day topography and archival geological profiles, a geophysical survey network for 2D and 3D microgravimetry surveys was designed. This included profile lines that were set out at an interval of approximately 100 m, and approximately 10 m spacing in a given line in the area of test field three. Furthermore, three detailed test sites were set out for 3D gravimetric surveys in a quasi-regular 6 × 6 m grid, including one site where a sinkhole was spotted in late 2019. The location of the profile lines and test sites is shown in Figure 5.

The gravimetric field surveys consisted of the establishment of 4917 measurement points and GNSS measurements of their location; 3900 points were located along the 2D profiles, and 1017 in the detailed test sites. For the microgravimetry surveys, two Scintrex CG-5 Autograv gravimeters were used. The gravimetric measurements, at scattered points, were made with the method of intermediate points. Thus, each measurement sequence began and ended with measurement on a base point with a known value of the acceleration of gravity g. To determine the accuracy of measurements, acceleration of gravity g values was determined twice at selected positions, distributed evenly over the entire area of study. The mean square error of a single measurement amounted to ±0.006 mGal. For data processing, Schlumberger Petrel (Grav&Mag plug in) and ZOND software were used.

The results of interpretation of the obtained geophysical data are presented in the form gravimetric maps in the Results and Discussion section. The prepared gravimetric maps highlighted the geostructural forms, and helped to define the range of possible land deformations that were further studied with the Electrical Resistivity Tomography (ERT) surveys.

3.3. Electrical Resistivity Tomography Surveys

The purpose of electrical surveys is to determine the subsurface resistivity distribution by making measurements on the ground surface. From these measurements, the true resistivity of the subsurface can be estimated. The ground resistivity is related to various geological parameters, such as mineral and fluid content, porosity, and degree of rock water saturation.

Resistivity measurements, or so-called resistivity sounding, are normally made by injecting a current into the ground through two current electrodes and measuring the resulting voltage difference at two potential electrodes. From the current I and voltage V values, an apparent resistivity ${\rho }_a\left[\mathsf{\Omega }m\right]$ is calculated:

(1)${\rho }_a=k\left(\frac{V}{I}\right),$

where k is the geometric factor, which depends on the arrangement of the electrodes. Using measurements of I and V for different k, we can estimate the resistivity changes in vertical direction. A more accurate model of the subsurface is a two-dimensional (2D) model, where the resistivity changes in the vertical direction as well as in the horizontal direction, along the survey line. Details of this method can be also found in [55,56,57], and example applications of the method in studies of post-mining areas are stated in the Introduction section.

In our case, ERT measurements were conducted in 15 profiles with a total length of 4413 m. The location of profile is shown in Figure 5. In the surveys, an asymmetric three-electrode scheme was used, with the spacing of the electrodes along the profile equal to 3 or 5 m. The scheme is presented in Figure 6. Data were recorded with the ARES II 850 W Automatic Resistivity System. The adopted methodology aimed for semi-detailed reconnaissance of the rock mass up to a depth of 70–80 m below surface level, in the central part of test field three, and detailed identification of test sites to a depth of 30–40 m below surface level. The collected data were processed by 2D inversion in profile lines. As a result, resistivity models of the ground were obtained. Electrical resistivity tomography is most often used in the 2D variant, which produces cross-sections showing the variability of resistivity along the measuring profile together with its depth. The results of ERT measurements are presented as profiles of apparent resistivity with geophysical interpretation.

4. Results and Discussion

The development of a 3D geospatial database of former mining activity allowed us to map the location and geometry of underground workings with respect to the present-day topography of the study area. The GIS-based geodatabase stores features representing individual vertical, horizontal, and inclined drifts, with information on their depth below the surface, heading, and geometry. The archival documentations, and thus the resulting database, cover the entire area of test field three. These were used to set out test field three and tests sites within its border.

Interpreting the gravimetry results, the areas at the greatest risk of post-mining deformations, related to the consolidation of the ground caused by mining, are associated with the occurrence of negative gravity anomalies resulting from the shortage of mass in the substratum.

The results of microgravimetry data interpretation (Figure 7) show that the Bouguer anomaly has a streaked, alternating pattern visible in the central and northern parts of the test field. From the center of the area to the south, there is a larger, relatively negative anomaly, reaching a minimum of −14.21 mGal at the southern border, near the anthropogenic lake above the excavated coal deposit. There are three streaks of relatively negative anomalies from the center in the NE part. These negative anomalies are associated with the lighter formations, also of anthropogenic origin, that fill the depressions. The positive anomalies of the SW-NE direction run from the center and continue beyond the northern and north-eastern boundaries of the area. The positive anomalies are probably caused by tertiary formations, with a higher bulk density in relation to the anthropogenic formations marked in the picture of Bouguer anomalies, in the form of relatively negative anomalous forms. The picture of the Bouguer anomaly in, test field three, is shown in Figure 7a. To focus on anomalies associated with shallow sources, separation of gravity Bouguer anomalies into residual components was done in the frequency domain, using the Butterworth high-pass filter. This filter rejects regional components resulting from deeper geological structures.

The image of residual anomalies is dominated by alternating positive and negative anomalous forms running from the central part of the area towards the north-east (Figure 7b). The two larger forms in the central part merge into one anomaly towards the SW (A). The anomaly near the NE border has three maxima, the largest of which reaches 0.28 mGal. In the southern part, near the anthropogenic lake (B), there are several small positive anomalies with values that do not exceed 0.15 mGal. There are three negative bands between the positive bands with a minimum value of −0.27 mGal (C). The sequences of positive values (in the central and north-eastern part) are contoured by distinct densities separating the areas of positive and negative anomalies. This is probably related to the difference in density between the heavier Miocene formations and lighter anthropogenic formations and brown coal outcrops. In the central and southern parts of test field three, the density contrast is slightly lower, and may indicate transformation of the original environment as a result of the brown coal mining and appearance of anthropogenic deposits (sand, gravel, peat, or silt). With respect to the particular test sites, the one located in the NE part of test field three is covered by a relatively negative residual anomaly, with a minimum of −0.055 mGal (D), located near the south-eastern border of the site. The positive anomalous forms are visible mainly at the edges, and partially in the central part, of this area. This pattern of anomalies is related to the post-mining character of the test site, i.e., the presence of depressions filled with anthropogenic formations of lower density than the surrounding Miocene forms. In the test site located in the SE part of test field three (E), the distribution of residual anomalies shows an extensive, relatively negative anomaly running from SW to NE. The lowest value of −0.05 mGal is located in the southern part of this site. The image of residual anomalies, as in the case of the previous test site, reflects the post-mining nature of the area, where, next to higher bulk density tertiary deposits, there are lighter anthropogenic formations. The last test site (F) shows the same character.

The 2D and 3D gravimetric modeling methods made it possible to quantify the image of the rock mass in the light of gravimetric data. The Bouguer anomaly maps confirmed the glaciotectonic character of the geological structure, emphasized the geostructural forms, and helped to define the extent of possible land deformations. The areas exposed to the threat of surface deformations are related to the consolidation of the ground disturbed by mining. In gravimetric images, these are associated with the occurrence of negative anomalies, and result from a shortage of mass. The process of tightening the goafs (underground workings) loosens the sediments above, and thus reduces their density in relation to the intact surroundings. These sites were further investigated with the Electrical Resistivity Tomography surveys. The results, in the form of eight apparent resistivity cross-sections beneath the ERT geophysical profiles, and their geophysical interpretation, indicate three distinct lithological complexes, spatially differing in resistivity and physical-mechanical properties (Table 1). The delineation of boundaries between individual complexes is based only on the values of the physical parameter, which is the resistivity of the geological medium. Their spatial extent, in an exemplary profile corresponding to the profile line marked in Figure 5, is shown in Figure 8.

Table 1

Lithological and physic-mechanical classification of resistance complexes in the study area.

Complex	Range of Resistivity [Ωm]	Filtration Coefficient [m/s]	Lithological Classification
loam	<40	Impermeable to poorly permeable (1 × 10−8–1 × 10−6)	Cohesive formations-clay formations (silts, clays, loams) and brown coal
sand-loam	40–200	poorly permeable to permeable(1 × 10−6–1 × 10−4)	sandy loams (values between 100 and 200 Ωm) indicate predominance of non-cohesive formations (sands, dusty sands, and gravel in the saturation zone).
sand	>200	Permeable (1 × 10−4–1 × 10−3)	Sands and gravels in the aeration zone; anthropogenic formations.

The ranges of resistivity values of coals versus the typical resistivity values of minerals and rocks (according to [58]).

Based on the adopted criteria, a schematic diagram of the hydrogeological cross-section, with marked extent of deformation zones caused by the exploitation or weathering of brown coal seams (Figure 9), was constructed. The discontinuity zones of poorly permeable and impermeable grounds constitute potential channels for the migration of non-cohesive formations into mining excavations zones. Based on the analysis of these results in ERT profiles, a map of subsidence risk areas was interpreted, and further geophysical studies provided additional geological information and identified locations of potential ground disturbance of anthropogenic (mining) origin (Figure 10).

The presented results provide information about conditions that can be transferred to other fields of underground mining in the study area, in order to map zones of potential land deformation (sinkhole) occurrence, such as in test field three (Figure 11). The location of this sinkhole is shown in Figure 5 and Figure 10. The feature was spotted in Autumn 2019 by local forest services, and has been under monitoring with a terrestrial laser scanner since the Spring of 2020. The results of the most recently processed scans (3.2021) indicate an average sinkhole depth of 4 m below surface level.

This is the first study using geophysical methods in this post-mining area, and probably the first one conducted in a structurally complex glaciotectonic environment. The only known description of mining-related deformations in a comparable area was prepared by [59], and concerns simple inventorying and description of different deformation types, e.g., subsidence basins and sinkholes in former brown coal mines located to the north of the study area. Different types of deformations were also briefly mentioned in [50]. A study of deformations caused by ceased shallow underground brown coal mining, conducted with the application of geophysical and geological surveys in different geological and tectonic settings, was described in [60]. There are numerous studies using different geophysical methods to investigate and map voids resulting from old mining activity, and aimed at identifying potential deformation risk areas. Notable examples are listed in the Introduction section. None of them, however, were conducted in a glaciotectonic area and concerned shallow underground mining of brown coal deposits in complicated conditions of mineral extraction.

5. Conclusions

The glaciotectonic Muskau Arch structure was disturbed by long-time mining of brown coal deposits. The effects of past mining in complicated geological conditions, such as secondary deformations, pose a threat to present day development of this area.

This study mapped, for the first time, the old underground workings in relation to present day topography, and identified zones of disturbed ground that are susceptible to the formation of sinkholes.

The results of gravity and ERT geophysical measurements revealed deep anthropogenic transformations of the natural geological conditions, and identified three distinct complexes in the ground.

Furthermore, a model of present-day hydrogeological conditions was proposed in the area of past underground mining. The natural system of separate aquifers there has been transformed, and forms one connected groundwater system.

These results relate the existing sinkholes with shallow underground workings and zones of anthropogenically disturbed ground. The periodic and seasonal fluctuation of the groundwater table causes subsequent watering and dewatering of upper parts of the ground and old mine workings, that activate land deformation processes.

The adopted methodology can be transferred to other parts of underground brown coal mining in the Muskau Arch area, as well as to other regions of shallow mining, and confirms the effectiveness of the adopted approach in determining subsidence risk zones in glaciotectonic areas.

Author Contributions

Conceptualization, J.B., E.W. and J.K.; methodology, J.B., E.W. and L.K.; software, J.B., P.T., J.W., M.B. and A.B.; validation, J.B., E.W., M.B. and M.W.; formal analysis, M.W., P.T., E.W. and L.K.; investigation, N.B., D.J. and M.B.; resources, J.K., L.K. and J.W.; data curation, P.T., N.B., A.B. and J.W.; writing—original draft preparation, J.B., E.W. and M.W.; writing—review and editing, A.B., N.B. and P.K.; visualization, J.B., A.B., J.K. and P.T.; supervision, J.B. and E.W.; project administration, J.B. and E.W.; funding acquisition, J.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financed by the National Science Centre financed OPUS-17 project, entitled “Genesis and course of anthropogenic and natural deformations of the terrain in post-mining areas of the former brown coal mine Babina”, no. 2019/33/B/ST10/02975.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors would like to express their gratitude to the “Lipinki” State Forest District Service for their assistance during the carrying out of the research.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ALOS	Advanced Land Observation Satellite
DEM	Digital Elevation Model
ERT	Electrical Resistivity Tomography
GIS	Geographic Information System
GNSS	Global Navigation Satellite Systems
GPR	Ground Penetrating Radar
InSAR	Interferometric Synthetic Aperture Radar
MIS	Marine Oxygen Isotope Stage
RI	Resistivity Imaging

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Table

Figure 1. General location of the study area: (a) location of the Muskau Arch (b) test field three; (c) extent of the Muskau Arch.

Figure 2. General ideological model of the deformation and erosional processes in the area of the transboundary Muskau Arch ([51] after [53]): (A) The approximately 430–510 m thick glacier dislocates the horizontal layered bed of tertiary and quaternary deposits by loading and horizontal push. The glacier acts down to a depth of about 270–290 m beneath the ice base. 1—terminal moraine; 2—fold and diapir structures; 3—scale structures; 4—level of moraine erosion. (B) The contemporary shape of the coal seam outcrop, depending on the type of glaciotectonic structures. 5—brown coal deposits associated with folds and diapirs; 6—brown coal deposits associated with scales; 7—present day level of terminal moraine erosion.

Figure 3. Area of test field three in relation to topography and documented underground workings; left, topographical background; right, orthophoto mosaic background.

Figure 4. An example of source documentation; a section of mining map depicting underground workings.

Figure 5. Location of geophysical surveys in relation to topography and documented underground workings in test field three.

Figure 6. Scheme of ERT measurement with a three-electrode system method.

Figure 7. Gravimetric map: (a) Gravimetric anomalies in the Bouguer reduction; (b) Residual anomalies.

Figure 8. Exemplary geoelectrical profile with geophysical interpretation.

Figure 9. Diagram of interpretation of the hydrogeological conditions along profile (shown in Figure 5) in test-field three.

Figure 10. Interpreted subsidence risk areas in zones of ground disturbed by mining activity.

Figure 11. Photo of the recent sinkhole in the study area.

References

1. Lin, P.M.; Peng, S.S.; Tsang, P. Dealing with Subsidence on Abandoned Mine Lands. Min. Eng.; 1990; 42, 1245.

2. Preuße, A.; Kateloe, H.-J.; Sroka, A. Subsidence and Uplift Prediction in German and Polish Hard Coal Mining. Proceedings of the 31st International Conference on Ground Control in Mining; Morgantown, WV, USA, 31 July–2 August 2012.

3. Contrucci, I.; Balland, C.; Kinscher, J.; Bennani, M.; Bigarré, P.; Bernard, P. Aseismic Mining Subsidence in an Abandoned Mine: Influence Factors and Consequences for Post-Mining Risk Management. Pure Appl. Geophys.; 2019; 176, pp. 801-825. [DOI: https://dx.doi.org/10.1007/s00024-018-2015-6]

4. Strzałkowski, P.; Ścigała, R. Assessment of Post-mining Terrain Suitability for Economic Use. Int. J. Environ. Sci. Technol.; 2020; 17, pp. 3143-3152. [DOI: https://dx.doi.org/10.1007/s13762-019-02617-8]

5. Carrión-Mero, P.; Montalván-Burbano, N.; Herrera-Narváez, G.; Morante-Carballo, F. Geodiversity and mining towards the development of Geotourism: A global perspective. Int. J. Des. Nat. Ecodynamics; 2021; 16, pp. 191-201. [DOI: https://dx.doi.org/10.18280/ijdne.160209]

6. Salmi, E.F.; Malinowska, A.; Hejmanowski, R. Investigating the Post-Mining Subsidence and the Long-Term Stability of Old Mining Excavations: Case of Cow Pasture Limestone Mine, West Midlands, UK. Bull. Eng. Geol. Environ.; 2020; 79, pp. 225-242. [DOI: https://dx.doi.org/10.1007/s10064-019-01575-2]

7. Padmanaban, R.; Bhowmik, A.; Cabral, P. A Remote Sensing Approach to Environmental Monitoring in a Reclaimed Mine Area. Int. J. Geo-Inf.; 2017; 6, 401. [DOI: https://dx.doi.org/10.3390/ijgi6120401]

8. Malinowska, A.; Witkowski, W.; Hejmanowski, R.; Chang, L.; van Leijen, F.; Hanssen, R. Sinkhole Occurrence Monitoring over Shallow Abandoned Coal Mines with Satellite-Based Persistent Scatterer Interferometry. Eng. Geol.; 2019; 262, 105336. [DOI: https://dx.doi.org/10.1016/j.enggeo.2019.105336]

9. De Vent, I.; Roest, H. Lagging Mining Damage in the Netherlands? Recent Signs of Soil Movement in the Zuid–Limburg Coal District. Proceedings of the XV International ISM Congress; Aachen, Germany, 16–20 September 2013; pp. 27-41.

10. Cuenca, M.C.; Hooper, A.J.; Hanssen, R.F. Surface Deformation Induced by Water Influx in the Abandoned Coal Mines in Limburg, The Netherlands Observed by Satellite Radar Interferometry. J. Appl. Geophys.; 2013; 88, pp. 1-11. [DOI: https://dx.doi.org/10.1016/j.jappgeo.2012.10.003]

11. Vervoort, A.; Declercq, P.-Y. Surface Movement above Old Coal Longwalls after Mine Closure. Int. J. Min. Sci. Technol.; 2017; 27, pp. 481-490. [DOI: https://dx.doi.org/10.1016/j.ijmst.2017.03.007]

12. Blachowski, J.; Kopeć, A.; Milczarek, W.; Owczarz, K. Evolution of Secondary Deformations Captured by Satellite Radar Interferometry: Case Study of an Abandoned Coal Basin in SW Poland. Sustainability; 2019; 11, 884. [DOI: https://dx.doi.org/10.3390/su11030884]

13. Malinowska, A.A.; Witkowski, W.T.; Guzy, A.; Hejmanowski, R. Satellite-Based Monitoring and Modeling of Ground Movements Caused by Water Rebound. Remote Sens.; 2020; 12, 1786. [DOI: https://dx.doi.org/10.3390/rs12111786]

14. Frank, O. Aspects of Surface and Environment Protection in German Mining Areas. Min. Sci. Tech.; 2009; 19, pp. 615-619. [DOI: https://dx.doi.org/10.1016/S1674-5264(09)60114-9]

15. Nyssen, J.; Vermeersch, D. Slope Aspect Affects Geomorphic Dynamics of Coal Mining Spoil Heaps in Belgium. Geomorphology; 2010; 123, pp. 109-121. [DOI: https://dx.doi.org/10.1016/j.geomorph.2010.07.004]

16. Adibee, N.; Osanloo, M.; Rahmanpour, M. Adverse Effects of Coal Mine Waste Dumps on the Environment and Their Management. Environ. Earth Sci.; 2013; 70, pp. 1581-1592. [DOI: https://dx.doi.org/10.1007/s12665-013-2243-0]

17. Abramowicz, A.; Rahmonov, O.; Chybiorz, R. Environmental Management and Landscape Transformation on Self-Heating Coal-Waste Dumps in the Upper Silesian Coal Basin. Land; 2020; 10, 23. [DOI: https://dx.doi.org/10.3390/land10010023]

18. Voigt, S.; Tetzlaff, A.; Zhang, J.; Kuenzer, C.; Zhukov, B.S.; Strunz, G.; Oertel, D.; Roth, A.; van Dijk, P.M.; Mehl, H. Integrating Satellite Remote Sensing Techniques for Detection and Analysis of Uncontrolled Coal Seam Fires in North China. Int. J. Coal Geol.; 2004; 59, pp. 121-136. [DOI: https://dx.doi.org/10.1016/j.coal.2003.12.013]

19. Erener, A. Remote Sensing of Vegetation Health for Reclaimed Areas of Seyitömer Open Cast Coal Mine. Int. J. Coal Geol.; 2011; 86, pp. 20-26. [DOI: https://dx.doi.org/10.1016/j.coal.2010.12.009]

20. Karan, S.K.; Samadder, S.R.; Maiti, S.K. Assessment of the Capability of Remote Sensing and GIS Techniques for Monitoring Reclamation Success in Coal Mine Degraded Lands. J. Environ. Manag.; 2016; 182, pp. 272-283. [DOI: https://dx.doi.org/10.1016/j.jenvman.2016.07.070] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27491028]

21. Vorovencii, I. Changes Detected in the Extent of Surface Mining and Reclamation Using Multitemporal Landsat Imagery: A Case Study of Jiu Valley, Romania. Environ. Monit. Assess.; 2021; 193, 30. [DOI: https://dx.doi.org/10.1007/s10661-020-08834-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33398530]

22. Chwastek, J. Possibility of the Impact of Old Mining Excavations in Lower Silesia on the Development of the Land Surface; Center for Forecasting Research, Wrocław University of Science and Technology: Wrocław, Poland, 1974; Volume 3, (In Polish)

23. Kratzsch, H. Mining Subsidence Engineering; Springer: Berlin/Heidelberg, Germany, 1983; ISBN 978-3-642-81925-4

24. Pilecki, Z.; Popiołek, E. Geodetic and Geophysical Identification of Collapsing Hazards. Work. Saf. Environ. Prot. Min.; 2010; 6, pp. 34-39. (In Polish)

25. Suh, J.; Choi, Y. Mapping Hazardous Mining-Induced Sinkhole Subsidence Using Unmanned Aerial Vehicle (Drone) Photogrammetry. Environ. Earth Sci.; 2017; 76, 144. [DOI: https://dx.doi.org/10.1007/s12665-017-6458-3]

26. Zheng, M.; Zhang, H.; Deng, K.; Du, S.; Wang, L. Analysis of Pre- and Post-Mine Closure Surface Deformations in Western Xuzhou Coalfield From 2006 to 2018. IEEE Access; 2019; 7, pp. 124158-124172. [DOI: https://dx.doi.org/10.1109/ACCESS.2019.2938846]

27. Muntean, A.; Mocanu, V.; Ambrosius, B. A GPS Study of Land Subsidence in the Petrosani (Romania) Coal Mining Area. Nat. Hazards; 2016; 80, pp. 797-810. [DOI: https://dx.doi.org/10.1007/s11069-015-1997-y]

28. Tichavský, R.; Jiránková, E.; Fabiánová, A. Dating of Mining-Induced Subsidence Based on a Combination of Dendrogeomorphic Methods and in Situ Monitoring. Engineering Geology. Eng. Geol.; 2020; 272, 105650. [DOI: https://dx.doi.org/10.1016/j.enggeo.2020.105650]

29. Akcin, H.; Kutoglu, H.S.; Kemaldere, H.; Deguchi, T.; Koksal, E. Monitoring Subsidence Effects in the Urban Area of Zonguldak Hardcoal Basin of Turkey by InSAR-GIS Integration. Nat. Hazards Earth Syst. Sci.; 2010; 10, pp. 1807-1814. [DOI: https://dx.doi.org/10.5194/nhess-10-1807-2010]

30. Lee, S.; Park, I. Application of Decision Tree Model for the Ground Subsidence Hazard Mapping near Abandoned Underground Coal Mines. J. Environ. Manag.; 2013; 127, pp. 166-176. [DOI: https://dx.doi.org/10.1016/j.jenvman.2013.04.010]

31. Oh, H.-J.; Lee, S. Assessment of Ground Subsidence Using GIS and the Weights-of-Evidence Model. Eng. Geol.; 2010; 115, pp. 36-48. [DOI: https://dx.doi.org/10.1016/j.enggeo.2010.06.015]

32. Sedlák, S.; Poljakovič, P. Particularities of Deformation Processes Solution with GIS Application for Mining Landscape Reclamation in East Slovakia. J. Geogr. Cartogr.; 2018; 4, 508. [DOI: https://dx.doi.org/10.24294/jgc.v1i4.508]

33. Guo, Q.; Meng, X.; Li, Y.; Lv, X.; Liu, C. A Prediction Model for the Surface Residual Subsidence in an Abandoned Goaf for Sustainable Development of Resource-Exhausted Cities. J. Clean. Prod.; 2021; 279, 123803. [DOI: https://dx.doi.org/10.1016/j.jclepro.2020.123803]

34. Akcin, H. A GIS-Based Building Risk Assessment for the Subsidence Due to Undercity Coal Mining Activities in Zonguldak, Turkey. Arab. J. Geosci.; 2021; 14, 376. [DOI: https://dx.doi.org/10.1007/s12517-021-06702-6]

35. Kotyrba, A. Application of Geophysical Methods in Geoengineering Research in Mining and Post-Mining Areas. Proceedings of the 3rd Conference Geology, Hydrogeology and Geophysics in Solving the Problems of Modern Mining, Energy and the Environment; Ustroń, Poland, June 2012; (In Polish)

36. Majer, E.; Sokołowska, M.; Frankowski, Z.; Barański, M.; Bestyński, Z.; Ostrowski, S.; Pasieczna, A.; Pietrzykowski, P.; Przyłucka, M.; Błachnio, O. et al. Principles of Geological and Engineering Documentation; Polish Geological Institute: Warszawa, Poland, 2018; ISBN 978-83-7863-774-5 (In Polish)

37. Johnson, W.J.; Snow, R.E.; Clark, J.C. Surface Geophysical Methods for the Detection of Underground Mine Workings. Proceedings of the Symposium on Geotechnical Methods for Mine Mapping Verifications; Charleston, WV, USA, 29 October 2002.

38. Shoemaker, M.L.; Anderson, N.L.; Hatheway, A.W.; Newton, T.E.; Shaw, A.E.; Baker, J.A.; Cardimona, S.; Webb, D.J.; Conley, J.L. Reflection Seismic and Ground-penetrating Radar Study of Previously Mined (Lead/Zinc) Ground, Joplin, Missouri. J. Environ. Eng. Geophys.; 1999; 4, pp. 105-112. [DOI: https://dx.doi.org/10.4133/JEEG4.2.105]

39. Soni, A.K.; Singh, K.K.K.; Prakash, A.; Singh, K.B.; Chakraboraty, A.K. Shallow Cover over Coal Mining: A Case Study of Subsidence at Kamptee Colliery, Nagpur, India. Bull. Eng. Geol. Environ.; 2007; 66, pp. 311-318. [DOI: https://dx.doi.org/10.1007/s10064-006-0072-z]

40. Pringle, J.K.; Stimpson, I.G.; Toon, S.M.; Caunt, S.; Lane, V.S.; Husband, C.R.; Jones, G.M.; Cassidy, N.J.; Styles, P. Geophysical Characterization of Derelict Coalmine Workings and Mineshaft Detection: A Case Study from Shrewsbury, United Kingdom. Near Surf. Geophys.; 2008; 6, pp. 185-194. [DOI: https://dx.doi.org/10.3997/1873-0604.2008014]

41. Bishop, I.; Styles, P.; Emsley, S.J.; Ferguson, N.S. The Detection of Cavities Using the Microgravity Technique: Case Histories from Mining and Karstic Environments. Geol. Soc. Lond. Eng. Geol. Spec. Publ.; 1997; 12, pp. 153-166. [DOI: https://dx.doi.org/10.1144/GSL.ENG.1997.012.01.13]

42. Pringle, J.K.; Styles, P.; Howell, C.P.; Branston, M.W.; Furner, R.; Toon, S.M. Long-Term Time-Lapse Microgravity and Geotechnical Monitoring of Relict Salt Mines, Marston, Cheshire, U.K. Geophysics; 2012; 77, pp. B287-B294. [DOI: https://dx.doi.org/10.1190/geo2011-0491.1]

43. Pilecki, Z. The Role of Geophysical Methods in the Estimation of Sinkhole Threat in the Post-Mining Areas of Shallow Exploitation in the Upper Silesian Coal Basin, Poland. Miner. Resour. Manag.; 2008; 24, pp. 27-40.

44. Kotyrba, B.; Schmidt, V. Combination of Seismic and Resistivity Tomography for the Detection of Abandoned Mine Workings in Münster/Westfalen, Germany: Improved Data Interpretation by Cluster Analysis. Near Surf. Geophys.; 2014; 12, pp. 415-426. [DOI: https://dx.doi.org/10.3997/1873-0604.2013056]

45. Tomecka-Suchoń, S.; Marcak, H. Interpretation of Ground Penetrating Radar Attributes in Identifying the Risk of Mining Subsidence. Arch. Min. Sci.; 2015; 60, pp. 645-656. [DOI: https://dx.doi.org/10.1515/amsc-2015-0042]

46. Isiaka, A.I.; Durrheim, R.J.; Manzi, M.S.D. High-Resolution Seismic Reflection Investigation of Subsidence and Sinkholes at an Abandoned Coal Mine Site in South Africa. Pure Appl. Geophys.; 2019; 176, pp. 1531-1548. [DOI: https://dx.doi.org/10.1007/s00024-018-2026-3]

47. Dimovski, S.; Stoyanov, N.; Tzankov, C.; Kisyov, A. A Geophysical Approach for Mapping of Abandoned Mining Workings and Unconsolidated Zones in Coal Mining Areas. J. Min. Geol. Sci.; 2017; 60, 6.

48. Porzucek, S.; Loj, M. Microgravity Survey to Detect Voids and Loosening Zones in the Vicinity of the Mine Shaft. Energies; 2021; 14, 3021. [DOI: https://dx.doi.org/10.3390/en14113021]

49. Doležalová, H. Suggestion of Complex Monitoring of Undermined Territory. Geophysics in Mining and Environmental Protection; Springer: Berlin/Heidelberg, Germany, 2011; pp. 81-89. ISBN 978-3-642-19096-4

50. Greinert, A. Lignite Mining and Reclamation of Post-Mining Areas in the Lubuskie Region; University of Zielona Góra Press: Zielona Góra, Poland, 2015; ISBN 978-83-937619-2-0 (In Polish)

51. Koźma, J. Analysis of the Landscape Evolution of the Polish Part of the Muskau Arch and Its Valorization in Terms of the Protection of Geological Heritage. Ph.D. Thesis; Polish Geological Institute—National Research Institute: Cracow, Poland, 2018; (In Polish)

52. Rascher, J.; Meier, J.; Kupetz, M. Der Geopark Muskauer Faltenbogen-Grundlagen, Stand, Perspektiven. Geowiss. Mitt. Thüringen. Beih.; 2000; 10, pp. 75-85. Weimar. 698. (In German)

53. Haracz, P.; Iwlew, B.; Jagiełło, K.; Koźma, J.; Maciantowicz, M. The European Muskau Arch. Geopark from Three Different Sides…; “Chroma” Printing House Krzysztof Raczkowski: Łęknica, Poland, 2012; ISBN 978-83-61315-15-5 (In Polish)

54. Fajklewicz, Z. Applied Gravimetry; University Scientific and Educational Publishers of the AGH University of Science and Technology: Cracow, Poland, 2007; ISBN 978-83-7464-032-9 (In Polish)

55. Loke, M. Electrical Imaging Surveys for Environmental and Engineering Studies; A Practical Guide to 2-D and 3-D Surveys. 2000; Available online: https://www.semanticscholar.org/paper/Electrical-Imaging-Surveys-for-Environmental-and-Loke/efa943ffb44f5a217e2d22255f7e334778932248 (accessed on 17 January 2022).

56. Pánek, T.; Hradecký, J.; Šilhán, K. Application of Electrical Resistivity Tomography (ERT) in the Study of Various Types of Slope Deformations in Anisotropic Bedrock: Case Studies from the Flysh Carpathians. Studia Geomorphol. Carpatho-Balc.; 2008; 42, pp. 57-73.

57. Kneisel, C.; Hauck, C. Multi-Method Geophysical Investigation of a Sporadic Permafrost Occurrence. Geophysical Applications in Geomorphology; Gebrüder Borntraeger: Stuttgart, Germany, 2003; pp. 145-159. ISBN 978-3-443-21132-5

58. Farag, K.S.I. Multi-Dimensional Resistivity Models of the Shallow Coal Seams at the Opencast Mine “Garzweiler I” (Northwest of Cologne) Inferred from Radiomagnetotelluric, Transient Electromagnetic and Laboratory Data. Ph.D. Thesis; Univesity of Cologne: Cologne, Germany, 2005.

59. Gontaszewska, A. Lignite Underground Extraction in Ziemia Lubuska—Former Mining, Current Problem. Min. Rev.; 2015; 71, pp. 1-8. (In Polish)

60. Kotyrba, A.; Frolik, A.; Kortas, Ł.; Siwek, S. Post-Mining Hazards in Areas of Former Underground Brown Coal Mines in the Region of Piła-Młyn (Kuyavian-Pomeranian Voivodeship). Min. Rev.; 2012; 68, pp. 58-66. (In Polish)