Preprint
Article

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

Altmetrics

Downloads

139

Views

192

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

16 October 2021

Posted:

19 October 2021

You are already at the latest version

Alerts
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: 
Subject: Environmental and Earth Sciences  -   Geophysics and Geology
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated