Improving Satellite-Retrieved Cloud Base Height with Ground-Based Cloud Radar Measurements

Cloud base height(CBH) is a crucial parameter for cloud radiative effect estimates, climate change simulations, and aviation guidance. However, due to the limited information on cloud vertical structures included in passive satellite radiometer observations, few operational satellite CBH products are currently available. This study presents a new method for retrieving CBH from satellite radiometers. The method first uses the combined measurements of satellite radiometers and ground-based cloud radars to develop a lookup table(LUT) of effective cloud water content(ECWC), representing the vertically varying cloud water content. This LUT allows for the conversion of cloud water path to cloud geometric thickness(CGT), enabling the estimation of CBH as the difference between cloud top height and CGT. Detailed comparative analysis of CBH estimates from the state-of-the-art ECWC LUT are conducted against four ground-based millimeter-wave cloud radar(MMCR) measurements, and results show that the mean bias(correlation coefficient) is0.18±1.79 km(0.73), which is lower(higher) than 0.23±2.11 km(0.67) as derived from the combined measurements of satellite radiometers and satellite radar-lidar(i.e., Cloud Sat and CALIPSO). Furthermore, the percentages of the CBH biases within 250 m increase by 5% to 10%, which varies by location. This indicates that the CBH estimates from our algorithm are more consistent with ground-based MMCR measurements. Therefore, this algorithm shows great potential for further improvement of the CBH retrievals as ground-based MMCR are being increasingly included in global surface meteorological observing networks, and the improved CBH retrievals will contribute to better cloud radiative effect estimates.

Time-varying baseline error correction method for ground-based micro-deformation monitoring radar

In recent years, ground-based micro-deformation monitoring radar has attracted much attention due to its excellent monitoring capability. By controlling the repeated campaigns of the radar antenna on a fixed track, ground-based micro-deformation monitoring radar can accomplish repeat-pass interferometry without a space baseline and thus obtain highprecision deformation data of a large scene at one time. However, it is difficult to guarantee absolute stable installation position in every campaign. If the installation position is unstable, the stability of the radar track will be affected randomly, resulting in time-varying baseline error. In this study, a correction method for this error is developed by analyzing the error distribution law while the spatial baseline is unknown. In practice, the error data are first identified by frequency components, then the data of each one-dimensional array(in azimuth direction or range direction) are grouped based on numerical distribution period, and finally the error is corrected by the nonlinear model established with each group.This method is verified with measured data from a slope in southern China, and the results show that the method can effectively correct the time-varying baseline error caused by rail instability and effectively improve the monitoring data accuracy of groundbased micro-deformation radar in short term and long term.

Constrained geometry analysis to resolve 3-D deformations from three ground-based radars

When multiple ground-based radars(GB-rads)are utilized together to resolve three-dimensional(3-D)deformations,the resolving accuracy is related with the measurement geometry constructed by these radars.This paper focuses on constrained geometry analysis to resolve 3-D deformations from three GB-rads.The geometric dilution of precision(GDOP)is utilized to evaluate 3-D deformation accuracy of a single target,and its theoretical equation is derived by building a simplified 3-D coordinate system.Then for a 3-D scene,its optimal accuracy problem is converted into determining the minimum value of an objective function with a boundary constraint.The genetic algorithm is utilized to solve this constrained optimization problem.Numerical simulations are made to validate the correctness of the theoretical analysis results.

An approach to wide-field imaging of linear rail ground-based SAR in high squint multi-angle mode

Ground-based synthetic aperture radar(GB-SAR) has been successfully applied to the ground deformation monitoring.However, due to the short length of the GB-SAR platform, the scope of observation is largely limited. The practical applications drive us to make improvements on the conventional linear rail GB-SAR system in order to achieve larger field imaging. First, a turntable is utilized to support the rotational movement of the radar.Next, a series of high-squint scanning is performed with multiple squint angles. Further, the high squint modulation phase of the echo data is eliminated. Then, a new multi-angle imaging method is performed in the wave number domain to expand the field of view. Simulation and real experiments verify the effectiveness of this method.

Comparative research and its significance of deformation measurements by technologies of laser real-time holographic interferome-try and radar differential interferometry

The principles and applications of laser real-time holographic interferometry （LRTHI） and radar differential interferometry （RDI） technologies are described in this paper, respectively. By using LRTHI, we can observe the deformation of samples under pressure in the lab and study the anomaly characteristics relating to different strain fields in different fracture-developing areas; while by using RDI, we can observe the landform and surface deformation. The results of deformation observed before and after the Ms=7.9 Mani earthquake （Tibet） and Ms=6.2 Shangyi-Zhangbei earthquake in China are obtained. It is pointed out that LRTHI and RDi are similar, which study the characteristics of anomalous deformation field by fringe variations for both of them. Therefore, the observation of deformation field in the seismogenic process, especially in the period impending an earthquake by RDI, and the comparative study in the lab by LRTHI are of great significance.

Ground Surface Ruptures and Near-Fault,Large-Scale Displacements Caused by the Wenchuan Ms8.0 Earthquake Derived from Pixel Offset Tracking on Synthetic Aperture Radar Images

The 12 May 2008 Wenchuan Ms8.0 earthquake produced surface displacements along the causative fault, the Yingxiu-Beichuan Fault, which are up to several meters near the fault. Because of the large gradient, satellite synthetic aperture radar （SAR） interferometric data are strongly incoherent; the usual SAR interferometry method does not allow such displacements to be measured. In the present study, we employed another approach, the technique based on pixel offset tracking, to solve this problem. The used image data of six tracks are from the Advanced Land Observing Satellite, Phased Array type L-band Synthetic Aperture Radar （ALOS/PALSAR） dataset of Japan. The results show that the entire surface rupture belt is 238 km long, extending almost linearly in a direction of 42°north-east. It is offset left laterally by a north-west-striking fault at Xiaoyudong, and turns at Gaochuan, where the rupture belt shifts toward the south by 5 km, largely keeping the original trend. In terms of the features of the rupture traces, the rupture belt can be divided into five sections and three types. Among them, the Beichuan-Chaping and Hongkou-Yingxiu sections are relatively complex, with large widths and variable traces along the trend. The Pingtong-Nanba and Qingping-Jingtang sections appear uniform, characterized by straight traces and small widths. West of Yingxiu, the rupture traces are not clear. North of the rupture belt, surface displacements are 2.95 m on average, mostly 2-3.5 m, with 7-9 m the maximum near Beichuan. South of the rupture belt, the average displacement is 1.75 m, dominated by 1-2 m, with 3-4 m at a few sites. In the north, the displacements in the radar line of sight are of subsidence, and in the south, they are uplifted, in accordance with a right-slip motion that moves the northern wall of the fault to the east, and the southern wall to the west, respectively. Along the Guanxian-Jiangyou Fault, there is a uplift zone in the radar line of sight, which is 66 km long, 1.5-6 km wide, and has vertical displacements of approximately 2 m, but no observable rupture traces.

Monitoring elevation change of glaciers on Geladandong Mountain using Tan DEM-X SAR interferometry

Glaciers play an important role in the climate system. The elevation change of a glacier is an important parameter in studies of glacier dynamics. Only a few ground-based measurements of high mountain glaciers are available due to their remoteness, high elevation, and complex topography. The acquisition from the German Tan DEM-X(Terra SAR-X add-on for Digital Elevation Measurement) SAR imaging configuration provides a reliable data sources for studying the elevation change of glaciers. In this study, the bistatic Tan DEM-X data that cover the Geladandong Mountain on the Tibetan Plateau were processed with SAR interferometry technique and the elevation changes of the mountain's glaciers during 2000–2014 were obtained. The results indicated that although distinct positive and negative elevation changes were found for different glacier tongues, the mean elevation change was about-0.14±0.26 m a-1. Geoscience Laser Altimeter System(GLAS) data were obtained for comparison and verification. The investigation using GLAS data demonstrated the efficacy of the proposed method in determining glacier elevation change. Thus, the presented approach is appropriate for monitoring glacier elevation change and it constitutes a valuable tool for studies of glacier dynamics.

Precise Three-Dimensional Deformation Retrieval in Large and Complex Deformation Areas via Integration of Offset-Based Unwrapping and Improved Multiple-Aperture SAR Interferometry:Application to the 2016 Kumamoto Earthquake

Conventional synthetic aperture radar(SAR)interferometry(InSAR)has been successfully used to precisely measure surface deformation in the line-of-sight(LOS)direction,while multiple-aperture SAR interferometry(MAI)has provided precise surface deformation in the along-track(AT)direction.Integration of the InSAR and MAI methods enables precise measurement of the two-dimensional(2D)deformation from an interferometric pair;recently,the integration of ascending and descending pairs has allowed the observation of precise three-dimensional(3D)deformation.Precise 3D deformation measurement has been applied to better understand geological events such as earthquakes and volcanic eruptions.The surface deformation related to the 2016 Kumamoto earthquake was large and complex near the fault line;hence,precise 3D deformation retrieval had not yet been attempted.The objectives of this study were to①perform a feasibility test of precise 3D deformation retrieval in large and complex deformation areas through the integration of offset-based unwrapped and improved multiple-aperture SAR interferograms and②observe the 3D deformation field related to the 2016 Kumamoto earthquake,even near the fault lines.Two ascending pairs and one descending the Advanced Land Observing Satellite-2(ALOS-2)Phased Array-type L-band Synthetic Aperture Radar-2(PALSAR-2)pair were used for the 3D deformation retrieval.Eleven in situ Global Positioning System(GPS)measurements were used to validate the 3D deformation measurement accuracy.The achieved accuracy was approximately 2.96,3.75,and 2.86 cm in the east,north,and up directions,respectively.The results show the feasibility of precise 3D deformation measured through the integration of the improved methods,even in a case of large and complex deformation.

Designing and implementation of variance-dependent Goldstein radar interferogram filtering

The variance-dependent Goldstein radar interferogram filter takes into account the information of both interferometric coherence and multilook factors,and can produce very consistent results for interferograms generated under a wide variety of multilook factors and with very different noise level.However,the filter is a bit complicated and its application is still very limited.We present the designing and implementation of the variance-dependent Goldstein radar interferogram filtering,emphasizing on the logic flow,the generation of look-up table,the determination of filtering parameter,and the handling of edge information loss.Experiments with real interferograms are provided to demonstrate the applications of the designed filtering.Comparisons with the result of the coherence-dependent Goldstein filter show that improvements from 18.4% to 36.9% are achieved when the variance-dependent filter is used,and the noisier the interferogram,the greater the improvement.

Temporal decorrelation model for the bistatic SAR interferometry

This paper develops a temporal decorrelation model for the bistatic synthetic aperture radar（BSAR） interferometry. The temporal baseline is one of the important decorrelation sources for the repeat-pass synthetic aperture radar（SAR） interferometry. The study of temporal decorrelation is challenging, especially for the bistatic configuration, since temporal decorrelation is related to the data acquisition geometry. To develop an appropriate theoretical model for BSAR interferometry, the existing models for monostatic SAR cases are extended, and the general BSAR geometry configuration is involved in the derivation. Therefore, the developed temporal decorrelation model can be seen as a general model.The validity of the theoretical model is supported by Monte Carlo simulations. Furthermore, the impacts of the system parameters and BSAR geometry configurations on the temporal decorrelation model are discussed briefly.

ADAPTIVE MODEL-BASED SCATTERING DECOMPOSITION OF POLARIMETRIC SAR INTERFEROMETRY

In this paper,a new decomposition method is proposed to solve the problems that vegetation component is overestimated and is not sensitive to directional scattering features with traditional polarimetric Synthetic Aperture Radar(SAR)decomposition.It uses a Polarimetric Interferometric Similarity Parameter(PISP)calculated from Polarimetric SAR Interferometry(PolInSAR)datasets to the scattering decomposition.The PISP is proposed to reveal the geometric sensitivity of SAR interferometry.It is defined by three optimized mechanisms obtained from PolInSAR datasets,therefore,it not only relates to the coherent scattering mechanism closely,but also sufficiently uses the phase and amplitude information.The PISP of building is high,and forest’s PISP is low.The proposed method uses the PISP as a judge condition to select different vegetation model adaptively.The decomposition results show the proposed method can effectively solve the vegetation ingredients overestimation problem.In addition,it is sensitive to the directional scattering.

利用InSAR形变反演滑坡滑面运动特征——以汶川县布瓦村滑坡为例

当前滑坡滑面运动特征获取多采用接触式手段,费时费力、成本高昂,且仅能获取稀疏点位滑面运动状态。因此,本文提出一种非接触式滑坡滑面运动特征获取方法,其基本思路是:基于位错理论建立滑坡表面形变与滑面运动间的关系模型,同时引入滑面运动平滑约束条件,最终以InSAR等技术获取的地表形变为约束,反演获得滑面运动特征。本文选择四川省汶川县布瓦村滑坡作为研究对象,以时序InSAR技术获取的滑坡坡面升降轨形变数据为约束,基于发展的理论方法反演获得布瓦村滑坡滑面几何与运动学参数,并结合野外勘查、Lidar观测以及模型残差较好地验证了反演结果的可靠性,其中反演升降轨模型残差仅3.1 mm·a^(-1)和3.4 mm·a^(-1)。反演结果显示,布瓦村滑坡滑面位于坡表以下平均~15.0 m深度处,滑面整体以沿坡向运动为主,最大运动量级达~80 mm·a^(-1),同时在坡体西侧中下部发现微量的垂直坡向运动。本文提出方法可方便、快速的获取空间连续的滑坡滑面运动特征,为滑坡灾害预警和失稳风险评估提供更为直接可靠的科学参考。

基于最优PS点获取方法的矿山工业广场沉降监测

为保障矿山地面和井下工作人员、工业广场周边居民以及运煤系统运行安全,实时监测地表变形情况,提出永久散射体(PS)点的优化获取方法,以提升PS-合成孔径雷达干涉测量(InSAR)技术在沉降监测中的适用性。首先,选择淮北市郊煤矿工业广场煤柱Ⅱ513工作面开采沉陷期间2020年10月31日—2022年8月22日共26景SAR卫星影像,运用PS-InSAR技术确定最优的PS点目标选取方法;然后,基于该方法获取该工作面回采过程中的开采损害保护区地面沉降速率、累计沉降量;最后,基于地面水准点实测数据验证PS-InSAR监测精度并分析研究区地表及建(构)筑物动态沉降情况。结果表明:相干系数和振幅离差指数双阈值法较适应于研究区沉降监测;研究区内可探测到的最大下沉速率为-26.5 mm/a,最大累计下沉值为-53.7 mm;同时,探测到研究区西北部存在岩溶塌陷地质灾害因素;利用水准数据验证精度,均方根误差(RMSE)仅为3.8 mm,决定系数达到0.91。

SBAS-InSAR技术融合CNN-LSTM模型的矿区开采沉陷监测与预测

针对传统矿区开采沉陷监测方法耗费人力财力和预测预警模型较少的问题,研究提出一种基于短基线集合成孔径雷达干涉测量(Small Baseline Subset-Interferometry Synthetic Aperture Radar,SBAS-InSAR)技术和卷积神经网络(Convolutional Neural Networks,CNN)与长短期记忆网络(Long Short-Term Memory,LSTM)相结合的矿区开采沉陷监测预测方法。首先,利用SBAS-InSAR技术对建新煤矿进行矿区开采沉陷监测,获取了该矿区的年平均沉降速率和累计沉降值。用GNSS监测数据与SBAS-InSAR结果进行对比验证,其拟合效果较好。其次,在此基础上利用CNN-LSTM模型预测后6期沉降数据,其结果与CNN和LSTM预测结果进行对比。研究显示,CNN-LSTM模型的平均绝对误差(S_(MAE))和均方根误差(S_(RMSE))比单一的CNN和LSTM分别至少降低了44.8%和40.6%,其决定系数均高于98%。最后,进一步预测前6期和中6期沉降数据,验证了CNN-LSTM预测模型在时间上的一致性。因此,SBAS-InSAR融合CNN-LSTM模型在类似矿山开采沉陷监测和预测中有较好的应用前景。

2023年积石山Ms6.2级地震InSAR同震形变探测与断层滑动分布反演

2023年12月18日,甘肃省临夏回族自治州积石山县发生了Ms6.2级地震,地震灾害损失显著高于同震级地震,查明其同震形变场与断层滑动分布特征,有助于揭示此次地震灾害损失严重的原因。利用Sentinel-1卫星升、降轨雷达影像,采用合成孔径雷达干涉技术获取了本次地震的同震形变场,进而基于Okada弹性位错模型,确定了本次地震的震源参数,并基于分布式滑动模型反演了本次地震断层面上的滑动分布。结果表明,积石山地震为逆冲型地震,升、降轨同震形变场沿视线向最大形变量分别为7.1 cm和7.8 cm,断层最大滑动量为0.31 m,主要集中在地下0~8 km。与2019年6月17日四川长宁Ms 6.0级地震对比发现,长宁地震地表最大形变量为8 cm,断层最大滑动量为0.38 m,均大于积石山地震,推测积石山地震灾害损失严重的主要原因为震区居民房屋相对集中、房屋抗震性能差和黄土区场地地震波放大效应等。

基于时序InSAR技术的中贵天然气管道天水市段沿线滑坡隐患识别与形变分析

中贵天然气管道天水市段沿线地处黄土高原,地貌类型多样、构造活跃、复杂地质条件下发育多种具有隐蔽性和潜伏性的地质灾害,其中滑坡灾害对中贵天然气管道的安全危害极大,因此,对中贵天然气管道天水段沿线滑坡隐患进行有效识别与分析具有重要意义。采用Sentinel-1A升降轨卫星数据,基于时序合成孔径雷达干涉测量(synthetic aperture radar interferometry,InSAR)对中贵天然气天水段沿线滑坡隐患进行了解译识别、现场复核、发育特征与典型滑坡形变分析。研究区共识别17处滑坡隐患点,现场复核最终确定13处,其余4处为人类工程活动区,其中常沟村和磨峪沟村滑坡形变受降雨影响。统计发现研究区管道沿线滑坡多发生在坡度40°~50°,高差400~600 m,东北坡向和距管道100~150 m范围内,且岩土体强度较低,地层以第四系全新统冲洪积层(Q_(4)^(a1+p1))为主。结果表明:联合升降轨的时序InSAR技术可以有效识别管道沿线滑坡隐患,为油气管线的安全运营及今后油气管道选线提供重要的科学依据和参考。

基于GAMMA软件二次开发的地表形变监测系统

InSAR技术具有大范围、高精度、全天时、全天候的优点,已被广泛地应用于地表形变监测与高程获取。GAMMA软件具有支持数据源丰富、处理功能全面、结果精度可靠和效率高等优势,是当前国际上常用的InSAR数据处理软件。但该软件采用命令行交互,操作命令繁杂、对专业要求程度高,不易于非专业用户使用。针对以上问题,本文利用界面友好、开源的QGIS对GAMMA进行二次开发。结合GAMMA提供的py_gamma库和QGIS提供的py_QGIS API接口,采用模块化设计、算法并行优化等技术,开发了一套界面化InSAR形变监测系统,并以西藏沃达地区106景哨兵1号数据进行处理试验,从软件性能和处理结果方面证明了基于QGIS和GAMMA二次开发的可行性和可靠性。相关内容可为命令行用户界面软件转图形化用户界面提供参考。

基于FL型地基雷达干涉测量技术的桥梁变形监测研究

地基雷达干涉测量技术可以在非接触状态下监测特大尺寸线状建(构)筑物的整体运营状态,且测量精度达到亚毫米级,是桥梁健康监测的重要手段。文中采用该技术FL模式监测特大跨径桥梁,获取全天候24 h桥梁形变演化数据,研究外界因素变化对桥梁变形规律的影响。结果表明,在自然环境中,桥梁变形与温度、气压呈负相关,与风速呈正相关,温度变化约1℃发生最大变形为11.44 mm。试验数据表明,自然状态下,变形最大值发生在边跨,而不是发生在跨中,是理论设计与复杂环境共同作用的结果。实测数据显示,地基雷达干涉测量系统可快速检测桥梁结构弯曲模态,尤其适合于大跨度桥梁体系在灾害后的紧急状态监测。

SIMULATION AND ANALYSIS ON OBSERVATIONAL ERRORS OF CLOUD INTENSITY AND STRUCTURE WITH TRMM PR AND GROUND-BASED RADAR

In order to explain theoretically the observational biases of reflectivity and structure of precipitation systems by TRMM Precipitation Radar (TRMM PR) and ground-based radar,the effects of wavelengths,incident direction of radar waves and radar beam width on the reflectivity observation are simulated.The results show that the error due to the different wavelength and incident direction of radar wave is within 2.0 dB,TRMM PR can observe a larger reflectivity than ground-based radar in echo center.TRMM PR smoothes the cloud structure,overestimates and underestimates reflectivity by 3-5 dB in strong and week echo areas,respectively.Beam width and long distance from TRMM PR to target cause it to overestimate the large echo area and area integrated rainfall amount,and to underestimate the averaged refleetivity.The theoretical results above can only explain part of observational facts,meaning that the comparison of observation results between TRMM PR and ground-based radar is complicated,the attenuation of radar wave within precipitation area is the main factor to affect the observed result.

采煤沉陷区立体化调查监测技术体系研究

采煤沉陷带来耕地受损、生态破坏和地灾隐患等一系列问题,开展采煤沉陷区调查监测,识别沉陷范围、预测沉陷边界、监测沉陷速率,有利于及时掌握沉陷地类、沉陷机理和稳沉状态,为耕地保护、生态修复及灾害防治等工作提供参考。本文通过合成孔径雷达干涉测量(Synthetic Aperture Radar Interferometry,InSAR)、微动勘探实测数据和沉陷区地层相近区域的地面沉降光纤孔监测数据研究皖北某煤矿沉陷区边界和采空区,结果表明:以10 mm为界线,利用InSAR可以快速划定沉降边界;通过微动勘探实测数据发现1000 m深度煤系地层的横波速度比正常地层低25~50 m/s,InSAR结果也显示该区域出现了地面沉降;分析收集到相近区域12组光纤孔监测数据,77.7~200 m深度段地层为主要变形区域,压缩量约为163.77 mm,数据频率能达到分钟级。由此可见,基于InSAR、微动和光纤孔监测三种方法实现多源异构数据耦合,能够在“高潜水位、中厚煤层、千米以深”的采煤沉陷区达到“由表、及里、实时”的立体化调查监测效果,构建立体化调查监测技术体系存在可行性。