3D mechanical modeling of the GPS velocity field along the North Anatolian fault

The North Anatolian fault (NAF) extends over 1500 km in a complex tectonic setting. In this region of the eastern Mediterranean, collision of the Arabian, African and Eurasian plates resulted in creation of mountain ranges (i.e. Zagros, Caucasus) and the westward extrusion of the Anatolian block. In this study we investigate the effects of crustal rheology on the long-term displacement rate along the NAF. Heat flow and geodetic data are used to constrain our mechanical model, built with the three-dimensional finite element code ADELI. The fault motion occurs on a material discontinuity of the model which is controlled by a Coulomb-type friction. The rheology of the lithosphere is composed of a frictional upper crust and a viscoelastic lower crust. The lithosphere is supported by a hydrostatic pressure at its base (representing the asthenospheric mantle). We model the long-term deformation of the surroundings of the NAF by adjusting the effective fault friction and also the geometry of the surface fault trace. To do so, we used a frictional range of 0.0^0.2 for the fault, and a viscosity varying between 10 19 and 10 21 Pa s. One of the most striking results of our rheological tests is that the upper part of the fault is locked if the friction exceeds 0.2. By comparing our results with geodetic measurements [McClusky et al., J. Geophys. Res. B 105 (2000) 5695^5719] and tectonic observations, we have defined a realistic model in which the displacement rate on the NAF reaches V17 mm/yr for a viscosity of 10 19 Pa s and a fault friction of 0.05. This strongly suggests that the NAF is a weak fault like the San Andreas fault in California. Adding topography with its corresponding crustal root does not induce gravity flow of Anatolia. Rather, it has the counter-intuitive effect of decreasing the westward Anatolian escape. We find a poor agreement between our calculated velocity field and what is observed with GPS in the Marmara and the Aegean regions. We suspect that the simple lithosphere model is responsible for this discrepancy. Taking into account the weaknesses of these deforming regions should allow us to build a more realistic model that would match ground observations more appropriately. On the other hand, our results fit well GPS measurements in central Anatolia, setting the basis of modeling crustal strain in Turkey. B 2003 Elsevier Science B.V. All rights reserved.

[1]  P. Bird.  Testing hypotheses on plate‐driving mechanisms with global lithosphere models including topography, thermal structure, and faults , 1998 .

[2]  R. Sabadini,et al.  The block‐like behavior of Anatolia envisaged in the modeled and geodetic strain rates , 2002 .

[3]  Michael Ortiz,et al.  Finite element analysis of strain localization in frictional materials , 1989 .

[4]  Rob Westaway,et al.  Present‐day kinematics of the Middle East and eastern Mediterranean , 1994 .

[5]  M. Toksöz,et al.  Source mechanism of the Bartin earthquake of September 3, 1968 in northwestern Turkey: Evidence for active thrust faulting at the southern Black Sea margin , 1986 .

[6]  Present-day crustal movements and the mechanics of cyclic deformation , 1990 .

[7]  Bertrand Meyer,et al.  Slip partitioning in the Sea of Marmara pull‐apart determined from GPS velocity vectors , 2002 .

[8]  J. Braun.  Three-dimensional numerical simulations of crustal-scale wrenching using a non-linear failure criterion , 1994 .

[9]  M. Toksöz,et al.  Tectonic consequences of the collision of the Arabian and Eurasian plates: Finite element models , 1983 .

[10]  Bertrand Meyer,et al.  Westward propagation of the North Anatolian fault into the northern Aegean: Timing and kinematics , 1999 .

[11]  N. N. Ambraseys,et al.  The Seismic Activity of the Marmara Sea Region over the Last 2000 Years , 2002 .

[12]  Ross D. Hartleb,et al.  The Surface Rupture and Slip Distribution of the 17 August 1999 İzmit Earthquake (M 7.4), North Anatolian Fault , 2002 .

[13]  P. Meijer,et al.  Temporal variation in the stress field of the Aegean Region , 1996 .

[14]  Sean D. Willett,et al.  Mechanical model for the tectonics of doubly vergent compressional orogens , 1993 .

[15]  L. Royden,et al.  Crustal rheology and faulting at strike‐slip plate boundaries: 2. Effects of lower crustal flow , 2000 .

[16]  Demitris Paradissis,et al.  Global Positioning System constraints on plate kinematics and dynamics in the eastern Mediterranean and Caucasus , 2000 .

[17]  O. M. Ilkişik.  REGIONAL HEAT-FLOW IN WESTERN ANATOLIA USING SILICA TEMPERATURE ESTIMATES FROM THERMAL SPRINGS , 1995 .

[18]  J. Dewey.  Extensional collapse of orogens , 1988 .

[19]  James Jackson,et al.  Active tectonics of the Alpine—Himalayan Belt between western Turkey and Pakistan , 1984 .

[20]  S. Kirby,et al.  Rheology of the lithosphere: Selected topics , 1987 .

[21]  M. Zoback,et al.  New Evidence on the State of Stress of the San Andreas Fault System , 1987, Science.

[22]  B. Hager,et al.  Viscoelastic deformation from North Anatolian Fault Zone earthquakes and the eastern Mediterranean GPS velocity field , 2002 .

[23]  D. McKenzie.  Active Tectonics of the Mediterranean Region , 1972 .

[24]  P Cundall,et al.  A MICROCOMPUTER PROGRAM FOR MODELLING LARGE-STRAIN PLASTICITY PROBLEMS. PROCEEDINGS OF THE SIXTH INTERNATIONAL CONFERENCE ON NUMERICAL METHODS IN GEOMECHANICS, 11-15 APRIL 1988, INNSBRUCK, AUSTRIA. VOLUMES 1 - 3 , 1988 .

[25]  Wayne Thatcher,et al.  Microplate versus continuum descriptions of active tectonic deformation , 1995 .

[26]  Bertrand Meyer,et al.  Asymmetric slip partitioning in the Sea of Marmara pull‐apart: a clue to propagation processes of the North Anatolian Fault? , 2001 .

[27]  Bertrand Meyer,et al.  The active Main Marmara Fault , 2001 .

[28]  L. Rybach,et al.  Geothermal reconnaissance of the Marmara Sea region (NW Turkey): surface heat flow density in an area of active continental extension , 1998 .

[29]  A. Barka,et al.  Slip distribution along the North Anatolian fault associated with the large earthquakes of the period 1939 to 1967 , 1996, Bulletin of the Seismological Society of America.

[30]  Haluk Ozener,et al.  Estimates of Seismic Potential in the Marmara Sea Region from Block Models of Secular Deformation Constrained by Global Positioning System Measurements , 2002 .

[31]  H. Houston,et al.  Orientation of the stress field surrounding the creeping section of the San Andreas Fault: Evidence for a narrow mechanically weak fault zone , 2001 .

[32]  A. Zanchi,et al.  Seismotectonics of western Anatolia: Regional stress orientation from geophysical and geological data , 1993 .

[33]  M. Boccaletti,et al.  Reproducing the velocity and stress fields in the aegean Region , 1997 .

[34]  W. Lettis,et al.  Influence of Releasing Step-Overs on Surface Fault Rupture and Fault Segmentation: Examples from the 17 August 1999 İzmit Earthquake on the North Anatolian Fault, Turkey , 2002 .

[35]  D. Albarello,et al.  Short and long term deformation patterns in the Aegean‐Anatolian Systems: Insights from space geodetic data (GPS) , 2001 .

[36]  Larry G. Mastin,et al.  Regional patterns of tectonic stress in Europe , 1992 .

[37]  A. Şengör,et al.  The North Anatolian transform fault: its age, offset and tectonic significance , 1979, Journal of the Geological Society.

[38]  H. Kahle,et al.  GPS and geologic estimates of the tectonic activity in the Marmara Sea region, NW Anatolia , 1997 .

[39]  J. Chéry,et al.  An integrated mechanical model of the San Andreas Fault in central and northern California , 2001 .

[40]  J. Chéry,et al.  A Physical Basis for Time Clustering of Large Earthquakes , 2001 .

[41]  A. Deschamps,et al.  Active deformation of the Corinth rift, Greece: Results from repeated Global Positioning System surveys between 1990 and 1995 , 2000 .

[42]  M. Gupta,et al.  Terrestrial heat flow and geothermal energy in Asia , 1995 .

[43]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[44]  D. Giardini,et al.  A geodynamic framework for eastern Mediterranean kinematics , 1998 .

[45]  J. Suppe,et al.  State of stress near the San Andreas fault: Implications for wrench tectonics , 1987 .

[46]  C. Froidevaux,et al.  Tectonics and topography for a lithosphere containing density heterogeneities , 1982 .

[47]  M. Richards,et al.  Finite element models of stress orientations in well‐developed strike‐slip fault zones: Implications for the distribution of lower crustal strain , 2001 .

[48]  J. Chéry,et al.  Study of plate deformation and stress in subduction processes using two‐dimensional numerical models , 1997 .

[49]  J. Baumgardner,et al.  Fault friction, regional stress, and crust‐mantle coupling in southern California from finite element models , 1984 .

[50]  P. Huchon,et al.  The deformation of aegea: A finite element study , 1982 .

[51]  New Interpretation as a Paleo,et al.  The North Anatolian Fault Zone: New Interpretation as a Paleo-Benioff Zone , 1975 .

[52]  E. Boschi,et al.  Numerical modelling of the Aegean–Anatolian region: geodynamical constraints from observed rheological heterogeneities , 2001 .