Array observations of an oscillating seismic sequence in the Reykjanes Peninsula, SW-Iceland, in December 2021&#160;
Authors: Blanck Hanna , Halldórsson Benedikt , Vogfjord Kristín 
Published: March 2022

&lt;div&gt;
&lt;p&gt;&lt;span&gt;In the evening hours of 21&amp;#160;December 2021, a seismic sequence started in south-central&amp;#160;Reykjanes peninsula in SW-Iceland. Both the science community and the civil protection agency were alarmed due to the proximity of this sequence&amp;#160;to the site of the 2021 Fagradalsfjall eruption (March &amp;#8211; September&amp;#160;2021), especially as&amp;#160;this&amp;#160;was the&amp;#160;most prominent&amp;#160;sequence since the end of the eruption and it&amp;#160;showed similar characteristic as the seismic activity that had been observed in the 3 weeks leading up to it.&amp;#160;In addition, the December&amp;#160;earthquake sequence was&amp;#160;located along a NE-SW striking alignment which, together with GPS and InSAR measurements, has been interpreted as a dike intrusion, which also was the origin of the March eruption.&amp;#160;We analyse the seismic activity using&amp;#160;a&amp;#160;small-aperture (D=1.7 km, d=0.5 km) urban seismic array, consisting of 5 Raspberry Shake 4D sensors (1 vertical geophone and 3&amp;#160;MEM&amp;#160;accelerometric components) located&amp;#160;in the nearby&amp;#160;municipality of Grindav&amp;#237;k&amp;#160;about 10 km WSW from the former eruption site.&amp;#160;During the first days of the seismic activity magnitudes reached up to M&lt;/span&gt;&lt;span&gt;L&lt;/span&gt;&lt;span&gt; 4.8 but on 30&amp;#160;December the activity subsided and then ceased,&amp;#160;with only few events reaching more than&amp;#160;M&lt;/span&gt;&lt;span&gt;L &lt;/span&gt;&lt;span&gt;2, which coincides with&amp;#160;the magnitude of completeness of the seismic array.&amp;#160;&lt;/span&gt;&lt;span&gt;&amp;#160;&lt;/span&gt;&lt;/p&gt;
&lt;/div&gt;&lt;div&gt;
&lt;p&gt;&lt;span&gt;We present the first insights into the spatiotemporal characteristics of the sequence provided by array processing of the most intense period of the sequence. To process the array data, we used the SeisComP module AUTOLAMBDA with&amp;#160;both the FK and PMCC (Progressive Multi-Channel Correlation) method to obtain back azimuth and slowness pairs of incoming waves.&amp;#160;During its first hours, the sequence&amp;#160;showed a systematic behaviour in the back azimuthal distribution of the incoming waves. Namely, over a repeated interval of a couple of hours the back azimuthal estimates increase steadily at a rate of 5 to 12&amp;#176;/h after which the source of the activity appears to drop back to the initial azimuthal values, and the cycle repeats. Over the following days, these bursts of oscillating activity become less frequent with relatively calm phases between. These periods of oscillating behaviour show that the seismic activity was systematically&amp;#160;migrating southwest to/from northeast and most likely is the signature of a pulsating&amp;#160;magma pressure front&amp;#160;in the dike itself. This behaviour is similar to some phases during the previous eruption&amp;#160;when lava was actively&amp;#160;erupting with hours of quiescence in between. These results show that the monitoring of automatic back azimuth and slowness estimates are a&amp;#160;useful tool in revealing&amp;#160;small-scale systematic behaviour of seismic sequences in the area&amp;#160;in real-time.&lt;/span&gt;&lt;span&gt;&amp;#160;&lt;/span&gt;&lt;/p&gt;
&lt;/div&gt;