High quality InSAR data linked to seasonal change in hydraulic head for an agricultural area in the San Luis Valley, Colorado

In the San Luis Valley (SLV), Colorado legislation passed in 2004 requires that hydraulic head levels in the confined aquifer system stay within the range experienced in the years 1978–2000. While some measurements of hydraulic head exist, greater spatial and temporal sampling would be very valuable in understanding the behavior of the system. Interferometric synthetic aperture radar (InSAR) data provide fine spatial resolution measurements of Earth surface deformation, which can be related to hydraulic head change in the confined aquifer system. However, change in cm-scale crop structure with time leads to signal decorrelation, resulting in low quality data. Here we apply small baseline subset (SBAS) analysis to InSAR data collected from 1992 to 2001. We are able to show high levels of correlation, denoting high quality data, in areas between the center pivot irrigation circles, where the lack of water results in little surface vegetation. At three well locations we see a seasonal variation in the InSAR data that mimics the hydraulic head data. We use measured values of the elastic skeletal storage coefficient to estimate hydraulic head from the InSAR data. In general the magnitude of estimated and measured head agree to within the calculated error. However, the errors are unacceptably large due to both errors in the InSAR data and uncertainty in the measured value of the elastic skeletal storage coefficient. We conclude that InSAR is capturing the seasonal head variation, but that further research is required to obtain accurate hydraulic head estimates from the InSAR deformation measurements

State-of-the-art in studies of glacial isostatic adjustment for the British Isles: a literature review

Understanding the effects of glacial isostatic adjustment (GIA) of the British Isles is essential for the assessment of past and future sea-level trends. GIA has been extensively examined in the literature, employing different research methods and observational data types. Geological evidence from palaeo-shorelines and undisturbed sedimentary deposits has been used to reconstruct long-term relative sea-level change since the Last Glacial Maximum. This information derived from sea-level index points has been employed to inform empirical isobase models of the uplift in Scotland using trend surface and Gaussian trend surface analysis, as well as to calibrate more theory-driven GIA models that rely on Earth mantle rheology and ice sheet history. Furthermore, current short-term rates of GIA-induced crustal motion during the past few decades have been measured using different geodetic techniques, mainly continuous GPS (CGPS) and absolute gravimetry (AG). AG-measurements are generally employed to increase the accuracy of the CGPS estimates. Synthetic aperture radar interferometry (InSAR) looks promising as a relatively new technique to measure crustal uplift in the northern parts of Great Britain, where the GIA-induced vertical land deformation has its highest rate. This literature review provides an in-depth comparison and discussion of the development of these different research approaches

The worsening impacts of land reclamation assessed with Sentinel-1: The Rize (Turkey) test case

Massive amounts of land are being reclaimed to build airports, new cities, ports, and highways. Hundreds of kilometers are added each year, as coastlines are extended further out to the sea. In this paper, this urbanization approach is monitored by Persistent Scatterer Interferometry (PSI) technique with Sentinel-1 SAR data. The study aims to explore this technology in order to support local authorities to detect and evaluate subtle terrain displacements. For this purpose, a large 3-years Sentinel-1 stack composed by 92 images acquired between 07/01/2015 to 27/01/2018 is employed and stacking techniques are chosen to assess ground motion. The test site of this study, Rize, Turkey, has been declared at high risk of collapse and radical solutions such as the relocation of the entire city in another area are been taken into consideration. A media fact-checking approach, i.e. evaluating national and international press releases on the test site, is considered for the paper and this work presents many findings in different areas of the city. For instance, alerts are confirmed by inspecting several buildings reported by the press. Critical infrastructures are monitored as well. Portions of the harbor show high displacement rates, up to 1 cm/year, proving reported warnings. Rural villages belonging to the same municipality are also investigated and a mountainous village affected by landslide is considered in the study. Sentinel-1 is demonstrated to be a suitable system to detect and monitor small changes or buildings and infrastructures for these scenarios. These changes may be highly indicative of imminent damage which can lead to the loss of the structural integrity and subsequent failure of the structure in the long-term. In Rize, only a few known motion-critical structures are monitored daily with in-situ technologies. SAR interferometry can assist to save expensive inspection and monitoring services, especially in highly critical cases such as the one studied in this paper

An accurate method to correct atmospheric phase delay for InSAR with the ERA5 global atmospheric model

Differential SAR Interferometry (DInSAR) has proven its unprecedented ability and merits of monitoring ground deformation on a large scale with centimeter to millimeter accuracy. However, atmospheric artifacts due to spatial and temporal variations of the atmospheric state often affect the reliability and accuracy of its results. The commonly-known Atmospheric Phase Screen (APS) appears in the interferograms as ghost fringes not related to either topography or deformation. Atmospheric artifact mitigation remains one of the biggest challenges to be addressed within the DInSAR community. State-of-the-art research works have revealed that atmospheric artifacts can be partially compensated with empirical models, point-wise GPS zenith path delay, and numerical weather prediction models. In this study, we implement an accurate and realistic computing strategy using atmospheric reanalysis ERA5 data to estimate atmospheric artifacts. With this approach, the Line-of-Sight (LOS) path along the satellite trajectory and the monitored points is considered, rather than estimating it from the zenith path delay. Compared with the zenith delay-based method, the key advantage is that it can avoid errors caused by any anisotropic atmospheric phenomena. The accurate method is validated with Sentinel-1 data in three different test sites: Tenerife island (Spain), Almería (Spain), and Crete island (Greece). The effectiveness and performance of the method to remove APS from interferograms is evaluated in the three test sites showing a great improvement with respect to the zenith-based approach.Peer ReviewedPostprint (published version

Imaging multi-age construction settlement behaviour by advanced SAR interferometry

This paper focuses on the application of Advanced Satellite Synthetic Aperture Radar Interferometry (A-DInSAR) to subsidence-related issues, with particular reference to ground settlements due to external loads. Beyond the stratigraphic setting and the geotechnical properties of the subsoil, other relevant boundary conditions strongly influence the reliability of remotely sensed data for quantitative analyses and risk mitigation purposes. Because most of the Persistent Scatterer Interferometry (PSI) measurement points (Persistent Scatterers, PSs) lie on structures and infrastructures, the foundation type and the age of a construction are key factors for a proper interpretation of the time series of ground displacements. To exemplify a methodological approach to evaluate these issues, this paper refers to an analysis carried out in the coastal/deltaic plain west of Rome (Rome and Fiumicino municipalities) affected by subsidence and related damages to structures. This region is characterized by a complex geological setting (alternation of recent deposits with low and high compressibilities) and has been subjected to different urbanisation phases starting in the late 1800s, with a strong acceleration in the last few decades. The results of A-DInSAR analyses conducted from 1992 to 2015 have been interpreted in light of high-resolution geological/geotechnical models, the age of the construction, and the types of foundations of the buildings on which the PSs are located. Collection, interpretation, and processing of geo-thematic data were fundamental to obtain high-resolution models; change detection analyses of the land cover allowed us to classify structures/infrastructures in terms of the construction period. Additional information was collected to define the types of foundations, i.e., shallow versus deep foundations. As a result, we found that only by filtering and partitioning the A-DInSAR datasets on the basis of the above-mentioned boundary conditions can the related time series be considered a proxy of the consolidation process governing the subsidence related to external loads as confirmed by a comparison with results from a physically based back analysis based on Terzaghi's theory. Therefore, if properly managed, the A-DInSAR data represents a powerful tool for capturing the evolutionary stage of the process for a single building and has potential for forecasting the behaviour of the terrain-foundation-structure combination

Advanced InSAR atmospheric correction: MERIS/MODIS combination and stacked water vapour models

A major source of error for repeat-pass Interferometric Synthetic Aperture Radar (InSAR) is the phase delay in radio signal propagation through the atmosphere (especially the part due to tropospheric water vapour). Based on experience with the Global Positioning System (GPS)/Moderate Resolution Imaging Spectroradiometer (MODIS) integrated model and the Medium Resolution Imaging Spectrometer (MERIS) correction model, two new advanced InSAR water vapour correction models are demonstrated using both MERIS and MODIS data: (1) the MERIS/MODIS combination correction model (MMCC); and (2) the MERIS/MODIS stacked correction model (MMSC). The applications of both the MMCC and MMSC models to ENVISAT Advanced Synthetic Aperture Radar (ASAR) data over the Southern California Integrated GPS Network (SCIGN) region showed a significant reduction in water vapour effects on ASAR interferograms, with the root mean square (RMS) differences between GPS- and InSAR-derived range changes in the line-of-sight (LOS) direction decreasing from ,10mm before correction to ,5mm after correction, which is similar to the GPS/MODIS integrated and MERIS correction models. It is expected that these two advanced water vapour correction models can expand the application of MERIS and MODIS data for InSAR atmospheric correction. A simple but effective approach has been developed to destripe Terra MODIS images contaminated by radiometric calibration errors. Another two limiting factors on the MMCC and MMSC models have also been investigated in this paper: (1) the impact of the time difference between MODIS and SAR data; and (2) the frequency of cloud-free conditions at the global scale

InSAR reveals land deformation at Guangzhou and Foshan, China between 2011 and 2017 with COSMO-SkyMed data

Subsidence from groundwater extraction and underground tunnel excavation has been known for more than a decade in Guangzhou and Foshan, but past studies have only monitored the subsidence patterns as far as 2011 using InSAR. In this study, the deformation occurring during the most recent time-period between 2011 and 2017 has been measured using COSMO-SkyMed (CSK) to understand if changes in temporal and spatial patterns of subsidence rates occurred. Using InSAR time-series analysis (TS-InSAR), we found that significant surface displacement rates occurred in the study area varying from -35 mm/year (subsidence) to 10 mm/year (uplift). The 2011-2017 TS-InSAR results were compared to two separate TS-InSAR analyses (2011-2013, and 2013-2017). Our CSK TS-InSAR results are in broad agreement with previous ENVISAT results and levelling data, strengthening our conclusion that localised subsidence phenomena occurs at different locations in Guangzhou and Foshan. A comparison between temporal and spatial patterns of deformations from our TS-InSAR measurements and different land use types in Guangzhou shows that there is no clear relationship between them. Many local scale deformation zones have been identified related to different phenomena. The majority of deformations is related to excessive groundwater extraction for agricultural and industrial purposes but subsidence in areas of subway construction also occurred. Furthermore, a detailed analysis on the sinkhole collapse in early 2018 has been conducted, suggesting that surface loading may be a controlling factor of the subsidence, especially along the road and highway. Roads and highways with similar subsidence phenomenon are identified. Continuous monitoring of the deforming areas identified by our analysis is important to measure the magnitude and spatial pattern of the evolving deformations in order to minimise the risk and hazards of land subsidence

Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain

Subsidence related to multiple natural and human-induced processes affects an increasing number of areas worldwide. Although this phenomenon may involve surface deformation with 3D displacement components, negative vertical movement, either progressive or episodic, tends to dominate. Over the last decades, differential SAR interferometry (DInSAR) has become a very useful remote sensing tool for accurately measuring the spatial and temporal evolution of surface displacements over broad areas. This work discusses the main advantages and limitations of addressing active subsidence phenomena by means of DInSAR techniques from an end-user point of view. Special attention is paid to the spatial and temporal resolution, the precision of the measurements, and the usefulness of the data. The presented analysis is focused on DInSAR results exploitation of various ground subsidence phenomena (groundwater withdrawal, soil compaction, mining subsidence, evaporite dissolution subsidence, and volcanic deformation) with different displacement patterns in a selection of subsidence areas in Spain. Finally, a cost comparative study is performed for the different techniques applied.The different research areas included in this paper has been supported by the projects: CGL2005-05500-C02, CGL2008-06426-C01-01/BTE, AYA2 010-17448, IPT-2011-1234-310000, TEC-2008-06764, ACOMP/2010/082, AGL2009-08931/AGR, 2012GA-LC-036, 2003-03-4.3-I-014, CGL2006-05415, BEST-2011/225, CGL2010-16775, TEC2011-28201, 2012GA-LC-021 and the Banting Postdoctoral Fellowship to PJG

Imaging Land Subsidence Induced by Groundwater Extraction in Beijing (China) Using Satellite Radar Interferometry

Beijing is one of the most water-stressed cities in the world. Due to over-exploitation of groundwater, the Beijing region has been suffering from land subsidence since 1935. In this study, the Small Baseline InSAR technique has been employed to process Envisat ASAR images acquired between 2003 and 2010 and TerraSAR-X stripmap images collected from 2010 to 2011 to investigate land subsidence in the Beijing region. The maximum subsidence is seen in the eastern part of Beijing with a rate greater than 100 mm/year. Comparisons between InSAR and GPS derived subsidence rates show an RMS difference of 2.94 mm/year with a mean of 2.41 ± 1.84 mm/year. In addition, a high correlation was observed between InSAR subsidence rate maps derived from two different datasets (i.e., Envisat and TerraSAR-X). These demonstrate once again that InSAR is a powerful tool for monitoring land subsidence. InSAR derived subsidence rate maps have allowed for a comprehensive spatio-temporal analysis to identify the main triggering factors of land subsidence. Some interesting relationships in terms of land subsidence were found with groundwater level, active faults, accumulated soft soil thickness and different aquifer types. Furthermore, a relationship with the distances to pumping wells was also recognized in this work.This work was supported by the National Natural Science Foundation of China under Grant 41201419 and a China Scholarship Council (CSC) scholarship to Mi Chen. Roberto Tomás was supported by the Ministry of Education, Culture and Sport through the project PRX14/00100. Part of this work is also supported by the Spanish Ministry of Economy and Competitiveness and EU FEDER funds under projects TIN2014-55413-C2-2-P, by the UK Natural Environmental Research Council (NERC) through the LICS and IRNHiC projects (ref. NE/K010794/1 and NE/N012151/1, respectively), the ESA-MOST DRAGON-3 projects (ref. 10607 and 10665) and the Open Fund from the Key Laboratory of Earth Fissures Geological Disaster, Ministry of Land and Resources (Geological Survey of Jiangsu Province)