Abstract
Volcanic eruptions at mid-ocean ridges are rarely witnessed due to their inaccessibility, and are therefore poorly understood. Shallow waters in the Red Sea allow the study of ocean ridge related volcanism observed close to sea level. On the 18th December 2011, Yemeni fishermen witnessed a volcanic eruption in the Southern Red Sea that led to the formation of Sholan Island. Previous research efforts to constrain the dynamics of the intrusion and subsequent eruption relied primarily on interferometric synthetic aperture radar (InSAR) methods, data for which were relatively sparse. Our study integrates InSAR analysis with seismic data from Eritrea, Yemen, and Saudi Arabia to provide additional insights into the transport of magma in the crust that fed the eruption. Twenty-three earthquakes of magnitude 2.1–3.9 were located using the Oct-tree sampling algorithm. The earthquakes propagated southeastward from near Sholan Island, mainly between December 12th and December 13th. The seismicity is interpreted as being induced by emplacement of a ∼12 km-long dike. Earthquake focal mechanisms are primarily normal faulting and suggest the seismicity was caused through a combination of dike propagation and inflation. We combine these observations with new deformation modeling to constrain the location and orientation of the dike. The best-fit dike orientation that satisfies both geodetic and seismic data is NNW-SSE, parallel to the overall strike of the Red Sea. Further, the timing of the seismicity suggests the volcanic activity began as a submarine eruption on the 13th December, which became a subaerial eruption on the 18th December when the island emerged from the beneath the sea. The new intrusion and eruption along the ridge suggests seafloor spreading is active in this region.
Introduction
Magma motion through the crust can induce earthquakes, the analysis of which is a useful tool to understand when, where, and how magma is transported. The understanding of magma intrusion and volcanism is key for understanding dynamic processes that occur at divergent boundaries (e.g., ; ; ; ). However, the understanding of the extensional processes at submarine ocean ridges in space and time is limited by the difficulty associated with observing seismicity and ground deformation in the deep oceans (e.g., ; ). Over the past 20 years, the Southern Red Sea has experienced six seismic swarms interpreted to be related to magma intrusions (). The Southern Red Sea has a large proportion of shallow waters, which provide a rare opportunity to study volcanic eruptions that occur above sea level (Figure 1). Three sub-aerial volcanic eruptions have occurred in the last decade: on Jebel at Tair island (2007) and two eruptions resulting in the formation of Sholan Island (2011) and Jadid Island (2013) (). All the islands lie on the rift axis of the Southern Red Sea. The eruptions of Jebel at Tair and Jadid lasted 8 and 2 months, respectively, and followed short duration (days) low magnitude (M > 4) seismic swarms thought to be caused by the intrusion of magma through the crust to the Earth’s surface (). The eruption of Jebel at Tair was subaerial, while the first half of the Jadid eruption was submarine ().
FIGURE 1
The eruption of Sholan Island occurred between Haycock and Rugged islands in the Zubair Archipelago and was first witnessed by Yemeni fishermen on the 18th December 2011 (Figure 2). The eruption started as submarine activity before erupting as a surtseyan eruption, and lasted until the 12th January 2012 (
FIGURE 2

The Zubair archipelago in the Southern Red Sea. Sholan and Jadid Island, erupted in 2011 and 2013, respectively, are shown in red. Proposed dike locations from
Preceding the Sholan eruption, three distinct seismic swarms occurred in April, June, and August 2011 and the Yemeni seismological network detected two earthquakes of magnitude 3.7 and 3.9 on the 13th December 2011, taken from the International Seismic Catalogue (ISC). In this study we investigate the temporal and spatial variability in seismic activity prior to the eruption of Sholan Island. Analysis of seismic data from local seismic networks in Eritrea, Yemen, and Saudi Arabia allows us to locate earthquakes surrounding the island which were not recorded by the ISC. The increased number of earthquakes that we detect, coupled with improved relative earthquake locations, allows us to place constraints on the location, timescales and direction of intrusion that fed the eruption. In addition, we use the seismic results to guide new deformation models created from InSAR data. Our data of land-based stations positioned relatively close to an ocean ridge place rare constraints on a seafloor-spreading episode and provide additional insights into the rifting mechanics of the Southern Red Sea.
Tectonic Setting
The Red Sea formed due to the divergence of the Nubian and Arabian plates (
South of 17°N, GPS measurements show the locus of extension splits into two branches: the continuation of the Red Sea ridge in the east, and the Danakil depression to the west (
The Zubair archipelago (at 15–15.2°N) is a 25 km by 10 km shallow platform [<100 m below sea level (bsl)] orientated parallel to the Red Sea ridge; it consists of ∼10 volcanic islands situated along the central axis of the Red Sea (
Data and Methods
We used seven broadband seismometers in Eritrea, Yemen, and Saudi Arabia (Figure 1). A Butterworth band pass filter was applied to the data to only allow frequencies between 1 and 10 Hz for the stations in Eritrea, and 1–8 Hz for the Yemen and Saudi Arabia stations due to a higher level of noise. Earthquakes were manually picked for both P and S waves, and events with a minimum of four arrival times at three stations were located with NonLinLoc, using the Oct-Tree Sampling Algorithm (
Magnitudes are calculated by measuring the maximum peak-to-peak amplitude and using the local magnitude scale of
Results
Using arrival times from a minimum of 3 and up to 7 seismic stations, we locate 28 individual earthquakes in NonLinLoc in the vicinity of the islands from the 5th December 2011 to the 12th January 2012. Latitude and longitude errors were calculated in NonLinLoc; these are statistical errors based on the probability density function of each earthquake’s location (Supplementary Figure 3). Five additional earthquakes were removed due to having error bars larger than 5 km, leaving 23 earthquakes that were used for final analysis that define a broadly NW-SE trend (Figure 3). Further information on the earthquakes location, errors and statistics can be found in Supplementary Table 1. The main seismic events are focused within the Zubair archipelago and occurred during the 12th and 13th December with just two events after this time period (Figure 4). The first earthquakes, on the 12th December occur to the north of Sholan island, beneath Haycock island, and then migrate southeast to Saba island with time (Figure 3). The computed errors in horizontal earthquake location are within ±5 km. The PDF scatter clouds produced by NonLinLoc can be seen in Figure 3.
FIGURE 3

The Sholan seismic sequence with earthquakes colored by time. Focal mechanisms show three normal faulting and two strike-slip faulting events. The black dashed line denotes the inferred dike location taken from earthquake hypocenters only; the red line denotes the final dike location taken from combining ground deformation and earthquake hypocenters.
FIGURE 4

Cumulative seismic moment release and number of events over December 2011 for the Sholan seismic sequence. Inset: histogram displaying the range of magnitudes for events in December 2011. Red line denotes the start of the subaerial eruption.
The calculated local magnitudes (ML) vary from 2.1 to 3.9, with the majority of earthquakes ranging from 2.5 to 3.5 (Figure 4). No temporal trend in ML with time can be seen over the 12th–13th December. However, a cluster of low magnitude events is observed at the end of the 12th December. Larger magnitudes are observed north of the island with ML 2–3 earthquakes further south (Figure 3). Cumulative seismic moment calculated over December shows the majority of the energy release occurring on the 12th–13th December, in agreement with the number of events. A gradual increase in energy release to 19 × 1014 Nm is observed over the 12th and 13th December, with little additional energy from the two events on the 16th and 27th (Figure 4). Due to quality of data, only five focal mechanisms were adequately constrained (Table 1 and Figure 3). The solutions show a combination of normal fault and strike-slip fault plane solutions (Figure 3).
Table 1
| ID | DATE | TIME | LAT. | LONG. | MAG. | STRIKE | DIP | RAKE |
|---|---|---|---|---|---|---|---|---|
| 1 | 12/12/11 | 16:59:27 | 15.182 | 42.092 | 3.27 | 127.4 | 44.9 | -35.3 |
| 2 | 12/12/11 | 20:45:20 | 15.164 | 42.126 | 2.99 | 117.0 | 54.9 | -36.3 |
| 3 | 13/12/11 | 01:22:12 | 15.188 | 42.109 | 3.92 | 323.5 | 50.1 | -86.1 |
| 4 | 13/12/11 | 05:01:49 | 15.190 | 42.095 | 3.92 | 322.1 | 50.1 | -86.1 |
| 5 | 13/12/11 | 13:02:31 | 15.113 | 42.164 | 3.47 | 322.1 | 50.1 | -86.7 |
Focal plane solutions of five events on 12th–13th December 2011.
Discussion
The locations of 23 earthquakes in the Zubair Archipelago shows a clustering of events that follow a NW-SE trend. The variation of event locations with time shows a progression toward the south originating from near Sholan island, with larger magnitude events occurring at either end of the sequence. Earthquakes are commonly induced at the leading edge of a propagating dike due to concentration of extensional stresses (e.g.,
The focal plane solutions show two focal mechanisms (1, 2) which occur on the 12th December in the early stages of dike propagation (Figure 3). Both these events display a large component of strike-slip deformation, with p-axes oriented 70° from the regional maximum horizontal compressive stress (σ1R). These two earthquakes are best explained by a model proposed by
In order to reconcile the observed ground deformation pattern, eruption and ground fracturing locations, and the seismicity we used rectangular dislocations in an elastic halfspace to model the dike opening, and calculated the expected surface deformation from several different possible dike locations. For the models we use a dike depth of 0–5 km, vertical dip, and thickness of 0.5 m. We show the modeled ground deformation patterns from
FIGURE 5

Interferograms produced by modeling the ground deformation pattern using different dike locations. (A) Original dike location produced by
In order to understand the migration of the earthquakes with time, a projection was taken along the length of the inferred dike and compared with the hypocenters of the earthquakes (Figure 6). Here we can clearly see that the southward dike propagation occurs between 17:00 on the 12th December and 02:00 on the 13th December, and by varying isotachs (lines of constant speed) we estimate that the dike propagated at a rate of 0.32 ms-1. This compares with a dike propagation event seen at Bár
arbunga in August 2014 which had dike propagation rates of up to 0.83 m/s (
FIGURE 6

The Sholan seismic sequence projected along the length of the proposed dike. The dike propagation phase is clearly seen with activity at either end of the dike continuing once propagation has ceased. The red dashed line shows the 0.32 ms-1 isotach indicating the speed of dike propagation.
The cessation of the dike to the south is likely due to the dike reaching a stress barrier such as caused by a topographic load (
Little seismicity is recorded after the 13th December, with only one further earthquake recorded prior to the first observation of the subaerial eruption on the 18th December. We suggest that the eruption began as a submarine eruption on the 13th December once the dike had ceased propagating to the south. We interpret this sequence of events using insights from subaerial fissure eruptions in Afar (
There is uncertainty in the nature of the Red Sea ridge between 15.5 and 14.5°N. A lack of magnetic anomalies suggests an absence in seafloor spreading in this region, however, bathymetry data suggest the ridge continues to 14.5°N (
Conclusion
We analyse a swarm of earthquakes in the Red Sea at 15–15.2°N and 42–42.3°E, associated with a subaerial volcanic eruption starting on the 18th December 2011 that resulted in the formation of Sholan Island. Intense earthquake activity on the 12th December forms a NW-SE trending cluster with the time evolution of seismicity showing a southward migration over a day period. The rapid migration of earthquakes, coupled with earthquakes focused at the lateral tips of the cluster after migration stops strongly suggests the earthquakes are caused by an intrusion that fed the eruption, broadly consistent with surface deformation patterns best modeled using the presence of a dike. The seismicity migration suggests that the dike propagated at speeds of 0.32 ms-1 and was approximately 12 km in length. Analysis of focal mechanisms suggests that the seismicity in the early stages of dike propagation was caused by stresses around the expanding dike. Once, propagation had halted seismicity is dominated by normal faulting suggesting dike-induced extension. We interpret dike position and a NNW-SEE orientation based on tested forward models of ground deformation pattern based on both geodetic and the new seismicity data. Observations from other dike intrusions in Afar, coupled with high solubility of SO2 in water, suggest that a submarine eruption may have begun on the 13th December once the dike had ceased propagating. The intrusion of the dike parallel to the orientation of the Red Sea suggests seafloor spreading is still active in this region, with the increased seismicity and magmatic activity likely related to a spreading event.
Statements
Author contributions
JE conducted seismology analysis, and led the interpretation and writing of the paper. FI-K contributed to seismology analysis, interpretation, and writing the paper. DK contributed to project design, interpretation, and writing the paper. JR contributed to the project design, interpretation, and writing the paper. SJ conducted the InSAR analysis and related modeling, and contributed to interpretation and writing the paper.
Funding
FI-K was funded through NERC studentship NE/L002531/1 and a grant to GSNOCS from Roy Franklin O.B.E and the ECLIPSE Program funded by the New Zealand Ministry of Business, Innovation and Employment. DK was supported by NERC grant NE/L013932, and grant number OSR-2015-CRG4-2643 from King Abdullah University of Science and Technology. JE was supported by NERC studentship NE/L002582/1.
Acknowledgments
We acknowledge the support of Martin Mai (KAUST) in facilitating acquisition of data from Farasan Island from the Saudi Geological Survey, who are also thanked. The remaining data was downloaded from the IRIS-DMC. We thank the reviewers EF and LDS, editor BT, and chief editor Valerio Acocella for constructive reviews.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2018.00141/full#supplementary-material
References
1
AhmedA.DoubreC.LeroyS.KassimM.KeirD.AbayazidA.et al (2016). Seafloor spreading event in western Gulf of Aden during the November 2010 - March 2011 period captured by regional seismic networks: evidence for diking events and interactions with a nascent transform zone.Geophys. J. Int.2051244–1266. 10.1093/gji/ggw068
2
AlmalkiK. A.BettsP. G.AilleresL. (2016). Incipient seafloor spreading segments: insights from the Red Sea.Geophys. Res. Lett.432709–2715. 10.1002/2016GL068069
3
BarnieT. D.KeirD.HamlingI.HofmannB.BelachewM.CarnS.et al (2016). A multidisciplinary study of the final episode of the Manda Hararo dike sequence, Ethiopia, and implications for trends in volcanism during the rifting cycle.Geol. Soc. Special Publ.420149–163. 10.1144/SP420.6
4
BarrancosJ.PadillaG.PadronE.HernándezP.CalvoD.MarquezA.et al (2012). Estimated CO2, SO2 and H2S emission to the atmosphere from the 2011 El Hierro submarine eruption (Canary Islands) on the basis of helicopter gas surveys.Geophys. Res. Abstr.14:12345.
5
BelachewM.EbingerC.CoteD.KeirD.RowlandJ.HammondJ.et al (2011). Comparison of dike intrusions in an incipient seafloor- spreading segment in Afar, Ethiopia: seismicity perspectives.J. Geophys. Res.116:B06405. 10.1029/2010JB007908
6
BerckhemerH.BaierB.BartlesenH.BehleA.BurkhardtH.GebrandeH.et al (1975). “Deep seismic soundings in the Afar region and on the highland of Ethiopia,” inAfar Depression of EthiopiaedsPilgerA.öslerA. R. (Stuttgart: Schweizerbart) 89–107.
7
BonattiE. (1985). Punctiform initiation of seafloor spreading in the Red Sea during transition from a continental to an oceanic rift.Nature31633–37. 10.1038/316033a0
8
BosworthW.HuchonP.McClayK. (2005). The Red Sea and Gulf of Aden basins.J. Afr. Earth Sci.43334–378. 10.1016/j.jafrearsci.2005.07.020
9
BrandsdottirB.EinarssonP. (1979). Seismic activity associated with the september 1977 deflation of the krafla central volcano in northeastern ice- land.J. Volcanol. Geothermal Res.6197–212. 10.1016/0377-0273(79)90001-5
10
ButterfieldD. A.NakamuraK. I.TakanoB.LilleyM. D.LuptonJ. E.ResingJ. A.et al (2011). High SO2 flux, sulfur accumulation, and gas fractionation at an erupting submarine volcano.Geology39803–806. 10.1130/G31901.1
11
CareyR.SouleS. A.MangaM.WhiteJ. D.McPhieJ.WysoczanskiR.et al (2018). The largest deep-ocean silicic volcanic eruption of the past century.Sci. Adv.4:e1701121. 10.1126/sciadv.1701121
12
EbingerC. J.van WijkJ.KeirD. (2013). The time scales of continental rifting: implications for global processes.Geol. Soc. Am. Special Papers500371–396. 10.1130/2013.2500(11)
13
EgloffF.RihmR.MakrisJ.IzzeldinY.BobsienM.MeierK.et al (1991). Contrasting structural styles of the eastern and western margins of the southern Red Sea: the 1988 SONNE experiment.Tectonophysics198329–353. 10.1016/0040-1951(91)90159-P
14
EinarssonP. (1991). “The Krafla rifting episode 1975–1989,” inNattura Myvatns (The Nature of Lake Myvatn)edsGardarssonA.EinarssonA. (Reykjavic: Icelandic Nature Science Society) 97–139.
15
FieldL.BarnieT.BlundyJ.BrookerR. A.KeirD.LewiE.et al (2012). Integrated field, satellite and petrological observations of the November 2010 eruption of Erta Ale.Bull. Volc.742251–2271. 10.1007/s00445-012-0660-7
16
GoitomB.OppenheimerC.HammondJ. O. S.GrandinR.BarnieT.DonovanA.et al (2015). First recorded eruption of Nabro volcano, Eritrea, 2011.Bull. Volc.77:85. 10.1007/s00445-015-0966-3
17
HamlynJ. E.KeirD.WrightT. J.NeubergJ. W.GoitomB.HammondJ. O. S.et al (2014). Seismicity and subsidence following the 2011 Nabro eruption, Eritrea: insights into the plumbing system of an off-rift volcano.J. Geophys. Res.1198267–8282. 10.1002/2014JB011395
18
HanksT. C.KanamoriH. (1979). A moment magnitude scale.J. Geophys. Res.842348–2350. 10.1029/JB084iB05p02348
19
Illsley-KempF.KeirD.BullJ. M.AyeleA.HammondJ. O.KendallJ.-M.et al (2017). Local earthquake magnitude scale and b-value for the Danakil region of northern Afar.Bull. Seis. Soc. Am.107521–531. 10.1785/0120150253
20
JónssonS.XuW. (2015). Volcanic Eruptions in the Southern Red Sea During 2007-2013. The Red Sea.Berlin: Springer175–186. 10.1007/978-3-662-45201-1_10
21
KanamoriH. (1977). The energy release in great earthquakes.J. Geophys. Res.822981–2987. 10.1029/JB082i020p02981
22
KeirD.HamlingI. J.AyeleA.CalaisE.EbingerC.WrightT. J.et al (2009). Evidence for focused magmatic accretion at segment centers from lateral dike injections captured beneath the Red Sea rift in Afar.Geology3759–62. 10.1130/G25147A.1
23
Le PichonX. T.GaulierJ.-M. (1988). The rotation of Arabia and the levant fault system.Tectonophysics153271–294. 10.1016/0040-1951(88)90020-0
24
LomaxA.VirieuxJ.VolantP.Berge-ThierryC. (2000). “Probabilistic earth- quake location in 3D and layered models,” inAdvances in Seismic Event Location, edsThurberC. H.RabinowitzN. (Berlin: Springer), 101–134. 10.1007/978-94-015-9536-0_5
25
McCluskyS.ReilingerR.OgubazghiG.AmlesonA.HealebB.VernantP.et al (2010). Kinematics of the southern Red Sea-Afar Triple Junction and implications for plate dynamics.Geophys. Res. Lett.37:L05301. 10.1029/2009GL041127
26
PagliC.MazzariniF.KeirD.RivaltaE.RooneyT. O. (2015). Introduction: anatomy of rifting: tectonics and magmatism in continental rifts, oceanic spreading centers, and transforms.Geosphere111256–1261. 10.1130/GES01082.1
27
PagliC.WrightT. J.EbingerC. J.YunS.-H.CannJ. R.BarnieT.et al (2012). Shallow axial magma chamber at the slow-spreading Erta Ale Ridge.Nat. Geosci.5284–288. 10.1038/ngeo1414
28
QinR.BuckW. R. (2008). Why meter-wide dikes at oceanic spreading centres?Earth Planet. Sci. Lett.265466–474. 10.1016/j.epsl.2007.10.044
29
RivaltaE.TaisneB.BungerA.KatzR. (2015). A review of mechanical models of dike propagation: schools of thought, results and future directions.Tectonophysics6381–42. 10.1016/j.tecto.2014.10.003
30
RomanD. C. (2005). Numerical models of volcanotectonic earthquake triggering on non-ideally oriented faults.Geophys. Res. Lett.32:4. 10.1029/2004GL021549
31
RomanD. C.CashmanK. V. (2006). The origin of volcano-tectonic earthquake swarms.Geology34457–460. 10.1130/G22269.1
32
RubinA. M.GillardD. (1998). Dike-induced earthquakes: theoretical considerations.J. Geophys. Res.10310017–10030. 10.1029/97JB03514
33
RuchJ.WangT.XuW.HenschM.JónssonS. (2016). Oblique rift opening revealed by reoccurring magma injection in central Iceland.Nat. Commun.7:12352. 10.1038/ncomms12352
34
SchettinoA.MacchiavelliC.PierantoniP. P.ZanoniD.RasulN. (2016). Recent kinematics of the tectonic plates surrounding the Red Sea and Gulf of Aden.Geophys. J. Int.207457–480. 10.1093/gji/ggw280
35
ScordilisE. M. (2006). Empirical global relations converting Ms and mb to moment magnitude.J. Seismol.10225–236. 10.1007/s10950-006-9012-4
36
SigmundssonF.HooperA.HreinsdóttirS.VogfjördK. S.ÓfeigssonB. G.HeimissonE. R.et al (2015). Segmented lateral dyke growth in a rifting event at Bararbunga volcanic system.Icel. Nat.517191–195. 10.1038/nature14111
37
SnokeJ. A. (2003). FOCMEC: focal mechanism determinations.Int. Handb. Earthq. Eng. Seismol.851629–1630. 10.1016/S0074-6142(03)80291-7
38
SouleS. A.FornariD. J.PerfitM. R.RubinK. H. (2007). New insights into mid-ocean ridge volcanic processes from the 2005-2006 eruption of the East Pacific Rise, 9 46 N-9 56 N.Geology351079–1082. 10.1130/G23924A.1
39
UrbaniS.AcocellaV.RivaltaE.CorbiF. (2018). Propagation and arrest of dikes under topography: models applied to the 2014 Bardarbunga (Iceland) rifting event.Geophys. Res. Lett.446692–6701. 10.1002/2017GL073130
40
WilcockW. S.TolstoyM.WaldhauserF.GarciaC.TanY. J.BohnenstiehlD. R.et al (2016). Seismic constraints on caldera dynamics from the 2015 axial seamount eruption.Science3541395–1399. 10.1126/science.aah5563
41
WrightT.EbingerC.BiggsJ.AyeleA.YirguG.KeirD.et al (2006). Magma-maintained rift segmentation at continental rupture in the 2005 Afar dyking episode.Nature442291–294. 10.1038/nature04978
42
WrightT. J.SigmundssonF.PagliC.BelachewM.HamlingI. J.BrandsdóttirB.et al (2012). Geophysical constraints on the dynamics of spreading centres from rifting episodes on land.Nat. Geosci.5242–250. 10.1038/ngeo1428
43
XuW.JónssonS. (2014). The 2007-8 volcanic eruption on Jebel at Tair island (Red Sea), observed by satellite radar and optical images.Bull. Volc.76795–808. 10.1007/s00445-014-0795-9
44
XuW.RuchJ.JónssonS. (2015). Birth of two volcanic islands in the southern Red Sea.Nat. Commun.6:7104. 10.1038/ncomms8104
Summary
Keywords
Red Sea, mid-ocean ridge, dike, seismicity, InSAR, eruption
Citation
Eyles JHW, Illsley-Kemp F, Keir D, Ruch J and Jónsson S (2018) Seismicity Associated With the Formation of a New Island in the Southern Red Sea. Front. Earth Sci. 6:141. doi: 10.3389/feart.2018.00141
Received
04 June 2018
Accepted
05 September 2018
Published
24 September 2018
Volume
6 - 2018
Edited by
Benoit Taisne, Nanyang Technological University, Singapore
Reviewed by
Luca De Siena, University of Aberdeen, United Kingdom; Eisuke Fujita, National Research Institute for Earth Science and Disaster Prevention, Japan
Updates

Check for updates
Copyright
© 2018 Eyles, Illsley-Kemp, Keir, Ruch and Jónsson.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Derek Keir, d.keir@soton.ac.uk
This article was submitted to Volcanology, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.