3D S-wave velocity imaging of Reykjanes Peninsula high-enthalpy geothermal fields with ambient-noise tomography

• Published in: Journal of Volcanology and Geothermal Research
• Main Authors: Martins, Joana E. (author), Weemstra, C. (author), Ruigrok, Elmer (author), Verdel, Arie (author), Jousset, Philippe (author), Hersir, Gylfi (author)
• Format: Article in Journal/Newspaper
• Language: English
• Published: 2019
• Subjects: ambient noise tomography (ANT); Volcano; Reservoir characterization; Seismic interferometry; Surface-waves; Velocity anomalies; Reykjanes Peninsula; magma chamber; magma plumbing; magma mush; Beamforming; Geothermal; high-enthalpy geothermal reservoir; monitoring; passive seismic; Reykjanes; Iceland
• Online Access: http://resolver.tudelft.nl/uuid:7c5cdb15-80fe-4dc6-8870-a45d705a73ef; https://doi.org/10.1016/j.jvolgeores.2019.106685

Tomographic imaging based on ambient seismic noise measurements has shown to be a powerful tool, especially in areas like Iceland, where the microseism illumination is excellent. In this paper, we produce a 3D S-wave tomographic image over the western Reykjanes Peninsula high-enthalpy geothermal fields and evaluate the reliability of the tomographic results for different resolutions through simulated and real data. We use 30 broadband stations operating for approximately one-and-a-half year and apply ambient noise seismic interferometry for each station-pair. This results in empirical Green's functions in which especially the ballistic surface waves (BSW) are well resolved. The retrieved BSW exhibit a high signal-to-noise ratio between 0.1 and 0.5 Hz, and the beamforming analysis indicates an apparent surface-wave velocity of 3 km/s over a broad azimuthal range. For the tomographic inversion, we invert the estimated phase velocities between all station pairs to frequency-dependent phase velocity maps in four different resolutions (1, 2, 3, and 4 km) using a Tikhonov regularisation. With the estimated regularisation parameter per frequency per resolution, we invert simulated data for checkerboard sensitivity tests per frequency for different combinations of velocity anomaly sizes and resolutions. Finally, after the inversion to depth, we detect S-wave velocity anomalies with variations between −15% and 15% with reference to an estimated average velocity using 1 km and 3 km of lateral resolutions and 1 km of vertical resolution. This study shows the potential of ambient noise tomography as complementary seismological tool for reservoir characterization. Invited Research Article Mathematical Geodesy and Positioning Applied Geophysics and Petrophysics