Coulomb Stress Perturbations Related to the Al Hoceima ( Morocco ) Earthquakes of 1994 and 2004

This paper exposes the results of the study of the stress perturbations caused by the Al Hoceima (Morocco) earthquakes of 1994 and 2004 by means of Coulomb modeling. Modeling was based on the compilation of all the studies carried out after both events, including seismological (location and depth of the main shocks and aftershocks), seismotectonic (source parameters, stress field), geodetic (GPS), tomographic and geological ones. It shows that the first earthquake is likely to have induced the second one when adopting appropriate epicenter locations, source and receiver fault planes. In detail, the model shows that motion along the N23E oriented Bousekkour-Aghbal fault in 1994 activated the eastern segment of a previously unknown NW-SE fault located at its southern end, which in turn originated the 2004 earthquake. The model also provides a suitable explanation for the distribution of the aftershock clusters. Finally, an attempt of prediction of the next event shows that it is likely to occur on NE-SW planes located to the NW (Bokkoya and offshore) and SE (reaching the Nekor fault) of the NW-SE fault, while Al Hoceima city should remain in a shadow zone.


Introduction
During the last decades, several research teams have investigated the possible transfer of stress after the occurrence of large/moderate earthquakes to other areas (faults), which may trigger similar events after a relatively short interval, e.g.[1] [2].The amount of stress transfer is commonly studied using Coulomb stress change (he-F.Medina reafter abbreviated to CSC) modeling; e.g.Coulomb 3 [3], which is a quantitative and graphical method for determining the change in stress within the region located around the main shock [4].CSC modeling has been applied with success in studies carried out on several seismic zones such as California [4] [5], Japan [6] and Turkey [7], and therefore it appears as an efficient tool for assessing the seismic hazard in a given area.
In northwest Africa, Morocco is located at the limit of the Nubian plate boundary along the Azores-Gibraltar Fault Zone (AGFZ) [8] (Figure 1) and is threatened by moderate and large earthquakes [9].The largest earthquakes that have been recorded instrumentally since 1934, date of installation of the first seismological station at Averroes Observatory [10], were the "twin" shocks of Al Hoceima on 26 May 1994 (Mw = 5.8 -6.0), and 24 February 2004 (Mw = 6.3 -6.5).The first caused two deaths and important material damage [11], while the second caused 628 fatalities and damaged 2539 constructions [12].
The occurrence of closely related large earthquakes in space and time arises two main questions: 1) did the first shock of 26 May 1994 trigger that of 24 February 2004, as already suggested by [13], who used the term "conjugate faults" for the almost-perpendicular rupture planes?And 2) if so, where should the next earthquake occur?In order to assess this issue, CSC modeling was carried out on the Al Hoceima area, as for the "twin" earthquakes (M = 5) which hit Rissani (south-eastern Morocco) on 23 and 30 October 1992 [14].After exposing the general tectonic and seismic characteristics of the area (Section 2), the available data for setting the main constraints for modeling CSC are carefully examined (Section 3).The possible relationship of the 1994 and 2004 earthquakes are then investigated (Section 4) using Coulomb stress change software, and the implications on the seismic hazard of the region are discussed (Section 5).

Regional Geodynamic Setting
Seismicity of the Ibero-Maghrebian area is related to the plate motion of the Nubian and Eurasian plates, which is divergent in the Azores region and right-lateral to convergent eastward e.g.[8] and references therein [9].Depending on the authors and the used methodology (magnetic anomalies and/or GPS data), the Nubia-Eurasia Eulerian pole (Eurasia fixed) is located between 18.2˚S and 21.6˚N and 15.9˚W and 30.8˚W [15] [16].In contrast to the more or less linear Azores-Gibraltar segment of the Nubia-Eurasia plate boundary (Figure 1), the Rif-Betics-Alboran area is a complex mosaic of micro-blocks undergoing variable motion.The Rif and Betics are symmetrical thrust belts undergoing compressional and/or strike-slip regime, while the western Alboran Sea is an extensional (back-arc) basin [17].The Al Hoceima area is located at the southern segment of a major NNE-SSW fault zone spanning the eastern Betics e.g.[18], the Alboran Sea [19] and the Central Rif [20].This shear Figure 1.Main structural features along the Azores-Gibraltar Fault Zone and tectonic and seismotectonic setting of Morocco.Bathymetry from NOAA (http://maps.ngdc.noaa.gov/viewers/bathymetry/). Atlantic structures from [8] [88].GCAW = Gulf of Cadiz accretionary wedge; TASZ = Trans-Alboran Shear Zone.Mean Africa (Nubia)-Eurasia plate motion adapted from [88] and references therein.zone, called the Trans-Alboran Shear Zone (TASZ) [21], may represent the eastern termination of the continental western Alboran basin and the eastward transition to the oceanic part (Figure 1).However, the tectonic framework is much more complex in the Alboran basin [22].The description of the detailed geological history of this area e.g.[17] is beyond the scope of our paper, but knowledge of its present-day kinematics is necessary for understanding its seismic activity and assessing the seismic hazard.
Analyses of GPS data indicate that the direction of convergence is WNW-ESE (Figure 1) with a rate of 4.3 -5.6 mm•yr -1 [9] and references in their Figure 5 [23]- [25].The regional model of [24] considers five blocks bounded by fault zones: 1) Iberia and the Atlantic segment north of the AGFZ; 2) Morocco and the Atlantic segment south of the AGFZ; 3) the Betics block; 4) the central Rif and western Alboran block, and 5) the eastern Alboran block.The faults bounding the Betics block are right-lateral, with the exception of the eastern fault, the TASZ, which is left-lateral.Southwards, the presumed faults bounding the central Rif block are left-lateral (the TASZ and the Rif front fault) with the exception of the western fault zone, which is right-lateral.Therefore, with respect to Nubia, a south-westward motion is observed in the central Rif and its foreland [20] [23] [24] [26] in relation to a SW-verging tectonic escape process observed from the surface [27] which may extend at depth [24].The boundary between both kinematic zones is located along the TASZ, the continental segment of which is located west of Al Hoceima [24].The amount of slip along this boundary is evaluated at 5.2 to 5.7 mm•yr -1 [24].Focal mechanism solutions and analysis of the subsequent state of stress in northern Morocco also indicate a NW-SE convergence [9] [28]- [36] in accordance with the kinematic models.

Local Geological Setting and Main Faults
The city of Al Hoceima is located in the central Rif chain, which runs parallel to the northern coast of Morocco (Figure 2).The most conspicuous structures consist of southward-verging thrust sheets which are, from the upper (northernmost) to the lower (southernmost): 1) the Bokkoya (Palaeozoic and Mesozoic-Cenozoic); 2) the Tiziren unit (Middle Jurassic to Early Cretaceous carbonates and flysch series); and 3) the Ketama metamorphic unit (Cretaceous flysch and limestones).Volcanic rocks are represented by the middle Miocene Ras-Tarf andesites [37] [38].Unconformable Miocene and younger deposits [39] are encountered in the lower Nekor graben, trending N-S, and in the Boudinar basin, east of Ras Tarf Cape [37] [40].Offshore, the main structures are the Plio-Quaternary Bokkoya basin, which extends until the Alboran ridge [41] and the lower Nekor basin which is the marine continuation of the onshore one [22] [41] [42].The NE-SW left-lateral Nekor fault is the major structure in the area e.g.[43]; however, it has remained seismically inactive since the onset of seismic recordings in Morocco [36] [44] [45].Other important faults are (Figure 2): 1) The Imzouren (NNW-SSE) and Trougout (N-S) normal faults, which respectively delimit the western and eastern boundaries of the lower Nekor graben; new studies show that the Trougout fault has an oblique (normal/strike-slip) slip rate of 0.73 -0.95 mm•yr −1 since 7 ka [46], and can be followed offshore for about 7 km with a surface offset [42]; authors [46] calculate a return period for major earthquakes of about 3 ka, which is a relatively long period; 2) The Jbel Hammam fault system (NNW-SSE), consisting of several normal faults with trace lengths of 20 km; 3) The Rouadi fault, a normal fault trending NNE-SSW over 15 km, the eastern block being downthrown; 4) The Bousekkour-Arhbal strike-slip fault (Figure 2), trending NNE-SSW, which crosscuts the whole Bokkoya unit with a sinistral displacement [47].The latter fault is thought to have been active during the 1994 seismic crisis [11], but radar interferometry and field studies suggest that the rupture did not reach the surface [48].However, a recent oceanographic survey [42] has shown that the fault can be followed offshore for about 4 km where it is expressed by an east-dipping normal fault plane (50 m offset of "basement") with a clear 5 m-surface fault scarp.Other offshore faults are those delimiting and affecting the Nekor graben [41], and especially the NNE-SSW striking, 4 km long, Bokkoya normal fault whose throw is about 50 m [42].

Deep Structure
Recent multidisciplinary studies e.g.[49] [50] show that the thickness of the crust in northern Morocco decreases northwards from 30 km in the Rif to 12 km in the central Alboran basin, whereas the depth of the lithosphere-asthenosphere boundary in the area decreases from 230 -250 km in the west, to about 70 km in the Alboran Sea, the isobaths running roughly N-S.Crustal thinning is accompanied by the development of a large thermal anomaly in the Alboran basin [51].Around Al Hoceima, the deep structure was mainly studied using seismic tomography [36] [52]- [54].The results obtained by these authors appear strongly related to the network configuration, and especially to the distance between recording stations; for instance, P-wave velocity maps obtained by [52] show in the layer 5 -15 km a large circular (~100 km diameter), high velocity zone south-east of Al Hoceima, which extends at depth in the layer 15 -30 km, contiguous to a high velocity zone.The boundary of both zones runs NNW-SSE across the lower Nekor graben.Authors [53] inferred low velocity zones within the layer 5 km beneath the Nekor graben and within the lower crust in the continental area, whereas the most important high velocity zone is located south of Al Hoceima city within the layer 10 km.Tomographic maps of [54] show a low-velocity zone at 15 km depth near Al Hoceima; however, the resolution is not appropriate for determining better orientations.Instead, the higher resolution tomography obtained by [36] shows better constrained velocity zones at 5 -15 km depth, with a determinable orientation.

Seismicity
Historically, the Al Hoceima area suffered numerous strong earthquakes as shown by the earliest earthquake catalogues [44] [55].The best documented events were recorded in 1522, 1624, 1790-1791, 1795 and 1800-1803 (M~6?), 1848, 1910-1914.Since the onset of the instrumental period, the early and following authors e.g.[44] [56]- [58] observed that seismicity is largely scattered on the scale of the Betic-Rif Arc.On the base of the data from NEIC, the epicenters of events with magnitudes > 3 from 1 st January 1960 to 3 October 2014 in the Al Hoceima area are shown in Figure 3.Most of the earthquakes are shallow, and their magnitudes are generally moderate, except for the 1994 and 2004 events.The epicenters are clearly aligned along a NNE-SSW trend parallel to the TASZ; however, on the regional scale, epicenters of some well-constrained micro-earthquakes and aftershocks of large events are clearly aligned along NNE-SSW, NNW-SSE and WNW-ESE trends [11] [36] [41] [59]- [61].The seismic flux (the seismic moment released per year and per unit area) calculated from a 104 year period is one of the largest in the Western Mediterranean, reaching more than 10 19 dyne•cm•km -2 •yr -1 [9].

Focal Mechanisms and Ground Surface Motion
Since the 1970's and until 2007, at least 100 focal mechanisms were determined for large and moderate earthquakes (M > 3.5) of the Al Hoceima area using first motion polarities of P-waves, waveform analysis and moment tensor inversion [62].Among these, a large number of solutions were determined by seismological agencies (IGN, IAG, MED, INGV) for the aftershocks of the 2004 earthquake.Most fault-plane solutions, including the largest aftershocks of the two major earthquakes, correspond to strike-slip faulting with a normal component and to normal faulting (Figure 4).A few solutions correspond to reverse faulting or to strike-slip motions with reverse component.The P axes are distributed in the NW and SE quadrants with variable plunge.The solutions obtained for interseismic periods [28] [45] [63] [64], also indicate that the dominant regime is strike-slip/normal, with a general lack of pure reverse mechanisms.The P axes are distributed along a NW-SE great circle, with largely variable plunges, whereas T axes are NE-SW to ESE-WNW with a shallow plunge [28] [65].
From GPS data listed by [23] [24], the general displacements with respect to fixed Africa calculated from the East and North components, are variable: east of the TASZ, permanent station MELI shows a SSE (N173) small motion (1.61 mm•yr −1 ) during the period 1999-2007, while nearby temporary station MDAR shows a NNE (N20) moderate motion (2.24 mm•yr −1 ) during the period 1999-2004.West of the trace of the TASZ, the motion at temporary stations BBFH (period 1999-2004) and KTMA (period 1999-2005) is SSW (N203 and N229 respectively), and the linear velocity is significant (3.31 and 3.7 mm•yr −1 respectively).

Principle
Coulomb stress change is best expressed by the equation e.g.[4]: where c ∆σ , ∆τ and N ∆σ are, respectively, the coseismic changes in Coulomb stress, shear stress and nor- mal stress, and ' µ is the effective coefficient of friction, e.g. [66]and references therein, which corresponds to ( ) where B is Skempton's coefficient, comprised between 0.5 and 0.9.Modeling was carried out on the region of coordinates 34.8˚N -35.7˚N × 3.3˚W -4.8˚W with the help of Coulomb 3.2 software [3].In a first step, we considered the fault related to the 1994 earthquake as the source fault and that of the 2004 shock as the receiver fault.In a second step, we considered the effects of the 2004 earthquake on the nearby region.CSC was resolved both on optimal strike-slip planes and on chosen particular planes parallel to the considered receiver faults.A discussion on this point can be found in [67].The following subsections expose in detail the model constraints which are necessary as input data to the construction of the model, such as the regional state of stress and the source parameters including the location of the epicenter, the nucleation depth, length and width of the rupture and aftershocks location.

Model Constraints
The model constraints are based on a large database collected after the 1994 and 2004 earthquakes by numerous international research teams that carried out studies on the source parameters of both earthquakes, using various methods such as classical seismic and micro-seismic surveys [11]

Regional State of Stress
The state of stress in the Al Hoceima and surrounding region, compiled in Table 1, was determined by several authors from background microseismicity [28], from the aftershocks of the seismic crises of 1994 [11] and 2004 [36] [69] and from general seismicity [28] [31]- [35].The determinations are mainly based on focal mechanism solutions using the right dihedra method [74] [75], and/or numerical stress inversion methods [76]- [83].The state of stress obtained for three sets of moderate to large earthquakes (isolated, 1994 and 2004) show similar orientation of axes, reflecting a strike-slip regime, with σ 1 and σ 3 horizontal and trending NNW-SSE (N140˚E to N156˚E) and ENE-WSW (N42˚E to N65˚E) respectively.For micro-earthquakes (shocks of magnitude 1 to 3), the state of stress is extensional in the Imzouren area and strike-slip elsewhere [28].As discussed in detail previously [11] [28] [32], large to moderate earthquakes are related to the regional state of stress driven by plate and micro-plate motions, while micro-earthquakes only reflect a local state of stress related to block adjustments.

The 26 May 1994 Earthquake
Characteristics and ground effects.The main shock of 26 May 1994, occurred at 8 h 27 min UTC.Various epicenter locations were proposed for the main shock (Figure 5), but all authors agree that it was located west of Al Hoceima, probably near-or offshore.The magnitude was Md = 5.6 (Mw = 6.0), and the hypocentral depth was 13 km.The values of the seismic moment Mo vary from 1.1 × 10 17 N•m (waveform analysis [70]) to 2.1 × 10 18 N•m (InSAR; [72]).As reported by [11], the ground effects observed after the main shock consisted of decametric-scale N40 to N50 en échelon cracks, and landslides along the coastal cliffs west of Al Hoceima and more inland.No unequivocal surface faults were observed during the survey.The isoseismal map elaborated by [11] based on EMS 1992 scale shows a maximum damage area (intensity VIII-IX EMS) trending NNE-SSW from 35.25˚N × 4˚W to 35.12˚N × 4.15˚W.
Focal mechanisms.At least 9 fault-plane solutions were determined by various seismological centers (USGS, CSEM and Harvard) and authors [11] [70] [72] for the main shock of 1994.Solutions obtained by moment tensor inversion (USGS, Harvard, CSEM, [84]) or waveform modeling [70] [72] show strike-slip motion with a slight normal component.The solution proposed by [11] on the base of first motions includes a reverse component.Reference [70] suggested that the main event consisted of two sub-events, the mechanism of sub-event 2 also showing a reverse component.Seven solutions determined for the aftershocks correspond to strike-slip/ normal faulting, except for the one of 3 June at 8 h 57 min, which shows a reverse component as for the main shock [11].The T-axes are almost-horizontal and are oriented NE-SW to ENE-WSW; the P-axes have a moderate plunge (60˚ maximum).
Location of fault rupture.In addition to the source parameter determinations from seismic data [70], relatively concordant radar interferometry analyses performed by [13] [72] suggest (Table 2) that the ground deformation can be best modeled by a N23 striking, 10 -16 km long, 10 -15 km wide fault dipping east, with a 0.6 -0.8 m left-lateral slip.The rupture area remains beneath the surface (depth ≥ 2 km).
Aftershock sequence.The DPG (Institut Scientifique) temporary network of 7 short-period (1 Hz) analogical stations installed by [11] from 27 May to 9 June 9 1994 recorded 512 events located onshore, whose highest magnitude was Md = 4.4.The 68 best constrained events correspond to 2 -18 km deep, NNE-SSW trending cluster of aftershocks, which is largely distributed over a 30 × 10 km region that matches the maximum damage area (Figure 6).The aftershocks appear southward of the epicenter of the main shock when the latter is taken to be offshore.These authors state that the shift may be due either to errors on the determination of the main shock epicenter, or to sensitivity threshold of the seismic stations, which would not record weak events located offshore.Among the vertical sections oriented at different angles to the NNE-SSW trending cluster, the one oriented WNW-ESE, perpendicular to the cluster, suggests that the foci are concentrated along a plane with a steep dip to the east-southeast [11] their figure 8.The aftershock sequence of the 1994 earthquake was also studied by [41], based on teleseisms recorded by the permanent digital stations of the CNRST (Morocco), and  the PDE of the USGS epicenters recorded from January 1994 to June 1995.The aftershock cluster appears to the east of that shown by [11].Authors [41] concluded that the observed shift to the west of PDE locations with respect to CNRST ones may be due to "an area of anomalous velocities, coupled with a relative dearth of seismic stations on the African continent for good azimuthal control".Authors [72] also studied the aftershock sequence of the 1994 shock.From ISC and PDE (USGS) files, only 9 well-constrained events were relocated; among which 4 are located along the NNE-SSW fault trace, whereas the other 5 are located either to the NW, to the SW and to the SSE.

The 24 February 2004 Earthquake
Characteristics and ground effects.The main shock of 24 February 2004 occurred at 2 h 27 min UTC.The epicenter of the earthquake was located southeast of Al Hoceima at Ait Kamra near Imzouren (35.28˚N × 3.99˚W).Its magnitude Md was 6.3 (Mw = 6.2), and the hypocentral depth was 6 km.The seismic moment Mo was estimated by [60] at 1.8 × 10 18 N m.From analysis of body wave forms at teleseismic distances, [60] found a complex rupture process along a NNE-SSW fault with bilateral rupture formed by 4 shallow sub-events, with rupture propagating towards the north.However, no numerical data were published to date.Authors [33] studied the rupture process using regional ASTM.Their results involve 2 sub-events separated by 3 sec, located on parallel faults striking N11, the second sub-event being located southwest of the first.The ground effects observed after the main shock consist of variable size (centimetric to kilometric) cracks oriented ENE-WSW to NNE-SSW, along a 20-km-wide corridor south of Al Hoceima, and metric to decametric landslides [48] [61] [68].The isoseismal map based on MKS scale shows a maximum damage area (intensity VIII-IX) trending NE-SW from 35.20˚N × 3.75˚W to 35.05˚N × 4.10˚W [68].
Focal mechanisms.Focal mechanisms determined for the main shock from moment tensor inversion (USGS, Harvard, LDG) show nodal planes oriented NNE-SSW and WNW-ESE with ESE and SSW steep dips respectively.Mechanisms of the strongest aftershocks, determined by IGN (Spain) correspond to strike-slip faulting with a reverse component mainly, and to reverse faulting (Figure 4).Mechanisms of the weakest aftershocks also correspond to strike-slip faulting with mainly a normal component [36] [69].P-axes concentrate around azimuth N330˚E.
Location of fault rupture.In order to find the location of the fault trace for this event, radar interferometry analyses were performed by several teams [13] [61] [71]- [73].Their results are summarized in Table 2. Authors [13] [71] indicate that the ground deformations are best modeled by a NW-SE curved-shape, 19 × 14 km fault dipping 87˚NE, with a maximum 2.7 m right-lateral slip without reaching the surface.Authors [72] also favor a NW-SE, 10 × 18 km fault dipping 87.4 NE, with a 1.4 right lateral slip.Authors [61] suggested a double rupture along two sub-vertical faults oriented N10 (9 × 11.5 km) and N312 (15 × 9 km), with displacements of 0.92 m (left-lateral) and 0.76 m (right-lateral) respectively.However, it is somewhat surprising that the determined NW-SE fault trace does not match the NW-SE alignment of numerous aftershock epicenters, but appears shifted to the south.Finally, [73] arrive to similar fault parameters using MTInSAR.
Aftershock sequence.Among the several aftershock surveys conducted using teleseismic, regional permanent and/or portable seismographs e.g.[12] [72], the most detailed survey carried out from 28 March to 10 April [36] [61] [69] clearly shows that the 650 aftershocks are 1 -16 km depth and aligned along NNE-SSW and NW-SE trends.The NNE-SSW cluster is located in the southward continuation of the suggested fault, whereas the NW-SE clusters are located parallel to the NW-SE plane.Other alignments appear to the south-east.Early aftershocks located with the help of the permanent seismological stations are aligned NW-SE but are also located to the NE of the proposed plane.
Surface ground motion.GPS data from [61] indicate that, taking into account the annual displacement with respect to fixed Africa (see Section 2.5), station BBFH, located 30 km to the west of Al Hoceima, underwent a 30 mm displacement toward the west, while station MDAR, located 40 km to the southeast of Al Hoceima, only showed a 5 mm displacement.Finally, comparison of two SPOT5 images shows that the co-seismic rupture could not have exceeded 5 -10 cm at the surface.

The 1994 Source Fault (Fault 1)
Epicenter location and focal depth.Based on the field survey of [11], we chose the epicenter located around 35.26˚N × 4.0˚W for the 1994 earthquake, which is situated on the chosen fault plane.The adopted depth is 7.5 km.
Fault parameters.For the 1994 shock (the source fault), we calculated CSC simplifying the parameters exposed by [72]: strike N23˚; dip 85˚(E); length 8 km; width 8 km (2 -10 km depth); mean displacement 0.7 m left-lateral (when uniform), which lead to a seismic scalar moment of about 1.4 × 10 25 dyn•cm (1.4 × 10 18 N•m) and a moment magnitude of ~6.0.In order to avoid in excess concentration of stress at the tips of the fault, the plane was tapered into 5 patches as recommended by the Coulomb 3 tutorial.
Stress orientation and magnitude.CSC calculations were carried out taking into account the regional state of stress, derived from our R4DT determination for the Al Hoceima area [28] [32]: σ 1 and σ 3 horizontal at 150˚E and N60˚E respectively.The stress magnitudes used by software Coulomb are 100 bar for σ 1 at the surface; 0 bar for σ 3 and 30 bar for σ 2 .No vertical gradient was taken into account.The CSC calculation depth is 7.5 km.
Fault strength parameters.We used the default values recommended by the Coulomb 3 software.These are 8 × 10 5 bar for Young modulus; 0.25 for Poisson's ratio, and 0.8 for the coefficient of friction.

The 2004 Source Fault (Fault 2)
Epicenter location and focal depth.For the 2004 shock, a source with coordinates 35.13˚N × 3.95˚W at a depth of 7.5 km was used.
Fault parameters.Fault parameters (Table 2) are also those calculated by [72]: strike N295˚; dip 87˚N; length 10 km; width 10 km (2 -12 km depth) and 1 m of right-lateral displacement, which lead to a seismic scalar moment of 3.2 × 10 25 dyn.cm and a moment magnitude Mw = 6.3.The fault trace was constrained by recent aftershock data published by [36] repository material.
Stress orientation and magnitudes and fault strength parameters are the same than those used in the previous case.

CSC after the 1994 Earthquake
Figure 6(A) illustrates the CSC distribution for optimal strike-slip planes after 0.7 m left-lateral slip on fault 1.It clearly shows an increase (≥5 bar, saturation value) in the magnitude of the stresses at the tips of the fault, and a decrease (≤-5 bar, saturation value) across the plane along a NE-SW oriented zone.Most of the epicenters of the aftershocks taken from [11] (for clarity, only the 80 first ones recorded from 27 to 30 May are shown) plot in the positive CSC zone oriented at a small anti-clockwise angle from the fault, but a large number plot in a negative CSC zone.The receiver fault undergoes positive CSC mainly at its centre.
In the case of CSC resolved on the N20E (=N200E) vertical planes oriented parallel to the main fault (Figure 6(B)), all aftershock epicenters fall within the large negative CSC area situated on both fault blocks.However, a small positive CSC zone affects the receiver fault at its eastern tip.
Finally, when CSC is resolved for specified N300E (N120E) striking dextral strike-slip planes (Figure 6(C)), parallel to the receiver fault, almost-all aftershocks fall in a negative CSC zone.The receiver fault also undergoes positive CSC all along its eastern segment.The aftershock epicenters determined by [36] (only the first 300 ones from 29 March to 6 April were plotted for clarity) fall within a Coulomb stress increase area (Figure 7(A)).However, the cluster at 35.15˚N × 4.05˚W falls in a shadow zone, which does not match the Coulomb stress change.No changes are observed when the CSC is resolved on planes which are parallel to the main fault (Figure 7(B)), and the cluster at 35.15˚N × 4.05˚W remains in the same CSC decrease zone.

CSC after the 2004 Earthquake
When CSC is resolved on sinistral planes striking N20E, no important changes are observed with respect to the previous panel (Figure 7C).Finally, CSC resolved on normal fault planes such as found by [36] at the eastern aftershock cluster (Figure 7(D)) shows that the aftershocks also fall in the increase zone, which is in conformity with aftershock triggering.

Robustness of CSC Modeling in Al Hoceima Area
As indicated in section 1, CSC modeling has already been successfully used for studying local and regional stress changes after major earthquakes, as well as for predicting the loci of occurrence of future ones in several seismic areas in the world.In the present study, CSC modeling was applied to the Al Hoceima region, which is an area of distributed moderate seismicity with relatively short faults.
From our results, CSC modeling was successful in helping understand several observations: 1) CSC after the 1994 earthquake may have led to a stress increase along the eastern segment of the WNW-ESE fault responsible for the 2004 earthquake; therefore, as for the Rissani twin earthquakes of 1992 [14], the model provides a consistent explanation for the space relationship between both events.
2) Aftershock distribution and reactivation of optimally oriented planes after both earthquakes are also successfully accounted for by CSC modeling.
However, although the geometrical parameters seem to be relatively well constrained on the base of the data taken from the numerous studies of the Al Hoceima earthquakes, there are intrinsic parameters which remain less constrained.This is the case of the variable input values suggested by the software such as the friction coefficient, the regional stress magnitude and its variability with depth, which are beyond the scope of the present article, focused on regional aspects.

Seismogenic Faults in the Al Hoceima Area
One of the major problems for constraining the relationship between earthquakes and faults in the Central Rif is that a few ruptures were observed during both events [11] [36] [48] [68], although rupture has been proved to reach the surface offshore [42], so the faults remained blind and their pattern was only delineated afterwards using focal mechanism solutions, aftershock distribution and radar interferometry.Indeed, the WNW-ESE fault responsible for the 2004 earthquake had no previous morphological or deep features which could have helped mapping it.
The post-seismic field survey of the 2004 earthquake [48] suggests that there seems to be some vertical partitioning of deformation, with thrust sheets at the surface and strike-slip faults at depth.However, this is not true everywhere, because it is obvious that some faults reach the surface, such as the Trougout-Bou Haddoud, Rouadi, Bousekkour-Aghbal, Bokkoya offshore and Boudinar normal faults which have been and still are potentially seismogenic [42] [46].
Finally, at depth, tomographic maps obtained by [54] show a low-velocity zone at 15 km depth near Al Ho-ceima; however, the resolution of the maps is not appropriate for determining accurate orientations.Instead, the higher resolution tomography obtained by [36] clearly shows that the aftershock hypocentres are located within high velocity zones at 5 -15 km depth, whose orientation is in conformity with the chosen faults, which may not be an artifact of the used methodology.

Seismic Hazard
In order to predict which fault planes are most loaded by positive CSC and may be reactivated, we plotted the main known faults of the Al Hoceima area and run the program with Coulomb stress resolved onto planes having the same strike than the major faults recognized in the area.

Figures 7(D)
and Figure 8 show CSC after the 2004 earthquake resolved onto the strikes of the regional fault planes located around the epicenter.The figures show that the faults around Al Hoceima city are not affected by an increase in stress, with the exception of the Nekor fault, the south-western segment of which is affected by an increase of about 7 bar.However, it should be emphasized that most stress was probably released through the very numerous aftershocks.
The time interval for the occurrence of the next event is for the moment difficult to determine because it depends on the tectonic load related to the plate motion of Nubia with respect to Eurasia, but also on the regional motion of blocks across the TASZ and the accommodation of motion by faults.Estimations should be possible when a larger and more precise GPS database will be available in the next years.
Anyway, CSC modeling is a tool which is, at a local scale (10,000 km 2 ), an efficient alternative method to probabilistic and deterministic approaches of seismic hazard in the area [86] [87].

Conclusions
Coulomb modeling of the stress perturbations caused by the Al Hoceima earthquakes of 1994 and 2004 based on the compilation of all the studies carried out after both events, including seismological (location and depth of the main shocks and aftershocks), seismotectonic (source parameters, stress field), geotedic (GPS), tomographic and geological ones, leads to the conclusions that: 1) The 1994 earthquake is likely to have induced the 2004 one when adopting the most accurate epicenter locations, source fault traces and receiver fault planes.Motion along the N23E oriented Bousekkour-Aghbal fault in 1994 activated the eastern segment of an unknown NW-SE fault located at its southern end, which originated the 2004 earthquake.
2) The model also provides a suitable explanation for the distribution of the aftershocks' clusters.
3) Finally, prediction of the next event shows that it is likely to occur on NE-SW planes located to the NW (Bokkoya and offshore) and SE (reaching the Nekor fault) of the NW-SE fault, while Al Hoceima city will remain in a shadow zone.

Figure 3 .
Figure 3. Epicenters of earthquakes with magnitudes > 3 recorded by NEIC (USGS) from 1 st January 1960 to 3 October 2014 in the Al Hoceima area.

Figure 4 .
Figure 4. Focal mechanisms of the Al Hoceima area from 1959 to 2006 after [65].

Figure 5 .
Figure 5. Location of the epicenter of the May 1994 shock according to different seismological agencies and fault traces used for Coulomb stress modeling.

Figure 7 (
Figure7(A)shows the Coulomb stress distribution after 1 m right-lateral slip on Fault 2 resolved on optimal strike-slip planes.An increase (≥5 bar, saturation value) in the magnitude of the stresses is observed at the tips of the fault, whereas a decrease (≤-5 bar, saturation value) can be seen across the plane along a NE-SW direction.The aftershock epicenters determined by[36] (only the first 300 ones from 29 March to 6 April were plotted for clarity) fall within a Coulomb stress increase area (Figure7(A)).However, the cluster at 35.15˚N × 4.05˚W falls in a shadow zone, which does not match the Coulomb stress change.No changes are observed when the CSC is resolved on planes which are parallel to the main fault (Figure7(B)), and the cluster at 35.15˚N × 4.05˚W remains in the same CSC decrease zone.When CSC is resolved on sinistral planes striking N20E, no important changes are observed with respect to the previous panel (Figure7C).Finally, CSC resolved on normal fault planes such as found by[36] at the

Figure 6 .
Figure 6.Results of CSC modeling of the 1994 earthquake: (A) On optimally oriented strike-slip planes; (B) On planes parallel to the main fault; (C) On planes oriented N300 (or N120).

Figure 7 .
Figure 7. Results of CSC modeling of the 2004 earthquake: (A) On optimally oriented strike-slip planes; (B) On planes parallel to the main fault; (C) On planes oriented N20; (D) On NNW-SSE oriented normal planes.

Figure 8 .
Figure 8. Results of CSC modeling of the 2004 earthquake on planes parallel to the Nekor fault.

Table 1 .
State of stress determined by various authors for the Al Hoceima area.

Table 2 .
Fault parameters determined by several authors for the Al Hoceima earthquakes of 1994 and 2004.Data in bold characters were those taken into account as input parameters for the CSC model.
* curved fault trace, ** determined graphically by the present author from figures.