articolo wasburn.pdf Key words paleoseismology – Altyn Tagh Fault – strike-slip faults – India-Eurasia collision 1. Introduction The sinistral Altyn Tagh Fault (ATF) travers- es the northern boundary of Tibet for > 1500 km and is a major structure in the India-Eurasia con- Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China Zack Washburn (1), J Ramón Arrowsmith (2), Guillaume Dupont-Nivet (3), Wang Xiao Feng (4), Zhang Yu Qiao (4) and Chen Zhengle (4) (1) Geomatrix Consultants Inc., Costa Mesa, CA, U.S.A. (2) Department of Geological Sciences, Arizona State University, Tempe, AZ, U.S.A. (3) Paleomagnetic Laboratory «Fort Hoofdijk», Faculty of Earth Sciences, Utrecht University, Netherlands (4) Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing, China Abstract Although the Altyn Tagh Fault (ATF) is thought to play a key role in accommodating India-Eurasian conver- gence, little is known about its earthquake history. Studies of this strike-slip fault are important for interpreta- tion of the role of faulting versus distributed deformation in the accommodation of the India- Eurasia collision. In addition, the > 1200 km long fault represents one of the most important and exemplary intracontinental strike-slip faults in the world. We mapped fault trace geometry and interpreted paleoseismic trench exposures to characterize the seismogenic behavior of the ATF. We identified 2 geometric segment boundaries in a 270 km long reach of the central ATF. These boundaries define the westernmost Wuzhunxiao, the Central Pingding, and the easternmost Xorxol (also written as Suekuli or Suo’erkuli) segments. In this paper, we present the results from the Camel paleoseismic site along the Xorxol Segment at 91.759°E, 38.919°N. There evidence for the last two earthquakes is clear and 14C dates from layers exposed in the excavation bracket their ages. The most recent earthquake occurred between 1456 and 1775 cal A.D. and the penultimate event was between 60 and 980 cal A.D. Combining the Camel interpretations with our published results for the central ATF, we conclude that mul- tiple earthquakes with shorter rupture lengths ( 50 km) rather than complete rupture of the Xorxol Segment bet- ter explain the paleoseismic data. We found 2-3 earthquakes in the last 2-3 kyr. When coupled with typical amounts of slip per event (5-10 m), the recurrence times are tentatively consistent with 1-2 cm/yr slip rates. This result favors models that consider the broader distribution of collisional deformation, rather than those with northward motion of India into Asia absorbed along a few faults bounding rigid blocks. tinental collision (e.g., Tapponnier and Molnar, 1977; Peltzer et al., 1989; Yin and Harrison, 2000; Tapponnier et al., 2001; fig. 1). Two com- peting end-member models have developed that describe how northward convergence of India is accommodated by Tibet and Southern Asia. The first model divides this region into a few rigid blocks bounded by major faults (i.e. Altyn Tagh and Kunlun faults, fig. 1) (e.g., Avouac and Tap- ponnier, 1993; Tapponnier, et al., 2001). The second model has convergence accommodat- ed by crustal thickening and along numerous minor and a few major faults (e.g., England and Molnar, 1997). Models of the first type are con- 1015 ANNALS OF GEOPHYSICS, VOL. 46, N. 5, October 2003 Mailing address: Dr. J Ramón Arrowsmith, Department of Geological Sciences, Arizona State University, Tempe, AZ 85287-1404 U.S.A.; e-mail: ramon.arrowsmith@asu.edu 1016 Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle China Tibet Altyn Tagh Faul t 80° 100° 120° 40° 20° Kunlun fault Red River Fault Himalayan Frontal ThrustsIndia Fig. 1. Overview map showing major structures and topography in the India-Eurasia continental collision zone. Inset shows location of fig. 2. Figure generated from GTOPO30 1 km digital elevation data set using GMT (Wessel and Smith, 1995). sistent with higher slip rates ( 3 cm/yr for the ATF) while lower slip rates are consistent with the second type because deformation is distrib- uted instead of localized along a small number of discrete fault planes. Global seismicity databases (CNSS) and the Chinese Catalog of Historic Strong Earthquakes (2300 B.C. to 1911 A.D.) record no major earthquakes along the ATF. However, we know that large earthquakes have occurred because of moletracks seen on satellite imagery, surface breaks shown by Ge et al. (1992), and limited field investigation by western scientists (Molnar et al., 1987). Despite the clear evidence of active faulting, slip rate determinations are discrepant. Over timescales of 10 s-100 s of kyr, slip rates are inferred to be 2-3 cm/yr. Peltzer et al. (1989) inferred 3 cm/yr rates in the area of 84°E-92°E based on an assumption of initial Holocene age for offset landforms. Recent work has focused on determining the slip rate by dating offset Quaternary landforms using cosmogenic meth- ods. Meriaux et al. (1999, 2000) determined 3 cm/yr between 85°E and 89°E and 2-3 cm/yr near 90°E. Merieaux et al. (2003) determined a minimum of 20 ± 3 mm/yr at 94°E. Chinese earthquake geology studies (mapping, paleoseis- mic trenching, and 14C) inferred a minimum rate of 5 mm/yr along the entire ATF (Ge et al., 1992). Studies of deformation in the India- Eur- asian collision over geodetic timescales (5-10 years; Bendick et al., 2000; Shen et al., 2001; Wang et al., 2001) indicate that distributed deformation in Tibet is an important process because faults such as the ATF have slip rates on the order of 10 mm/yr (9 ± 5 mm/yr, Ben- dick et al., 2000; 9 ± 2 mm/yr, Shen et al., 2001). Holt et al. (2000) also used a rate of 10 ± 8 mm/yr along the ATF in their analysis of the velocity field in Asia. Two major strike-slip earthquakes have recently occurred in Northern Tibet. The 1997 Mw 7.6 Manyi earthquake ruptured for a dis- tance of about 170 km with peak left-lateral slip of about 7 m, apparently reactivating a Quaternary fault (fig. 2; Peltzer et al., 1999). The 2001 Mw 7.8 Kokoxili or Central Kunlun earthquake had a 300-400 km long rupture zone and sinistral slip up to 16.3 m (Lin et al., 2002; Klinger et al., 2003). It ruptured the eastern portion of the Kunlun Fault (van der Woerd et al., 2002). These earthquakes draw attention to the style of earthquake rupture along the great strike-slip faults of Asia and also to their response to the geodynamic setting. Because the historic record in this region is only com- plete for about the last 100 years, paleoseismo- logical methods are important for investigating the timing of large ruptures like these and those along neighboring faults. Earthquake displacements, surface rupture lengths, and recurrence intervals along major strike-slip faults such as the ATF have implica- tions for which geodynamic models better explain continental deformation and for region- al earthquake sequences in continental settings. Therefore, we conducted a study of the earth- quake geology of the ATF to characterize the seismogenic behavior of this major strike-slip fault. In this paper, we first highlight the geo- Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China 80 85 90 95 35 40 Kunlun Fault Tibet 6 7 8 White < 2500 years BP SSB data historic earthquake M scale 5 6 7 8 CNSS data Black > 2500 years BP 1997 Mw 7.6 Manyi earthquake Tarim Basin QaidamBasin Qimen Tagh Milan Ruoqiang NA TF Huatago Lapequan Alty n Ta gh F ault 1924 M>7 earthquakes 1951 ~M 6 earthquake Altyn Tagh Altyn Tagh Fault Akato Tagh Minfeng Subei Sulamu Tagh Tula 2001 Mw 7.8 Kokoxili earthquake ? Fig. 2. Overview map showing the Central Altyn Tagh Fault (ATF), location of study area (rectangular box, fig. 3), instrumental seismicity (Council of the National Seismic System (CNSS); http://quake.geo.berkeley.edu/cnss/cata- log-search.html), historic earthquakes (Chinese Catalog of Historic Strong Earthquakes 2300 B.C. to A.D., 1991), paleoearthquakes (State Seismological Bureau (SSB); Ge et al., 1992), and major towns. The Ge et al. (1992) data are plotted at the reported locations with their inferred magnitude and color-coded by age range. The events most relevant for this study are < 2500 years BP (white). The low seismic activity is unexpected because the ATF is thought to accommodate a significant portion of India-Asia convergence. The historical record also shows fewer earthquakes than anticipated. North Altyn Tagh Fault (NATF) mapping from Cowgill et al. (2000). Manyi earthquake surface rupture from Peltzer et al. (1999) and Kokoxili earthquake location and approximate surface rupture are from van der Woerd et al. (2002) and Lin et al. (2002). Base figure generated from GTOPO30 1 km digital elevation data set using GMT (Wessel and Smith, 1995). 1017 metric segmentation of the central ATF by di- viding it into three sections (westernmost Wuzhun- xiao, the central Pingding, and easternmost Xorxol-also written as Suekuli or Suo’erkuli). Secondly, we present new paleoseismic results from the Xorxol Segment at 91.759°E, 38.919°N at the Camel paleoseismic site. Finally, we re- view results from Washburn et al. (2001) and Washburn (2001), to develop earthquake rup- ture scenarios for the central ATF. The geomet- ric boundaries and 14C dates defining age ran- ges for paleoearthquakes imply that one large earthquake rupturing all 3 segments (270 km) is unlikely. Furthermore, multiple earthquakes with shorter rupture lengths ( 50 km) rather than complete rupture of the Xorxol Segment better explain the paleoseismic data. When coupled with typical amounts of slip per event (5-10 m), the recurrence times are consistent with 1-2 cm/yr slip rates. 2. Geomorphic and geologic framework The study area is located in Xinjiang, Qinghai and Gansu provinces of Western China between 89.619°E, 38.302°N and 92.523°E, 39.138°N (figs. 1, 2 and 3). This is an arid re- gion of alluvial fan and playa environments around 3000 m in elevation with neighboring glaciated peaks rising up to 5000 m. In the western part of the study area, the ATF passes through a 4000 m high, axial valley south of Liweiqiming Shan and descends across alluvial fans and saline lakes of the Wuzhunxiao area (figs. 2 and 3). From here, the ATF traverses bedrock and small perched alluvial basins (< 0.5 km2) within an isolated mountain range, mark- ed by the 4800 m Pingding Shan. From Gobi- ling, the ATF descends into the alluvium filled strike-slip rift of Xorxol Valley. This straight and distinct tectonic landscape continues east for Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle 39.00 E, 91.75 N 38.50 E, 89.75 N 38.50 E, 91.00 N 39.00 E, 92.50 N Liweiqiming Shan Kulesayi to Milan Magnei Lake Wuzhunxiao Arkateng Dayang Shan Pingding Shan Huatago (Huatougou) old road to Milan Ku Shui He Gobiling Camel site Figure 19 Figure 18 Bitter Seato Milan to Lapequa n Wuzhunxia o Segment Pingding Segment Xorxol Seg ment Xorx ol Va lley Figure 4 Explanation strike-slip fault Active normal fault Qm normal fault Debris flow Jeep road active trace (well located, inferred) Qm fault trace Tectonic lineament Interpreted trace River Fig. 3. Fault trace map of Central Altyn Tagh Fault (ATF) showing geometry of most recent and Late Qua- ternary rupture traces. Major roads, location of geographic landmarks, and towns are also shown for reference. This portion of the ATF is predominantly left lateral and has well-developed strike-slip topography (e.g., Xorxol rift valley). However, it has not produced instrumentally recorded earthquakes > M 6. The active trace traverses alluvial flats around Lake Wuzhunxiao and the Bitter Sea and bedrock in the Arkateng and Pingding mountain ranges. Two geometric boundaries are identified: the 7 km wide Arkateng Restraining Bend and the 4 km wide Pingding releasing step. These geometric boundaries de fine the westernmost Wuzhunxiao, the Central Pingding, and the easternmost Xorxol segments. The boundaries serve as a framework for interpreting offset data and esti- mating rupture extent because such discontinuities commonly act as rupture barriers. 1018 over 110 km (figs. 2 and 3). Proterozoic intru- sive rocks and slivers of Paleozoic-Mesozoic granitoid, ultramafic, marine carbonate, and ophiolitic rocks comprise the mountains in the study area (Cowgill et al., 2003). Thick Neoge- ne boulder gravels to silts encircle the Pro- terozoic-Mesozoic mountains and fill adjacent basins. Intense folding of these strata in the Qaidam Basin and throughout the field area indicate a contractional regime (regionally) fol- lowing deposition of the Neogene units. Relatively thin Quaternary gravels and silts (pied- mont deposits and primary and reworked loess) overlie the Neogene strata. 3. Fault trace mapping and paleoseismic site identification In the field, we mapped the fault traces that comprise the individual faults that have had recent rupture based on geomorphic age of the scarps, disruptions of Holocene deposits and landforms. The active trace designation gener- ally refers to surface breaks created by the last 1-3 earthquakes. We have identified the active traces and most Quaternary faults, but may not show all of the early Quaternary faults because we are primarily concerned with the geometry and earthquake history for the last few kyr. In fig. 3, faults labeled Qm are similar in age to or slightly younger than surfaces and deposits of Late Pleistocene and Early Holocene inferred age. These were mapped as unit Qm and are the most extensive of the Quaternary deposits and form broad floors that depositionally abut adja- cent hill sides and often de fine a sharp break in slope. They generally have a few to 10 s of m incision and may be covered with varnished clasts. Tectonic lineaments refer to suspicious geomorphic alignments that are most likely faults older than Qm. Our mapping approach was to inspect the principal active traces as identified on CORONA imagery. We mapped the main fault traces, the Quaternary units discussed above, and occa- sionally obtained bedding or other structural attitudes. We typically covered a 5 km wide zone centered over and parallel to the active trace to identify relevant structures and geolog- ic relations. We made many small excavations with shovels to clarify stratigraphic relation- ships and expose bedding surfaces. These small- er excavations helped add confidence to our geo- morphic and stratigraphic-based paleoseismic site decisions. We logged natural and hand-dug exposures, and documented offset (both vertical and horizontal) landforms and strata. While mapping the fault trace, we looked for sites where landforms indicated that steady accumulation of distinctive sediments would preserve evidence for past ruptures. In a num- ber of excavations at such sites, broken and tilt- ed beds overlain by continuous beds and upward terminations of fault traces identify past earth- quakes in excavations. Disruption seen at the same stratigraphic level at several locations and on both trench walls qualifies as a distinct event (after Grant and Sieh, 1994). Organic matter within the stratigraphy was dated using 14C and fine silt was dated by infrared stimulated lumi- nescence (IRSL) (Forman, 1999). We use the 2 calibrated sample ages and OxCal v3.5 (Ram- sey, 2000) for timing constraints. We also use a Monte Carlo program that uses ordering con- straints to «trim» probability distributions of ra- diocarbon samples and then yields probabilities for events (Hilley and Young, 2003). This approach allows for more objective event corre- lation between paleoseismic sites. 4. Fault trace geometry A small-scale representation of the active traces that we mapped is shown in fig. 3. For convenience, we refer to locations along the mapped traces in terms of distance from the western end of the study area. Washburn et al. (2001) described the western 150 km of the area shown in fig. 3 and established the geo- metric segmentation. Our subsequent map- ping added another 120 km toward the east (Washburn, 2001). A view along the strike of the active trace in fig. 3 shows a prominent right bend (65 to 85 km) and left step (97 km) in the generally straight Central Altyn Tagh Fault. We followed the conventions of dePolo et al. (1991) and Knuepfer (1989) and identi- fied these discontinuities as geometric bound- 1019 Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China aries which in turn define the westernmost Wuz- hunxiao Segment, the Central Pingding Seg- ment, and the Eastern Xorxol Segment. These segments are not necessarily rupture segments; however, this approach is useful for characterizing fault zones because geometric boundaries may be locations of significant change in slip or moment release (e.g., Wesnousky, 1989; Harris and Day, 1993, 1999). Historical rupture data show that rup- ture endpoints commonly coincide with geometric boundaries and these boundaries appear to retard rupture (Knuepfer, 1989). The major geometric boundary in the study area is the Pingding Releasing Step. Here, the active trace jumps from the southern side of Pingding Shan 4 km north to the active trace that continues to the Xorxol Valley (fig. 3; Wash- burn, 2001; Washburn et al., 2001). The 3-4 m vertical offset of Holocene-age fan on the sout- heast side of Pingding Shan, a prominent scarp on the east face of Pingding Shan, and the pres- ence of recent landslides indicate that the nor- mal fault that links the bounding strike-slip faults ruptured recently. Our mapping stops at 268 km in a straight and simple portion of the Xorxol Segment (fig. 3). This 070° striking simple geometry continues eastward, but just off the eastern edge of fig. 3, around 300 km, a prominent northern strand de- velops (fig. 2). We infer that the Xorxol geomet- ric segment is at least 200 km long and remark- ably straight (Washburn, 2001). 5. Xorxol Segment paleoseismic investigation Overview Seismogenic properties such as fault straight- ness and maximum displacement during the last earthquakes at specific sites along the central ATF were presented above and in Washburn et al. (2001) and Washburn (2001). The next steps in characterizing the seismogenic behav- ior of the Central Altyn Tagh Fault are determi- nations of recent earthquake timing and esti- mates of magnitude. In order to establish the timing of past earthquakes, we made 10 exca- Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle Table I. 14C data for the Camel paleoseismic site. All analyses are standard AMS dates from CAMS facility at LLNL. .13C values are the assumed values according to Stuiver and Polach (1977). Values measured for the material itself are given with a single decimal place. The quoted age is in 14C years using the Libby half life of 5568 years and following the conventions of Stuiver and Polach (1977). 14C concentration is given as fraction modern, D14C, and conventional 14C age. Field sample CAMS # .13C Fraction ± D14C ± 14.C age ± 2 calibrated name modern age (*) (cal A.D.) C 1 69687 – 25 0.8923 0.0035 – 107.7 3.5 920 40 1024-1209 C 7 69688 – 25 0.9592 0.005 – 40.8 5 330 50 1456-1650 C 8 69689 – 25 0.9855 0.0048 – 14.5 4.8 120 40 1676-1939 C 10A 69690 – 25 0.7921 0.0034 – 207.9 3.4 1870 40 60-240 C 16 69691 – 25 0.8912 0.0043 – 108.8 4.3 930 40 1022-1191 C 20 69692 – 25 0.8639 0.0047 – 136.1 4.7 1180 50 717-980 C 22 69693 – 25 0.9803 0.0039 – 19.7 3.9 160 40 1662-1950 C 24 69694 – 25 0.7376 0.0042 – 262.4 4.2 2440 50 762-403 B.C. C 25 69695 – 25 0.9874 0.0039 – 12.6 3.9 100 40 1678-1937 C 26 69696 – 25 0.8929 0.0045 – 107.1 4.5 910 50 1022-1221 C 32 69697 – 25 0.6212 0.003 – 378.8 3 3820 40 2404-2140 B.C. C 45 69698 – 25 0.5745 0.0027 – 425.5 2.7 4450 40 3337-2924 B.C. SL rabbit 69699 – 25 1.1245 0.0054 124.5 5.4 > Modern > 1950 SL camel 69700 – 25 1.2074 0.0058 207.4 5.8 > Modern > 1950 SL stick 69701 – 25 1.1898 0.0052 189.8 5.2 > Modern > 1950 (*) Determined with Calib 4.2 (Stuiver et al., 1998). Does not include curve intercepts comprising < 2% of the probability distribution. 1020 1021 Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China 195 197.5 200 21o 19o 10o ~30o 12o Playa Hill 10o N O trench site N Explanation strike-slip fault Qm normal fault active trace Qm fault trace jeep road Along-fault distance (km) Fig. 4. Fault trace map showing main strike-slip rupture cutting Playa Hill, the eastern end of the playa, small alluvial fans at the base of north facing escarpment, and the Camel paleoseismic site. Note the 058°-striking secondary strike- slip fault, which intersects the main trace at (N) and dies out within 2 km along strike to the southwest. The Qm nor- mal faults located along the piedmont and range front to the south have 1-8 m high degraded scarps indicating they have not ruptured recently. The only exception is a slightly sharper scarp seen at (O). However its more distal position may explain this morphologic difference. Base image is a rectified CORONA photograph. See fig. 3 for location of this site. vations somewhat evenly spaced along the field area (fig. 3). Interpretations of logs of these excavations are presented in Washburn (2001). In addition, early results from 1998 and 1999 excavations are presented in Washburn et al. (2001). Here, we present observations and infe- rences from the Camel site where the record of repeated ground rupture was best preserved (figs. 3 and 5). All of the 14C data are presented in table I. 5.1. Camel site setting The Camel excavation has an understand- able relationship with neighboring faults, re- veals a clear stratigraphy, and preserves abun- dant organic material. This site yields the most robust timing constraints for recent earthquakes of all our paleoseismic sites. A secondary 058° striking, strike-slip fault intersects the active trace one kilometer west of the Camel trench site and normal scarps cut the Qm piedmont 2.5 km south of this site (fig. 4). These secondary structures could allow for an incomplete earth- quake record in excavations across the active trace. Fortunately, degraded surface rupture implies that the 058° strands probably did not rupture during the last earthquake and the nor- mal faults also probably have not ruptured dur- ing the last several earthquakes. Additionally, it is unlikely that the 058° system would rupture independently because of its close proximity to the active trace. In support of this assumption, the portion closest to the active trace has sharp- er scarps (shown in red on fig. 4) and preserves geomorphic offsets (moderately reliable 3.5 m), while farther west there is no evidence for re- cent surface rupture. Therefore, we hypothesize that the last earthquake on the main trace acti- vated the only the easternmost part of the 058° fault. Furthermore, because the site emphasizes earthquake timing, a broader rupture zone should provide better and redundant constraints than a narrow one. See Washburn (2001) for more detailed geometric observations and further 1022 Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle 0 25 m N Camel trench siterectified kite photos N contour interval 1.0 m relative to arbitrary datum 0 25 m trench fault trace topographic breakline channel strandline A B C C B B B A B B strandline 1000 1005 1010 1015 1020 a b Fig. 5a,b. a) Rectified kite aerial photograph mosaic showing the geomorphic setting and geometry of surface rupture traces at the Camel trench site. Note the playa strandline at the north end of the trench, small dunes just south of the strandline on the east side of the trench, Neogene hills at the southern end of the trench, and allu- vial fans at the base of these hills. b) Topographic map quantifying the landforms at this site. The letters (A, B, C) on both figures identify corresponding fault traces discussed in the text. review of secondary structures at the Camel site. Four geomorphic processes interact to cre- ate an environment that is conducive to paleo- seismology (fig. 5a,b). On the northeast side of the trench, small dunes are visible. They create distinctive cross bedded deposits comprised of well rounded, medium to coarse-grained sand (fig. 9). Small alluvial fans on the west side of the trench deposit muddy sand and fine gravels units seen in the trench. To the south, low hills of indurated Neogene silts comprise the «base- ment» units exposed in the southern part of the trench. The playa edge contributes fine silt and muddy strata to the northern end of the trench and introduces organic material at regular inter- vals. Individual pieces of wood and dung are sporadically exposed throughout the trench stra- ta. In addition, playa strandlines are also preserved in numerous places in the strata. A modern strandline is seen in fig. 5a,b; it formed by shallow inundation of the playa and subsequent rafting of organic material towards the leeward end of the playa. Regular seasonal winds blow to the east and have subsequently concentrated woody material, flower parts, cam- el dung, and rabbit dung in < 1 m wide by 5 cm high by 10 s of m long deposits ringing the playa edge. Radiocarbon dating of the modern strandline allows us to evaluate the residence time of these deposits and therefore possibili- ty for detrital ages from the trench strandlines (table I). Although 2 unconformities were identified within the trench, their planar con- tact implies that erosion was not significant. Furthermore, the presence of strandlines at different stratigraphic levels in the trench gen- erally implies that nearly continual deposition has been occurring at this site. This is key to preservation of a complete earthquake record. The intermixing of these various sedimentary deposits results in distinctive, interpretable, and datable stratigraphy at the Camel trench, which in turn allows us to determine earth- quake timing at this site (figs. 6, 7a,b and 9). The active trace near the trench sites con- sists of two prominent fault traces, 12-15 m apart that have been almost entirely degraded in places (figs. 4 and 5a,b). A possible third trace breaks along the Neogene hill front. We exca- vated a 34 m long trench across all three of these traces to minimize the potential for miss- ing an event. The trench stratigraphy closely corresponds with the geomorphic processes at this site (figs. 6, 7a,b and 9). Playa and dune sediments thicken toward the northern end of the trench, while coarser sand and gravel (allu- vial fan) units and Neogene bedrock dominate the southern end. We identified 6 fault zones with recurrent disruption of the stratigraphy and number these starting with fault zone 1 in the north end of the trench. We also number strati- graphic units upward and show their contact relationships in fig. 7a,b. Projection of pre- served portions of active trace indicate that the northern trace (A on fig. 5a,b) corresponds with fault zones 1 and 2 and the southern trace (B on fig. 5a,b) forms fault zones 3 and 4 (fig. 9). 5.2. Camel Most Recent Earthquake (MRE ) Faults in zones 1-4 cut units 100 through 130 and sometimes 140, however fractures in fault zones 1 and 2 cut higher units (fig. 9). These frac- tures show no offsets across them and the often show downward terminations at medial strati- graphic levels. These characteristics and proximi- ty to the playa show the fractures are probably desiccation-mud cracks and non-tectonic. From fault zones 1 to 4, the silts compris- ing unit 140 grade to coarse sand and fine gravel and erode the top of dune unit 130 (figs. 7a,b and 9). On the east wall, fault zone 1 shows faults terminating in the middle of unit 140, while zone 2 shows faults capped by unit 140. In portions of fault zones 3 and 4, unit 140 is missing. The west wall shows the same general trend of thinning and eventual disap- pearance of unit 140 south of fault zone 4. We interpret unit 140 as a time transgressive unit that was deposited before and after the last earthquake based on these cross-cutting rela- tionships and because the faults seen in zones 1-4 ruptured the surface based on the fault traces on either side of the trench. Unit 140 is probably a sheet flow deposit coming from runoff-driven sediment transport on the small alluvial fan. The base of 150 is the event hori- zon for the last earthquake in these fault zones. 1023 Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China On the east wall of zone 5, a single fault surface cuts unit 130 and is covered by unit 160 (fig. 9). Unit 150 pinched out north of this point, so these cross cutting relationships are similar to those in fault zones 1-4. On the west side, unit 123 is disrupted by the same fault, however a < 3 cm wide fissure implies higher rupture through unit 127. In fault zone 6, a sin- gle fault cuts to the same level (unit 127). South of the zone 6, numerous faults cut Neogene base- ment and strongly indurated Quaternary gravels. Compared to fault zones 1-4, the smaller fault zones and amount of disruption suggest that fault zones 5 and 6 slipped only a small amount and may be reactivated older principal rupture zones. The consistent level of upward fault ter- minations indicates that the base of 150 is the event horizon for the last earthquake. Units 130 and 140 represent the ground surface at this time. The southern end of the trench has thinner units and thus does not record the surface rup- tures as well as the lower gradient northern end. Note the gently convex warp in the strata on east wall of fault zone 3. This is apparently a buried moletrack formed during the last earthquake and subsequent erosion beveled the top off. Below it gently folded strata are strata deformed prior to the last earthquake. Erosion over the top of fault zone 6 makes cor- relation of fracture events there difficult. However, fractures do terminate at the base of a unit correlat- ed to 140 in the northern portion of the trench. 5.3. Camel penultimate earthquake Different levels of induration and an irregu- lar contact sometimes showing centimeter deep rills identify an unconformity between units 105 and 100. The west side of fault zone 2 shows folded and faulted strata of unit 100 overlain by undeformed silts and sands of unit 105. This relationship is also seen on the other side of this fault zone and indicates the penultimate rupture 1024 Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle 170 160 150 140 MRE 130 127 125 123 115 Fig. 6. Photograph of fault zone 1 on east wall of Camel trench (fig. 9). Note the metric tape measure on the trench wall for scale. Photograph shows, from top to bottom, units 170, 160, 150, 140, 130 (thick crossbedded sand), 127, 125, 123, and 115. Note that the faults have upward terminations between the top of unit 130 and the base of unit 150. This relationship shows rupture resulting from the last earthquake at the Camel site. 1025 Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China (fig. 9). Additional evidence for the penultimate event is seen in heavily deformed strata over- lain by less deformed strata in fault zone 3. Both sides of fault zone 4 show upward fracture terminations and increased local erosion at the 100-105 unit contact. Fault zone 5 shows a clear north side up sense of offset below the unconfor- mity. Therefore unit 100 is the youngest unit cut by the penultimate event and the unconformity is the penultimate event horizon. As noted above, the MRE and the penultimate events produced local highs (moletracks) that are espe- cially well manifest on the east side of fault zone 3, but also evident in fault zones 4 and 5. Each time they formed, they were erosionally beveled. N2 x x x x x C22 1662-1950 A.D. C25 1678-1937 A.D. C7 1456-1650 A.D. C26 1022-1221 A.D. C16 1022-1191 A.D. C20 717-980 A.D. C1 1024-1193 A.D. C10 60-240 A.D. C24 762-403 B.C. C32 2404-2140 B.C. C45 3337-2924 B.C. uncf uncf uncf 170 160 150 140 130 125 123 127 115 110 105 100 90 85 80 C8 1676-1939 A.D. younger than MRE older than MRE/younger than pen.x RADIOCARBON SAMPLES older than penultimate event uncf unconformity 115 unit number 4000 cal B.C. 2000 cal B.C. cal BC/cal A.D. 2000 cal A.D. C22 160 ± 40BP C25 100 ± 40BP C8 120 ± 40BP C7 330 ± 50BP C26 910 ± 50BP C16 930 ± 40BP C1 920 ± 40BP C20 1180 ± 50BP C10 1870 ± 40BP C24 2440 ± 50BP C32 3820 ± 40BP C45 4450 ± 40BP Calibrated dates Radiocarbon dates and calibrated curves from OxCal 3.5Unit # Fig. 7a,b. a) Stratigraphic column of Camel trench showing relative unit thickness, unit number, major uncon- formities, and stratigraphic location of radiocarbon samples. Diamonds show dates above the MRE, Xs show dates older than the MRE, and Os are older than the penultimate event. Weak evidence in the trench near fault zone 2 suggests an event at the level of unit 125. b) Probability distributions for calendar dates determined from the radiocarbon age for each sample using OxCal v3.5 (Ramsey, 2000). a b 5.4. Other possible earthquakes in the Camel trench There is limited evidence for a possible event occurring before the MRE, but after the penultimate event. Two 50° south dipping fault traces apparently terminate in unit 125 and 130 on the east wall of fault zone 2 (fig. 9). The die out of the tip in unit 130 is most like- ly due to dilatation of the sand unit. The base of unit 127 could be a rupture horizon that lies between the MRE and penultimate event hori- zons. However, we do not see similar relation- ships on the west wall. There, the two faults have tips within unit 130; therefore it is likely that the tipline of the faults plunge eastward. A fault trace splays into three strands that are apparently overlain by continuous bedding of unit 130 on the east wall of fault zone 3 (fig. 9). This potential rupture horizon could also represent an intermediate earthquake. The cor- responding fault on the west wall ruptures to a higher level, supporting the idea of upward die out or plunging tip line (fig. 9). We also see a fault trace that ruptures to an intermediate level in the southern part of fault zone 4. On the east wall this fault is capped by unit 110, but breaks to the top of 110 on the west wall (fig. 9). This situation is similar to fault zone 3 and could also represent an intermediate event. However, the different rupture levels could also result from a plunging tip line of an MRE fault trace. If the intermediate level terminations seen in fault zones 2, 3 and 4 represent an intermedi- ate event, then it produced only small amounts of disruption of the trench strata. In contrast, the MRE and penultimate events produced intense disruption of the strata along the entire > 30 m length of the trench. Therefore, the lim- ited deformation, seen as only 3 faults, indi- cate that this intermediate event was different and had lower magnitude than the MRE and penultimate events. We favor an interpretation of upward die out and/or plunging tip line of MRE fractures because the faulting relation- ships are not the same on both trench walls and the much smaller degree of deformation com- pared to ruptures seen in other trenches in the study area. More substantial evidence for a pre-penulti- mate event is seen 1 m north of and in the bot- tom of fault zone 5. Here, faulted and tilted strata are overlain by flat lying bedded sand and fine gravel of unit 85. These relationships are suggestive of an earthquake before the penulti- mate event, however the trench was not deep enough to confirm this event. Fault zone 6 preserves evidence of the fail- ure of large blocks of the local Neogene base- ment indurated silts into mostly eolian sands. The blocks are strongly fractured and in combi- nation with the windblown sands represent a slightly different environment than the eolian sand and alluvial fan muddy sands and gravels to the north, further making correlation of e- vents with the rest of the trench difficult. 5.5. Camel 14C dating We collected 57 samples of organic mate- rial from the trench and dated 12, as well as 3 from a modern strandline. Because much of the organic material in the Camel trench was represented in modern strandlines, and often deposited as strandlines, we dated the modern ones to assess their ability to define the depo- sitional age. The 3 strandline samples are all modern. Thus, the 14C dates in the trench prob- ably do not include a significant pre-deposi- tional age. The 14C data are presented in table I and the probability distributions from calibra- tions using OxCal (Ramsey, 2000) are shown in fig. 7a,b. An 8 by 5 by 3 mm piece of wood with intact bark (C8) from unit 150 provides an upper bound of 1676-1939 cal A.D., while a rabbit dropping (C7) yields a lower bound of 1456-1650, for the last earthquake in this area. Likewise, an elongate rabbit dropping (C20) places an upper bounds of 766-980 cal A.D. and 10 small pieces of woody debris were mixed (C10) to provide a lower bounds of 60-240 cal A.D. for the penultimate event. The oldest strata exposed in the trench are be- tween 3337-2924 cal B.C. based on a radio- carbon date from piece of wood (C45). Because of the complete stratigraphic a- greement within all 12 samples, we used OxCal 1026 Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle 1027 Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China to refine the constraints on event timing (Ram- sey, 2000). This program uses stratigraph- ic ordering constraints to eliminate overlap- ping age probabilities, which refines our earthquake timing brackets. We used 1950 A.D. as an upper event boundary because the world wide seismic network was able to detect moderate sized earthquakes by this time. OxCal is most useful for samples that are closely spaced stratigraphically and yield similar radiocarbon ages. C22 is located near the top of unit 160 while C25 in the lower part of this unit and C8 is located in unit 150. These samples all yield similar calendar ages (table I), however their stratigraphic position allows OxCal to cut off overlapping portions of their probabil- ity calendar age distributions. This is useful be- cause this part of the calibration curve results in a wide age range (e.g., 1676-1939 for C8). After running OxCal, C8 yields a narrow up- per bound of 1676-1775 for the Camel MRE. The lower bound of 1456-1650 cal. A.D. remains unchanged because the next lowest sample (C26 1022-1221 cal A.D.) is much older and therefore probability distributions do not overlap. The lower bound on the pe- nultimate earthquake is also unchanged for the same reason. Above the penultimate event horizon, unit 125 contains 2 radiocarbon sam- ples, however they are spaced over an 8 m horizontal distance and sit at the approximate same stratigraphic level. Therefore we cannot apply ordering constraints to these samples, although C26 in unit 127 does shift their age distributions 30 years older. Figure 8a,b shows the earthquake timing constraints at a number of sites along the cen- tral ATF. The ages displayed in fig. 8a are cal- ibrated but otherwise not refined with strati- 0 50 100 150 200 250distance (km) C al en d ar Y ea rs -2500 -1500 -500 500 1500 K ul es ay i 2- 6- 17 2- 6- 18 C am elExplanation 14 MRE x x x most recent event younger than MRE older than MRE/younger than pen. older than penultimate event earthquake from Active Altun Fault monograph (Ge et al.) 2 sigma calibrated age of C (solid line) and IRSL (dotted line) samples older than earliest event nosym. -2500 -1500 -500 500 1500 Camel Probability 0.01 0.0 C al en d ar Y ea rs Explanation MRE Penultimate event Fig. 8a,b. a) Earthquake timing constraints from stratigraphic relationships and 14C and Infrared Stim- ulated Luminescence (IRSL) samples (from this study and Washburn et al., 2001). Note the active trace map above plot A for reference. These are the calibrated ages without OxCal or other refinements from stratigraphic ordering. b) Monte Carlo proba- bility plot for the Camel paleoseismic site earthquake timing data (Hilley and Young, 2003). The area under the curves represents the probability distribution of the age of past earthquakes. Sharper peaks represent well-constrained events and show the refinement of the Camel event ages. a b Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle 1028 Fig. 9. Trench logs for the Camel paleoseismic site at 198 km. Most faults terminate at the upper blue line, the event horizon for the last earthquake at this site. The penultimate earthquake is identified by upward fault ter- minations and folded bedding overlain by continuous strata, which is marked by the lower blue line. This penul- timate event horizon coincides with a small unconformity (yellow line) identified by different levels of indura- tion and an irregular contact with cm deep rills in places. Stratigraphic column of Camel trench at lower left shows relative unit thickness, unit number, major unconformities, and stratigraphic location of radiocarbon sam- ples. See discussion in text for further interpretation and timing constraints on earthquakes seen in this trench. N2 x x x x x C22 1662-1950 A.D. C25 1678-1937 A.D. C7 1456-1650 A.D. C26 1022-1221 A.D. C16 1022-1191 A.D. C20 717-980 A.D. C1 1024-1193 A.D. C10 60-240 A.D. C24 762-403 B.C. C32 2404-2140 B.C. C45 3337-2924 B.C. uncf uncf uncf 170 160 150 140 130 125 123 127 115 110 105 100 90 85 80 C8 1676-1939 A.D. Stratigraphic column Unit # Sample calibrated date younger than MRE older than MRE/younger than pen. older than penultimate event x Radiocarbon samples uncf unconformity 115 unit number Key Explanation matrix supported massive sand gravel clast supported gravel bedded sand eolian sand silt sandy silt prominent strand line event horizon unconformity 14C sample fault, fracture Neogene bedrock east wall C32 (2404-2360, 2354-2140 cal B.C.) C26 (1022-1221 cal A.D.) C25 (1803-1907, 1678-1744 cal A.D.) C16 (1022-1191 cal A.D.) C1 (1024-1193 cal A.D.) west wall fault zone 1fault zone 2 fault zone 3 fault zone 4 fault zone 5 fault zone 6 C24 (762-678, 671-608, 599-403 cal B.C.) C7 (1456-1650 cal A.D.) C8 (1802-1939, 1676-1765 cal A.D.) C20 (766-980 cal A.D ) C22 (1662-1711, 1717-1886, 1911-1950 cal A.D.) C45 (3337- 3208, 3194-3150, 3140-3008, 2987-2924 cal B.C.) C10 (60-240 cal A.D.) 0 5m prominent strand line 14C sample Explanation event horizon unconformity fault, fracture matrix supported clast supported bedded sand massive sand gravel gravel eolian sand silt sandy silt Neogene bedrock 120 130 SOUTH NORTH Unit numbersh Unit numbers 170 160 155 140 130 123 115 105 100 170 160 130 115 110 150 140 90 120 130 100 105 170 160 155 140 130 123 115 105 100 Unit numbers 80 90 100 110 160 170 Fig. 6. (continued) 1029 Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China graphic ordering constraints using OxCal. Hil- ley and Young (2003) have developed a 14C age refinement method that is similar to OxCal and produces probability curves for events with a Monte Carlo-based analysis of 14C age calibration and stratigraphic ordering. We show the probability curves for the last two well-defined Camel events determined using the Hilley and Young (2003) method in fig. 8b. 6. Discussion 6.1. Earthquake ruptures along the Central Altyn Tagh Fault Correlation of earthquakes between paleo- seismic sites can lead to misinterpretation be- cause not all strands in wide fault zones may be activated during an earthquake. We have elimi- nated most of the potential for this scenario by thorough evaluation of paleoseismic sites and excavation across only the sharpest, most active traces. However, «dogtail» rupture terminations in the middle of geometric segments and may not rupture through all excavations in question (e.g., Ward, 1997). Therefore trenches showing overlapping geochronologic earthquake con- straints that are only 10 s of kilometers apart may not record the same earthquake. We there- fore utilize other observations to help clarify the rupture history. We know the peak surface displacement at a given site from the offset data (Washburn, 2001; Washburn et al., 2001) and we have con- strained earthquake timing at a number of paleo- seismic sites. With these data, we can determine crude recurrence intervals, because we have only 1-3 events in each trench. We can also make magnitude estimates using the displace- ment data. Better estimates of moment magni- tude (Mw) can be achieved by using Surface Rupture Length (SRL) (Wells and Coppersmith, 1994). We estimate surface rupture length by utilizing geometric constraints, offset distribu- tion, surface rupture morphology, timing con- straints from the trench sites. Figure 8a,b shows that the last earthquake at Kulesayi, 2-6-17, 2-6-18, and the Camel site all occurred at the approximate same time, be- tween 1215-1775 cal A.D. (data from Washburn et al., 2001, and this paper). However, it is unlikely that one large earthquake ruptured the entire field area between 1215-1775 cal A.D. In such a case, the MRE would have jumped two geometric boundaries and ruptured from the Kulesayi trench (23 km) to at least the Camel trench (199 km) and likely farther based on the offset distribution and sharp surface rup- ture that continues to > km 240. This scenario yields a SRL of > 220 km and Mw > 7.8. However, consideration of details of the age control from the excavations and the geometric boundaries makes it unlikely, and the relatively high slip im- plied ( 10 m) is not consistent with our infer- red offset distribution (maximum offset 7 m; Washburn, 2001; Washburn et al., 2001). It is possible that the 4 km Pingding left step barred rupture. The Arkateng Restraining Bend and associated complexities around Dayang Shan are also potential rupture barriers fig. 3. Furthermore, the discontinuous and degraded surface rupture around these boundaries argues for smaller earthquakes that did not rupture the whole field area. Therefore we separate Wuzhunxiao and Pingding events from Xorxol events based on surface rupture morphology and geometric boundaries. Consequently, the last earthquake at Kulesayi occurred between 1215 and 1750 A.D. and may have produced the 3-5 m offsets on the west side of Lake Wuzhunxiao (fig. 8a,b; Washburn et al., 2001). Along the Xorxol Segment, the timing con- straints imply that the last earthquake ruptured its western portion before the eastern portion broke. The last earthquake in the west ruptur ed 1270 and 1430 cal A.D. (from Bitter Sea published results of Washburn et al., 2001; fig. 8a,b). This rupture may go from western end of the Xorxol Segment to around the Bitter Sea ( 50 km SRL) with peak displacement of 7 m (Washburn, 2001; Washburn et al., 2001). The younger Camel MRE occurred between 1456 and 1775 cal A.D. (fig. 8a,b). It ruptured from somewhere east of the Bitter Sea through the Camel trench and likely farther based on the rounded and semi continuous moletracks be- tween 247 km to the edge of the study area (268 km, fig. 3). The degraded nature of the surface rupture suggests that this reach is either the eastern terminus of the Xorxol MRE or that ero- sion rates, which are detectably higher in this area, led to the more degraded surface rupture. Higher erosion rates are plausible because of the position of the Eastern Xorxol Valley relative to the «mouth» of the Tarim Basin only 7 km north of the playa at 265 km on fig. 3. The mouth allows more frequent storms and intense daily winds to blow through Eastern Xorxol Valley and more animal access and potential degradation. In summary, the last earthquakes along the Central Altyn Tagh Fault: 1) occurred between 1215 and 1775 cal A.D.; 2) had maximum displacement of 5 m on the Wuzhunxiao and Pingding segments and 7 m on the Xorxol Se- gment (Washburn, 2001; Washburn et al., 2001), and 3) had surface rupture lengths between 50 and > 150 km with corresponding magnitudes of low Mw 7 s to high Mw 7 s, respectively. 6.2. Implications for slip rate based on the size and frequency of earthquakes in the last 3 kyr Although the earthquake data do not tightly constrain the surface rupture length or Mw of the MRE, the recurrence times from trenches show good agreement. We find 2 earthquakes in the last 0.8-2.0 kyr at the Camel site, 3 events in the last 2.4-3.0 kyr at the Bitter Sea, and 2 Kulesayi events in the last 0.8-2.2 kyr (this paper; Washburn, 2001; Washburn et al., 2001). Although we have a limited number of events for robust recurrence calculations, taking averages from these data give repeat times of 0.7 ± 0.3 kyr for the Camel site, 0.9 ± 0.1 kyr for the Bitter Sea trenches, and 0.75 ± 0.35 kyr for the Kule- sayi site. These frequencies are consistent with the 2 earthquakes in the past 3 kyr reported by Ge et al. (1992) in the study area and are gen- erally similar to the 800-1000 year recurrence times estimated for the Kunlun Fault (van der Woerd et al., 1998). Given the limited number of earthquakes in this study, we emphasize the tentative results derived from relating earthquake recurrence to slip rates. However, systematic underinterpre- tation or rupture location variability must be invoked for significantly higher recurrence rates or slip/event. If a large earthquake oc- curred in the near future, our recurrence rates would be maxima and any inferred slip rate would increase. To compare these paleo-earthquake data with other findings along the ATF we cast slip rate in terms of average recurrence intervals by using 5 and 10 m of slip per event because these values bracket the peak offsets recorded for the MRE(s) (Washburn et al., 2001). Using: 1) 10 mm/yr determined from repeat GPS surveys (Bendick et al., 2000; Shen et al., 2001) yields recurrence intervals of 0.5 and 1.0 kyr for 5 and 10 m of slip per event, respectively; 2) 3 cm/yr determined by large scale post glacial offsets and offset terraces and moraines in areas west of the study area (Peltzer et al. 1989; Meriaux et al., 2000) yields recurrence times of 0.2 and 0.3 kyr. Figure 10 shows a compilation of the various slip rates that have been inferred for the central ATF. Comparison of the different recur- rence intervals show that the paleoearthquake data most closely agree with a 1-2 cm/yr slip rate because we would expect to see 2-3 times the number of earthquakes in our paleoseismic investigations and the faulting and sedimentary record in them is not consistent with so many missing earthquakes. An important issue in slip rate comparisons is the observed eastward decrease in slip rate along the ATF (fig. 10). Peltzer et al. (1989) in- ferred 3 cm/yr rates in the area of 84°-92°E. Meriaux et al. (1999, 2000) determined 3 cm/yr between 85° and 89°E and 2-3 cm/yr near 90°E. Ryerson et al. (2003) report a well constrain- ed 18 mm/yr at the Pingding Shan (90.5°E). Meriaux et al. (2003) determined a minimum of 20 ± 3 mm/yr at 94°E. Ryerson et al. (2003) also report 12.5-15 mm/yr at 95°E. Meyer et al. (1996) concluded that slip rates east of 96° were 4 ± 2 mm/yr, consistent with decreasing slip rate along the ATF and accommodation of the slip gradient by thrusting along the south side in the Qaidam Basin (also suggested by Peltzer et al., 1989; Chen et al., 1999). Meyer et al. (1998) observed that the thrusts branch from the ATF and are thus kinematically linked; and van der Woerd et al. (2001) showed the decrease in slip rate and geometry of the fault junction at the Altyn Tagh Fault-Tanghenan Shan Thrust triple junction. 1030 Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle This study covers the portion of the ATF from 89.6°E to 92.5°E. It is within the area covered by the higher rates in Peltzer et al. (1989) and spans the Meriaux et al. (1999) and Ryerson et al. (2003) site near the Pingding Shan (18 mm/yr). Our best earthquake recur- rence data come from the Camel site is at about 91.75°E, so we may begin to observe a decrease in slip and thus recurrence rate due to interac- tion between the ATF and the adjacent active thrusts of the Qaidam Basin (for example, note that the active northwest striking Yousha-Shan thrust fault of Chen et al. (1999) approaches the ATF near the Gobiling). Although the re- sults of Ryerson et al. (2003) and our inferred slip rate from near the Pingding Shan overlap (fig. 10), the 20 ± 3 mm/yr minimum rate at 94°E (Meriaux et al., 2003) implies a higher earthquake recurrence rate than we inferred for the Xorxol Segment. The geodetic investigations are not suffi- ciently dense to image variations in slip rate along the central ATF. Given their great breadth (100 s of km perpendicular to the ATF), near fault locking and seismic cycle variation in velocity cannot explain the slip rate discrepancy (fig. 10; Thatcher, 1990; Shen et al., 2001). Shen et al. (2001) also point out that because there is no major convergence parallel to the ATF on its south side across the Qimen Tagh 1031 Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China 84 86 88 90 92 94 96 98 0 0.5 1 1.5 2 2.5 3 3.5 4 Longitude (°E) A T F p a ra lle l s lip r a te ( cm ) Ge et al. (1992) Shen et al. (2001) Bendick et al. (2000) Meriaux et al. (1999, 2000) Peltzer et al. (1989) Meriaux et al. (2003) Ryerson et al. (2003) This study Ryerson et al. (2003) Meyer et al. (1996) Fig. 10. Inferred Altyn Tagh Fault-parallel slip rates from long-term studies (blue), geodesy (red), and earth- quake geology (black) compiled from the literature and this study. The eastward decrease in slip rate indicated by the long-term studies is consistent with the activity of numerous thrust faults splaying from and to the south of the ATF (e.g., Meyer et al., 1996, 1998; Tapponnier et al., 2001; although not consistent with the geodetic results from Shen et al., 2001). In general, the geodetic results show lower slip rates relative to the long-term studies. Our earthquake geology results span the two other methods in both the time scale of measurements (a few kyr versus decadal and > 10 kyr) and the resulting inferred rates. Mountains at the southern margin of the Qaidam Basin, the slip rate along the ATF must not decrease significantly between 80° and 92°E. One explanation for the discrepancy in recur- rence times from GPS site velocities and the earthquake record and geology is that a major earthquake along the central ATF is overdue. Another intriguing explanation is that the long term loading rate along the ATF has decreased in the past 100 kyr and strain is now accommodat- ed by other structures in Central Asia. For exam- ple, the faults that produced four M 8 earth- quakes in Mongolia in the past century released significant strain in an area of relatively lower inferred long term deformation rate (Baljinnyam et al., 1993). Earthquakes, such as the 1997 Mw 7.6 Manyi earthquake that ruptured 170 km along the Kunlun Fault (Peltzer et al., 1999) and the 2001 Mw 7.8 Kokoxili earthquake (van der Woerd et al., 2002; Lin et al., 2002) may indi- cate that the western portion of the Kunlun Fault and other central and north Tibetan strike-slip faults may be more active than previously thought. Finally, shortening within the Qimen Tagh and other contractional structures in Northern Tibet may now accommodate a higher percentage of India-Eurasia convergence (Yin and Harrison, 2000; Cowgill, personal commu- nication). 7. Conclusions In this paper, we have characterized the fault trace geometry and paleoseismology of the Central Altyn Tagh Fault, emphasizing its 230 km long and remarkably straight Xorxol geometric segment. Trace geometry, timing of paleoearthquakes, and inferences of rupture magnitudes contribute to the global catalogue of strike-slip fault systems and provide informa- tion about its geodynamic role in the India- Eurasian collision. The trace geometry of the central ATF is dominated by the Arkateng Restraining Bend, Pingding Releasing Step, and the notably straight Xorxol Segment. The geometric bound- aries are probable earthquake rupture bound- aries based on surface rupture continuity, distri- bution of small geomorphic offsets, and histor- ical data. Major Holocene deformation along the Altyn Tagh Fault System is confined to a narrow zone less than 3 km wide. The geometric boundaries and timing con- straints probably preclude one large earthquake rupturing all 3 segments (270 km). Further- more, multiple earthquakes with shorter rupture lengths ( 50 km) rather than complete rupture of the Xorxol Segment better explain the paleo- seismic data. We found 2-3 earthquakes in the last 2-3 kyr. When coupled with typical a- mounts of slip per event (5-10 m), the recur- rence times are tentatively consistent with 1-2 cm/yr slip rates. This result favors models that consider the broader distribution of collisional deformation, rather than those with northward motion of India into Asia absorbed along a few faults bounding rigid blocks. Acknowledgements We thank P. Molnar, R. Bilham, P. Tap- ponnier, M. Hamburger, R. Bürgmann, A. Yin, S. Graham, and J. Freymueller for early reviews of this research, our colleagues at the Institute of Geomechanics in Beijing, and S. Holloway, M. Baillie, S. Selkirk, E. Catlos, C. Schrein, and S. Zhang for help in the lab, field, and the office. Constructive reviews by K. Okamura and F. Galadini and the guidance of Daniela Pantosti for Annals of Geophysics were also very helpful. R. Ryerson arranged our 14C dat- ing at Lawrence Livermore National Labora- tory/University of California Center for Acce- lerator Mass Spectrometry where we were guided by G. Sietz, M. Caffee, J. Southon, and their colleagues. G. Hilley provided impor- tant advice on interpretation of 14C data. The Northern California Earthquake Data Center (NCEDC), and the member networks of the Council of the National Seismic System (CNSS) provided the seismicity data shown in fig. 2. We are particularly grateful to An Yin for inviting us to join the project and to Eric Cowgill for introducing us to field work in Xinjiang. This work is part of a collaborative research effort supported by the U.S. National Science Foundation Continental Dynamics Program (grant EAR-9725780). 1032 Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle REFERENCES AVOUAC, P.J. and P. TAPPONNIER (1993): Kinematic model of active deformation in Central Asia, Geophys. Res. Lett., 20, 895-898. BALJINNYAM, I., A. BAYASGALAN, B.A. BORISOV, A. CISTERNAS, M.G. DEM’YANOVICH, L. GANBAATAR, V.M. KOCHETOV, R.A. KURUSHIN, P. MOLNAR, H. PHILIP and Y.Y. VASHILOV (1993): Ruptures of major earthquakes and active deformation in Mongolia and its surround- ings, Geol. Soc. Am. Mem., 181, pp. 62. BENDICK, R., R. BILHAM, J. FREYMUELLER, K. LARSON, G. PELTZER and G. YIN (2000): Geodetic evidence for a low slip rate on the Altyn Tagh Fault and constraints on the deformation of Asia, Nature, 404, 69-72. CHEN, W-P., C-Y. CHEN and J.L. NÁBELEK (1999): Present- day deformation of the Qaidam basin with implications for intra-continental tectonics, Tectonophysics, 305, 165-181. COWGILL, E., A. YIN, X. WANG and Q. ZHANG (2000): Is the North Altyn Fault part of a strike-slip duplex along the Altyn Tagh Fault System?, Geology, 28, 255-258. COWGILL, E., A. YIN, J.R. ARROWSMITH, W. XIAOFENG and S. ZHANG (2003): The Akato Tagh bend along the Altyn Tagh Fault, NW Tibet, 1. Cenozoic structure, smooth- ing by vertical-axis rotation, and the effect of topo- graphic stresses on borderland faulting, Geol. Soc. Am. (in review). DEPOLO, C.M., D. CLARK, B. SLEMMONS and A. RAMELLI (1991): Historical surface faulting in the Western Basin and Range Province, Western North America: impli- cations for fault segmentation, J. Struct. Geol., 13, 122-136. ENGLAND, P.C. and P. MOLNAR (1997): The field of crustal velocity in Asia calculated from Quaternary rates of slip on faults, Geophys. J. Int., 130, 551-582. FORMAN, S.L. (1999): Infrared and red stimulated lumines- cence dating of Late Quaternary nearshore sediments from Spitsbergen, Svalbard, Arctic, Antarct. Alpine Res., 31, 34-49. GE, S. and 34 OTHERS (1992): Active Altun Fault Zone Monograph (State Seismological Bureau, Seismological Press, Beijing, China), 1-319. GRANT, L.B. and K. SIEH (1994): Paleoseismic evidence of clustered earthquakes on the San Andreas Fault in the Carrizo Plain, California, J. Geophys. Res., 99, 6819- 6841. HARRIS, R.A. and S.M. DAY (1993): Dynamics and fault interaction: parallel strike-slip faults, J. Geophys. Res., 98, 4461-4472. HARRIS, R.A. and S.M. DAY (1999): Dynamic 3D simula- tions of earthquakes on en échelon faults, Geophys. Res. Lett., 206, 2089-2092. HILLEY, G.E. and J.J. YOUNG (2003): Determining even magnitude and timing from paleoseismic and geomor- phic data along the Central San Andreas Fault, California, Bull. Seismol. Soc. Am. (in review). HOLT, W.E., N. CHAMOT-ROOKE, X. LE PICHON, A.J. HAINES, B. SHEN-TU and J. REN (2000): Velocity field in Asia inferred from Quaternary fault slip rates and Global Positioning System observations, J. Geophys. Res., 105, 19,185-19,209. KLINGER, Y., J. VAN DER WOERD, P. TAPPONNIER, X.W. XU, G. KING, W.B. CHEN, W.T. MA and D. BOWMAN (2003): Detailed strip map of the Kokoxili earthquake rupture (MW 7.8, 14/11/01) from space, Geophys. Res. Abstr., 5, 08487. KNUEPFER, P.L.K. (1989): Implications of the characteris- tics of end-points of historical surface fault ruptures for the nature of fault segmentation, in Fault Segmentation and Controls of Rupture Initiation and Termination, edited by D. SCHWARTZ and R. SIBSON, U.S. Geol. Surv. Open-File Rep. 89-315, 193-228. LIN, A., B. FU, J. GUO, Q. ZENG, G. DANG, W. HE and Y. ZHAO (2002): Co-seismic strike-slip and rupture length produce by the 2001 Ms 8.1 Central Kunlun earth- quake, Science, 296, 2015-2017. MÉRIAUX, A.S., P. TAPPONNIER, F.J. RYERSON, J. VAN DER WOERD, C. LASSERRE, X. XU, R. FINKEL and M. CAFFEE (1999): Large-scale strain patterns, great earthquake breaks, and Late Pleistocene slip-rate along the Altyn Tagh Fault (China), European Union of Geosciences Conference Abstracts, EUG 10, Terra Abstracts, 11, p. 55. MÉRIAUX, A.S., F.J. RYERSON, P. TAPPONNIER, J. VAN DER WOERD, R.C. FINKEL, M.W. CAFFEE, C. LASSERRE, X. XU, H. LI AND Z. XU (2000): Fast extrusion of the Tibet Plateau: a 3 cm/yr, 100 kyr slip-rate on the Altyn Tagh Fault, Eos Trans. Am. Geophys. Un., 81 (48), p. 1137. MÉRIAUX, A.S., P. TAPPONNIER, F.J. RYERSON, X. XU, J. VAN DER WOERD, G. KING, R.C. FINKEL, LI HAIBING, M.W. CAFFEE and X. ZHIQIN (2003): Post-glacial slip rate on the Aksay segment of the Northern Altyn Tagh Fault, derived from cosmogenic radionuclide dating of morphological offset features, Geophys. Res. Abstr., 5, 08062. MEYER, B., P. TAPPONNIER, Y. GAUDEMER, G. PELTZER, G. SHUNMIN and C. ZHITAI (1996): Rate of left-lateral movement along the easternmost segment of the Altyn Tagh Fault, east of 96°E (China), Geophys. J. Int., 124, 29-44. MEYER, B., P. TAPPONNIER, L. BOURJOT, F. MÉTEVIER, Y. GAUDEMER, G. PELTZER, G. SHUNMIN and C. ZHITAI (1998): Crustal thickening in Gansu-Qinghai, lithos- pheric mantle subduction, and oblique, strike-slip con- trolled growth of the Tibet plateau, Geophys. J. Int., 135, 1-47. MOLNAR, P., C. BURCHFIEL, L. K’UANGYI and Z. ZIYUN (1987): Geomorphic evidence for active faulting in the Altyn Tagh and Northern Tibet and qualitative esti- mates of its contribution to the convergence of India and Eurasia, Geology, 15, 249-253. PELTZER, G., P. TAPPONNIER and R. ARMIJO (1989): Magnitude of Late Quaternary left-lateral displace- ments along the north edge of Tibet, Science, 246, 1285-1289. PELTZER, G., F. CRAMPÉ and G. KING (1999): Evidence of nonlinear elasticity of the crust from the Mw 7.6 Manyi (Tibet) earthquake, Science, 286, 272-276. RAMSEY, C.B. (2000): OxCal v3.5 software that uses strati- graphic constraints to trim calendar age probability dis- tributions, http://www.rlaha.ox.ac.uk/oxcal/oxcal.htm. RYERSON F.J., A-S. MÉRIAUX, P. TAPPONNIER, J. VAN DER WOERD, XU XIWEI, LI HAIBING and XU ZHIQIN (2003): Northeastwards decrease in the Late Pleistocene- Holocene slip-rate and propagation of the Altyn Tagh Fault (China), Geophys. Res. Abs., 5, 08070. 1033 Paleoseismology of the Xorxol Segment of the Central Altyn Tagh Fault, Xinjiang, China SHEN, Z-K., M. WANG, Y. LI, D.D. JACKSON, A. YIN, D. DONG and P. FANG (2001): Crustal deformation along the Altyn Tagh Fault System, Western China, from GPS, J. Geophys. Res., 106, 30,607-30,621. STUIVER, M. and H.A. POLACH (1977): Discussion, report- ing of 14C data, Radiocarbon, 19, 355-363. STUIVER, M., P.J. REIMER and E. BARD (1998): INTCAL98 radiocarbon age calibration, 24.000-0 cal BP, Radio- carbon, 40, 1041-1081. TAPPONNIER, P. and P. MOLNAR (1977): Active faulting and tectonics in China, J. Geophys. Res., 82, 2905- 2930. TAPPONNIER, P., X. ZHIQIN, F. ROGER, B. MEYER, N. ARNAUD, G. WITTLINGER and Y. JINGSUI (2001): Oblique stepwise rise and growth of the Tibet plateau, Science, 294, 1671-1677. THATCHER, W. (1990): Present-day crustal movements and the mechanics of cyclic deformation, in The San Andreas Fault System, California, edited by R.E. WALLACE, USGS Prof. Paper 1515, 189-205. VAN DER WOERD, J., F.J. RYERSON, P. TAPPONNIER, Y. GAUDEMER, R. FINKEL, A.S. MERIAUX, M.W. CAFFEE, Z. GUOGUANG and H. QUNLU (1998): Holocene left-slip rate determined by cosmogenic surface dating on the Xidatan Segment of the Kunlun Fault (Qinghai, China), Geology, 26, 695-698. VAN DER WOERD, J., X. XU, H. LI, P. TAPPONIER, B. MEYER, F.J. RYERSON, A.S. MÉRIAUX and X. ZHIQIN (2001): Rapid active thrusting along the northwestern range front of the Tanghe Nan Shan (Western Gansu, China), J. Geophys. Res., 106, 30,475-30,504. VAN DER WOERD, J., A. MÉRIAUX, Y. KLINGER, F.J. RYERSON, Y. GAUDEMER and P. TAPPONNIER (2002): The 14 No- vember 2001, Mw = 7.8 Kokoxili earthquake in Nort- hern Tibet (Quinghai Province, China), Seismol. Res. Lett., 73, 125-135. WANG, Q., P-Z. ZHANG, J.T. FREYMULLER, R. BILHAM, K.M. LARSON, X. LAI, X. YOU, Z. NIU, J. WU, Y. LI, J. LIU, Z. YANG and Q. CHEN (2001): Present-day crustal defor- mation in China constrained by Global Positioning System measurements, Science, 294, 574-577. WARD, S.N. (1997): Dogtails versus rainbows: synthetic earth- quake rupture models as an aid in interpreting geological data, Bull. Seismol. Soc. Am., 87, 1422-1441. WASHBURN, Z. (2001): Quaternary tectonics and earthquake geology of the Central Altyn Tagh Fault, Xinjiang, China: implications for tectonic process along the northern margin of Tibet, M.S. Thesis, Tempe, Arizona State University, pp. 179. WASHBURN, Z., J.R. ARROWSMITH, S. FORMAN, E. COWGILL, X.F. WANG, Y. ZHANG and C. ZHENGLE (2001): Late Holocene earthquake history of the Central Altyn Tagh Fault, China, Geology, 29, 1051-1054. WELLS, D.L. and K. COPPERSMITH (1994): New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement, Seismol. Soc. Am. Bull., 84, 974-1002. WESNOUSKY, S.G. (1989): Seismicity and the structural evo- lution of strike-slip faults, in Fault Segmentation and Controls of Rupture Initiation and Termination, edited by D. SCHWARTZ and R. SIBSON, U.S. Geol. Surv. Open- File Rep. 89-315, 409-431. WESSEL, P. and W.H.F. SMITH (1995): New version of the Generic Mapping Tools released, Eos, Trans. Am. Geo- phys. Un., 76 (33), p. 329. YIN, A. and M. HARRISON (2000): Geologic evolution of the Himalayan-Tibetan orogen, Annu. Rev. Earth Planet. Sci., 28, 211-280. 1034 Zack Washburn, J Ramón Arrowsmith, Guillaume Dupont-Nivet, Wang Xiao Feng, Zhang Yu Qiao and Chen Zhengle