Preprint Article Version 1 Preserved in Portico This version is not peer-reviewed

On the Potential of 3D Transdimensional Surface Wave Tomography for Geothermal Prospecting of the Reykjanes peninsula

Version 1 : Received: 16 October 2021 / Approved: 19 October 2021 / Online: 19 October 2021 (08:23:56 CEST)

A peer-reviewed article of this Preprint also exists.

Rahimi Dalkhani, A.; Zhang, X.; Weemstra, C. On the Potential of 3D Transdimensional Surface Wave Tomography for Geothermal Prospecting of the Reykjanes Peninsula. Remote Sens. 2021, 13, 4929. Rahimi Dalkhani, A.; Zhang, X.; Weemstra, C. On the Potential of 3D Transdimensional Surface Wave Tomography for Geothermal Prospecting of the Reykjanes Peninsula. Remote Sens. 2021, 13, 4929.

Journal reference: Remote Sens. 2021, 13, 4929
DOI: 10.3390/rs13234929

Abstract

Seismic travel time tomography using surface waves is an effective tool for three-dimensional crustal imaging. Historically, these surface waves are the result of active seismic sources or earthquakes. More recently, however, also surface waves retrieved through the application of seismic interferometry are exploited. Conventionally, two-step inversion algorithms are employed to solve the tomographic inverse problem. That is, a first inversion results in frequency-dependent, two-dimensional maps of phase velocity, which then serve as input for a series of independent, one-dimensional frequency-to-depth inversions. As such, a two-dimensional grid of localized depth-dependent velocity profiles are obtained. Stitching these separate profiles together subsequently yields a three-dimensional velocity model. Relatively recently, a one-step three-dimensional non-linear tomographic algorithm has been proposed. The algorithm is rooted in a Bayesian framework using Markov chains with reversible jumps, and is referred to as transdimensional tomography. Specifically, the three-dimensional velocity field is parameterized by means of a polyhedral Voronoi tessellation. In this study, we investigate the potential of this algorithm for the purpose of recovering the three-dimensional surface-wave-velocity structure from ambient noise recorded on and around the Reykjanes Peninsula, southwest Iceland. To that end, we design a number of synthetic tests that take into account the station configuration of the Reykjanes seismic network. We find that the algorithm is able to recover the 3D velocity structure at various scales in areas where station density is high. In addition, we find that the standard deviation on the recovered velocities is low in those regions. At the same time, the velocity structure is less well recovered in parts of the peninsula sampled by fewer stations. This implies that the algorithm successfully adapts model resolution to the density of rays. Also, it adapts model resolution to the amount of noise on the travel times. Because the algorithm is computationally demanding, we modify the algorithm such that computational costs are reduced while sufficiently preserving non-linearity. We conclude that the algorithm can now be applied adequately to travel times extracted from (time-averaged) station-station cross correlations by the Reykjanes seismic network.

Keywords

Seismic interferometry; Transdimensional tomography; Surface wave dispersion; probabilistic inversion; Markov chain Monte Carlo

Subject

EARTH SCIENCES, Geophysics

Comments (0)

We encourage comments and feedback from a broad range of readers. See criteria for comments and our diversity statement.

Leave a public comment
Send a private comment to the author(s)
Views 0
Downloads 0
Comments 0
Metrics 0


×
Alerts
Notify me about updates to this article or when a peer-reviewed version is published.
We use cookies on our website to ensure you get the best experience.
Read more about our cookies here.