This dissertation explores subduction zone-related deformation both on short time scales in the form of subduction zone earthquakes and over larger time and geographical scales in the form of subduction rollback or detachment of the subducting lithosphere. The study presented here is composed of two parts. First, we analyzed the source-rupture processes of the April 1, 2007 Solomon Islands Earthquake (Mw=8.1) using a body-wave inversion technique. Our analysis indicated that the earthquake ruptured approximately 240 km of the southeast Pacific subduction zone in two sub-events.In the second part of this study, we used shear-wave splitting analysis to investigate the effects of the subducting African lithosphere on the upper-mantle flow field beneath the Anatolian Plate in the Eastern Mediterranean region. Our shear-wave splitting results are consistent with relatively uniform southwest-directed flow towards the actively southwestward-retreating Aegean slab. Based on spatial variations in observed delay times we identified varying flow speeds beneath Anatolia and we attribute this variation to the differential retreat rates of the Aegean and the Cyprean trenches.Finally, we used teleseismic P-wave travel-time tomography to image the geometry of the subducting African lithosphere beneath the Anatolia region. Our tomograms show that the subducting African lithosphere is partitioned into at least two segments along the Cyprean and the Aegean trenches. We observed a gap between the two segments through which hot asthenosphere ascends beneath the volcanic fields of western Anatolia. Our results show that the Cyprean slab is steeper than the Aegean slab. We inferred that this steep geometry, in part, controls the flow regime of asthenosphere beneath Anatolia causing variations in flow speeds inferred from shear-wave splitting analysis.

Flow control of transitional shear flows is investigated by means of numerical simulations. The attenuation of three-dimensional wavepackets of Tollmien-Schlichting (TS) and streaks in the boundary layer is obtained using active control in combination with localised sensors and actuators distributed near the rigid wall. Due to the dimensions of the discretized Navier-Stokes operator, reduced-order models are identified, preserving the dynamics between the inputs and the outputs of the system. Balanced realizations of the system are computed using balanced truncation and system identification. We demonstrate that the energy growth of the perturbations is substantially and efficiently mitigated, using relatively few sensors and actuators. The robustness of the controller is analysed by varying the number of actuators and sensors, the Reynolds number, the pressure gradient and by investigating the nonlinear, transitional case. We show that delay of the transition from laminar to turbulent flow can be achieved despite the fully linear approach. This configuration can be reproduced in experiments, due to the localisation of sensing and actuation devices. The closed-loop system has been investigated for the corresponding twodimensional case by using full-dimensional optimal controllers computed by solving an iterative optimisation based on the Lagrangian approach. This strategy allows to compare the results achieved using open-loop model reduction with model-free controllers. Finally, a parametric analysis of the actuators/ sensors placement is carried-out to deepen the understanding of the inherent dynamics of the closed-loop. The distinction among two different classes of controllers – feedforward and feedback controllers - is highlighted. A second shear flow, a confined turbulent jet, is investigated using particle image velocimetry (PIV) measurements. Proper orthogonal decomposition (POD) modes and Koopman modes via dynamic mode decomposition (DMD) are computed and analysed for understanding the main features of the flow. The frequencies related to the dominating mechanisms are identified; the most energetic structures show temporal periodicity. / <p>QC 20130207</p>

Two types of shear flows are investigated in this thesis; numerical simulations are performed for the analysis and control of the perturbation arising in a boundary layer over a flat plate, whereas PIV measurements are analysed for the investigation of a confined turbulent jet. Modal structures of the flows are identified: the aim is to understand the flow phenomena and to identify reduced-order models for the feedback control design. The attenuation of three-dimensional wavepackets of streaks and Tollmien-Schlichting (TS) waves in the boundary layer is obtained using feedback control based on arrays of spatially localized sensors and actuators distributed near the rigid wall. In order to tackle the difficulties arising due to the dimension of the discretized Navier-Stokes operator, a reduced-order model is identified, preserving the dynamics between the inputs and the outputs; to this end, approximate balanced truncation is used. Thus, control theory tools can be easily handled using the low-order model. We demonstrate that the energy growth of both TS wavepackets and streak-packets is substantially and efficiently mitigated, using relatively few sensors and actuators. The robustness of the controller is investigated by varying the number of actuators and ensors, the Reynolds number and the pressure gradient. The configuration can be possibly reproduced in experiments, due to the localization of sensing and actuation devices. A complete analysis of a confined turbulent jet is carried out using timeresolved PIV measurements. Proper orthogonal decomposition (POD) modes and Koopman modes are computed and analysed for understanding the main features of the flow. The frequencies related to the dominating mechanisms are identified; the most energetic structures show temporal periodicity. / QC 20110214

The Alpine Fault is a major plate boundary structure, which accommodates up to 50-80% of the total plate boundary motion across the South Island of New Zealand. The fault has not ruptured historically although limited off-fault shaking records and on-fault dating suggest large to great (~ Mw 8) earthquakes (every ~100-480 years; most recently in 1717), making it potentially one of the largest onshore sources of seismic hazard in New Zealand. The central section of the Alpine Fault, which bounds the highest elevations in the Southern Alps, is one of the most poorly characterised sections along the fault. On-fault earthquake timing in addition to the amount of dextral slip during major earthquakes was unknown along a 200-km-long section of the central Alpine Fault, while the amount of co-seismic hanging wall uplift was poorly known, prior to the present work. In this thesis I address these knowledge gaps through a combination of light detection and ranging (lidar), field, and stratigraphic mapping along with sample dating to constrain earthquake timing, style of faulting, and hanging wall rock uplift rates. Using lidar data coupled with field mapping I delineated the main trace of the Alpine Fault at Gaunt Creek as a north-striking fault scarp that was excavated and logged; this is part of a 2-km-wide restraining bend dominated by low-angle thrust faulting and without the clear strike-slip displacements that are present nearby (<5 km distant along strike in both directions). Where exposed in this scarp, the fault-zone is characterized by a distinct 5-50 cm thick clay fault-gouge layer juxtaposing hanging wall bedrock (mylonites and cataclasites) over unconsolidated late-Holocene footwall colluvium. An unfaulted peat at the base of the scarp is buried by post-most recent event (MRE) alluvium and yields a radiocarbon age of A.D. 1710–1930, consistent with sparse on-fault data, validating earlier off-fault records that suggest a 1717 MRE with a moment magnitude of Mw 8.1 ± 0.1, based on the 380-km-long surface rupture. Lidar and field mapping also enabled the identification and measurement of short (<30 m), previously unrecognized dextral offsets along the central section of the Alpine Fault. Single-event displacements of 7.5 ± 1 m for the 1717 earthquake and cumulative displacements of 12.9 ± 2 m and 22 ± 2.7 m for earlier ruptures can be binned into 7.1 ± 2.1 m increments of repeated dextral (uniform) slip along the central Alpine Fault. A comparison of these offsets with the local paleoseismic record and known plate kinematics suggests that the central Alpine Fault earthquakes in the past 1.1 ka may have: (i) bimodal character, with major surface ruptures (!Mw 7.9) every 270 ± 70 years (e.g. the 1717 event) and with moderate to large earthquakes (!Mw 7) occurring between these ruptures (e.g. the 1600 event); or (ii) that some shaking data may record earthquakes on other faults. If (i) is true, the uniform slip model (USM) perhaps best represents central Alpine Fault earthquake recurrence, and argues against the applicability of the characteristic earthquake model (CEM) there. Alternatively, if (ii) is true, perhaps the fault is “characteristic” and some shaking records proximal to plate boundary faults do not necessarily reflect plate-boundary surface ruptures. Paleoseismic and slip data suggest that (i) is the most plausible interpretation, which has implications for the understanding of major plate-boundary faults worldwide. Field mapping, geological characterisation, geophysical mapping, and optically stimulated luminescence (OSL) dating of on-fault hanging wall sediments were used to better constrain the geometry and kinematics of Holocene deformation along the rangefront of the Southern Alps at the Alpine Fault near the Whataroa River. The fault here is dextral-reverse, although primarily strike-slip with clear fault traces cutting across older surfaces of varying elevations. Deformational bulges are observed along these traces that are likely thrust-bounded. A terrace of Whataroa River sediments was found on the hanging wall of the Alpine Fault approximately ~ 55-75 m (when considering uncertainties) above the floodplain of the Whataroa River. OSL ages for a hanging wall sediments of 10.9 ± 1.0 ka for the aforementioned terrace, 2.8 ± 0.3 ka for Whataroa River terrace deposits in a deformational bulge, and 11.1 ± 1.2 ka for a rangefront derived fan indicate Holocene aggradation along the rangefront and hanging wall uplift rates of 6.0 ± 1.1 mm/yr. The sub-horizontal, laterally continuous, and planar-bedded Whataroa-sourced terrace deposits suggest that the adjacent bounding faults are steeply-dipping faults without geometries in the shallow subsurface that would tend to cause sedimentary bed rotation and tilting. Using data from the approximately 100-m deep pilot DFDP boreholes together with lidar and field mapping, I present a review of the Quaternary geology, geomorphology, and structure of the fault at Gaunt Creek, and estimate new minimum Late-Pleistocene hanging wall rock uplift rates of 5.7 ± 1.0 mm/yr to 6.3 ± 1.1 mm/yr (without considering local erosion) that suggest that the Southern Alps are in a dynamic steady state here. GPS-derived “interseismic” vertical uplift rates are < 1 mm/yr at the Alpine Fault, so the majority of rock uplift at the rangefront happens during episodic major earthquakes, confirming with on-fault data that slip occurs coseismically. Notably the uplift rates from both Mint and Gaunt Creek are consistent between the two sites although the primary style of faulting at the surface is different between the two sites, suggesting consistent coseisimc uplift of the Southern Alps rangefront along the Alpine Fault in major earthquakes. This thesis collected new on-fault datasets that confirm earlier inferences of plate-boundary fault behaviour. This study of the high-uplift central section of the Alpine Fault provides the first on-fault evidence for the MRE (i.e. 1717) and repeated of dextral slip during the MRE and previous events as well as new hanging wall uplift data which suggests that the majority of rangefront uplift occurs in earthquakes along the Alpine Fault. Because the fault has not ruptured for ~300 years, it poses a significant seismic hazard to southern New Zealand.

Campaign GPS measurements during 2001-2003 in the transition between subduction and collision of the Banda arc reveal how strain is partitioned away from the trench and distributed to other parts of the arc-trench system. Genrich, et. al. (1996) conducted a GPS campaign (1992-1994) throughout the Eastern Sunda and Banda arcs that demonstrated partial accretion of the arc to the Australian plate. We reoccupied many of the sites from this earlier study and 7 additional stations, 3 of which are new benchmarks. Our study shortened many baselines and extended the observation epoch to ten years for many key stations. The resulting GPS velocity field for the active Banda arc-continent collision reveals: 1) several mostly rigid crustal blocks exist in the transition from subduction to collision, 2) relative to an Asian reference frame, most of these blocks move in the same direction as the Australian lower plate, but at different rates, 3) block boundaries may exist between the islands of Lombok and Komodo, Flores and Sumba, Savu and West Timor, and between Timor and Darwin, 4) the Timor Trough may account for at least 20 mm/yr of motion between Timor and Darwin, 5) a major transverse fault off the coast of West Timor separates the Savu/Flores/Sumba block from the Timor/Wetar Block and may account for variations in movement in Rote, 6) the Flores thrust moves the eastern Sunda arc north relative to Asia by decreasing amounts to the west, 7) the back-arc Wetar Thrust system takes up the majority of plate convergence between Australia and Asia, and 8) fault boundaries are not found between many blocks, such as various islands of the Sunda arc and forearc with different amounts of motion.

Temporally developing direct numerical simulations (T-DNS) are performed and validated for bypass transition of a zero pressure gradient flat plate boundary layer to understand the interplay between pressure-strain terms and flow instability mechanisms, and to propose and validate a phenomenological hypothesis for the identification of a robust transition onset marker for use in transition-sensitive Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations. Results show that transition initiates at a location where the slow pressure-strain term becomes more dominant than the rapid term in the pre-transitional boundary layer region. The slow pressure strain term is responsible for the transfer of turbulence energy from the streamwise component to other components while the rapid pressure strain term counteracts with the slow term in the pre-transitional regime before transition onset akin to a shear sheltering like effect. The relative magnitudes of the slow and rapid terms thus provide a basis for the development of physically meaningful large-scale parameters that can be used as a transition onset marker for Reynolds averaged Navier-Stokes RANS simulations.

L’archipel des Açores dans l’océan Atlantique est édifiées sur un épais plateau océanique, à proximité de la jonction triple entre les plaques Nord-américaine (Na), Nubienne (Nu) et Eurasienne (Eu). La formation du plateau et l’origine du volcanisme ont été le plus souvent attribués à la présence d’une instabilité mantellique. Cependant, la répartition et la morphologie des édifices volcaniques semblent avoir été grandement influencés par la déformation régionale liée à la migration de la frontière de plaque (Eu/Nu). En effet, la frontière serait passée d’une faille transformante aujourd’hui inactive, la zone de fracture est des Açores (EAFZ), à un rift ultra lent actif appelé le Rift de Terceira (TR).Lors de ce travail, nous utilisons le volcanisme comme marqueur de la déformation régionale. Nous nous intéressons particulièrement aux îles de S. Miguel et Graciosa, qui sont localisées à l’intérieur du TR, et à Santa Maria, une île volcanique éteinte qui se situe entre la EAFZ et le TR. De par leur position, ces trois îles constituent donc des cibles particulièrement appropriées afin d’étudier l’architecture et l’évolution de la frontière de plaque Eu/Nu durant les dernier Millions d’années. A partir de nouvelles données géomorphologiques, stratigraphiques, géochronologiques et tectoniques, couplées aux données bathymétriques et géophysiques disponibles, nous reconstruisons les étapes successives de construction et de démantèlement de ces îles puis discutons de leur signification géodynamique. Ces données sont ensuite complétées par des expériences de mécanique des fluides afin d’investiguer les liens possibles entre un panache mantellique, la migration de la frontière de plaque sur plusieurs échelles d’espace et de temps.Les résultats montrent que les édifices localisés dans le TR se construisent via des pulses volcaniques courts (<100 kyr) et relativement synchrones, séparés par des épisodes d’effondrements catastrophiques. Nous proposons qu’une telle évolution reflète des épisodes brefs et intenses de déformation régionale le long de la frontière de plaque active. La distribution des marqueurs tectoniques ainsi que leurs orientations N110 et N150 dans la partie Est de S. Miguel, nous conduit à proposer que l’extension oblique du TR est principalement accommodée par les failles bordières majeures du rift. Nous identifions une nouvelle tendance tectonique orientée N50° qui pourrait représenter des failles transformantes accommodant les variations d’obliquité du TR. L’activité de île de Santa Maria est ici datée entre 5.7 et 2.8 Ma. S. Maria a été façonnée par plusieurs effondrements sectoriels catastrophiques, le plus probablement déclenchés par les mouvements tectoniques régionaux. Nous identifions également une nouvelle structure de type graben reliant les îles de S. Maria et S. Jorge plus loin au NW. La forme de ce graben est semblable au TR et est située entre l’ancienne et la nouvelle frontière Eu/Nu. Nous interprétons ce graben comme un ancien rift transitionnel et donc comme une ancienne frontière de plaque Eu/Nu. A partir de nos données géochronologiques, nous proposons que la partie Est de ce rift transitionnel aurait migré vers la partie Est du TR entre 2.8 et 1.7 Ma.La migration de la frontière Eu/Nu a été interprétées par Vogt and Jung (2004) comme résultant de sauts successifs vers le NE de l’axe du Rift afin de maintenir sa position au dessus d’un point chaud fixe. Nos expériences de mécanique des fluides suggèrent que l’archipel des Açores, comme celui des Canaries, du Cap Vert, de Madère ainsi que les volcans sous marins de Great Meteor sont la signature en surface d’un groupe d’instabilités mantellique prenant naissance et remontant à partir du sommet d’un dôme thermochimique situé dans le manteau inférieur. De plus, Ces panaches secondaires pourraient être suffisamment faibles pour adapter leurs mouvements aux équilibres de forces pré-existants, notamment la structure et la morphologie de la lithosphère. / The Azores archipelago in the Atlantic comprises nine volcanic islands which developed on a thick oceanic plateau close to the Triple Junction between the North American (Na), the Nubian (Nu), and the Eurasian (Eu) lithospheric plates. The formation of the plateau and the origin of the volcanism remain controversial, but have been generally attributed to a plume-like mantle instability. However, the distribution of the volcanic edifices east of the Mid-Atlantic Ridge (MAR) appears greatly influenced by regional deformation associated with the northward migration of the Eu/Nu plate boundary from an extinct old transform fault, the East Azores Fracture Zone (EAFZ), up to the presently active ultra-slow Terceira Rift (TR). In this thesis, we use the volcanism as a marker for regional deformation. We especially focus on S. Miguel and Graciosa, which are located within the TR, and on S. Maria, an old volcanically extinct island located between the EAFZ and the TR. These three islands thus constitute particularly suitable targets to track the architecture and the evolution of the Eu/Nu plate boundary during the last few Myr. From new geomorphological, stratigraphic, geochronologic, structural/tectonic data, and existing bathymetric and geophysical data, we reconstruct the successive stages of growth and destruction of the islands, and discuss their geodynamic meaning. These data are then complemented by fluid dynamic modelling using laboratory experiments to examine the possible links between mantle instability, plate boundary migration and the development of the volcanism on various spatial and temporal scales.The new results on the islands show that the edifices located within the TR grew through short (<100 kyr) and partly synchronous volcanic pulses, separated by catastrophic sector collapses. We propose that such evolution reflects brief and intense episodes of regional deformation along the still active Eu/Nu plate boundary. The distribution of tectonic markers and the recognition of N110 and N150 tectonic structures in eastern S. Miguel leads us to propose that oblique extension in the TR is mainly accommodated by the master faults of the rift, and that the TR is presently not the locus of appreciable sea-floor spreading. Furthermore, we identify a new N050 trend, which may represent transform faults accommodating the variation in obliquity of the TR. The activity of S. Maria is here dated between 5.7 and 2.8 Ma. Like the recent islands, S. Maria experienced catastrophic flank collapses, most probably triggered by regional tectonics. We identify a new graben structure linking Santa Maria to the island of S. Jorge further NW. The shape of this graben is similar to the TR and it is located between the EAFZ and the current plate boundary. We interpret this graben as a former transient rift, and therefore an old Eu/Nu plate boundary. From the new data, we propose that the eastern part of the transient rift migrated to the eastern part of the TR between 2.8 Ma and 1.7 Ma.The overall migration of the Eu/Nu plate boundary to the north and the creation of the Azores plateau has been interpreted by Vogt and Jung (2004) as resulting from successive NE jumps of the rift axis to maintain its position over a fixed ‘hotspot’. Our fluid mechanics experiments suggest that the Azores, as Canary, Cape Verde, Madeira Islands and Great Meteor seamounts might be the surface signature of a cluster of mantle instabilities rising from the top of a large thermochemical dome located in the lower mantle. However, such secondary plumes present a strong time-dependence 5-40 Myr time scale. Moreover, they could be sufficiently weak to adapt their motions to the pre-existing force balances and morphology of the lithosphere. We therefore present a scenario of the Azores area evolution combining a triple junction and decompression melting buoyant material (i.e. such in volatiles and/or temperature) under a thickening lithosphere.

The upper 8-12 km of the Alpine Fault, South Island, New Zealand, accommodates relative Australia-Pacific plate boundary motion through coseismic slip accompanying large-magnitude earthquakes. Earthquakes occur due to frictional instabilities on faults, and their nucleation, propagation, and arrest is governed by tectonic forces and fault zone properties. A multi-disciplinary dataset is presented on the lithological, microstructural, mineralogical, geochemical, hydrological, and frictional properties of Alpine Fault rocks collected from natural fault exposures and from Deep Fault Drilling Project (DFDP-1) drillcore. Results quantify and describe the physical and chemical processes that affect seismicity and slip accommodation. Oblique dextral motion on the central Alpine Fault in the last 5-8 Myr has exhumed garnet-oligoclase facies mylonitic fault rocks from depths of up to 35 km. During the last phase of exhumation, brittle deformation of these mylonites, accompanied by fluid infiltration, has resulted in complex mineralogical and lithological variations in the fault rocks. Petrophysical, geochemical, and lithological data reveal that the fault comprises a central alteration zone of protocataclasites, foliated and nonfoliated cataclasites, and fault gouges bounded by a damage zone containing fractured ultramylonites and mylonites. Mineralogical results suggest that at least two stages of chemical alteration have occurred. At, or near, the brittle-to-ductile transition (c. >320 °C), metasomatic alteration reactions resulted in plagioclase and feldspar replacement by muscovite and sausserite, and biotite (phlogopite), hornblende (actinolite) and/or epidote replacement by chlorite (clinochlore). At lower temperatures (c. >120°C), primary minerals were altered to kaolinite, smectite and pyrite, or kaolinite, smectite, Fe-hydroxide (goethite) and carbonate, depending on redox conditions. Ultramylonites, nonfoliated and foliated cataclasites, and gouges in the hanging wall and footwall contain the high-temperature phyllosilicates chlorite and white mica (muscovite/illite). Brown principal slip zone (PSZ) gouges contain the low-temperature phyllosilicates kaolinite and smecite, and goethite and carbonate cements. The frictional and hydrological properties of saturated intact samples of central Alpine Fault surface-outcrop gouges and cataclasites were investigated in room temperature experiments conducted at 30-33 MPa effective normal stress (σn') using a double-direct shear configuration and controlled pore fluid pressure in a triaxial pressure vessel. Surface-outcrop samples from Gaunt Creek, location of DFDP-1, displayed, with increasing distance (up to 50 cm) from the contact with footwall fluvioglacial gravels: (1) an increase in fault normal permeability (k = 7.45 x 10⁻²⁰ m² to k = 1.15 x 10⁻¹⁶ m²), (2) a transition from frictionally weak (μ=0.44) fault gouge to frictionally strong (μ=0.50’0.55) cataclasite, (3) a change in friction rate dependence (a–b) from solely velocity strengthening to velocity strengthening and weakening, and (4) an increase in the rate of frictional healing. The frictional and hydrological properties of saturated intact samples of southern Alpine Fault surface-outcrop gouges were also investigated in room temperature double-direct shear experiments conducted at σn'= 6-31 MPa. Three complete cross-sections logged from outcrops of the southern Alpine Fault at Martyr River, McKenzie Creek, and Hokuri Creek show that dextral-normal slip is localized to a single 1-12 m-thick fault core comprising impermeable (k=10⁻²⁰ to 10⁻²² m²), frictionally weak (μ=0.12 – 0.37), velocity-strengthening, illite-chlorite and trioctahedral smectite (saponite)-chlorite-lizardite fault gouges. In low velocity room temperature experiments, Alpine Fault gouges tested have behaviours associated with aseismic creep. In a triaxial compression apparatus, the frictional properties of PSZ gouge samples recovered from DFDP-1 drillcore at 90 and 128 m depths were tested at temperatures up to T=350°C and effective normal stresses up to σn'=156 MPa to constrain the fault's strength and stability under conditions representative of the seismogenic crust. The chlorite/white mica-bearing DFDP-1A blue gouge is frictionally strong (μ=0.61–0.76) across a range of experimental conditions (T=70–350°C, σn'=31.2–156 MPa) and undergoes a stability transition from velocity strengthening to velocity weakening as T increases past 210°C, σn'=31.2–156 MPa. The coefficient of friction of smecite-bearing DFDP-1B brown gouge increases from μ=0.49 to μ=0.74 with increasing temperature and pressure (T=70–210°C, σn'=31.2–93.6 MPa) and it undergoes a transition from velocity strengthening to velocity weakening as T increases past 140°C, σn'=62.4 MPa. In low velocity hydrothermal experiments, Alpine Fault gouges have behaviours associated with potentially unstable, seismic slip at temperatures ≥140°C, depending on mineralogy. High-velocity (v=1 m/s), low normal stress (σn=1 MPa) friction experiments conducted on a rotary shear apparatus showed that the peak coefficient of friction (μp) of Alpine Fault cataclasites and fault gouges was consistently high (mean μp=0.69±0.06) in room-dry experiments. Variations in fault rock mineralogy and permeability were more apparent in experiments conducted with pore fluid, wherein the peak coefficient of friction of the cataclasites (mean μp=0.64±0.04) was higher than the fault gouges (mean μp=0.24±0.16). All fault rocks exhibited very low steady state coefficients of friction (μss) (room-dry mean μss=0.18±0.04; saturated mean μss=0.10±0.04). Three high-velocity experiments conducted on saturated smectite-bearing principal slip zone (PSZ) fault gouges had the lowest peak friction coefficients (μp=0.13-0.18), lowest steady state friction coefficients (μss=0.02-0.10), and lowest breakdown work values (WB=0.07-0.11 MJ/m²) of all the experiments performed. Lower strength (μ < c. 0.62) velocity-strengthening fault rocks comprising a realistically heterogeneous fault plane represent barrier(s) to rupture propagation. A wide range of gouges and cataclasites exhibited very low steady state friction coefficients in high-velocity friction experiments. However, earthquake rupture nucleation in frictionally strong (μ ≥ c. 0.62), velocity-weakening material provides the acceleration necessary to overcome the low-velocity rupture propagation barrier(s) posed by velocity-strengthening gouges and cataclasites. Mohr-Coulomb theory stipulates that sufficient shear stress must be resolved on the Alpine Fault, or pore fluid pressure must be sufficiently high, for earthquakes to nucleate in strong, unstable fault materials. A three-dimensional stress analysis was conducted using the average orientation of the central and southern Alpine Fault, the experimentally determined coefficient of friction of velocity-weakening DFDP-1A blue gouge, and the seismologically determined stress tensor and stress shape ratio(s). Results reveal that for a coefficient of friction of μ ≥ c. 0.62, the Alpine Fault is unfavourably oriented to severely misoriented for frictional slip.

Tsunami hazards were greatly underestimated along the coasts of countries bordering the northeastern Indian Ocean until the occurrence of the 26 December 2004, Mw 9.2 earthquake and its ensuing tsunami. Sourced off the coast of northern Sumatra, on the plate boundary between the Indo-Australian and Eurasian plates, the rupture of the 2004 earthquake propagated ~1300 km northward. The magnitude of this earthquake and the reach of its tsunami exceeded all known precedents, based on instrumental and historic records. The coseismic deformational and post-tsunami depositional features facilitated opportunities to conduct tsunami geology studies along the coasts of countries bordering the Indian Ocean. Several questions are being posed, the answers of which have implications for tsunami hazard assessment. How did this plate boundary behave prior to and after the great earthquake? Was the 2004 earthquake the first of its kind on the Sumatra-Andaman plate boundary? If it had a predecessor, when did it occur and was it a true predecessor in terms of its rupture dimensions and tsunamigenic potential? What types of depositional evidence are preserved and how can we use them to develop the history of past tsunamigenic earthquakes? Researchers are exploring the affected regions and using the imprints left by the 2004 event, to address these questions. There are two components to this study: one, a seismotectonic analysis of the region from the perspective of plate driving forces and their relative roles in the interseismic and post-seismic phases. This study uses global data catalogs like the NEIC PDE (National Earthquake Information Centre Preliminary Determination of Epicenters) and the Global Centroid Moment Tensor (CMT) solutions for earthquake source parameters to understand the along-strike variations in seismicity patterns before and after the 2004 earthquake. The 2004 experience was unprecedented in South Asia. Unaffected by tsunami hazards in the past, tsunami geology is a nascent field for most South Asian researchers. Very little background field data is available on the deformational features of great earthquakes along this plate boundary and the depositional characteristics of extreme coastal surges, such as tsunamis and storms. Where do we begin our search for evidence of past tsunamigenic earthquakes? How best can we use the 2004 tsunami and its deposits as a proxy? What problems are encountered in the interpretations? This thesis addresses these questions in part and presents observations from the Andaman Islands (the ~400 km, northern segment of the Sumatra-Andaman subduction zone) and the southeast coast of India, towards developing a reliable database of tsunami geology for 2004-type events. The premise is that regions affected by the 2004 earthquake are more likely to conserve signatures from older events. Based on the stratigraphic context of the proxy and quality of age estimates, this work presents evidence for past earthquake related deformation and tsunami deposition. In this work we use deformational and depositional features from the Andaman Islands, falling within the 2004 rupture zone and from one location on the Tamil Nadu coast of India (Kaveripattinam). From a perceptive understanding of the features related to tectonic deformation of the Sumatra-Andaman subduction zone, we have selected the Andaman segment that demonstrates explicit evidence for deformation and tsunami deposition through geomorphological and stratigraphic features, which are key to our exploration. A gist of each chapter is given below. The introduction (chapter 1) presents the background, motivation and scope of this work and the organization of this thesis, also summarizing the contents of each chapter. Chapter 2 provides a review of literature on subduction zone earthquakes and updates on tsunami geology, to place this study in the global context. The next two chapters discuss the seismotectonics of the Sumatra-Andaman plate boundary, the important earthquakes and their source processes. In chapter 3 we discuss the Andaman segment (from 10–15° N), characterized by relatively lower level seismicity, but distinctive, as it falls within the northern limit of the 2004 rupture. The deformational and depositional features here are better exposed due to availability of land straddling the hinge line separating the areas of 2004 uplift and subsidence. Here, the pre-2004 earthquakes used to occur along a gently dipping subducting slab, up to a depth of about 40 km. Post-2004, the earthquakes moved up-dip, extending also to the outer-rise and outer-ridge regions, expressing post-earthquake relaxation [Andrade and Rajendran, 2011]. The southern Nicobar segment (5–10° N) differs from the Andaman segment in its style of deformation and seismic productivity. The decreasing obliquity of convergence, the likely influence of a subducting ocean ridge on the subducting plate and the character of the subducting oceanic plate make this segment distinctly different. In chapter 4 we present an analysis of its seismotectonic environment based on the well-constrained focal mechanisms of historic and recent earthquakes. We report that left-lateral strike-slip faulting on near N-S oriented faults control the deformation and the style of faulting is consistent to ~80 km within the subducting slab [Rajendran, K. et al., 2011]. The 11 April 2012 sequence of earthquakes on the subducting oceanic plate, between the Sumatra Trench and the Ninety East Ridge are the more recent among the oceanic intraplate earthquakes that demonstrate the reactivation of N-S oriented fossil fractures. The limited availability of land and the 2004 coseismic deformation dominated by subsidence, followed by prolonged waterlogging makes exploration difficult in the Nicobar segment. Thus, we focus on the Andaman Islands for deformational and depositional evidence, using observations that can be corroborated through multiple proxies and depositional environments that are not prone to other coastal surges, such as cyclones and storms. The criteria for selection of sites, evaluation of deposits and determination of limiting ages are discussed in chapters 5 through 9. In chapter 5 we discuss different types of coastal environments and their response to high-energy sea surges. We also give a brief review of the comparative analyses of storm and tsunami deposits, a highly debated issue and then discuss important characteristics of these two deposits, using examples from the 2004 tsunami and the 2011 Thane cyclone that affected parts of the Tamil Nadu coast. An important component of tsunami geology is the ability to identify and select datable material from tsunami deposits and chose an appropriate method for dating (chapter 6). The types of material used vary from peat layers, peat-rich soil, gastropod shells, wood, charcoal, organic remains such as bones, coral fragments, pottery sherds and buried soil. Techniques such as AMS Carbon-14 and Thermoluminescence are commonly used with appropriate calibrations and corrections. In addition to the dates generated in this study (based on wood and shell dates) we use some previous dates from the entire stretch of the rupture within the Indian Territory and assign a relative grading to these ages, based on the quality criterion evolved in this study. We believe that this is the first attempt to segregate age data obtained from coastal deposits, and assign them a specific quality grading based on their environment of deposition and the type of material dated. Chapter 7 presents results of our investigations in the Andaman Islands, which cover ~30% of the rupture area. A coseismically subsided mangrove from Rangachanga (Port Blair, east coast of South Andaman) led us to a former subsidence during AD 770–1040, which we believe is the most convincing evidence for a previous tectonic event. Data based on inland deposits of coral and organic debris yielded a younger age in the range of AD 1480–1660. Both these dates fall in the age brackets reported from other regions of this plate boundary (mainly Sumatra) as well as distant shores of Sri Lanka, Thailand and mainland India. To understand the nature of distant deposits, we present observations from Kaveripattinam, an ancient port city on the east coast of India, where a high-energy sea surge deposit, found 1 km inland is attributed to a paleotsunami. The inland location of this archeological site at an elevation of 2 m and characteristics of the deposit that help discriminate it from typical storm deposition provide clinching evidence in favor of a 1000-year old regional tsunami (chapter 8). In chapter 9 we discuss the results of our study. We evaluate the nature of deformation/deposition and the calibrated age data in the context of their environments. Ages based on the organic material associated with coral debris (at Hut Bay and Interview Island) and the remains of mangrove roots, 1 m below the present ground level (at Port Blair) are considered as reliable estimates, due to their sheltered inland location and the in situ root horizon used for dating. Age data from Kaveripattinam is also considered reliable, based on its inland location beyond the reach of storm surges, sediment characteristics typical of tsunami deposition and ages based on multiple methods and samples. The age data based on the sites presented in this thesis are more conclusive about the 800 to 1100 AD and 1250 to 1450 AD tsunamis, and the former is represented from regions closer to the 2004 source as well as distant shores reached by its tsunami. Chapter 10 presents our conclusions and the scope for future studies. We present this as the first study of its kind in the northeastern Bay of Bengal, wherein the coseismic vertical coastal deformation features along an interplate subduction boundary and a variety of tsunami deposits are used to categorize depositional environments and ages of paleoearthquakes and tsunamis. To our knowledge, this is the first study of its kind where the effects of a recent tsunami have been used to evaluate paleodeposits based on their respective environments of occurrence. Our results have implications for tsunami geology studies in coastal regions prone to tsunami hazard.

Les rides de Barracuda et de Tiburon sont deux reliefs sous-marins situés dans la partie ouest de l'océan Atlantique, là où la lithosphère océanique des plaques Amérique du Nord (NAM) et Amérique du Sud (SAM) est entraînée par subduction sous la plaque Caraïbe, formant l'arc volcanique des Petites Antilles et le prisme d’accrétion de Barbade. Le processus et la période de soulèvement conduisant au relief actuel de ces rides (qui semblent être un marqueur important dans l'histoire géodynamique de la région) sont sujets à débat depuis des décennies.L’interprétation de nouvelles données de sismique réflexion et de bathymétrie multifaisceaux acquises à travers les rides de Barracuda et de Tiburon (campagne Antiplac, 2007 ) a permis de dater les périodes de soulèvements des rides et réaliser des reconstructions paléogéographiques incluant les flux sédimentaires majeurs, depuis le Crétacé jusqu’ à l’Actuel.L’analyse structurale révèle des phases de réactivations tardives d’anciennes zones de fractures dans un contexte transpressif, conduisant aux surrections des rides de Tiburon et de Barracuda.Les processus géologiques possibles impliqués dans la formation des rides de Barracuda et de Tiburon coïncident avec les modèles cinématiques récents décrivant les mouvements relatifs entre les plaques NAM et SAM, le long de la limite de plaque diffuse.Ces résultats permettent de mieux définir la limite de plaque entre NAM et SAM. Elle est nécessairement hétérogène exploitant les zones de faiblesses dans la lithosphère que sont les zones de fracture. Au sein de cette limite de plaque la lithosphère serait donc fragmentée. / The Barracuda Ridge and the Tiburon Rise, two oceanic-basement ridges, lie in the western Atlantic Ocean, where oceanic lithosphere of the North American (NAM) and South American (SAM) plates is subducted beneath the Caribbean plate, creating the Lesser Antilles volcanic arc and the Barbados Ridge accretionary complex. The process and the timing of the uplift leading to the present day morphologies of the Tiburon and Barracuda ridges, that seem to be key markers in the geodynamic history of the region, has remained a matter of debate for decades.From the analysis of new multibeam and seismic reflection profiles acquired in 2007 (Antiplac crusie) DSDP-ODP boreholes available, we provide new information on the timing of the formation of the Barracuda Ridge and Tiburon Rise in their present-day configurations. We propose paleogeographic reconstructions with the main sediments fluxes deposited in the area of the Barracuda and Tiburon ridges from the Late Cretaceous to present. Structural analysis shows reactivation of fracture zones in a transpressive setting leading to the uplifts of the Barracuda and Tiburon Ridges.The location of the Barracuda Ridge and the Tiburon Rise and the timing of the uplift fit well with recent global plate kinematic models describing the movements of NAM relative to SAM along a diffuse plate boundaryThis NAM-SAM plate boundary zone, therefore must most certainly be heterogeneous in nature, exploiting weaknesses in the lithosphere provided by fracture zones where mechanically advantageous, but forming new boundary segments elsewhere, to transfer motion between reactivated segments of the fracture zones.