# First recorded eruption of Nabro volcano, Eritrea, 2011

Berhe Goitom, Clive Oppenheimer, James Hammond, Raphael Grandin, Talfan Barnie, Amy Donovan, Ghebrebrhan Ogubazghi, Ermias Yohannes, Goitom Goitom, J- Michael Kendall, et al.

► To cite this version:

Berhe Goitom, Clive Oppenheimer, James Hammond, Raphael Grandin, Talfan Barnie, et al. First recorded eruption of Nabro volcano, Eritrea, 2011. Bulletin of Volcanology, 2015, 77 (10), 10.1007/s00445-015-0966-3. hal-02114238

HAL Id: hal-02114238

https://uca.hal.science/hal-02114238

Submitted on 15 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# First recorded eruption of Nabro volcano, Eritrea, 2011

Berhe Goitom<fnref n="1" />,<fnref n="2" /> · Clive Oppenheimer<fnref n="3" /> · James O. S. Hammond<fnref n="4" /> · Raphaël Grandin<fnref n="5" /> ·

Talfan Barnie<fnref n="6" /> · Amy Donovan<fnref n="3" /> · Ghebrebrhan Ogubazghi<fnref n="2" /> · Ermias Yohannes<fnref n="7" /> ·

Goitom Kibrom<fnref n="7" /> · J- Michael Kendall<fnref n="1" /> · Simon A. Carn<fnref n="8" /> · David Fee<fnref n="9" /> ·

Christine Sealing<fnref n="8" /> · Derek Keir<fnref n="10" /> · Atalay Ayele<fnref n="11" /> · Jon Blundy<fnref n="1" /> · Joanna Hamlyn<fnref n="12" /> ·

Tim Wright<fnref n="12" /> · Seife Berhe<fnref n="13" />

Received: 23 March 2015 / Accepted: 10 August 2015 / Published online: 7 September 2015
© The Author(s) 2015. This article is published with open access at Springerlink.com

**Abstract** We present a synthesis of diverse observations of the first recorded eruption of Nabro volcano, Eritrea, which began on 12 June 2011. While no monitoring of the volcano was in effect at the time, it has been possible to reconstruct the nature and evolution of the eruption through analysis of regional seismological and infrasound data and satellite remote sensing data, supplemented by petrological analysis of erupted products and brief field surveys. The event is notable for the comparative rarity of recorded historical eruptions in the region and of caldera systems in general, for the prodigious quantity of SO₂ emitted into the atmosphere and the significant human impacts that ensued notwithstanding the low population density of the Afar region. It is also relevant in understanding the broader magmatic and tectonic significance of the volcanic massif of which Nabro forms a part and which strikes obliquely to the principal rifting directions in the Red Sea and northern Afar. The whole-rock compositions of

the erupted lavas and tephra range from trachybasaltic to trachybasaltic andesite, and crystal-hosted melt inclusions contain up to 3,000 ppm of sulphur by weight. The eruption was preceded by significant seismicity, detected by regional networks of sensors and accompanied by sustained tremor. Substantial infrasound was recorded at distances of hundreds to thousands of kilometres from the vent, beginning at the onset of the eruption and continuing for weeks. Analysis of ground deformation suggests the eruption was fed by a shallow, NW-SE-trending dike, which is consistent with field and satellite observations of vent distributions. Despite lack of prior planning and preparedness for volcanic events in the country, rapid coordination of the emergency response mitigated the human costs of the eruption.

**Keywords** Nabro · InSAR · Seismicity · Afar · Danakil · Volcano monitoring · Satellite remote sensing

Editorial responsibility: G. Giordano

<figure><img src="image_2.png" /></figure>

<footnote>Berhe Goitom
bg12363@bristol.ac.uk</footnote>

<footnote n="1">¹ School of Earth Sciences, University of Bristol, Queens Road, Bristol BS8 1RJ, UK</footnote>

<footnote n="2">² Department of Earth Sciences, Eritrea Institute of Technology, PO Box 12676, Asmara, Eritrea</footnote>

<footnote n="3">³ Department of Geography, Downing Place, Cambridge CB2 3EN, UK</footnote>

<footnote n="4">⁴ Department of Earth Science and Engineering, Imperial College, London SW7 2AZ, UK</footnote>

<footnote n="5">⁵ Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Univ Paris Diderot CNRS, 75005 Paris, France</footnote>

<footnote n="6">⁶ Laboratoire Magmas et Volcans, Université Blaise Pascal, Clermont Ferrand, France</footnote>

<footnote n="7">⁷ Department of Mines, Eritrea Geological Surveys, PO Box 272, Asmara, Eritrea</footnote>

<footnote n="8">⁸ Department of Geological and Mining Engineering and Sciences, Michigan Technological University, 1400 Townsend Dr, Houghton, MI 49931, USA</footnote>

<footnote n="9">⁹ Wilson Infrasound Observatories, Alaska Volcano Observatory, Geophysical Institute, University of Alaska Fairbanks, Fairbanks, AK, USA</footnote>

<footnote n="10">¹⁰ National Oceanography Centre Southampton, University of Southampton, Southampton SO14 3ZH, UK</footnote>

<footnote n="11">¹¹ Institute of Geophysics, Space Science and Astronomy, Addis Ababa University, Addis Ababa, Ethiopia</footnote>

<footnote n="12">¹² COMET, School of Earth and Environment, University of Leeds, Leeds, UK</footnote>

<footnote n="13">¹³ Global Resources Development Consultants, Asmara, Eritrea</footnote>
## Introduction

On 12 June 2011, Nabro volcano was the site of an eruption that had significant societal and environmental consequences. It is the first eruption of Nabro on record. Previously, the only securely recorded historical eruption in Eritrea was that of Dubbi volcano in 1861, located just over 30 km from the Red Sea coast (Wiart and Oppenheimer 2000; Wiart et al. 2000). That event reportedly claimed 175 lives despite the sparse population of the Southern Red Sea region of Eritrea, likely as a result of inundation of settled areas by pyroclastic density currents (Wiart and Oppenheimer 2000). Ash clouds and tephra fall associated with the 1861 eruption also disrupted navigation in the Red Sea.

Nabro volcano, some 25 km southwest from Dubbi, is part of a much larger massif referred to as the Bidu Volcanic Complex (Wiart and Oppenheimer 2005; Fig. 1), which is comprised of two calderas (Nabro and, just across the international border in Ethiopia, Mallahle). Nabro's 2011 eruption began with very little warning—at the time, there were no seismic or other monitoring networks operating in Eritrea. As we show here, seismometers in Ethiopia, Yemen and Djibouti did record seismicity associated with the volcano, but they did not provide information of operational value at the time of the eruption. However, felt earthquakes on the volcano, occurring over several hours before the eruption, did prompt a rapid evacuation of settlements, notably those within Nabro's caldera. This likely saved many lives. Nevertheless, seven people were killed and three others injured during the eruption. About 12,000 people were ultimately displaced and cared for in temporary camps in the region (Solomon 2012).

In broad terms, Nabro is sited in an extensional setting in the geographical area of Afar, close to the Mesozoic crustal block of the Danakil Alps. It reaches a maximum elevation of over 2,200 m above sea level and has an 8-km-wide summit caldera with associated ignimbrites (Wiart and Oppenheimer 2005). The axis of the Bidu Volcanic Massif, i.e. the alignment of volcanic centres, bears NE–SW, in contrast to the NW–SE trend of the Red Sea and axial volcanic ranges in northern Afar. The Nabro eruption occurred during a period of apparently heightened tectonic and volcanic activity in the Afar/Red Sea region (Yirgu et al. 2014; Jónsson and Xu 2015), which began with the Manda Hararo rifting episode and Dabbahu eruption in Ethiopia (in 2005; e.g. Grandin et al. 2009). Further, there have been eruptions of Dalaffilla (in 2008; Keir et al. 2013) and Erta 'Ale (2010; Field et al. 2012) on the Erta 'Ale range in Ethiopia and Jebel al Tair (2007) and Zubair (2011–13) in the southern Red Sea (Xu et al. 2015).

The 2011 Nabro eruption began shortly before midnight local time on 12 June following several hours of ground shaking. The preceding and accompanying seismicity is the first of note in global catalogues for this part of Afar. Remarkably, the eruption produced the largest stratospheric aerosol

perturbation since the 1991 Pinatubo eruption (Fromm et al. 2014). This reflected the substantial SO₂ yield, estimated at around 4.5 Tg for the first 15 days of the eruption by Theys et al. (2013), and penetration of the plume into the lower stratosphere (Fromm et al. 2014; Clarisse et al. 2014). It has been argued that Nabro's is one of a number of eruptions in the twenty-first century to have had a detectable effect on near-global troposphere temperature, sea surface temperature and precipitation, thereby contributing to the apparent 'hiatus' in tropospheric warming since 1998 (Santer et al. 2015).

We present here a preliminary synthesis of the nature and mechanisms of the eruption based on multiple observations, obtained via means of satellite remote sensing, regional seismic and infrasound data, fieldwork and petrology. In particular, we aim to reach a first-order characterisation of the eruption in respect of the geometries and kinematics of associated intrusive activity, the compositions of lavas and tephra, the regional context of extensional tectonics and the hazards and humanitarian response. We hope to stimulate further work on the eruption and, more widely, on volcanism in the region. Our study is also relevant to any future work that might consider the geothermal potential of Nabro volcano.

## Methods

We outline here the principal datasets and techniques used in this multi-disciplinary study. In addition to those described below, we have incorporated datasets obtained from operational Ozone Monitoring Instrument (OMI) SO₂ retrievals (Yang et al. 2007; Sealing 2013) and a Comprehensive Nuclear Test-Ban Treaty Organisation (CTBTO) infrasound station (IS32) in Nairobi, 1,708 km from Nabro (Fee et al. 2013). Additionally, we report in the Overview of the petrology of the 2011 lavas and tephra section a digest of conventional XRF analyses of the whole-rock compositions of samples of lava and tephra erupted at Nabro in 2011.

## Seismology

At the time of the eruption of Nabro in 2011, there were no operational seismometers in Eritrea (triggered permanent stations in the country were out of service) nor was there any routine surveillance of the volcano. The first broadband instruments in the country were deployed (serendipitously) on 23 June 2011, 11 days after the eruption began. This deployment of six sensors was supported by a collaborative project between several of the authors, with the aim of studying regional mantle structure (Hammond et al. 2013). Fortunately, in the months preceding the eruption, seismometers were operating elsewhere in Afar, and we obtained data from ten stations in Ethiopia (Hammond et al. 2013; Barnie et al. 2015b), eight in Djibouti (for which we received arrival time picks) and one in
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Seismic stations used in our study (inverted blue triangles and red diamonds indicate temporary and permanent stations, respectively) and regional seismicity (red dots) from the USGS catalogue for the period</figcaption></figure>

1960–2011. The green triangle locates Nabro. Inset map shows the volcanoes in the Bidu Volcanic Complex (coverage of the red box in the larger map)

Yemen. Combined, the data from these networks enable us to characterise the seismicity before, during and after the Nabro eruption (Fig. 1). We focus our analysis on events spanning the period from 23 February to 17 September 2011 (3.5 months before the eruption to 3 months after). We note that in August 2011, eight broadband stations were deployed around the volcano (Hamlyn et al. 2014). These provided a valuable 14-month dataset that is especially suited to studies of the relaxation of the system after the eruption. We did not include these data in the present analysis.

We used HYPOINVERSE-2000 (Klein 2002) to locate events from measurements of P- and S-wave arrival times, drawing on a 1D velocity model previously used for investigations in the Afar region (Keir et al. 2009). This model is based on the crustal structure deduced from wide-angle controlled-source seismology (Makris and Ginzburg 1987). We located earthquakes using a minimum of five stations. Typical

horizontal errors are <10 km, but depths are only poorly constrained. When we added data from the more proximal Eritrean stations (available from 23 June 2011), errors are reduced, but locations are systematically shifted ~2 km to the NE. We also relocated events by double-difference methods (Waldhauser and Ellsworth 2000), reducing relative location errors (relative to the master event errors) to between 2 and 50 m, with only six events having relative location errors exceeding 20 m. Local magnitudes were calculated using amplitudes measured on the horizontal components of simulated Wood Anderson displacement seismographs and the distance correction for the Main Ethiopian rift (Keir et al. 2006, 2011). Seismic moment release was estimated using empirical relationships between $M_L$, $m_b$ and $M_o$ (e.g. Kanamori, 1977; Scordilis 2006) as previously applied in a study of Dallol volcano (Nobile et al. 2012). We used Real-time Seismic Amplitude Measurement (RSAM) to follow
temporal changes in seismicity (Endo and Murray 1991). This
was calculated from average amplitude of a 10-min moving
window for the vertical component at selected stations.

Spaceborne synthetic aperture radar imagery

Data

We have analysed spaceborne synthetic aperture radar (SAR)
imagery to characterise surface changes associated with the
Nabro eruption. We examined backscatter images acquired
through time to track the evolution of lava flows and tephra
deposits, and of large-scale topographic changes, and to mea-
sure surface motion based on interferometric processing and
image correlation. Data were acquired from the X-band SAR
sensors TerraSAR-X (operated by DLR) and COSMO-
SkyMed (operated by ASI). The satellites provide SAR im-
ages with nominal spatial resolutions of 2 m for TerraSAR-X
('Stripmap mode') and about 1 m for COSMO-SkyMed
('Spotlight' mode).

Our differential analyses require identical acquisition ge-
ometry for the SAR images, to avoid geometric distortions
and artefacts arising from differences in illumination direc-
tions. Pairwise combinations are thus restricted to intervals
of 11 days for TerraSAR-X and 9 days for COSMO-SkyMed.

Mapping volcanic deposits

X-band SAR images can map geological structures down to the decametric scale. However, geometric distortions associated with complex topography complicate their interpretation. At large scale (from tens of pixels across), the average amplitude of the SAR images is mostly sensitive to the local illumination angle and therefore to slope. Due to the inclined line-of-sight direction of SAR images (usually about 30° from the vertical), systematic geometric distortions and radiometric variations are evident in scenes with topographic relief: while slopes facing towards the sensor appear bright and compressed, slopes inclined away from the satellite are darker and stretched. SAR amplitude images thus provide a qualitative characterisation of topographic changes (Fig. 2).

At the pixel scale, SAR amplitude also depends on back-
scattering characteristics of the surface. The latter are deter-
mined by the dielectric properties of the surface and surface
roughness at the scale of the radar wavelength (about 1 cm in
our case). We differentiated the SAR amplitude images ac-
quired before and after the eruption to identify the location
and nature of deposits (Fig. 2). The newly emplaced lava
flows generally appear bright, likely due to high surface
roughness. In contrast, decreased backscattered amplitude in-
dicates ‘smoothing’ of the topography by tephra deposits. A
significant decrease of the SAR amplitude can be seen in a
20×20-km² area to the SW of Nabro.

Ground deformation analysis

We have analysed ground deformation associated with the
eruption in two ways. First, the line-of-sight component of
ground displacement is retrieved via interferometric SAR
(InSAR) analysis of pairs of SAR images spanning the erup-
tion. We used the NSBAS software, a modified version of
ROI_PAC adapted to processing of time series SAR images
(Rosen et al. 2004; Doin et al. 2011). The InSAR method
requires preservation of surface properties of the ground—it
thus fails where tephra deposits or lava flows have accumu-
lated between image pairs. Thus, we only obtain meaningful
InSAR measurements beyond Nabro’s caldera rim. InSAR is
sensitive only to the component of ground motion that is di-
rected towards or away from the satellite, in the line-of-sight
(LOS) direction, about 30° from nadir. Accordingly, we com-
puted interferograms on descending and ascending overpasses
(looking from ESE and WSW, respectively) to evaluate the
importance of horizontal and vertical components of ground
motion (Fig. 3). For SAR data captured after the eruption
began, we selected scenes that were acquired within a short
time of each other to ensure that the InSAR pairs captured
comparable deformation fields. We focus here on two pairs
of SAR images acquired by TerraSAR-X: 11 June 2010 and 1
July 2011 for the descending pass and 14 May 2010 and 6
July 2011 for the ascending pass. An analysis of ground de-
formation spanning the period July 2011 to September 2012
(after the main phase of the eruption) is provided by Hamlyn
et al. (2014).

In addition to InSAR, we performed a sub-pixel correlation
of SAR amplitude images using the ‘ampcor’ routine of
ROI_PAC (Fig. 3). This resolves two components of the
ground motion along the native coordinates of the SAR im-
ages: horizontal motion along the satellite track (the ‘azimuth
offset’) and motion towards or away from the satellite LOS, in
the slant direction (the ‘range offset’). Displacements exceed-
ing around 10 % of the pixel size can be resolved with this
method, yielding a detection threshold of 10–20 cm, suitable
for measuring large deformation (Grandin et al. 2009). The
component of displacement sensed by sub-pixel correlation in
the slant range direction is the same as for InSAR but
does not saturate in zones of large deformation.
Furthermore, in contrast to InSAR, sub-pixel correlation
can provide meaningful measurements in tephra-covered
areas so long as underlying landforms are still apparent.
However, it is difficult to distinguish actual motion from
apparent motion due to the addition of material at the
surface. Such spurious effects are likely to particularly
contaminate the range offsets, as the slant range compo-
nent is nearly vertical. This should only suggest apparent
uplift. In contrast, azimuth (horizontal) offsets, when mea-
surable, should be little biased by resurfacing by tephra
deposits.
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Left panel shows a differential analysis of backscattered SAR amplitude for the Nabro area based on analysis of two TerraSAR-X scenes from 11 June 2010 (pre-eruption) and 1 July 2011. The normalised difference of the amplitude images is computed using the formula: $\text{Delta\_amp} = 2 \times \frac{(\text{amp}_1 - \text{amp}_2)}{(\text{amp}_1 + \text{amp}_2)}$. Red areas indicate decreased SAR amplitude and topographic smoothing (at the cm</figcaption></figure>

scale), while blue areas correspond to increased SAR amplitude and increased surface roughness. The upper-right panel shows the Nabro caldera and 2011 lava flow. The lower-right panel shows the corresponding topography. Coloured triangles indicate different generations of eruptive vents identified from time series analysis of SAR and thermal imagery (Fig. 10)

Finally, using two satellite overpass directions (ascending and descending), four independent components of the ground motion can be resolved (motion in the slant range direction is determined simultaneously by InSAR and sub-pixel correlation). These four independent components permit, after solving an over-determined set of linear equations for each pixel, retrieval of the east-west, north-south and up-down components of ground displacement (Fig. 4).

## Satellite optical imagery

We use satellite images from several sources to track the eruption by its thermal emissions and geomorphologic changes on the ground. These include the following four instruments: the Spinning Enhanced Visible and Infrared Imager (SEVIRI), a multispectral imager aboard Eumetsat's Meteosat geostationary weather satellites; the Moderate-resolution imaging spectroradiometer (MODIS), a multispectral imager aboard NASA's Aqua and Terra satellites; the Advanced Very High-Resolution Radiometer (AVHRR), another multispectral imager aboard various NASA and Eumetsat satellites; and the Advanced Land Imager (ALI), aboard NASA's Earth Observing 1 (EO-1) satellite. SEVIRI has a spatial resolution of about 3 km at nadir, with 96 acquisitions per day; MODIS

has a spatial resolution of 1 km at nadir with approximately four acquisitions per day; AVHRR has a spatial resolution of 1 km with a few acquisitions per day; and ALI has a spatial resolution of 30 m with acquisitions typically every few days to weeks when tasked to view an event, either manually or automatically.

We interpret the MODIS, AVHRR and ALI images visually. They contain useful spatial information but are acquired infrequently. In contrast, the SEVIRI images enable near-continuous monitoring albeit at comparatively coarse spatial resolution. We extract the volcanogenic thermal signal from the SEVIRI dataset following Barnie and Oppenheimer (2015a) as follows. Independent component analysis (ICA) is applied to the time series of radiances associated with each non-saturated band 4 (3.9 µm) pixel, in a 9×9 pixel widow surrounding the eruption, to find statistically independent contributing sources. The magnitude of the contribution of each source to each pixel is found, again excluding saturated hotspots. Volcanic sources are identified by a high temperature event (HTE) threshold (the product of the absolute skewness of the source and per-pixel contributions). The per-pixel volcanogenic radiances are then found by summing the outer product of the source and contributions for all sources above the threshold. In this way, the image cube is decomposed into
Ascending pass (2010/05/14–2011/07/06)

<figure><img src="image_5.png" /></figure>

**Fig. 3** Comparison of observed (*Data*), predicted (*Model*) and residual (*Residual*) deformation for the various components of ground motion deduced from ascending (**left panel**) and descending (**right panel**) SAR data. The **top row** corresponds to InSAR, the **middle row** to range offsets and the **lower row** to azimuth offsets. The **thick-coloured arrows** show the direction of sensed ground displacement: **red** for line-of-sight (LOS) and **blue** for along-track (AZ). InSAR and range offsets are sensitive to motion only in the LOS direction. The LOS vector makes an angle of 29° with the vertical for the ascending pass (34° for the descending pass). For clarity, only the horizontal component of the LOS vector is shown by the **red arrow**. In contrast, azimuth offsets only provide information on the horizontal component of ground motion in the direction of the satellite track (AZ). In each panel, motion in the direction indicated by the **thick-coloured arrow** has negative sign (**cool hues**), whereas motion in the opposite direction has positive sign (**warm hues**). Note the colour scale is common to all components. **Thin arrows** represent the modelled horizontal component of ground motion derived from the elastic inversion. The **shaded area** SW of Nabro received substantial tephra fallout and was excluded in the inversion. The **blue rectangle** represents the surface projection of the normal fault deduced from the model (Fig. 14). The **red rectangle** shows the surface projection of the dike, and the **orange circle** the surface projection of the pressure source. See Fig. 12 for a perspective view

Descending pass (2010/06/11–2011/07/01)

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

**Fig. 4** As preceding figure but for ground motion resolved in east–west, north–south and up–down components. The **left panel** corresponds to a reconstruction of the 3D components of ground motion by combining the six available components of ground motion, among which only four yield independent directions. The **colours** represent the vertical displacement, while **arrows** indicate the horizontal component
a number of radiance sources, the non-volcanic sources excluded, and the image cube reassembled with only non-saturated volcanic sources. For events with many saturated pixels, the correction may underestimate the true radiance but not as substantially as the truncated saturated values. Estimating the contribution of cloud contamination to any given pixel is challenging, and interpretations of the time series are necessarily qualitative: for instance, a drop in radiance could result from either occlusion of hot material or declining lava effusion. We note also that the vent and lava flow were frequently obscured by the volcanic plume.

## Fieldwork

Fieldwork was conducted between 8 and 15 October 2011 by a subset of the authors (BG, CO, JOSH, EY and GK), based first at Afambo and later at Wade. Given the limited time, the extent of the lava flow fields and tephra fall, and the proximity of the international border, a comprehensive sampling strategy was impractical. However, lavas and tephra were sampled over a wide area within and beyond the caldera.

Meaningful estimates of lava flow thicknesses were not possible given the dimensions of the flow field. The flow surfaces were well above ambient temperatures in places and were still emitting fumes. The longest flow, which had traveled beyond the caldera, had overrun cultivated and settled areas. It was sampled towards its front. The vent region was approached from the northern part of the caldera. Within the caldera, three shorter flows were sampled.

A second field mission was conducted in October 2012 to decommission the temporary seismic network. It afforded time for a reconnaissance of the vent region from the southern part of the caldera. The two missions provide some relevant observations but were not focused on obtaining specific ground truth with which to validate interpretations based on our analysis of satellite and other monitoring data.

## The eruption: observations and models

In this section, we report key observations (seismology, satellite remote sensing and ground surveys) and explanatory models for intrusive and eruptive processes involved in the 2011 eruption.

### Seismicity

The area around Nabro has hitherto appeared to be seismically quiet. No major earthquake was reported by international or regional centres (USGS or ISC) prior to the eruption, with the exception of a $4.5-M_W$ event on 31 March 2011 (i.e. just over 10 weeks before the eruption). This was the first earthquake ever reported by the USGS for the area (Fig. 1). The epicentral

coordinates given by the USGS for this event are 13.17° N, 41.83° E (horizontal error > 20 km), i.e. about 25 km southeast of Nabro, with an origin time of 18:33:38 UTC.

On 15 May 2011, a scientific team from Asmara visited the area (Ogubazghi et al. 2014) to investigate the effects of this earthquake and to interview eyewitnesses. It was found that most houses within the Nabro caldera were damaged, with the severity ranging from cracked walls to partial collapse (Fig. 5). The event had also induced landslides and, according to local people, killed many domesticated animals including sheep and goats. Small tensional cracks were still present on the ground (Ogubazghi et al. 2014). The length of these cracks reached up to 10 m but, reportedly, they had extended several tens of metres immediately after the event. The cracks trended NNW-SSE, parallel to the direction of plate opening determined from plate kinematic and moment tensor analyses (e.g. McClusky et al. 2010).

While residents of Siroru, a village situated on the southwest floor of Nabro's caldera, clearly felt the shaking, people in the vicinity of Maebele and Wade (close to the USGS-located epicentre) were not awoken by the event. Moreover, damaged houses were observed around Siroru but not in Wade or Maebele. A maximum intensity of VII was estimated for the Siroru area and of less than IV in Wade and Maebele (Ogubazghi et al. 2014). Based on these observations, the epicentre was relocated to 13.33° N, 41.68° E, placing it within the caldera (Ogubazghi et al. 2014). This is close to our relocation of the earthquake, using data from the stations in Ethiopia and Djibouti, to 13.35° N, 41.69° E (horizontal error of ±1.2 km), with a magnitude $M_L$ 4.8 and time of 18:33:37 UTC on 31 March 2011 (yellow star in Fig. 6a).

From analysis of the regional data, we detected just six other events in the Nabro region before 12 June, with $M_L$ ranging between 2.8 and 3.9. We did not identify any earthquakes in the days preceding the eruption (the last one was on 29 May) earthquakes with magnitudes below our detection threshold of $M_L \sim 2.1$ may well have occurred.

The eruption was preceded by 5 h of earthquake swarms, with the first event at 15.37 UTC on 12 June 2011. Two further notable events occurred at 20:32 and 21:03 UTC, the earlier of the two likely very close to the onset of the eruption (Fig. 7).

The most seismically active days were 12 and 17 June, with 79 events ($>M_L$ 2.0) on 12 June and 28 events on 17 June (totalling 55 % of the recorded events). Local magnitudes during this period vary between $M_L$ 2.7 to 5.9, with ten events exceeding $M_L$ 5. The largest event ($M_L$ 5.9) occurred on 17 June; the second largest, at $M_L$ 5.8, took place on 12 June. On 12 June, seismicity occurred beneath both Nabro and Mallahle volcanoes. In contrast, on 17 June, while there was little seismicity at Nabro, Mallahle remained seismically active (Fig. 6b). This suggests a linkage between the two volcanoes, either via their magmatic systems or associated with stress
field changes. This shift in seismicity to the SW precedes the availability of data from the Eritrean stations, so it is not an artefact of changing station coverage. Mallahle remained seismically active in the months after the eruption, as seen from analysis of data from the local seismic network installed in August 2011 (Hamlyn et al. 2014).

While only ten earthquakes were recorded between 12 and 17 June, RSAM analysis of seismograms from the regional stations indicates sustained tremor (Fig. 8a). Spectrograms for the closest stations, KOZE and FAME, reveal a signal dominated by energy at frequency bands <2 Hz (Fig. 8b). The tremor was very high on the first day of the eruption but reduced during the next 2 days. On 16 June, the tremor almost stopped, but it picked up abruptly with the onset of a second period of intense seismicity on 17 June. After this, tremor persisted until 14 July with many fluctuations (Fig. 8c). The tremor was observed at many stations with different azimuths precluding interpretation of the observed signal as a path effect (e.g. Cote et al. 2010). For the period 18 June to 17 September, we located 93 events (43 % of the total events in our analysis). Events were detected on most days, with 3 days marked by earthquakes sized >$M_L$ 4.

Looking at the double-difference relative locations, we identify a NE–SW trend (i.e. parallel to the Nabro–Mallahle axis) in epicentres beneath Mallahle on 17 June (Fig. 6b). In contrast, after 17 June, the distribution of epicentres on Mallahle follows more of a NW–SE (Red Sea) trend (stretching as far as Sork ‘Ale volcano; Fig. 6c).

## Observations from spaceborne optical and SAR imagers

The Nabro caldera is up to 8 km across (north to south). In its elevated centre, prior to the 2011 eruption, there were two craters, one about 400 m deep and 2 km across and the other, partly recessed into its SW rim, up to 1.4 km across and over 200 m deep. This smaller crater, whose floor was about 100 m deeper than its neighbour and which we refer to as the SW pit, proved to be the main focus of the 2011 eruption. On its floor were the remnants of a small low-aspect-ratio crater about 140 m across. The inner walls of the SW pit, especially on its east side, were draped with what appear from satellite images to have been scoriaceous lavas, perhaps resulting from drain-back of lava. These lavas are less evident on the western wall of the SW pit, which showed signs of substantial hydrothermal alteration.

The time series of SEVIRI images indicates that the eruption started between 20:27 and 20:42 UTC on 12 June, i.e. shortly before midnight local time as confirmed by eyewitnesses. This is evident in Fig. 9a, b, which shows a measure of thermal radiation for two pixels, one representing the final lava flow front and the other the vent region inside the caldera. The signal is ‘spiky’ due to the intermittent presence of thick plume and cloud that obscure the hot material at the surface.

Figure 9c shows the first SEVIRI image to register the eruption, showing an incipient thermal anomaly and plume, placing eruption onset between the time of acquisition of this image and its predecessor (20:42 and 20:27 UTC, respectively). While the SEVIRI images usefully reveal the structure and dispersal of the plume, the thermal anomalies on the ground are obscured until 15 June.

The first high-spatial resolution image we have is an ALI scene acquired at 07:32 UTC on 14 June, i.e. about 35 h after the start of the eruption. A plume rich in condensed steam drifts to the west. It covers most of the caldera, hiding the vent area from view. The next spatially detailed scene, a COSMO-SkyMed SAR image acquired on 18 June, 5 days after the eruption onset, reveals substantial changes compared with a pre-eruption scene—(i) the innermost pit crater has been largely infilled by a tephra cone, (ii) several lava flows have issued from points on and around this cone, (iii) the longest of these flows has breached the open SW margin of the caldera and then travelled NW and (iv) a discontinuous alignment of small pits has opened up along a NW–SE trend, reaching about 2 km from the vent region. Inspection of the pre-eruption topography indicates that the main lava flow (labeled (c) in Fig. 10) was sourced at the lowest point on the SW pit rim (about 1,515 m above sea level), on its western side.

To follow events of the first days of the eruption, we have to turn to lower-spatial resolution sensors including MODIS and AVHRR. Figure 9d shows thermal images from these instruments covering the period from 14 to 17 June. The interval can be subdivided into three periods. In period 1, the plume dominates the SEVIRI signal, and the MODIS/AVHRR images show a small thermal anomaly confined to the crater. In period 2, the plume margins develop thermal anomalies at 3.9 µm, which could be due to forward scattering of infrared radiation from the advancing lava flow. Minor

<figure><img src="image_8.png" /><figcaption>Fig. 5 Building damaged by the 31 March 2011 earthquake. Better constructed walls, as seen in the background, sustained only cracking</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Fig. 6** a Seismicity around Nabro for the period 23 February to 17 September 2011. Locations for the 31 March event are shown by the blue star (USGS location based on distant stations of the Global Seismic Network) and yellow star (our location based on stations in Ethiopia and Djibouti). Colour codes display temporal evolution, events</figcaption></figure>

with horizontal errors exceeding 2 km are *coloured grey* and symbol size represents magnitude (*M*<sub>L</sub>) classes (see legend). Plots (b) and (c) show double-difference relative locations for events occurring between 12 and 17 June and between 18 June and 17 September 2011, respectively

thermal anomalies are also registered in the lava flow front pixel indicating that it had covered most of the distance by the evening of 15 June. During period 3 (from the evening of 16 June), the plume has weakened, revealing the thermal anomaly at the vent and showing the lava flow has reached its maximum length by this point.

Returning to the high-spatial resolution imagery, we have two more SAR scenes acquired on 20 and 21 June before the next optical images from ALI were recorded (Fig. 9e). Dramatic topographic changes occurred in the vicinity of Nabro's former SW crater. As a result, images acquired syn-

or post-eruption that are geo-coded with a pre-eruption DEM are severely distorted in this area (see Fig. 2). Therefore, we analyse changes in the native geometry of the SAR system. The alternating ascending and descending passes of the satellite further complicate this comparison (successive passes in the same direction are acquired every 11 days by TerraSAR-X). The most prominent change apparent in SAR imagery is the complete filling of the SW crater, already seen in the first image acquired after the eruption onset (20 June). Several new scoria cones are also identified. We calculate that the SW pit was filled from 1,350 m elevation (pre-eruptive crater floor) to
<figure><img src="image_10.png" /><figcaption>**Fig. 7** Time series of a daily totals of seismic events located using the Ethiopian and Djibouti stations and b corresponding seismic moment release (SMR)</figcaption></figure>

at least 1,550 m elevation (highest altitude of overflow of the crater rim), requiring a volume of lava and tephra exceeding 0.1 km³.

The ALI scenes acquired between 07:18 UTC on 24 June and 07:27 UTC on 2 July are quite spectacular (Fig. 9e), notably in the shortwave infrared (SWIR) channels. These indicate intense thermal radiation around the vent area located close to the highest point on the SW crater (on its NNE rim) and along a number of lava flows. In all these images, the source of the main lava flow remains ‘hot’ in SWIR channels, though the intensity of the anomaly has diminished by 29 June. The front of this flow appears still active on 24 June. Video taken by one of us (SB) near the flow front reveals widespread exposure of incandescent lava along the flow margins, in addition to flames leaping from engulfed trees. The strong thermal anomaly centred on the new tephra cone is suggestive of a perched lava lake filling a flared vent.

The main lava flow (labelled (c) in Fig. 10) initially travelled SW from the vent region until reaching the breached part of the horseshoe-shaped caldera. At this point, it turned NW, reaching an overall flow length of about 17.5 km. In addition, both SAR images and ALI infrared scenes reveal successive generations of smaller lava flows emplaced to the S, E and N of the vent region. With the exception of a short lava flow fed from a series of aligned vents to the SW of the old crater (flow (a) in Fig. 10), all these lava flows appear to originate from

four distinct sources at the four cardinal points of the SW pit (flows (b), (c), (d), (e) in Fig. 10). The first lava flows were sourced at points around the W (c), N (b) and S ((a) and (d)) rims of the SW crater. Vents located to the north of the rim ((b) and (e)) sourced lava flows that were confined to a pre-existing pit crater. Lava supply from vent (b) had ceased by 20 June. In contrast, the flows emitted by vent (d) were emplaced in two phases: a short flow was emitted before 20 June and was subsequently covered by a longer flow that remained active until at least 1 July. Finally, onset of effusive activity at vent (e) seems to have occurred after 20 June and continued until at least 12 July. We note that vent (c), which sourced the main lava flow, marked the lowest point of the pre-existing crater rim (1,515 m elevation). Vent (e), the second longest-lived, was located at 1,525 m elevation; vent (d), the third longest-lived, at 1,545 m elevation; and vent (b) stands at an altitude of 1,585 m. This pattern may suggest a progressive lowering of the magma column.

It is not until 29 June, when there remains only a weak ashy plume drifting to the south, that the ground surface south and southeast of the tephra cone becomes unobscured. A smaller lava flow has advanced east from the SW pit rim and then diverged with one branch heading north where it ponded in the adjacent crater and another moving south down the flank of intracaldera edifice. High temperatures are evident along levées of the main lava flow, where incandescent lava was
<figure><img src="image_11.png" /><figcaption>**Fig. 8** a Raw vertical component seismogram from station KOZE, 120 km from Nabro volcano. b Spectrogram from station KOZE. c Composite plot showing RSAM at stations KOZE (blue) and FAME (yellow), SO₂ emission rate computed from OMI data (red) and infrasound signal recorded at the Nairobi station (green). SO₂ emission rates were derived using a plume traverse technique (e.g. Theys et al. 2013), assuming SO₂ advection in a constant wind field, with wind speed derived from a trajectory model. The magnitude and temporal variation of our OMI-based SO₂ fluxes are broadly consistent with analyses of other satellite SO₂ measurements for Nabro (Theys et al. 2013). The dashed vertical line indicates the eruption onset (20:27–20:42) estimated from SEVERI images</figcaption></figure>

likely more exposed. By 2 July, the eastern lava flow had advanced further into the adjacent crater. Particularly striking in the 29 June ALI SWIR bands is the revelation that the aligned cinder cones seen in the SAR images from 18 June onwards coincide with intense thermal anomalies. The location and orientation of these vents suggest the possibility of their source being a NW–SE-striking dike (Fig. 11). These vents do not appear to have contributed significantly to lava eruption after 20 June.

Notably, one of these extends from the SE rim of the SW crater, which appears to have been substantially infilled by lava and tephra. The main flow appears to be at its maximum extent, though the ‘spillway’ at the vent that sources it remains ‘hot’, suggesting sustained lava effusion, in the shortwave infrared channels.

Comparing the SAR images acquired on 18 June and 4 July indicates that the main flow only advanced a further 400–500 m in the intervening period, reaching 750 m elevation,
<figure><img src="image_12.png" /><figcaption>**Fig. 9** Spaceborne optical observations of the Nabro eruption. a A 3.9-µm image of the saturation-corrected volcanic radiance extracted from SEVIRI at 15:12 UTC on 16 June, superimposed on GDEM topography for comparison. Low-to-high radiance is scaled from *blue-to-red*, and the final extent of the lava flow is outlined in *red*. Pixels covering the vent area and lava flow front are outlined in *blue* and *purple*, respectively. The time series of volcanogenic radiance for these pixels are shown in **b**. Note that the radiance time series are 'spiky' due to presence of cloud and plume. c The SEVIRI image that signals the start of the eruption (and providing the isolated spike in radiance seen late on 12 June in **b**, before the vent is obscured by plume). d Low-resolution images tracking the eruption onset, with bands around 4, 10 and 12 µm as RGB, respectively,</figcaption></figure>

and a total length of approximately 17.5 km. The ALI images show that between 24 June and 2 July, the main lava flow advanced only ~140 m at its front, though it likely continued to thicken and widen in places along its length during this period. Between 17 and 23 July, the SAR images reveal drainage of lava from the main flow's channel, quite close to the source and within the caldera. There has been a minor advance of the intra-caldera flow to the SE. As early as 20 June, scarps are apparent around the tephra cone, where the original rim of the pit crater is being expressed or exposed. Much more of the new tephra cone apparently collapsed later in July, though between 23 July and 8 August, no further substantive changes are evident in the SAR imagery.

By 28 September, the SWIR channels of the ALI image indicate anomalous radiation only in association with hot spots around the vent region and along the rim of the original pit crater. The lava flows appear 'cold' at these wavelengths.

in each case. Saturation of AVHRR channels results in loss of a strip of pixels, leading to the banding effect seen in some images. During period 1, the thermal anomaly is restricted to the SW pit in the Nabro caldera, while during period 2, a weak thermal anomaly around the plume margins may result from forward scattering of thermal emission from the advancing lava flow indicating breaching of the SW pit. Period 2 is associated with weak thermal anomalies in the SEVIRI time series seen in **b**. Period 3 follows the first break in the plume, revealing that lava has reached almost its final extent by this point. e High-resolution shortwave infrared images of the vent region during the later stages of the eruption. *Red polygon* marks the outline of the original SW pit

This pattern persists in the 6 October ALI scene but, 10 days later, the anomalies around the vent region in the SWIR bands look less intense. This is consistent with observations made by several of us during fieldwork carried out on the volcano between 8 and 15 October. We observed incandescent vents in the vent area from the caldera rim on the night of 12 October and a weak plume drifting across the caldera during the daytime.

A prominent signature of the initial flows emitted from vents (a) and (c) (Fig. 10) is that their texture appears to be smoother in the first image acquired after the start of the eruption (20 June) than seen subsequently. The smooth texture of early lava flows seen in the SAR images could reflect lower surface roughnesses of the lavas or, more likely, tephra fallout from the early phases of the eruption (the effects of which are readily apparent elsewhere within the caldera; Fig. 2). This would suggest that the output of ash diminished substantially
TerraSAR-X (ascending)

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /><figcaption>**Fig. 10** Time series of SAR images focused on the source region of the flows, with interpretations of lava flow emplacement and topographic changes. Upper and lower panels show ascending and descending geometry, respectively (in the native geometry of the sensor to avoid artefacts caused by erroneous orthorectification due to syn-eruptive topographic changes; see Fig. 3). The images have been rotated to</figcaption></figure>

facilitate interpretation. In each SAR image, the far range is at the bottom, and the near range is at the top. Note significant distortions and radiometric changes associated with steep slopes. The central panel shows lava flow emplacement through time, with the dates of the images used for interpretations (dates indicated in grey were unavailable)

after 20 June, as later lava flows retain their ‘blocky’ appearance thereafter.

In the SAR images acquired after 20 June, we see that lava eruption from sources (d), (c) and (e) (Fig. 10) continued up to early to mid-July. This is consistent with the volcanic thermal

radiance registered in the pixel situated above the vent, which persisted until around midday on 16 July, before decaying rapidly to low values before midnight (Fig. 9b), marking the cessation of intense eruptive activity. The radiant signal from the flow front decayed significantly by early July, suggesting

<figure><img src="image_18.png" /><figcaption>**Fig. 11** Photograph of part of the fissure system taken in 2012. View approximately to the W—a small cone, formed in the 2011 eruption, is apparent on the left and the front of lava effused from vent (e) (Fig. 10) apparent on the right. Behind is the modified cone of the SW pit crater</figcaption></figure>
waning advance. Compared with the 12 July SAR image, the 17 July scene reveals a set of ring faults or fractures associated with subsidence of the newly emplaced central cone. The geometries of these fault traces recapitulate the trace of the former SW pit rim, suggesting localisation of brittle deformation along the interface between the previous crater walls and the new infill of tephra and lava. Lava eruption seems to have ceased in conjunction with this substantial collapse of the material, although minor changes prevailed near vent (c), possibly arising from drainage of lava already within the channel or slope instability.

Both the SAR and optical images clearly reveal the fallout from the eruption, fanning out in a 130° sector from the vent area (spanning SSE to WNW). In the SAR images, this is apparent from a decrease in backscattered intensity (area in red in Fig. 2). This reflects a smoothing of the topography at the scale of the radar wavelength (about 1 cm) and decrease of the energy backscattered towards the antenna. The most significant change in radar intensity occurred before the first acquisition (20 June), which suggests that the bulk of the ash was deposited in the first 6 days of the eruption. The area of tephra fallout evident in the SAR imagery covers an area of about 700 km². It is not possible to assess the thickness of tephra on the ground solely from the available satellite images.

## Satellite geodesy

The decomposition of ground displacement into horizontal (east-west and north-south) and vertical components derived from SAR imagery reveals a complex pattern of ground motion (Fig. 4). At the large scale, the horizontal displacement field indicates an overall NE-SW extension, concurrent with NW-SE compression. Maximum horizontal displacements reach 0.75–1 m beyond Nabro’s caldera rim. The vertical component of ground displacement indicates widespread subsidence, reaching an average of 0.2 m to the N, E and SE of Nabro. Localised subsidence, reaching up to 0.5 m, is apparent to the NW of the caldera rim within an area elongated along a NW-SE axis. Displacements retrieved within the caldera and to the SW of the vent region are particularly complex. The most prominent feature is localised uplift reaching about 0.75 m within a circular area, 5 km in diameter. The vents lie on the edge of this zone. Further south, apparent subsidence of Mallahle caldera exceeds 0.5 m. Horizontal displacements are directed to the SW and exceeded 1 m. This area of complex deformation coincides with a decorrelation of the InSAR data (Fig. 4) and a decrease of the backscattered SAR amplitude (Fig. 2). This suggests that thick tephra deposits cover the area and that the retrieved ground motion is at least partly an artefact. We therefore neglect this area in the quantitative analysis of ground deformation that follows.

For a first-order interpretation of the deformation source, we performed a series of forward elastic models to reproduce

the main features of the ground displacement field in well-resolved sectors of the volcano and its environs. These were based on models from Okada (1985) and Mogi (1958) to predict surface deformation induced by rectangular dislocations or isotropic point sources, respectively. We tested various combinations of faults, dikes and pressure sources. We found three elementary sources sufficient to explain the pattern and amplitude of ground deformation. Our preferred model includes a NW-SE dike, responsible for the NE-SW horizontal dilation pattern observed to the NE of the volcano. This is consistent with the orientation of vents observed within the crater (Observations from spaceborne optical and SAR imagers section). A negative pressure source beneath the vent area explains the overall subsidence signal, as well as NW-SE contraction. The net NE-SW dilation results from the predominance of dike-induced dilation over Mogi-induced contraction in the direction perpendicular to the dike. In contrast, the dike does not produce significant horizontal motion beyond its tips, so that the effect of the negative pressure source accounts for the NW-SE-oriented pattern of contraction observed to the SE and to the NW of the volcano. Finally, a shallow NW-SE-striking, NE-dipping normal fault reproduces the localised subsidence seen to the NW of Nabro.

We then used these three sources of deformation to initialise a nonlinear inversion of the surface deformation that takes the six displacement fields as an input (InSAR, SAR range offset and SAR azimuth offset, for ascending and descending passes). We followed Tarantola and Valette (1982) to optimise the geometry (location, size, dip, depth) and kinematics (slip magnitude and direction or magnitude of pressure change) of the three elementary sources. Not all the parameters could be constrained due to incomplete spatial coverage of the data and various tradeoffs. The best-fit model is shown in perspective view in Fig. 12. This preferred solution includes a negative pressure source equivalent to a volume loss of 0.07 km³, situated at 6 km depth. This depth coincides with the depth of a deflating pressure source deduced from analysis of ground deformation derived from InSAR time series (spanning 1 July 2011 to 5 October 2012), and with a cluster of micro-earthquake sources 5 km below sea level, located using the local seismic network operating from 31 August 2011 (Hamlyn et al. 2014). The modelled dike is vertical, strikes at 320° N, is 8-km-long and extends from 0.3 km down to 4.3 km depth. It corresponds to a volume inflation of 0.04 km³, which is comparable to the volume lost from the underlying pressure source. As noted, the strike direction of the dike is roughly parallel to the alignment of the vents identified from high-resolution imagery (Observations from spaceborne optical and SAR imagers section). This trend is perpendicular to the direction of regional spreading associated with motion of the Arabian and Nubian plates.

Lastly, the normal fault strikes at 320° N, slips by 0.7 m, has a dip angle of 55° and extends from 0.8 km down to
<figure><img src="image_19.png" /><figcaption>**Fig. 12** Model deduced from elastic inversion of the SAR-derived surface displacement field. The *left panel* shows the surface projection of the three elementary sources used in the model in map view. The *centre*</figcaption></figure>

and right panels show the same sources in a 3D perspective view, looking from the SE and SW, respectively

3.2 km depth. The seismic moment associated with this fault is $4.3 \times 10^{17}$ N m, which is equivalent to an $M_w$ 5.8 earthquake (assuming a shear modulus of $3.3 \times 10^{10}$ Pa). This is comparable to the sum of the two most energetic earthquakes that occurred on the day of peak seismic activity on 12 June, at 20:32 and 21:03 UTC, with magnitudes $M_w$ 5.6 and $M_w$ 5.4, respectively ($2.91 \times 10^{17}$ and $1.42 \times 10^{17}$ N m according to gCMT). The focal mechanisms of these two earthquakes indicate normal faulting striking 280–290° N, albeit with a substantial non-double-couple component (Fig. 13). Given the uncertainties in the analysis, this is broadly consistent with the 320° N direction determined from our interferometric SAR analysis.

Overall, the elastic modelling suggests that the Nabro eruption was associated with a shallow dike intrusion. This dike may have linked a deeper source to the surface eruption. The onset of seismicity suggests the dike intruded on 12 June. It triggered slip on one or several normal faults running parallel to the dike at shallower depth, as observed elsewhere in regions undergoing local or regional extension (e.g. Rubin and Pollard 1988; Rowland et al. 2007). It is evident from our InSAR analysis that faulting breached the surface on the NW flank of Nabro. Figure 14 shows the co-eruptive displacement along a profile perpendicular to the inferred normal fault. Motion away from the satellite to the NE of the surface trace of the fault clearly points towards a normal fault dipping to the NE with most of the vertical displacement concentrated in the hanging wall. Fault geometry in the vent region cannot be readily identified due to resurfacing by lava flows and tephra deposits.

Dike intrusions are common during volcanic eruptions. However, the broad scale of fringes in the interferograms indicates that the dike opened permanently with some magma retained in the crust (i.e. not all was erupted). This inflation of the dike suggests a response to the regional extensional stress field and a degree of tectonic influence on the eruption. However, the relative importance of tectonics versus magmatic processes in nucleating, feeding and maintaining the

inferred dike remains unclear. It is interesting that the initial dike injection beneath Nabro followed a NW–SE, Red Sea trend (Satellite geodesy section), while on 17 June, the seismicity shifted to Mallahle and defined a NE–SW trend (Fig. 6b). Seismicity beneath Mallahle then switched back to more of a NW–SE, Red Sea trend, stretching as far as Sork ‘Ale volcano (Fig. 6c). This may indicate a perturbation of the stress regime in which the direction of maximum horizontal stress changed. A closer study of earthquake source mechanisms and seismic anisotropy may provide further evidence and insights.

The residual displacement field determined after subtracting the best-fit model reveals several other notable features. Beyond the masked area (shaded zones in Figs. 3 and 4), the residual horizontal and vertical displacements are generally small. In contrast, a residual uplift, reaching about 1 m, is seen to the SW of the vent region. This is likely the result of tephra accumulation in the first days of the eruption. A distinct negative anomaly, corresponding to subsidence, is also apparent within the caldera of Mallahle volcano. This residual signal might correspond to deflation of a reservoir below Mallahle, induced by activity of the neighbouring magmatic system. This activation of Mallahle is seen also in the seismicity, beginning 17 June (Seismicity section) and continuing in late 2011 (Hamlyn et al. 2014).

While our model captures the main features of the ground deformation field, we emphasise that it is non-unique for several reasons: (i) the retrieved surface displacements do not span exactly the same observation intervals, (ii) spatial coverage of the data is patchy due to the tephra fall deposits, (iii) sub-pixel correlation results are fairly noisy around the volcano, likely due to slope effects and relatively adverse surface conditions for X-band SAR imaging (compared with the Afar plains), (iv) the deformation field is particularly complex and (v) the elastic inversion problem is inherently non-unique. Thus, alternative solutions based on a different set of assumptions could be proposed. For instance, Hamlyn et al. (2013)
<figure><img src="image_20.png" /><figcaption>**Fig. 13** Moment tensors for all events with moment magnitudes ($M_w$) exceeding 5.0 and detected near Nabro during June 2011 and recorded in the Global Centroid Moment Tensor (gCMT) catalogue (Dziewonski et al. 1981; Ekström and G.M. and Dziewonski, A.M. 2012). The plots</figcaption></figure>

show the full moment tensor of each event, with *grey lines* showing the fault planes of the double-couple component of the moment tensor. Strike, dip and rake of the two nodal planes are indicated

explained the first-order contribution to the observed deformation as a NNE-SSW-trending, left lateral, strike-slip fault located at shallow depth beneath Nabro volcano, consistent with the focal mechanism of one of the most energetic earthquakes of the eruption episode (09:16 UTC, 17 June 2011, $M_w$ 5.6).

## Synthesis of observations

Figure 8c integrates observations of the Nairobi infrasound station, RSAM and SO₂ emission rates computed from available OMI satellite data. Several features stand out in the time series: (i) the persistence of seismic tremor and its infrasonic

<figure><img src="image_21.png" /><figcaption>**Fig. 14** Descending (**left**) and ascending (**right**) interferograms showing the line-of-sight component of ground displacement around Nabro volcano. Blue-to-purple shading corresponds to motion away from the satellite. The thick black arrows show the line-of-sight direction. The red</figcaption></figure>

lines indicate the profiles across an inferred normal fault (surface trace in white with bars indicating the hanging wall), shown in the corresponding lower panels
equivalent (i.e. continuous vibration of the atmosphere) dur-
ing and after the main eruptive stage, (ii) approximate corre-
lations between the temporal trends in the three parameters,
(iii) a pronounced drop in RSAM to near zero around 17 June
and (iv) plateau-like temporal signals in RSAM prior to 17
June but more fluctuating trends afterwards. Most of the lo-
cated events occurred on 12 and 17 June. The largest earth-
quakes also occurred on these days.

The SEVIRI images bracket the eruption onset to between approximately 20:27 and 20:42 UTC on 12 June—the later scene shows a pronounced thermal spike over the caldera. The prominent earthquake at 20:32 UTC could signify the final dike advance that breaks the surface. The infrasound recorded at the closest station in Djibouti (264 km from Nabro) picks up abruptly at 21:44 UTC (Fee et al. 2013). For an acoustic travel time of circa 20 min from Nabro, this indicates that explosive activity either began or picked up substantially at 21:24 UTC. This may indicate that the first 50 min or so of the eruption saw predominantly effusive activity at one or more vents, although the SEVIRI images indicate a plume that rapidly gains altitude between 20:42 and 20:57 UTC.

Optical and SAR observations indicate that the main lava flow had advanced to close to its maximum extent within about 3 days of the eruption onset on 12 June. In addition, most of the tephra was emplaced within about 6 days of the onset. The persistent tremor between 12 and 16 June suggests sustained high lava effusion rates during this period. On the other hand, infrasonic signals point to continuous but fluctuating explosive activity, consistent with the evidence for airborne ash and tephra accumulation on the ground. Unique features in the infrasound signal suggest the volcanic jet was gas-rich and supersonic (Fee et al. 2013). SO₂ emissions were strong during this period, likely with efficient decoupling of gas from erupting magma at the vent(s). The comparative quiescence in activity around 16/17 June is puzzling but may correspond to a choking in magma supply to the vent(s).

The abrupt resumption of tremor on 17 June and accompa-
nying seismicity suggests renewed magma supply.
Interestingly, most of the seismicity is located beneath
Mallahle volcano, suggesting either triggering by changes in
the local stress field (e.g. Walter et al. 2014) or potentially
interconnection between magma reservoirs (e.g. Albino and
Sigmundsson 2014). The former explanation seems more like-
ly based on the observation that post-eruption seismicity (31
August–7 October 2011) beneath Mallahle had significantly
lower *b* values than that associated with Nabro (Hamlyn et al.
2014). Optical imagery indicates vigorous activity at the vent
region, likely characterised by a lava lake sourcing a steam-
rich, low-altitude plume. SO₂ was still detected by OMI but
burdens and emission rates were reduced, suggesting lower
lava effusion rates and consistent with generally low
explosivity at the vent(s). Explosive activity on 18 and 21/
22 June is apparent in the Nairobi infrasound data and

coincides with a pulse in RSAM and SO₂ output. By 29
June, the optical shortwave infrared images (ALI) indicate
diminished lava effusion. While SO₂ emissions persisted, out-
puts were further reduced.

Overview of the petrology of the 2011 lavas and tephra

The lava erupted is trachybasaltic to trachybasaltic andesite in composition—between 48–53 wt% SiO₂ and 5–7 wt% total alkalis, following a mildly alkali-enriched trend relative to rocks from the axial volcanic ranges in Afar (e.g. Barberi et al. 1972; Fig. 15). Whole-rock MgO contents are low (3–5 wt%). There are abundant phenocrysts of olivine, clinopyroxene, plagioclase and magnetite, with rare ilmenite, chromite and fluorapatite. Microlites of olivine, clinopyroxene, plagioclase and Fe-Ti oxides are present in the ground mass. Some tephra samples are glassy with abundant plagioclase microlites, while others are more crystalline. Phenocrysts show both normal and reverse zoning, indicating a complex plumbing system. This is rare in Afar and may be the result of the presence of thicker crust in the proximity of the Danakil Alps.

Olivine phenocrysts are typically compositionally zoned with fayalite-rich rims and forsteritic cores (Fo about 85). Some contain inclusions of chromite. There are also a few smaller crystals that are reverse-zoned. Clinopyroxenes exhibit diverse textural features. One population has no discernible zoning, while the majority of phenocrysts are normally zoned. A few, however, show slight reverse zoning. Plagioclase is abundant and shows normal zoning, with cores about An80 and rims about An55. Magnetite and ilmenite are both present and are relatively homogeneous. There are two apatite populations, one containing very little sulphur and the other exhibiting some zoning in sulphur.

The lava is highly crystalline and shows evidence of along-
flow fractionation. The variation in the tephra textures may
represent discrete batches of magma from depth. There is
localised contamination from lithics of rhyolitic composition,
most probably introduced in the upper regions of the conduit.
Textural and compositional features suggest that some of the
magma was reheated prior to the eruption, while other pheno-
crysts appear to have been transported rapidly from depth. It is
likely that there was at least one recharge event prior to the
2011 eruption, which would be consistent with seismic data.

Clinopyroxene-hosted melt inclusions contain up to ~3, 100 ppm S, whereas matrix glass contains very low sulphur. A separate sulphide phase is present in some samples and as inclusions in phenocrysts. It contains up to 3 % Cu and 1.3 % Ni. Taking the SO₂ mass released into the atmosphere as 4.5 Tg (Theys et al. 2013), a crystallinity of 60 % and the aforementioned value as indicative of melt sulphur content, the dense rock volume of the eruption amounts to approximately 0.47 km³ (magnitude 5.1, given by [log₁₀ (mass in kg)]
<figure><img src="image_22.png" /><figcaption>**Fig. 15** Total alkalis versus silica plot for eruptive products from the 2011 eruption</figcaption></figure>

−7]). This is roughly compatible with the scale of emplaced lava flows (approximately 18 km² in area for an assumed mean thickness of 20 m) and tephra deposits (5 m mean thickness inside the caldera, 0.1 m mean thickness beyond). We hope in the future to be able to obtain further tandem SAR data products to enable a more accurate determination of the lava and tephra volumes. This will be essential for a robust treatment of the volatile budget.

## Impacts of the eruption and emergency response

The Danakil region is generally dry and sparsely populated. However, Nabro, which rises to more than 2,200 m above sea level, engenders a microclimate that favours denser settlement. Several villages are, or were, located within the caldera and on the flanks of the volcano. Siroru, the main settlement within the caldera, had a population of around 3,000 prior to the eruption. The inhabitants of the region are mostly pastoralists keeping herds of sheep, goats, camels, donkeys and cattle. Despite the limited number of recorded earthquakes, the seismic hazard in the region had been previously recognised (e.g. Midzi et al. 1999). The eruption of Dubbi in 1861 (Wiart and Oppenheimer 2000) is also a testament to the volcanic hazard and is remembered by communities living in the area today. Nevertheless, prior to the 2011 eruption, a national risk mitigation strategy for geophysical threats was

lacking. There was no monitoring of any kind on Nabro nor was there a formalised warning system in place.

We believe that the experience of the 4.5 M<sub>w</sub> event on 31 March 2011 may have served to sensitise and increase preparedness of the population in advance of the eruption. Almost all the houses in the affected area that were built with non-mortared stone walls suffered severe damage or collapse; those built with mortared walls suffered varying degrees of cracking (Fig. 5). During the visit of the scientific team in May 2011, village and district administrators and community leaders were briefed on what to do in the event of future earthquakes, namely to vacate their homes and avoid structures that could collapse and areas susceptible to landslides.

Village and district administrators in Eritrea are mostly war veterans well-versed in aspects of emergency management and in coping with situations with limited means at their disposal. Their experience was likely decisive in the immediate local management of the unfolding emergency on 12 June 2011. An army detachment based close to the local administration building in Siroru maintained radio contact with the Southern Red Sea Region Administration (SRSRA) in the regional capital, Assab (about 120 km away, as the crow flies). Together with other nearby units, they helped to evacuate people as the crisis developed on the night of 12/13 June 2011. Clearly, people living within the caldera, on the volcano flanks and in range of the tephra fallout, were most affected by the eruption and its associated seismicity. In
particular, Siroru was destroyed by a combination of seismic shaking, lava flows and burial in tephra. While some people refused to be evacuated, by and large the rescue efforts were conducted with the full cooperation and consent of the affected people.

Another factor that may have contributed to the rapid response is that one of us (JOSH) happened to be in Asmara at the time of the eruption. He received an automated e-mail via the USGS National Earthquake Information Center's 'Preliminary Determination of Epicenters' service concerning the first larger earthquake on 12 June (15:33:12 UTC, i.e. 5 h before the eruption). This information was relayed via another member of our team (GO) to the SRSRA.

The first meeting of the SRSRA committee tasked with coordinating the response to the unfolding humanitarian crisis took place at 08:00 local time on 13 June 2011 in Assab. This resulted in the dispatch of a medical team and provision of emergency relief (food supplies and shelters, Solomon, 2012). Other government agencies provided longer term assistance to displaced people in respect of health care and schooling.

The SRSRA classified the displaced persons into two categories. Those in category 1 came from villages that were destroyed. Those in category 2 were those whose villages were threatened by the eruption including the effects of fumes and tephra fallout. According to Solomon (2012), during the first 2 days of the eruption, the estimated numbers in category 1 were 700–800 families (about 2,500–3,000 individuals). The corresponding values for category 2 were about 2,250 families (about 9,000 individuals). On 17 June 2011, when the eruption intensified, most category 2 evacuees were reclassified as category 1. In total, about 11,780 people were affected. Seven people were killed and three were injured during the eruption (Solomon 2012).

The losses in household goods were substantial (Solomon, 2012). In previously settled areas of the caldera floor south of the vent region, tephra accumulated to depths of several metres, burying or partially burying many structures. Two villages were completely destroyed, and recorded losses include 19,839 goats and sheep, 460 cattle, 834 camels and 142 donkeys. Total economic losses were estimated to be in the region of US$3 million (Solomon 2012). The contamination of water wells in the area, at least on the Eritrean side, appears to have been limited (Solomon 2012).

Many of the people displaced by the eruption were unable to return to their homelands; substantial numbers have been permanently resettled elsewhere. Resettlement can have major consequences for livelihoods, communities and culture and can expose people to new hazards (Usamah and Haynes 2012). Resettlement programmes are arguably most effective when the affected communities are closely involved in decision-making. According to Solomon (2012), this appears to have been the case following the Nabro eruption.

Information on impacts on the Ethiopian side of the border is limited, but the caldera of Mallahle volcano was inhabited and received significant fallout of tephra from the eruption plumes. We were told that the affected communities received assistance in Eritrea (although the two countries do not maintain diplomatic relations). Communications received from staff at a mining operation in Dallol, 170 km to the northwest in Ethiopian Danakil (and below sea level) suggested grounding of a portion of the gas and aerosol emissions, causing discomfort to workers. Respiratory stress may have been exacerbated by the extreme temperatures experienced in this part of northern Afar in the summer.

In terms of civil aviation issues, the region falls under the responsibility of the Toulouse Volcanic Ash Advisory Centre. The volcanic clouds resulted in disruptions to aviation in the region with international flights in East Africa and the Middle East particularly affected. In the first days of the eruption, there was confusion in the outside world as to whether it was Dubbi or Nabro that was erupting.

## Concluding remarks and perspectives

The 2011 Nabro eruption offers a valuable opportunity to develop our understanding of unrest and eruptive activity of caldera systems, the local interactions between tectonics and volcanism and between neighbouring volcanoes and the origins and significance of the off-axis volcanic ranges in the wider Afar region. We summarise key observations and preliminary interpretations as follows:

1. The June eruption was preceded by a damaging earthquake (the 31 March 2011 event) and further felt seismicity and earthquake swarms up to 5 h beforehand. (We heard anecdotally from a local guide that the vigour of steam vents inside the caldera had increased in advance of the eruption. These vents were apparently used by people living within the caldera for therapeutic purposes but had heated up to the point where they were no longer comfortable.) These signs indicate that the eruption was presaged by sensible phenomena.

2. The eruption began between 20:27 and 20:42 UTC, according to satellite images, and possibly close to 20:32 UTC, the time of a large earthquake. It may have been predominantly effusive for the first 50 min or so, then sufficiently explosive to be recorded at distant infrasound stations. A 17.5-km-long lava flow was observed on 16 June 2011, exceeding 1 km in width in places. Explosive activity generated significant tephra clouds and a substantial quantity of SO₂ was released to the atmosphere. Compositions of erupted lavas and tephra range from trachybasalt to basaltic trachyandesite. (At present, the eruption magnitude is poorly constrained, but it could be calculated in future if additional SAR data are acquired for construction of pre- and post-eruption digital topography.)
3. Approximately 12,000 people were displaced as a result of the eruption, and seven fatalities were recorded. Considering that more than 3,000 people lived within the Nabro caldera itself, the toll is remarkably low. This reflects a rapid, largely spontaneous evacuation of settlements close to the vent region, prompted by premonitory ground shaking and the first signs of eruption.

4. Our analysis of regional seismic records for the period 23 February–17 September 2011 revealed 217 events within a 25–40 km radius of Nabro. Seven of these events occurred before the eruption. Their local magnitudes range between 2.1 and 5.9, with 76 events exceeding $M_L$ 4 and eight events larger than $M_L$ 5. While the located events identify temporal variation of seismicity beneath Nabro volcano in the lead up to the eruption, we note that our detection threshold is quite high, about $M_L$ 2.1. Double-difference relocations of the events reveal clustering beneath both Nabro and Mallahle volcanoes and hint at alignments of epicentres, along a NE–SW trend, notably on 17 June 2011 (parallel to the Nabro-Mallahle axis), and, subsequently, along a NW–SE trend (parallel to the Red Sea).

5. Geodetic modelling suggests that the eruption was associated with intrusion of a shallow dike that triggered slip on parallel normal faults. The orientation of this dike is roughly NW–SE, parallel to the trend of an eruptive fissure formed during the eruption and to the Red Sea.

6. In general terms, temporal trends in recorded seismicity, $SO_2$ output and infrasound signals are correlated, with several pulses and fluctuations. The onset of the eruption was marked by high seismicity and infrasound followed by strong $SO_2$ emissions to the atmosphere. Seismic and infrasonic tremor persisted for about a month. A lull occurred on 16 June, but seismicity increased abruptly on 17 June.

Our study of the Nabro eruption has highlighted the importance of diverse Earth observation techniques for monitoring volcanoes in comparatively remote regions, especially where local ground-based sensor networks are limited or lacking. It also exemplifies the precarious status of risk management where it is not underpinned by operational monitoring.

**Acknowledgments** This research was funded by the NERC (project no. NE/J012297/1 ‘Mechanisms and implications of the 2011 eruption of Nabro volcano, Eritrea’) with contributions from JOSH’s NERC fellowship (NE/1020342/1) and supported with access to the NERC Ion Microprobe Facility. Additional funding was received in the form of grant NE/L013932/1 held by DK and the support of several of the authors by the NERC-funded Centre for Observation and Modelling of Earthquakes, Volcanoes and Tectonics (COMET). We acknowledge funding for BG from the EPSRC and the School of Earth Sciences at the University of Bristol; for AD by the Leverhulme Trust and Isaac Newton Trust; for SAC from NASA through the Aura Science Team (grant NNX11AF42G) and MEaSURES (grant NNX13AF50G) programmes; and DF from NSF grant EAR-1113294. TerraSAR-X data were provided by DLR under grant award GEO1206. We acknowledge the cooperation we received

from the Eritrean government, Southern and Northern Red Sea Administrations, local sub-zones and village administrations. We thank the Department of Mines, Ministry of Energy and Mines for their continued support throughout the project. Special thanks go to Zerai Berhe, Mebrahtu Fisseha, Michael Eyob, Ahmed Mohammed, Kibrom Nerayo, Asresehey Ogbatsien (plausible future presenter of *Top Gear*) Andemichael Solomon and Isaac Tuum. We thank Alem Kibreab for his vital help in facilitating the fieldwork. We are very grateful to the two referees (Cindy Ebinger and John Pallister) and the handling editor (Guido Giordamo) for their beneficial and enthusiastic comments on the original manuscript.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## References

Albino F, Sigmundsson F (2014) Stress transfer between magma bodies: influence of intrusions prior to 2010 eruptions at Eyjafjallajökull volcano, Iceland. J Geophys Res Solid Earth 119(4):2964–2975

Barberi F, Tazieff H, Varet J (1972) Volcanism in the Afar depression: its tectonic and magmatic significance. In: Girdler RW (ed) East African rifts, vol 15, Tectonophysics., pp 19–29

Barnie T, Oppenheimer C (2015) Extracting high temperature event radiance from satellite images and correcting for saturation using independent component analysis. Remote Sens Environ 158:56–68

Barnie TD, Keir D, Hamling I, Hofmann B, Belachew M, Carn S, Eastwell D, Hammond JOS, Ayele A, Oppenheimer C. and Wright T (2015b) A multidisciplinary study of the final episode of the Manda Hararo dyke sequence, Ethiopia, and implications for trends in volcanism during the rifting cycle. Geological Society London Special Publications, 420 (doi:10.1144/SP420.6)

Clarisse L, Coheur P-F, Theys N, Hurtmans D, Clerbaux C (2014) The 2011 Nabro eruption, a $SO_2$ plume height analysis using IASI measurements. Atmos Chem Phys 14:3095–3111. doi:10.5194/acp-14-3095-2014

Cote DM, Belachew M, Quillen AC, Ebinger CJ, Keir D, Ayele A, Wright T (2010) Low-frequency hybrid earthquakes near a magma chamber in Afara: quantifying path effects. Bull Seismol Soc Am 100(5A):1892–1903

Doin MP, Guilloso S, Jolivet R, Lasserre C, Lodge F, Ducret G, Grandin R (2011) Presentation of the small baseline NSBAS processing chain on a case example: the Etna deformation monitoring from 2003 to 2010 using Envisat data, In Proceedings of the European Space Agency Symposium “Fringe”. Frascati, Italy

Dziewonski AM, Chou T-A, Woodhouse JH (1981) Determination of earthquake source parameters from waveform data for studies of global and regional seismicity. J Geophys Res 86(2825–2852): 1981. doi:10.1029/JB086iB04p02825

Ekström G, Nettles M, Dziewonski AM (2012) The global CMT project 2004-2010: centroid-moment tensors for 13,017 earthquakes. Phys Earth Planet Inter 200–201(1–9):2012. doi:10.1016/j.pepi.2012.04.002

Endo ET, Murray T (1991) Real-time seismic amplitude measurement (RSAM): a volcano monitoring and prediction tool. Bull Volcanol 53(7):533–545

Fee D, Matoza RS, Gee KL, Neilsen TB, Ogden DE (2013) Infrasonic crackle and supersonic jet noise from the eruption of Nabro Volcano, Eritrea. Geophys Res Lett 40:4199–4203. doi:10.1002/grl.50827
Field L, Barnie T, Blundy J, Brooker RA, Keir D, Lewi E, Saunders K (2012) Integrated field, satellite and petrological observations of the November 2010 eruption of Erta Ale. Bull Volcanol 74(10):2251–2271

Fromm M, Kablick G III, Nedoluha G, Carboni E, Grainger R, Campbell J, Lewis J (2014) Correcting the record of volcanic stratospheric aerosol impact: Nabro and Sarychev Peak. J Geophys Res Atmos 119:10,343–10,364. doi:10.1002/2014JD021507

Grandin R et al (2009) September 2005 Manda Hararo-Dabbahu rifting event, Afar (Ethiopia): constraints provided by geodetic data. J Geophys Res 114:B08404. doi:10.1029/2008JB005843

Hamiel Y, Grandin R, Baer G, Wright TJ, Doubre C, Hamlyn J, Peltzer G (2013) InSAR investigation of crustal deformation associated with the 2011 eruption of the Nabro volcano, Eritrea. In AGU Fall Meeting Abstracts (Vol. 1, p. 0955)

Hamlyn J, Keir D, Wright TJ, Neuberg J, Goitom B, Hammond JO, Pagli C, Oppenheimer C, Kendall J-M, Grandin R (2014) Seismicity and subsidence following the 2011 Nabro eruption, Eritrea: insights into the plumbing system of an off-rift volcano. J Geophys Res Solid Earth 119:8267–8282. doi:10.1002/2014JB011395

Hammond JOS, Kendall J-M, Stauart GW, Ebibger CJ, Bastow ID, Keir D, Ayele A, Belachew M, Goitom B, Ogubazghi G, Wright TJ (2013) Mantle upwelling and initiation of rift segmentation beneath the Afar Depression. Geology 41:635–638

Jónsson S, & Xu W (2015). Volcanic eruptions in the southern Red Sea during 2007–2013. In *The Red Sea* (pp. 175–186). Springer Berlin Heidelberg

Kanamori H (1977) The energy release in great earthquakes. J Geophys Res 82:2981–2987. doi:10.1029/JB082i020p02981

Keir D, Stuart GW, Jackson A, Ayele A (2006) Local earthquake magnitude scale and seismicity rate for the Ethiopian rift. Bull Seismol Soc Am 96(6):2221–2230

Keir D, Hamling IJ, Ayele A, Calais E, Ebinger C, Wright TJ, Jacques E, Mohamed K, Hammond JOS, Belachew M, Baker E, Rowland JV, Lewi E, Bennati L (2009) Evidence for focused magmatic accretion at segment centers from lateral dike injections captured beneath the Red Sea rift in Afar. Geology 37(1):59–62. doi:10.1130/G25147A.1

Keir D, Pagli C, Bastow IaI.n D, Ayele A (2011) The magma-assisted removal of Arabia in Afar: evidence from dike injection in the Ethiopian rift captured using InSAR and seismicity. Tectonics, 30, (2), TC2008. (doi:10.1029/2010TC002785)

Keir D, Bastow ID, Pagli C, Chambers EL (2013) The development of extension and magmatism in the Red Sea rift of Afar. Tectonophysics 607:98–114. doi:10.1016/j.tecto.2012.10.015

Klein, F.W., (2002). User's guide to HYPOLINVERSE-2000, a Fortran program to solve for earthquake locations and magnitudes. Open File Report 02-171, US Geological Survey

Makris J, Ginzburg A (1987) The Afar depression—transition between continental rifting and sea-floor spreading. Tectonophysics 141(1–3):199–214

McClusky S et al (2010) Kinematics of the southern Red Sea–Afar Triple Junction and implications for plate dynamics. Geophys Res Lett 37: L05301. doi:10.1029/2009GL041127

Midzi V, Hlatywayo DJ, Chapola LS, Kebede F, Atakan K, Lombe DK, Turyomurugyendo G, Tugume FA (1999) Seismic hazard assessment in eastern and Southern Africa. Ann Geofis 42:1067–1083

Mogi K (1958) Relations between the eruptions of various volcanoes and the deformations of the ground surfaces around them. Bull Earthq Res Inst 36:99–134

Nobile, A., Pagli, C., Keir, D., Wright, T. J., Ayele, A., Ruch, J., and Acocella, V. (2012) Dike-fault interaction during the 2004 Dallol intrusion at the northern edge of the Erta Ale Ridge (Afar, Ethiopia). Geophysical Research Letters, 39, (19), L19305. (doi: 10.1029/2012GL053152).

Ogubazghi G, Goitom B, Kibreab A, King J (2014) Relocation of the 31 March 2011 earthquake in the Nabro Volcanic Range of Southern Eritrea. Eritrean J Sci Eng 1:93–102

Okada Y (1985) Surface deformation due to shear and tensile faults in a half-space. Bull Seismol Soc Am 75(4):1135–1154

Rosen PA, Hensley S, Peltzer G, Simons M (2004) Updated repeat orbit interferometry package released. EOS Trans Am Geophys Union 85(5):47–47

Rowland J, Baker E, Ebinger CJ, Keir D, Kidane T, Biggs J, Wright TJ (2007) Fault growth at a nascent slow-spreading ridge: 2005 Dabbahu rifting episode, Afar. Geophys J Int 171(3):1226–1246

Rubin AM, Pollard DD (1988) Dike-induced faulting in rift zones of Iceland and Afar. Geology 16(5):413–417

Santer, B. D., et al. (2015) Observed multivariable signals of late 20th and early 21st century volcanic activity. Geophys. Res. Lett., 42, 500–509, doi:10.1002/2014GL062366 DOI:10.1002/2014GL062366#_blank#Link to external resource: 10.1002/2014GL062366

Scordilis EM (2006) Empirical global relations converting Ms and mb to moment magnitude. J Seismol 10:225–236. doi:10.1007/s10950-006-9012-4

Sealing, C. R. (2013) Characterizing the first historic eruption of Nabro, Eritrea: insights from thermal and UV remote sensing. MS thesis, Michigan Technological University

Solomon A (2012) Summary of activities of the Southern Red Sea Administration during the Nabro eruption. Department of Infrastructure of the Southern Red Sea Administration, Assab

Tarantola A, Valette B (1982) Generalized nonlinear inverse problems solved using the least squares criterion. Rev Geophys 20(2):219–232

Theys N, Campion R, Clarisse L, Brenot H, van Gent J, Dils B, Corradini S, Merucci L, Coheur P-F, Van Roozendael M, Hurtmans D, Clerbaux C, Tait S, Ferrucci F (2013) Volcanic SO₂ fluxes derived from satellite data: a survey using OMI, GOME-2, IASI and MODIS. Atmos Chem Phys 13:5945–5968. doi:10.5194/acp-13-5945-2013

Usamah M, Haynes K (2012) An examination of the resettlement program at Mayon Volcano: what can we learn for sustainable volcanic risk reduction? Bull Volcanol 74(4):839–859

Waldhauser F, Ellsworth WL (2000) A double difference earthquake location algorithm: method and application to the Northern Hayward Fault, California. Bull Seismol Soc Am 90(6):1353–1368

Walter TR, Shirzaei M, Manconi A, Solaro G, Pepe A, Manzo M, Sansosti E (2014) Possible coupling of Campi Flegrei and Vesuvius as revealed by InSAR time series, correlation analysis and time dependent modeling. J Volcanol Geotherm Res 280:104–110

Wiart PAM, Oppenheimer C (2000) Largest known historic eruption in Africa: Dubbi volcano, Eritrea, 1861. Geology 28:291–294

Wiart P, Oppenheimer C (2005) Large magnitude silicic volcanism in north Afar: the Nabro Volcanic Range and Ma'alalta volcano. Bull Volcanol 67:99–115

Wiart PAM, Oppenheimer C, Francis P (2000) Eruptive history of Dubbi volcano, northeast Afar (Eritrea), revealed by optical and SAR image interpretation. Int J Remote Sens 21:911–936

Xu W, Ruch J, Jonsson S (2015) Birth of two volcanic islands in the southern Red Sea. Nat Commun 6:7104. doi:10.1038/ncomms8104

Yang K, Krotkov NA, Krueger AJ, Cam SA, Bhartia PK, Levelt PF (2007) Retrieval of large volcanic SO2 columns from the Aura Ozone Monitoring Instrument: comparison and limitations. *Journal of Geophysical Research: Atmospheres* (1984–2012), 112(D24)

Yirgu G, Ferguson DJ, Barnie TD, Oppenheimer C (2014) Recent volcanic eruptions in the Afar rift, northeastern Africa, and implications for volcanic risk management in the region. In: Ismail-Zadeh A, Fucugauchi JU, Kijko A, Takeuchi K, Zaliapin I (eds) Extreme natural hazards, disaster risks and societal implications. Cambridge University Press, Cambridge, p 200–213