Mapping land subsidence over the eastern Beijing city using satellite radar interferometry*

Beijing City has suffered from groundwater-induced subsidence since the late 1930s and the over-exploration of groundwater could lead to subsidence as much as ?12.0 cm?yr^?1. Previous studies on the ground deformation at Beijing City mainly focused on the period before the year of 2014 when a mega-engineering project was launched to reduce water shortage in Beijing. To study the most recent ground deformation, 19 L-band ALOS-1 PALSAR images (June 2007–January 2011), 24 C-band Sentinel-1 SAR images (June 2015–November 2016) together with 9 ALOS-2 PALSAR acquisitions (September 2014–February 2017) were analysed in this work. Levelling measurements were exploited to verify the ALOS-1-based time series InSAR (TS-InSAR) result while Sentinel-1 and ALOS-2 result were cross-verified with each other. Furthermore, the whole study area was divided into four sub-zones, and the result indicated that the subsidence rates over five townships, Cuigezhuan, Jinzhan, Liyuan, Songzhuang and Yanjiao were accelerating and more attentions should be paid. On the contrary, the town centre of Douge Zhuang township experienced a decreasing trend between these two temporal-periods. Additionally, the time series measurements with respect to five selected measurement points and the profile line along the subsidence hot spots were analysed.     

宽幅InSAR大地测量学与大尺度形变监测方法

随着人们对大尺度地形信息与地表环境变化监测需求的提升,以及哨兵1号（Sentinel-1A）、大地2号（ALOS-2）等合成孔径雷达（synthetic aperture radar,SAR）卫星的宽幅模式数据不断获取,宽幅雷达干涉测量（SAR interferometry,InSAR）技术已成为大尺度地形测绘、地球动力学（地震、火山、滑坡等）与人工地物结构健康监测等领域的研究热点。分析了两类宽幅SAR数据,即扫描（ScanningSAR,ScanSAR）模式与逐行扫描地形观测（terrain observation by progressive scans SAR,TOPSAR）模式开展干涉测量的主要限制条件与解决方法,探讨了宽幅InSAR形变监测关键误差估计与改正方法、时间序列分析技术与方位向位移观测技术,并给出2008年矩震级Mw 7.1新疆于田地震同震、震后形变监测应用。随着宽幅SAR数据的不断积累,宽幅InSAR大地测量学有望得到深入发展与应用。     

Monitoring active open-pit mine stability in the Rhenish coalfields of Germany using a coherence-based SBAS method

With the recent progress in synthetic aperture radar (SAR) technology, especially the new generation of SAR satellites (Sentinel-1 and TerraSAR-X), our ability to assess slope stability in open-pit mines has significantly improved. The main objective of this work is to map ground displacement and slope instability over three open-pit mines, namely, Hambach, Garzweiler and Inden, in the Rhenish coalfields of Germany to provide long-term monitoring solutions for open-pit mining operations and their surroundings. Three SAR datasets, including Sentinel-1A data in ascending and descending orbits and TerraSAR-X data in a descending orbit, were processed by a modified small baseline subset (SBAS) algorithm, called coherence-based SBAS, to retrieve ground displacement related to the three open-pit mines and their surroundings. Despite the continuously changing topography over these active open-pit mines, the small perpendicular baselines of both Sentinel-1A and TerraSAR-X data were not affected by DEM errors and hence could yield accurate estimates of surface displacement. Significant land subsidence was observed over reclaimed areas, with rates exceeding 500 mm/yr, 380 mm/yr, and 310 mm/yr for the Hambach, Garzweiler and Inden mine, respectively. The compaction process of waste materials is the main contributor to land subsidence. Land uplift was found over the areas near the active working parts of the mines, which was probably due to excavation activities. Horizontal displacement retrieved from the combination of ascending and descending data was analysed, revealing an eastward movement with a maximum rate of ∼120 mm/yr on the western flank and a westward movement with a maximum rate of ∼ 60 mm/yr on the eastern flank of the pit. Former open-pit mines Fortuna-Garsdorf and Berghein in the eastern part of Rhenish coalfields, already reclaimed for agriculture, also show subsidence, at locations reaching 150 mm/yr. The interferometric results were compared, whenever possible, with groundwater information to analyse the possible reasons for ground deformation over the mines and their surroundings.     

InSAR to support sustainable urbanization over compacting aquifers: The case of Toluca Valley,Mexico

This paper illustrates how InSAR alone can be used to delineate potential ground fractures related to aquifer system compaction. An InSAR-derived ground fracturing map of the Toluca Valley, Mexico, is produced and validated through a field campaign. The results are of great interest to support sustainable urbanization and show that InSAR processing of open-access Synthetic Aperture Radar (SAR) data from the Sentinel-1 satellites can lead to reliable and cost-effective products directly usable by cities to help decision-making.The Toluca Valley Aquifer (TVA) sustains the water needs of two million inhabitants living within the valley, a growing industry, an intensively irrigated agricultural area, and 38% of the water needs of the megalopolis of Mexico City, located 40 km east of the valley. Ensuring water sustainability, infrastructure integrity, along with supporting the important economic and demographic growth of the region, is a major challenge for water managers and urban developers. This paper presents a long-term analysis of ground fracturing by interpreting 13 years of InSAR-derived ground displacement measurements. Small Baseline Subset (SBAS) and Persistent Scatterer Interferometry (PSI) techniques are applied over three SAR datasets totalling 93 acquisitions from Envisat, Radarsat-2, and Sentinel-1A satellites and covering the period from 2003 to 2016.From 2003 to 2016, groundwater level declines of up to 1.6 m/yr, land subsidence up to 77 mm/yr, and major infrastructure damages are observed. Groundwater level data show highly variable seasonal responses according to their connectivity to recharge areas. However, the trend of groundwater levels consistently range from −0.5 to −1.5 m/yr regardless of the well location and depth. By analysing the horizontal gradients of vertical land subsidence, we provide a potential ground fracture map to assist in future urban development planning in the Toluca Valley.     

利用SBAS时序技术识别地热资源开采引起的地表形变

本文以地热资源丰富的雄安新区为研究对象,利用SBAS技术对2015年7月至2019年2月的42幅Sentinel-1区域卫星影像和2016年8月至2018年12月的10幅ALOS-2区域卫星影像进行处理,以此获取雄安新区的时间序列地表形变。研究结果表明,雄安新区整体地表沉降严重,且沉降区主要集中在雄县和容城县,其中雄县沉降最为严重,沉降中心累计沉降量超过330 mm,容城县沉降中心累计下沉量超过120 mm。结合地热井位置和地热开采历史,由时间序列分析表明,大量开采地热资源会引起地表沉降现象。     

沧州地区地面沉降现状Sentinel-1A/B时序InSAR监测与分析

获取了22景IW模式的Sentinel-1A/B数据,分别采用PS-InSAR和SBAS-InSAR两种时序技术得到了沧州地区2015年11月至2016年11月期间地表平均沉降速率及累计沉降量。经分析,两种时序监测结果具有较高的一致性,其结果表明：献县、沧县、沧州市区地面表现为回升趋势;青县、沧县东部、南皮县西部、东光县等地区地面沉降量较大,其中东光县何庄村累计下沉量达47 mm,而引起这些区域地面沉降的主要原因为农业生产等带来的地下水开采。本文对基于Sentinel-1A/B SAR数据的时序分析方法在城市沉降监测方面的应用进行了探讨,为有关部门提供了重要的参考数据。     

Mapping wetland characteristics using temporally dense Sentinel-1 and Sentinel-2 data: A case study in the St. Lucia wetlands,South Africa

Wetlands have been determined as one of the most valuable ecosystems on Earth and are currently being lost at alarming rates. Large-scale monitoring of wetlands is of high importance, but also challenging. The Sentinel-1 and -2 satellite missions for the first time provide radar and optical data at high spatial and temporal detail, and with this a unique opportunity for more accurate wetland mapping from space arises. Recent studies already used Sentinel-1 and -2 data to map specific wetland types or characteristics, but for comprehensive wetland characterisations the potential of the data has not been researched yet. The aim of our research was to study the use of the high-resolution and temporally dense Sentinel-1 and -2 data for wetland mapping in multiple levels of characterisation. The use of the data was assessed by applying Random Forests for multiple classification levels including general wetland delineation, wetland vegetation types and surface water dynamics. The results for the St. Lucia wetlands in South Africa showed that combining Sentinel-1 and -2 led to significantly higher classification accuracies than for using the systems separately. Accuracies were relatively poor for classifications in high-vegetated wetlands, as subcanopy flooding could not be detected with Sentinel-1’s C-band sensors operating in VV/VH mode. When excluding high-vegetated areas, overall accuracies were reached of 88.5% for general wetland delineation, 90.7% for mapping wetland vegetation types and 87.1% for mapping surface water dynamics. Sentinel-2 was particularly of value for general wetland delineation, while Sentinel-1 showed more value for mapping wetland vegetation types. Overlaid maps of all classification levels obtained overall accuracies of 69.1% and 76.4% for classifying ten and seven wetland classes respectively.     

哨兵-2号光学影像地表形变监测:以2016年Mw7.8新西兰凯库拉地震为例

光学影像已经广泛地应用于地表形变监测研究。哨兵-2号光学影像作为一种新的对地观测数据,具有重访周期短、空间分辨率高、影像覆盖范围大以及数据免费等优点,因此该数据在地表形变监测上有广泛的应用潜力。本文以COSI-Corr软件包为数据处理平台,基于亚像素的频率域相关性匹配技术,处理多时相的哨兵-2号数据获取地表形变。本文选取2016年M_w7.8新西兰凯库拉（Kaikoura）地震覆盖区域为例,对哨兵-2号影像地表形变中存在的各种系统误差源进行了系统分析和改正,并提出了改进的均值相减法去除形变场中的卫星姿态角误差。另外,还对哨兵-2号4个10 m空间分辨率的波段（Band 2/3/4/8）中可用于地表形变监测的最佳波段进行了分析,统计结果显示Band8的地表形变监测效果最好。最后,利用哨兵-2号光学影像获取了2016年11月14日M_w7.8新西兰凯库拉地震的同震形变场;分析了沿地震主要断层的滑移分布,结果表明最大水平滑移量达10 m;并与同期Landsat8全色影像的同震形变监测结果进行对比分析,结果表明哨兵-2号结果精度更高。本文的研究成果可以为哨兵-2号光学影像的应用提供参考。     

卫星高度计定标现状

定标是卫星高度计数据精度的重要保障,随着高度计卫星HY-2A的发射及后续卫星组网规划,中国将获取长时间序列的自主高度计观测资料,定标对数据精度和长期一致性的重要性日益凸显。总结了卫星高度计定标常用技术方法的国内外研究进展和现状,阐述了中国自主海上定标场的建设和应用情况,重点对青岛千里岩定标场的大地基准测量、地壳沉降监测及HY-2A等多颗卫星高度计的定标结果进行总结分析,并对规划建设中的珠海万山海上综合定标场和中国沿海定标场网做了介绍。此外利用GNSS水汽反演技术对星载微波辐射计观测的大气湿延迟开展了精度检核实验,得到了Jason-2卫星2010年—2016年微波辐射计大气湿延迟观测精度,证明了利用中国沿海GNSS连续运行站标定星载微波辐射计大气湿延迟的可行性,对于充分了解和认识卫星高度计定标的研究现状和发展趋势有一定的借鉴和指导意义。     

星载加速度计增强北斗自主定轨性能

卫星自主定轨是提高全球卫星导航系统(GNSS)可靠性、稳健性、完整性和生存能力的重要保证。新一代的北斗卫星已可以进行星间链路测距,从而达到提高卫星全球跟踪能力以及实现整个卫星导航系统的自主定轨。然而由于卫星运行会受到多种摄动力的影响,如果不能对这些摄动力进行精密的改正,在没有地面或其他天体提供绝对约束的条件下,导航系统会随着自主定轨时间的延长出现星座整体旋转。卫星所受摄动力分为保守力和非保守力两部分:对于保守力,如地球非球形摄动、潮汐摄动、太阳月球和其他三体引力,现在已有的力学模型可以很精确地进行改正;而非保守力(如太阳光压摄动),则难以用精确的模型进行改正,因此成为影响卫星定轨精度的主要因素。星载加速度计可以高精度地测量非保守力,并已成功应用于重力卫星(CHAMP、GRACE、GOCE)的重力场反演与大气研究中。本文研究主要探讨采用星上加速度计提高北斗卫星自主定轨精度和延长自主定轨时长的可行性。利用模拟的卫星轨道和星间链路数据,以及现有的星载加速度计误差模型,对北斗卫星系统分别使用星间链路数据和星间链路与加速度计组合数据,进行自主定轨与精度评定。计算结果表明,使用星间链路与星载加速度计数据进行自主定轨,较单纯使用星间链路数据精度具有明显改进。在模拟的星间测距观测数据具有0.33m随机噪声以及分米级系统误差,自主定轨两个月的情况下,联合使用加速度计数据的自主定轨IGSO和MEO卫星精度为分米级,而仅使用星间链路数据的定轨精度约为3~6m,比使用加速度计精度低一个量级。     

Coseismic surface deformation of the 2014 Napa earthquake mapped by Sentinel-1A SAR and accuracy assessment with COSMO-SkyMed and GPS data as cross validation

ABSTRACT

The new land observation satellite Sentinel-1A was launched on 25 April 2014 with a C-band synthetic aperture radar (SAR) sensor, which has the significant enhancements in terms of revisit period and high resolution. The Mw 6.1 Napa, California earthquake occurring on 24 August 2014, almost 4 months after the launch, is the first moderate earthquake imaged by the Sentinel-1A. This provides an opportunity to map the coseismic deformation of the event and evaluate the potential of Sentinel-1A SAR for earthquake study. Two techniques including the interferometric SAR (InSAR) and pixel offset-tracking (PO) are, respectively, employed to map the surface deformation along the radar line of sight (LOS), azimuth and slant-range directions. The cross comparison between Sentinel-1A InSAR LOS deformation and GPS observations indicates good agreement with an accuracy of ～2.6?mm. We further estimate the earthquake source model with the external COSMO-SkyMed InSAR and GPS data as constraints, and forward calculate the surface deformation as cross validation with the Sentinel-1A observations. The comparison between the observed and modeled deformation shows that the Sentinel-1A measurement accuracy can achieve 1.6?cm for InSAR technique along LOS direction, and 6.3 and 6.7?cm for PO along azimuth and range directions, respectively.     

High-resolution mapping of forest canopy height using machine learning by coupling ICESat-2 LiDAR with Sentinel-1, Sentinel-2 and Landsat-8 data

Forest canopy height is an important indicator of forest carbon storage, productivity, and biodiversity. The present study showed the first attempt to develop a machine-learning workflow to map the spatial pattern of the forest canopy height in a mountainous region in the northeast China by coupling the recently available canopy height (H_canopy) footprint product from ICESat-2 with the Sentinel-1 and Sentinel-2 satellite data. The ICESat-2 H_canopy was initially validated by the high-resolution canopy height from airborne LiDAR data at different spatial scales. Performance comparisons were conducted between two machine-learning models – deep learning (DL) model and random forest (RF) model, and between the Sentinel and Landsat-8 satellites. Results showed that the ICESat-2 H_canopy showed the highest correlation with the airborne LiDAR canopy height at a spatial scale of 250 m with a Pearson’s correlation coefficient (R) of 0.82 and a mean bias of -1.46 m, providing important evidence on the reliability of the ICESat-2 vegetation height product from the case in China’s forest. Both DL and RF models obtained satisfactory accuracy on the upscaling of ICESat-2 H_canopy assisted by Sentinel satellite co-variables with an R-value between the observed and predicted H_canopy equalling 0.78 and 0.68, respectively. Compared to Sentinel satellites, Landsat-8 showed relatively weaker performance in H_canopy prediction, suggesting that the addition of the backscattering coefficients from Sentinel-1 and the red-edge related variables from Sentinel-2 could positively contribute to the prediction of forest canopy height. To our knowledge, few studies have demonstrated large-scale vegetation height mapping in a resolution ≤ 250 m based on the newly available satellites (ICESat-2, Sentinel-1 and Sentinel-2) and DL regression model, particularly in the forest areas in China. Thus, the present work provided a timely and important supplementary to the applications of these new earth observation tools.     

基于少量控制点的Radarsat-2影像快速几何纠正技术研究

Radarsat-2卫星依据其搭载的GPS接收机可实现3倍中误差小于60m的精确实时定轨。由此本文提出依据其影像元数据信息实现快速几何纠正的方法,该方法利用少量的几个控制点来消除Radarsat-2影像与待纠正参考系间的系统误差。从而实现Radarsat-2影像的快速几何纠正。本文并依据SAR斜距成像原理的纠正公式和实地采集的GPS数据,验证了元数据中所提供RPF模型的内部精度和外部符合精度。通过实验验证了本快速纠正技术可以获得中误差小于2个像素的平面几何纠正精度。     

HJ-1光学卫星遥感应用前景分析

HJ-1星座光学卫星作为我国环境与灾害监测预报小卫星星座的首发卫星,其运行与应用状况对于后续星座的发展具有重要 的理论和实践指导意义。本研究从HJ-1光学卫星的分辨率（时间、空间和光谱）、光谱及幅宽等数据特征入手,综合考虑遥感应用 数据要求,全面分析评价HJ-1光学卫星数据可用性。研究表明,两颗光学卫星可基本满足我国及周边国家环境监测与灾害管理方面 的主要应用需求,并在其它领域也具有广阔的应用前景。     

多源星载SAR地形干涉测量精度分析

星载InSAR技术具有高效率、高精度获取全球DEM数据及其他增值产品的优势,是地形测绘领域的研究热点。本文综合利用国内外多源星载SAR影像数据,开展地形干涉测量试验及精度对比分析,旨在为全球测绘提供技术参考。现有的X/C/L波段卫星SAR系统,如COSMO-SkyMed、GF-3、ALOS-2获取的干涉数据集,在青海典型复杂地形实验区均成功获取了DEM数据产品。分析结果表明,在3种典型数据源中,基于COSMO-SkyMed干涉测量DEM的精度与细节质量相对较高,GF-3干涉结果次之,ALOS-2数据也实现了较好的地形测图精度。相关结果从侧面论证了国产GF-3数据具备在空间基线合适的条件下获得高精度DEM数据产品的潜力。     

Water body mapping method with HJ-1A/B satellite imagery

This paper proposes an integrated water body mapping method with HJ-1A/B satellite imagery, the CCD (charge coupled device) data of the Chinese environmental satellites that were launched on September 6th, 2008. It combines the difference between NDVI and NDWI (NDVI–NDWI) with SLOPE and near-infrared (NIR) band. The NDVI–NDWI index is used to enhance the contrast between water bodies and the surrounding surface features; the topographic SLOPE is used to eliminate the mountain shadow; and the NIR band is used to reduce the effects of artificial construction land. The objectives are evaluating the potential of the HJ-1A/B imagery on water body monitoring, and proposing ideally mapping method. The test study results indicated that the NDVI–NDWI index is superior to the single index of NDVI and NDWI to enhance the contrast between water bodies and the rest of the features. On the basis of the accurately mapped water bodies in the HJ-1A/B CCD images of the study area, we conclude that the HJ-1A/B multi-spectral satellite images is an ideal data source for high spatial and temporal resolution water bodies monitoring. And the integrated water body mapping method is suitable for the applications of HJ-1A/B multi-spectral satellite images in this field.     

利用轨道参数修正的无控制点星载SAR图像几何校正方法

使用距离多普勒模型进行SAR图像几何校正时,卫星轨道误差、系统成像参数误差和DEM高程的误差会影响几何校正精度。本文提出了一种基于轨道参数修正的星载SAR图像几何校正方法。首先利用多项式对卫星轨道进行参数化,然后使用模拟SAR图像与真实SAR图像进行匹配得到控制点来修正轨道参数,最后利用修正后的参数进行几何精校正,从而提高几何校正精度。该方法无需地面控制点,适用于不易于人工测量获取地面控制点地区的SAR图像几何校正,与基于模拟SAR图像匹配并使用多项式改正的几何校正方法相比,本文方法具有更高的精度。使用Radarsat-2图像进行试验,并使用地面实测GPS控制点验证了本方法的有效性。     

天绘一号03星定位精度初步评估

天绘一号卫星是我国第一颗传输型立体测绘卫星,主要用于无地面控制点条件下高精度定位及1∶5万比例尺地形图测绘。01、02星分别于2010年8月24日及2012年5月6日成功发射,03星于2015年11月26日发射,目前3颗星在轨正常运行。01星定位精度经过系统检测后,定位精度为10.3 m/5.7 m(平面/高程),02星定位精度与01星精度相当。本文首先简要介绍了天绘一号卫星的总体概况及其地面处理系统中的关键技术,重点对03星的无控定位精度进行初步试验评估。通过国内外3条航线的试验发现,03星较02星精度有较大提高,在无地面控制点条件下可实现平面7.2 m、高程2.6 m的定位精度。     

北斗二号对北斗三号伪距单点定位精度影响分析

针对当前BDS-2/BDS-3组合定位性能,基于国际GNSS服务(IGS)跟踪站对BDS-2/BDS-3的组合定位实测数据,分析了北斗二号(BDS-2)不同类型卫星对北斗三号(BDS-3)卫星B1I、B3I与B1I/B3I双频组合伪距单点定位精度的影响. 经研究发现,BDS-2不同类型星座卫星都能有效改善BDS-3卫星可见数与卫星空间几何构型,地球静止轨道(GEO)卫星、倾斜地球同步轨道(IGSO)卫星、中圆地球轨道(MEO)卫星使BDS-3单频和双频伪距单点定位精度的提升在20%以内,而全部BDS-2卫星使BDS-3单频和双频伪距单点定位精度的提升在30%以内,整体提升关系为:BDS-2全部卫星>IGSO卫星>GEO卫星>MEO卫星.      

时序InSAR技术在地球环境监测及其资源管理中的应用:以交城-清徐地区为例

利用ALOS-1（2007-2010）、Sentinel-1A（2017-2018）存档数据对山西交城-清徐地区的地面沉降进行监测。结合小基线和永久散射体技术优点,在增加时间采样密度的基础上利用二维线性回归分析得到研究区域的形变速率和时间序列。经同期GPS观测结果校核表明:交城-清徐地区持续发生地面沉降,但山区和平原区域形变的空间分布、量级不同,引起形变原因也不同。平原地带沉降空间分布受构造断裂控制,断裂带两侧呈现明显的差异性形变,且最大沉降速率为－200 mm/a,沉降的主要原因是地下水超采,但经治理后地面沉降灾害有所缓解,表现为沉降速率小于－30 mm/a。山区主要由于矿产资源的长期开采,沉降中心不断向南移动,最大形变速率为－462 mm/a。