The Al Hoceima Mw 6.4 earthquake of 24 February 2004 and its aftershocks sequence
Graphical abstract
Introduction
The region of Al Hoceima, Morocco, was shaken by a strong Mw 6.4 earthquake (CSEM, USGS, Harvard) on February 24th, 2004 at 2h27mn (UTC) (Jabour et al., 2004). It is the largest earthquake since at least two centuries (Fig. 1; El Mrabet, 2005a, El Mrabet, 2005b, Pelaez et al., 2007) in the Al Hoceima region considered as the most seismically active of Morocco (Cherkaoui et al., 1990, El Alami et al., 1998, Bezzeghoud and Buforn, 1999). An earthquake of magnitude Mw 6.0 occurred on 2nd May 1994 in the same area about 10 km to the west-northwest of the 2004 event in the morning of a spring day, when most of the local people were in the fields and only few casualties were to deplore. The 2004 earthquake occurred in the night, in the winter season, and caused more than 800 fatalities. While the magnitude of the 2004 earthquake was larger, this difference alone does not explain the difference in human loss. Traditional houses made of stones and adobe and heavy roofs suffered considerable destructions while the modern building in Al Hoceima experienced weak damages. However, recent buildings suffered heavy damages, as far as complete collapse in Imzouren, some 12 km south-east of Al Hoceima (Fig. 2) resulting mainly from bad design (Murphy Corella, 2005) though without significant site effect (Goula et al., 2005). An accelerometer installed close to a dam on the Nekor river, about 10 km southeast from the epicentre recorded a maximum horizontal acceleration of 0.24 g, a relatively high value (CNRST, 2004).
Despite several geophysical studies since the earthquake (Stich et al., 2005, Cakir et al., 2006, Akoglu et al., 2006, Biggs et al., 2006, Tahayt et al., 2009), the source of the 2004 event remains debated and puzzling (Galindo-Zaldivar et al., 2009). Indeed, similarly to the 1994 event, the 2004 event did not rupture to the surface and occurred along a fault that has no cumulative surface expression, to the contrary of well expressed normal faults in the Nekor basin (Trougout and Ajdir faults) and to the west (Rouadi fault), or such as the Nekor strike-slip fault (Morel, 1988, Van der Woerd et al., 2005) (Fig. 2, Fig. 3). Surface deformations (open cracks, ground ruptures, landslides) were observed mainly along a 10 km-wide strip running from Ajdir to Beni Abdellah, in a roughly NNE-SSW direction (Fig. 2, Fig. 4; Ait Brahim et al., 2004), but modelling of radar interferograms (InSAR) suggest a WNW-ESE fault plane (Cakir et al., 2006, Akoglu et al., 2006, Tahayt et al., 2009, Biggs et al., 2006). In addition, the aftershock clouds extend along two almost perpendicular directions similarly to the directions of the nodal planes of the main shock focal mechanism (this study; Harvard catalogue; USGS) with the main shock at the intersection of the two clouds (Fig. 2b) making any attempt to link the seismicity to a preferential fault plane difficult.
The Rif region of northern Morocco belongs to the complex plate boundary between western Africa and western Europe, which extends for more than 1500 km E-W from Portugal to Tunisia and 1000 km N-S from France to south Morocco, and comprises several collisional belts, the Pyrenees, the Betics, the Rif, the Middle Atlas, the Atlas, and the Tellian Atlas (Fig. 1). North of the Rif, the Alboran Sea, characterized by a thinned continental crust, resulted from extension in a back arc setting during the Oligo-Miocene (Watts et al., 1993, Seber et al., 1996). The knowledge of the present day tectonics in the Rif has strong implications in the understanding of strain pattern in the western Mediterranean region (Ait Brahim and Chotin, 1984, Ait Brahim et al., 1990, Morel and Meghraoui, 1996, Lopez Casado et al., 2001). While the convergence between western Africa and western Europe can be assessed by plate tectonic models or GPS velocity fields (DeMets et al., 1990, DeMets et al., 1994, Nocquet and Calais, 2003, Nocquet, 2012), the manner and where the strain is taken into account remains debated, mostly because the rates of deformation are slow and need to be averaged over long time intervals.
Both earthquakes of 1994 and 2004 are characterized by blind ruptures that have not propagated to the surface making their seismotectonic characterizations challenging. The challenge is re-enforced as these two earthquakes are the largest known in this part of the Rif since several centuries (El Mrabet, 2005a, Pelaez et al., 2007). It has been proposed that the ruptures remained located below a structural level made of shallow dipping décollement, remnants of the Miocene Rif nappes (Galindo-Zaldivar et al., 2009, Negro et al., 2007, Mattauer, 1964). However, the nearby occurrence of active faults shaping the landscape of the lower Nekor basin as well as the well-known Nekor and Jebbha strike-slip-faults remains to be explained (Leblanc and Olivier, 1984).
To better document the deformation of the Al Hoceima region, we present here an aftershocks study of the 2004 event together with a review of recent seismological and geological data.
The main shock of February 24, 2004 was followed by a large number of aftershocks, some of them with magnitude greater than 5 between February 25th and 28th, and on March 20th (Table S3). The main shock and the aftershocks were located independently by the national seismological network of Spain and Morocco and by the European-Mediterranean Seismological Center. The analysis of the aftershocks we present here results from 15 days of recording with a temporary network installed in the epicentral area from March 29 to April 11 2004 (Dorbath et al., 2005; Table 1). The main shock and the largest aftershocks were relocated. The P and SH waveforms recorded at teleseismic distances were then inverted. We then compare these results with the main tectonic features of the region.
Section snippets
Tectonic summary
The Al Hoceima region is part of the Rif mountain belt, which, with the Betics, form the west vergent Gibraltar arc of lower Miocene age (Andrieux et al., 1971, Tapponnier, 1977, Gutscher et al., 2002, Frizon de Lamotte et al., 1991; Fig. 1). Several models have been proposed to explain the simultaneous uplift of the Betics and Rif ranges and the subsidence of the Alboran Sea in between during the Miocene implying various subduction geometries and roll-back mechanisms (Platt and Vissers, 1989,
Data and processing
The post-seismic field study resulted from a common initiative between Centre National de la Recherche Scientifique et Technique (CNRST, Rabat), Ecole et Observatoire des Sciences de la Terre (EOST, Strasbourg) and Institut de Physique du Globe de Paris (IPG, Paris). Unfortunately, the fieldwork began quite late, more than 4 weeks after the occurrence of the main shock. In a first step, 17 stations were installed. After few days, 3 of them were moved to improve the coverage of the aftershocks
Relocation of the main shock and early aftershocks
To gain more precision on the aftershock sequence we relocated the main shock and aftershocks that occurred before the installation of our temporary network through a master event technique. We used stations of the Spanish national network (IGN) and the national network of Morocco (CNRST). The arrival times of P and S waves were taken from the respective bulletins. We selected a dozen of Spanish stations to have the best azimuthal coverage as possible to the north and to the east, where the
Modelling of the main shock from P and SH waves, teleseismic inversion
Teleseismic P and SH waveforms were inverted to obtain an independent determination of the main shock focal mechanism, seismic moment, source depth, as well as an estimation of source dimension. Broadband seismograms used here were recorded by the IRIS and GEOSCOPE networks and data processing includes deconvolution from the instrument response, integration to obtain displacement, equalization to a common magnification and epicentral distance, and bandpass filtering from 0.01 Hz to 0.8 Hz (P
Seismic tomography
The travel-times of the best-located aftershocks have been used for a tomographic inversion using tomoDD code (Zhang and Thurber, 2003). From the 676 selected events, we obtained about 11,400 absolute travel-times (∼38% being S waves), and constructed more than 50,100 differential travel-times with inter event distances of less than 10 km.
We tested several inversion grids and found stable results, i.e., a velocity model independent of the sampling of the volume, with a grid spacing of 4 km
Complexity of the seismicity
The source fault of the 2004 Al Hoceima earthquake remains uncertain despite many efforts in seismological, geodetic and geologic data analysis. We present the first results from a local seismological aftershock survey (Dorbath et al., 2005). From March 29 to April 11, 2004, our seismological field survey, about 1 month after the event, recorded a set of more than 600 aftershocks in the epicentral area (Fig. 2, Fig. 6). At first order, the aftershocks are distributed along two directions with
Conclusions
The studies on the late aftershocks recorded by the temporary seismological network should be summarized in the following points:
- (1)
The distribution of the aftershocks does not allow deciphering which of the two branches displayed is the actual fault where the main shock occurred. Some observations favour the N20E direction: (a) The aftershock cluster is densest in the N10-15E direction and is about 20 km-long. (b) This direction is in agreement with the left-lateral slipping plane from the main
Acknowledgements
We are grateful for support from Ministère de l’Intérieur, France, via the “Service des Missions et Conférences à l’étranger” EGIDE for travel and field expenses during the field survey in 2004. The field survey has been made possible thanks to the support and welcome of the Wilaya of Al Hoceima, the Wali and Gouverneur of Al Hoceima, M. Mohamed Afoud, the Wali of Taza, M. El Hassan Benameur and the Wali of Taounate, M. Mohamed Amghouz. We thank the Général de Brigade Abdelkrim El Yakoubi, de
References (74)
- et al.
The 1994–2004 Al Hoceima (Morocco) earthquake sequence: conjugate fault ruptures deduced from InSAR
Earth Planet. Sci. Lett.
(2006) - et al.
Sur un modèle explicatif de l’arc de Gibraltar
Earth Planet. Sci. Lett.
(1971) - et al.
Seismotectonics of the Ibero-Maghrebian region
Tectonophysics
(1995) - et al.
The 2012 Sumatra great earthquake sequence
Earth Planet. Sci. Lett.
(2012) - et al.
Are the seismological and geological observations of the Al Hoceima (Morocco, Rif) 2004 earthquake (M = 6.3) contradictory?
Tectonophysics
(2009) - et al.
Microearthquake seismicity and fault plane solutions around the Nékor strike-slip fault, Morocco
Earth Planet. Sci. Lett.
(1993) - et al.
New GPS constraints on active deformation along the Africa–Iberia plate boundary
Earth Planet. Sci. Lett.
(2011) - et al.
Role of strike-slip faults in the Betic-Rifian orogeny
Tectonophysics
(1984) Notes on the Neogene basin history of the western Alboran Sea and its implications for the tectonic evolution of the Rif-Betic orogenic belt
J. Afr. Earth Sci.
(1992)Present-day kinematics of the Mediterranean: a comprehensive overview of GPS results
Tectonophysics
(2012)