Search

Abstract Mine operational safety is an important aspect of maintaining the operational continuity of a mining area. In this study, we used the InSAR time series to analyze land surface changes using the ICOPS (improved combined scatterers with optimized point scatters) method. This ICOPS method combines persistent scatterers (PS) with distributed scatterers (DS) to increase surface deformation analysis’s spatial coverage and quality. One of the improvements of this study is the use of machine learning in postprocessing, based on convolutional neural networks, to increase the reliability of results. This study used data from the Sentinel-1 SAR C-band satellite during the 2016–2022 observation period at the Musan mine, North Korea. In the InSAR surface deformation time analysis, the maximum average rate of land subsidence was approximately > 15.00 cm per year, with total surface deformation of 170 cm and 70 cm for the eastern dumping area and the western dumping area, respectively. Analyzing the mechanism of land surface changes also involved evaluating the geological conditions in the Musan mining area. Our research findings show that combining machine learning and statistical methods has great potential to enhance the understanding of mine surface deformation.

We analyzed ground subsidence at the coal ash disposal sites of Stanton Energy Center, a power plant located in Orlando, Florida, USA, by applying 157 Sentinel-1 SAR images obtained between May 2017 and December 2022 in ascending orbit to the PSInSAR technique. A LiDAR DEM with 1 m posting was used for the DInSAR and StaMPS processing for PSInSAR. The results showed significant ground subsidence on the area where solar panels were installed on top of the coal ash landfill. The coal ash landfill was divided into three sites (A, B, and C) according to the landfill sequence. The spatially averaged PSInSAR showed subsidence rates of 7.3 mm/year, 6.2 mm/year, and 8.8 mm/year in sites A, B, and C, respectively. In particular, relatively newly deposited sites A and B showed a decreasing trend in subsidence rate with higher quadratic components in regression function, indicating a stabilization of the subsidence. On the other hand, the oldest site C exhibited the highest (and a relatively constant) subsidence rate, suggesting that the settlement occurred earlier and is now at a constant rate. It is also suspected that new dumping activity near C might have caused a higher subsidence rate than in sites A and B. No subsidence occurred at other solar panel installations on the ground outside the landfill, suggesting that the subsidence was caused by the gravitational compaction of the landfill materials rather than by the instability of the solar facilities. Comparison of PSInSAR results with lower resolution DEMs, such as SRTM and Copernicus DEM, showed range errors of the PS positions proportional to the height deviation from LiDAR DEM, highlighting the importance of accurate DEMs for the time-series analysis of SAR data.

The Okgye limestone mine, which is the largest open-pit limestone mine located in a mountainous area in Korea, suffered a collapse in 2012 that claimed four casualties. Restoration work on the rocky mined-out slopes, as well as mining and dumping activities, are still in progress. Monitoring slope stability is important to prevent the sudden collapse of slopes, which can be efficiently performed by satellite-based interferometric synthetic aperture radar (InSAR) techniques. Firstly, we obtained elevation changes using InSAR-generated Copernicus 30 m DEM in 2014 and an SRTM 1Sec DEM in 2000, through which the area was roughly classified into the mining area, tailings storage area, and the mined-out area. A time series of 12-day coherence images produced by Sentinel-1B SAR were averaged annually to produce an RGB-composite image to observe the change in mining activities during 2018, 2019, and 2020. We found many persistent scatterers (PS) when observing the ground displacement, both in the ascending and descending orbits, from which we decomposed this into the vertical and east components. The largest displacement of 63.6 mm/year was observed during 2019 and 2020 in the tailings storage area in the direction of the dumping slope. For the rocky outcrops and the transmission tower, we found a seasonal oscillation, which can be interpreted as the thermal expansion of limestone and iron. This paper demonstrated that the surface stability and deformation of open-pit mines could be effectively monitored by combining InSAR DEM, coherence, and PSInSAR techniques.

In this study, elevation change and surface morphology of CookE2, one of the most active subglacial lakes in East Antarctica, were analyzed by using Differential Interferometric Synthetic Aperture Radar (DInSAR) and a newly adapted Time-Segmented Persistent Scatterer Interferometric Synthetic Aperture Radar (TS-PSInSAR) techniques. Firstly, several DInSAR pairs were used to study the surface morphology of the subglacial lake during the rapid discharge event in 2007 and the subsequent recharge in 2010 by using ALOS PALSAR data and the continuous recharge from 2018 to 2020 by using Sentinel-1 SAR data. For time-series observation from 2018 to 2020, however, simple integration of DInSAR deviates largely from the satellite altimeter data because errors from the horizontal flow of the surrounding ice field or atmospheric phase accumulate. Conventional PSInSAR deviates from the altimeter data if the LOS displacement exceeds 300 mm, i.e., approximately 1/4 of the slant range resolution of the Sentinel-1 SAR in Interferometric Wide-swath (IW) mode, during the time window. Therefore, a series of Time-Segmented PSInSAR with a 4-month time window could accurately distinguish 1.10 ± 0.01 m/year of highly linear (R2 = 0.99) surface rise rate of CookE2 and 0.63 m/year of horizontal deformation rate of the surrounding ice field from 2018 to 2020.

Ground-based synthetic aperture radar (GB-SAR) is a useful tool to simulate advanced SAR systems with its flexibility on RF system and SAR configuration. This paper reports an indoor experiment of bistatic/multistatic GB-SAR operated in Ku-band with two antennae: one antenna was stationary on the ground and the other was moving along a linear rail. Multiple bistatic GB-SAR images were taken with various stationary antenna positions, and then averaged to simulate a multistatic GB-SAR configuration composed of a moving Tx antenna along a rail and multiple stationary Rx antennae with various viewing angles. This configuration simulates the use of a spaceborne/airborne SAR system as a transmitting antenna and multiple ground-based stationary antennae as receiving antennae to obtain omni-directional scattering images. This SAR geometry with one-stationary and one-moving antennae configuration was analyzed and a time-domain SAR focusing algorithm was adjusted to this geometry. Being stationary for one antenna, the Doppler rate was analyzed to be half of the monostatic case, and the azimuth resolution was doubled. Image quality was enhanced by identifying and reducing azimuth ambiguity. By averaging multiple bistatic images from various stationary antenna positions, a multistatic GB-SAR image was achieved to have better image swath and reduced speckle noise.

Regional changes in the flow velocity of Antarctic glaciers can affect the ice sheet mass balance and formation of surface crevasses. The velocity anomaly of a glacier can be detected using the Double-Differential Interferometric Synthetic Aperture Radar (DDInSAR) technique that removes the constant displacement in two Differential Interferometric SAR (DInSAR) images at different times and shows only the temporally variable displacement. In this study, two circular-shaped ice-velocity anomalies in Campbell Glacier, East Antarctica, were analyzed by using 13 DDInSAR images generated from COSMO-SkyMED one-day tandem DInSAR images in 2010–2011. The topography of the ice surface and ice bed were obtained from the helicopter-borne Ice Penetrating Radar (IPR) surveys in 2016–2017. Denoted as A and B, the velocity anomalies were in circular shapes with radii of ~800 m, located 14.7 km (A) and 11.3 km (B) upstream from the grounding line of the Campbell Glacier. Velocity anomalies were up to ~1 cm/day for A and ~5 cm/day for B. To investigate the cause of the two velocity anomalies, the ice surface and bed profiles derived from the IPR survey crossing the anomalies were analyzed. The two anomalies lay over a bed hill along the glacial valley where stick-slip and pressure melting can occur, resulting in temporal variation of ice velocity. The bright radar reflection and flat hydraulic head at the ice bed of A observed in the IPR-derived radargram strongly suggested the existence of basal water in a form of reservoir or film, which caused smaller friction and the reduced variation of stick-slip motion compared to B. Crevasses began to appear at B due to tensile stress at the top of the hill and the fast flow downstream. The sporadic shift of the location of anomalies suggests complex pressure melting and transportation of the basal water over the bed hill.

In this study, time-series of Sentinel-1A/B Interferometric Synthetic Aperture Radar (InSAR) coherence images were used to monitor the mining activity of Musan open-pit mine, the largest iron mine in North Korea. First, the subtraction of SRTM DEM (2000) from TanDEM-X DEM (2010–2015) has identified two major accumulation areas, one in the east (+112.33 m) and the other in the west (+84.03 m), and a major excavation area (−42.54 m) at the center of the mine. A total of 89 high-quality coherence images with a 12-day baseline from 2015 to 2020 were converted to the normalized difference activity index (NDAI), a newly developed activity indicator robust to spatial and temporal decorrelation. An RGB composite of annually averaged NDAI maps (red for 2019, green for 2018, and blue for 2017) showed that overall activity has diminished since 2018. Dumping slopes were categorized into shrinking, expanding, or transitional, according to the color pattern. Migration and expansion of excavation sites were also found on the pit floor. Time series of 12-day NDAI graphs revealed the date of activities with monthly accuracy. It is believed that NDAI with continuous acquisition of Sentinel-1A/B data can provide detailed monitoring of various types of activities in open-pit mines especially with limited in situ data.

Recent advancement of satellite synthetic aperture radar (SAR) techniques require more sophisticated system configurations such as the use of bistatic antennas or multi-frequencies. A ground-based experiment is a cost-effective and efficient way to evaluate those new configurations especially in the early stage of the system development. In this paper, a ground-based synthetic aperture radar (GB-SAR) system was constructed and operated in a bistatic mode at Ku-band where a receiving antenna (Rx) follows a transmitting antenna (Tx) separated by a baseline B. A new bistatic GB-SAR focusing algorithm was developed by modifying a conventional range-Doppler algorithm (RDA), and its performance has been evaluated by comparing the results with those from a back-projection algorithm (BPA). The results showed good performance of RDA at far range approaching nominal resolutions of 9.4 cm in range and 4.5 cm in azimuth, but limited quality at near range due to the approximation used in RDA. Signals from three trihedral corner reflectors (CR) reduced with increasing B, showing a typical bidirectional scattering behavior of CR. This GB-SAR system will be a testbed for new SAR imaging configurations with variations in antenna positions and target properties.

The international symposium on remote sensing 2018 (ISRS 2018) was held in Pyeongchang, Korea, 9−11 May 2018 [...]

This paper presents the development and experiment of three-dimensional image formation by using a ground-based tomographic synthetic aperture radar (GB-TomoSAR) system. GB-TomoSAR formulates two-dimensional synthetic aperture by the motion of antennae, both in azimuth and vertical directions. After range compression, three-dimensional image focusing is performed by applying Deramp-FFT (Fast Fourier Transform) algorithms, both in azimuth and vertical directions. Geometric and radiometric calibrations were applied to make an image cube, which is then projected into range-azimuth and range-vertical cross-sections for visualization. An experiment with a C-band GB-TomoSAR system with a scan length of 2.49 m and 1.86 m in azimuth and vertical-direction, respectively, shows distinctive three-dimensional radar backscattering of stable buildings and roads with resolutions similar to the theoretical values. Unstable objects such as trees and moving cars generate severe noise due to decorrelation during the eight-hour image-acquisition time.