# [SIT28-P10] Preliminary study using harmonic decomposition of receiver functions to image non-isotropic characteristics of Japan subduction zone

Keywords:Receiver Function, Harmonic decomposition, Anisotropy, Subduction zones

Conventional receiver function methods assume horizontal geometry and isotropy for velocity discontinuity analysis. However, in subduction zones, the isotropy assumption unlikely holds. In case of anisotropic velocity or dipping velocity discontinuity P-to-S receiver function show variation by back azimuth (Shiomi & Park, 2008). We employed the harmonic decomposition method (Bianchi et al. 2010, Agostinetti & Miller, 2014) to extract non-isotropic component from receiver functions to image the Pacific plate subduction under Japan. The harmonic decomposition gives five components (isotropic, cos(kθ), sin(kθ) terms for k=1, 2) from linear matrix inversion using radial and transverse receiver functions. The preliminary analysis using data from three Hi-net stations (ANIH, IHEH, KZMH) located along the 40N latitude with varying distance from trench shows following features: (1) Within first harmonics, the EW (sin(θ)) component is larger than the NS (cos(θ)) component at timing of the oceanic Moho phase. This well reproduces westward dipping of the Pacific slab. (2) In upper most crust (0-1 s), amplitude of k=2 harmonics is larger than k=1 harmonics, which implies large horizontal symmetric axis anisotropy. (3) Between continental Moho and top of subducting oceanic crust, large k=1 and k=2 harmonics are observed in mantle wedge . (4) Below subducting oceanic crust, both k=1 and k=2 harmonic components decrease consistently for all three stations, but locally large k=1 harmonics appear. Signature of previously reported hydrated mantle above subducting oceanic crust (Kawakatsu & Watada, 2007) is observed in station ANIH. At later positive peak, large amplitude of k=1 and k=2 harmonics is observed, which might indicate existence of dipping structure having horizontal symmetric anisotropy beneath. Our results show a possibility of applying the harmonic decomposition method to image non-isotropic component of subduction zones using receiver functions.