Search

Katla is one of the most active volcanoes in Iceland and is characterised by persistent seismicity. It is partly covered by the Mýrdalsjökull glacier and its historic activity is dominated by phreatomagmatic eruptions within the caldera associated with catastrophic glacial floods. In July 2011 a sudden jökulhlaup was released from the glacier, associated with tremor, elevated seismicity inside the caldera and a new cluster of seismicity on the south flank. This was likely caused by a hydrothermal or magmatic event, possibly a small subglacial eruption. Similar unrests occurred in 1955 and 1999. We have identified changes of the seismicity pattern coinciding with the 2011 unrest, suggesting a modification in the volcanic system. It may be speculated that if the persistent seismicity at Katla is an indication of a pressurized magma system ready to erupt, small events like those of 1955, 1999 and 2011 may trigger larger eruptions in the future. We have also conducted a pilot study of the geology of the southern flank, where the new seismicity is recorded, and identified sources for flank eruptions in the recent eruptive history of Katla. These include rhyolitic domes and surtseyan craters. Therefore, a wide range of volcanic processes have to be taken into account as possible source for the new seismicity and volcanic hazard.

SUMMARYThe vibration modes of the ground have been described both in the 1-D and 2-D case. The 1-D resonance is found on geological structures whose aspect ratio is low, that is on layers with a lateral width much larger than their thickness. A typical example is that of a horizontal soft sediment layer overlying hard bedrock. In this case, the 1-D resonance frequency, traditionally detected by means of the microtremor H/V (horizontal to vertical spectral ratio) technique, depends on the bedrock depth and on the shear wave velocity of the resonating cover layer. The H/V technique is thus used both to map the resonance frequencies in seismic microzonation studies and for stratigraphic imaging. When 2-D resonance occurs, generally on deep and narrow valleys, the whole sedimentary infill vibrates at the same frequency and stratigraphic imaging can no longer be performed by means of the 1-D resonance equation. Understanding the 1-D or 2-D resonance nature of a site is therefore mandatory to avoid wrong stratigraphic and dynamic interpretations, which is in turn extremely relevant for seismic site response assessment. In this paper, we suggest a procedure to address this issue using single-station approaches, which are much more common compared to the multistation synchronized approach presented by research teams in earlier descriptions of the 2-D resonances. We apply the procedure to the Bolzano sedimentary basin in Northern Italy, which lies at the junction of three valleys, for which we observed respectively 1-D-only, 1-D and 2-D, and 2-D-only resonances. We conclude by proposing a workflow scheme to conduct experimental measurements and data analysis in order to assess the 1-D or 2-D resonance nature of a site using a single-station approach.

A 23-hour tremor burst was recorded on July 8–9th 2011 at the Katla subglacial volcano, one of the most active and hazardous volcanoes in Iceland. This was associated with deepening of cauldrons on the ice cap and a glacial flood that caused damage to infrastructure. Increased earthquake activity within the caldera started a few days before and lasted for months afterwards and new seismic activity started on the southern flank. No visible eruption broke the ice and the question arose as to whether this episode relates to a minor subglacial eruption with the tremor being generated by volcanic processes, or by the flood. The tremor signal consisted of bursts with varying amplitude and duration. We have identified and described three different tremor phases, based on amplitude and frequency features. A tremor phase associated with the flood was recorded only at stations closest to the river that flooded, correlating in time with rising water level observed at gauging stations. Using back-projection of double cross-correlations, two other phases have been located near the active ice cauldrons and are interpreted to be caused by volcanic or hydrothermal processes. The greatly increased seismicity and evidence of rapid melting of the glacier may be explained by a minor sub-glacial eruption. A less plausible interpretation is that the tremor was generated by hydrothermal boiling and/or explosions with no magma involved. This may have been induced by pressure drop triggered by the release of water when the glacial flood started. All interpretations require an increase of heat released by the volcano.

We introduce a back-projection method to locate tremor sources using products of cross-correlation envelopes of time series between seismic stations. For a given subset of n stations, we calculate ...

A double-correlation method is introduced to locate tremor sources based on stacks of complex, doubly-correlated tremor records of multiple triplets of seismographs back projected to hypothetical source locations in a geographic grid. Peaks in the resulting stack of moduli are inferred source locations. The stack of the moduli is a robust measure of energy radiated from a point source or point sources even when the velocity information is imprecise. Application to real data shows how double correlation focuses the source mapping compared to the common single correlation approach. Synthetic tests demonstrate the robustness of the method and its resolution limitations which are controlled by the station geometry, the finite frequency of the signal, the quality of the used velocity information and noise level. Both random noise and signal or noise correlated at time shifts that are inconsistent with the assumed velocity structure can be effectively suppressed. Assuming a surface-wave velocity, we can constrain the source location even if the surface-wave component does not dominate. The method can also in principle be used with body waves in three dimensions, although this requires more data and seismographs placed near the source for depth resolution.

Relative locationmethods are commonly used to precisely locate earthquake clusters consisting of similar waveforms. Repeating waveforms are often recorded at volcanoes, where, however, the crust structure is expected to contain strong heterogeneities and therefore the 1-D velocity model assumption that is made in most location strategies is not likely to describe reality. A peculiar cluster of repeating low-frequency seismic events was recorded on the south flank of Katla volcano (Iceland) from 2011. As the hypocentres are located at the rim of the glacier, the seismicity may be due to volcanic or glacial processes. Information on the size and shape of the cluster may help constraining the source process. The extreme similarity of waveforms points to a very small spatial distribution of hypocentres. In order to extract meaningful information about size and shape of the cluster, we minimize uncertainty by optimizing the cross-correlation measurements and relative-location process.With a synthetic test we determine the best parameters for differential-time measurements and estimate their uncertainties, specifically for each waveform. We design a location strategy to work without a pre-defined velocity model, by formulating and inverting the problem to seek changes in both location and slowness, thus accounting for azimuth, take-off angles and velocity deviations from a 1-D model. We solve the inversion explicitly in order to propagate data errors through the calculation. With this approach we are able to resolve a source volume few tens of metres wide in horizontal directions and around 100 metres in depth. There is no suggestion that the hypocentres lie on a single fault plane and the depth distribution indicates that their source is unlikely to be related to glacial processes as the ice thickness is not expected to exceed few tens of metres in the source area. Our method is designed for a very small source region, allowing us to assume a constant slowness for the whole cluster and to include the effects of 3-D heterogeneity such as refraction. Similar circumstances may arise in other volcanic regions with a high level of heterogeneity and where densely clustered earthquakes are often recorded.

Katla is a threatening volcano in Iceland, partly covered by the M\'yrdalsj\"okull ice cap. The volcano has a large caldera with several active geothermal areas. A peculiar cluster of long-period seismic events started on Katla's south flank in July 2011, during an unrest episode in the caldera that culminated in a glacier outburst. The seismic events were tightly clustered at shallow depth in the Gvendarfell area, 4 km south of the caldera, under a small glacier stream on the southern margin of M\'yrdalsj\"okull. No seismic events were known to have occurred in this area before. The most striking feature of this seismic cluster is its temporal pattern, characterized by regular intervals between repeating seismic events, modulated by a seasonal variation. Remarkable is also the stability of both the time and waveform features over a long time period, around 3.5 years. No comparable examples have been found in the literature. Both volcanic and glacial processes can produce similar waveforms and therefore have to be considered as potential seismic sources. Discerning between these two causes is critical for monitoring glacier-clad volcanoes and has been controversial at Katla. For this new seismic cluster on the south flank we regard volcano-related processes as more likely than glacial ones for the following reasons: 1) the seismic activity started during an unrest episode involving sudden melting of the glacier and a j\"okulhlaup; 2) the glacier stream is small and stagnant; 3) the seismicity remains regular and stable for years; 4) there is no apparent correlation with short-term weather changes, such as rain storms. We suggest that a small, shallow hydrothermal system was activated on Katla's south flank in 2011, either by a minor magmatic injection or by changes of permeability in a local crack system.