We studied the seismic velocity structure beneath the Krafla central volcano, NE Iceland, by performing 3‐D tomographic inversions of 1453 earthquakes recorded by a temporary local seismic network between 2009 and 2012. The seismicity is concentrated primarily around the Leirhnjúkur geothermal field near the center of the Krafla caldera. To obtain robust velocity models, we incorporated active seismic data from previous surveys. The Krafla central volcano has a relatively complex velocity structure with higher P wave velocities (Vp) underneath regions of higher topographic relief and two distinct low‐Vp anomalies beneath the Leirhnjúkur geothermal field. The latter match well with two attenuating bodies inferred from S wave shadows during the Krafla rifting episode of 1974–1985. Within the Leirhnjúkur geothermalreservoir, we resolved a shallow (−0.5 to 0.5 km below sea level; bsl) region with low‐Vp/Vs values and a deeper (0.5–1.5 km bsl) high‐Vp/Vs zone. We interpret the difference in the velocity ratios of the two zones to be caused by higher rock porosities and crack densities in the shallow region and lower porosities and crack densities in the deeper region. A strong low‐Vp/Vs anomaly underlies these zones, where a superheated steam zone within felsic rock overlies rhyolitic melt.

We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective QS/QP ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small.

Seismicity was monitored beneath the Krafla central volcano, NE Iceland, between 2009 and 2012 during a period of volcanic quiescence, when most earthquakes occurred within the shallow geothermal field. The highest concentration of earthquakes is located close to the rock‐melt transition zone as the Iceland Deep Drilling Project‐1 (IDDP‐1) wellbore suggests and decays quickly at greater depths. We recorded multiple swarms of microearthquakes, which coincide often with periods of changes in geothermal field operations, and found that about one third of the total number of earthquakes are repeating events. The event size distribution, evaluated within the central caldera, indicates average crustal values with b = 0.79 ± 0.04. No significant spatial b value contrasts are resolved within the geothermal field nor in the vicinity of the drilled melt. Besides the seismicity analysis, focal mechanisms are calculated for 342 events. Most of these short‐period events have source radiation patterns consistent with double‐couple (DC) mechanisms. A few events are attributed to non‐shear‐faulting mechanisms with geothermal fluids likely playing an important role in their source processes. Diverse faulting styles are inferred from DC events, but normal faulting prevails in the central caldera. The best fitting compressional and tensional axes of DC mechanisms are interpreted in terms of the principal stress or deformation rate orientations across the plate boundary rift. Maximum compressive stress directions are near‐vertically aligned in different study volumes, as expected in an extensional tectonic setting. Beneath the natural geothermal fields, the least compressive stress axis is found to align with the regional spreading direction. In the main geothermal field both horizontal stresses appear to have similar magnitudes causing a diversity of focal mechanisms.

No abstract