利用Sentinel-1A数据监测地表形变案例分析

雷达卫星结合InSAR技术已广泛应用于高精度地表形变监测领域。本文选取2017年九寨沟地震为研究案例,利用Sentinel-1A地震前后的单视复数影像,基于D-InSAR技术获取该次地震的同震形变场。结果显示:震中西北侧表现出相对均匀的下沉现象,沉降漏斗区雷达视线向最大沉降量达25.1 cm;东南侧呈现不均匀抬升状态,地表破碎较为明显,最大抬升量为11.6 cm。研究表明基于Sentinel-1A数据的D-InSAR技术可以为地震形变场的定量分析提供一种快速有效的手段,为阐释地震发震机理及评估受灾情况提供必要的数据支撑,具有广阔的应用前景。     

哨兵1号和高分三号SAR数据海浪谱反演精度评估

Ocean wave is one of the important marine dynamic phenomenon that affect human activities. At present, the main observation means include buoy observation, marine numerical prediction model, and microwave remote sensing observation. However, we cannot conduct large-scale observation by buoy, and the marine numerical prediction model’s result is not measured data. Spectrometers and altimeters in microwave remote sensing instruments can also measure spectral parameters. However, SAR, which has a higher resolution, can provide 2D sea surface information. The Sentinel-1 satellite of ESA and GF-3 satellite independently developed by China are now in orbit, and numerous teams are working to retrieve wave parameters from SAR data of these two satellites. In this work, we compared the wave parameter inversion accuracy of Sentinel-1 SAR Interferometric Wide Swath model and GF-3 SAR strip model based on wave spectrum, which provides a reference for the wide application of GF-3 SAR data. The sea states according to the ERA-5 data of ECMWF are divided into three categories: low, moderate, and high sea states. The sea areas of Hormuz and Malacca Straits of the maritime Silk Road in the Indian Ocean and the coastal waters of the Pacific and Atlantic Ocean are selected as the study areas. Meanwhile, the SAR data of Sentinel-1 and GF-3 satellites under different sea states are selected as the data source. The MPI method is used to retrieve the wave spectrum and wave parameters using the E spectrum as the initial guess. Subsequently, the SAR data inversion results of the two satellites under different sea states are compared with the ERA-5 and buoy wave data. The inversion accuracy of the wave parameters can be verified by calculating the values of the Root Mean Square Error (RMSE) and Scatter Index (SI), and the inversion accuracy of the wave parameters of the two satellites under different sea conditions can be compared. The RMSEs of significant wave height (Hs) retrieved by GF-3 SAR under low, moderate, and high sea conditions are 0.30, 0.34, and 0.48 m, and those of mean wave period (Tm) are 1.02, 0.99, and 0.95 s, respectively, compared with the ERA-5 data. In addition, the RMSE of Hs retrieved by Sentinel-1 SAR under low, moderate, and high sea conditions are 0.30, 0.29, and 0.33 m, respectively, and the RMSEs of Tm are 0.94, 0.51, and 0.64 s, respectively. The RMSEs of Hs and Tm under different sea conditions retrieved by GF-3 SAR are 0.38 m and 0.99 s, and those of Hs and Tm retrieved by Sentinel-1 SAR are 0.31 m and 0.70 s, respectively, compared with the ERA-5 data. The RMSEs of the retrieved Hs and Tm of GF-3 satellite are 0.42 m and 0.94 s, and those of the retrieved Hs and Tm of Sentinel-1 are 0.40 m and 0.91 s, respectively, compared with the buoy data. The SAR wave parameter inversion of Sentinel-1 and GF-3 SAR based on the wave spectrum shows that the inversion results of the two satellites meet the index requirements in this field, and the accuracy of the inversion results of wave spectrum is the same. The strip mode SAR data of GF-3 satellite, China’s first self-developed SAR satellite, has broad prospects in marine research fields. © 2023 National Remote Sensing Bulletin. All rights reserved.     

Sentinel-1A数据的四川宜宾Ms6.0同震形变分析

为了获取2019年6月17日发生的四川宜宾Ms6.0地震引起的地表形变情况,该文利用欧空局宽幅模式的高分辨率新型Sentinel-1A卫星获取了此次地震的第一对同震干涉像对数据,使用D-InSAR技术获取宜宾市长宁县地区的同震形变场。结果显示,本次地震在震中西北方向分别形成了1个明显的沉降区和抬升区,在雷达视线方向上的最大沉降量为7.9 cm,最大抬升量为8.1 cm。通过与同一时间内的GPS高程测量形变量相比,D-InSAR解算的地表形变量与GPS监测点形变量基本一致,均不超过3 mm,表明了本文的D-InSAR形变解算结果的可靠性,体现了新型Sentinel-1A雷达卫星在地震形变监测领域有着很高的应用价值和潜力。     

Sentinel-1 observations of the 2016 Menyuan earthquake: A buried reverse event linked to the left-lateral Haiyuan fault

Knowledge on the interaction of active structures is essential to understand mechanics of continental deformation and estimate the earthquake potential in complex tectonic settings. Here we use Sentinel-1A radar imagery to investigate coseismic deformation associated with the 2016 Menyuan (Qinghai) earthquake, which occurred in the vicinity of the left-lateral Haiyuan fault. The ascending and descending interferograms indicate thrust-dominated slip, with the maximum line-of-sight displacements of 58 and 68 mm, respectively. The InSAR observations fit well with the uniform-slip dislocation models except for a larger slip-to-width ratio than that predicted by the empirical scaling law. We suggest that geometric complexities near the Leng Long Ling restraining bend confine rupture propagation, resulting in high slip occurred within a small area and much higher stress drop than global estimates. Although InSAR observations cannot distinguish the primary plane, we prefer the west-dipping solution considering aftershocks distribution and the general tectonic context. Both InSAR modelling and aftershock locations indicate that the rupture plane linked to the Haiyuan fault at 10 km depth, a typical seismogenic depth in Tibet. We suggest that the earthquake more likely occurred on a secondary branch at a restraining bend of the Haiyuan fault, even though we cannot completely rule out the possibility of it being on a splay of the North Qilian Shan thrusts.     

Assessing the suitability of data from Sentinel-1A and 2A for crop classification

Sentinel-1A C-SAR and Sentinel-2A MultiSpectral Instrument (MSI) provide data applicable to the remote identification of crop type. In this study, six crop types (beans, beetroot, grass, maize, potato, and winter wheat) were identified using five C-SAR images and one MSI image acquired during the 2016 growing season. To assess the potential for accurate crop classification with existing supervised learning models, the four different approaches namely kernel-based extreme learning machine (KELM), multilayer feedforward neural networks, random forests, and support vector machine were compared. Algorithm hyperparameters were tuned using Bayesian optimization. Overall, KELM yielded the highest performance, achieving an overall classification accuracy of 96.8%. Evaluation of the sensitivity of classification models and relative importance of data types using data-based sensitivity analysis showed that the set of VV polarization data acquired on 24 July (Sentinel-1A) and band 4 data (Sentinel-2A) had the greatest potential for use in crop classification.     

多视角InSAR数据解算2017两伊地震三维同震形变场

2017两伊地震是自1900年以来发生在扎格罗斯山脉的最大地震,为了研究此次地震引起的同震形变场,利用覆盖同一地区的3对Sentinel-1A升降轨数据分别进行两通差分DInSAR处理,得到了研究区3个视线向的地表同震形变场,通过直接解算法重建了研究区的三维同震形变场。试验表明：3种视角的升降轨视线向上升与沉降总体趋势基本一致;联合多个视角的观测结果可以实现三维形变场的重建;根据地表视线向和三维同震形变的特征以及地质构造背景推测了发震断层很有可能为扎格罗斯山前断层。     

时序Sentinel-1A数据支持的长江中下游汛情动态监测

合成孔径雷达（SAR）因其对地观测全天候、全天时优势,成为多云多雨天气限制下洪水动态监测中不可或缺的数据来源之一。由于GEE（Google Earth Engine）云计算平台的兴起和短重访Sentinel-1数据的可获取性,洪水监测与灾害评估目前正面向动态化、广域化快速发展。顾及洪水淹没区土地覆盖变化的复杂性和发生时间的不确定性,基于时序Sentinel-1A卫星数据提出了针对大尺度范围、连续长期的汛情自动检测及动态监测方法。该方法首先,利用图像二值化分割时序SAR数据实现水体时空分布粗制图,逐像素计算时间序列中被识别为水体候选点的频率。然后,利用Sentinel-2光学影像对精度较粗的初期SAR水体提取结果进行校正,得到精细的水体分布图。最后,针对不同频率区间的淹没特点,采用差异化的时序异常检测策略识别淹没范围：对低频覆水区利用欧氏距离检测时序断点,以提取扰动强度大、淹没时间短的洪涝灾害区;对高频覆水区利用标准分数（Z-Score）检测时序断点,以提取季节性水体覆盖区。在GEE平台上利用该方法,实现了2020-05—10长江中下游地区全域洪水淹没范围时空信息的自动、快速、有效监测,揭示了不同区域汛情发展模式的差异性。本文提出的洪水快速监测方法对大尺度下的汛情动态监测、灾害定量评估和快速预警响应具有重要的现实意义。