Unrest at Domuyo Volcano, Argentina, detected by geophysical and geodetic data and morphometric analysis

New volcanic unrest has been detected in the Domuyo Volcanic Center (DVC), to the east of the Andes Southern Volcanic Zone in Argentina. To better understand this activity, we investigated new seismic monitoring data, gravimetric and magnetic campaign data, and interferometric synthetic aperture radar (InSAR) deformation maps, and we derived an image of the magma plumbing system and the likely source of the unrest episode. Seismic events recorded during 2017-2018 nucleate beneath the southwestern flank of the DVC. Ground deformation maps derived from InSAR processing of Sentinel-1 data exhibit an inflation area exceeding 300 km2, from 2014 to at least March 2018, which can be explained by an inflating sill model located 7 km deep. The Bouguer anomaly reveals a negative density contrast of ~35 km wavelength, which is spatially coincident with the InSAR pattern. Our 3D density modeling suggests a body approximately 4-6 km deep with a density contrast of -550 kg/m3. Therefore, the geophysical and geodetic data allow identification of the plumbing system that is subject to inflation at these shallow crustal depths. We compared the presence and dimensions of the inferred doming area to the drainage patterns of the area, which support long-established incremental uplift according to morphometric analysis. Future studies will allow us to investigate further whether the new unrest is hydrothermal or magmatic in origin.Fil: Astort, Ana. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Estudios Andinos "Don Pablo Groeber". Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Estudios Andinos "Don Pablo Groeber"; ArgentinaFil: Walter, Thomas R. German Research Centre for Geosciences; AlemaniaFil: Ruiz, Francisco. Universidad Nacional de San Juan. Facultad de Ciencias Exactas, Físicas y Naturales. Instituto Geofísico Sismológico Volponi; ArgentinaFil: Sagripanti, Lucía. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Estudios Andinos "Don Pablo Groeber". Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Estudios Andinos "Don Pablo Groeber"; ArgentinaFil: Nacif, Andres Antonio. Universidad Nacional de San Juan. Facultad de Ciencias Exactas, Físicas y Naturales. Instituto Geofísico Sismológico Volponi; ArgentinaFil: Acosta, Gemma. Universidad Nacional de San Juan. Facultad de Ciencias Exactas, Físicas y Naturales. Instituto Geofísico Sismológico Volponi; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Folguera Telichevsky, Andres. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Estudios Andinos "Don Pablo Groeber". Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Estudios Andinos "Don Pablo Groeber"; Argentin

The 2010 MW 6.8 Yushu (Qinghai, China) earthquake: constraints provided by InSAR and body wave seismology

By combining observations from satellite radar, body wave seismology and optical imagery, we have determined the fault segmentation and sequence of ruptures for the 2010 Mw 6.8 Yushu (China) earthquake. We have mapped the fault trace using displacements from SAR image matching, interferometric phase and coherence, and 2.5 m SPOT-5 satellite images. Modeling the event as an elastic dislocation with three segments fitted to the fault trace suggests that the southeast and northwest segments are near vertical, with the central segment dipping 70° to the southwest; slip occurs mainly in the upper 10 km, with a maximum slip of 1.5 m at a depth of 4 km on the southeastern segment. The maximum slip in the top 1 km (i.e., near surface) is up to 1.2 m, and inferred locations of significant surface rupture are consistent with displacements from SAR image matching and field observations. The radar interferograms show rupture over a distance of almost 80 km, much larger than initial seismological and field estimates of the length of the fault. Part of this difference can be attributed to slip on the northwestern segment of the fault being due to an Mw 6.1 aftershock two hours after the main event. The remaining difference can be explained by a non-uniform slip distribution with much of the moment release occurring at depths of less than 10 km. The rupture on the central and southeastern segments of the fault in the main shock propagated at a speed of 2.5 km/s southeastward toward the town of Yushu located at the end of this segment, accounting for the considerable building damage. Strain accumulation since the last earthquake on the fault segment beyond Yushu is equivalent to an Mw 6.5 earthquake

An Automatic Mosaicking Algorithm for the Generation of a Large-Scale Forest Height Map Using Spaceborne Repeat-Pass InSAR Correlation Magnitude

This paper describes an automatic mosaicking algorithm for creating large-scale mosaic maps of forest height. In contrast to existing mosaicking approaches through using SAR backscatter power and/or InSAR phase, this paper utilizes the forest height estimates that are inverted from spaceborne repeat-pass cross-pol InSAR correlation magnitude. By using repeat-pass InSAR correlation measurements that are dominated by temporal decorrelation, it has been shown that a simplified inversion approach can be utilized to create a height-sensitive measure over the whole interferometric scene, where two scene-wide fitting parameters are able to characterize the mean behavior of the random motion and dielectric changes of the volume scatterers within the scene. In order to combine these single-scene results into a mosaic, a matrix formulation is used with nonlinear least squares and observations in adjacent-scene overlap areas to create a self-consistent estimate of forest height over the larger region. This automated mosaicking method has the benefit of suppressing the global fitting error and, thus, mitigating the “wallpapering” problem in the manual mosaicking process. The algorithm is validated over the U.S. state of Maine by using InSAR correlation magnitude data from ALOS/PALSAR and comparing the inverted forest height with Laser Vegetation Imaging Sensor (LVIS) height and National Biomass and Carbon Dataset (NBCD) basal area weighted (BAW) height. This paper serves as a companion work to previously demonstrated results, the combination of which is meant to be an observational prototype for NASA’s DESDynI-R (now called NISAR) and JAXA’s ALOS-2 satellite missions

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

Long-term monitoring of geodynamic surface deformation using SAR interferometry

Thesis (Ph.D.) University of Alaska Fairbanks, 2014Synthetic Aperture Radar Interferometry (InSAR) is a powerful tool to measure surface deformation and is well suited for surveying active volcanoes using historical and existing satellites. However, the value and applicability of InSAR for geodynamic monitoring problems is limited by the influence of temporal decorrelation and electromagnetic path delay variations in the atmosphere, both of which reduce the sensitivity and accuracy of the technique. The aim of this PhD thesis research is: how to optimize the quantity and quality of deformation signals extracted from InSAR stacks that contain only a low number of images in order to facilitate volcano monitoring and the study of their geophysical signatures. In particular, the focus is on methods of mitigating atmospheric artifacts in interferograms by combining time-series InSAR techniques and external atmospheric delay maps derived by Numerical Weather Prediction (NWP) models. In the first chapter of the thesis, the potential of the NWP Weather Research & Forecasting (WRF) model for InSAR data correction has been studied extensively. Forecasted atmospheric delays derived from operational High Resolution Rapid Refresh for the Alaska region (HRRRAK) products have been compared to radiosonding measurements in the first chapter. The result suggests that the HRRR-AK operational products are a good data source for correcting atmospheric delays in spaceborne geodetic radar observations, if the geophysical signal to be observed is larger than 20 mm. In the second chapter, an advanced method for integrating NWP products into the time series InSAR workflow is developed. The efficiency of the algorithm is tested via simulated data experiments, which demonstrate the method outperforms other more conventional methods. In Chapter 3, a geophysical case study is performed by applying the developed algorithm to the active volcanoes of Unimak Island Alaska (Westdahl, Fisher and Shishaldin) for long term volcano deformation monitoring. The volcano source location at Westdahl is determined to be approx. 7 km below sea level and approx. 3.5 km north of the Westdahl peak. This study demonstrates that Fisher caldera has had continuous subsidence over more than 10 years and there is no evident deformation signal around Shishaldin peak.Chapter 1. Performance of the High Resolution Atmospheric Model HRRR-AK for Correcting Geodetic Observations from Spaceborne Radars -- Chapter 2. Robust atmospheric filtering of InSAR data based on numerical weather prediction models -- Chapter 3. Subtle motion long term monitoring of Unimak Island from 2003 to 2010 by advanced time series SAR interferometry -- Chapter 4. Conclusion and future work

Satellite SAR Interferometry for Earth’s Crust Deformation Monitoring and Geological Phenomena Analysis

Synthetic aperture radar interferometry (InSAR) and the related processing techniques provide a unique tool for the quantitative measurement of the Earth’s surface deformation associated with certain geophysical processes (such as volcanic eruptions, landslides and earthquakes), thus making possible long-term monitoring of surface deformation and analysis of relevant geodynamic phenomena. This chapter provides an application-oriented perspective on the spaceborne InSAR technology with emphasis on subsequent geophysical investigations. First, the fundamentals of radar interferometry and differential interferometry, as well as error sources, are briefly introduced. Emphasis is then placed on the realistic simulation of the underlying geophysics processes, thus offering an unfolded perspective on both analytical and numerical approaches for modeling deformation sources. Finally, various experimental investigations conducted by acquiring SAR multitemporal observations on areas subject to deformation processes of particular geological interest are presented and discussed

Mapping ground instability in areas of geotechnical infrastructure using satellite InSAR and Small UAV Surveying: a case study in Northern Ireland

Satellite Interferometric Synthetic Aperture Radar (InSAR), geological data and Small Unmanned Aerial Vehicle (SUAV) surveying was used to enhance our understanding of ground movement at five areas of interest in Northern Ireland. In total 68 ERS-1/2 images 1992–2000 were processed with the Small Baseline Subset (SBAS) InSAR technique to derive the baseline ground instability scenario of key areas of interest for five stakeholders: TransportNI, Northern Ireland Railways, Department for the Economy, Arup, and Belfast City Council. These stakeholders require monitoring of ground deformation across either their geotechnical infrastructure (i.e., embankments, cuttings, engineered fills and earth retaining structures) or assessment of subsidence risk as a result of abandoned mine workings, using the most efficient, cost-effective methods, with a view to minimising and managing risk to their businesses. The InSAR results provided an overview of the extent and magnitude of ground deformation for a 3000 km2 region, including the key sites of the disused salt mines in Carrickfergus, the Belfast–Bangor railway line, Throne Bend and Ligoniel Park in Belfast, Straidkilly and Garron Point along the Antrim Coast Road, plus other urbanised areas in and around Belfast. Tailored SUAV campaigns with a X8 airframe and generation of very high resolution ortho-photographs and a 3D surface model via the Structure from Motion (SfM) approach at Maiden Mount salt mine collapse in Carrickfergus in 2016 and 2017 also demonstrate the benefits of very high resolution surveying technologies to detect localised deformation and indicators of ground instabilit

Estimation of Forest Height Using Spaceborne Repeat-Pass L-Band InSAR Correlation Magnitude over the US State of Maine

This paper describes a novel, simple and efficient approach to estimate forest height over a wide region utilizing spaceborne repeat-pass InSAR correlation magnitude data at L-band. We start from a semi-empirical modification of the RVoG model that characterizes repeat-pass InSAR correlation with large temporal baselines (e.g., 46 days for ALOS) by taking account of the temporal change effect of dielectric fluctuation and random motion of scatterers. By assuming (1) the temporal change parameters and forest backscatter profile/extinction coefficient follow some mean behavior across each inteferogram; (2) there is minimal ground scattering contribution for HV-polarization; and (3) the vertical wavenumber is small, a simplified inversion approach is developed to link the observed HV-polarized InSAR correlation magnitude to forest height and validated using ALOS/PALSAR repeat-pass observations against LVIS lidar heights over the Howland Research Forest in central Maine, US (with RMSE \u3c 4 m at a resolution of 32 hectares). The model parameters derived from this supervised regression are used as the basis for propagating the estimates of forest height to available interferometric pairs for the entire state of Maine, thus creating a state-mosaic map of forest height. The present approach described here serves as an alternative and complementary tool for other PolInSAR inversion techniques when full-polarization data may not be available. This work is also meant to be an observational prototype for NASA’s DESDynI-R (now called NISAR) and JAXA’s ALOS-2 satellite missions

Error estimation in multitemporal InSAR deformation time series, with application to Lanzarote, Canary Islands

Interferometric Synthetic Aperture Radar (InSAR) is a reliable technique for measuring crustal deformation. However, despite its long application in geophysical problems, its error estimation has been largely overlooked. Currently, the largest problem with InSAR is still the atmospheric propagation errors, which is why multitemporal interferometric techniques have been successfully developed using a series of interferograms. However, none of the standard multitemporal interferometric techniques, namely PS or SB (Persistent Scatterers and Small Baselines, respectively) provide an estimate of their precision. Here, we present a method to compute reliable estimates of the precision of the deformation time series. We implement it for the SB multitemporal interferometric technique (a favorable technique for natural terrains, the most usual target of geophysical applications). We describe the method that uses a properly weighted scheme that allows us to compute estimates for all interferogram pixels, enhanced by a Montecarlo resampling technique that properly propagates the interferogram errors (variance-covariances) into the unknown parameters (estimated errors for the displacements). We apply the multitemporal error estimation method to Lanzarote Island (Canary Islands), where no active magmatic activity has been reported in the last decades. We detect deformation around Timanfaya volcano (lengthening of line-of-sight ∼ subsidence), where the last eruption in 1730–1736 occurred. Deformation closely follows the surface temperature anomalies indicating that magma crystallization (cooling and contraction) of the 300-year shallow magmatic body under Timanfaya volcano is still ongoing.Peer reviewe