Detection of small earthquakes along the Pacific-Antarctic Ridge from T-waves recorded by abyssal ocean-bottom observatories

A number of seismoacoustic T-wave events were observed between January 2003 and January 2004 by broadband ocean-bottom seismometers installed on the French Polynesia seafloor at depths of 4,000?C5,000?m, well below the conjugate depth of the SOFAR channel. Using T-wave arrival times, we located 89 T-wave events along the Pacific-Antarctic Ridge. Among these, 63 events were not detected by land-based seismic observations of the United States Geological Survey (USGS), which was nearly twice the number of earthquakes reported by the USGS in the area during the observation period. We used a simple method to estimate earthquake magnitude from T-wave energy. The magnitudes of the events unidentified by the USGS ranged from 2.3 to 5.5, whereas magnitudes of the earthquakes reported by the USGS ranged from 4.1 to 6.2. Our study suggests that T-wave observations with abyssal ocean-bottom seismometers can improve the detection of small earthquakes and help our understanding of the seismotectonics of remote oceanic areas.     

APPLICATION OF MEM SPECTRAL ESTIMATION TO MU RADAR OBSERVATION

Preliminary results of the wind velocity estimation using the Maximum Entropy Method (MEM) to MUradar observation data sets are presented. The comparison of the results from the periodogram method and theMEM shows that the MEM estimation is reliable, and has higher accuracy, resolution and detectability than theestimation from periodogram method. The high accuracy power spectrum obtained by the MEM is veryuseful to studying the atmospheric turbulence structure. However. the MEM needs the longer computingtime for obtaining the high accuracy spectrum. Particularly, the estimation of MEM will bring serious devia-tion at lower signal-to-noise ratio.     

Iterative deconvolution of complex body waves from great earthquakes—the Tokachi-Oki earthquake of 1968

A method of body-wave inversion is developed in an attempt to extract the information about asperities or barriers in a fault zone. A sequence of point sources, each being characterized with the seismic moment, the onset time and the location, are iteratively derived from observed records at multi-stations, where the two-dimensional extent of the source location is taken into account. A modification is made of the iterative method of Kikuchi and Kanamori on the formulation of inversion procedure to facilitate the computation.Using this method, we analyse long period P waves of the Tokachi-Oki earthquake of 1968 (M_w = 8.2) and obtain several significant subevents with time durations of ～ 10 s. Their spatio-temporal distribution shows that the rupture process consists of three characteristic stages: (A) a stage of introductory rupture, (B) a stage of main rupture and (C) a stage of aftershocks. The main rupture takes place in the form of clustering around a few sites of the fault plane. The largest subevent occurs in the northwestern corner. The stress drop associated with this event is estimated to be ～ 200 bars, one order of magnitude higher than the stress drop averaged over the entire fault plane. The sum of the seismic moments of the individual subevents amounts to 2.3 × 10²⁸ dyn. cm which approximately coincides with the one estimated from the analysis of long-period surface waves. This implies that the source of the Tokachi-Oki earthquake consists of several major subevents with time durations of ～ 10 s in addition to other minor subevents.     

Full-correlation analysis of turbulent scattering layers in the mesosphere observed by the MU radar

We have applied a full-correlation analysis technique to the echo power fluctuations observed by the MU radar (35°N, 136°E), and analyzed the horizontal structure of the scattering pattern in the mesosphere as well as their horizontal motions. The velocity of the scattering pattern did not agree with the background wind velocity, but was associated with the horizontal propagating direction of a saturated inertia gravity wave identified in the wind field. The length of the long axis of the characteristic ellipse of the scattering pattern was approximately 50 km, and the direction was almost perpendicular to the propagating direction of the wave. The correlation time of the scattering pattern was approximately 700 s, which is much longer than the lifetime of the isolated turbulence itself. This implies that the observed scattering pattern is associated with a region where the saturated inertia gravity wave generates turbulence.     

Diurnal variations of the planetary boundary layer observed with anL-band clear-air doppler radar

Based on continuous observations of the planetary boundary layer (PBL) with anL-band (1357.5 MHz) boundary-layer radar (BLR) at a hilly location in Japan, we have discovered that on clear days, a thin enhanced echo layer corresponding to the top of the PBL (or mixed layer) appeared at about 500 m height in the morning and ascended to about 1500 m in the afternoon. Strong upward velocities were observed below the echo layer (or inside the PBL), reaching 1500 m in the afternoon.     

Topography of the mantle discontinuities beneath the South Pacific superswell as inferred from broadband waveforms on seafloor

We analyzed seismic waveforms recorded by the broadband ocean bottom seismographs deployed in the South Pacific superswell to determine the depths of the mantle discontinuities using a receiver function method. We estimated the thickness of the mantle transition zone (MTZ) to be 245 km on average beneath the superswell region, which is close to the global average. The MTZ is found to be thinned locally beneath the Society hot spot by 30 km. Temperature anomalies computed from the thinned transition zone and the Clapeyron slope of the olivine phase transformations are 150–200 K beneath the Society hot spot. Previous studies, using land-based data, suggested the presence of a hot MTZ beneath the Pitcairn hot spot. The locally hot transition zone beneath the hot spots is attributed to narrow mantle plumes probably rising from the lower mantle. The normal average thickness obtained from the present study indicates that there is no broad upwelling across the transition zone beneath the entire superswell region at present.     

The height of the maximum ionospheric electron density over the MU radar

Ionospheric F₂-layer peak height h_mF₂ variations, as measured over 1986–1995 by the MU radar (34.85°N, 136.1°E) and as calculated with a theoretical model, are discussed. The diurnal variations of the measured peak height for different seasons and levels of solar activity are compared with those estimated from ionosonde M3000F₂ and IRI predictions. Also given are the measured ion drift velocities and meridional neutral winds needed to understand the dynamic behavior of the F₂-layer. It is found that: (1) h_mF₂ is generally higher during periods of the solar maximum than during periods of the solar minimum, and higher in summer than in winter; (2) for the solar maximum, h_mF₂ drops markedly in the morning and in the afternoon, while, for the solar minimum, the h_mF₂ minimum occurs in the morning during summer and usually in the afternoon during winter. In general, the measured h_mF₂ is well reproduced by our model when we use the observed drift velocities and plasma temperatures as inputs. Our modeling study shows that the neutral wind contributes strongly to the diurnal variation of h_mF₂ in winter by lowering the ionization layer by day, particularly for the solar maximum; it also helps to enlarge the day–night difference of h_mF₂ in summer. The northward electromagnetic drifts that usually cancel the neutral wind effect have only a minor effect for the location of the MU radar. Other features of the observed h_mF₂ variations, e.g., the solar maximum–minimum difference, the summer–winter difference, and the morning and afternoon drops, are explained by the basic processes of O⁺ production, loss and diffusion, as influenced by the atomic oxygen concentration and neutral and plasma temperatures.     

Signature of remnant slabs in the North Pacific from P-wave tomography

A 3-D ray-tracing technique was used in a global tomographic inversion in order to obtain tomographic images of the North Pacific. The data reported by the Geophysical Survey of Russia (1955–1997) were used together with the catalogues of the International Seismological Center (1964–1991) and the US Geological Survey National Earthquake Information Center (1991–1998), and the recompiled catalogue was reprocessed. The final data set, used for following the inversion, contained 523 430 summary ray paths. The whole of the Earth's mantle was parametrized by cells of 2° × 2° and 19 layers. The large and sparse system of observation equations was solved using an iterative LSQR algorithm.
A subhorizontal high-velocity anomaly is revealed just above the 660 km discontinuity beneath the Aleutian subduction zone. This high-velocity feature is observed at latitudes of up to ~70°N and is interpreted as a remnant of the subducted Kula plate, which disappeared through ridge subduction at about 48 Ma. A further positive velocity perturbation feature can be identified beneath the Chukotka peninsula and Okhotsk Sea, extending from ~300 to ~660 km depth and then either extending further down to ~800 km (Chukotka) or deflecting along the 660 km discontinuity (Okhotsk Sea). This high-velocity anomaly is interpreted as a remnant slab of the Okhotsk plate accreted to Siberia at ~55 Ma.     

Incessant excitation of the Earth's free oscillations: global comparison of superconducting gravimeter records

Records of superconducting gravimeters (SGs) at Canberra (Australia), Esashi (Japan), Metsähovi (Finland) and Syowa Stations (Antarctica) were analyzed to search for further evidence of background free oscillations of the Earth. Spectrograms for 1-year period and averaged power spectra for seismically quiet periods were obtained for each of the stations. Anomalous features of the oscillations observed at Syowa Station, such as an apparent seasonal variation and a high intensity at frequencies between 3 and 4 mHz, were absent at the other SG stations. Among the SG stations used in this study, the background free oscillations were detected most consistently and distinctly at Canberra, where the noise level was comparable to that at the IDA quietest station, while that at Syowa Station was close to the critical limit for detecting the oscillations. The background free oscillations provide a good reference to evaluate the noise level in the milliHertz band.