Fault identification for buried strike-slip earthquakes using InSAR: The 1994 and 2004 Al Hoceima, Morocco earthquakes

SUMMARY The 1994 M w 6.0 and 2004 M w 6.5 Al Hoceima earthquakes are the largest to have occurred in Morocco for 100 yr, and give valuable insight into the poorly understood tectonics of the area. Bodywave modelling indicates the earthquakes occurred on near-vertical, strike-slip faults with the nodal planes oriented NW‐SE and NE‐SW. Distinguishing between the primary fault plane and auxiliary planes, using either geodetic or seismic data, is difficult due to the spatial symmetry in deformation fields and radiation pattern of moderately sized, buried, strikeslip earthquakes. Preliminary studies, using aftershock locations and surface observations, have been unable to identify the orientation of the primary fault plane for either earthquake conclusively. We use radar interferometry and aftershock relocation of the earthquake sequence to resolve the ambiguity. For the 2004 earthquake, inverting the interferograms for a uniform slip model based either of the two potential nodal planes results in similar misfits to the data. However, the NE‐SW best-fit fault plane has an unrealistically high fault slip-to-length ratio and we conclude the NW‐SE striking nodal plane is the primary fault plane and slip was right lateral. We carry out tests on synthetic data for a buried strike-slip earthquake in which the orientation of the fault plane is known a priori. Independent of geometry, missing data, and correlated noise, models produced assuming the auxiliary plane to be the fault plane have very high fault slip-to-length ratios. The 1994 earthquake had a smaller magnitude and comparisons of model misfits and slip-to-length ratios do not conclusively indicate which of the nodal planes is the primary fault plane. Nonetheless, the InSAR data provides valuable information by improving the accuracy of the earthquake location by an order of magnitude. We carry out a multiple event relocation of the entire earthquake sequence, including aftershocks, making use of the absolute locations for the 1994 and 2004 main shocks from our InSAR study. The aftershock locations are consistent with a NW‐SE orientated fault plane in 2004 and suggests that the 1994 earthquake occurred on a NE‐SW fault; perpendicular to the fault which ruptured in 2004. Previous tectonic models of the area proposed a bookshelf model of block rotation with NNE‐SSW left-lateral faults. This model requires modification to accommodate the observation of right-lateral slip on a NW‐SE fault plane for the 2004 earthquake and we prefer to interpret the fault orientations as due to a zone of distributed shear with a right-lateral fault striking at ∼115 ◦ and conjugate, clockwise rotating, left-lateral faults striking at ∼25 ◦ .

[1]  J. Morales,et al.  Source analysis of the Mw 6.3 2004 Al Hoceima earthquake (Morocco) using regional apparent source time functions , 2005 .

[2]  T. Wright,et al.  The 1998 Aiquile, Bolivia earthquake: A seismically active fault revealed with InSAR , 2005 .

[3]  Rowena B. Lohman,et al.  Some thoughts on the use of InSAR data to constrain models of surface deformation: Noise structure and data downsampling , 2005 .

[4]  Rongjiang Wang,et al.  Implications of the 2003 Bingöl Earthquake for the Interaction between the North and East Anatolian Faults , 2004 .

[5]  Ziyadin Cakir,et al.  Coastal uplift and thrust faulting associated with the Mw = 6.8 Zemmouri (Algeria) earthquake of 21 May, 2003 , 2004 .

[6]  J. Jackson,et al.  Active faulting and seismicity of the Dasht-e-Bayaz region, eastern Iran , 2004 .

[7]  P. Rosen,et al.  Updated repeat orbit interferometry package released , 2004 .

[8]  J. Nocquet,et al.  Evidence for a post-3.16-Ma change in Nubia Eurasia North America plate motions ? , 2003 .

[9]  Simon McClusky,et al.  GPS constraints on Africa (Nubia) and Arabia plate motions , 2003 .

[10]  Zhong Lu,et al.  Source model for the Mw 6.7, 23 October 2002, Nenana Mountain Earthquake (Alaska) from InSAR , 2003 .

[11]  Rowena B. Lohman,et al.  Location and mechanism of the Little Skull Mountain earthquake as constrained by satellite radar interferometry and seismic waveform modeling , 2002 .

[12]  L. Rivera,et al.  Coseismic Deformation from the 1999 Mw 7.1 Hector Mine, California, Earthquake as Inferred from InSAR and GPS Observations , 2002 .

[13]  H. Zebker,et al.  Fault Slip Distribution of the 1999 Mw 7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements , 2002 .

[14]  S Das,et al.  Simultaneous Rupture Along Two Conjugate Planes of the Wharton Basin Earthquake , 2001, Science.

[15]  K. Priestley,et al.  A re‐assessment of focal depth distributions in southern Iran, the Tien Shan and northern India: do earthquakes really occur in the continental mantle? , 2000 .

[16]  R. Stein The role of stress transfer in earthquake occurrence , 1999, Nature.

[17]  Peter J. Clarke,et al.  Source parameters of the 1 October 1995 Dinar (Turkey) earthquake from SAR interferometry and seismic bodywave modelling , 1999 .

[18]  Mourad Bezzeghoud,et al.  Source parameters of the 1992 Melilla (Spain, MW = 4.8), 1994 Alhoceima (Morocco, MW = 5.8), and 1994 Mascara (Algeria, MW = 5.7) earthquakes and seismotectonic implications , 1999, Bulletin of the Seismological Society of America.

[19]  James Jackson,et al.  The 1997 May 10 Zirkuh (Qa'enat) earthquake (Mw 7.2): faulting along the Sistan suture zone of eastern Iran , 1999 .

[20]  M. Ramdani,et al.  The Al Hoceima earthquake of May 26, 1994 and its aftershocks: a seismotectonic study , 1998 .

[21]  C. Werner,et al.  Radar interferogram filtering for geophysical applications , 1998 .

[22]  A. D. Miller,et al.  Non‐double‐couple earthquakes 1. Theory , 1998 .

[23]  E. Engdahl,et al.  Global teleseismic earthquake relocation with improved travel times and procedures for depth determination , 1998, Bulletin of the Seismological Society of America.

[24]  R. Bro,et al.  A fast non‐negativity‐constrained least squares algorithm , 1997 .

[25]  F. Gomez,et al.  An integrated geophysical investigation of recent seismicity in the Al-Hoceima region of north Morocco , 1997, Bulletin of the Seismological Society of America.

[26]  E. R. Engdahl,et al.  Constraints on seismic velocities in the Earth from traveltimes , 1995 .

[27]  Hiromichi Tsuji,et al.  Coseismic crustal deformation from the 1994 Hokkaido‐Toho‐Oki Earthquake Monitored by a nationwide continuous GPS array in Japan , 1995 .

[28]  D. Hatzfeld,et al.  Microearthquake seismicity and fault plane solutions around the Nékor strike-slip fault, Morocco , 1993 .

[29]  P. Segall,et al.  Comparison of various inversion techniques as applied to the determination of a geophysical deformation model for the 1983 Borah Peak earthquake , 1992, Bulletin of the Seismological Society of America.

[30]  Masaya Matsuura,et al.  Geodetic data inversion using a Bayesian information criterion for spatial distribution of fault slip , 1992 .

[31]  P. Molnar,et al.  Fault plane solutions of earthquakes and active tectonics of the Tibetan Plateau and its margins , 1989 .

[32]  K. Hudnut,et al.  Cross-fault triggering in the November 1987 Superstition Hills earthquake sequence , 1989 .

[33]  C. Werner,et al.  Satellite radar interferometry: Two-dimensional phase unwrapping , 1988 .

[34]  Robert McCaffrey,et al.  SYN3: A Program for Inversion of Teleseismic Body Wave Forms on Microcomputers , 1988 .

[35]  C. Morley Origin of a major cross-element zone: Moroccan Rif , 1987 .

[36]  Stuart A. Spikin Estimation of earthquake source parameters by the inversion of waveform data: Global seismicity, 1981-1983 , 1986 .

[37]  S. Sipkin Interpretation of non‐double‐couple earthquake mechanisms derived from moment tensor inversion , 1986 .

[38]  Y. Okada Surface deformation due to shear and tensile faults in a half-space , 1985 .

[39]  P. Olivier,et al.  Role of strike-slip faults in the Betic-Rifian orogeny , 1984 .

[40]  S. Sipkin,et al.  Estimation of earthquake source parameters by the inversion of waveform data: synthetic waveforms , 1982 .

[41]  J. Jackson,et al.  Relations between surface deformation, fault geometry, seismicity, and rupture characteristics during the El Asnam (Algeria) earthquake of 10 October 1980 , 1981 .

[42]  Christopher H. Scholz,et al.  Off-fault aftershock clusters caused by shear stress increase? , 1981 .

[43]  Thomas H. Jordan,et al.  Teleseismic location techniques and their application to earthquake clusters in the South-Central Pacific , 1981 .

[44]  John H. Woodhouse,et al.  Determination of earthquake source parameters from waveform data for studies of global and regional seismicity , 1981 .

[45]  J. M. Robb,et al.  Evolution of the continental margin of southern Spain and the Alboran Sea , 1980 .

[46]  G. Backus,et al.  Uniqueness in the inversion of inaccurate gross Earth data , 1970, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[47]  T. J. Farrugia,et al.  The Shuttle Radar Topography Mission-Data Validation and Applications , 2006 .

[48]  James Jackson,et al.  The 1994 Sefidabeh (eastern Iran) earthquakes revisited: new evidence from satellite radar interferometry and carbonate dating about the growth of an active fold above a blind thrust fault , 2006 .

[49]  Akira Nishitani,et al.  Geodetic data inversion using ABIC to estimate slip history during one earthquake cycle with viscoelastic slip-response functions , 2004 .

[50]  Mustapha Meghraoui,et al.  Tectonique plio-quaternaire de la chaine tello-rifaine et de la mer d'Alboran; une zone complexe de convergence continent-continent , 1996 .

[51]  P. Chotin,et al.  Failles actives dans le Rif central et oriental (Maroc) , 1990 .

[52]  Ian Parsons,et al.  Surface deformation due to shear and tensile faults in a half-space , 1986 .

[53]  NábÄ lek,et al.  Determination of earthquake source parameters from inversion of body waves , 1984 .

[54]  W. Menke Geophysical data analysis : discrete inverse theory , 1984 .

[55]  Victor E. Camp,et al.  The Sistan suture zone of eastern Iran , 1983 .