Search

The oldest oceanic basin (160–180 Ma) in the western Pacific is the birthplace of the Pacific Plate and is thus essential for understanding the formation and evolution of the oceanic plate. However, the upper mantle structure beneath the region has not been thoroughly investigated because of the remoteness and difficulties of long‐term in situ seismic measurements at the ocean bottom. From 2018 to 2019, the Oldest‐1 experiment on the oldest seafloor was conducted as part of the international Pacific Array initiative. We present the first three‐dimensional P‐wave velocity structure down to a depth of 350 km based on the relative travel time residuals of teleseismic earthquakes recorded by 11 broadband ocean‐bottom seismometers operated during the Oldest‐1 experiment. Our result shows a fast P‐wave velocity anomaly (VP perturbation of 2%–4% faster than average) at a depth of 95–185 km beneath the northeast of the study area. This structure is interpreted as evidence of dry, viscous, and rigid materials at depths below the lithosphere. Two slow anomalies (VP perturbation of 2%–4% slower than average) are seen beneath the southwestern and eastern (the oldest seafloor >170 Ma) parts of the array site. The low‐velocity zones are found at depths of 95–305 km. The observed velocity structures can be indicative of plume activities that affected the upper mantle as the Pacific Plate migrated over hotspots from the southeast. Alternatively, the observed velocity features may provide seismic evidence for small‐scale sublithospheric convection. Key Points: One‐year ocean‐bottom geophysical investigation on the oldest Pacific provides seismic mantle structure of the regionDetailed 3‐D mantle structure implies complex evolution process of Pacific PlateOur model implies thermochemical modification of the upper mantle by plume interaction or small‐scale convection [ABSTRACT FROM AUTHOR]

Inversion for seismological structures of the oceanic lithosphere–asthenosphere system is important to understand the mechanisms of plate tectonics. Previous models of the oceanic upper mantle have been primarily obtained via global tomography using surface waveforms. However, besides scarcity of waveform data in the oceanic regions, difficulties in fitting phases for shorter-period components in the previous global tomography have yielded resultant models that possess poor resolutions above |$\sim \, 50$| km depth. Recent developments of broad-band ocean-bottom seismometer (BBOBS) arrays provide larger amount of seismic data with epicentral distances of <20°. In this study, we develop an appropriate method to fully utilize the information contained in the shorter-period components of BBOBS arrays. We first fit the envelopes without phase information to analyse the shorter-period components (8–60 s) which are generally unavailable in the conventional phase fitting. We then use the resultant model as our initial model for waveform inversion of the longer periods (12.5–200 s) to fit the phase, which allows us to infer a continuous structure model from the crust to the asthenosphere. We demonstrate the validity of this combined envelope-fitting and waveform inversion method by analysing the waveform data from a BBOBS array that was deployed in the Northwestern Pacific and has recorded events in the vicinity of the Japan Trench to obtain the average velocity structure between the event and station arrays. We independently resolve the crustal compressional and shear wave velocities, and thickness by analysing the envelopes, which minimizes biases in the subsequent waveform inversion. We also find that the waveform inversion improves the resolution in the asthenosphere. Our results suggest that further extension of this method should improve our knowledge of the oceanic lithosphere–asthenosphere system. [ABSTRACT FROM AUTHOR]

A study of the seismograms from a seismic array in the Japanese islands offers evidence for 920-km deep seismic discontinuity in the Earth's mantle below the Tonga subduction zone. The discontinuity may be attributed to a change in the mantle's chemical content or the signature of the garnet layer of the subducted zone. The discontinuity also exists below subduction areas in the Flores and Japan seas.

Strong long‐period seismic signals at periods around 25 and 18 s appear in the ambient noise cross‐correlation functions recorded by an array of ocean bottom seismometers (OBSs) in the western‐central Pacific. The signal amplitude varies from time to time, and the apparent travel times of the signals are typically smaller than those expected for the Rayleigh waves propagating along the great circle connecting the station pairs. From the cross‐correlation functions, the signal sources are located in the Vanuatu Arc. Local data analysis suggests the signals originate from two different sources possibly located at depths of ~0–1 km below the sea level beneath the active cones of Ambrym volcano. Plain Language Summary: Spatially localized sources that persistently generate long‐period seismic waves at periods longer than 10 s are reported worldwide. For some sources, the excitation mechanisms have been well‐investigated. A volcano in Japan generates such waves that travel more than thousands of kilometers. For others, however, the excitation origins are still poorly understood. These include the one observed in the Vanuatu Arc. A temporal observation network of seismometers installed on the seafloor in the western‐central Pacific observed persistent waves peaking at around two different periods. From the data analysis, we find the signals originate from an active volcano, Ambrym volcano, in the Vanuatu Arc. Besides, we find that those seismic waves are most likely originated from two different sources possibly located at shallow depths beneath the volcano. Key Points: Persistent long‐period (25 and 18 s) seismic signals appear in ambient noise cross‐correlations of OBS data in the western‐central PacificThe signals are originated from Ambrym volcano in the Vanuatu ArcSource locations of the two signals may be different, but both are possibly located at shallow depths right beneath the active cones [ABSTRACT FROM AUTHOR]

We found that SH diff phases generated by earthquakes in the Fiji–Tonga, recorded in India, are accompanied by secondary pulses. We interpreted them as a consequence of multipathing of S waves caused by the Pacific large low-shear-velocity province (LLSVP). We analysed the differential traveltimes between SH diff and the secondary pulse, together with the absolute SH diff arrival times, to constrain the thickness and velocity perturbations in the western end of the Pacific LLSVP. Our preferred model shows a lateral variation in the thickness of the LLSVP; the southern part reveals a thicker (300 km) low-velocity region compared to the northern part (200 km). However, the velocity perturbations of the LLSVP appear to be comparable (⁠|$-1.5\hbox{ per cent}$|⁠). The results are consistent with a scenario that the LLSVP is a chemically distinct pile with significant surface topography. [ABSTRACT FROM AUTHOR]

We conducted broadband dispersion survey by deploying two arrays of broadband ocean bottom seismometers in the northwestern Pacific Ocean at seafloor ages of 130 and 140 Ma. By combining ambient noise and teleseismic surface wave analyses, dispersion curves of Rayleigh waves were obtained at a period range of 5–100 s and then used to invert for one‐dimensional isotropic and azimuthally anisotropic βV (VSV) profiles beneath each array. The obtained profiles show ~2% difference in isotropic βV in the low‐velocity zone (LVZ) at a depth range of 80–150 km in spite of the small difference in seafloor ages and the horizontal distance of ~1,000 km. Forward dispersion‐curve calculation for thermal models indicates that simple cooling models cannot explain the observed difference and an additional mechanism, such as sublithospheric small‐scale convection, is required. In addition, the fastest azimuths of azimuthal anisotropy in the LVZ significantly deviate from the current plate motion direction. We infer that these observations are consistent with the presence of small‐scale convection beneath the study area. As for azimuthal anisotropy in the Lid, the peak‐to‐peak intensity is 3–4% at the depth from Moho to ~40 km. The fastest direction is almost perpendicular to magnetic lineation in area A at 130 Ma and oblique to magnetic lineations in area B at 140 Ma, suggesting complex mantle flow beneath the infant Pacific Plate surrounded by three ridge axes. The intensity of azimuthal anisotropy in the LVZ is ~2%, indicating that radial anisotropy is stronger than azimuthal anisotropy therein. Key Points: Broadband ocean bottom seismic observation is conducted in the NW Pacific Ocean at two sites with seafloor ages of 130 and 140 MaSignificant seismic velocity variation in the LVZ of two areas is revealed that may indicate the presence of small‐scale convectionAzimuthal anisotropy depth profiles indicate that the fast directions do not necessarily coincide with the ancient spreading directions [ABSTRACT FROM AUTHOR]

We investigate slab and mantle structure of the NE Japan subduction zone from P wave azimuthal and radial anisotropy using travel time tomography. Trench normal E-W-trending azimuthal anisotropy (AA) and radial anisotropy (RA) with VPV > VPH are found in the mantle wedge, which supports the existence of small-scale convection in the mantle wedge with flow-induced LPO of mantle minerals. In the subducting Pacific slab, trench parallel N-S-trending AA and RA with VPH > VPV are obtained. Considering the effect of dip of the subducting slab on apparent anisotropy, we suggest that both characteristics can be explained by the presence of laminar structure, in addition to AA frozen-in in the subducting plate prior to subduction. [ABSTRACT FROM AUTHOR]

Abstract: We report anomalous core phase PKPbc‐PKPdf differential travel times relative to 1‐D spherically symmetric model with a uniformly anisotropic inner core recorded by stations in Alaska for South Sandwich Islands (SSI) earthquakes. The data sample the inner core for the polar paths, as well as the lowermost mantle beneath Alaska. Our major observations are the following: (1) fractional travel time residuals of PKPbc‐PKPdf increase rapidly within 2 °, (2) a clear shift of the residual pattern could be seen for earthquakes with different locations, and (3) the residuals show systematic lateral variation: at the northern part, no rapid increase of residual can be seen. A structural boundary with a P wave velocity contrast of about 3% at the lowermost mantle beneath East Alaska is invoked to explain the observation, and the required strength of anisotropy in the quasi‐western hemisphere of the inner core might be reduced if those anomalous data are excluded from analysis. [ABSTRACT FROM AUTHOR]

Properties of a new parameter, ηκ, that is recently introduced by Kawakatsu et al. for transverse isotropy are examined. It is illustrated that the parameter nicely characterizes the incidence angle dependence of bodywave phase velocities for vertical transverse isotropy models that share the same P- and S-wave anisotropy. When existing models of upper-mantle radial anisotropy are compared in terms of this new parameter, PREM shows a distinct property. Within the anisotropic layer of PREM (a depth range of 24.4-220 km), ηκ < 1 in the upper half and ηκ > 1 in the lower half. If ηκ > 1, anisotropy cannot be attributed to a layering of homogeneous isotropic layers, and thus requires the presence of intrinsic anisotropy. [ABSTRACT FROM AUTHOR]

We have developed a real-time moment tensor monitoring system (RMT) which takes advantage of a grid-based moment tensor inversion technique and real-time broad-band seismic recordings to automatically monitor earthquake activities in the vicinity of Taiwan. The centroid moment tensor (CMT) inversion technique and a grid search scheme are applied to obtain the information of earthquake source parameters, including the event origin time, hypocentral location, moment magnitude and focal mechanism. All of these source parameters can be determined simultaneously within 117 s after the occurrence of an earthquake. The monitoring area involves the entire Taiwan Island and the offshore region, which covers the area of 119.3°E to 123.0°E and 21.0°N to 26.0°N, with a depth from 6 to 136 km. A 3-D grid system is implemented in the monitoring area with a uniform horizontal interval of 0.1° and a vertical interval of 10 km. The inversion procedure is based on a 1-D Green's function database calculated by the frequency–wavenumber (fk) method. We compare our results with the Central Weather Bureau (CWB) catalogue data for earthquakes occurred between 2010 and 2012. The average differences between event origin time and hypocentral location are less than 2 s and 10 km, respectively. The focal mechanisms determined by RMT are also comparable with the Broadband Array in Taiwan for Seismology (BATS) CMT solutions. These results indicate that the RMT system is realizable and efficient to monitor local seismic activities. In addition, the time needed to obtain all the point source parameters is reduced substantially compared to routine earthquake reports. By connecting RMT with a real-time online earthquake simulation (ROS) system, all the source parameters will be forwarded to the ROS to make the real-time earthquake simulation feasible. The RMT has operated offline (2010–2011) and online (since January 2012 to present) at the Institute of Earth Sciences (IES), Academia Sinica (http://rmt.earth.sinica.edu.tw). The long-term goal of this system is to provide real-time source information for rapid seismic hazard assessment during large earthquakes. [ABSTRACT FROM AUTHOR]

The rapid detection and characterization of megathrust earthquakes is a difficult task given their large rupture zone and duration. These events produce very strong ground vibrations in the near field that can cause weak motion instruments to clip, and they are also capable of generating large-scale tsunamis. The 2011 M 9 Tohoku-oki earthquake that occurred offshore Japan is one member of a series of great earthquakes for which extended geophysical observations are available. Here, we test an automated scanning algorithm for great earthquakes using continuous very long-period (100–200 s) seismic records from K-NET strong-motion seismograms of the earthquake. By continuously performing the cross-correlation of data and Green's functions (GFs) in a moment tensor analysis, we show that the algorithm automatically detects, locates and determines source parameters including the moment magnitude and mechanism of the great Tohoku-oki earthquake within 8 min of its origin time. The method does not saturate. We also show that quasi-finite-source GFs, which take into account the effects of a finite-source, in a single-point source moment tensor algorithm better fit the data, especially in the near-field. We show that this technique allows the correct characterization of the earthquake using a limited number of stations. This can yield information usable for tsunami early warning. [ABSTRACT FROM AUTHOR]

The vibration of a fluid-filled crack is considered to be one of the most plausible source mechanisms for the long-period events and volcanic tremors occurring around volcanoes. As a tool for the quantitative interpretation of source process of such volcanic seismic signals, we propose a method to numerically simulate the dynamic response of a fluid-filled crack. In this method, we formulate the motions of the fluid inside and the elastic solid outside of the crack, using boundary integrals in the frequency domain and solve the dynamic interactions between the fluid and the elastic solid using the point collocation method. The present method is more efficient compared with the time-domain finite difference method, which has been used in simulations of a fluid-filled crack and enables us to study the dynamics of a fluid-filled crack over a wide range of physical parameters. The method also allows us direct calculation of the attenuation quality factor of the crack resonance, which is an indispensable parameter for estimating the properties of the fluid inside the crack. The method is also designed to be flexible to many applications, which may be encountered in volcano seismology, and thus, extensions of the method to more complicated problems are promising. [ABSTRACT FROM AUTHOR]

Examines the broad-band recordings of deep events for the presence of signals from the lower-mantle discontinuities. Use of the novel technique; Important of the discontinuities for geodynamics, geochemistry and mineral physics; Details on the seismograms findings.

The time reversal operation in seismic source estimation is considered. We show that the time reversal operation, equally the adjoint operation, for seismic source imaging gives an approximate solution to more conventional seismic source inverse problem through the ‘happy approximation’ underlined by Claerbout. Practical applications of such methods in a long-period range to monitor earth's activities in realtime are also discussed. [ABSTRACT FROM AUTHOR]