Counterclockwise exhumation of a hot orogen: The Paleoproterozoic ultrahigh-temperature granulites in the North China Craton

Regional metamorphic belts provide important constraints on the plate tectonic architecture of orogens. We report here a detailed petrologic examination of the sapphirine-bearing ultra-high temperature (UHT) granulites from the Jining Complex within the Khondalite Belt of the North China Craton (NCC). These granulites carry diagnostic UHT assemblages and their microstructures provide robust evidence to trace the prograde, peak and retrograde metamorphic evolution. The P–T conditions of the granulites estimated from XMgGrt(Mg/Fe + Mg) − XMgSpr isopleth calculations indicate temperature above 970 °C and pressures close to 7 kbar. We present phase diagrams based on thermodynamic computations to evaluate the mineral assemblages and microstructures and trace the metamorphic trajectory of the rocks. The evolution from Spl–Qtz–Ilm–Crd–Grt–Sil to Spr–Qtz–Crd–Opx–Ilm marks the prograde stage. The Spl–Qtz assemblage appears on the low-pressure side of the P–T space with Spr–Qtz stable at the high-pressure side, possibly representing an increase in pressure corresponding to compression. The spectacular development of sapphirine rims around spinel enclosed in quartz supports this inference. An evaluation of the key UHT assemblages based on model proportion calculation suggests a counterclockwise P–T path. With few exceptions, granulite-facies rocks developed along collisional metamorphic zones have generally been characterized by clockwise exhumation trajectories. Recent evaluation of the P–T paths of metamorphic rocks developed within collisional orogens indicates that in many cases the exhumation trajectories follow the model subduction geotherm, in accordance with a tectonic model in which the metamorphic rocks are subducted and exhumed along a plate boundary. The timing of UHT metamorphism in the NCC (c. 1.92 Ga) coincides with the assembly of the NCC within the Paleoproterozoic Columbia supercontinent, a process that would have involved subduction of passive margins sediments and closure of the intervening ocean. Thus, the counterclockwise P–T path obtained in this study correlates well with a tectonic model involving subduction and final collisional suturing, with the UHT granulites representing the core of the hot or ultra-hot orogen developed during Columbia amalgamation.

Retrogressed eclogite (20 kbar, 1020 °C) from the Neoproterozoic Palghat–Cauvery suture zone, southern India

The Palghat–Cauvery suture zone in southern India separates Archaean crustal blocks to the north and the Proterozoic Madurai block to the south. Here we present the first detailed study of a partially retrogressed eclogite (from within the Sittampundi anorthositic complex in the suture zone) that occurs as a 20-cm wide layer in a garnet gabbro layer in anorthosite. The eclogite largely consists of an assemblage of coexisting porphyroblasts of almandine–pyrope garnet and augitic clinopyroxene. However, a few garnets contain inclusions of omphacite. Rims and symplectites composed of Na–Ca amphibole and plagioclase form a retrograde assemblage. Petrographic analysis and calculated phase equilibria indicate that garnet–omphacite–rutile–melt was the peak metamorphic assemblage and that it formed at ca. 20 kbar and above 1000 °C. The eclogite was exhumed on a very tight hairpin-type, anticlockwise P–T path, which we relate to subduction and exhumation in the Palghat–Cauvery suture zone. The REE composition of the minerals suggests a basaltic oceanic crustal protolith metamorphosed in a subduction regime. Geological–structural relations combined with geophysical data from the Palghat–Cauvery suture zone suggest that the eclogite facies metamorphism was related to formation of the suture zone. Closure of the Mozambique Ocean led to development of the suture zone and to its western extension in the Betsimisaraka suture of Madagascar.

Neoarchean sanukitoid series in the Karelian Province, Finland

Sanukitoid series intrusions can be found throughout the Archean Karelian Province of the Fennoscandian shield. All sanukitoids share the same controversial elemental characteristics: they have high content of incompatible elements such as K, Ba, and Sr as well as high content of the compatible elements Mg, Cr, and Ni, and high Mg#. This composition is explained by an enriched mantle wedge origin in a Neoarchean subduction setting. This study concentrates on sanukitoid intrusions and tonalite-trondhjemite-granodiorite series (TTGs) from Finnish part of the Karelian Province. The collected rock samples have been studied in the field and under microscope as well as for their whole-rock (including isotopes) and mineral compositions. The new data together with previously published analyses help us to better understand the petrogenesis, tectonic setting and reworking of the Archean rock units. TTGs from the Karelian Province form a voluminous series of granitoids and reworked migmatites. This study divides TTG series into two subgroups based on their elemental composition: low-HREE (heavy rare earth element) TTGs and high-HREE TTGs indicating pressure differences in their source. Sanukitoid series is a minor, divergent group of intrusions. These intrusions are variable sized, and the texture varies from even-grained to K-feldspar porphyritic. The elemental composition differentiates sanukitoids from more voluminous TTG groups, the SiO2 in sanukitoids varies to include series of gabbro, diorite, and granodiorite. U Pb age determinations from sanukitoid series show temporally limited emplacement between ~ 2745 2715 Ma after the main crust forming period in the area. Hafnium, neodymium, common lead, and oxygene isotopes indicate well homogenized characteristics. Recycled crust has made a variable, yet minor, contribution to sanukitoids, as evidenced by oxygene isotopes and inherited zircon cores. A proposed tectonic setting for the formation of the sanukitoid series is slab breakoff of oceanic lithosphere in subduction setting, with sanukitoids deriving from an enriched mantle wedge. The proposed setting explains some of the peculiar features of sanukitoids, such as their temporally limited occurrence and controversial elemental composition. Sanukitoids would occur after cessation of the regional growth of Archean crust, and they could be derived from mantle wedge previously enriched by melts and fluids from oceanic crust and sediments. A subsequent event during the Paleoproterozoic Svecofennian orogeny at ~1.9 Ga affected the appearance and microstructures of the rocks as well as caused redistribution of lead between minerals and whole rock. However, the deformation was not able to obliterate the original geochemical characteristics of these sanukitoids.

Paleoseismological evidence of surface faulting along the northeastern Himalayan front, India: Timing, size, and spatial extent of great earthquakes

The similar to 2500 km long Himalayan arc has experienced three large to great earthquakes of M-w 7.8 to 8.4 during the past century, but none produced surface rupture. Paleoseismic studies have been conducted during the last decade to begin understanding the timing, size, rupture extent, return period, and mechanics of the faulting associated with the occurrence of large surface rupturing earthquakes along the similar to 2500 km long Himalayan Frontal Thrust (HFT) system of India and Nepal. The previous studies have been limited to about nine sites along the western two-thirds of the HFT extending through northwest India and along the southern border of Nepal. We present here the results of paleoseismic investigations at three additional sites further to the northeast along the HFT within the Indian states of West Bengal and Assam. The three sites reside between the meizoseismal areas of the 1934 Bihar-Nepal and 1950 Assam earthquakes. The two westernmost of the sites, near the village of Chalsa and near the Nameri Tiger Preserve, show that offsets during the last surface rupture event were at minimum of about 14 m and 12 m, respectively. Limits on the ages of surface rupture at Chalsa (site A) and Nameri (site B), though broad, allow the possibility that the two sites record the same great historical rupture reported in Nepal around A.D. 1100. The correlation between the two sites is supported by the observation that the large displacements as recorded at Chalsa and Nameri would most likely be associated with rupture lengths of hundreds of kilometers or more and are on the same order as reported for a surface rupture earthquake reported in Nepal around A.D. 1100. Assuming the offsets observed at Chalsa and Nameri occurred synchronously with reported offsets in Nepal, the rupture length of the event would approach 700 to 800 km. The easternmost site is located within Harmutty Tea Estate (site C) at the edges of the 1950 Assam earthquake meizoseismal area. Here the most recent event offset is relatively much smaller (<2.5 m), and radiocarbon dating shows it to have occurred after A.D. 1100 (after about A.D. 1270). The location of the site near the edge of the meizoseismal region of the 1950 Assam earthquake and the relatively lesser offset allows speculation that the displacement records the 1950 M-w 8.4 Assam earthquake. Scatter in radiocarbon ages on detrital charcoal has not resulted in a firm bracket on the timing of events observed in the trenches. Nonetheless, the observations collected here, when taken together, suggest that the largest of thrust earthquakes along the Himalayan arc have rupture lengths and displacements of similar scale to the largest that have occurred historically along the world's subduction zones.

Intraplate Response to the Great 2004 Sumatra-Andaman Earthquake: A Study from the Andaman Segment

The 1300-km rupture of the 2004 interplate earthquake terminated at around 15 degrees N, in the northernmost segment of the Andaman-Nicobar subduction zone. This part of the plate boundary is noted for its generally lower level seismicity, compared with the southern segments. Based on the Global Centroid Moment Tensor (CMT) and National Earthquake Information Center (NEIC) data, most of the earthquakes of M-w >= 4.5 prior to 2004 were associated with the Andaman Spreading Ridge (ASR), and a few events were located within the forearc basin. The 2004 event was followed by an upward migration of hypocenters along the subducting plate, and the Andaman segment experienced a surge of aftershock activity. The continuing extensional faulting events, including the most recent earthquake (10 August 2009; M-w 7.5) in the northern end of the 2004 rupture, suggest the reduction of compressional strain associated with the interplate event. The style of faulting of the intraplate events before and after a great plate boundary earthquake reflects the relative influences of the plate-driving forces. Here we discuss the pattern of earthquakes in the Andaman segment before and after the 2004 event to appraise the spatial and temporal relation between large interplate thrust events and intraplate deformation. This study suggests that faulting mechanisms in the outer-ridge and outer-rise regions could be indicative of the maturity of interplate seismic cycles.

Multiple garnet growth in garnet-kyanite-staurolite gneiss, Pangong metamorphic complex, Ladakh Himalaya: New constraints on tectonic setting

Garnet-kyanite-staurolite gneiss in the Pangong complex, Ladakh Himalaya, contains porphyroblastic euhedral garnets, blades of kyanite and resorbed staurolite surrounded by a fine-grained muscovite-biotite matrix associated with a leucogranite layer. Sillimanite is absent. The gneiss contains two generations of garnet in cores and rims that represent two stages of metamorphism. Garnet cores are extremely rich in Mn (X(Sps) = 0.35-038) and poor in Fe (X(Alm) = 0.40-0.45), whereas rims are relatively Mn-poor (X(Sps) =0.07-0.08), and rich in Fe (X(Alm), = 0.75-0.77). We suggest that garnet cores formed during prograde metamorphism in a subduction zone followed by abrupt exhumation, during early collision of the Ladakh arc and Karakoram block. The subsequent India-Asia continental collision subducted the metamorphic rocks to a mid-crustal level, where the garnet rims overgrew the Mn-rich cores at ca. 680 degrees C and ca. 8.5 kbar. PT calculations were estimated from phase diagrams calculated using a calculated bulk chemical composition in the Mn-NCKFMASHT system for the garnet-kyanite-staurolite-bearing assemblage. Muscovites from the metamorphic rocks and associated leucogranites have consistent K-Ar ages (ca. 10 Ma), closely related to activation of the Karakoram fault in the Pangong metamorphic complex. These ages indicate the contemporaneity of the exhumation of the metamorphic rocks and the cooling of the leucogranites. (C) 2011 Elsevier B.V. All rights reserved.

Constraining large earthquakes along the Andaman trench using deepwater turbidites: prospects and challenges

The 2004 Sumatra-Andaman earthquake was unprecedented in terms of its magnitude (M-w 9.2), rupture length along the plate boundary (1300 km) and size of the resultant tsunami. Since 2004, efforts are being made to improve the understanding of the seismic hazard in the Sumatra-Andaman subduction zone in terms of recurrence patterns of major earthquakes and tsunamis. It is reasonable to assume that previous earthquake events in the Myanmar Andaman segment must be preserved in the geological record in the form of seismo-turbidite sequences. Here we present the prospects of conducting deep ocean palaeoseismicity investigations in order to refine the quantification of the recurrence pattern of large subduction-zone earthquakes along the Andaman-Myanmar arc. Our participation in the Sagar Kanya cruise SK-273 (in June 2010) was to test the efficacy of such a survey. The primary mission of the cruise, along a short length (300 km) of the Sumatra Andaman subduction front was to collect bathymetric data of the ocean floor trenchward of the Andaman Islands. The agenda of our piggyback survey was to fix potential coring sites that might preserve seismo-turbidite deposits. In this article we present the possibilities and challenges of such an exercise and our first-hand experience of such a preliminary survey. This account will help future researchers with similar scientific objectives who would want to survey the deep ocean archives of this region for evidence of extreme events like major earthquakes.

Ages and relative sizes of pre-2004 tsunamis in the Bay of Bengal inferred from geologic evidence in the Andaman and Nicobar Islands

Geologic evidence along the northern part of the 2004 Aceh-Andaman rupture suggests that this region generated as many as five tsunamis in the prior 2000years. We identify this evidence by drawing analogy with geologic records of land-level change and the tsunami in 2004 from the Andaman and Nicobar Islands (A&N). These analogs include subsided mangrove swamps, uplifted coral terraces, liquefaction, and organic soils coated by sand and coral rubble. The pre-2004 evidence varies in potency, and materials dated provide limiting ages on inferred tsunamis. The earliest tsunamis occurred between the second and sixth centuries A.D., evidenced by coral debris of the southern Car Nicobar Island. A subsequent tsunami, probably in the range A.D. 770-1040, is inferred from deposits both in A&N and on the Indian subcontinent. It is the strongest candidate for a 2004-caliber earthquake in the past 2000years. A&N also contain tsunami deposits from A.D. 1250 to 1450 that probably match those previously reported from Sumatra and Thailand, and which likely date to the 1390s or 1450s if correlated with well-dated coral uplift offshore Sumatra. Thus, age data from A&N suggest that within the uncertainties in estimating relative sizes of paleo-earthquakes and tsunamis, the 1000year interval can be divided in half by the earthquake or earthquakes of A.D. 1250-1450 of magnitude >8.0 and consequent tsunamis. Unlike the transoceanic tsunamis generated by full or partial rupture of the subduction interface, the A&N geology further provides evidence for the smaller-sized historical tsunamis of 1762 and 1881, which may have been damaging locally.

Neoarchean continental growth through arc magmatism in the Nilgiri Block, southern India

The formation and growth of continental crust in the Archean have been evaluated through models of subduction-accretion and mantle plume. The Nilgiri Block in southern India exposes exhumed Neoarchean lower crust, uplifted to heights of 2500 m above sea level along the north western margin of the Peninsula. Major lithologies in this block include charnockite with or without garnet, anorthosite-gabbro suite, pyroxenite, amphibolite and hornblende-biotite gneiss (TTG). All these rock types are closely associated as an arc magmatic suite, with diffuse boundaries and coeval nature. The charnockite and hornblende-biotite gneisses (TTG) show SiO2 content varying from 64 to 73 wt.%. The hornblende-biotite gneisses (TTG) are high-Al type with Al2O3 >15 wt.% whereas the charnockites show Al2O3 <15 wt.%. The composition of charnockite is mainly magnesian and calcic to calc-alkaline. The mafic-ultramafic rocks show composition close to that of tholeiitic series. The low values of K(2)o (<3 wt.%), (K/Rb)/K2O (<500), Zr/Ti, and trace element ratios like (La/Yb)n/(Sr/Y), (Y/Nb), (Y + Nb)/Rb, (Y+Ta)/Rb, Yb/Ta indicate a volcanic arc signature for these rocks. The geochemical signature is consistent with arc magmatic rocks generated through oceanic plate subduction. The primitive mantle normalized trace element patterns of these rocks display enrichment in large ion lithophile elements (LILE) and comparable high field strength elements (HFSE) in charnockite and hornblende-biotite gneisses (TTG) consistent with subduction-related origin. Primitive mantle normalized REE pattern displays an enrichment in LREE in the chamockite and hornblende-biotite gneisses (TTG) as compared to a flat pattern for the mafic rocks. The chondrite normalized REE patterns of zircons of all the rock types reveal cores with high HREE formed at ca. 2700 Ma and rims with low HREE formed at 2500-2450 Ma. Log-transformed La/Th-Nb/Th-Sm/Th-Yb/Th discrimination diagram for the mafic and ultramafic rocks from Nilgiri displays a transition from mid-oceanic ridge basalt (MORB) to island arc basalt (IAB) suggesting a MORB source. The U-Pb zircon data from the charnockites, mafic granulites and hornblende-biotite gneisses (TTG) presented in our study show that the magma generation during subduction and accretion events in this block occurred at 2700-2500 Ma. Together with the recent report on Neoarchean supra-subduction zone ophiolite suite at its southern margin, the Nilgiri Block provides one of the best examples for continental growth through vertical stacking and lateral accretion in a subduction environment during the Neoarchean. (c) 2014 Elsevier B.V. All rights reserved.

Evolution of high-pressure mafic granulites and pelitic gneisses from NE Madagascar: Tectonic implications

The occurrence of high-pressure mafic-ultramafic bodies within major shear zones is one of the indicators of paleo-subduction. In mafic granulites of the Andriamena complex (north-eastern Madagascar) we document unusual textures including garnet-clinopyroxene-quartz coronas that formed after the breakdown of orthopyroxene-plagioclase-ilmenite. Textural evidence and isochemical phase diagram calculations in the Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O-TiO2 system indicate a pressure-temperature (P-T) evolution from an isothermal (780 degrees C) pressure up to c. 24 kbar to decompression and cooling. Such a P-T trajectory is typically attained in a subduction zone setting where a gabbroic/ultramafic complex is subducted and later exhumed to the present crustal level during oceanic closure and final continental collision. The present results suggest that the presence of such deeply subducted rocks of the Andriamena complex is related to formation of the Betsimisaraka suture. LA-ICPMS U-Pb zircon dating of pelitic gneisses from the Betsimisaraka suture yields low Th/U ratios and protolith ages ranging from 2535 to 2625 Ma. A granitic gneiss from the Alaotra complex yields a zircon crystallization age of ca. 818 Ma and Th/U ratios vary from 1.08 to 2.09. K-Ar dating of muscovite and biotite from biotite-kyanite-sillimanite gneiss and garnet-biotite gneiss yields age of 486 +/- 9 Ma and 459 +/- 9 Ma respectively. We have estimated regional crustal thicknesses in NE Madagascar using a flexural inversion technique, which indicates the presence of an anomalously thick crust (c. 43 km) beneath the Antananarivo block. This result is consistent with the present concept that subduction beneath the Antananarivo block resulted in a more competent and thicker crust. The textural data, thermodynamic model, and geophysical evidence together provide a new insight to the subduction history, crustal thickening and evolution of the high-pressure Andriamena complex and its link to the terminal formation of the Betsimisaraka suture in north-eastern Madagascar. (C) 2015 Elsevier B.V. All rights reserved.

Archean tectonics and crustal evolution of the Biligiri Rangan Block, southern India

The Southern Granulite Terrain in India is a collage of crustal blocks ranging in age from Archean to Neoproterozoic. This study investigate the tectonic evolution of one of the northernmost block- the Biligiri Block (BRB) through a multidisciplinary approach involving field investigation, petrographic studies, LA-ICPMS zircon U-Pb geochronology, Hf isotopic analyses, metamorphic P-T phase diagram computations, and crustal thickness modeling. The garnet bearing quartzofeldspathic gneiss from the central BRB preserve Mesoarchean magmatic zircons with ages between 3207 and 2806 Ma and positive epsilon Hf value (+2.7) which possibly indicates vestiges of a Mesoarchean primitive continental crust. The occurrence of quartzite-iron formation intercalation as well as ultramafic lenses along the western boundary of the BRB is interpreted to indicate that the Kollegal structural lineament is a possible paleo-suture. Phase diagram computation of a metagabbro from the southwestern periphery of the Kollegal suture zone reveals high-pressure (similar to 18.5 kbar) and medium-temperature (similar to 840 degrees C) metamorphism, likely during eastward subduction of the Western Dharwar oceanic crust beneath the Mesoarchean BRB. In the model presented here, slab subduction, melting and underplating processes generated arc magmatism and subsequent charnockitization within the BRB between ca. 2650 Ma and ca. 2498 Ma. These results thus reveal Meso- to Neoarchean tectonic evolution of the BRB. The spatial variation of crustal thickness, derived from flexure inversion technique, provides additional constraints on the tectonic linkage of the BRB with its surrounding terrains. In conjunction with published data, the Moyar and the Kollegal suture zones are considered to mark the trace of ocean closure along which the Nilgiri and Biligiri Rangan Blocks accreted on to the Western Dharwar Craton. (C) 2016 Elsevier B.V. All rights reserved.

Determination of source finiteness and depth of large earthquakes

In this thesis, a method to retrieve the source finiteness, depth of faulting, and the mechanisms of large earthquakes from long-period surface waves is developed and applied to several recent large events.

In Chapter 1, source finiteness parameters of eleven large earthquakes were determined from long-period Rayleigh waves recorded at IDA and GDSN stations. The basic data set is the seismic spectra of periods from 150 to 300 sec. Two simple models of source finiteness are studied. The first model is a point source with finite duration. In the determination of the duration or source-process times, we used Furumoto's phase method and a linear inversion method, in which we simultaneously inverted the spectra and determined the source-process time that minimizes the error in the inversion. These two methods yielded consistent results. The second model is the finite fault model. Source finiteness of large shallow earthquakes with rupture on a fault plane with a large aspect ratio was modeled with the source-finiteness function introduced by Ben-Menahem. The spectra were inverted to find the extent and direction of the rupture of the earthquake that minimize the error in the inversion. This method is applied to the 1977 Sumbawa, Indonesia, 1979 Colombia-Ecuador, 1983 Akita-Oki, Japan, 1985 Valparaiso, Chile, and 1985 Michoacan, Mexico earthquakes. The method yielded results consistent with the rupture extent inferred from the aftershock area of these earthquakes.

In Chapter 2, the depths and source mechanisms of nine large shallow earthquakes were determined. We inverted the data set of complex source spectra for a moment tensor (linear) or a double couple (nonlinear). By solving a least-squares problem, we obtained the centroid depth or the extent of the distributed source for each earthquake. The depths and source mechanisms of large shallow earthquakes determined from long-period Rayleigh waves depend on the models of source finiteness, wave propagation, and the excitation. We tested various models of the source finiteness, Q, the group velocity, and the excitation in the determination of earthquake depths.

The depth estimates obtained using the Q model of Dziewonski and Steim (1982) and the excitation functions computed for the average ocean model of Regan and Anderson (1984) are considered most reasonable. Dziewonski and Steim's Q model represents a good global average of Q determined over a period range of the Rayleigh waves used in this study. Since most of the earthquakes studied here occurred in subduction zones Regan and Anderson's average ocean model is considered most appropriate.

Our depth estimates are in general consistent with the Harvard CMT solutions. The centroid depths and their 90 % confidence intervals (numbers in the parentheses) determined by the Student's t test are: Colombia-Ecuador earthquake (12 December 1979), d = 11 km, (9, 24) km; Santa Cruz Is. earthquake (17 July 1980), d = 36 km, (18, 46) km; Samoa earthquake (1 September 1981), d = 15 km, (9, 26) km; Playa Azul, Mexico earthquake (25 October 1981), d = 41 km, (28, 49) km; El Salvador earthquake (19 June 1982), d = 49 km, (41, 55) km; New Ireland earthquake (18 March 1983), d = 75 km, (72, 79) km; Chagos Bank earthquake (30 November 1983), d = 31 km, (16, 41) km; Valparaiso, Chile earthquake (3 March 1985), d = 44 km, (15, 54) km; Michoacan, Mexico earthquake (19 September 1985), d = 24 km, (12, 34) km.

In Chapter 3, the vertical extent of faulting of the 1983 Akita-Oki, and 1977 Sumbawa, Indonesia earthquakes are determined from fundamental and overtone Rayleigh waves. Using fundamental Rayleigh waves, the depths are determined from the moment tensor inversion and fault inversion. The observed overtone Rayleigh waves are compared to the synthetic overtone seismograms to estimate the depth of faulting of these earthquakes. The depths obtained from overtone Rayleigh waves are consistent with the depths determined from fundamental Rayleigh waves for the two earthquakes. Appendix B gives the observed seismograms of fundamental and overtone Rayleigh waves for eleven large earthquakes.

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

Characterization of diverse megathrust fault behavior related to seismic supercycles, Mentawai Islands, Sumatra

Long paleoseismic histories are necessary for understanding the full range of behavior of faults, as the most destructive events often have recurrence intervals longer than local recorded history. The Sunda megathrust, the interface along which the Australian plate subducts beneath Southeast Asia, provides an ideal natural laboratory for determining a detailed paleoseismic history over many seismic cycles. The outer-arc islands above the seismogenic portion of the megathrust cyclically rise and subside in response to processes on the underlying megathrust, providing uncommonly good illumination of megathrust behavior. Furthermore, the growth histories of coral microatolls, which record tectonic uplift and subsidence via relative sea level, can be used to investigate the detailed coseismic and interseismic deformation patterns. One particularly interesting area is the Mentawai segment of the megathrust, which has been shown to characteristically fail in a series of ruptures over decades, rather than a single end-to-end rupture. This behavior has been termed a seismic “supercycle.” Prior to the current rupture sequence, which began in 2007, the segment previously ruptured during the 14th century, the late 16th to late 17th century, and most recently during historical earthquakes in 1797 and 1833. In this study, we examine each of these previous supercycles in turn.

First, we expand upon previous analysis of the 1797–1833 rupture sequence with a comprehensive review of previously published coral microatoll data and the addition of a significant amount of new data. We present detailed maps of coseismic uplift during the two great earthquakes and of interseismic deformation during the periods 1755–1833 and 1950–1997 and models of the corresponding slip and coupling on the underlying megathrust. We derive magnitudes of Mw 8.7–9.0 for the two historical earthquakes, and determine that the 1797 earthquake fundamentally changed the state of coupling on the fault for decades afterward. We conclude that while major earthquakes generally do not involve rupture of the entire Mentawai segment, they undoubtedly influence the progression of subsequent ruptures, even beyond their own rupture area. This concept is of vital importance for monitoring and forecasting the progression of the modern rupture sequence.

Turning our attention to the 14th century, we present evidence of a shallow slip event in approximately A.D. 1314, which preceded the “conventional” megathrust rupture sequence. We calculate a suite of slip models, slightly deeper and/or larger than the 2010 Pagai Islands earthquake, that are consistent with the large amount of subsidence recorded at our study site. Sea-level records from older coral microatolls suggest that these events occur at least once every millennium, but likely far less frequently than their great downdip neighbors. The revelation that shallow slip events are important contributors to the seismic cycle of the Mentawai segment further complicates our understanding of this subduction megathrust and our assessment of the region’s exposure to seismic and tsunami hazards.

Finally, we present an outline of the complex intervening rupture sequence that took place in the 16th and 17th centuries, which involved at least five distinct uplift events. We conclude that each of the supercycles had unique features, and all of the types of fault behavior we observe are consistent with highly heterogeneous frictional properties of the megathrust beneath the south-central Mentawai Islands. We conclude that the heterogeneous distribution of asperities produces terminations and overlap zones between fault ruptures, resulting in the seismic “supercycle” phenomenon.

Exploiting seismic waveforms of ambient noise and earthquakes

In this thesis, I apply detailed waveform modeling to study noise correlations in different environments, and earthquake waveforms for source parameters and velocity structure.

Green's functions from ambient noise correlations have primarily been used for travel-time measurement. In Part I of this thesis, by detailed waveform modeling of noise correlation functions, I retrieve both surface waves and crustal body waves from noise, and use them in improving earthquake centroid locations and regional crustal structures. I also present examples in which the noise correlations do not yield Green's functions, yet the results are still interesting and useful after case-by-case analyses, including non-uniform distribution of noise sources, spurious velocity changes, and noise correlations on the Amery Ice Shelf.

In Part II of this thesis, I study teleseismic body waves of earthquakes for source parameters or near-source structure. With the dense modern global network and improved methodologies, I obtain high-resolution earthquake locations, focal mechanisms and rupture processes, which provide critical insights to earthquake faulting processes in shallow and deep parts of subduction zones. Waveform modeling of relatively simple subduction zone events also displays new constraints on the structure of subducted slabs.

In summary, behind my approaches to the relatively independent problems, the philosophy is to bring observational insights from seismic waveforms in critical and simple ways.