The excessive extraction and recharge of groundwater lead to long-time seasonal land subsidence in Los Angeles, USA, and especially in the Santa Ana basin. The rate of land subsidence in the Santa Ana basin has been rising, which could pose a danger to infrastructure and human lives. However, the most recent research on land surface deformation in the area was conducted using the traditional parameter estimation method, resulting in little understanding of the regional spatiotemporal characteristics. The parametric method consists of a least square linear inversion, using the pre-defined mathematical geometric or geophysical theoretical models to describe groundwater deformation, and it requires precise external environmental variables and accurate geophysical parameters, which are more difficult to implement. In this study, multitemporal InSAR-derived deformation time series are analyzed by using 69 descending C-band Sentinel-1A SAR scenes acquired from 2015 to 2018. A method based on independent component analysis (ICA) is applied to characterize the spatial pattern and temporal evolution of land subsidence in the Santa Ana basin. The results reveal two different spatial and temporal deformation patterns in the basin. First, a widespread seasonal deformation is identified by the first component, related to annual seasonal groundwater level changes, and the overall deformation shows a concentrated spatial pattern. The second component captures a long-term signal with a large-scale spatial pattern. For quantitative assessment, the obtained deformation time series are compared with the GNSS data, validating an accuracy of millimeters. We further calculate the cross-correlation coefficient and the elastic skeletal storage coefficient from the ICA-derived seasonal deformation and groundwater level, which reveals that the deformation responds quickly (i.e., a lag of 8 days) to the change in groundwater and the Santa Ana aquifer retains almost the same elasticity for at least 15 years. Quantifying the spatiotemporal characteristics of the deformation in the Santa Ana basin can provide a reference for the monitoring and managing of groundwater.

The Xi’an region of China has been suffering from groundwater depletion, ground fissure hazards, and surface subsidence for a long time. Due to the complex tectonics and frequent human and natural activities, land deformation in the region is aggravated, posing a threat to infrastructure and human life. This study adopted the multi-orbit and multi-temporal InSAR technology to measure multi-dimensional displacements and time-series displacements in Xi’an City. Through the multi-dimensional deformation verification, it was found that the control of groundwater flow direction by ground fissures is the cause of horizontal deformation. On the contrary, the flow direction of groundwater from west to east was inferred using multi-dimensional deformation. Further analysis was performed by calculating the deformation gradient of the cumulative deformation to obtain differential land subsidence and angular distortions, and it was quantitatively determined that the threshold for the generation of ground fissures caused by differential subsidence is 1/500. Then, through the mutual verification of the time series data and the groundwater level, a positive correlation was obtained. However, due to the inconsistent geological conditions and soil layers at the monitoring positions of Well 2 and Well 3, the lag time was 64 days and 4 days, respectively. Finally, the relationship between the surface deformation and the groundwater in the sustained uplift areas was explored. The Well 1 groundwater-level data with a monitoring period of 22 years and the corresponding monitoring points’ time series data were modeled; it was concluded that, in the future, the groundwater level will continue to rise and surface deformation will mainly increase, without a slowing trend. Therefore, research on the impact of surface uplift on infrastructure should be strengthened. By quantifying the relationship between land subsidence, ground fissures, and the groundwater level in Xi’an, the results of this study provide a reference for groundwater monitoring and management.

Groundwater resources are crucial to socio-economic development and the ecosystem, and over-extraction can cause the groundwater level to drop, deplete reserves, and trigger geological hazards like land subsidence. The North China Plain (NCP) has experienced both subsidence and groundwater depletion due to over-extraction in the past 70 years. In this study, we used MT-InSAR technology and ascending C-band Sentinel-1 SAR data from 2017 to 2023 to study land deformation in the junction area of Shijiazhuang–Baoding–Cangzhou–Hengshui. We identified multiple subsidence funnels with a maximum rate exceeding −150 mm/year and a total deformation surpassing 600 mm. Seasonal decomposition methods accurately separated seasonal signals in the time-series deformation and groundwater level data. An exponential function model applied to long-term deformation showed no significant decrease in subsidence in severely affected areas. By modeling seasonal deformation and seasonal groundwater levels, we determined the elastic skeletal storage coefficients (S_ke) to be in the range of 1.02 × 10⁻³~6.53 × 10⁻³ in subsidence areas. We obtained the spatiotemporal evolution of the total groundwater storage (TGWS), irreversible ground storage (IGWS), and recoverable ground storage (RGWS). The TGWS and IGWS decreased annually while the RGWS increased, which is attributable to the implementation of the South-to-North Water Diversion Project (SNWDP) and the issuance of groundwater withdrawal policies in the NCP.

Monitoring the actual surface deformation caused by sub-surface fluid activities is crucial for comprehending the sub-surface fluid mutation process. Incorporation of InSAR observations in a prior volume strain model (VSM) can facilitate the estimation of sub-surface fluid volumetric changes and three-dimensional surface deformation. However, the model parameters, i.e. the thickness and depth of the sub-surface fluid reservoir, are assumed to be known a priori. Therefore, it is impossible to acquire the precise model parameters in practice. Therefore, we proposed a novel approach for adaptive estimation of the model parameters by utilizing the black hole strategy method integrated with the particle swarm optimization algorithm. First, the performance of this method was tested through simulation experiments. The results suggest that the proposed method can precisely determine the model parameters and three-dimensional(3-D) deformations of the sub-surface fluids. Finally, the proposed model was tested using actual data from the Shale gas area in Texas and the Kilauea Caldera volcano in Hawaii, USA. For long-term, slow small deformations or sudden large deformations caused by the sub-surface fluids, the proposed algorithm can obtain corresponding three-dimensional deformations, especially in the presence of incoherent regions, which can still obtain a complete three-dimensional deformation field.

Time-series PolSAR are capable for continuous change monitoring of natural resources and urban land-covers regardless of weather and lighting conditions. However, in the big SAR data era, the scarcity of labelled PolSAR samples poses new challenge to the traditional change detection methods. To reduce the dependence on labeled samples and ensure the efficiency of long time-series PolSAR interpretation, an unsupervised and pseudo label-based change detection method is proposed. First, the similarity maps of time-series PolSAR are gauged by three selected polarimetric distance measurements (PDM), which are suitable for PolSAR distribution characteristics and have the potential to reflect PolSAR changes. Second, the high-confidence changed pseudo-samples are selected based on the similarity maps, and the unchanged pseudo-samples are selected based on the non-similarity maps. Thirdly, the limited selected pseudo-samples (changed and unchanged) and multi-dimensional features are used to train the ConvLSTM network for change detection, and the input features include the T₃ coherence matrix elements of time-series PolSAR and the aforementioned PDMs. Finally, the change detection results based on pseudo-samples and the ConvLSTM network can be obtained, without additional manual labels. Adequate experiments are conducted on Radarsat-2, UAVSAR full-polarized, and Sentinel-1 dual-polarized datasets, achieving improved unsupervised change detection accuracy at 89.59–93.24%.

Permanent Scatterer Interferometric Synthetic Aperture Radar (PSInSAR), with millimeter-level accuracy and full-resolution capabilities, is essential for monitoring urban deformation. With the advancement of SAR sensors in spatial and temporal resolution and the expansion of wide-swath observation capabilities, the number of permanent scatterers (PSs) in high-density urban areas has surged exponentially. To address these computational and memory challenges in high-density urban PSInSAR processing, this paper proposes an efficient method for integrating regular grid partitioning and connected component constraints. First, adaptive dynamic regular grid partitioning was employed to divide monitoring areas into sub-blocks, balancing memory usage and computational efficiency. Second, a weighted least squares adjustment model using common PS points in overlapping regions eliminated systematic inter-sub-block biases, ensuring global consistency. A graph-based connected component constraint mechanism was introduced to resolve multi-component segmentation issues within sub-blocks to preserve discontinuous PS information. Experiments on TerraSAR-X data covering Fuzhou, China (590 km²), demonstrated that the method processed 1.4 × 10⁷ PS points under 32 GB memory constraints, where it achieved a 25-fold efficiency improvement over traditional global PSInSAR. The deformation rates and elevation residuals exhibited high consistency with conventional methods (correlation coefficient ≥ 0.98). This method effectively addresses the issues of memory overflow, connectivity loss between sub-blocks, and cumulative merging errors in large-scale PS networks. It provides an efficient solution for wide-area millimeter-scale deformation monitoring in high-density urban areas, supporting applications such as geohazard early warning and urban infrastructure safety assessment.