Abstract
Ground deformation can cause serious environmental issues such as infrastructure damage, ground compaction, and reducing the ground capacity to store water. Mashhad, as one of the largest ...and most populated cities in the Middle East, has been suffering from extreme subsidence. In the last decade, some researchers have been interested in measuring land subsidence rates in the Mashhad valley by InSAR techniques. However, most of those studies were based on inaccurate measurements introducing uncertainties in the resulting subsidence rates. These researches used a small number of EnviSat data with long perpendicular and inhomogeneous temporal baseline. This paper seeks to determine the subsidence rate in urban areas of Mashhad in recent years, the threat that was neglected by the city managers and decision-makers. For this purpose, the Persistent Scatterer InSAR technique was applied in the study area using two time-series of descending and ascending Sentinel-1A acquisitions between 2014 and 2017. The results demonstrated the maximum line-of-sight deformation rate of 14.6 cm/year and maximum vertical deformation (subsidence) rate about 19.1 cm/year which could have irreversible consequences. The results were assessed and validated using piezometric data, GPS stations, and geotechnical properties. This assessment confirms that the main reason for subsidence in the interested area is groundwater over-extraction. Also, investigation of geotechnical properties shows that thick fine-grained layers in the northwest of the city could strongly affect the results. At the end of this paper, a new simplified method was proposed to estimate specific storage in special cases to predict the subsidence rate.
Abstract
The Ardabil plain, with an approximate area of 1097.2 km
2
in northwestern Iran, has experienced land subsidence due to intensive groundwater withdrawal and long seasons of drought in recent ...years. Different techniques have been used to investigate and evaluate subsidence in this region including: Global Positioning Systems (GPS), Levelling, and Geotechnical methods. These methods are typically expensive, time-consuming, and identify only a small fraction of the areas prone to subsidence. This study employs an Interferometric Synthetic Aperture Radar (InSAR) technique to measure the long-term subsidence of the plain. An open-source SAR interferometry time series analysis package, LiCSBAS, that integrates with the automated Sentinel-1 InSAR processor (COMET-LiCSAR) is used to analyze Sentinel-1 satellite images from October 2014 to January 2021. Processing of Sentinel-1 images shows that the Ardabil plain has been facing rapid subsidence due to groundwater pumping and reduced rainfall, especially between May 2018 to January 2019. The maximum subsidence rate was 45 mm/yr, measured at the southeastern part of the plain. While providing significant advantages (less processing time and disk space) over other InSAR processing packages, implementation of the LiCSBAS processing package and its accuracy for land subsidence measurements at different scales needs further evaluation. This study provides a procedure for evaluating its efficiency and accuracy for land subsidence measurements by comparing its measurements with the results of the GMTSAR and geotechnical numerical modeling. The results of geotechnical numerical modeling showed land subsidence with an average annual rate of 38 mm between 2006 and 2020, which was close to measurements using the InSAR technique. Comparison of the subsidence measurements of the Ardabil plain using the LiCSBAS package with results obtained from other techniques shows that LiCSBAS is able to accurately detect land deformation at large scales (~ km). However, they may not be optimized for more local deformations such as infrastructure monitoring.
Space-borne Synthetic Aperture Radar (SAR) Interferometry (InSAR) is now a key geophysical tool for surface deformation studies. The European Commission’s Sentinel-1 Constellation began acquiring ...data systematically in late 2014. The data, which are free and open access, have global coverage at moderate resolution with a 6 or 12-day revisit, enabling researchers to investigate large-scale surface deformation systematically through time. However, full exploitation of the potential of Sentinel-1 requires specific processing approaches as well as the efficient use of modern computing and data storage facilities. Here we present Looking Into Continents from Space with Synthetic Aperture Radar (LiCSAR), an operational system built for large-scale interferometric processing of Sentinel-1 data. LiCSAR is designed to automatically produce geocoded wrapped and unwrapped interferograms and coherence estimates, for large regions, at 0.001° resolution (WGS-84 coordinate system). The products are continuously updated at a frequency depending on prioritised regions (monthly, weekly or live update strategy). The products are open and freely accessible and downloadable through an online portal. We describe the algorithms, processing, and storage solutions implemented in LiCSAR, and show several case studies that use LiCSAR products to measure tectonic and volcanic deformation. We aim to accelerate the uptake of InSAR data by researchers as well as non-expert users by mass producing interferograms and derived products.
A new algorithm for forest height estimation based on dual polarimetric interferometric SAR data is presented in this study. The main objective is to consider the efficiency of the dual-polarization ...data compared to the full polarimetric images with respect to forest height retrieval. Accordingly, the forest height estimation based on the random volume over the ground model is examined using a geometrical procedure named the three-stage method. An exhaustive search polarization optimization technique is also applied to improve the results by employing the efficiency of all the polarization bases based on the four-dimensional lexicographic PolInSAR vector. The repeat-pass experimental SAR (ESAR) images, which include both L- and P-band full polarimetric data, are employed for the accuracy assessment of the dual PolInSAR data and the newly proposed method for forest height estimation. The experimental results on the L-band PolInSAR data show the ability of the dual PolInSAR data for forest height estimation with an average root mean square error (RMSE) of 4.97 m against Lidar data based on the conventional three-stage method. Additionally, the proposed method results in an accuracy of 2.95 m for forest height estimation, indicating its high potential for tree height retrieval.
An InSAR‐GNSS Velocity Field for Iran Watson, Andrew R.; Elliott, John R.; Lazecký, Milan ...
Geophysical research letters,
28 May 2024, Letnik:
51, Številka:
10
Journal Article
Recenzirano
Odprti dostop
We present average ground‐surface velocities and strain rates for the 1.7 million km2 area of Iran, from the joint inversion of InSAR‐derived displacements and GNSS data. We generate interferograms ...from 7 years of Sentinel‐1 radar acquisitions, correct for tropospheric noise using the GACOS system, estimate average velocities using LiCSBAS time‐series analysis, tie this into a Eurasia‐fixed reference frame, and perform a decomposition to estimate East and Vertical velocities at 500 m spacing. Our InSAR‐GNSS velocity fields reveal predominantly diffuse crustal deformation, with localized interseismic strain accumulation along the North Tabriz, Main Kopet Dagh, Main Recent, Sharoud, and Doruneh faults. We observe signals associated with recent groundwater subsidence, co‐ and postseismic deformation, active salt diaprism, and sediment motion. We derive high‐resolution strain rate estimates on a country‐ and fault‐scale, and discuss the difficulties of mapping diffuse strain rates in areas with abundant non‐tectonic and anthropogenic signals.
Plain Language Summary
Across the entire country of Iran, the ongoing convergence of the Arabian and Eurasian tectonic plates at about two cm every year is deforming the Earth's crust, producing earthquakes in this process of continental collision. Accurate measurements of how the ground is moving today are critical to understanding both the county‐scale deformation, and the local‐scale earthquake hazard, where the crust is deforming quickly and building up strain that is likely to be released in future major earthquakes. We combine multiple series of satellite radar images with GPS velocities to estimate East and vertical ground motion across all of Iran, at a higher level of detail than previous GPS‐only studies have been able to capture. Our velocity fields show a complex mix of ground motion signals, from crustal deformation on a country‐scale, to rapid land subsidence caused by the extraction of groundwater from aquifers. Some of the major faults are clearly building up strain for future earthquakes, but other regions are deforming much more diffusely, making it difficult to ascertain the locations of any future seismic hazard and ground shaking.
Key Points
We generate high‐resolution East and vertical velocity fields for Iran using Sentinel‐1 InSAR and GNSS observations
Regional deformation is diffuse, with interseismic strain localized onto the Doruneh, Main Kopet Dagh, North Tabriz, and Sharoud faults
Iran contains a wealth of time‐varying, short‐wavelength signals associated with groundwater extraction, salt diaprism, and sediment motion
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.
The unexpected collapse of land surface due to subsidence is one of the most significant geohazards that threatens human life and infrastructure. Kabudrahang and Famenin are two Iranian plains ...experiencing several sinkholes due to the characteristics of the underground soil layers and extreme groundwater depletion. In this study, space-based Synthetic Aperture Radar images are used to investigate the ground displacement behavior to examine the feasibility of Sentinel-1 data in detecting precursory deformation proceeding before the sinkhole formation. The selected sinkhole occurred in August 2018 in the vicinity of Kerdabad village in Hamedan province with a 40 m diameter and depth of ~40 m. Time series of the European constellation Sentinel-1 data, spanning from January 2015 to August 2018, is analyzed, and the results revealed a 3 cm annual subsidence (–3 cm/year) along with the line-of-sight direction. Time-series analysis demonstrated that the driving mechanism of the sinkhole formation had a gradual process. Displacement of persistent scatterers (PSs) near the cave area had an acceleration by approaching the sinkhole formation date. In contrast, other areas that are far from the cave area show linear subsidence behavior over time. Additionally, the one-kilometer deformation profile over the cave area indicates a high subsidence rate precisely at the location where the sinkhole was formed later on 20 August 2018.
Land cover classification is one of the most important applications of polarimetric SAR images, especially in urban areas. There are numerous features that can be extracted from these images, hence ...feature selection plays an important role in PolSAR image classification. In this study, three main steps are used to address this task: (1) feature extraction in the form of three categories, namely original data features, decomposition features, and SAR discriminators; (2) feature selection in the framework of the single and multi-objective optimization; and (3) image classification using the best subset of features. In single objective methods, we employ genetic algorithms (GAs) and support vector machines (SVMs) or multi-layer perceptron (MLP) neural network in order to maximize classification accuracy. Then a new method is proposed to perform an efficient land cover classification of the San Francisco Bay urban area based on the multi-objective optimization approach. The objectives are to minimize the error of classification and the number of selected PolSAR parameters. The experimental results on Radarsat-2 fine-quad data show that the proposed method outperforms the single objective approaches tested against it, while saving computational complexity. Finally, we show that the our method has a better performance than the SVM with full set of features and the Wishart classifier which is based on the covariance matrix.
In recent years, subsidence is considered one of the major environmental challenges in semiarid and arid regions of Iran. The main reason of land subsidence in Iran is injudicious groundwater ...extraction for agricultural and industrial activities. This study focuses on monitoring the displacement level of the ground in the west of Tehran which is highly affected by subsidence of Shahriar and Tehran plains. To support this study, Sentinel-1 images in both ascending (28) and descending (35) paths were used. The maximum subsidence rates of 16 cm/year and 15 cm/year were estimated in the southwest of Tehran and Shahriar Plain, respectively. Three subsidence centers were discovered in Shahriar Plain which is located in an urban area. The maximum LOS velocity was calculated about 5 cm/year in Fath Highway. The maximum rate of subsidence was estimated approximately 20-70 mm/year in the vicinity of Qods City. In Shahriar Plain, accuracy of InSAR results was assessed using cross-comparison analysis of the ascending and descending paths. However, in the southwest of Tehran, the vertical velocity obtained from Sentinel-1 data is validated using the GPS observations. The validations confirmed the accuracy of the PSInSAR algorithm and the InSAR results.
This paper proposes a new method for forest height estimation using single-baseline single frequency polarimetric synthetic aperture radar interferometry (PolInSAR) data. The new algorithm estimates ...the forest height based on the random volume over the ground with a volume temporal decorrelation (RVoG+VTD) model. We approach the problem using a four-stage geometrical method without the need for any prior information. In order to decrease the number of unknown parameters in the RVoG+VTD model, the mean extinction coefficient is estimated in an independent procedure. In this respect, the suggested algorithm estimates the mean extinction coefficient as a function of a geometrical index based on the signal penetration in the volume layer. As a result, the proposed four-stage algorithm can be used for forest height estimation using the repeat pass PolInSAR data, affected by temporal decorrelation, without the need for any auxiliary data. The suggested algorithm was applied to the PolInSAR data of the European Space Agency (ESA), BioSAR 2007 campaign. For the performance analysis of the proposed approach, repeat pass experimental SAR (ESAR) L-band data, acquired over the Remningstorp test site in Southern Sweden, is employed. The experimental result shows that the four-stage method estimates the volume height with an average root mean square error (RMSE) of 2.47 m against LiDAR heights. It presents a significant improvement of forest height accuracy, i.e., 5.42 m, compared to the three-stage method result, which ignores the temporal decorrelation effect.