Monitoring ground subsidence induced by salt mining in the city of Tuzla (Bosnia and Herzegovina)

Salt mining induced ground subsidence is a major hazard in the city of Tuzla (Northeastern sector of Bosnia and Herzegovina) and its surroundings since 1950, when solution mining of salt deposits by boreholes began. An analysis of the large (and never before processed) amount of topographical data collected during two periods: from 1956 to the Balkan War, and from 1992 to 2003 has been made. The analysis reveals a cumulative subsidence as high as 12 m during the whole period, causing damage to buildings and infrastructures within an area that includes a large portion of the historical town. Human-induced subsidence, (with rates up to 40 cm/year in the most developed area), has been investigated to recognize the areas affected by the sinking phenomenon and to produce a subsidence hazard. The time series of topographical observations have been enlarged by conducting new surveys in the urban area by modern space-geodesy methodologies, such as static relative GPS (Global Positioning System) and high resolution satellite imageries. The GPS monitoring started in 2004 and detected a decrease in the subsidence rates to 20 cm/year related to the reduction of salt exploitation. There is close correlation between the average subsidence rate and the annual amount of salt extracted.     

Monitoring of coal mining subsidence in peri-urban area of Zonguldak city (NW Turkey) with persistent scatterer interferometry using ALOS-PALSAR

This paper presents the results of spatiotemporal monitoring of surface subsidence over a mining area in Zonguldak Province of Turkey using Synthetic Aperture Radar (SAR) data, providing maps of subsidence rates in the radar line of sight direction. A total of 18 SAR images, acquired between January 2007 and June 2010 by the Japanese Advanced Land Observing Satellite, have been used to map the surface displacements using the Persistent Scatterer Interferometry technique. The use of Phased Array Type L-band Synthetic Aperture Radar data has proved useful for avoiding signal decorrelation and estimating surface deformation in the heavily vegetated study region. The technique enables the monitoring of continuous small displacements over a large area. Our findings present that many Persistent Scatterers were located on the vegetation cover. The results reveal areas of ground surface subsidence up to 44 mm/year that are well correlated with the underground coal mining galleries particularly in the Gelik region where the Karadon mining galleries are present.     

多基线距DInSAR反演河南省永城煤矿区地表形变

永城是河南省重要的煤炭生产基地,煤矿开采后形成的塌陷区中有70%以上是耕地,已严重威胁到当地居民的生产生活秩序以及生命财产安全,迫切需要采用快速高效的技术手段获取塌陷区的信息来为相关部门提供决策依据;研究中选用2004～2006年的6景ASAR数据,采用基于相干目标的多基线距DInSAR技术对河南永城煤矿区的地表形变进行反演;研究中采用的地表形变反演技术克服了传统DInSAR技术失相干的问题,反演得出的地表形变的范围和程度与调查结果相吻合,表明该项技术具有很强的实用性和推广性;最后分析反演误差产生的原因并提出今后可采用L波段以及多平台联合监测来提高煤矿区地表形变的精度。     

Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain

Subsidence related to multiple natural and human-induced processes affects an increasing number of areas worldwide. Although this phenomenon may involve surface deformation with 3D displacement components, negative vertical movement, either progressive or episodic, tends to dominate. Over the last decades, differential SAR interferometry (DInSAR) has become a very useful remote sensing tool for accurately measuring the spatial and temporal evolution of surface displacements over broad areas. This work discusses the main advantages and limitations of addressing active subsidence phenomena by means of DInSAR techniques from an end-user point of view. Special attention is paid to the spatial and temporal resolution, the precision of the measurements, and the usefulness of the data. The presented analysis is focused on DInSAR results exploitation of various ground subsidence phenomena (groundwater withdrawal, soil compaction, mining subsidence, evaporite dissolution subsidence, and volcanic deformation) with different displacement patterns in a selection of subsidence areas in Spain. Finally, a cost comparative study is performed for the different techniques applied.     

Joint analysis of SAR interferometry and electrical resistivity tomography surveys for investigating ground deformation: the case-study of Satriano di Lucania (Potenza, Italy)

We have jointly applied microwave remote sensing imaging and ground-based geophysical methodologies for investigating ground deformation. In particular, the Differential Interferometric Synthetic Aperture Radar (DInSAR) technique referred to as Small BAseline Subset (SBAS) algorithm and the Electrical Resistivity Tomography (ERT) method have been applied. The former is a tool for investigating large areas, allowing us to detect deformation phenomena and to analyse their spatial patterns and temporal evolution; the latter permits us to reconstruct, with high spatial resolution, the geometry of near-surface geological structures and to estimate the electrical resistivity of outcropping terrains. The starting point of our analysis has been a DInSAR survey carried out on an area of about 400 km² centred around the town of Potenza (Basilicata region, Italy), for which SAR data acquired by the ERS-1/2 sensors in the time interval 1992–2000 were available. The DInSAR analysis allowed us to detect an unknown ground deformation effect which involves the urban area of Satriano di Lucania, located close to Potenza. An on-site ERT survey was then carried out in order to characterize in detail the physical properties of the terrains involved in the main deformation pattern detected via the DInSAR processing. The integration of the DInSAR and ERT measurements permitted us to investigate areas where the deformation have space-variant characteristics and allowed us to formulate a hypothesis on the origin of the detected displacement processes that, although not fully conclusive, is consistent with the overall DInSAR and ERT analysis.     

高应力区露天转地下开采岩体移动规律

以金川龙首矿为例,基于2003―2008年长达12期的GPS实测数据,按照时空分布比较分析总体变形趋势,得出整个矿区的地表沉降特征、平面位移特征、行线（勘探线）随时间变化的沉降和平面位移特征。据此分析了以水平构造应力为主导的高应力区露天转地下开采引起岩体的变形规律,并与该矿闭坑之初的变形规律比较,通过理论分析的方法对产生差异的原因进行初步解释,并用数值模拟的方法加以验证,其研究结果对认识露天转地下开采引起的岩体移动、破坏特征,指导矿区安全生产有积极意义     

蔚县矿区地面沉陷InSAR多维形变监测

河北蔚县是我国华北地区最大的地下采煤区之一,该地区长期存在采矿塌陷灾害,不仅威胁采矿安全,而且严重破坏当地生态环境。本文基于合成孔径雷达干涉（InSAR）这一新型空间对地观测技术,采用61景Sentinel-1A/B干涉宽幅（Interferometric Wide swath,IW）模式数据进行矿区形变观测,获取整个矿区在2017—2018年间的地表形变空间分布特征,并对矿区地表的沉降量级及面积进行详细的统计分析。此外,采用融合多轨道SAR数据的多维形变时序估计方法,对西细庄矿数据进行东西向和垂向的二维形变分解,获取该矿的二维形变时间序列。结果表明:除南留庄井田外,其余三大井田在监测期间均存在不同程度的地面沉陷灾害;整个矿区年沉陷速率超过-10 cm/a的区域达到了2.16 km2;受成像几何影响,不同轨道数据获取的形变结果存在一定差异;西细庄矿以垂向形变为主,伴随明显的东西向水平形变。研究结果为蔚县矿区地面沉陷监测与煤矿安全开采提供数据参考。      

Using differential SAR interferometry to map land subsidence: a case study in the Pingtung Plain of SW Taiwan

Synthetic aperture radar (SAR) interferometry (InSAR) is a geodetic tool widely applied in the studies of earth-surface deformation. This technique has the benefits of high spatial resolution and centimetre-scale accuracy. Differential SAR interferometry (DInSAR) is used to measure ground deformation with repeat-pass SAR images. This study applied DInSAR and persistent scatterers InSAR (PSInSAR) for detecting land subsidence in the Pingtung Plain, southern Taiwan, between 1995 and 2000. In recent years, serious land subsidence occurred along coastal regions of Taiwan as a consequence of over-pumping of underground water. Results of this study revealed that the critical subsidence region is located on the coast near the estuary of Linpien River. It is also found that subsidence was significantly higher during the dry season than the wet season. The maximum annual subsidence rate of the dry season is up to −11.51 cm/year in critical subsidence region and the vertical land movement rate is much slower during the wet season. The average subsidence rates in wet and dry seasons are −0.31 and −3.37 cm/year, respectively. As a result, the subsidence rate in dry seasons is about 3 cm larger than in wet seasons.     

一矿山开采崩落区的岩体变形及崩落角确定

矿山地下开采过程中,由于采出了大量的矿体,上部岩体失去支撑,破坏了原有的平衡,因此岩体就破裂、变形,地表上则表现为塌陷、裂缝以及肉眼难以分辨的地面变形,矿坑内则表现为片帮、冒顶、底鼓等。     

基于InSAR技术的缓倾煤层开采诱发顺层岩体地表变形模式研究

贵州贞丰县某煤矿开采煤层以向斜缓倾的三叠系上统火把冲组（T3h）为主,与贵州省大部分煤矿开采的背斜反倾煤层不同,其采矿活动诱发的地面沉降和滑坡风险亦表现出不同的变形破坏模式(背斜反倾煤层易诱发倾倒崩塌、顺层缓倾煤层易诱发地面塌陷与滑坡)。论文利用升、降轨观测的共15期3 m空间分辨率L波段PALSAR -2 SAR为数据源,开展了多期地表变形D -InSAR测量,确定出变形发生的位置、范围与滞后时间。经实地调查验证,InSAR解算结果较好地吻合了矿区开采范围和地表破坏情况,证实了InSAR在煤矿区识别时序性地表形变的准确性。进而分解计算了地表三维变形,并通过与地下开采范围和过程的相关性分析,深化了对该地区缓倾煤层地下开采诱发的顺层滑坡变形模式的认识:（1）InSAR可以识别计算出采矿区地表变形的范围与沉降量,矿区变形在干涉影像中表现为以采空区地表为中心向四周扩散的圆环状变形条纹;（2）地表变形区域覆盖地下采空区上方及附近地表区域,根据地表变形情况与地下采空区范围计算出该地区上山边界角约70°、下山边界角约58°;（3）地下采空与地表沉降变形存在约30 d的时间滞后;（4）顺层地下采空引发的地表水平移动方向受地层产状、地表坡向共同作用,水平向为沿层面的顺层滑移与向沉降中心汇聚的合成运动结果;（5）沿层面的顺层滑移与地表坡度因素叠加造成采空区地表上山侧岩石受拉产生拉裂缝,下山侧则易产生塌陷坑及裂缝。     

城市连拱隧道施工地表沉降分析及预测研究

地表沉降的研究工作大都集中在地铁隧道盾构法施工的情况,而对于矿山法城市双连拱隧道施工引起地表沉降的研究几乎没有。本文以重庆彭家花园双连拱隧道为工程背景,就矿山法双连拱隧道施工引起的地表变形进行研究。在不同施工方法条件下,针对不同埋深、断面大小进行回归分析,从理论上弄清地表沉降的规律,根据实测资料,提出矿山法施工条件下地表沉降的预测公式。并结合三维数值模拟验证了回归效果     

Multi-sensor DInSAR applied to the spatiotemporal evolution analysis of ground surface deformation in Cerro Prieto basin,Baja California,Mexico, for the 1993–2014 period

The combined effects of active tectonics and anthropogenic activities, primarily geothermal resources exploitation for electricity production in Cerro Prieto geothermal field, influence the ground surface deformation in Cerro Prieto basin, Baja California, Mexico. In this study, a large set of multi-sensor C-band SAR images have been employed to reconstruct the spatiotemporal evolution of aseismic ground surface deformation that has affected Cerro Prieto basin from 1993 to 2014. Conventional DInSAR together with the interferograms stacking procedure was applied. The results showed that the study area presented considerable surface deformation (mainly subsidence) during the entire time of the investigation. The main changes in rate and pattern of surface deformation have a good correlation in time and space with the changes in production in the Cerro Prieto geothermal field. Comparison of LOS displacement maps from different viewing geometries, and decomposition (where possible) of LOS displacement into vertical and horizontal (east–west) components, revealed considerable horizontal displacement which mostly reflects the ground movement at and beyond the margin of the subsidence basin toward the areas of highest subsidence rates. In addition, the validation of the DInSAR results by comparing them against measurements from leveling surveys was performed, confirming the high reliably of satellite interferometry for the ground surface deformation rate mapping in the study area.     

Vulnerability to ground deformation phenomena in the city of Tuzla (BiH): a GIS and multicriteria approach

Ground deformation phenomena affecting the town of Tuzla (Bosnia Herzegovina), as consequence of underground salt deposits exploitation, have been surveyed since 60?years ago. The growing of cavities beneath the town and successive surface deformations have induced several hazard factors such as ground subsidence, deep and shallow fracturing and intense groundwater table oscillations. Fortunately, this processes had no catastrophic consequences such as sinkholes developing. Previous studies analyzed the historical database of topographic and piezometric data and highlighted a cumulative subsidence rate up to 12 meters over a period spanning from the 1956 to the 2003. Fractures arose as obvious result of ground deformation and caused damages and demolitions to thousands of buildings, with almost 15,000 people evacuated. Nowadays, in order to diminish the deformation processes, a strong reduction in brine withdrawal has been introduced by the local authorities. This caused an uplift of the water table which actually makes the spatial pattern of deformation more complex. In fact, in the most urbanized area, the sinking phenomena are now turning into uplift displacements. The presence of such a hazard required an urban vulnerability assessment able to support local authorities in the forthcoming planning procedure. The vulnerability assessment of urban elements over a wide area, potentially affected by deformation phenomena, has been accomplished by an integration between GIS solutions and a multicriteria approach based on ELECTRE-TRI methodology. The vulnerability procedure is able to define a discrete map where the 50?×?50?m width elementary cells report the most vulnerable parts of the town. In particular, the ELECTRE-TRI methodology is used to combine each vulnerable urban element with the others, expressed as vector data in a geographical dataset. The multicriteria approach is performed outside the GIS environment by the ELECTRE-TRI 2.0 package and requires a preliminary rasterization of the used dataset.     

基于InSAR三维分解技术的贵州省贞丰某煤矿地下采空区反演

贵州是我国主要煤炭产地之一, 成本与地方条件的限制使贵州贞丰县山区出现大量的非法小煤窑, 其大多开采一年、半年便废弃煤矿, 调查其地下开采范围对矿区生态修复、土地资源再利用、煤矿越界开采监测等有一定意义。本文利用升、降轨观测共16期3 m空间分辨率的L波段PALSAR-2影像为数据源, 在贵州省贞丰县采用InSAR技术对煤矿采空区开展了短期动态地表沉降监测, 结合研究区地质条件、地下开采条件、开采时间、范围等信息, 计算了该地区煤矿开采沉陷影响角的规律与变形时间滞后关系; 而后利用上述参数及地表变形信息, 在其附近一处废弃煤矿, 对其地下采空范围及开采时间进行反演:(1)研究区煤炭开采走向影响角为83°、上山影响角75°、下山影响角80°, 地表变形时间与地下开采时间差约1个月; (2)计算反演出研究区废弃煤矿地下开采范围为380 m×150 m; (3)在反演煤矿采区, 对比FLAC3D数值模拟计算结果和野外现场调查, 验证了InSAR计算结果基本准确, 验证了此方法对煤矿采区反演的合理性与可靠性。     

Prediction of ground subsidence in Samcheok City,Korea using artificial neural networks and GIS

This study shows the construction of a hazard map for presumptive ground subsidence around abandoned underground coal mines (AUCMs) at Samcheok City in Korea using an artificial neural network, with a geographic information system (GIS). To evaluate the factors governing ground subsidence, an image database was constructed from a topographical map, geological map, mining tunnel map, global positioning system (GPS) data, land use map, digital elevation model (DEM) data, and borehole data. An attribute database was also constructed by employing field investigations and reinforcement working reports for the existing ground subsidence areas at the study site. Seven major factors controlling ground subsidence were determined from the probability analysis of the existing ground subsidence area. Depth of drift from the mining tunnel map, DEM and slope gradient obtained from the topographical map, groundwater level and permeability from borehole data, geology and land use. These factors were employed by with artificial neural networks to analyze ground subsidence hazard. Each factor’s weight was determined by the back-propagation training method. Then the ground subsidence hazard indices were calculated using the trained back-propagation weights, and the ground subsidence hazard map was created by GIS. Ground subsidence locations were used to verify results of the ground subsidence hazard map and the verification results showed 96.06% accuracy. The verification results exhibited sufficient agreement between the presumptive hazard map and the existing data on ground subsidence area. An erratum to this article can be found at     

PS技术及其在地表形变监测中的应用现状与发展

合成孔径雷达干涉测量技术（InSAR）为地表形变监测提供了极具应用潜力的手段,具有大面积、高空间分辨率、全天候及成本低的优点。但由于大气条件变化、地表覆被等时间空间去相干的影响,其精度和普适性受到极大限制。 近几年发展的永久散射体（PS）技术在传统差分干涉测量（DInSAR）中引入时间维,分析长时间内保持稳定的像元集相位变化,获得毫米级的地表形变测量精度,同时有效地解决了时间空间去相关和大气非均质性影响的问题,目前在滑坡、地面沉降和地质灾害监测等领域得到了广泛的应用。PS技术具有高精度、高时间分辨率、能极大提高影像利用率的优点;但只适用较小区域、需要大量影像、且不适于分析快速突变的地表形变。为克服PS应用中的问题,近年来出现了三角反射体技术、多平台PS技术及相关性像元分析（CPT）技术,使PS技术应用具有更广泛的适应性。     

露天与地下联合开采引起岩层移动规律的模型试验研究

非法开采对矿山开采稳定性具有重要影响,以某铁矿为例,进行了二维物理相似模型试验研究,分析得出了在考虑和不考虑非法开采两种情况下的地表变形规律和矿体围岩移动规律。试验结果表明,从-270 m采至-550 m水平期间,地表沉降和围岩变形呈缓慢增加趋势,但当从-550 m水平继续向下开采时,由于矿体倾角变得平缓,地表沉降、深部围岩变形都急剧增加,非法开采会造成矿体上部岩体的应力重分布,虽没有明显增大地表的沉降范围,但增加了开采点附近的地表沉降量,而且大大增大了地表岩层的水平移动。     

复杂地质条件下矿山地下开采地表变形规律的研究

地质条件是矿山开采岩层移动的基础,复杂的地质条件常使得地表移动角、沉陷角预测不准,导致灾难性后果。以典型复杂地质矿山程潮铁矿东区为例,以多年地表变形监测结果为依据,对在其特殊地质条件下地表变形和岩层移动的机制进行了分析研究。结果表明,地下采矿是矿区岩体破坏的直接诱因,地下水产生的静水压力和动水压力强化了岩体的变形与破坏;矿区岩体中断层、节理等特殊地质结构是岩体破坏的基础;构造应力场的存在是围岩向采空区产生较大移动变形的主要影响因素;在采矿之初,地表塌陷主要是由地下水疏干引起的;地表大规模塌陷形成以后,地下采空区的扩大是引起地表塌陷的主要原因,但矿区特殊的地质条件对地表塌陷范围的扩展速度有重要影响。     

Crustal deformation in the Central Ionian Islands (Greece): Results from DGPS and DInSAR analyses (1995–2006)

Ground deformation studies based on Differential GPS (DGPS) and Differential Interferometric SAR (DInSAR) analyses have been conducted in the seismically active area of the Central Ionian Islands. Local GPS networks were installed in Cephallonia (2001) and Zakynthos (2005). The Cephallonian network has been remeasured five times and Zakynthos' once as of July 2006. The studies have yielded detailed information regarding both local and regional deformations that are occurring in the area.For Lefkas Island, DInSAR analysis (March to September 2003) revealed 56 mm of uplift in the central and western parts and is attributed to the August 2003 earthquake (M_w = 6.3) that occurred offshore to the west. Synthetic DInSAR modelling of the magnitude and extent of deformation is consistent with the seismologically deduced parameters for the ruptured segment along the Lefkas Transform Fault. Subsidence (< 28 mm) along the northern part of the island is attributed to local conditions unrelated to the earthquake. For Zakynthos Island, large-magnitude earthquakes that occurred offshore to the south in October 2005 and April 2006 most likely contributed to the observed deformation as deduced from DGPS measurements for an encompassing period (August 2005 to July 2006). The largest amount of horizontal deformation occurred in the south, where its western part moved in a W–NW direction, while the eastern part moved towards the NE, with magnitudes ranging from 15 to 26 mm. The southern part of the island uplifted a maximum of 65 mm whereas the north subsided from 12 to 28 mm.For Cephallonia Island, DInSAR analysis (1995 to 1998) indicated ground deformation up to 28 mm located in small sections of the island. Further interferometric analysis for the period 2003 to 2004, encompassing the occurrence of the Lefkas earthquake in August 2003, indicated 28 mm of uplift in the northern part, while during the next two years (2004 to 2005), further uplift of at least 56 mm had taken place at the western and northern part of the island.DGPS measurements for the period 2001 to 2006 revealed a clockwise rotation of the island with respect to a centrally located station on Aenos Mt. The horizontal component of deformation generally ranged from 6 to 34 mm, with the largest values at the western and northern parts of the island. Considering the vertical deformation, two periods are distinguished. The first one (2001 to 2003) is consistent with anticipated motions associated with the main geological and tectonic features of the island. The second one (2003 to 2006) has been tentatively attributed to dilatancy in which relatively small uplift (20–40 mm) occurred along the southern and southeastern parts of the island, while larger magnitudes (> 50 mm) happened at the western part (Paliki Peninsula). These large magnitudes of uplift over an extended area (> 50 km), in conjunction with an accelerated Benioff strain determined from the analysis of the seismicity in the broader region, are consistent with dilatancy. This effect commenced some time after 2003 and is probably centered in the area between Zakynthos and Cephallonia. If this interpretation is correct, it may foreshadow the occurrence of a very strong earthquake(s) sometime during 2007 to 2008 in the above designated region.     

太原市地面沉降的现状及发展变化趋势

地面沉降作为一种广泛的地质灾害,不仅对地面设施造成很大的破坏,而且也会使地下环境发生变化。从太原市地下水开采历史与地面沉降演变历史总体上相符可以看出,太原市地面沉降是由于过量开采中深层孔隙地下水引起地下水位急剧下降造成的。随着2003年太原市采取关井压采及引黄入并等一系列举措以来,地下水位呈现回升的趋势,地面沉降速率随之减缓。