High resolution Vp and Vs imaging for geothermal exploration in the Ilan Plain, northeastern Taiwan by joint inversion of body-wave and surface-wave data

The Ilan Plain in northeastern Taiwan located at the southern termination of the Okinawa Trough, which is formed by the back-arc extension of the Ryukyu subduction zone. Owing to this extensional tectonic setting, the Ilan Plain shows distributed hot springs with high geothermal gradient in several places, such as the Chingshui, Jiaoxi, Longde-Lize and Hongchailin areas. These are potential locations for geothermal exploration which require a thorough investigation to understand their geothermal system and to evaluate their potential in geothermal power generation. Seismic velocity provides key information to determine properties of a geothermal system including the distribution of fluids, fractures, and porosity. This study obtained seismic velocity images of the Ilan Plain and the surrounding area by joint inversion of body wave travel times of local earthquakes and Rayleigh wave phase velocities from ambient noise correlations recorded by a temporary dense geophone array in the Ilan Plain and the Formosa Array in northeastern Taiwan. The joint inversion takes advantage of the complementary sensitivities of both data types to better constrain both shallow and deep structures. In addition, we also lift the flat earth surface assumption commonly assumed in the surface wave tomography to avoid possible biased velocity structure due to topography by adapting the surface wave sensitivity kernel to the surface topography in the inversion. The obtained image at the 1-km depth shows low Vp and Vs corresponding to unconsolidated Quaternary sediments. The deeper images at 2- to 5-km depths show a low-velocity corridor extending from the Guishan volcanic island northeastern offshore to the southwest corner of the Ilan Plain, implying that the heat source of the geothermal systems in the Ilan Plain may be related to magmatic intrusions of backarc rifting. Locally, in the shallow region of the Chingshui and Jiaoxi area and the deeper region of the Hongchailin and Longde-Lize area, the high Vp/Vs ratio with little increase in Vp and significant decrease in Vs could indicate a saturated rock with high porosity, based on the rock physics principles in geothermal environment summarized in previous study. The low Vp, Vs, and Vp/Vs ratio in the shallow region of the Hongchailin area could indicate vapor-rich high temperature environment. Our results demonstrate the potential of new seismic imaging integrating body- and surface-wave data for geothermal exploration.