Search

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]

Abstract: Oceanic plates subducted at trenches penetrate into the deep mantle, and encounter a structural boundary at a depth of 410km where olivine, the dominant element of mantle rocks, transforms into a higher density form wadsleyite. This transformation may be delayed within the coldest core of subducting plates (slabs) due to kinetic effects, and it has been suggested that metastable olivine may persist deeper than 410km. Using high density seismic array data in Japan, we show the direct image of the structure corresponding to this metastable olivine wedge (MOW) beneath southwest Japan. Numerical simulation of a subducting slab, including the kinetic effect of water (H2O) on the olivine–wadsleyite transformation, indicates that the presence of the imaged MOW requires an insignificant amount of water (less than 100wt. ppm) be present in the slab mantle, thus a deep dry cold slab. We infer that the transportation of water into the deep mantle occurs along the top surface of the subducting slab, but no significant amount within the slab itself. We also demonstrate that a numerical simulation including the kinetics of 660-km phase transformation can reconcile the observed deep depression of the 660-km discontinuity with a gentle Clapeyron slope. [Copyright &y& Elsevier]

Abstract: The petrological consequences of deep subhorizontal deflection (“stagnation”) of subducting slabs should affect both apparent seismic velocity structures and slab morphology. We construct kinematic thermal models of stagnant slabs and perform thermodynamic modeling of the consequent perturbation of high-pressure phase transitions in mantle minerals, focusing upon Japan as our study area. We calculate associated thermo-petrological buoyancy forces and bending moments which (along with other factors such as viscosity variations and rollback dynamics) may contribute to slab deformation. We consider effects of variations in depth of stagnation, post-stagnation dip angle, phase transition sharpness, transition triplication due to multiple intersection of geotherms with phase boundaries, and potential persistence of metastable phases due to kinetic hindrance. We also estimate seismic velocity anomalies, as might be imaged by seismic tomography, and corresponding seismic velocity gradients, as might be imaged by receiver-function analysis. We find that buoyant bending moment gradients of petrological origin at the base of the transition zone may contribute to slab stagnation. Such buoyancy forces vary with the depth at which stagnation occurs, so that slabs may seek an equilibrium slab stagnation depth. Metastable phase bending moment gradients further enhance slab stagnation, but they thermally decay after ∼600–700km of horizontal travel, potentially allowing stagnant slabs to descend into the lower mantle. Stagnant slabs superimpose zones of negative velocity gradient onto a depressed 660-km seismic discontinuity, affecting the seismological visibility of such features. Seismologically resolvable details should depend upon both stagnation depth and the nature of the imaging technique (travel-time tomography vs. boundary-interaction phases). While seismic tomography appears to yield images of stagnant slabs, discontinuity topography beneath Japan resolved by migrated receiver functions appears to be consistent with slab penetration of the transition zone. However, model slabs which bottom around ∼780–810km and then bend upwards by a few degrees can match both the tomographic and receiver-function images. [Copyright &y& Elsevier]

We report seismic evidence for the transportation of water into the deep mantle in the subduction zone beneath northeastern Japan. Our data indicate that water is released from the hydrated oceanic crust at shallow depths (≤100 kilometers) and then forms a channel of hydrated mantle material on top of the subducting plate that is the pathway for water into the deep mantle. Our result provides direct evidence that shows how water is transported from the ocean to the deep mantle in a cold subduction zone environment. [ABSTRACT FROM AUTHOR]

Discusses the discovery of a discontinuity in the mantle beneath the subduction zones in Tonga and the Japan and Flores seas. Data analysis using of the J-array of Japan; Identification of 920-kilometer discontinuity; Speculation as a global feature; Significance.

High frequency energy radiation of large earthquakes is a key to evaluating shaking damage and is an important source characteristic for understanding rupture dynamics. We developed a new inversion method, Image Deconvolution Back-Projection (IDBP) to retrieve high frequency energy radiation of seismic sources by linear inversion of observed images from a back-projection approach. The observed back-projection image for multiple sources is considered as a convolution of the image of the true radiated energy and the array response for a point source. The array response that spreads energy both in space and time is evaluated by using data of a smaller reference earthquake that can be assumed to be a point source. The synthetic test of the method shows that the spatial and temporal resolution of the source is much better than that for the conventional back-projection method. We applied this new method to the 2001 M w 7.8 Kunlun earthquake using data recorded by Hi-net in Japan. The new method resolves a sharp image of the high frequency energy radiation with a significant portion of supershear rupture. [ABSTRACT FROM AUTHOR]

Abstract: At Aso volcano, Kyushu, Japan, several different types of volcanic tremor have been observed for many years. One of them is the continuous tremor, the ground vibration which has dominant frequency between 3 and 10 Hz and has approximately constant amplitudes without any clear beginning and ending. We observed the continuous tremor at Aso using short period seismometer arrays for 3 days in 1999. We locate the source of the continuous tremor by seismic array data processing. We use the semblance coefficients in order to estimate the arrival azimuth and apparent slowness by grid search. The epicenters of the continuous tremor are located around the currently active crater, and the source depths are likely to be shallower than 600 m. We find that the estimated epicenter clearly migrates synchronized with the change in the tremor amplitude. The migration often occurs periodically with a period of about 80 s, but aperiodic occurrence of the migration is also often seen. In both cases, the epicenter is located southeastward (northwestward) when the amplitude is larger (smaller). We propose that there are two or more independent tremor sources with fixed locations, and that their amplitudes modulate either aperiodically or periodically with periods nearly 80 s. The tremor signals from those sources are mixed at the arrays, and the estimated epicenter parameters vary according to which signal dominates the seismograms. The simplest model is that there are two point sources, one at west and the other at south of the crater (we call the two sources as “NW source” and “SE source”, respectively), and the amplitude of the SE source changes with time. Consequently, the SE source dominates the seismogram when the observed amplitude is larger, whereas the NW source dominates when the amplitude is smaller. We generate synthetic seismograms, and apply the location technique to them to verify the validity of the two point source model and to search the locations of two sources which can explain the observed synchronization between the amplitude and the apparent epicenter location. We find that the distance between the two sources needs to be more than 400 m to agree with the observation. We also analyze the seismic array data observed in 2001, and infer that the NW source in 1999 may be identical to the tremor source of the 2001 data. [Copyright &y& Elsevier]

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]

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]

Abstract: In order to investigate the continuous volcanic tremor of Aso Volcano in Japan, we performed a series of temporary short-period seismic array observations near the Nakadake first crater (the active crater) during five years from 1999 to 2003. We deployed in all of the temporary observations a seismic array at the same location about 700 m west of the active crater, in order to investigate long-term changes in the tremor activity. In 1999 and 2003, another array was simultaneously deployed at a different location 700 m north of the crater to help locate the tremor sources. We developed a frequency domain semblance method and applied it to the waveform data of the frequency range where the continuous tremor is dominant (3–6 Hz). We measured arrival azimuths and slownesses of the continuous tremor signals as functions of frequency, which are then used to locate the epicenters of the tremor signals corresponding to the principal peaks of the power spectra. For the observations in 1999 and 2002, the continuous tremor amplitudes are relatively small, and the slowness of the tremor signal observed at the west array takes a local minimum (0.5 to 0.6 s/km) near the frequency (∼4.7 Hz for 1999 and ∼4.8 Hz for 2002) which corresponds to the highest spectral peak. This implies that body waves dominate the tremor signals at the west array around the frequency. The tremor epicenters corresponding to 4.7 Hz for the observation in 1999 are located at the west rim of the currently active crater. While the surface crater activity of Aso remains low and the tremor activity is not clearly linked with the surface activity until early 2003, a close link between the tremor and crater activity appears in the middle of 2003, when a small phreatic eruption occurred a month before the array observation (July 10, 2003). Tremor signals of the observation in 2003 show a large spectral peak (4.2 Hz) where the slowness measured for the west array is very large (1.1 s/km), clearly suggesting that surface waves are dominant. The epicenter is again located at the western rim of the active crater. We interpret these observations as follows: in 1999 and 2002 when the surface activities of Aso were low, the continuous tremor excitation was deep and inactive. In the middle of 2003 when Aso Volcano became active with a series of phreatic eruptions a shallower tremor source was activated, possibly masking the deeper sources. This shallowing of the dominant tremor source could be due to the increase in the volcanic gas flow rate triggered by the phreatic eruptions. [Copyright &y& Elsevier]