%0 Article %J Geophysical Journal International %D 2020 %T Upper mantle discontinuities beneath Australia from trans-dimensional Bayesian inversions using multi-mode surface waves and receiver functions %A Taira, T. %A Yoshizawa, K. %1 10.1093/gji/ggaa442 %K Mineos %X Radially anisotropic S wave structures under the permanent seismic stations in Australia are reconstructed using multi-mode surface waves and receiver functions in a framework of the Bayesian inference. We have developed a fully nonlinear method of joint inversions incorporating P receiver functions and multi-mode Rayleigh and Love waves, based on the trans-dimensional Hierarchical Bayesian formulation. The method allows us to estimate a probabilistic Earth model taking account of the complexity and uncertainty of Earth structure, by treating the model parameters and data errors as unknowns. The Parallel Tempering algorithm is employed for the effective parameter search based on the reversible-jump Markov Chain Monte Carlo method. The use of higher modes enables us to enhance the sensitivity to the depth below the continental asthenosphere. Synthetic experiments indicate the importance of higher-mode surface waves for the better recovery of radial anisotropy in the whole depth range of the upper mantle. The method is applied to five Global Seismographic Network stations in Australia. While the S-wave models in eastern Australia show shallow Lithosphere-Asthenosphere Boundary (LAB) above 100Â~km depth, those in central and western Australia exhibit both Mid-Lithosphere Discontinuities (MLDs) and LAB. Also, seismic velocity jumps equivalent to the Lehmann Discontinuity (L-D) are found in all seismic stations. The L-D under the Australian continents is found at around 200 -- 300Â~km depth, depending on locations. Radial anisotropy in the depth range between LAB and L-D tends to show faster SH anomalies, which may indicate the effects of horizontal shear underneath the fast-moving Australian plate.