Geodetic investigation of triggered slip below Fiordland

<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>

Geodetic investigation of triggered slip below Fiordland

<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>

Tsunami Hazard Assessment Along Chinese Coast Using Scaling Relations Developed for Tsunami Prediction

Potential tsunamis in the western Pacific Ocean pose great threats to the Chinese coastal areas. Among all possible tsunami source regions, the Manila subduction zone draws the most attention and there have been many research works on the tsunami hazards in the South China Sea. In this study, we evaluate the tsunami hazard along the Chinese coast by investigating more potential sources, including the subduction zones of Manila, Ryukyu, Nankai, Izu–Bonin and Mariana. Two tsunami scenarios are considered for each subduction zone, a worst scenario of earthquake magnitude 9.0 and a scenario of largest earthquake magnitude known in history in this zone. Earthquake source parameters are calculated using scaling relations that have been shown to be suitable for tsunami generation. Our results show that for the Chinese coast, tsunami hazards from the Manila and Ryukyu subduction zones are severe in the worst scenarios, and tsunami hazards from the Nankai, Izu–Bonin and Mariana subduction zones are mild. Using the largest earthquake magnitude in history, tsunami hazards from all the investigated subduction zones are almost negligible. Through a sensitivity test on earthquake magnitude, we find that earthquakes of magnitude of 8.5 or larger in the Manila and Ryukyu subduction zones cause severe tsunami hazard along the Chinese coast with wave amplitude over 2 m.

Parallel Dislocation Model Implementation for Earthquake Source Parameter Estimation on Multi-Threaded GPU

Graphics processing units (GPUs) have been in the spotlight in various fields because they can process a massive amount of computation at a relatively low price. This research proposes a performance acceleration framework applied to Monte Carlo method-based earthquake source parameter estimation using multi-threaded compute unified device architecture (CUDA) GPU. The Monte Carlo method takes an exhaustive computational burden because iterative nonlinear optimization is performed more than 1000 times. To alleviate this problem, we parallelize the rectangular dislocation model, i.e., the Okada model, since the model consists of independent point-wise computations and takes up most of the time in the nonlinear optimization. Adjusting the degree of common subexpression elimination, thread block size, and constant caching, we obtained the best CUDA optimization configuration that achieves 134.94×, 14.00×, and 2.99× speedups over sequential CPU, 16-threads CPU, and baseline CUDA GPU implementation from the 1000×1000 mesh size, respectively. Then, we evaluated the performance and correctness of four different line search algorithms for the limited memory Broyden–Fletcher–Goldfarb–Shanno with boundaries (L-BFGS-B) optimization in the real earthquake dataset. The results demonstrated Armijo line search to be the most efficient one among the algorithms. The visualization results with the best-fit parameters finally derived by the proposed framework confirm that our framework also approximates the earthquake source parameters with an excellent agreement with the geodetic data, i.e., at most 0.5 cm root-mean-square-error (RMSE) of residual displacement.

Modeling of earthquake source parameters on December 12, 2018 (08:49:56,16; 36,4478° N; 140,5788° E; H = 62,0 km; Mw = 4,3, Japan)

Modeling of earthquake source parameters, such as the orientation of the fault plane and the direction of the fault slip, is important for understanding the physics of earthquake source processes, determining the stress-strain state of the geological medium and seismic hazard estimation. For modeling source parameters of the earthquake on December 12, 2018 at 08:49:56,16 (UTC) in Japan (36,4478° N, 140,5788° E, Northern Ibaraki Pref region) at a depth of 62 km with a magnitude of Mw = 4.3, the waveforms inversion was used to determine seismic moment tensor and representation it through a focal mechanism. The earthquake source is considered as a point source of seismic waves which propagate in a medium represented by a set of horizontally homogeneous elastic layers. An algorithm for determining seismic tensor components based on the forward problem solved by the matrix method, and using the generalized inverse solution, selecting only direct waves is applied. The input data for determining seismic moment components are data of only direct P waves selected from the observed records at six seismic stations of the Japanese local network NIED F-net: TSK, YMZ, ASI, ONS, SBT, KSK. The seismic moment tensor components were determined through waveform inversion using the matrix method. The obtained results, presented through a focal mechanism, are compared to the results obtained by the National Research Institute of Earth Sciences and Resistance to Natural Disasters (NIED CMT solutions). As a result of focal mechanisms comparison, it is concluded that the proposed algorithm for determining seismic moment tensor components can be used if it is impossible to use another method, or requires some refinement for another method. This approach is especially relevant for regions with low seismicity and insufficient number of stations. In addition, this method reduces the effects of an inaccurate medium model, because direct waves are much less distorted than reflected and converted, and that increases the accuracy and reliability of the method.

Strong ground motion simulation of the 27 May 2006 Yogyakarta earthquake using Empirical Green’s Function

On May 27, 2006, 05:54 am local time, a moderate crustal earthquake of magnitude Mw 6,3 struck the Yogyakarta province, especially in the Bantul regency in the south part of the province. The earthquake damaged or destroyed more than 400,000 houses and buildings and caused more than 5,700 people killed. Several earthquake stations recorded the ground vibration caused by the mainshock very well, except at the stations closest to the earthquake source, namely YOGI in Gamping, West of Yogyakarta, which experienced saturation due to significant vibration. Therefore, information about the maximum ground acceleration near the source is yet not known. We model the ground vibrations near the earthquake source using a stochastic Green’s Function approach to obtain information about the ground motions’ maximum amplitude. The earthquake source parameters we referred to is the moment tensor solution from the Harvard Moment Tensor. The calculations show that the amplitude is consistent with observations recorded at the BJI Banjarnegara (0.04g) and YOGI Yogyakarta (0.3g).

Identification of Ransiki fault segment in South Manokwari Regency, West Papua Province, Indonesia based on analysis of a high-resolution of global gravity field: Implications on the Earthquake Source Parameters

The province of West Papua in Indonesia is an area crossed by three major faults, including Sorong, Koor, and Ransiki, leading to the collision of Australia, the Pacific, and Eurasia. In the past, there have been strong and damaging earthquakes on these faults, manly Ransiki fault in the South Manokwari regency. Identification of the Ransiki fault segment was conducted by geological subsurface modeling using the earth gravity field of the Global Gravity Map (GGM) based on satellite measurements implicates for earthquake source parameters. The GGM is seen as a solution for the unavailability of direct measurements in the region. The gravity field analysis begins with data reduction using SRTM2gravity as modern terrain correction to obtain a complete Bouguer anomaly. Furthermore, the gravity gradient approach through vertical and horizontal gradients, analytical signal, and the tilt angle are applied to emphasize a contact or fault structures that are not visible using a 2D fast Fourier transform. Overall, the gravity gradient analysis obtained results that were compatible with the alignment of the Ransiki fault segment which direction of the northwest to south. The gravity inversion produces a geological subsurface model that clearly shows the Ransiki fault segment, associated with a low rock density distribution, thought to the Befoor formation and quaternary sediments, located between high-density rocks correlated to Arfak volcanic rocks as a basement.

Strain to ground motion conversion of distributed acoustic sensing data for earthquake magnitude and stress drop determination

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' 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 earthquakes, including those capable of generating tsunami.