Static Ground Displacement for an Induced Earthquake Recorded on Broadband Seismometers

2020 ◽  
Vol 110 (5) ◽  
pp. 2216-2224
Author(s):  
Megan Zecevic ◽  
Thomas S. Eyre ◽  
David W. Eaton

ABSTRACT Using geodetic methods, significant static ground deformation has been observed for many large natural earthquakes. Some of the largest earthquakes induced by hydraulic-fracturing operations have been observed in the Western Canada Sedimentary Basin; however, because of the size and depths of these events, the associated static ground deformations have not yet been observed using traditional geodetic techniques. A seismic processing technique, developed for small volcano-seismic events, has the potential to resolve micrometer-scale static displacements using broadband seismic data. In this study, we test this processing method using vertical component broadband recordings of an Mw 4.1 event acquired at four nearby broadband seismometers. Estimated static displacements at the four stations are compared with the theoretical surface displacement field for a dislocation on a finite rectangular source within a homogeneous, elastic half-space. The theoretical displacements have the same polarities as the measured displacements across the seismic network and have similar amplitudes for three of the four stations. However, one station yielded unstable results, which shows that care must be taken when using this method. These results suggest that this processing method has potential for obtaining surface deformation for small to moderate-sized earthquakes using broadband data.

2019 ◽  
Vol 220 (1) ◽  
pp. 585-597 ◽  
Author(s):  
Maria C Araya ◽  
Juliet Biggs

SUMMARY Tectonic slivers form in the overriding plate in regions of oblique subduction. The inner boundaries of the sliver are often poorly defined and can consist of well-defined faults, rotating blocks or diffuse fault systems, which pass through or near the volcanic arc. The Guanacaste Volcanic Arc Sliver (GVAS) as defined by Montero et al., is a segment of the Central American Forearc Sliver, whose inner boundary is the ∼87-km-long Haciendas-Chiripa fault system (HCFS), which is located ∼10 km behind the volcanic arc and consists of strike slip faults and pull apart steps. We characterize the current ground motion on this boundary by combining earthquake locations and focal mechanisms of the 2016 Bijagua earthquake sequence, with the surface ground deformation obtained from Interferometric Synthetic Aperture Radar (InSAR) images from the ALOS-2 satellite. The coseismic stack of interferograms show ∼6 cm of displacement towards the line of sight of the satellite between the Caño Negro fault and the Upala fault, indicating uplift or SE horizontal surface displacement. The largest recorded earthquake of the sequence was Mw 5.4, and the observed deformation is one of the smallest earthquakes yet detected by InSAR in the Central American region. Forward and inverse models show the surface deformation can be partially explained by slip on a single fault, but it can be better explained by slip along two faults linked at depth. The best-fitting model consists of 0.33 m of right lateral slip on the Caño Negro fault and 0.35 m of reverse slip on the Upala fault, forming a positive flower structure. As no reverse seismicity was recorded, we infer the slip on the Upala fault occurred aseismically. Observations of the Bijagua earthquake sequence suggests the forearc sliver boundary is a complex and diffuse fault system. There are localized zones of transpression and transtension and areas where there is no surface expression suggesting the fault system is not yet mature. Although aseismic slip is common on subduction interfaces and mature strike-slip faults, this is the first study to document aseismic slip on a continental tectonic sliver boundary fault.


Author(s):  
M. Ali ◽  
M. I. Shahzad ◽  
M. Nazeer ◽  
J. H. Kazmi

Earthquake cause ground deformation in sedimented surface areas like Pasni and that is a hazard. Such earthquake induced ground displacements can seriously damage building structures. On 7 February 2017, an earthquake with 6.3 magnitudes strike near to Pasni. We have successfully distinguished widely spread ground displacements for the Pasni earthquake by using InSAR-based analysis with Sentinel-1 satellite C-band data. The maps of surface displacement field resulting from the earthquake are generated. Sentinel-1 Wide Swath data acquired from 9 December 2016 to 28 February 2017 was used to generate displacement map. The interferogram revealed the area of deformation. The comparison map of interferometric vertical displacement in different time period was treated as an evidence of deformation caused by earthquake. Profile graphs of interferogram were created to estimate the vertical displacement range and trend. Pasni lies in strong earthquake magnitude effected area. The major surface deformation areas are divided into different zones based on significance of deformation. The average displacement in Pasni is estimated about 250 mm. Maximum pasni area is uplifted by earthquake and maximum uplifting occurs was about 1200 mm. Some of areas was subsidized like the areas near to shoreline and maximum subsidence was estimated about 1500 mm. Pasni is facing many problems due to increasing sea water intrusion under prevailing climatic change where land deformation due to a strong earthquake can augment its vulnerability.


2019 ◽  
Vol 11 (17) ◽  
pp. 2027 ◽  
Author(s):  
Chengsheng Yang ◽  
Dongxiao Zhang ◽  
Chaoying Zhao ◽  
Bingquan Han ◽  
Ruiqi Sun ◽  
...  

Fluid extraction or injection into underground reservoirs may cause ground deformation, manifested as subsidence or uplift. Excessive deformation may threaten the infrastructure of the oilfield and its surroundings and may even induce earthquakes. Therefore, the monitoring of surface deformation caused by oil production activities is important for the safe production of oilfields and safety assessments of surrounding infrastructure. Karamay oilfield is one of the major oil and gas fields in China. In this study, we take the Karamay oilfield in Xinjiang as a case study to detect surface deformation caused by subsurface fluid injection. Sentinel-1A images of 32 ascending (Path 114) and 34 descending (Path 165) tracks spanning March 2017 to August 2018, were used to derive vertical and horizontal deformation over Karamay oilfield using the MSBAS-InSAR method. The observed two-dimensional deformation exhibited significant vertical and east-west deformation in this region. The maximum uplift and horizontal velocity was approximately 160 mm/year and 65 mm/year, respectively. We also modeled one of the typical deformation zones using a dislocation model in a homogenous elastic half-space.


2021 ◽  
Vol 13 (3) ◽  
pp. 406
Author(s):  
Mehdi Darvishi ◽  
Georgia Destouni ◽  
Saeid Aminjafari ◽  
Fernando Jaramillo

Changes in subsurface water resources might alter the surrounding ground by generating subsidence or uplift, depending on geological and hydrogeological site characteristics. Improved understanding of the relationships between surface water storage and ground deformation is important for design and maintenance of hydraulic facilities and ground stability. Here, we construct one of the longest series of Interferometric Synthetic Aperture Radar (InSAR) to date, over twenty-five years, to study the relationships between water level changes and ground surface deformation in the surroundings of Lake Mead, United States, and at the site of the Hoover Dam. We use the Small Baseline Subset (SBAS) and Permanent scatterer interferometry (PSI) techniques over 177 SAR data, encompassing different SAR sensors including ERS1/2, Envisat, ALOS (PALSAR), and Sentinel-1(S1). We perform a cross-sensor examination of the relationship between water level changes and ground displacement. We found a negative relationship between water level change and ground deformation around the reservoir that was consistent across all sensors. The negative relationship was evident from the long-term changes in water level and deformation occurring from 1995 to 2014, and also from the intra-annual oscillations of the later period, 2014 to 2019, both around the reservoir and at the dam. These results suggest an elastic response of the ground surface to changes in water storage in the reservoir, both at the dam site and around the reservoir. Our study illustrates how InSAR-derived ground deformations can be consistent in time across sensors, showing the potential of detecting longer time-series of ground deformation.


2018 ◽  
Vol 89 (6) ◽  
pp. 2214-2226
Author(s):  
Qingjun Meng ◽  
Sidao Ni ◽  
Aizhi Guo ◽  
Yong Zhou

ABSTRACT The Mw 7.8 Kaikōura earthquake on 13 November 2016 is one of the most complex events ever recorded, with surface rupture found on more than a dozen faults. Within about 10 minutes after the mainshock, an Mw 5.8 event occurred and caused an 8 cm static displacement at high‐rate Global Positioning System (GPS) station KAIK, which was not accounted for in previous mainshock studies. In this article, we focus on the Mw 5.8 aftershock including (1) relocating the hypocenter using the hypo2000 method, (2) conducting a grid search for its point‐source mechanism and centroid location using seismic waveforms at four nearby stations, (3) inverting finite‐fault models of this event based on grid‐searched fault mechanism, and (4) calculating the surface ground deformation and estimating the deformation in the line of sight (LoS) directions of the ascending and descending Advanced Land Observation Satellite‐2 (ALOS‐2). Although we are not able to resolve the ruptured fault of the Mw 5.8 aftershock because of limited observation data, we estimate that this event can generate 10–20 cm ground surface displacement and affect the ground displacement observed on the Interferometric Synthetic Aperture Radar (InSAR) data near the Kaikōura Peninsular.


2018 ◽  
Vol 10 (12) ◽  
pp. 1901 ◽  
Author(s):  
Emanuela Valerio ◽  
Pietro Tizzani ◽  
Eugenio Carminati ◽  
Carlo Doglioni ◽  
Susi Pepe ◽  
...  

We investigate the Mw 6.5 Norcia (Central Italy) earthquake by exploiting seismological data, DInSAR measurements, and a numerical modelling approach. In particular, we first retrieve the vertical component (uplift and subsidence) of the displacements affecting the hangingwall and the footwall blocks of the seismogenic faults identified, at depth, through the hypocenters distribution analysis. To do this, we combine the DInSAR measurements obtained from coseismic SAR data pairs collected by the ALOS-2 sensor from ascending and descending orbits. The achieved vertical deformation map displays three main deformation patterns: (i) a major subsidence that reaches the maximum value of about 98 cm near the epicentral zones nearby the town of Norcia; (ii) two smaller uplift lobes that affect both the hangingwall (reaching maximum values of about 14 cm) and the footwall blocks (reaching maximum values of about 10 cm). Starting from this evidence, we compute the rock volumes affected by uplift and subsidence phenomena, highlighting that those involved by the retrieved subsidence are characterized by significantly higher deformation values than those affected by uplift (about 14 times). In order to provide a possible interpretation of this volumetric asymmetry, we extend our analysis by applying a 2D numerical modelling approach based on the finite element method, implemented in a structural-mechanic framework, and exploiting the available geological and seismological data, and the ground deformation measurements retrieved from the multi-orbit ALOS-2 DInSAR analysis. In this case, we consider two different scenarios: the first one based on a single SW-dipping fault, the latter on a main SW-dipping fault and an antithetic zone. In this context, the model characterized by the occurrence of an antithetic zone presents the retrieved best fit coseismic surface deformation pattern. This result allows us to interpret the subsidence and uplift phenomena caused by the Mw 6.5 Norcia earthquake as the result of the gravitational sliding of the hangingwall along the main fault plane and the frictional force acting in the opposite direction, consistently with the double couple fault plane mechanism.


2020 ◽  
Vol 12 (7) ◽  
pp. 1189 ◽  
Author(s):  
Pietro Mastro ◽  
Carmine Serio ◽  
Guido Masiello ◽  
Antonio Pepe

This work presents an overview of the multiple aperture synthetic aperture radar interferometric (MAI) technique, which is primarily used to measure the along-track components of the Earth’s surface deformation, by investigating its capabilities and potential applications. Such a method is widely used to monitor the time evolution of ground surface changes in areas with large deformations (e.g., due to glaciers movements or seismic episodes), permitting one to discriminate the three-dimensional (up–down, east–west, north–south) components of the Earth’s surface displacements. The MAI technique relies on the spectral diversity (SD) method, which consists of splitting the azimuth (range) Synthetic Aperture RADAR (SAR) signal spectrum into separate sub-bands to get an estimate of the surface displacement along the azimuth (sensor line-of-sight (LOS)) direction. Moreover, the SD techniques are also used to correct the atmospheric phase screen (APS) artefacts (e.g., the ionospheric and water vapor phase distortion effects) that corrupt surface displacement time-series obtained by currently available multi-temporal InSAR (MT-InSAR) tools. More recently, the SD methods have also been exploited for the fine co-registration of SAR data acquired with the Terrain Observation with Progressive Scans (TOPS) mode. This work is primarily devoted to illustrating the underlying rationale and effectiveness of the MAI and SD techniques as well as their applications. In addition, we present an innovative method to combine complementary information of the ground deformation collected from multi-orbit/multi-track satellite observations. In particular, the presented technique complements the recently developed Minimum Acceleration combination (MinA) method with MAI-driven azimuthal ground deformation measurements to obtain the time-series of the 3-D components of the deformation in areas affected by large deformation episodes. Experimental results encompass several case studies. The validity and relevance of the presented approaches are clearly demonstrated in the context of geospatial analyses.


2019 ◽  
Vol 93 (12) ◽  
pp. 2651-2660 ◽  
Author(s):  
Sergey Samsonov

AbstractThe previously presented Multidimensional Small Baseline Subset (MSBAS-2D) technique computes two-dimensional (2D), east and vertical, ground deformation time series from two or more ascending and descending Differential Interferometric Synthetic Aperture Radar (DInSAR) data sets by assuming that the contribution of the north deformation component is negligible. DInSAR data sets can be acquired with different temporal and spatial resolutions, viewing geometries and wavelengths. The MSBAS-2D technique has previously been used for mapping deformation due to mining, urban development, carbon sequestration, permafrost aggradation and pingo growth, and volcanic activities. In the case of glacier ice flow, the north deformation component is often too large to be negligible. Historically, the surface-parallel flow (SPF) constraint was used to compute the static three-dimensional (3D) velocity field at various glaciers. A novel MSBAS-3D technique has been developed for computing 3D deformation time series where the SPF constraint is utilized. This technique is used for mapping 3D deformation at the Barnes Ice Cap, Baffin Island, Nunavut, Canada, during January–March 2015, and the MSBAS-2D and MSBAS-3D solutions are compared. The MSBAS-3D technique can be used for studying glacier ice flow at other glaciers and other surface deformation processes with large north deformation component, such as landslides. The software implementation of MSBAS-3D technique can be downloaded from http://insar.ca/.


1990 ◽  
Vol 12 (1) ◽  
pp. 58-70 ◽  
Author(s):  
Lin Xin Yao ◽  
James A. Zagzebski ◽  
Ernest L. Madsen

In previous work, we demonstrated that accurate backscatter coefficient measurements are obtained with a data reduction method that explicitly accounts for experimental factors involved in recording echo data. An alternative, relative processing method for determining the backscatter coefficient and the attenuation coefficient is presented here. This method involves comparison of echo data from a sample with data recorded from a reference phantom whose backscatter and attenuation coefficients are known. A time domain processing technique is used to extract depth and frequency dependent signal ratios for the sample and the reference phantom. The attenuation coefficient and backscatter coefficient of the sample are found from these ratios. The method is tested using tissue-mimicking phantoms with known scattering and attenuation properties.


Author(s):  
Kuo Ding ◽  
Hui Li

Over the past several years, a metal mine by block caving method has experienced a long-term and progressive surface deformation and fracturing, and then we start our investigation based on this background. The location of surface rupture was based on a series of mapping activities and the deformation data was collected by GPS from 2013 to 2016. In this paper, emphasis was put on the analysis of the fissures, deformation and stress of surface subsidence. Results reveal the diversity magnitude and structural features of surface deformation and ground fissures. In addition, the time dependent behavior is comprehended and the subsidence zone reflects different types of time-displacement curve – regressive phase, steady phase and progressive phase, all these achievements indicate the complexity and diversity of the subsidence zone. On the other hand, stress calculation which inspired from the mechanical model of the cracking of hole wall is carried out, it is meaningful to understand the relation between fracture features, displacement vectors and horizontal stress.


Sign in / Sign up

Export Citation Format

Share Document