Between 1993 and 1998, the Hengill volcanic region in southwest Iceland was subject to a significant volcano-tectonic event which caused a total of 8 cm uplift and induced >90 thousand earthquakes, ...most of magnitude smaller than 2 ML and a small intrusion has been proposed as its cause. In this paper, we study fault activation around the uplift center to analyze whether the Coulomb stress change induced by such an intrusion would cause the observed crustal response. The locations of a subset of 12,000 micro-earthquakes are significantly improved by a double-difference relocation algorithm and the use of cross-correlation of waveforms to refine phase arrival times. The vast majority of seismic activity is scattered, but twenty-five (25) rectangular prisms, consisting of 11 to 363 earthquakes and 200 to 3000 m long, were identified and are interpreted as faults. Most of these faults are sub-vertical (within 5° from vertical) and only two faults dip at angles of 62° and 65°, respectively. Focal mechanism analysis of individual events was applied to estimate slip direction on these faults. Simulation of the Coulomb stress changes reveals that the stress changes would have been able to induce the observed fault movement on 22 of 25 faults. Discrepancies on the remaining faults are likely caused by uncertainties in the rake directions. Our results emphasize that the seismicity caused by a volume change in the crust is controlled by the interplay between the local stress field and the induced stresses.
•Relative relocations of earhtquakes in Hengill, SW-Iceland•Joint interpretation of locations and focal mechanisms to map faults•Twenty-three of twenty-five faults are near-vertical strike-slip faults•Modeled Coulomb stress changes account for movement on faults
In 1991, a new seismic monitoring network named SIL was started in Iceland with a digital seismic system and automatic operation. The system is equipped with software that reports the automatic ...location and magnitude of earthquakes, usually within 1–2 min of their occurrence. Normally, automatic locations are manually checked and re-estimated with corrected phase picks, but locations are subject to random errors and systematic biases. In this article, we consider the quality of the catalogue and produce a revised catalogue for South Iceland, the area with the highest seismic risk in Iceland. We explore the effects of filtering events using some common recommendations based on network geometry and station spacing and, as an alternative, filtering based on a multivariate analysis that identifies outliers in the hypocentre error distribution. We identify and remove quarry blasts, and we re-estimate the magnitude of many events. This revised catalogue which we consider to be filtered, cleaned, and corrected should be valuable for building future seismicity models and for assessing seismic hazard and risk. We present a comparative seismicity analysis using the original and revised catalogues: we report characteristics of South Iceland seismicity in terms of
b
value
and magnitude of completeness. Our work demonstrates the importance of carefully checking an earthquake catalogue before proceeding with seismicity analysis.
Glacier runoff and melt from volcanic and geothermal activity accumulates in glacier dammed lakes in glaciated areas around the world. These lakes eventually drain, creating hazardous subglacial ...floods that are usually only confirmed after they exit the glacier and reach local river systems, which can be many tens of kilometres from the flood source. Once in the river systems, they travel rapidly to populated areas. Such delayed detection represents a potentially lethal shortcoming in early-warning. Here we demonstrate how to advance early-warning potential through the analysis of four such floods in a glaciated region of Iceland. By comparing exceptional multidisciplinary hydrological, GPS and seismic ground vibration (tremor) data, we show that array analysis of seismic tremor can be used for early location and tracking of the subglacial flood front. Furthermore the timing and size of the impending flood can be estimated, prior to it entering the river system. Advanced warnings of between 20 to 34 hours are achieved for large (peak discharge of more than 3000 m
/s, accumulation time of ~ 5.25 years) to small floods (peak discharges from 210 to 380 m
/s, accumulation times of ~ 1.3 years) respectively.
SUMMARY
The reliable forecasting of seismic sequences following a main shock has practical implications because effective post-event response is crucial in earthquake-stricken regions, aftershocks ...can progressively cause increased damage and compound economic losses. In the South Iceland Seismic Zone (SISZ), one of two large transform zones in Iceland where earthquake hazard is the highest, an intense seismic sequence took place during 17–24 June 2000, starting with a ${M}_{\rm{w}}$ 6.4 main shock on 17 June 2000, followed by another ${M}_{\rm{w}}$ 6.5 main shock four days later and on a different fault. Both earthquakes caused considerable damage and incurred heavy economic losses. They were immediately followed by intense aftershock activity on the causative faults and triggered earthquakes as far as 80 km away along the transform zone. To investigate the feasibility of forecasting the progression of such complex sequences, we calibrated a spatio-temporal epidemic-type aftershock sequence (ETAS) clustering model to the June 2000 seismic sequence in the framework of Bayesian statistics. Short-term seismicity forecasts were carried out for various forecasting intervals and compared with the observations, the first generated a few hours after the first main shock and followed by daily forecasts. The reliability of the early forecasts was seen to depend on the initial model parameters. By using an adaptive parameter inference approach where the posteriors from each preceding forecasting interval served as informative priors for the next, the fast convergence of the parametric values was ensured. As a result, the 16–84 percentile range of the forecasted number of events captured the actual number of observed events in all daily forecasts, and the model exhibited a strong spatial forecasting ability, even only a few hours after the main shock, and over all subsequent daily forecasts. We present the spatio-temporal ETAS parameters for the June 2000 sequence as ideal candidates of prior estimates for future operational earthquake forecasting of other Icelandic aftershock sequences. Past seismic sequences need to be analysed retrospectively to confirm the stability of the parameters of this study, effectively enable the application of the Bayesian ETAS model as an operational earthquake forecasting system for aftershocks in Iceland.
Recent decades have highlighted the increasing need to connect and strengthen the volcanology community at European level. Indeed, research in the volcanology field is highly qualified in Europe and ...the volcano monitoring infrastructures have achieved valuable know-how, becoming the state-of-the-art in the world. However, the lack of common good practices in sciences and technologies, missing standards, as well as a significant fragmentation of the community requires coordination to move forward and guarantee a trans-national harmonisation. The European Plate Observing System (EPOS) represented the first opportunity to initiate this process of coordination by encouraging the creation of a European volcanological scientific infrastructure for data and service sharing. During the preparation and the design of EPOS, the volcanology community identified the objectives and the needs of the community building, the services to be provided and the work plan to implement the infrastructure. To achieve this aim, the contribution from three European projects FUTUREVOLC, MED-SUV and EUROVOLC was essential. This paper presents the main steps performed during the last years for building the community and implementing the infrastructure. This paper also describes the strategic choices and actions taken to realise the infrastructure such as the establishment of the Volcano Observation Thematic Core Service (TCS), whose structure and activity are described.
Strong earthquakes occurring frequently in south Iceland are always characterized by intense aftershock activity. The aftershock sequences cause progressive damages to structures and considerable ...psychological strain on inhabitants as the threat of upcoming large earthquakes looms. In this study, we apply a Bayesian epidemiological spatiotemporal aftershock sequence (ETAS) clustering model to forecast the aftershock sequence that commenced with the damaging 29 May 2008 Mw6.3 Ölfus earthquake. The forecast is carried out manually but in a framework of a provisional operational earthquake forecasting (OEF) system and based on an informative set of ETAS parameters derived from the strong June 2000 sequence in the region. We show that the short-term forecasts are unreliable at the early stages of the ongoing seismic sequence due to incomplete aftershocks collected to feed the Bayesian updating procedure of the ETAS model parameters, mainly during the first hours after the Ölfus mainshock. In contrast, issuing 24-h long forecasts, starting 6–8 h after the mainshock gives the best forecasting results. We allow the daily forecasts to adapt to the ongoing 2008 sequence using the adaptive prior estimation method, i.e., by assigning an informative prior distribution to ETAS parameters for the subsequent forecasting interval based on the preceding ETAS parameter posterior distributions. The parametric posteriors do not show narrow distributions, thus indicating inconsistent forecasting performance, both in stark contrast to the June 2000 sequence. Regardless, our results provide a quantitative estimate of the reliability of the forecasts which is a prerequisite for evaluating their usefulness to civil protection.
•Calibrating a Bayesian ETAS model to forecast an intense seismic sequence in Iceland•Short-term aftershock forecasting of the damaging 29 May 2008 M6.3 earthquake•Comparison of retrospective forecasting of June 2000 and May 2008 seismic sequences•Daily forecasts using informative and adaptive ETAS priors over an ongoing sequence
We recorded high-frequency (>10Hz) harmonic tremor with spectral gliding at Hekla Volcano in Iceland. Particle motion plots indicated a shallow tremor source. We observed up to two overtones beneath ...our Nyquist frequency of 50Hz and could resolve a source of closely spaced pulses of very short duration (0.03-0.1s) on zoomed seismograms. Volcanic tremor with fundamental frequencies above 5Hz, frequency gliding and/or repetitive sources similar to our observations were observed on different volcanoes around the world. However, this frequency content, duration and occurrence of volcano-related tremor was not observed in the last 35years of seismic observations at Hekla. Detailed analysis reveals that this tremor was related to helicopters passing the volcano. This study relates the GPS track of a helicopter with seismic recordings of the helicopter at various distances. We show the effect the distance, number of rotor blades and velocity of the helicopter has on the observed up and down glidings at up to 40km distance. We highlight similarities and differences between volcano-related and helicopter tremor in order to help avoid misinterpretations.
•Helicopter-generated tremor consists of regularly repeating pulses.•Fundamental frequencies above 10Hz and overtones at integer multiples are observed.•Doppler Effect causes the frequency gliding and is used to derive the helicopter speed.•Effects of distance, speed and wind direction is revealed using a GPS track.•Signals could be misinterpreted as volcanic tremor on poorly monitored volcanoes.
Crust at many divergent plate boundaries forms primarily by the injection of vertical sheet-like dykes, some tens of kilometres long. Previous models of rifting events indicate either lateral dyke ...growth away from a feeding source, with propagation rates decreasing as the dyke lengthens, or magma flowing vertically into dykes from an underlying source, with the role of topography on the evolution of lateral dykes not clear. Here we show how a recent segmented dyke intrusion in the Bárðarbunga volcanic system grew laterally for more than 45 kilometres at a variable rate, with topography influencing the direction of propagation. Barriers at the ends of each segment were overcome by the build-up of pressure in the dyke end; then a new segment formed and dyke lengthening temporarily peaked. The dyke evolution, which occurred primarily over 14 days, was revealed by propagating seismicity, ground deformation mapped by Global Positioning System (GPS), interferometric analysis of satellite radar images (InSAR), and graben formation. The strike of the dyke segments varies from an initially radial direction away from the Bárðarbunga caldera, towards alignment with that expected from regional stress at the distal end. A model minimizing the combined strain and gravitational potential energy explains the propagation path. Dyke opening and seismicity focused at the most distal segment at any given time, and were simultaneous with magma source deflation and slow collapse at the Bárðarbunga caldera, accompanied by a series of magnitude M > 5 earthquakes. Dyke growth was slowed down by an effusive fissure eruption near the end of the dyke. Lateral dyke growth with segment barrier breaking by pressure build-up in the dyke distal end explains how focused upwelling of magma under central volcanoes is effectively redistributed over long distances to create new upper crust at divergent plate boundaries.