2014 GSA Annual Meeting in Vancouver, British Columbia (19–22 October 2014)

Paper No. 181-8
Presentation Time: 10:10 AM

TESTING TECTONIC MODELS USING TIME-VARYING CRUSTAL THICKNESS VARIATIONS ACROSS SOUTHERN TIBET


HARRISON, T. Mark1, WIELICKI, Matthew2, LOVERA, Oscar M.3, DEPAOLO, D.J.4, ZHU, Di-Cheng5, ZHAO, Zhidan5, YIN, An6 and MO, Xuanxue7, (1)Department of Earth, Planetary and Space Sciences, UCLA, Los Angeles, CA 90095-1567, (2)Department of Geological Sciences, University of Alabama, 201 7th Avenue, Tuscaloosa, AL 35406, (3)Dept. of Earth, Planetary and Space Sciences, Univ. of California, Los Angeles, CA 90095-1567, (4)Earth and Planetary Science, Univ California - Berkeley, 301 McCone Hall, Berkeley, CA 94720-4767, (5)School of Earth Science and Mineral Resources, China University of Geosciences, 29 Xueyuan Road, Beijing, 100083, China, (6)Department of Earth, Planetary, and Space Sciences, University of California, Los Angeles, Los Angeles, CA 90095, (7)Earth Science and Resources, China University of Geosciences, Beijing, 29 Xueyuan Road, Beijing, 100083, China

More than 30 years after the classic tectonic models for the Tibetan-Himalayan orogen were formulated, no clear test has been devised to select between a pure-shear thickening history and one involving pre-collisional topography and crustal extrusion. That these models imply vastly different crustal thickness histories underscores how little we know of this parameter. From isotopic data, we previously estimated a crustal thickness profile at the initiation of collision (~50 Ma) ranging from 20 km at the ITS to 45 km thick 130 km to the N. Since S. Tibet has a uniform ~85 km thickness today (inc. ~30 km of Indian crust) and experienced little post-50 Ma upper crustal deformation, we inferred that ~20 km was added by tectonic accretion and ~15 km from juvenile magma since collision. To assess pre-collisional variations, we undertook a traverse across S. Tibet at 92oE. Zircons from 82 granitoids (28 to 230 Ma) were analyzed for U-Pb, δ18O, εHf, and trace elements; results show the high degree of variability expected from RAFC processes. Zircon εHf data show a similar distribution to whole rock εNd permitting its use as a crustal thickness proxy. Zircons from 195 to 50 Ma granitoids near the ITS all show +ve εHf whereas the northernmost samples (100 to 230 Ma) range from -7 to -22. Although RAFC modeling is often limited by its multivariate nature relative to available geochemical controls, zircon Txlln data permit new insights into these processes. We developed an RAFC model featuring a parameterized T-crystallization relationship coupled to a zircon growth model that predicts zircon εHf as a function of Txlln. The model shows that recharge yields a zircon population that records the full spectrum of εHf in the system whereas no recharge yields a much narrower distribution. Thus zircons may have the potential to preserve a record of the magnitude of recharge. Insights gained from modelling reinforce our view that the relationship between assimilation and geothermal structure can be used to estimate past crustal thickness of convergent margins. Interpreted through this model, our data show that the southern margin of Asia remained thin (<25 km) until collision began whereas the northern terrane was continuously thick from 200 Ma. This observation reinforces the growing view that Tibet hosted significant topography prior to collision.