Strength variation along the Altyn Tagh and the Kunlun fault, northern Tibetan plateau, inferred from 3D mechanical modeling

The Altyn Tagh (ATF) and the Kunlun (KLF) fault distribute around the northern Tibetan plateau from west to east about 2000 km and 1200 km in length, and deform predominately with left-lateral strike-slip motion. Previous geological and geodetic observations suggested that over at least 800-km length of the faults, the slip rate averaged on active deformation period is quite uniform, for the ATF being about 9-10 mm/yr and the KLF about 10-12mm/yr. Strike-slip deformation of these faults is undoubtedly result from regional loading by ongoing collision between the India and the Eurasia continent. Whereas, dense GPS measurements show that along the central Tibetan plateau from west to east, the GPS velocity field changes greatly both on magnitude and on direction, suggesting that tectonic loading to the ATF and the KLF could be changed along their strike directions. To investigate how a non-uniform tectonic loading condition as documented by the GPS velocity field could cause a relatively uniform slip rate of the two active faults, we built a three-dimensional viscoelastic finite element model, in which motion of the strike-slip fault is governed by frictional strength. Given a reasonable bound of model parameters, we at first test the numerical calculation with uniform frictional coefficient of the faults. At this condition, the predicted slip rate is inevitably largest near center of the faults and gradually decreasing to the fault ends. To better fitting the observed uniform slip rate along the faults over 1000km length, variation of fault strength along the ATF and the KLF must be invoked. By testing numerous models, an optimum result was obtained, among which the frictional coefficient for the ATF is varied from 0.02 to 0.12 between 820E and 1000E with its maximum at 840E, and for the KLF from 0.02 to 0.10 with its maximum between 950E and 970E. This means that the strength of the two large-scale strike-slip faults exists significant difference along their strikes. We believe that the predicted fault pattern could play an important role on partitioning strain aside the fault, together on determination of potential rupture during an earthquake.