Search

Circular synthetic aperture radar (CSAR) possesses the capability of multi-angle observation, breaking through the geometric observation constraints of traditional strip SAR and holding the potential for three-dimensional imaging. Its sub-wavelength level of planar resolution, resulting from a long synthetic aperture, makes CSAR highly valuable in the field of high-precision mapping. However, the motion geometry of CSAR is more intricate compared to traditional strip SAR, demanding high precision from navigation systems. The accumulation of errors over the long synthetic aperture time cannot be overlooked. CSAR exhibits significant coupling between the range and azimuth directions, making traditional motion compensation methods based on linear SAR unsuitable for direct application in CSAR. The dynamic nature of flight, with its continuous changes in attitude, introduces a significant deformation error between the non-rigidly connected Inertial Measurement Unit (IMU) and the Global Positioning System (GPS). This deformation error makes it difficult to accurately obtain radar position information, resulting in imaging defocus. The research in this article uncovers a correlation between the deformation error and radial acceleration. Leveraging this insight, we propose utilizing radial acceleration to estimate residual motion errors. This paper delves into the analysis of Position and Orientation System (POS) errors, presenting a novel high-resolution CSAR motion compensation method based on airborne platform acceleration information. Once the system deformation parameters are calibrated using point targets, the deformation error can be directly calculated and compensated based on the acceleration information, ultimately resulting in the generation of a high-resolution image. In this paper, the effectiveness of the method is verified with airborne flight test data. This method can compensate for the deformation error and effectively improve the peak sidelobe ratio and integral sidelobe ratio of the target, thus improving image quality. The introduction of acceleration information provides new means and methods for high-resolution CSAR imaging. [ABSTRACT FROM AUTHOR]

Affected by geological factors, the scale of surface deformation in a hilly semi-desertification mining area varies. Meanwhile, there is certain dense vegetation on the ground, so it is difficult to construct a high-precision and complete image of a subsidence basin by using a single monitoring method, and hence the laws of the deformation and inversion of mining parameters cannot be known. Therefore, we firstly propose conducting collaborative monitoring by using InSAR (Interferometric Synthetic Aperture Radar), UAV (unmanned aerial vehicle), and 3DTLS (three-dimensional terrestrial laser scanning). The time-series complete surface subsidence basin is constructed by fusing heterogeneous data. In this paper, SBAS-InSAR (Small Baseline Subset) technology, which has the characteristics of reducing the time and space discorrelation, is used to obtain the small-scale deformation of the subsidence basin, oblique photogrammetry and 3D-TLS with strong penetrating power are used to obtain the anomaly and large-scale deformation, and the local polynomial interpolation based on the weight of heterogeneous data is used to construct a complete and high-precision subsidence basin. Compared with GNSS (Global Navigation Satellite System) monitoring data, the mean square errors of 1.442 m, 0.090 m, 0.072 m are obtained. The root mean square error of the high-precision image of the subsidence basin data is 0.040 m, accounting for 1.4% of the maximum subsidence value. The high-precision image of complete subsidence basin data can provide reliable support for the study of surface subsidence law and mining parameter inversion. [ABSTRACT FROM AUTHOR]

The continuous downward movement exhibited by the Lashagou landslide group in recent years poses a significant threat to the safety of both vehicles and pedestrians traversing the highway G310. By integrating geomorphological interpretation using multi-temporal optical images, interferometric synthetic aperture radar (InSAR) measurements, and continuous global navigation satellite system (GNSS) observations, this paper traced the formation period of the Lashagou landslide group, and explored its kinematic behavior under external drivers such as rainfall and snowmelt. The results indicate that the formation period can be specifically categorized into three periods: before, during, and after the construction of highway G310. The construction of highway G310 is the direct cause and prerequisite for the formation of the Lashagou landslide group, whereas summer precipitation and spring snowmelt are the external driving factors contributing to its continuous downward movement. Additionally, both the long-term seasonal downslope movement and transient acceleration events are strongly controlled by rainfall, and there is a time lag of approximately 1–2 days between the transient acceleration and heavy rainfall events. This study highlights the benefits of leveraging multi-source remote sensing data to investigate slow-moving landslides, which is advantageous for the implementation of effective control and engineering intervention to mitigate potential landslide disasters. [ABSTRACT FROM AUTHOR]

As the demand for precise navigation of aircraft increases in modern society, researching high-precision, high-autonomy navigation systems is both theoretically valuable and practically significant. Because the inertial navigation system (INS) has systematic and random errors, its output information diverges. Therefore, it is necessary to combine them with other navigation systems for real-time compensation and correction of these errors. The SAR matching positioning and navigation system uses synthetic aperture radar (SAR) image matching for platform positioning and compensates for the drift caused by errors in the inertial measurement unit (IMU). Images obtained by SAR are matched with digital landmark data, and the platform's position is calculated based on the SAR imaging geometry. However, SAR matching positioning faces challenges due to seasonal variations in SAR images, the need for typical landmarks for matching, and the lack of elevation information in two-dimensional SAR image matching. This paper proposes an airborne platform positioning method based on interferometric SAR (InSAR) interferogram matching. InSAR interferograms contain terrain elevation information, are less affected by seasonal changes, and provide higher positioning accuracy and robustness. By matching real-time InSAR-processed interferograms with simulated interferograms using a digital elevation model (DEM), three-dimensional position information about the matching points has been obtained. Subsequently, a three-dimensional positioning model for the platform has bene established using the unit line-of-sight vector decomposition method. In actual flight experiments using an FMCW Ku-band Interferometric SAR system, the proposed platform positioning framework demonstrated its ability to achieve precise positioning in the absence of signals from the global navigation satellite system (GNSS). [ABSTRACT FROM AUTHOR]

To explore the degree of constraint by Global Navigation Satellite System (GNSS) and Interferometric Synthetic Aperture Radar (InSAR) data on the Maduo earthquake within a layered earth model structure and to gain an insight into the seismogenic mechanism and the seismic risk in the surrounding area, this study employs D-InSAR technology to acquire the InSAR co-seismic deformation field of the Maduo earthquake on 22 May 2021. Utilizing both GNSS and InSAR data, the inversions constrained by single and joint data are conducted and compared to determine the co-seismic slip model and fault plane stress distribution of the Maduo earthquake. Additionally, this paper calculates the Coulomb stress changes induced by 14 M ≥ 7 strong earthquakes, considering co-seismic effects, post-seismic viscoelastic relaxation, and inter-seismic tectonic stress loading, on 19 fault segments within the Bayan Har block research area (96°E~106°E, 29°N~36°N) since 1900. The findings are as follows: (1) The maximum line-of-sight (LOS) deformation was approximately 0.9 m. The joint inversion rupture was primarily located in the Dongcao Along Lake section (~98.6°E), aligning with previous research outcomes. (2) The cumulative Coulomb stress at the Maduo earthquake's source location was −0.1333 MPa, while the inter-seismic stress loading amounted to 0.0745 MPa. The East Kunlun Fault, Maduo–Gande Fault, Ganzi–Yushu Fault, and Dari Fault C exhibited considerable stress loading, warranting attention due to heightened seismic risk. (3) Based on three different co-seismic slip models, the stress disturbance results caused by the Maduo earthquake to the surrounding area and fault did not differ significantly. After the earthquake, the seismogenic fault still has high seismic risk. [ABSTRACT FROM AUTHOR]

When airborne synthetic aperture radar (SAR) encounters long-time global navigation satellite system (GNSS) denial, the system cannot eliminate inertial navigation system (INS) accumulated drift. Platform positioning technology based on SAR image-matching is one of the important auxiliary navigation methods. This paper proposes a three-step positioning method for an airborne SAR platform, which can achieve the robust and high-precision estimation of platform position and velocity. Firstly, the motion model of the airborne SAR platform is established and a nonlinear overdetermined equation set of SAR Range-Doppler based on the ground-control points set obtained by SAR image-matching is constructed. Then, to overcome the ill-conditioned problem generated by the singular Jacobian matrix when solving the equations directly, a three-step robust and high-precision estimation of platform position and velocity is achieved through singular value decomposition and equation decoupling. Furthermore, the error transfer model of systematic and random platform positioning errors is derived. Finally, a set of semi-physical simulation experiments of airborne SAR is conducted to verify the effectiveness of the positioning method and the accuracy of the error model presented in this paper. [ABSTRACT FROM AUTHOR]

The hybrid interferometric synthetic aperture radar system is a combination of an along-track configuration and cross-track configuration. Based on linear ocean wave theory, an ocean wave inversion algorithm for a hybrid interferometric synthetic aperture radar system is proposed in this work. Using the interferometric synthetic aperture radar images acquired by the TerraSAR-X and TanDEM-X satellites and the interferometric synthetic aperture radar images acquired by an airborne interferometric radar altimeter with a certain degree of squint, the profile of ocean waves and the corresponding orbital velocities were retrieved by combining the new inversion algorithm with the cross-spectra. Meanwhile, key parameters of ocean waves, such as the significant wave height, significant wave orbital velocity, propagation direction, and wavelength of the dominant waves, were also extracted from the ocean wave spectra retrieved in this study. In order to evaluate the reliability of the new inversion algorithm, the retrieved significant wave heights were compared with those provided by the European Centre for Medium-Range Weather Forecasts and measured by a Global Navigation Satellite System buoy. The results showed that for the ocean waves retrieved from the spaceborne hybrid interferometric synthetic aperture radar images, the differences between the retrieved significant wave heights of the four subareas selected in this paper and those provided by European Centre for Medium-Range Weather Forecasts were approximately 0.01, –0.17, –0.55, and –0.37 m, respectively, and for the ocean waves retrieved from the airborne interferometric radar altimeter images, the differences between the retrieved significant wave heights corresponding to the M920 and M3120 images used in this paper and those measured by the Global Navigation Satellite System buoy were approximately –0.05 and –0.09, respectively. Therefore, the method proposed in this work could retrieve the ocean wave spectra well when the velocity bunching had a small influence; however, as the nonlinear influence of the velocity bunching increased, the difference between the significant wave heights retrieved using this method and provided by the European Centre for Medium-Range Weather Forecasts also increased. [ABSTRACT FROM AUTHOR]

The landslide monitoring method that uses the Sentinel 1A Interferogram Synthetic Aperture Radar (InSAR) through the Stanford Method for Persistent Scatterers (StaMPS) method is a complimentary but complex procedure without exact guidelines. Hence, this paper delivered a parametric test by examining the optimal settings of the Sentinel 1A Persistent Scatterer (PS)- and Small Baseline Subset (SBAS)-InSAR using the StaMPS compared to the Global Navigation Satellite Systems (GNSS) in landslide cases. This study first revealed parameters with the suggested values, such as amplitude dispersion used to describe amplitude stability, ranging from 0.47 to 0.48 for PS and equal to or more than 0.6 for SBAS in WuWanZai, Ali Mt. The study further examined the suggested values for other factors, including the following: unwrap grid size to re-estimate the size of the grid; unwrap gold n win as the Goldstein filtering window to reduce the noise; and unwrap time win as the smoothing window (in days) for estimating phase noise distributions between neighboring pixels. Furthermore, the study substantiated the recommended settings in the Woda and Shadong landslide cases with the GNSS, inferring that the SBAS has adequate feasibility in practical landslide monitoring. [ABSTRACT FROM AUTHOR]

Atmospheric delay is the primary error in ground-based synthetic aperture radar (GBSAR). The existing compensation methods include the external meteorological data correction method, the polynomial fitting method, and the persistent scatterers SAR interferometry (PSInSAR) calibration method. Combining the polynomial fitting and the persistent scatterers targets is the most popular method of GBSAR atmospheric delay compensation. However, the displacement component of the coherent target is always ignored in the atmospheric delay compensation, which is unpractical. A joint estimation method of ground displacement and atmospheric model parameters is developed in this paper. The displacement component is determined by the spatial and temporal features of the objects. The atmospheric delay component is regarded as a systematic error represented by a quadratic polynomial related to distance. The result is resolved by the least-square method. Compared to the existing method, the root-mean-square error (RMSE) of the proposed method had a significant improvement in the validation experiment. In the real in situ experiment, the time series obtained by the GBSAR had a similar trend to that acquired by the Global Positioning System (GPS) receiver. It is indicated that the proposed method can lead to a better deformation estimation by taking the deformation component into account in the atmospheric compensation. [ABSTRACT FROM AUTHOR]

Landslides along the Three Gorges Reservoir in China pose a threat to coastal residents and waterway safety. To reduce false positive misjudgments caused by a sudden local change in the landslide deformation curve, in this paper, we propose an effective method for predicting the deformation trend of reservoir bank landslides. We take reservoir bank landslides in the Wanzhou District of the Three Gorges Reservoir area as the research object. The Time-Series Interferometric Synthetic Aperture Radar (InSAR) method and 62 Sentinel-1A images from 2018 to 2022 were selected for landslide deformation monitoring, and the Hurst index was calculated to characterize the deformation trend. Furthermore, we propose a method for predicting the deformation trend based on the statistical distribution of deformation rates and the physical significance of the Hurst index. After the field survey and Global Positioning System (GPS) verification, the Time-Series InSAR results are shown to be reliable. We take the Sifangbei landslide as a representative case to analyze the validation results. It is found that the determined Sifangbei landslide deformation trend is consistent with the conclusions for the region. In addition, the deformation trend of a reservoir bank slope has obvious spatial and temporal differences. Changes in the reservoir water level and concentrated rainfall play roles similar to those of catalysts. The proposed method, involving the combination of Time-Series InSAR and the Hurst index, can effectively monitor deformation and predict the stability trend of reservoir bank landslides. The presented research results provide new ideas and solutions for landslide prevention and risk mitigation in reservoir areas. [ABSTRACT FROM AUTHOR]

Global navigation satellite system (GNSS) array antenna receivers are effective for suppressing wideband jamming. However, its anti-jamming performance decreases sharply when the number of wideband interference surpasses the number of array elements. Since a large number of jammers are often used in navigation countermeasures, it is crucial to keep array antenna receivers available in such conditions. Aiming at this issue, two main tasks were performed in this research and are presented in this paper: Firstly, the direction sensitivity of the sup-freedom anti-jamming performance is revealed and an anti-jamming method for array antenna receivers based on antenna rotation is proposed. Secondly, in order to determine the optimal rotation angle rapidly, a variable-step iteration algorithm based on gradient descent is proposed. Theoretical analysis and simulation show the effectiveness of the proposed anti-jamming method and the efficiency of the implementation algorithm. In a typical airborne scenario with a maximum azimuth difference of 90°, the anti-jamming ability of the proposed method improved by 21~26 dB and 5~10 dB for arrays adopting the PI (power inversion) and MVDR (minimum variance distortionless response) algorithms, respectively. The iterative efficiency improved by 78.35–99.63% in comparison with a traversal of 0.1° search resolution. The proposed method and algorithm are not limited to airborne scenarios and they might be influential to anti-jamming algorithms in the data domain. [ABSTRACT FROM AUTHOR]

With the rapid development of interferometric synthetic aperture radar (InSAR) measurement technology, its measurement accuracy requirements are increasing. Atmospheric delay errors must be corrected, especially in the case of crustal deformation monitoring, the 20% variation of tropospheric water vapor among InSAR pairs generally produces range from 10 cm to 14 cm deformation errors. Such errors can be of the same magnitude as the annual changes in crustal deformation, or even greater, masking crustal deformation information and seriously affecting the results of crustal deformation monitoring. Therefore, in order to obtain a more accurate InSAR atmospheric delay correction model, this paper calculated and integrated atmospheric delays that were estimated by different sources, including the 37 pressure levels of the fifth generation of the European Centre for Medium-Range Weather Forecasts (ECMWF)) numerical weather prediction model, ECMWF Reanalysis v5 (ERA5), and Global Navigation Satellite System (GNSS) measurement data from the crustal movement observation network of China, based on the variance component estimation (VCE) weighting method. The results showed that the integrated model, based on the VCE method, is better than the generic atmospheric correction online service (GACOS) model for InSAR measuring of crustal deformation. The precision in monitoring crustal deformations was improved by approximately 5 mm, the correlation coefficient of atmospheric delay errors and crustal deformations improved from 0.287 to 0.347, and accuracy improved by approximately 25%. However, the improvement in accuracy was limited because of system error decoherence that was induced by atmospheric noise caused by abundant vegetation or snow cover. Therefore, in order to achieve more accurate results, we recommend the adoption of the multi-source integrated atmospheric delay correction model, based on the VCE method, for InSAR high-precision measuring of crustal deformation and seismic activities. [ABSTRACT FROM AUTHOR]

Airborne aircrafts are dependent on the Global Navigation Satellite System (GNSS), which is susceptible to interference due to the satellite base-station and cooperative communication. Synthetic aperture radar altimeter (SARAL) provides the ability to measure the topographic terrain for matching with Digital Elevation Model (DEM) to achieve positioning without relying on GNSS. However, due to the near-vertical coupling in the delay-Doppler map (DDM), the similarity of DDMs of adjacent apertures is high, and the probability of successful matching is low. To this end, a novel neural network of terrain matching and aircraft positioning is proposed based on the airborne SARAL imagery. The model-driven terrain matching and aircraft positioning network (TMP-Net) is capable of realizing aircraft positioning by utilizing the real-time DDMs to match with the DEM-based DDM references, which are generated by a point-by-point coupling mechanism between the airborne routine and ground terrain DEM. Specifically, the training dataset is established by a numerical simulation method based on a semi-analytical model. Therefore, DEM-based DDM references can be generated by forward deduction when only regional DEM can be obtained. In addition to the model-based DDM generation, feature extraction, and similarity measurement, an aircraft positioning module is added. Three different positioning methods are designed to achieve the aircraft positioning, where three-point weighting exhibits the best performance in terms of positioning accuracy. Due to the fact that both the weighted triplet loss and softmax loss are employed in a cooperative manner, the matching accuracy can be improved and the positioning error can be reduced. Finally, both simulated and measured airborne datasets are used to validate the effectiveness of the proposed algorithm. Quantitative and qualitative evaluations show the superiority. [ABSTRACT FROM AUTHOR]

The GNSS-based bistatic SAR (GNSS-BSAR) system has emerged as a hotspot due to its low power consumption, nice concealment, and worldwide reach. However, the weak landing power density of the GNSS signal often necessitates prolonged integration to achieve an adequate signal-to-noise ratio (SNR). In this case, the effects of the receiver's time-frequency errors and atmospheric disturbances are significant and cannot be ignored. Therefore, we propose an ionospheric scintillation compensation-based imaging scheme for dual-channel GNSS-BSAR system. This strategy first extracts the reference phase, which contains the ionospheric phase scintillation and other errors. Subsequently, the azimuth phase of the target is divided into difference phase and reference phase. We apply the two-step phase compensation to eliminate Doppler phase errors, thus achieving precise focusing of SAR images. Three sets of experiments using the GPS L5 signal as the illuminator were conducted, coherently processing a 1.5 km by 0.8 km scene about 300 s. The comparative results show that the proposed method exhibited better focusing performance, avoiding the practical challenges encountered by traditional autofocus algorithms. Additionally, ionospheric phase scintillation extracted at different times of the day suggest diurnal variations, preliminary illustrating the potential of this technology for ionospheric-related studies. [ABSTRACT FROM AUTHOR]

The Cyprus Continuously Operating Natural Hazards Monitoring and Prevention System, abbreviated CyCLOPS, is a national strategic research infrastructure devoted to systematically studying geohazards in Cyprus and the Eastern Mediterranean, Middle East, and North Africa (EMMENA) region. Amongst others, CyCLOPS comprises six permanent sites, each housing a Tier-1 GNSS reference station co-located with two calibration-grade corner reflectors (CRs). The latter are strategically positioned to account for both the ascending and descending tracks of SAR satellite missions, including the ESA's Sentinel-1. As of June 2021, CyCLOPS has reached full operational capacity and plays a crucial role in monitoring the geodynamic regime within the southeastern Mediterranean area. Additionally, it actively tracks landslides occurring in the western part of Cyprus. Although CyCLOPS primarily concentrates on geohazard monitoring, its infrastructure is also configured to facilitate the radiometric calibration and geometric validation of Synthetic Aperture Radar (SAR) imagery. Consequently, this study evaluates the performance of Sentinel-1A SAR by exploiting the CyCLOPS network to determine key parameters including spatial resolution, sidelobe levels, Radar Cross-Section (RCS), Signal-to-Clutter Ratio (SCR), phase stability, and localization accuracy, through Point Target Analysis (PTA). The findings reveal the effectiveness of the CyCLOPS infrastructure to maintain high-quality radiometric parameters in SAR imagery, with consistent spatial resolution, controlled sidelobe levels, and reliable RCS and SCR values that closely adhere to theoretical expectations. With over two years of operational data, these findings enhance the understanding of Sentinel-1 SAR product quality and affirm CyCLOPS infrastructure's reliability. [ABSTRACT FROM AUTHOR]

Deformation along a subway rail network is related to the safe operation of the subway and the stability of construction facilities on the surface, making long-term deformation monitoring imperative. Long-term monitoring of surface deformation along the subway network and statistical analysis of the overall deformation situation are lacking in China. Therefore, targeting 35 Chinese cities whose subway mileage exceeds 50 km, we extracted their surface deformation along subway networks between 2018 and 2022, using spaceborne synthetic aperture radar (SAR) interferometry (InSAR) technology and Sentinel-1 satellite data. We verified the results with the continuous global navigation satellite system (GNSS) stations' data and found that the root mean square error (RMSE) of the InSAR results was 3.75 mm/year. Statistical analysis showed that ground subsidence along the subways was more prominent in Beijing, Tianjin, and other areas in the North China Plain, namely Kunming (which is dominated by karst landforms), as well as Shanghai, Guangzhou, Qingdao, and other coastal cities. In addition, an analysis revealed that the severity of surface subsidence correlated positively with a city's gross domestic product (GDP) with a Pearson correlation of 0.787, since the higher the GDP, the more frequent the construction and maintenance of subway, and the more commuters there are, which in turn exacerbates the disturbance to the surface. Additionally, the type of land cover also affects the ground deformation. Our findings provide a reference for constructing, operating, and maintaining the urban subway systems in China. [ABSTRACT FROM AUTHOR]

Interferometric Synthetic Aperture Radar (InSAR) is one of the most sufficient technologies to provide global digital elevation modeling (DEM). It unwraps the interferometric phase to provide observations for phase-to-height conversion. However, the phase gradient has great influence on the phase unwrapping quality, thereby affecting the topographic mapping accuracy. Multi-look processing can improve the reliability of phase unwrapping by reducing the noise phase gradient. Nevertheless, it reduces the spatial resolution while increasing the height phase gradient, thus lowering the reliability of height values. In this paper, we propose a multi-look iteration algorithm to suppress the noise and maintain the reliability of height values. First, we use a large number of looks (NL) to suppress the noise phase gradient and obtain a coarse DEM. Then taking this coarse DEM as a reference, we remove the topographic phase from the interferogram with a small NL, which will reduce height phase gradient and ensure the accuracy of phase unwrapping. Finally, we obtain DEM products with high precision and fine resolution. We validate the proposed algorithm using both simulated and real data, and obtain DEM products in a greatly undulated region using TanDEM-X data. Results show that the proposed method is capable of providing DEM with resolution of 4 m and accuracy of 1.73 m. [ABSTRACT FROM AUTHOR]

Interferometric synthetic aperture radar (InSAR) enables us to obtain precipitable water vapor (PWV) maps with high spatial resolution through the phase difference caused by refraction in the atmosphere. Although previous studies have evaluated the error level of InSARPWV observations, they validated it only with C-band InSARPWV observations. Since ionospheric disturbance seriously contaminates the InSAR phase in the case of the lower-frequency SAR system, it is necessary for a PWV error level evaluation correcting the ionospheric effect appropriately if we use lower-frequency SAR systems, such as the Advanced Land Observing Satellite-2 (ALOS-2). In this paper, we evaluated the error level of the L-band InSARPWV observation obtained from ALOS-2 data covering four areas in Japan. We compared the InSAR observations with global navigation satellite system (GNSS) atmospheric observations and estimated the L-band InSARPWV error value by utilizing the error propagation theory. As a result, the L-band InSARPWV absolute error reached 2.83 mm, which was comparable to traditional PWV observations. Moreover, we investigated the impacts of the seasonality, the interferometric coherence, and the height dependence on the PWV observation accuracy in InSAR. [ABSTRACT FROM AUTHOR]

The traditional leveling, total station, and global navigation satellite system (GNSS) and the new differential interferometric synthetic aperture radar (DInSAR) and terrestrial laser scanning (TLS) systems have their own advantages and limitations in the deformation monitoring of mining areas. It is difficult to obtain accurate deformation information only using single-source measurement data. In this study, we propose an LOS deformation correction method for DInSAR in mining areas by fusing ground data without control points. Based on free space data, small deformations at the edges of mining influence areas accurately obtained using DInSAR. By combining leveling/GNSS and TLS methods, it was possible to obtain large deformations in central areas without the need for control points located outside the mining influence range. For overcoming the non-uniform coordinates of the "space–ground" data and the limited overlap of the effective measurement ranges, the subsidence prediction model was employed to assist in its fusion. In addition, in LOS deformation correction, we retained the non-full cycle phase of DInSAR and replaced the full cycle phase with the one from the data fusion. Engineering experiments have shown that the correction results preserve the differences in the LOS deformations at the edge areas of the mine influence range, and they recover the lost LOS deformations at the center areas. Using the difference in the LOS deformation before and after correction as the verification indicator, the maximum absolute value of the errors after correction was 143 mm, which was approximately 6.4% of the maximum LOS deformation. In addition, there were still two errors that were large (−112 mm and −89 mm, respectively), and the absolute values of errors were not more than 75 mm. For all errors, the mean absolute value was 36 mm. Compared with 399 mm before correction, the error was reduced by 91%. This study provides technical support and theoretical reference for deformation monitoring and control in mining areas. [ABSTRACT FROM AUTHOR]

Spaceborne interferometric synthetic aperture radar (InSAR) techniques are important for landslide detection and monitoring; however, several limitations and uncertainties, such as the unique north–south flying direction and side-look radar observing geometry, currently limit the ability of InSAR to credibly detect landslides, especially those related to high and steep slopes. Here, we conducted experimental and statistical analysis on the feasibility of time-series InSAR monitoring for steep slopes using ascending and descending SAR images. First, the theoretical (TGNSS), practical (PGNSS), and terrain (Hterrain) (T-P-H) indices for sensitivity evaluations of the slope displacement monitoring results from time-series InSAR were proposed for slope monitoring. Subsequently, two experimental and statistical studies were conducted for the cases with and without Global Navigation Satellite System (GNSS) monitoring data. Our experimental results of two high and steep open-pit mines showed that the defined theoretical and practical sensitivity indices can quantitatively evaluate the feasibility of ascending and descending InSAR observations in steep-slope deformation monitoring with GNSS data, and the terrain sensitivity index can qualitatively evaluate the feasibility of landslide monitoring results from ascending and descending Sentinel-1 satellite data without GNSS data. We further demonstrate the generalizability of these proposed indices using four landslide cases with both public GNSS and InSAR monitoring data and 119 landslide cases with only InSAR monitoring data. The statistical results indicated that greater indices correlated with higher reliability of the monitoring results, suggesting that these novel indices have wide suitability and applicability. This study can help to improve the practice of slope deformation monitoring using spaceborne InSAR, especially for high and steep slopes. [ABSTRACT FROM AUTHOR]

Excessive exploitation of the groundwater has resulted in obvious three-dimensional (3D) deformation features on the surface of the Beijing Plain. This paper, by combining Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS) technologies, has obtained time-series information of the 3D surface deformation in the Beijing Plain, analyzing its spatial distribution characteristics. On this basis, the relationship between different controlling factors with the 3D deformation of the surface has been analyzed as well. The following results are obtained: (1) From 2013 to 2018, the land subsidence, which generally showed the trend of slowing down, was mainly concentrated in the eastern, northern, and southern regions of Beijing Plain, with multiple subsidence centers. (2) Under the International Terrestrial Reference Frame 2005 (ITRF2005), the horizontal direction of all GPS points in the plain is basically the same, with the dominant movement direction being NE112.5°~NE113.8°. Under the Eurasian reference frame, the horizontal movement rate of GPS points significantly decreases. The movement rate and direction of each point are not characteristic of overall trend activity. (3) The distribution and extent of the 3D surface deformation in the Beijing Plain are controlled by the basement structure. Part of the subsided area corresponds to a Quaternary depression formed at the junction of active faults disrupting the area. Similarly, the distribution of horizontal deformation in the E-W and N-S directions of the plain is controlled by the regional basement structure comprising major faults bounding horizontal deformation. (4) Groundwater exploitation is the main cause of the 3D surface deformation in the Beijing Plain. The groundwater funnels of the second and third confined aquifer are in suitable agreement with the land subsidence. The horizontal movement in the Beijing Plain is either directed toward the center of the groundwater or the land subsidence funnel, and the deformation is directed from areas with higher to areas with lower groundwater levels. [ABSTRACT FROM AUTHOR]

Global Navigation Satellite System Reflectometry Bistatic Synthetic Aperture Radar (GNSS-R BSAR) is becoming more and more important in remote sensing because of its low power, low mass, low cost, and real-time global coverage capability. The Back Projection Algorithm (BPA) was usually selected as the GNSS-R BSAR imaging algorithm because it can process echo signals of complex geometric configurations. However, the huge computational cost is a challenge for its application in GNSS-R BSAR. Graphics Processing Units (GPU) provides an efficient computing platform for GNSS-R BSAR processing. In this paper, a solution accelerating the BPA of GNSS-R BSAR using GPU is proposed to improve imaging efficiency, and a matching pre-processing program was proposed to synchronize direct and echo signals to improve imaging quality. To process hundreds of gigabytes of data collected by a long-time synthetic aperture in fixed station mode, a stream processing structure was used to process such a large amount of data to solve the problem of limited GPU memory. In the improvement of the imaging efficiency, the imaging task is divided into pre-processing and BPA, which are performed in the Central Processing Unit (CPU) and GPU, respectively, and a pixel-oriented parallel processing method in back projection is adopted to avoid memory access conflicts caused by excessive data volume. The improved BPA with the long synthetic aperture time is verified through the simulation of and experimenting on the GPS-L5 signal. The results show that the proposed accelerating solution is capable of taking approximately 128.04 s, which is 156 times lower than pure CPU framework for producing a size of 600 m × 600 m image with 1800 s synthetic aperture time; in addition, the same imaging quality with the existing processing solution can be retained. [ABSTRACT FROM AUTHOR]

In the last thirty years, Synthetic Aperture Radar interferometry (InSAR) and the Global Navigation Satellite System (GNSS) have become fundamental space geodetic techniques for mapping surface deformations due to tectonic movements. One major limiting factor to those techniques is the effect of the troposphere, as surface velocities are of the order of a few mm yr−1, and high accuracy (to mm level) is required. The troposphere introduces a path delay in the microwave signal, which, in the case of GNSS Precise Point Positioning (PPP), can nowadays be partly removed with the use of specialized mapping functions. Moreover, tropospheric stratification and short wavelength spatial turbulences produce an additive noise to the low amplitude ground deformations calculated by the (multitemporal) InSAR methodology. InSAR atmospheric phase delay corrections are much more challenging, as opposed to GNSS PPP, due to the single pass geometry and the gridded nature of the acquired data. Thus, the precise knowledge of the tropospheric parameters along the propagation medium is extremely useful for the estimation and correction of the atmospheric phase delay. In this context, the PaTrop experiment aims to maximize the potential of using a high-resolution Limited-Area Model for the calculation and removal of the tropospheric noise from InSAR data, by following a synergistic approach and integrating all the latest advances in the fields of remote sensing meteorology (GNSS and InSAR) and Numerical Weather Forecasting (WRF). In the first phase of the experiment, presented in the current paper, we investigate the extent to which a high-resolution 1 km WRF weather re-analysis can produce detailed tropospheric delay maps of the required accuracy, by coupling its output (in terms of Zenith Total Delay or ZTD) with the vertical delay component in GNSS measurements. The model is initially operated with varying parameterization, with GNSS measurements providing a benchmark of real atmospheric conditions. Subsequently, the final WRF daily re-analysis run covers an extended period of one year, based on the optimum model parameterization scheme demonstrated by the parametric analysis. The two datasets (predicted and observed) are compared and statistically evaluated, in order to investigate the extent to which meteorological parameters that affect ZTD can be simulated accurately by the model under different weather conditions. Results demonstrate a strong correlation between predicted and observed ZTDs at the 19 GNSS stations throughout the year (R ranges from 0.91 to 0.93), with an average mean bias (MB) of –19.2 mm, indicating that the model tends to slightly underestimate the tropospheric ZTD as compared to the GNSS derived values. With respect to the seasonal component, model performance is better during the autumn period (October–December), followed by spring (April–June). Setting the acceptable bias range at ±23 mm (equal to the amplitude of one Sentinel-1 C-band phase cycle when projected to the zenithal distance), it is demonstrated that the model produces satisfactory results, with a percentage of ZTD values within the bias margin ranging from 57% in summer to 63% in autumn. [ABSTRACT FROM AUTHOR]

This paper discusses a full interferometry processing chain based on dual-orbit Sentinel-1A and Sentinel-1B (S1) synthetic aperture radar data and a combination of open-source routines from the Sentinel Application Platform (SNAP), Stanford Method for Persistent Scatterers (StaMPS), and additional routines introduced by the authors. These are used to provide vertical and East-West horizontal velocity maps over a study area in the south-western sector of the Po Plain (Italy) where land subsidence is recognized. The processing of long time series of displacements from a cluster of continuous global navigation satellite system stations is used to provide a global reference frame for line-of-sight–projected velocities and to validate velocity maps after the decomposition analysis. We thus introduce the main theoretical aspects related to error propagation analysis for the proposed methodology and provide the level of uncertainty of the validation analysis at relevant points. The combined SNAP–StaMPS workflow is shown to be a reliable tool for S1 data processing. Based on the validation procedure, the workflow allows decomposed velocity maps to be obtained with an accuracy of 2 mm/yr with expected uncertainty levels lower than 2 mm/yr. Slant-oriented and decomposed velocity maps provide new insights into the ground deformation phenomena that affect the study area arising from a combination of natural and anthropogenic sources. [ABSTRACT FROM AUTHOR]

The paper presents a multi-physics investigation of the ground movements related to the cyclical and seasonal injection and withdrawal of natural gas in/from a depleted reservoir located in the Po Plain area, Italy. Interferometric Synthetic Aperture Radar (InSAr) data (from 2003) and Global Navigation Satellite System (GNSS) data (from 2008) provided a full and coherent panorama of almost two decades of ground movement in the monitored area (more extended than the field boundary). The analysis of the acquired millimetric-scale movements together with the detailed geological analysis, both at reservoir and at regional scale, represents the focal point for understanding the investigated phenomena. Based on this information, a fully integrated and multidisciplinary geological, fluid-flow and geomechanical numerical modeling approach was developed to reproduce the main geometrical and structural features of the involved formations together with the poromechanics processes induced by the storage operations. The main achievement of the adopted methodology is a deep knowledge of the system and the involved processes, which is mandatory for the safety of the urbanized areas and the effective management of the underground resources. [ABSTRACT FROM AUTHOR]

David Glacier and Drygalski Ice Tongue are massive glaciers in Victoria Land, Antarctica. The ice from the East Antarctic Ice Sheet is drained through the former, and then discharged into the western Ross Sea through the latter. David Drygalski is the largest outlet glacier in Northern Victoria Land, floating kilometers out to sea. The floating and grounded part of the David Glacier are the main focus of this article. During the XXI Italian Antarctic Expedition (2005–2006), within the framework of the National Antarctic Research Programme (PNRA), two GNSS stations were installed at different points: the first close to the grounding line of David Glacier, and the second approximately 40 km downstream of the first one. Simultaneous data logging was performed by both GNSS stations for 24 days. In the latest data processing, the kinematic PPP technique was adopted to evaluate the dominant diurnal components and the very small semi-diurnal variations in ice motion induced by the ocean tide and the mean ice flow rates of both GNSS stations. Comparison of the GNSS time series with predicted ocean tide calculated from harmonic coefficients of the nearest tide gauge stations, installed at Cape Roberts and Mario Zucchelli Station, highlight different local response of the glacier to ocean tide, with a minor amplitude of vertical motion at a point partially anchored at the bedrock close to the grounding line. During low tide, the velocity of the ice flow reaches its daily maximum, in accordance with the direction of seawater outflow from the fjord into the ocean, while the greatest daily tidal excursion generates an increase in the horizontal ice flow velocity. With the aim to extend the analysis in spatial terms, five COSMO-SkyMED Stripmap scenes were processed. The comparison of the co-registered offset tracking rates, obtained from SAR images, with the GNSS estimation shows good agreement. [ABSTRACT FROM AUTHOR]

Tropospheric delay error must be reduced during interferometric synthetic aperture radar (InSAR) measurement. Depending on different geographical environments, an appropriate correction method should be selected to improve the accuracy of InSAR deformation monitoring. In this study, surface deformation monitoring was conducted in a high mountain gorge region in Yunnan Province, China, using Sentinel-1A images of ascending and descending tracks. The tropospheric delay in the InSAR interferogram was corrected using the Linear, Generic Atmospheric Correction Online Service for InSAR (GACOS) and ERA-5 meteorological reanalysis data (ERA5) methods. The correction effect was evaluated by combining phase standard deviation, semi-variance function, elevation correlation, and global navigation satellite system (GNSS) deformation monitoring results. The mean value of the phase standard deviation (Aver) of the linear correction interferogram and the threshold value (sill) of the semi-variogram were reduced by –20.98% and –41%, respectively, while the accuracy of the InSAR deformation points near the GNSS site was increased by 58%. The results showed that the three methods reduced the tropospheric delay error of InSAR deformation monitoring by different degrees in low-latitude mountains and valleys. Linear correction was the best at alleviating the tropospheric delay, followed by GACOS, while ERA5 had poor correction stability. [ABSTRACT FROM AUTHOR]

The differential interferometric synthetic aperture radar (DInSAR) method is based on phase variation between the complex value of pixels of timely separated scenes in interferometric SAR pairs. This phase variation has five components: surface topography, curvature of planet's surface, terrain displacement, volume scatterers, and atmospheric propagation effects. The terrain displacement is the main product of the DInSAR method, while the last two effects are unpredictable and bring inaccuracy into the terrain displacement measurements. In this work, the propagation conditions in the troposphere and ionosphere were studied during two DInSAR measurements examining the Zagreb 22 March 2020 earthquake, with terrain raising of up to +3 cm at the epicenter. For the troposphere, the vertical profile of the modified refraction index, which incorporates local curvature change with height, was reconstructed using aerological balloon probing data. Ionospheric conditions were determined based on total electron content (TEC) calculated from the Croatian positioning system (CROPOS) and global navigation satellite system (GNSS) reference stations' measurements. One of the DInSAR measurements was conducted in unfavorable tropospheric refractive conditions, which resulted in an overall bias of −2 cm. The variability of propagation conditions indicates the need for examining the atmospheric propagation effects when calculating terrain displacements using the DInSAR method. The results of DInSAR indicate slight displacements, comparable with the amplitude of atmospheric variations, and should therefore be approached with caution. [ABSTRACT FROM AUTHOR]

Current studies of global navigation satellite systems (GNSS)-based bistatic synthetic aperture radar (GNSS-SAR) is focused on static objects on land. However, moving target imaging is also very significant for modern SAR systems. Imaging a moving target has two main problems. One is the unknown range cell migration; the other is the motion parameter estimation, such as the target's velocity. This paper proposes a moving target imaging formation algorithm for GNSS-SAR. First, an approximate bistatic range history is derived to describe the phase variation of the target signal along the azimuth time. Then, a keystone transform is employed to correct the range cell migration. To address the motion parameter estimation, a chirp rate estimation method based on short-time Fourier transform and random sample consensus is proposed with high processing efficiency and robust estimation errors in low signal-to-noise ratio scenes. The estimated chirp rate can calculate the target's velocity. Finally, azimuth compression derivation is performed to accomplish GNSS-SAR imaging. A maritime experimental campaign is conducted to validate the effectiveness of the proposed algorithm. The two cargo ships in the SAR images have good accordance with the ground truth in terms of the target-to-receiver vertical distances along the range and the ships' length along the cross-range. [ABSTRACT FROM AUTHOR]

The wettest 2018 long rainy season (March to May (MAM)) resulted in daily intensive rainfall events in East Africa that have seriously affected the environment and economy in many countries. Land subsidence is one of the environmental disasters that has occurred due to the long rainy season in Kenya for many years. However, it has received limited scientific attention. In this paper, we incorporate hydrological (soil moisture active–passive (SMAP) and loading models) and geodetic data (global positioning system (GPS) and interferometric synthetic aperture radar (InSAR)) to study hydrological changes and their associated subsidence potential in Kenya. Results show that widespread subsidence of more than 20 mm was associated with the MAM season in Kenya during 2018, based on SBAS InSAR measurements. The high values of land subsidence were well correlated with the areas of intense flooding during the MAM season. The widespread subsidence during the wet season has implications for the stability of the earth's surface during the season rather than creating the possibility of potential stresses along active faults. These stresses may trigger seismicity that is expected to pose risks to urban features. The results of the current study can help governmental authorities to adopt proper urban planning that avoids or minimizes the risks of land subsidence in the areas of sinkholes. [ABSTRACT FROM AUTHOR]

This paper deals with the integration of deformation rates derived from Synthetic Aperture Radar Interferometry (InSAR) and Global Navigation Satellite System (GNSS) data. The proposed approach relies on knowledge of the variance/covariance of both InSAR and GNSS measurements so that they may be combined accounting for the spectral properties of their errors, hence preserving all spatial frequencies of the deformation detected by the two techniques. The variance/covariance description of the output product is also provided. A performance analysis is carried out on realistic simulated scenarios in order to show the boundaries of the technique. The proposed approach is finally applied to real data. Five Sentinel-1A/B stacks acquired over two different areas of interest are processed and discussed. The first example is a merged deformation map of the northern part of the Netherlands for both ascending and descending geometries. The second example shows the deformation at the junction between the North and East Anatolian Fault using three consecutive descending stacks. [ABSTRACT FROM AUTHOR]

The ML 5.8 earthquake that hit the island of Crete on 27 September 2021 is analysed with InSAR (Interferometry from Synthetic Aperture Radar) and GNSS (Global Navigation Satellite System) data. The purpose of this work is to create a model with sufficient detail for the geophysical processes that take place in several kilometres below the earth's surface and improve our ability to observe active tectonic processes using geodetic and seismic data. InSAR coseismic displacements maps show negative values along the LOS of ~18 cm for the ascending orbit and ~20 cm for the descending one. Similarly, the GNSS data of three permanent stations were used in PPK (Post Processing Kinematic) mode to (i) estimate the coseismic shifts, highlighting the same range of values as the InSAR, (ii) model the deformation of the ground associated with the main shock, and (iii) validate InSAR results by combining GNSS and InSAR data. This allowed us to constrain the geometric characteristics of the seismogenic fault and the slip distribution on it. Our model, which stands on a joint inversion of the InSAR and GNSS data, highlights a major rupture surface striking 214°, dipping 50° NW and extending at depth from 2.5 km down to 12 km. The kinematics is almost dip-slip normal (rake −106°), while a maximum slip of ~1.0 m occurred at a depth of ca. 6 km. The crucial though indirect role of inherited tectonic structures affecting the seismogenic crustal volume is also discussed suggesting their influence on the surrounding stress field and their capacity to dynamically merge distinct fault segments. [ABSTRACT FROM AUTHOR]

The paper aims to detect the main changes that occurred in the area surrounding the Mario Zucchelli Station (MZS) through analysis of multi-temporal remote sensing integrated by geophysical measurements. Specific attention was directed at realizing an integrated geomorphological study of the Boulder Clay Glacier, a partially debris-covered glacier belonging to the Northern Foothills (Victoria Land, Antarctica). This area was recently chosen as the location for the construction of a new semi-permanent gravel runway for MZS logistical airfreight operations. Photogrammetric analysis was performed by comparing three historical aerial photogrammetric surveys (carried out in 1956, 1985, and 1993) and Very High Resolution (VHR) GeoEye-1 satellite stereo-image coverage acquired in 2012. The comparison of geo-referenced orthophoto-mosaics allowed the main changes occurring in some particular areas along the coast nearby MZS to be established. Concerning the study of the Boulder Clay Glacier, it has to be considered that glaciers and moraines are not steady-state systems by definition. Several remote sensing and geophysical investigations were carried out with the main aim of determining the general assessment of this glacier: Ground Penetrating Radar (GPR); Geodetic Global Positioning System (GPS) network; multi-temporal satellite Synthetic Aperture Radar (SAR) interferometry. The analysis of Boulder Clay Glacier moraine pointed out a deformation of less than 74 mm y−1 in a time span of 56 years, value that agrees with velocity and deformation data observed by GPS and InSAR methods. The presence of unexpected brine ponds at the ice/bedrock interface and the deformation pattern observed in the central part of the moraine has to be monitored and studied, especially under the long-term maintenance of the future runway. [ABSTRACT FROM AUTHOR]

In the last years, the advanced synthetic aperture radar (SAR) interferometry (InSAR) has proven its effectiveness in the assessment of ground gotion with millimetric accuracy. Its integrated use with traditional (in-situ) topographic height determination techniques, such as geometric leveling and Global Navigation Satellite System (GNSS), is consolidated in underground fluids extraction areas for detecting and monitoring land subsidence. Nevertheless, the lack of a specific standardized methodology does not allow for evaluating different results obtained from different types of analysis. Starting from the description of two independent estimations of land subsidence in the Agosta (Comacchio, Italy) area, where an environmental impact assessment procedure was carried out following a request for gas exploitation, this paper points out the need for a standardized methodology, focused on the in-situ calibration of InSAR data. This last purpose requires an adequately dense and homogeneous reference GNSS network. The in progress initiatives, at European and national level, aiming at providing a Copernicus Ground Motion service could offer the opportunity to structure a reliable and dedicated GNSS network, starting from the large amount of stations run by different institutions already existing in Italy. [ABSTRACT FROM AUTHOR]

The use of global navigation satellite system reflectometry (GNSS-R) measurements for classification of inundated wetlands is presented. With the launch of NASA's Cyclone Global Navigation Satellite System (CYGNSS) mission, space-borne GNSS-R measurements have become available over ocean and land. CYGNSS covers latitudes between ±38°, providing measurements over tropical ecosystems and benefiting new studies of wetland inundation dynamics. The GNSS-R signal over inundated wetlands is driven mainly by coherent scattering associated with the presence of surface water, producing strong forward scattering and a distinctive bistatic scattering signature. This paper presents a methodology used to classify inundation in tropical wetlands using observables derived from GNSS-R measurements and ancillary data. The methodology employs a multiple decision tree randomized (MDTR) algorithm for classification and wetland inundation maps derived from the phased-array L-band synthetic aperture radar (PALSAR-2) as reference for training and validation. The development of an innovative GNSS-R wetland classification methodology is aimed to advance mapping of global wetland distribution and dynamics, which is critical for improved estimates of natural methane production. The results obtained in this manuscript demonstrate the ability of GNSS-R signals to detect inundation under dense vegetation over the Pacaya-Samiria Natural Reserve, a tropical wetland complex located in the Peruvian Amazon. Classification results report an accuracy of 69% for regions of inundated vegetation, 87% for open water regions, and 99% for non-inundated areas. Misclassification of inundated vegetation, primarily as non-inundated area, is likely related to the combination of two factors: partial inundation within the GNSS-R scattering area, and signal attenuation from dense overstory vegetation, resulting in a low signal. [ABSTRACT FROM AUTHOR]

Los Angeles has experienced ground deformations during the past decades. These ground displacements can be destructive for infrastructure and can reduce the land capacity for groundwater storage. Therefore, this paper seeks to evaluate the existing ground displacement patterns along a new metro tunnel in Los Angeles, known as the Sepulveda Transit Corridor. The goal is to find the most crucial areas suffering from subsidence or uplift and to enhance the previous reports in this metropolitan area. For this purpose, we applied a Persistent Scatterer Interferometric Synthetic Aperture Radar using 29 Sentinel-1A acquisitions from June 2017 to May 2018 to estimate the deformation rate. The assessment procedure demonstrated a high rate of subsidence in the Inglewood field that is near the study area of the Sepulveda Transit Corridor with a maximum deformation rate of 30 mm/yr. Finally, data derived from in situ instruments as groundwater level variations, GPS observations, and soil properties were collected and analyzed to interpret the results. Investigation of geotechnical boreholes indicates layers of fine-grained soils in some parts of the area and this observation confirms the necessity of more detailed geotechnical investigations for future constructions in the region. Results of investigating line-of-sight displacement rates showed asymmetric subsidence along the corridor and hence we proposed a new framework to evaluate the asymmetric subsidence index that can help the designers and decision makers of the project to consider solutions to control the current subsidence. [ABSTRACT FROM AUTHOR]

In this paper, aiming at the limitation of persistence scatterers (PS) points selection, a new method for selecting PS points has been introduced based on the average coherence coefficient, amplitude dispersion index, estimated signal-to-noise ratio and displacement standard deviation of multiple threshold optimization. The stability and quality of this method are better than that of a single model. In addition, an atmospheric correction model has also been proposed to estimate the atmospheric effects on Ground-based synthetic aperture radar (GBSAR) observations. After comparing the monitoring results before and after correction, we clearly found that the results are in good agreement with the actual observations after applying the proposed atmospheric correction approach. [ABSTRACT FROM AUTHOR]

Several methods allow accurate measurement of terrain surface motions. Global navigation satellite systems (GNSSes) and interferometry with synthetic aperture radar (InSAR) stand out in terms of measurement accuracy among them. In principle, both methods make it possible to evaluate a three-dimensional vector of the motion of points on the terrain surface. In this work, we dealt with the evaluation of motions in the up–down (U–D) and east–west direction (E–W) over underground gas storage (UGS) from InSAR. One crucial step in breaking down PSInSAR line of sight (LOS) measurements to U–D and E–W components is getting time series derived from individual tracks to the same time frame. This is usually performed by interpolation, but we used an innovative approach: we analyzed individual time series using the Lomb–Scargle periodogram (LSP), which is suitable for periodic noisy and irregularly sampled data; we selected the most significant period, created LSP models, and used them instead of the original time series. Then, it was possible to derive time series values for any arbitrary time step. To validate the results, we installed one GNSS receiver in the Tvrdonice UGS test area to perform independent measurements. The results show a good agreement in the evaluation of motions by both methods. The correlation coefficient between horizontal components from both PSInSAR and GNSS was 0.95 in the case of the E–W component, with an RMSE of 1.75 mm; for U–D they were 0.78 and 2.35 mm, respectively. In addition to comparing the motions in the U–D and E–W directions, we also created a comparison by converting GNSS measurements to a line of sight of the Sentinel-1 satellite to evaluate the conformity of InSAR and GNSS measurements. Based on descending track, the correlation coefficient between LOS from both methods is, on average, 0.97, with an RMSE of 2.70 mm. [ABSTRACT FROM AUTHOR]

The fusion of global navigation satellite system (GNSS) and interferometric synthetic aperture radar (InSAR) deformation data can leverage the advantages of GNSS high temporal resolution and InSAR high spatial resolution, and obtain more abundant deformation data for constraints on geophysical structural and mechanical parameters. Existing studies seldom consider the spatial heterogeneity of largescale deformation data, which easily leads to obvious spatial aggregation of errors in the results of fusion. Here, we propose a novel multiresolution segmentation fusion (MRSF) method that uses a multiresolution segmentation algorithm to automatically classify the spatial heterogeneity of InSAR deformation data with similar deformation characteristics. We applied the MRSF method to the fusion of GNSS and InSAR deformation data covering the central valley aquifer system (CVAS) in southern California to verify its precision and robustness. Results show that the MRSF method can accurately reflect spatiotemporal evolution characteristics of displacement data and reliably estimate deformation for the times and locations of missing data. We then tested this method for geophysical parameter estimation by constructing three different sets of data, including dense GNSS sites, sparse GNSS sites, and sparse GNSS sites fused with InSAR data using MRSF, to invert the slip distribution of the Cascadia subduction zone. Results show that the inverted slip of the fused InSAR and GNSS data is comparable to that of the dense GNSS sites. Therefore, the MRSF method can obtain deformation results with high precision and high spatiotemporal resolution and effectively compensate for the lack of data caused by sparse GNSS sites during the geophysical inversion process. [ABSTRACT FROM AUTHOR]

A long-strip differential interferometric synthetic aperture radar (DInSAR) measurement based on multi-frame image mosaicking is currently the realizable approach to measure large-scale ground deformation. As the spatial range of the mosaicked images increases, the nonlinear variation of ground ocean tidal loading (OTL) displacements is more significant, and using plane fitting to remove the large-scale errors will produce large tidal displacement residuals in a region with a complex coastline. To conveniently evaluate the ground tidal effect on mosaic DInSAR interferograms along the west coast of the U.S., a three-dimensional ground OTL displacements grid is generated by integrating tidal constituents' estimation of the GPS reference station network and global/regional ocean tidal models. Meanwhile, a solid earth tide (SET) model based on IERS conventions is used to estimate the high-precision SET displacements. Experimental results show that the OTL and SET in a long-strip interferogram can reach 77.5 mm, which corresponds to a 19.3% displacement component. Furthermore, the traditional bilinear ramp fitting methods will cause 7.2~20.3 mm residual tidal displacement in the mosaicked interferograms, and the integrated tidal constituents displacements calculation method can accurately eliminate the tendency of tidal displacement in the long-strip interferograms. [ABSTRACT FROM AUTHOR]

This study proposes a new set of processing procedures based on the strain model and the Kalman filter (SM-Kalman) to obtain high-precision three-dimensional surface deformation time series from interferometric synthetic aperture radar (InSAR) and global navigation satellite system (GNSS) data. Implementing the Kalman filter requires the establishment of state and observation equations. In the time domain, the state equation is generated by fitting the pre-existing deformation time series based on a deformation model containing linear and seasonal terms. In the space domain, the observation equation is established with the assistance of the strain model to realize the spatial combination of InSAR and GNSS observation data at each moment. Benefiting from the application of the Kalman filter, InSAR and GNSS data at different moments can be synchronized. The time and measurement update steps are performed dynamically to generate a 3-D deformation time series with high precision and a high resolution in the temporal and spatial domains. Sentinel-1 SAR and GNSS datasets in the Los Angeles area are used to verify the effectiveness of the proposed method. The datasets include twenty-seven ascending track SAR images, thirty-four descending track SAR images and the daily time series of forty-eight GNSS stations from January 2016 to November 2018. The experimental result demonstrates that the proposed SM-Kalman method can produce high-precision deformation results at the millimeter level and provide two types of 3-D deformation time series with the same temporal resolution as InSAR or GNSS observations according to the needs of users. The new method achieves a high degree of temporal and spatial fusion of GNSS and InSAR data. [ABSTRACT FROM AUTHOR]

This work addresses a methodology based on the interferometric synthetic aperture radar (InSAR) applied to analyze and monitor ground-motion phenomena induced by underground mining activities in the Legnica-Glogow copper district, south-western Poland. The adopted technique employs an InSAR processing chain that exploits a stack of Sentinel-1 synthetic aperture radar (SAR) images using a small baseline multitemporal approach. Interferograms with small temporal baselines are first selected, then their network is optimized and reduced to eliminate noisy data, in order to mitigate the effect of decorrelation sources related to seasonal phenomena, i.e., snow and vegetation growth, and to the radar acquisition geometry. The atmospheric disturbance is mitigated using a spatio-temporal filter based on the nonequispaced fast Fourier transform. The estimated displacement maps and time series show the effect of both linear and impulsive ground motion and are validated against global navigation satellite system (GNSS) measurements. In this context, a significant threat to the built environment is represented by seismic tremors triggered by underground mining activities, which are analyzed using the proposed method to integrate the information gathered by in situ seismometer devices. [ABSTRACT FROM AUTHOR]

Despite the growing number of remote-sensing products from satellite sensors, mapping of the combined spatial distribution and temporal variability of inundation in tropical wetlands remains challenging. An emerging innovative approach is offered by Global Navigation Satellite System reflectometry (GNSS-R), a concept that takes advantage of GNSS-transmitting satellites and independent radar receivers to provide bistatic radar observations of Earth's surface with large-scale coverage. The objective of this paper is to assess the capability of spaceborne GNSS reflections to characterize surface inundation dynamics in a complex wetlands environment in the Peruvian Amazon with respect to current state-of-the-art methods. This study examines contemporaneous ALOS2 PALSAR-2 L-band imaging radar, CYGNSS GNSS reflections, and ground measurements to assess associated advantages and challenges to mapping inundation dynamics, particularly in regions under dense tropical forest canopies. Three derivatives of CYGNSS Delay-Doppler maps (1) peak signal-to-noise ratio (SNR), (2) leading edge slope, and (3) trailing edge slope, demonstrated statistically significant logarithmic relationships with estimated flooded area percentages determined from SAR, with SNR exhibiting the strongest association. Aggregated Delay-Doppler maps SNR time series data examined for inundated regions undetected by imaging radar suggests GNSS-R exhibits a potentially greater sensitivity to inundation state beneath dense forest canopies relative to SAR. Results demonstrate the capability for mapping extent and dynamic wetlands ecosystems in complex tropical landscapes, alone or in combination with other remote-sensing techniques such as those based on imaging radar, contributing to enhanced mapping of these regions. However, several aspects of GNSS-R observations such as noise level, spatial resolution, and signal coherence need to be further examined. [ABSTRACT FROM AUTHOR]

Motion error is one of the most serious problems in airborne synthetic aperture radar (SAR) data processing. For a smoothly distributing backscatter scene or a seriously speed-varying velocity platform, the autofocusing performances of conventional algorithms, e.g., map-drift (MD) or phase gradient autofocus (PGA) are limited by their estimators. In this paper, combining the trajectories measured by global position system (GPS) and inertial navigation system (INS), we propose a novel error compensation method for varying accelerated airborne SAR based on the best linear unbiased estimation (BLUE). The proposed compensating method is particularly intended for varying acceleration SAR or homogeneous backscatter scenes, the processing procedures and computational cost of which are much simpler and lower than those of MD and PGA algorithms. [ABSTRACT FROM AUTHOR]

Ice velocity is one of the products associated with the Ice Sheets Essential Climate Variable. This paper describes the intercomparison and validation of ice-velocity measurements carried out by several international research groups within the European Space Agency Greenland Ice Sheet Climate Change Initiative project, based on space-borne Synthetic Aperture Radar (SAR) data. The goal of this activity was to survey the best SAR-based measurement and error characterization approaches currently in practice. To this end, four experiments were carried out, related to different processing techniques and scenarios, namely differential SAR interferometry, multi aperture SAR interferometry and offset-tracking of incoherent as well as of partially-coherent data. For each task, participants were provided with common datasets covering areas located on the Greenland ice-sheet margin and asked to provide mean velocity maps, quality characterization and a description of processing algorithms and parameters. The results were then intercompared and validated against GPS data, revealing in several cases significant differences in terms of coverage and accuracy. The algorithmic steps and parameters influencing the coverage, accuracy and spatial resolution of the measurements are discussed in detail for each technique, as well as the consistency between quality parameters and validation results. This allows several recommendations to be formulated, in particular concerning procedures which can reduce the impact of analyst decisions, and which are often found to be the cause of sub-optimal algorithm performance. [ABSTRACT FROM AUTHOR]

The geometric accuracy of synthetic aperture radar (SAR) data is usually derived from level-1 products using accurately surveyed corner reflector positions. This paper introduces a novel approach that derives the range delay and azimuth shift from acquired SAR raw data (level-0 products). Therefore, the propagation path is completely retrieved from SAR pulse transmission up to the reception of the point target’s backscatter. The procedure includes simple pulse compression in range and azimuth instead of full SAR data processing. By applying this method, the geometric accuracy of ESA’s Sentinel-1 SAR satellites (Sentinel-1A and Sentinel-1B) is derived for each satellite overpass by using corner reflectors with precisely surveyed GPS positions. The results show that the azimuth bias of about 2 m found in level-1 products for Stripmap acquisitions is reduced to about 15 cm. This indicates an artificial bias arising from operational Sentinel-1 SAR data processing. The remaining range bias of about 1.0 m, observed in L0-products, is interpreted as the offset between the SAR antenna phase center and the spacecraft’s center of gravity. The relative pixel localization accuracy derived with the proposed method is about 12 cm for the evaluated acquisitions. Compared to the full processed level-1 SAR data products, this accuracy is similar in the range direction, but, for the azimuth direction, it is improved by about 50% with the proposed method. [ABSTRACT FROM AUTHOR]

The aim of this paper is to assess the potential use of data recorded by the Global Navigation Satellite System Reflectometry (GNSS-R) Cyclone Global Navigation Satellite System (CYGNSS) constellation to characterize desert surface roughness. The study is applied over the Sahara, the largest non-polar desert in the world. This is based on a spatio-temporal analysis of variations in Cyclone Global Navigation Satellite System (CYGNSS) data, expressed as changes in reflectivity (Γ). In general, the reflectivity of each type of land surface (reliefs, dunes, etc.) encountered at the studied site is found to have a high temporal stability. A grid of CYGNSS Γ measurements has been developed, at the relatively fine resolution of 0.03° × 0.03°, and the resulting map of average reflectivity, computed over a 2.5-year period, illustrates the potential of CYGNSS data for the characterization of the main types of desert land surface (dunes, reliefs, etc.). A discussion of the relationship between aerodynamic or geometric roughness and CYGNSS reflectivity is proposed. A high correlation is observed between these roughness parameters and reflectivity. The behaviors of the GNSS-R reflectivity and the Advanced Land Observing Satellite-2 (ALOS-2) Synthetic Aperture Radar (SAR) backscattering coefficient are compared and found to be strongly correlated. An aerodynamic roughness (Z0) map of the Sahara is proposed, using four distinct classes of terrain roughness. [ABSTRACT FROM AUTHOR]

Ground-based radar interferometry, which can be specifically classified as ground-based synthetic aperture radar (GB-SAR) and ground-based real aperture radar (GB-RAR), was applied to monitor the Liusha Peninsula landslide and Baishazhou Yangtze River Bridge. The GB-SAR technique enabled us to obtain the daily displacement evolution of the landslide, with a maximum cumulative displacement of 20 mm in the 13-day observation period. The virtual reality-based panoramic technology (VRP) was introduced to illustrate the displacement evolutions intuitively and facilitate the following web-based panoramic image browsing. We applied GB-RAR to extract the operational modes of the large bridge and compared them with the global positioning system (GPS) measurement. Through full-scale test and time-frequency result analysis from two totally different monitoring methods, this paper emphasized the 3-D display potentiality by combining the GB-SAR results with VRP, and focused on the detection of multi-order resonance frequencies, as well as the configure improvement of ground-based radars in bridge health monitoring. [ABSTRACT FROM AUTHOR]

Global Navigation Satellite System Reflectometry (GNSS-R) based Bi-static Synthetic Aperture Radar (BSAR) is becoming more and more important in remote sensing, given its low power, low mass, low cost, and real-time global coverage capability. Due to its complex configuration, the imaging for GNSS-R BSAR is usually based on the Back-Projection Algorithm (BPA), which is very time consuming. In this paper, an efficient and general imaging algorithm for GNSS-R BSAR is presented. A Two Step Range Cell Migration (TSRCM) correction is firstly applied. The first step roughly compensates the RCM and Doppler phase caused by the motion of the transmitter, which simplifies the SAR data into the quasi-mono-static case. The second step removes the residual RCM caused by the motion of the receiver using the modified frequency scaling algorithm. Then, a cubic phase perturbation operation is introduced to equalize the Doppler frequency modulation rate along the same range cell. Finally, azimuth phase compensation and geometric correction are completed to obtain the focused SAR image. A simulation and experiment are conducted to demonstrate the feasibility of the proposed algorithm, showing that the proposed algorithm is more efficient than the BPA, without causing significant degradation in imaging quality. [ABSTRACT FROM AUTHOR]

In this study, joint inversions of Synthetic Aperture Radar (SAR) and Global Position System (GPS) measurements are used to investigate the source parameters of four Mw > 5 events of the 2016-2017 Central Italy earthquake sequence. The results show that the four events are all associated with a normal fault striking northwest-southeast and dipping southwest. The observations, in all cases, are consistent with slip on a rupture plane, with strike in the range of 157˚ to 164˚ and dip in the range of 39˚ to 44˚ that penetrates the uppermost crust to a depth of 0 to 8 km. The primary characteristics of these four events are that the 24 August 2016 Mw 6.2 Amatrice earthquake had pronounced heterogeneity of the slip distribution marked by two main slip patches, the 26 October 2016 Mw 6.1 Visso earthquake had a concentrated slip at 3-6 km, and the predominant slip of the 30 October 2016 Mw 6.6 Norcia earthquake occurred on the fault with a peak magnitude of 2.5 m at a depth of 0-6 km, suggesting that the rupture may have reached the surface, and the 18 January 2017 Mw 5.7 Campotosto earthquake had a large area of sliding at depth 3-9 km. The positive static stress changes on the fault planes of the latter three events demonstrate that the 24 August 2016 Amatrice earthquake may have triggered a cascading failure of earthquakes along the complex normal fault system in Central Italy. [ABSTRACT FROM AUTHOR]