JOURNAL OF GEOPHYSICAL RESEARCH Supporting Information for Interseismic Coupling on the Main Frontal Thrust DOI:./ V. L. Stevens and JP Avouac, 4 5 6 7 8 9 Contents of this file. Tables S to S.. Figures S to S7. Additional Supporting Information (Files uploaded separately). Dataset S - GPS Data used in the Inversion.. Dataset S - Electronic data for coupling map of Figure.. Dataset S - Electronic data for stress rate map of Figure b. Corresponding author: V. L. Stevens, Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 95, USA. (vstevens@caltech.edu) Geological and Planetary sciences, California Institute of Technology, Pasadena, California, USA. Bullard Laboratories, Cambridge University, Cambridge, Cambridgeshire UK.
X -V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST Introduction Here we show extra figures relating to data, model set-up, resolution and residuals. References 4 5 6 7 8 9 4 5 7 9 Ader, T., J.-P. Avouac, J. Liu-Zeng, H. Lyon-Caen, L. Bollinger, J. Galetzka, J. Genrich, M. Thomas, K. Chanard, S. N. Sapkota, S. Rajaure, P. Shrestha, L. Ding, and M. Flouzat (), Convergence rate across the nepal himalaya and interseismic coupling on the main himalayan thrust: Implications for seismic hazard, J. Geophys. Res., 7 (B4). Berthet, T., J.-F. Ritz, M. Ferry, P. Pelgay, R. Cattin, D. Drukpa, R. Braucher, and G. Hetnyi (4), Active tectonics of the eastern himalaya: New constraints from the first tectonic geomorphology study in southern bhutan, Geology, 4 (5), 47 4, doi:./g56.. Bollinger, L., S. N. Sapkota, P. Tapponnier, Y. Klinger, M. Rizza, J. Van der Woerd, D. R. Tiwari, R. Pandey, A. Bitri, and S. Bes de Berc (4), Estimating the return times of great himalayan earthquakes in eastern nepal: Evidence from the patu and bardibas strands of the main frontal thrust, Journal of Geophysical Research: Solid Earth, 9 (9), 7 76, doi:./4jb97. Burgess, W. P., A. Yin, C. S. Dubey, Z.-K. Shen, and T. K. Kelty (), Holocene shortening across the main frontal thrust zone in the eastern himalaya, Earth and Planetary Science Letters, 5758 (), 5 67, doi:http://dx.doi.org/.6/j.epsl..9.4. Colchen, M. (999), The thakkholamustang graben in nepal and the late cenozoic extension in the higher himalayas, Journal of Asian Earth Sciences, 7 (56), 68 7,
V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST X - 5 6 7 8 9 4 4 4 4 44 45 46 47 48 49 5 5 5 5 doi:http://dx.doi.org/.6/s67-9(99)7-. Gan, W., P. Zhang, Z.-K. Shen, Z. Niu, M. Wang, Y. Wan, D. Zhou, and J. Cheng (7), Present-day crustal motion within the tibetan plateau inferred from gps measurements, J. Geophys. Res., (B8), B8,46. Kundu, B., R. K. Yadav, B. S. Bali, S. Chowdhury, and V. K. Gahalaut (4), Oblique convergence and slip partitioning in the nw himalaya: Implications from gps measurements, Tectonics, (), 4, doi:./4tc6. Lav, J., and J. P. Avouac (), Active folding of fluvial terraces across the siwaliks hills, himalayas of central nepal, Journal of Geophysical Research: Solid Earth, 5 (B), 575 577, doi:.9/999jb99. Mugnier, J., P. Huyghe, P. Leturmy, and F. Jouanne (), Episodicity and rates of thrust-sheet motion in the himalayas (western nepal), AAPG Mem., 8, 4. Parkash, B., R. Rathor, P. Pati, R. Jakhmola, and S. S. (), Convergence rates along the himalayan frontal thrust inferred from terraces at chandidevi temple hill, hardwar, northwestern himalaya, Curr Sci,, 4 4. Styron, R., M. Taylor, and K. Okoronkwo (), Database of active structures from the indo-asian collision, Eos, Transactions American Geophysical Union, 9 (), 8 8, doi:.9/eo. Thakur, V., M. Joshi, D. Sahoo, N. Suresh, R. Jayangondapermal, and A. Singh (4), Partitioning of convergence in northwest sub-himalaya: estimation of late quaternary uplift and convergence rates across the kangra reentrant, north india, International Journal of Earth Sciences, (4), 7 56, doi:.7/s5-4-6-7.
X - 4 V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST 54 55 56 57 58 59 6 6 Vernant, P., R. Bilham, W. Szeliga, D. Drupka, S. Kalita, A. K. Bhattacharyya, V. K. Gaur, P. Pelgay, R. Cattin, and T. Berthet (4), Clockwise rotation of the brahmaputra valley relative to india: Tectonic convergence in the eastern himalaya, naga hills, and shillong plateau, Journal of Geophysical Research: Solid Earth, 9 (8), 6558 657, doi:./4jb96. Wesnousky, S. G., S. Kumar, R. Mohindra, and V. C. Thakur (999), Uplift and convergence along the himalayan frontal thrust of india, Tectonics, 8 (6), 967 976, doi:.9/999tc9.
V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST X - 5 Table S. Long term slip rates from geomorphic studies and modelled here Lon Slip Rate (cm/yr) Uncertainty Modelled Modelled Reference Slip Rate (cm/yr) Uncertainty 76...8 4.6.5 Thakur et al. [4] 77.5.9. 7..7 Wesnousky et al. [999] 78..8. 7..7 Parkash et al. [] 8. 9 6. 9.. Mugnier et al. [] 85..5 9.7.6 Lav and Avouac [] 85.8 8. 6 9.7.6 Bollinger et al. [4] 9..8 8.8 7.4. Berthet et al. [4] 9.5 6... Burgess et al. []
X - 6 V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST Table S. Rotation Poles with respect to India for the Sikkim/Shillong and Bhutan/Assam Blocks of this study, and the Shillong and Assam Blocks in the study of Vernant et al. [4]. Block Lon Lat Rotation Rate Study Assam 88..6 -. this study Shillong 87. 5.6 -.6 this study Assam 87.8.8 -. Vernant et al. [4] Shillong 88.8.4 -.5 Vernant et al. [4] Latitude ( N) 9 Continuous GPS mm/yr 7 Campaign GPS mm/yr InSAR data -,, + mm LOS/yr Leveling data,, 6 mm uplift/yr 5 74 76 78 8 8 84 86 88 9 9 94 96 Longitude ( E) Figure S. Location of data used in the inversion. Arrows show velocities relative to India (as defined by Ader et al. []) determined from continuous (green) and campaign (black) GPS data. Ellipses show the uncertainties at the 67% confidence level. The blue circles are InSAR data. The magenta circles are leveling data. The thick red line is the simplified geometry of the MHT used in the model. Country borders are marked in grey. July 7, 5, :6pm
Latitude ( N) V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST X - 7 9 SHILLONG ASSAM KC 7 TH DS YG 5 74 76 78 8 8 84 86 88 9 9 94 96 Longitude ( E) Figure S. The model setup. The red line shows the simplified surface trace of the fault, and the grid shows the discretization used in this model. The dashed black lines are the boundaries between different regions of uniform long term slip rate. The blue lines show the edges of the two blocks for which rotation poles are found. The red dashed lines show the rough location of major grabens. KC = Kaurik Changdu rift [Kundu et al., 4]. TH = Thakkola graben [Colchen, 999]. DS = Dingjie-Sehnza fault zone [Gan et al., 7]. YG = Yadong-Gulu rift [Gan et al., 7]. Thinner red/brown lines show faults from Styron et al. [].
X - 8 V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST. ±.7 Coupling.5 8.5 ±.8. ±..9 ±.9 Latitude (ºN) 9.4 ±.4 7.6 ±.9 Continuous GPS Campaign GPS Model 74 Figure S. 76 78 8 8 ASSAM SHILLONG INDIA mm/yr 84 86 Longitude (ºE) 88 9 9 94 Same as Figure, except with the poles of Vernant et al. [4] imposed and not solved for. Coupling.5 Latitude (ºN) 74 Figure S4. 76 78 8 8 84 86 Longitude (ºE) 88 9 9 94 96 Same as Figure a, except with uniform weighting of the Laplacian. The coupling model is thus not regulated by the seismicity. July 7, 5, :6pm
V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST X - 9 9 mm Latitude (ºN) 84 86 - - - Continuous GPS Campaign GPS mm/yr 74 76 78 8 8 84 86 88 9 9 94 96 Longitude (ºE) Figure S5. Residuals to the GPS data, with uncertainy ellipses. Campaign measurements are in red, and continuous measurements in green. The inset shows the residuals of the InSAR and leveling data.
X - V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST 5 5 Resolution, km 4 6 8 Latitude (ºN) Continuous GPS station Campaign GPS station Leveling point InSAR point 74 76 78 8 8 84 86 Longitude (ºE) 88 9 9 94 96 Figure S6. Resolution on each patch of the fault. The resolution here is the characteristic size of the smallest inhomogeneities of coupling which could in principle be resolved, given the spatial distribution and uncertainties of the measurements. The figure is saturated at a resolution of 8 km, as above this value we assume there is no resolution on the corresponding patch, and so slip on this patch is just determined by slip on the neighbouring patches. See Ader et al. [] for details of the calculation. The location of data points used to find the resolution are also shown. July 7, 5, :6pm
Elevation, km Modeled Uplift Rate, mm/yr Elevation, km Modeled Uplift Rate, mm/yr Elevation, km Modeled Uplift Rate, mm/yr Elevation, km Modeled Uplift Rate, mm/yr V L STEVENS AND JP AVOUAC: INTERSEISMIC COUPLING ON THE MAIN HIMALAYAN THRUST X - 6 75 E 6 8 E 5 5 4 4 -.5.5.5 Latitude ( N) -.5 9 9.5.5 Latitude ( N) 6 85 E 6 9 E 5 5 4 4 -.5 7 7.5.5 9 Latitude ( N) -.5 7 7.5.5 9 Latitude ( N) Figure S7. Smoothed elevation (red) compared with modelled uplift rates (blue) at different longitudes. Profile locations are shown in Figure 5 of the main text.