Resolving co- and early post-seismic slip variations of the 2021 MW 7.4 Madoi earthquake in east Bayan Har block with a block-wide distributed deformation mode from satellite synthetic aperture radar data

On 21 May 2021(UTC),an MW 7.4 earthquake jolted the east Bayan Har block in the Tibetan Plateau.The earthquake received widespread attention as it is the largest event in the Tibetan Plateau and its surroundings since the 2008 Wenchuan earthquake,and especially in proximity to the seismic gaps on the east Kunlun fault.Here we use satellite interferometric synthetic aperture radar data and subpixel offset observations along the range directions to characterize the coseismic deformation of the earthquake.Range offset displacements depict clear surface ruptures with a total length of~170 km involving two possible activated fault segments in the earthquake.Coseismic modeling results indicate that the earthquake was dominated by left-lateral strike-slip motions of up to 7 m within the top 12 km of the crust.The well-resolved slip variations are characterized by five major slip patches along strike and 64%of shallow slip deficit,suggesting a young seismogenic structure.Spatial-temporal changes of the postseismic deformation are mapped from early 6-day and 24-day InSAR observations,and are well explained by time-dependent afterslip models.Analysis of Global Navigation Satellite System(GNSS)velocity profiles and strain rates suggests that the eastward extrusion of plateau is diffusely distributed across the east Bayan Har block,but exhibits significant lateral heterogeneities,as evidenced by magnetotelluric observations.The block-wide distributed deformation of the east Bayan Har block along with the significant co-and post-seismic stress loadings from the Madoi earthquake imply high seismic risks along regional faults,especially the Tuosuo Lake and Maqên-Maqu segments of the Kunlun fault that are known as seismic gaps.

Application of Interferometric Spaceborne Synthetic Aperture Radar to Morphotectonic Research in the Bengcuo Area of Tibet

Using ELLS-1/2 radar satellite data, we generated a Digital Elevation Model by Interferometric Synthetic Aperture Radar, and the ERS-1/2 DEM that we generated is obviously superior to both the 1/25 ten thousand-scale DEM of the National Fundamental Geographic Information System of China and the 90-m spatial resolution＇s SRTM DEM that America published in 2004 when it showed the characteristics of tiny structure relief. By analyzing the relief characteristics of the Bengcuo fault zone based on ERS-1/2 DEM, we find that the relief on the connection location of the Bengcuo and Pengcuo fault zones has complex characteristics. A structure relief that is similar to the Pengcuo fault zone crosses through the Dazi-Dasha fault on the the Bengcuo fault zone, while the Dazi-Dasha fault crosses through a gully at this place. This indicates that the Dazi-Dasha fault has been active at this place recently. At the same time, the Naka-Naduiduo fault is severed by the gully which was cut through by the Dazi-Dasha fault. Therefore, the Naka-Naduiduo fault was formed earlier than the Dazi-Dasha fault.

Effect of digital elevation models on monitoring slope displacements in open-pit mine by differential interferometry synthetic aperture radar

Displacement monitoring in open-pit mines is one of the important tasks for safe management of mining processes.Differential interferometric synthetic aperture radar(DInSAR),mounted on an artificial satellite,has the potential to be a cost-effective method for monitoring surface displacements over extensive areas,such as open-pit mines.DInSAR requires the ground surface elevation data in the process of its analysis as a digital elevation model(DEM).However,since the topography of the ground surface in open-pit mines changes largely due to excavations,measurement errors can occur due to insufficient information on the elevation of mining areas.In this paper,effect of different elevation models on the accuracy of the displacement monitoring results by DInSAR is investigated at a limestone quarry.In addition,validity of the DInSAR results using an appropriate DEM is examined by comparing them with the results obtained by global positioning system(GPS)monitoring conducted for three years at the same limestone quarry.It is found that the uncertainty of DEMs induces large errors in the displacement monitoring results if the baseline length of the satellites between the master and the slave data is longer than a few hundred meters.Comparing the monitoring results of DInSAR and GPS,the root mean square error(RMSE)of the discrepancy between the two sets of results is less than 10 mm if an appropriate DEM,considering the excavation processes,is used.It is proven that DInSAR can be applied for monitoring the displacements of mine slopes with centimeter-level accuracy.

A synthetic aperture radar sea surface distribution estimation by n-order Bézier curve and its application in ship detection

To dates,most ship detection approaches for single-pol synthetic aperture radar（SAR） imagery try to ensure a constant false-alarm rate（CFAR）.A high performance ship detector relies on two key components：an accurate estimation to a sea surface distribution and a fine designed CFAR algorithm.First,a novel nonparametric sea surface distribution estimation method is developed based on n-order Bézier curve.To estimate the sea surface distribution using n-order Bézier curve,an explicit analytical solution is derived based on a least square optimization,and the optimal selection also is presented to two essential parameters,the order n of Bézier curve and the number m of sample points.Next,to validate the ship detection performance of the estimated sea surface distribution,the estimated sea surface distribution by n-order Bézier curve is combined with a cell averaging CFAR（CA-CFAR）.To eliminate the possible interfering ship targets in background window,an improved automatic censoring method is applied.Comprehensive experiments prove that in terms of sea surface estimation performance,the proposed method is as good as a traditional nonparametric Parzen window kernel method,and in most cases,outperforms two widely used parametric methods,K and G0 models.In terms of computation speed,a major advantage of the proposed estimation method is the time consuming only depended on the number m of sample points while independent of imagery size,which makes it can achieve a significant speed improvement to the Parzen window kernel method,and in some cases,it is even faster than two parametric methods.In terms of ship detection performance,the experiments show that the ship detector which constructed by the proposed sea surface distribution model and the given CA-CFAR algorithm has wide adaptability to different SAR sensors,resolutions and sea surface homogeneities and obtains a leading performance on the test dataset.

A Review on Applications of Imaging Synthetic Aperture Radar with a Special Focus on Cryospheric Studies

The cryosphere is the frozen part of the Earth’s system. Snow and ice are the main constituents of the cryosphere and may be found in different states, such as snow, freshwater ice, sea ice, perma-frost, and continental ice masses in the form of glaciers and ice sheets. The present review mainly deals with state-of-the-art applications of synthetic aperture radar (SAR) with a special emphasize on cryospheric information extraction. SAR is the most important active microwave remote sensing (RS) instrument for ice monitoring, which provides high-resolution images of the Earth’s surface. SAR is an ideal sensor in RS technology, which works in all-weather and day and night conditions to provide useful unprecedented information, especially in the cryospheric regions which are almost inaccessible areas on Earth. This paper addresses the technological evolution of SAR and its applications in studying the various components of the cryosphere. The arrival of SAR radically changed the capabilities of information extraction related to ice type, new ice formation, and ice thickness. SAR applications can be divided into two broad classes-polarimetric applications and interferometric applications. Polarimetric SAR has been effectively used for mapping calving fronts, crevasses, surface structures, sea ice, detection of icebergs, etc. The paper also summarizes both the operational and climate change research by using SAR for sea ice parameter detection. Digital elevation model (DEM) generation and glacier velocity mapping are the two most important applications used in cryosphere using SAR interferometry or interferometric SAR (InSAR). Space-borne InSAR techniques for measuring ice flow velocity and topography have developed rapidly over the last decade. InSAR is capable of measuring ice motion that has radically changed the science of glaciers and ice sheets. Measurement of temperate glacier velocities and surface characteristics by using airborne and space-borne interferometric satellite images have been the significant application in glaciology and cryospheric studies. Space-borne InSAR has contributed to major evolution in many research areas of glaciological study by measuring ice-stream flow velocity, improving understanding of ice-shelf processes, yielding velocity for flux-gate based mass-balance assessment, and mapping flow of mountain glaciers. The present review summarizes the salient development of SAR applications in cryosphere and glaciology.

Effect of orbital errors on the geosynchronous circular synthetic aperture radar imaging and interferometric processing

The geosynchronous circular synthetic aperture radar (GEOCSAR) is an innovative SAR system,which can produce high resolution three-dimensional (3D) images and has the potential to provide 3D deformation measurement.With an orbit altitude of approximately 36 000 km,the orbit motion and orbit disturbance effects of GEOCSAR behave differently from those of the conventional spaceborne SAR.In this paper,we analyze the effects of orbit errors on GEOCSAR imaging and interferometric processing.First,we present the GEOCSAR imaging geometry and the orbit errors model based on perturbation analysis.Then,we give the GEOCSAR signal formulation based on imaging geometry,and analyze the effect of the orbit error on the output focused signal.By interferometric processing on the 3D reconstructed images,the relationship between satellite orbit errors and the interferometric phase is deduced.Simulations demonstrate the effects of orbit errors on the GEOCSAR images,interferograms,and the deformations.The conclusions are that the required relative accuracy of orbit estimation should be at centimeter level for GEOCSAR imaging at L-band,and that millimeter-scale accuracy is needed for GEOCSAR interferometric processing.

Formation flying orbit design for the distributed synthetic aperture radar satellite

Formation flying orbit design is one of the key technologies for system design and performance analysis of the distributed SAR satellites. The approximately analytic solution of the passive stable formation flying orbit elements is explored based on the expansion form of Kepler's equation. A new method of orbital parameters design for three-dimensional formation flying SAR satellites is presented, and the precision of the orbital elements is analyzed. Formation flying orbit elements are calculated for the L-Band distributed SAR satellites using the formulas deduced in this paper. The accuracy of the orbital elements is validated by the computer simulation results presented in this paper.

GF-3 data real-time processing method based on multi-satellite distributed data processing system

Due to the limited scenes that synthetic aperture radar(SAR)satellites can detect,the full-track utilization rate is not high.Because of the computing and storage limitation of one satellite,it is difficult to process large amounts of data of spaceborne synthetic aperture radars.It is proposed to use a new method of networked satellite data processing for improving the efficiency of data processing.A multi-satellite distributed SAR real-time processing method based on Chirp Scaling(CS)imaging algorithm is studied in this paper,and a distributed data processing system is built with field programmable gate array(FPGA)chips as the kernel.Different from the traditional CS algorithm processing,the system divides data processing into three stages.The computing tasks are reasonably allocated to different data processing units(i.e.,satellites)in each stage.The method effectively saves computing and storage resources of satellites,improves the utilization rate of a single satellite,and shortens the data processing time.Gaofen-3(GF-3)satellite SAR raw data is processed by the system,with the performance of the method verified.

AN ALGORITHM FOR INTERFEROMETRIC SAR DATA PROCESSING

In this paper, an algorithm of generating INSAR unwrapped phase image from SAR single-look complex images is presented. Besides the general processing technique, this article focuses on the methods of flat-earth phase removal, phase noise reduction and phase unwrapping. The availability is tested by the results of processing ERS-1/2 SAR images.

Frequency Domain Filtering SAR Interferometric Phase Noise Using the Amended Matrix Pencil Model

Interferometric phase filtering is one of the key steps in interferometricsynthetic aperture radar (InSAR/SAR). However, the ideal filtering results are difficult toobtain due to dense fringe and low coherence regions. Moreover, the InSAR/SAR datarange is relatively large, so the efficiency of interferential phase filtering is one of themajor problems. In this letter, we proposed an interferometric phase filtering methodbased on an amended matrix pencil and linear window mean filter. The combination ofthe matrix pencil and the linear mean filter are introduced to the interferometric phasefiltering for the first time. First, the interferometric signal is analyzed, and theinterferometric phase filtering is transformed into a local frequency estimation problem.Then, the local frequency is estimated using an amended matrix pencil at a window. Thelocal frequency can represent terrain changes, thus suggesting that the frequency can beaccurately estimated even in dense fringe regions. Finally, the local frequency is filteredby using a linear window mean filter, and the filtered phase is recovered. The proposedmethod is calculated by some matrices. Therefore, the computational complexity isreduced, and the efficiency of the interferometric phase filtering is improved.Experiments are conducted with simulated and real InSAR data. The proposed methodexhibits a better filtering effect and an ideal efficiency as compared with the traditionalfiltering method.

Phase synchronization processing method for alternating bistatic mode in distributed SAR

In the distributed synthetic aperture radar （SAR）, the alternating bistatic mode can perform phase reference without a synchronization link between two satellites compared with the pulsed alternate synchronization method. The key of the phase synchronization processing is to extract the oscillator phase differences from the bistatic echoes. A signal model of phase synchronization in the alternating bistatic mode is presented. The phase synchronization processing method is then studied. To reduce the phase errors introduced by SAR imaging, a sub-aperture processing method is proposed. To generalize the sub-aperture processing method, an echo-domain processing method using correlation of bistatic echoes is proposed. Finally, the residual phase errors of the both proposed processing methods are analyzed. Simulation experiments validate the proposed phase synchronization processing method and its phase error analysis results.

A statistical distribution texton feature for synthetic aperture radar image classification

We propose a novel statistical distribution texton(s-texton) feature for synthetic aperture radar(SAR) image classification. Motivated by the traditional texton feature, the framework of texture analysis, and the importance of statistical distribution in SAR images, the s-texton feature is developed based on the idea that parameter estimation of the statistical distribution can replace the filtering operation in the traditional texture analysis of SAR images. In the process of extracting the s-texton feature, several strategies are adopted, including pre-processing, spatial gridding, parameter estimation, texton clustering, and histogram statistics. Experimental results on Terra SAR data demonstrate the effectiveness of the proposed s-texton feature.

Three-dimensional bistatic interferometric ISA Rimaging

An approach based on interferometry technique is proposed for three-dimensional （3D） bistatic inverse synthetic aperture radar （ISAR） imaging. It is converted to a monostatic problem by using the theory that a bistatic radar equals a monostatic radar located on the bisector of bistatic an- gle. Then, interferometric phases extracted from a pair of cross shaped antennas are used to esti- mate the height and associated rotational velocity. Finally, numerical simulations are provided to e- valuate this method.

利用InSAR技术反演甘肃积石山M_(S)6.2地震同震破裂模型及形变时间序列

使用升降轨Sentinel-1A卫星数据,利用差分干涉技术获取甘肃积石山M_(S)6.2地震高精度同震形变场,并以此为约束反演同震破裂模型,然后采用SBAS-InSAR技术获取震中附近区域近6个月的累积形变时间序列。结果表明,本次地震地表形变以抬升为主,升降轨数据视线向最大同震形变量分别约为6.8 cm和7.3 cm,破裂主要集中在7~15 km深度范围内,最大滑动量约为0.55 m;震后2个月内,震中近场及西北侧的抬升形变呈扩大趋势。综合分析认为,此次积石山地震发震断层走向NW、倾向NE,隶属于拉脊山南缘断裂带东南侧的隐伏分支断层;震后的大面积突发形变可能与孔隙压力释放有关,地表稳定状态有待进一步观察。

Coseismic slip distribution of 2009 L'Aquila earthquake derived from InSAR and GPS data

To better understand the mechanism of the Mw6.3 L＇Aquila （Central Italy） earthquake occurred in 2009, global positioning system （GPS） and interferometric synthetic aperture radar （InSAR） data were used to derive the coseismic slip distribution of the earthquake fault. Firstly, based on the homogeneous elastic half-space model, the fault geometric parameters were solved by the genetic algorithm. The best fitting model shows that the fault is a 13.7 km×14.1 km rectangular fault, in 139.3° strike direction and 50.2° southwest-dipping. Secondly, fixing the optimal fault geometric parameters, the fault plane was extended and discretized into 16× 16 patches, each with a size of 1 kmx 1 krn, and the non-uniform slip distribution of the fault was inverted by the steepest descent method with an appropriate smoothing ratio based on the layered crustal structure model. The preferred solution shows that the fault is mainly a normal fault with slight right-lateral strike slip, the maximum slip of 1.01 m is located in the depth of 8.28 km, the average rake is -100.9°, and the total geodetic moment is about 3.34× 1018 N.m （Mw 6.28）. The results are much closer than previous studies in comparison with the seismological estimation. These demonstrate that the coseismic fault slip distribution of the L＇Aauila earthauake inverted by the crustal model considering layered characters is reliable.

一种天基双基地InSAR系统的相位同步方法及验证情况

以德国的TanDEM-X双星编队干涉合成孔径雷达(InSAR)卫星系统为背景,提出了基于点频连续波同步信号双向对传,实现编队卫星InSAR系统相位同步的方法。文中给出了该同步方法实现的框架,对该相位同步方法的原理进行了理论推导,得出了补偿相位差的表达式,并对影响同步精度的因素进行了初步的分析。该方法可用于天基双/多基地合成孔径雷达系统相位同步链路的设计,仿真分析的结果验证了该方法的理论可行性。此外,通过搭建实物系统对系统的工程可实现精度进行了验证,结果表明获取的载波频率相位同步精度与倍频器件的频率精度和一致性密切相关。该相位同步方法为与TanDEM-X相近InSAR系统的设计提供参考。

MCJ-UNet:一种双/多通道联合InSAR相位解缠网络

干涉合成孔径雷达(InSAR)可实现地表高程的高效获取,在地形测绘中应用广泛。双/多通道InSAR技术可借助不同通道(基线、频点)的高程模糊度差异,解决相位欠采样问题,完成高程陡变区域的干涉相位解缠,实现InSAR技术在测绘困难区域的有效应用。该文即面向高效高精度相位解缠需求,利用深度学习这一有力工具,结合不同通道的相位特征及相互约束关系,提出了一种双/多通道联合干涉相位解缠网络:Multi-Channel-Joint-UNet(MCJ-UNet)。该网络的构建以双通道(双频、双基线)InSAR为基本观测构型,并可实现向多通道构型的扩展,其构建的核心思路主要包括3点:首先,将干涉相位解缠中的模糊数估计问题转化为语义分割问题,并采用UNet网络完成分割处理;其次,引入挤压激励模块(SE)动态调整信息权重,以增强网络不同通道对其所需信息的感知能力;最后,利用多通道联合约束下的相位残差优化损失函数,实现网络调谐。此外,为避免语义分割结果的边缘细节误差对解缠效果的影响,该文还提出了一种基于多通道联合约束的解缠误差自修正方法,以保证解缠质量。模拟地形仿真数据、真实地形仿真数据以及TerraSAR-X实测数据验证了所提方法的有效性。

融合长短基线的干涉图解缠错误自动改正

合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)技术因具有全天时、全天候、大范围的独特优势,已被广泛应用于地球科学研究、自然灾害监测和环境管理等领域。然而,受地表失相干、干涉条纹密集等因素影响,InSAR数据处理关键的步骤——干涉图相位解缠常出现错误,进而影响高程反演和地表形变监测精度。为此,文章提出了一种融合长短基线干涉相位闭合环的解缠错误自动探测与改正方法,并利用湖南省-江西省高植被覆盖区的Sentinel-1雷达影像开展了算法验证。实验结果表明,该方法可有效地改正干涉相位解缠错误,这对InSAR技术提升高程反演和形变监测精度具有重要意义。

基于InSAR技术的大跨桥梁温度变形监测研究

以国内某高速铁路钢拱桥为研究对象,选取2017—2018年期间59幅C波段Senti⁃nel-1号雷达卫星影像,利用PS-InSAR技术处理影像获得桥梁的视线向(Line of Sight,LOS)位移,根据SAR成像空间几何关系解算出支座的纵向位移.研究结果表明:支座纵向位移的时空特性与实际桥梁结构相符合,验证了PS-InSAR技术观测桥梁结构位移的可行性.建立支座纵向位移与温度的线性相关模型,并与结构健康监测系统的实测结果进行对比.两者吻合良好,相对误差控制在10%以内,验证了PS-InSAR测量桥梁结构位移的可靠性.利用有限元模拟温度作用下桥梁支座的位移变化,并与PS-InSAR位移时间序列进行对比.两者趋势基本一致,LOS向位移误差在[-10,10]mm,验证了PS-InSAR测量桥梁结构位移的准确性.

基于时序InSAR技术的山区高速公路挖方边坡监测

针对传统人工水准监测风险高、难度大的问题,以本桓高速公路挖方边坡为例,提出一种利用SBAS-InSAR与PS-InSAR相结合的技术对挖方边坡进行形变监测。首先,采用SBAS-InSAR技术处理2023年4月—2024年1月获取的22景Sentinel-1A数据,生成小基线差分干涉对,反演得出边坡形变值与形变趋势。然后,选取相同监测时间的水准点测量结果与PS点进行验证,并通过灰色关联度分析降雨量与边坡形变的关联度。监测结果显示,边坡沉降速率集中在10.15~30.27 mm/a,沉降分析结果显示此路堑边坡整体较为稳定,时序InSAR相对水准测量具有全天时、全天候、全维度的特点,并具有可回溯性,表明时序InSAR技术对山区高速公路挖方边坡监测具有可行性与可靠性,确保高速公路挖方边坡安全稳定。