scholarly journals Geodetic investigation of triggered slip below Fiordland

2021 ◽  
Author(s):  
◽  
Danielle Lindsay

<p>Secretary Island, at the head of Doubtful Sound in Fiordland, has been seismically active in past 30 years, with earthquakes larger than M w 6.5: the 1989 Doubtful Sound, 1993 Secretary Island, and 2003 Fiordland earthquakes. These events were approximately coincident with the 17° bend in the strike of the young, obliquely-converging, and steeply dipping Puysegur Subduction Zone. This section of the plate interface also has a history of triggered slip: the 1989 earthquake is inferred to have triggered the 1993 earthquake and, further north at George Sound, triggered afterslip was reported following the 2009 Dusky Sound earthquake. We have used L-band (23.6 cm-wavelength) Synthetic Aperture Radar (SAR) data from the ALOS1 and ALOS2 satellites, and C-band (5.5 cm-wavelength) SAR data from Sentinel 1A/B satellites, to test the hypothesis that triggered slip also occurred in the vicinity of Secretary Island following the 2007 George Sound, 2009 Dusky Sound and 2016 Kaikōura earthquakes. SAR images were aligned, interfered, filtered, and unwrapped using GMTSAR processing tools. Long-wavelength ionosphere noise was removed by inverting for the best-fitting linear plane, and we assumed a linear function of height to remove short-wavelength atmospheric noise. Small Baseline Subset (SBAS) timeseries analysis indicated a localised deformation signal centred on Secretary Island following the Dusky Sound earthquake. A re-analysis was undertaken of the co- and post-seismic deformation caused by the Dusky Sound earthquake so that any surface deformation centred on Secretary Island could be isolated. Campaign and continuous Global Positioning System (GPS) data were simultaneously inverted with co- and post-seismic interferograms using a statistical Bayesian modelling approach to determine the optimal Dusky Sound earthquake source parameters. Limitations arising from orbital drift, the frequency of SAR acquisitions and the observation geometry hindered our ability to constrain the timing, magnitude and location of reactivated slip from a source similar to the 2003 Secretary Island earthquake. Our findings indicate that slip was not triggered following either the 2007 George Sound earthquake or 2016 Kaikōura earthquake. However, we cannot rule out triggered slip near Secretary Island following the 2009 Dusky Sound earthquake. Any such slip likely occurred on an area of c. 350 km² (c. 15 km updip of the Secretary Island epicentre) with an average slip of 1–3 m, producing motion away from the satellite of c. 25 mm at Secretary Island.</p>

2021 ◽  
Author(s):  
◽  
Danielle Lindsay

<p>Secretary Island, at the head of Doubtful Sound in Fiordland, has been seismically active in past 30 years, with earthquakes larger than M w 6.5: the 1989 Doubtful Sound, 1993 Secretary Island, and 2003 Fiordland earthquakes. These events were approximately coincident with the 17° bend in the strike of the young, obliquely-converging, and steeply dipping Puysegur Subduction Zone. This section of the plate interface also has a history of triggered slip: the 1989 earthquake is inferred to have triggered the 1993 earthquake and, further north at George Sound, triggered afterslip was reported following the 2009 Dusky Sound earthquake. We have used L-band (23.6 cm-wavelength) Synthetic Aperture Radar (SAR) data from the ALOS1 and ALOS2 satellites, and C-band (5.5 cm-wavelength) SAR data from Sentinel 1A/B satellites, to test the hypothesis that triggered slip also occurred in the vicinity of Secretary Island following the 2007 George Sound, 2009 Dusky Sound and 2016 Kaikōura earthquakes. SAR images were aligned, interfered, filtered, and unwrapped using GMTSAR processing tools. Long-wavelength ionosphere noise was removed by inverting for the best-fitting linear plane, and we assumed a linear function of height to remove short-wavelength atmospheric noise. Small Baseline Subset (SBAS) timeseries analysis indicated a localised deformation signal centred on Secretary Island following the Dusky Sound earthquake. A re-analysis was undertaken of the co- and post-seismic deformation caused by the Dusky Sound earthquake so that any surface deformation centred on Secretary Island could be isolated. Campaign and continuous Global Positioning System (GPS) data were simultaneously inverted with co- and post-seismic interferograms using a statistical Bayesian modelling approach to determine the optimal Dusky Sound earthquake source parameters. Limitations arising from orbital drift, the frequency of SAR acquisitions and the observation geometry hindered our ability to constrain the timing, magnitude and location of reactivated slip from a source similar to the 2003 Secretary Island earthquake. Our findings indicate that slip was not triggered following either the 2007 George Sound earthquake or 2016 Kaikōura earthquake. However, we cannot rule out triggered slip near Secretary Island following the 2009 Dusky Sound earthquake. Any such slip likely occurred on an area of c. 350 km² (c. 15 km updip of the Secretary Island epicentre) with an average slip of 1–3 m, producing motion away from the satellite of c. 25 mm at Secretary Island.</p>


Teknik ◽  
2019 ◽  
Vol 39 (2) ◽  
pp. 126
Author(s):  
Arliandy Pratama Arbad ◽  
Wataru Takeuchi ◽  
Yosuke Aoki ◽  
Achmad Ardy ◽  
Mutiara Jamilah

Penginderaan jauh kini memainkan peranan penting dalam pengamatan perilaku gunung api. Penelitian ini bertujuan untuk mengamati deformasi permukaan Gunung Bromo, yang terletak di Jawa bagian Timur, Indonesia, yang masuk dalam rangkaian sistem volkanik di Taman Nasional Bukit Tengger Semeru (TNBTS). Penggunaan algoritma SAR Interferometry (InSAR) yang disebut sebagai pendekatan Small Baseline Subset (SBAS) memungkinkan perancangan peta kecepatan deformasi rata-rata dan and peta time series displacement di wilayah kajian. Teknik SBAS yang biasa menghasilkan rangkaian observasi tahap interferometrik. Ini tercatat sebagai kombinasi linear dari nilai fase SAR  scene untuk setiap pixel secara tersendiri. Analisis yang dilakukan terutama berdasarkan 22 data SAR data yang diperoleh melalui sensor ALOS/PALSAR selama kurun waktu 2007–2011. Beberapa penelitian menunjukkan bahwa kemampuan analisis InSAR dalam menyelidiki siklus gunung api, terutama Gunung Bromo yang memiliki karakteristik erupsi stratovolcano dalam satu hingga lima tahun. Analisis hasil memperlihatkan adanya kemajuan dari kajian sebelumnya akan InSAR wilayah tersebut, yang lebih fokus  kepada deformasi yang berpengaruh kepada kaldera. Hal ini menunjukkan bahwa penelitian ini bisa diimplementasikan pada manajemen risiko atau manajemen infrastruktur


Author(s):  
Riccardo Lanari ◽  
Manuela Bonano ◽  
Sabatino Buonanno ◽  
Francesco Casu ◽  
Claudio De Luca ◽  
...  

&lt;p&gt;The Sentinel-1 constellation of the Copernicus Program already represents a big revolution within the Earth Observation (EO) scenario. This result is mainly due to the capability of this constellation to acquire huge volumes of SAR data all over the globe, with a wide spatial coverage, a short revisit time (12 or 6 days in the case of one or two operating satellites, respectively), and a free and open access data policy. In particular, the availability of such a large amount of SAR data acquired through the TOPS mode, characterized by a short &amp;#8220;orbital tube&amp;#8221; (with a 200m nominal diameter) and a specific design for ensuring differential SAR interferometry (DInSAR) applications, has opened the possibility to investigate Earth surface deformation phenomena at unprecedented spatial scale and with a high temporal rate.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;Among several advanced DInSAR algorithms, a widely used approach is the Small BAseline Subset (SBAS) technique, which has already proven its effectiveness to investigate surface displacements with centimeter- to millimeter-level accuracy in different scenarios. Moreover, a parallel algorithmic solution for the SBAS approach, referred to as Parallel Small BAseline Subset (P-SBAS), has been recently developed. This approach permits to generate, in an automatic and unsupervised way, advanced DInSAR products by taking full benefit from parallel computing architectures, such as cluster, grid and, above all, cloud computing infrastructures.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;In this work we present the results of a DInSAR experiment, based on the P-SBAS approach, carried out at the European scale. In particular, we exploited the entire available Sentinel-1 dataset collected through the TOPS acquisition mode between March 2015 and September 2018 from descending orbits over large part of Europe. Moreover, the overall analysis wasbcarried out by using the Copernicus Data and Information Access Services (DIAS) and, in particular, those provided by the ONDA DIAS platform, which was selected through a public tender. This activity, carried out as stress test of the EPOSAR service included in the Satellite Data Thematic Core Service of the EPOS infrastructure, permitted to investigate the DIAS capacity to operationally serve systematic and automatic DInSAR processing services, such as the one based on the P-SBAS approach.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;Our experiment was successfully completed, allowing the retrieval of the deformation time-series of the overall investigated area with the final products having the main characteristics summarized in the following:&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;ul&gt;&lt;li&gt;Exploited Sentinel-1 data: ~72.000&lt;/li&gt; &lt;li&gt;Covered Area: ~4.500.000 km&lt;sup&gt;2&lt;/sup&gt;&lt;/li&gt; &lt;li&gt;Coherent (multilook) SAR pixels: ~120.000.000&lt;/li&gt; &lt;li&gt;Final products pixel dimension: ~80 m&lt;/li&gt; &lt;li&gt;Time elapsed: ~6 months&lt;/li&gt; &lt;/ul&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;The presented discussion will highlight the main pros and cons of the exploited solution for such wide area DInSAR experiment. Moreover, the analysis of the achieved results will also show the high quality of the retrieved DInSAR results, that can be of interest for the Solid Earth scientific community, and the potentially positive impact of the presented solution for what concerns the future development of the European Ground Motion Service.&lt;/p&gt;&lt;p&gt;This work is supported by: the 2019-2021 IREA-CNR and Italian Civil Protection Department agreement; the H2020 EPOS-SP project (GA 871121); the I-AMICA (PONa3_00363) project; and the IREA-CNR/DGSUNMIG agreement.&lt;/p&gt;


2021 ◽  
Author(s):  
Steven J. Gibbons ◽  
Stefano Lorito ◽  
Marc de la Asunción ◽  
Manuela Volpe ◽  
Jacopo Selva ◽  
...  

&lt;p&gt;Advances in GPU-based High-Performance Computing (HPC) facilities, combined with improvements in GPU-optimized shallow water models for tsunami inundation, allow us to perform large numbers of numerical simulations of earthquake-generated tsunamis on high-resolution numerical grids. Large numbers of simulations are necessary to investigate the multi-dimensional parameter space that defines the tsunami hazard, including situations where the tsunami is generated outside major tectonic structures, where fault geometry is uncertain and can take widely different orientations. With over 1500 numerical simulations, we perform suites of systematic parameter searches to investigate the sensitivity of inundation at the towns of Catania and Siracusa on Sicily to changes both in the earthquake source parameters and in the specification of the Manning friction coefficient. The inundation is modelled using the GPU-based Tsunami-HySEA code on a system of nested topo-bathymetric grids with a finest spatial resolution of 10 meters. We consider tsunamigenesis by large earthquakes with uniform slip where the location, focal depth, fault dimensions and slip, together with the angles of strike, dip, and rake, are defined by the standard Okada parameters. We consider sources both close to the shore, in which significant co-seismic deformation occurs, and offshore, where co-seismic deformation is negligible. For the offshore earthquake sources, we see systematic and intuitive changes in the inundation with changes in strike, dip, rake, and depth. For the near-shore sources, the dependency is far more complicated and co-seismic deformation becomes significant in determining the inundation. The sensitivity studies provide clear guidelines as to the necessary resolution for source discretization for Probabilistic Tsunami Hazard Analysis, with a need for a far finer discretization of local sources than for more distant sources. For a small number of earthquake sources, we study systematically the inundation as a function of the Manning Friction Coefficient. The sensitivity of the inundation to this parameter varies greatly for different earthquake sources and topo-bathymetry at the coastline of interest. An understanding of all these dependencies is needed to better understand the consequences of tsunamigenic earthquake models with more complex geometries, and in quantifying the epistemic uncertainty in the tsunami hazard.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;This work is partially funded by the European Union&amp;#8217;s Horizon 2020 Research and Innovation Program under grant agreement No 823844 (ChEESE Center of Excellence, www.cheese-coe.eu). Computational resources made available through Sigma2/UNINETT on Saga at NTNU, Trondheim, Norway (in project nn5008k) and through PRACE on Marconi-100 at CINECA, Rome, Italy (through PRACE grant Pra21_5386/TsuHazAP).&lt;/p&gt;


2020 ◽  
Vol 14 (4) ◽  
pp. 405-412
Author(s):  
Endra Gunawan ◽  
Takuya Nishimura ◽  
Susilo Susilo ◽  
Sri Widiyantoro ◽  
Nanang T. Puspito ◽  
...  

AbstractOn 6 December 2016 at 22:03 UTC, a devastating magnitude 6-class strike-slip earthquake occurred along an unidentified and unmapped fault in Pidie Jaya, northern Sumatra. We analysed the possible fault using continuous Global Positioning System (GPS) observation available in the region. In our investigation, we searched for the fault source parameters of the north- and south-dipping left-lateral faults and the west- and east-dipping right-lateral faults. We identified that the fault responsible for the earthquake was located offshore, with a southwest-northeast direction. We also computed the Coulomb failure stress and compared the result with the distribution of the aftershocks. In this study, we demonstrated that the result of the geological field survey conducted soon after the mainshock was attributed to the secondary effects of ground shaking and near-surface deformation, and not surface faulting. The newly identified offshore fault proposed by this study calls for further investigation of the corresponding submarine morphological attributes in this particular region.


Sensors ◽  
2020 ◽  
Vol 20 (16) ◽  
pp. 4464
Author(s):  
Jing Wang ◽  
Chao Wang ◽  
Hong Zhang ◽  
Yixian Tang ◽  
Xuefei Zhang ◽  
...  

The dynamic changes of the thawing and freezing processes of the active layer cause seasonal subsidence and uplift over a large area on the Qinghai–Tibet Plateau due to ongoing climate warming. To analyze and investigate the seasonal freeze–thaw process of the active layer, we employ the new small baseline subset (NSBAS) technique based on a piecewise displacement model, including seasonal deformation, as well as linear and residual deformation trends, to retrieve the surface deformation of the Beiluhe basin. We collect 35 Sentinel-1 images with a 12 days revisit time and 9 TerraSAR-X images with less-than two month revisit time from 2018 to 2019 to analyze the type of the amplitude of seasonal oscillation of different ground targets on the Beiluhe basin in detail. The Sentinel-1 results show that the amplitude of seasonal deformation is between −62.50 mm and 11.50 mm, and the linear deformation rate ranges from −24.50 mm/yr to 5.00 mm/yr (2018–2019) in the study area. The deformation trends in the Qinghai–Tibet Railway (QTR) and Qinghai–Tibet Highway (QTH) regions are stable, ranging from −18.00 mm to 6 mm. The InSAR results of Sentinel-1 and TerraSAR-X data show that seasonal deformation trends are consistent, exhibiting good correlations 0.78 and 0.84, and the seasonal and linear deformation rates of different ground targets are clearly different on the Beiluhe basin. Additionally, there are different time lags between the maximum freezing uplift or thawing subsidence and the maximum or minimum temperature for the different ground target areas. The deformation values of the alpine meadow and floodplain areas are higher compared with the alpine desert and barren areas, and the time lags of the freezing and thawing periods based on the Sentinel-1 results are longest in the alpine desert area, that is, 86 days and 65 days, respectively. Our research has important reference significance for the seasonal dynamic monitoring of different types of seasonal deformation and the extensive investigations of permafrost in Qinghai Tibet Plateau.


2021 ◽  
Author(s):  
Itzhak Lior ◽  
Anthony Sladen ◽  
Diego Mercerat ◽  
Jean-Paul Ampuero ◽  
Diane Rivet ◽  
...  

&lt;p&gt;The use of Distributed Acoustic Sensing (DAS) presents unique advantages for earthquake monitoring compared with standard seismic networks: spatially dense measurements adapted for harsh environments and designed for remote operation. However, the ability to determine earthquake source parameters using DAS is yet to be fully established. In particular, resolving the magnitude and stress drop, is a fundamental objective for seismic monitoring and earthquake early warning. To apply existing methods for source parameter estimation to DAS signals, they must first be converted from strain to ground motions. This conversion can be achieved using the waves&amp;#8217; apparent phase velocity, which varies for different seismic phases ranging from fast body-waves to slow surface- and scattered-waves. To facilitate this conversion and improve its reliability, an algorithm for slowness determination is presented, based on the local slant-stack transform. This approach yields a unique slowness value at each time instance of a DAS time-series. The ability to convert strain-rate signals to ground accelerations is validated using simulated data and applied to several earthquakes recorded by dark fibers of three ocean-bottom telecommunication cables in the Mediterranean Sea. The conversion emphasizes fast body-waves compared to slow scattered-waves and ambient noise, and is robust even in the presence of correlated noise and varying wave propagation directions. Good agreement is found between source parameters determined using converted DAS waveforms and on-land seismometers for both P- and S-wave records. The demonstrated ability to resolve source parameters using P-waves on horizontal ocean-bottom fibers is key for the implementation of DAS based earthquake early warning, which will significantly improve hazard mitigation capabilities for offshore and tsunami earthquakes.&lt;/p&gt;


2021 ◽  
Author(s):  
Simone Cesca ◽  
Carla Valenzuela Malebrán ◽  
José Ángel López-Comino ◽  
Timothy Davis ◽  
Carlos Tassara ◽  
...  

&lt;p&gt; A complex seismic sequence took place in 2014 at the Juan Fern&amp;#225;ndez microplate, a small microplate located between Pacific, Nazca and Antarctica plates. Despite the remoteness of the study region and the lack of local data, we were able to resolve earthquake source parameters and to reconstruct the complex seismic sequence, by using modern waveform-based seismological techniques. The sequence started with an exceptional Mw 7.1-6.7 thrust &amp;#8211; strike slip earthquake doublet, the first subevent being the largest earthquake ever recorded in the region and one of the few rare thrust earthquakes in a region otherwise characterized by normal faulting and strike slip earthquakes. The joint analysis of seismicity and focal mechanisms suggest the activation of E-W and NE-SW faults or of an internal curved pseudofault, which is formed in response to the microplate rotation, with alternation of thrust and strike-slip earthquakes. Seismicity migrated Northward in its final phase, towards the microplate edge, where a second doublet with uneven focal mechanisms occurred. The sequence rupture kinematics is well explained by Coulomb stress changes imparted by the first subevent. Our analysis show that compressional stresses, which have been mapped at the northern boundary of the microplate, but never accompanied by large thrust earthquakes, can be accommodated by the rare occurrence of large, impulsive, shallow thrust earthquakes, with a considerable tsunamigenic potential.&lt;/p&gt;


2021 ◽  
Author(s):  
Figen Eskikoy ◽  
Semih Ergintav ◽  
Uğur Dogan ◽  
Seda Özarpacı ◽  
Alpay Özdemir ◽  
...  

&lt;p&gt;On 2020 October 30, an M&lt;sub&gt;w&lt;/sub&gt;6.9 earthquake struck offshore&amp;#160;Samos Island. Severe structural damages were observed in Greek Islands and city of Izmir (Turkey). 114 people lost their lives and more than a thousand people were injured in Turkey. The earthquake triggered local tsunami. Significant seismic activity occurred in this region following the earthquake and ~1800 aftershocks (M&gt;1) were recorded by KOERI within the first three days. In this study, we analyze the slip distribution and aftershocks of the 2020 earthquake.&lt;/p&gt;&lt;p&gt;For the aftershock relocations, the continuous waveforms were collected from NOA, Disaster and Emergency Management Authority of Turkey (AFAD) and KOERI networks. The database &amp;#160;&amp;#160;was created based on merged catalogs from AFAD and KOERI. For estimating optimized aftershock location distribution, the P and S phases of the aftershocks are picked manually and relocated with double difference algorithm. In addition, source mechanisms of aftershocks M&gt;4 are obtained from regional body and surface waveforms.&lt;/p&gt;&lt;p&gt;The surface deformation of the earthquake was obtained from both descending and ascending orbits of the Sentinel-1 A/B and ALOS2 satellites. Since the rupture zone is beneath the Gulf of Kusadas&amp;#305;, earthquake related deformation in the interferograms can only be observed on the northern part of the Samos Island. We processed all possible pairs chose the image pairs with the lowest noise level.&lt;/p&gt;&lt;p&gt;In this study, we used 25 continuous GPS stations which are compiled from TUSAGA-Aktif in Turkey and NOANET in Greece. In addition to continuous GPS data, on 2020 November 1, GPS survey was initiated and the earthquake deformation was measured on 10 GNSS campaign sites (TUTGA), along onshore of Turkey.&lt;/p&gt;&lt;p&gt;The aim of this study is to estimate the spatial and temporal rupture evolution of the earthquake from geodetic data jointly with near field displacement waveforms. To do so, we use the Bayesian Earthquake Analysis Tool (BEAT).&lt;/p&gt;&lt;p&gt;As a first step of the study, rectangular source parameters were estimated by using GPS data. In order to estimate the slip distribution, we used both ascending and descending tracks of Sentinel-1 data, ALOS2 and GPS displacements. In our preliminary geodetic data based finite fault model, we used the results of focal mechanism and GPS data inversion solutions for the initial fault plane parameters. The slip distribution results indicate that earthquake rupture is ~35 km long and the maximum slip is ~2 m normal slip along a north dipping fault plane. This EW trending, ~45&amp;#176; north dipping normal faulting system consistent with this tectonic regime in the region. This seismically active area is part of a N-S extensional regime and controlled primarily by normal fault systems.&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Acknowledgements&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;This work is supported by the Turkish Directorate of Strategy and Budget under the TAM Project number&amp;#160;2007K12-873.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document