Search

S-net is a large-scale ocean bottom (OB) network in the Japan Trench area, consisting of inline-type 150 observatories equipped with seismometers and pressure gauges. Among them, 41 observatories have been buried about one meter beneath the seafloors in the shallow water regions (water depth <1500 m). We analyzed the strong-motion data recorded at the S-net sites from earthquakes with magnitudes 3.5 < Mw ≤ 7 and focal depths < 70 km to understand the site amplification effect on the recorded motions. We used the spectral inversion technique and obtained some fundamental properties of the earthquake source spectra, path attenuation, and site factors from the horizontal-component S-wave portions of the recordings. We obtained that the source spectra followed the ω - 2 source model generally well, and the estimated magnitudes were mostly within ± 0.3 magnitude units of the catalog magnitudes. The stress drops increased systematically with the focal depths, and the values for the shallow earthquakes in the Pacific Plate were higher than those for the interplate earthquakes with comparable focal depths. The path-averaged quality factors were generally frequency-dependent and were somewhat larger than those in the past studies. The peak site frequencies ranged between about 0.2 and 10 Hz, while the peak amplification factors ranged between 10 and 50. Even though the peak frequencies and amplification factors differed from site to site, the peak frequencies were mostly higher than 2 Hz at the outer trench stations, while they were lower than 2 Hz at many inner trench stations. The amplification factors at a few OB sites in the shallow water regions were comparable with the theoretical ones computed from the 1-D subsea model. The amplification factors at intermediate frequencies (~ 0.3 to 2 Hz) generally increased with P-wave travel time in the sediments estimated from the multi-channel seismic survey. At about 20% of the sites (mainly at the unburied stations), spurious site spectra were recognized at frequencies higher than about 4 Hz. If the site spectra between 4 and 10 Hz are required, using only the X-component records is recommended. [ABSTRACT FROM AUTHOR]

Density perturbations accompanying seismic waves are expected to generate prompt gravity perturbations preceding the arrival of P-waves. Vallée et al. (Science 358:1164–1168, 2017, https://doi.org/10.1126/science.aao0746) reported the detection of such pre-P-wave signals in broadband seismograms during the 2011 Tohoku-oki earthquake. Kimura et al. (Earth Planets Space 71:27, 2019, https://doi.org/10.1186/s40623-019-1006-x) considered that their detection involved some uncertain points, including a concern regarding their signal processing procedure. Specifically, to remove the instrumental response, Vallée et al. (2017) applied acausal deconvolution to the seismograms truncated at the P-wave arrivals. Generally, acausal deconvolution produces artifacts at the edge of the time window. However, they did not present quantitative assessment whether the detected signals were artifacts due to the signal processing. To avoid this concern, Kimura et al. (2019) employed another procedure that eliminated acausal processes, resulting in the detection of a pre-P-wave signal with a statistical significance of 7σ in stacked broadband seismograms. Subsequently, Vallée et al. (Earth Planets Space 71:51, 2019, https://doi.org/10.1186/s40623-019-1030-x) commented that the procedure employed by Kimura et al. (2019) for the signal detection was inappropriate because it dismissed the low-frequency components of data. Although we admit the loss of low-frequency components in the data in Kimura et al. (2019), Vallée et al. (2019) have not yet provided a full account of the validity of their own procedure. Here, we assessed the validity of the procedure employed by Vallée et al. (2017) by quantitatively evaluating the magnitude of the acausal artifacts. First, we investigated how the input acceleration waveform, having an ideal signal-like shape, was distorted by their procedure. Their acausal deconvolution indeed generated a large-amplitude terminal artifact; however, it was removed by the causal band-pass filtering performed after the deconvolution and consequently became negligible. Next, we constrained the maximum amplitude of the artifact due to the noise in a seismogram and showed that it was sufficiently small compared to the reported signal amplitudes. These results suggest that the signal waveforms seen after their procedure were not artifacts but were representing the input acceleration with sufficient accuracy. Namely, their procedure well functions as a detection method for pre-P-wave signals. In the context of this validation, we replied to the comments of Vallée et al. (2019). [ABSTRACT FROM AUTHOR]

Dynamic earthquake rupture causes mass redistribution around the fault, and the emitted propagating seismic waves are accompanied by bulk density perturbations. Both processes cause transient gravity changes prior to the arrival of P-waves. Such pre-P gravity signals have been detected in previous studies of several large earthquakes. However, the detections were limited to the vertical component of the signal owing to the high noise level in the horizontal records. In this study, we analyzed dense tiltmeter array data in Japan to search for the horizontal components of the signal from the 2011 Mw 9.1 Tohoku-Oki earthquake. Based on the synthetic waveforms computed for a realistic Earth model, we stacked the horizontal records and identified a signal that evidently exceeded the noise level. We further performed a waveform inversion analysis to estimate the source parameters. The horizontal tiltmeter data, combined with the vertical component of the broadband seismometer array data, yielded a constraint on the dip angle and magnitude of the earthquake in the ranges of 11.5°–15.3° and 8.75°–8.92°, respectively. Our results indicate that the analysis of the three components of the pre-P gravity signal avoids the intrinsic trade-off problem between the dip angle and seismic moment in determining the source mechanism of shallow earthquakes. Pre-P gravity signals open a new observation window for earthquake source studies. [ABSTRACT FROM AUTHOR]

Density perturbations accompanying seismic waves are expected to generate prompt gravity perturbations preceding the arrival of P-waves. Vallée et al. (Science 358:1164–1168, 2017, https://doi.org/10.1126/science.aao0746) reported the detection of such pre-P-wave signals in broadband seismograms during the 2011 Tohoku-oki earthquake. Kimura et al. (Earth Planets Space 71:27, 2019, https://doi.org/10.1186/s40623-019-1006-x) considered that their detection involved some uncertain points, including a concern regarding their signal processing procedure. Specifically, to remove the instrumental response, Vallée et al. (2017) applied acausal deconvolution to the seismograms truncated at the P-wave arrivals. Generally, acausal deconvolution produces artifacts at the edge of the time window. However, they did not present quantitative assessment whether the detected signals were artifacts due to the signal processing. To avoid this concern, Kimura et al. (2019) employed another procedure that eliminated acausal processes, resulting in the detection of a pre-P-wave signal with a statistical significance of 7σ in stacked broadband seismograms. Subsequently, Vallée et al. (Earth Planets Space 71:51, 2019, https://doi.org/10.1186/s40623-019-1030-x) commented that the procedure employed by Kimura et al. (2019) for the signal detection was inappropriate because it dismissed the low-frequency components of data. Although we admit the loss of low-frequency components in the data in Kimura et al. (2019), Vallée et al. (2019) have not yet provided a full account of the validity of their own procedure. Here, we assessed the validity of the procedure employed by Vallée et al. (2017) by quantitatively evaluating the magnitude of the acausal artifacts. First, we investigated how the input acceleration waveform, having an ideal signal-like shape, was distorted by their procedure. Their acausal deconvolution indeed generated a large-amplitude terminal artifact; however, it was removed by the causal band-pass filtering performed after the deconvolution and consequently became negligible. Next, we constrained the maximum amplitude of the artifact due to the noise in a seismogram and showed that it was sufficiently small compared to the reported signal amplitudes. These results suggest that the signal waveforms seen after their procedure were not artifacts but were representing the input acceleration with sufficient accuracy. Namely, their procedure well functions as a detection method for pre-P-wave signals. In the context of this validation, we replied to the comments of Vallée et al. (2019). [ABSTRACT FROM AUTHOR]

Inter-station cross-correlation functions estimated using continuous ambient noise or microtremor records were used to extract the seismic wave propagation characteristics of the Osaka sedimentary basin, Japan. Temporary continuous observations were conducted at 15 sites in the Osaka basin between 2011 and 2013. The data were analyzed using seismic interferometry. The target period range was 2-8 s. Cross-correlations between all of the possible station pairs were calculated and stacked to produce a year-long data set, and Rayleigh wave signals in the vertical and radial components and Love wave signals in the transverse component were identified from the results. Simulation of inter-station Green's functions using the finite difference method was conducted to check the performance of the current three-dimensional velocity structure model. The measured time lag between the observed and theoretical Green's functions was less than 2 s for most station pairs, which is less than the wave period of interest in the target frequency range. Group velocity tomography was applied to group delay times estimated by means of multiple filter analysis. The estimated group velocities for longer periods of 5-8 s exhibited spatial variation within the basin, which is consistent with the bedrock depth distribution; however, the group velocities for shorter periods of 2-3 s were almost constant over the studied area. The waveform and group velocity information obtained by seismic interferometry analysis can be useful for future reconstruction of a three-dimensional velocity structure model in the Osaka basin. [ABSTRACT FROM AUTHOR]

We analyzed S-wave horizontal-to-vertical (S-H/V) spectral ratios at six ocean bottom seismograph (OBS) sites of K-NET located in the Sagami Bay area of Japan for nonlinear site responses. The degree of nonlinearity was computed by comparing the S-H/V spectral ratios for strong motions (PGA ≥ 20 cm/s) with those for weak motions (PGA < 20 cm/s). Our analyses, which showed that the weak-motion S-H/V spectral ratios differ from site to site, indicate that the underlying site geology is not uniform at the OBS sites. It was found that the threshold PGA causing a nonlinear site response is generally different from site to site. Recordings having horizontal PGAs greater than about 50-150 cm/s display clear signatures of nonlinear site effects, i.e., the shift of predominant frequencies to lower ones and/or the decrease in high-frequency spectral ratios. We also found that the degree of nonlinearity is generally larger at the OBS sites due to the smaller threshold motions that cause a nonlinear site response compared with the available data at land sites. The above findings suggest the possibility of a widespread nonlinear site response at the OBS sites for offshore earthquakes with a large magnitude. However, frequencies lower than about 2 Hz are not affected by the nonlinear site response in the analyzed data ranges (PGA < 467 cm/s). These results indicate the need for careful utilization of recorded strong motions at OBS sites for applications such as real-time ground motion predictions as front detections. [ABSTRACT FROM AUTHOR]

The nationwide strong-motion seismograph network of K-NET and KiK-net in Japan successfully recorded the strong ground motions of the 2016 Kumamoto earthquake sequence, which show the several notable characteristics. For the first large earthquake with a JMA magnitude of 6.5 (21:26, April 14, 2016, JST), the large strong motions are concentrated near the epicenter and the strong-motion attenuations are well predicted by the empirical relation for crustal earthquakes with a moment magnitude of 6.1. For the largest earthquake of the sequence with a JMA magnitude of 7.3 (01:25, April 16, 2016, JST), the large peak ground accelerations and velocities extend from the epicentral area to the northeast direction. The attenuation feature of peak ground accelerations generally follows the empirical relation, whereas that for velocities deviates from the empirical relation for stations with the epicentral distance of greater than 200 km, which can be attributed to the large Love wave having a dominant period around 10 s. The large accelerations were observed at stations even in Oita region, more than 70 km northeast from the epicenter. They are attributed to the local induced earthquake in Oita region, whose moment magnitude is estimated to be 5.5 by matching the amplitudes of the corresponding phases with the empirical attenuation relation. The real-time strong-motion observation has a potential for contributing to the mitigation of the ongoing earthquake disasters. We test a methodology to forecast the regions to be exposed to the large shaking in real time, which has been developed based on the fact that the neighboring stations are already shaken, for the largest event of the Kumamoto earthquakes, and demonstrate that it is simple but effective to quickly make warning. We also shows that the interpolation of the strong motions in real time is feasible, which will be utilized for the real-time forecast of ground motions based on the observed shakings. [ABSTRACT FROM AUTHOR]

Strong ground motions from the 2015 Mw 7.8 Gorkha, Nepal, earthquake and its eight aftershocks recorded by a strong-motion seismograph at Kantipath (KATNP), Kathmandu, were analyzed to assess the ground-motion characteristics and site effects at this location. Remarkably large elastic pseudo-velocity responses exceeding 300 cm/s at 5 % critical damping were calculated for the horizontal components of the mainshock recordings at peak periods of 4-5 s. Conversely, the short-period ground motions of the mainshock were relatively weak despite the proximity of the site to the source fault. The horizontal components of all large-magnitude (Mw ≥ 6.3) aftershock recordings showed peak pseudo-velocity responses at periods of 3-4 s. Ground-motion prediction equations (GMPEs) describing the Nepal Himalaya region have not yet been developed. A comparison of the observational data with GMPEs for Japan showed that with the exception of the peak ground acceleration (PGA) of the mainshock, the observed PGAs and peak ground velocities at the KATNP site are generally well described by the GMPEs for crustal and plate interface events. A comparison of the horizontal-to-vertical ( H/ V) spectral ratios for the S-waves of the mainshock and aftershock recordings suggested that the KATNP site experienced a considerable nonlinear site response, which resulted in the reduced amplitudes of short-period ground motions. The GMPEs were found to underestimate the response values at the peak periods (approximately 4-5 s) of the large-magnitude events. The deep subsurface velocity model of the Kathmandu basin has not been well investigated. Therefore, a one-dimensional velocity model was constructed for the deep sediments beneath the recording station based on an analysis of the H/ V spectral ratios for S-wave coda from aftershock recordings, and it was revealed that the basin sediments strongly amplified the long-period components of the ground motions of the mainshock and large-magnitude aftershocks. [ABSTRACT FROM AUTHOR]

National Research Institute for Earth Science and Disaster Resilience (NIED) integrated the land observation networks established since the 1995 Kobe earthquake with the seafloor observation networks established since the 2011 Tohoku earthquake and tsunami as MOWLAS (Monitoring of Waves on Land and Seafloor) in November 2017. The purpose of MOWLAS is to provide comprehensive, accurate, and rapid observation and monitoring of earthquake, tsunami, and volcano events throughout Japan and its offshore areas. MOWLAS data are widely utilized for long-term earthquake forecasting, the monitoring of current seismic activity, seismic and tsunami hazard assessments, earthquake early warning, tsunami warning, and earthquake engineering, as well as earthquake science. Ocean bottom observations provide an extension of observations to areas where no people are living and have the advantage of increasing lead time of earthquake early warning and tsunami warning. The application of recent technology advancements to real-time observations as well as the processing of MOWLAS data has contributed to the direct disaster mitigation of ongoing earthquakes. These observations are fundamental for both science and disaster resilience, and thus it is necessary to continue ceaseless operation and maintenance. [ABSTRACT FROM AUTHOR]

The observed peak ground accelerations and peak ground velocities (PGVs) of the 2018 Mw 6.6 Hokkaido eastern Iburi earthquake generally followed the median values from ground motion prediction equations with reasonable errors at fault distances ≥ 50 km. However, at smaller distances, the equations significantly underestimated the peak ground motions, and it was eminent for PGVs. A comparison of surface-to-borehole spectral ratios of S-waves during the mainshock and other events revealed that the sites at smaller distances experienced various degrees of nonlinear site response. The two most widely known characteristics of nonlinear site response are the weakening of higher-frequency components and shifting of predominant frequencies to lower ones in comparison with the linear site response. At one of the sites that recorded the largest intensity of 7 in JMA scale of 0–7, the latter nonlinear effect was so dominant that the ground motions around the new predominant frequency got intensified by one order of magnitude in comparison with that during the weak-motions. Two sites, which were closely located, recorded vertical peak ground accelerations exceeding 1 g for the up-going motions. The recordings showed asymmetric waveforms and amplitudes characteristics of the nonlinear site response in extreme vertical ground motions recorded during a few earthquakes in the past. Few sites having lower vertical peak ground accelerations were also suspected of being experienced nonlinear site response on vertical motions. These findings suggest taking a cautious approach to enumerate the reduction in amplification at higher frequencies using the single-station horizontal-to-vertical (H-to-V) spectral ratio technique. However, we found that the H-to-V technique was still useful to detect nonlinearity. Finally, an ad hoc equation was derived to correct the nonlinear site amplification in predicting horizontal PGVs with respect to one of the most widely used attenuation models in Japan. The results indicated that the effect was much stronger for a larger input motion than that for a proportional change in the Vs30 values. [ABSTRACT FROM AUTHOR]

Earthquake ruptures cause mass redistribution, which is expected to induce transient gravity perturbations simultaneously at all distances from the source before the arrival of P-waves. A recent research paper reported the detection of such prompt gravity signals from the 2011 Tohoku-Oki earthquake by comparing observed acceleration waveforms and model simulations. The 11 observed waveforms presented in that paper recorded in East Asia shared a similar trend above the background seismic noise and were in good agreement with the simulations. However, the signal detection was less quantitative because the significance of the observed signals was not discussed and the waveforms at other stations in the region were not shown. In this study, similar trends were not observed in most of the data recorded near the stations used in the aforementioned study, suggesting that the reported signals were only local noises. We thus took a different approach to identify the prompt signals. We optimized the multi-channel data recorded by superconducting gravimeters, broadband seismometers, and tiltmeters. Though no signal was identified in the single-trace records, the stacked trace of the broadband seismometer array data in Japan showed a clear signal above the reduced noise level. The signal amplitude was 0.25 nm/s2 for an average distance of 987 km from the event hypocenter. This detection was confirmed with a statistical significance of 7σ, where σ is the standard deviation of the amplitude of the background noise. This result provided the first constraint on the amplitude of the observed prompt signals and may serve as a reference in the detection of prompt signals in future earthquakes. [ABSTRACT FROM AUTHOR]

We analyzed strong-motion and broadband recordings of two moderate-magnitude earthquakes that occurred in the Nankai Trough. The first event was the 2004 Mw 6.5 southeast off-Kii peninsula earthquake, an aftershock event inside the Philippine Sea Plate near the Nankai Trough axis. The second event was the 2016 Mw 5.8 southeast off-Mie Prefecture earthquake, an independent event in the rupture area of the 1944 Mw ~ 8 Tonankai earthquake. The centroid depths were 11 and 14 km for the 2004 and 2016 events, respectively. Despite a large difference in the moment magnitude between the two events, the JMA magnitude (Mj) was 6.5 for both the events. We found that the short-period ground motions (e.g., response spectra at periods < 1 s) as well as the much longer-period ground motions (> 20 s) for the 2016 event scaled generally well with the moment magnitude of the event. In contrast, the ground motions from the 2016 event were comparable to those for the larger-moment-magnitude 2004 event at equal distances at periods of about 2-20 s in wide areas and the observed acceleration response spectra at those periods were noticeably underestimated for the 2016 event by the ground motion prediction equation (GMPE) that employs Mw. An examination of the existing subsurface velocity model suggested that the difference in the relative location of the two events with respect to the thick accretionary prism of low seismic velocity most probably caused the comparable amplitude of the seismic waves at those periods. As a result, we posit that the values of Mj are equal for the two events because Mj is estimated using the displacement amplitude of ground motions at periods smaller than about 6 s. On the other hand, GMPE employing Mj generally described the observed data well. The results suggested that the plate boundary earthquakes in the Nankai Trough may excite strong long-period ground motions of engineering importance, and these ground motions appear to be better explained by Mj than by Mw in GMPEs for moderate-magnitude earthquakes in the Nankai Trough subduction zone. [ABSTRACT FROM AUTHOR]