Finite element modelling of surface wave transmission across regions of subduction

Abstract

The propagation of Rayleigh waves across the regions of subduction of Japan and of the North Island of New Zealand is examined in the frequency domain by the finite element method. Study of the Japan region of subduction ESE of the island of Honshu has shown that, for perfectly elastic 2D models, although Love waves decrease in velocity across the region, Rayleigh waves of periods of 35-60s increase slightly in velocity. The material properties of the predominantly crustal material being subducted appear to slow the Love waves, whereas the Rayleigh waves of periods of from 35 to 60 s tend to retain their oceanic velocities, which are higher than the velocities of Rayleigh waves of those periods for the island of Honshu. Amplitude calculations have shown that, at periods of 60 s down to 10 s, Love waves are increasingly forward scattered into higher Love modes, whereas, at a period of 20s, 98 per cent of the energy of Rayleigh waves is transmitted as the Rayleigh fundamental mode. This result, together with that of Drake for the continental margin at Berkeley, California, confirms the practice of using Rayleigh waves at a period of approximately 20 s measured at seismographic stations near continental margins and near regions of subduction to estimate the magnitude of teleseisms. Study of the New Zealand region of subduction SE of Lake Taupo on the North Island of New Zealand has shown similar results to those for the Japan region of subduction. Although Love waves decrease in velocity across the region, Rayleigh waves of periods of 60 s down to 20 s increase slightly in velocity. Amplitude calculations have shown that at periods of 60 s down to 15 s, Love waves are increasingly forward scattered into higher Love modes. On the other hand, at a period of 20 s, 96 per cent of the energy of Rayleigh waves is transmitted as the Rayleigh fundamental mode. At a period of 10 s, in contrast to the results for the Japan region of subduction and for the continental margin at Berkeley, 73 per cent of the energy of Love waves is transmitted as the Love fundamental mode and 97 per cent of the energy of Rayleigh waves is transmitted as the Rayleigh fundamental mode. Amplitude and phase changes of Love and Rayleigh waves at continental margins and in regions of subduction need to be allowed for when waveform inversions of waves which have traversed these regions are made.

Approximation of surface wave mode conversion at a passive continental margin by a mode-matching technique

“…With decreasing depth of the lateral heterogeneity, f R min and f L min increase (Fig. This was ¢rst noted by Lysmer & Drake (1971) and Drake (1972) for a frequency of 50 Hz and later con¢rmed for passive margins by Gregersen (1978), Drake & Bolt (1980), Schlue (1981), Badal & Seron (1987) and Fitas & Mendes-Victor (1992) and for subduction zones by Drake & Bolt (1989). The percentage of energy lost by the fundamental mode varies from less than 1 per cent for a lateral heterogeneity con¢ned to the asthenosphere to 87 per cent for a lateral heterogeneity in the upper crust (Table 1).…”

Beaufort Ambient Seismo-Acoustics Beneath Ice Cover (BASIC)

“…or coupling between Rayleigh modes at me continental shelf. This problem has been extensively studied, but only at frequencies below the microseism peak (Drake and Bolt, 1989). It may be feasible to calculate coupling coefficients for down slope propagation of Rayleigh waves .…”

Guided wave propagation across sharp lateral heterogeneities: the complete wavefield at plane vertical discontinuities

“…There already exist various publications on related problems, e.g. (arbitrary selection) (Alsop (1966), Drake & Bolt (1989), Its & Yanovskaya (1985), McGarr & Alsop (1967), Malischewsky (1987) andVaccari et al (1989). In all papers, the boundary conditions at the vertical interfaces are not satisfied exactly because of the neglect of either body waves, non-propagating modes or conversion between Love and Rayleigh modes.…”