Constraining the active tectonic deformation of Andaman-Nicobar Arc in the background of December 26,2004 Great Sumatra-Andaman Earthquake

The Andaman-Nicobar Islands in the Bay of Bengal lies in a zone where the Indian plate subducts beneath the Burmese microplate, and therefore forms a belt of frequent earthquakes. Few efforts, not withstanding the available historical and instrumental data were not effectively used before the Mw 9.3 Sumatra-Andaman earthquake to draw any inference on the spatial and temporal distribution of large subduction zone earthquakes in this region. An attempt to constrain the active crustal deformation of the Andaman-Nicobar arc in the background of the December 26, 2004 Great Sumatra-Andaman megathrust earthquake is made here, thereby presenting a unique data set representing the pre-seismic convergence and co-seismic displacement.Understanding the mechanisms of the subduction zone earthquakes is both challenging sCientifically and important for assessing the related earthquake hazards. In many subduction zones, thrust earthquakes may have characteristic patterns in space and time. However, the mechanism of mega events still remains largely unresolved.Large subduction zone earthquakes are usually associated with high amplitude co-seismic deformation above the plate boundary megathrust and the elastic relaxation of the fore-arc. These are expressed as vertical changes in land level with the up-dip part of the rupture surface uplifted and the areas above the down-dip edge subsided. One of the most characteristic pattern associated with the inter-seismic era is that the deformation is in an opposite sense that of co-seismic period.This work was started in 2002 to understand the tectonic deformation along the Andaman-Nicobar arc using seismological, geological and geodetic data. The occurrence of the 2004 megathrust earthquake gave a new dimension to this study, by providing an opportunity to examine the co-seismic deformation associated with the greatest earthquake to have occurred since the advent of Global Positioning System (GPS) and broadband seismometry. The major objectives of this study are to assess the pre-seismic stress regimes, to determine the pre-seismic convergence rate, to analyze and interpret the pattern of co-seismic displacement and slip on various segments and to look out for any possible recurrence interval for megathrust event occurrence for Andaman-Nicobar subduction zone. This thesis is arranged in six chapters with further subdivisions dealing all the above aspects.

Crustal growth in the 3.4-2.7Ga São José de Campestre Massif, Borborema Province, NE Brazil

The Archean (3.45-2.70Ga) rocks of the São José do Campestre Massif (SJCM) in the Borborema Province (NE Brazil) make up a small area (~6000km2) and are composed of granitoids and metasupracrustal rocks that define a complex magmatic and deformational history. The massif provides the opportunity to study mantle- and crustal-derived magmas generated since the Palaeoarchean. The orthogneisses of the SJCM are composed of: (1) tonalite to granodiorite with diorite enclaves (Bom Jesus gneiss, 3412±8Ma; TDM Nd model ages from 4.1 to 3.5Ga and negative epsilon Nd values); (2) biotite and ferroan-diopside monzogranite (Presidente Juscelino complex, 3356±21Ma and 3251±44Ma; TDM model ages range from 4.1 to 3.4Ga and epsilon Nd values that are slightly positive to negative); (3) hornblende tonalite to granodiorite (Brejinho complex, 3333±77Ma and 3187±8Ma; dominantly positive epsilon Nd values and TDM ages from 3.6 to 3.2Ga); (4) biotite monzogranite (São Pedro do Potengi gneiss, 3120±22Ma; TDM =3.5Ga; negative epsilon Nd value); (5) ferroan-diopside-grossular anorthosite and metagabbro (Senador Elói de Souza complex, 3033±3Ma); and (6) quartz diorite to syenogranite (São José do Campestre complex; 2685±9Ma and 2655±4Ma; negative epsilon Nd values and TDM ages from 3.9 to 3.3Ga). The orthogneisses are subalkaline to faintly alkaline, magnesian to ferroan, M- and I-type granitoids that follow either the K-enrichment or the trondhjemite trends. Each group has a subset with REE characteristics similar to Archean TTG and another that is analogous to Phanerozoic granitoids. They have negative Ta-Nb and Ti anomalies and have trace element contents of granitoids from subduction zones. Geochemical and Nd isotope data suggest that subducted oceanic crust and a depleted and metasomatised mantle wedge both acted as the magma sources. We propose a convergent tectonic model in which hybridisation of the upper mantle occurs through interactions with adakitic or trondhjemitic melts and recycling of earlier crust. The results imply that both the subducted oceanic crust and the mantle wedge played major roles in continent formation throughout successive episodes of arc accretion in Palaeo- and Mesoarchean times. The Archean rocks of the SJCM shares some similarities with the Pilbara, Kaapvaal, West African, and São Francisco cratons. However, the most reliable comparisons with the SJCM are with the neighbouring basement of the Nigeria and Cameroon shields. © 2012 Elsevier B.V.

Geochemical constraints on blueschist- and amphibolite-facies rocks of the Central Cordillera of Colombia: the Andean Barragan region

Subduction zones are one of the most characteristic features of planet Earth. Convergent plate junctions exert enormous influence on the formation and recycling of continental crust, and they are also responsible for major mineral resources and earthquakes, which are of crucial importance for society. A subduction-related geologic unit containing high-pressure rocks occurs in the Barragan area (Valle del Cauca Department) on the western flank of the Central Cordillera of the Colombian Andes. Blueschists and amphibolites, serpentinized meta-ultramafic rocks, graphite-chlorite-muscovite-quartz schists, protocataclasites, and graphite-chlorite-andalusite-andesine-garnet-muscovite +/- titanite schists are exposed in this region. In spite of the petrotectonic importance of blueschists, the high-pressure metamorphism of the Central Cordillera of Colombia has been rarely studied. New geochemical data indicate that protoliths of the blueschist- and amphibolite-facies rocks possessed normal mid-ocean ridge basalt bulk compositions. Ar-40/Ar-39 geochronology for a metapelite rock associated with the blueschists shows a plateau age of similar to 120 million years. We suggest that high-P/T conditions were present from similar to 150 to 125 Ma, depending on the model of generation and exhumation considered.

The rupture process of recent tsunamigenic earthquakes by geophysical data inversion

Subduction zones are the favorite places to generate tsunamigenic earthquakes, where friction between oceanic and continental plates causes the occurrence of a strong seismicity. The topics and the methodologies discussed in this thesis are focussed to the understanding of the rupture process of the seismic sources of great earthquakes that generate tsunamis. The tsunamigenesis is controlled by several kinematical characteristic of the parent earthquake, as the focal mechanism, the depth of the rupture, the slip distribution along the fault area and by the mechanical properties of the source zone. Each of these factors plays a fundamental role in the tsunami generation. Therefore, inferring the source parameters of tsunamigenic earthquakes is crucial to understand the generation of the consequent tsunami and so to mitigate the risk along the coasts. The typical way to proceed when we want to gather information regarding the source process is to have recourse to the inversion of geophysical data that are available. Tsunami data, moreover, are useful to constrain the portion of the fault area that extends offshore, generally close to the trench that, on the contrary, other kinds of data are not able to constrain. In this thesis I have discussed the rupture process of some recent tsunamigenic events, as inferred by means of an inverse method. I have presented the 2003 Tokachi-Oki (Japan) earthquake (Mw 8.1). In this study the slip distribution on the fault has been inferred by inverting tsunami waveform, GPS, and bottom-pressure data. The joint inversion of tsunami and geodetic data has revealed a much better constrain for the slip distribution on the fault rather than the separate inversions of single datasets. Then we have studied the earthquake occurred on 2007 in southern Sumatra (Mw 8.4). By inverting several tsunami waveforms, both in the near and in the far field, we have determined the slip distribution and the mean rupture velocity along the causative fault. Since the largest patch of slip was concentrated on the deepest part of the fault, this is the likely reason for the small tsunami waves that followed the earthquake, pointing out how much the depth of the rupture plays a crucial role in controlling the tsunamigenesis. Finally, we have presented a new rupture model for the great 2004 Sumatra earthquake (Mw 9.2). We have performed the joint inversion of tsunami waveform, GPS and satellite altimetry data, to infer the slip distribution, the slip direction, and the rupture velocity on the fault. Furthermore, in this work we have presented a novel method to estimate, in a self-consistent way, the average rigidity of the source zone. The estimation of the source zone rigidity is important since it may play a significant role in the tsunami generation and, particularly for slow earthquakes, a low rigidity value is sometimes necessary to explain how a relatively low seismic moment earthquake may generate significant tsunamis; this latter point may be relevant for explaining the mechanics of the tsunami earthquakes, one of the open issues in present day seismology. The investigation of these tsunamigenic earthquakes has underlined the importance to use a joint inversion of different geophysical data to determine the rupture characteristics. The results shown here have important implications for the implementation of new tsunami warning systems – particularly in the near-field – the improvement of the current ones, and furthermore for the planning of the inundation maps for tsunami-hazard assessment along the coastal area.

Numerical models of trench migration for lateral heterogeneous subducting plates

The aim of this Thesis is to investigate the effect of heterogeneities within the subducting plate on the dynamics of subduction. In particular, I study the motion of the trench for oceanic and continental subduction, first, separately, and, then, together in the same system to understand how they interact. The understanding of these features is fundamental to reconstruct the evolution of complex subduction zones, such as the Central Mediterranean. For this purpose, I developed 2D and 3D numerical models of oceanic and continental subduction where the rheological, geometrical and compositional properties of the plates are varied. In these models, the trench and the overriding plate move self-consistently as a function of the dynamics of the system. The effect of continental subduction on trench migration is largely investigated. Results from a parametric study showed that despite different rheological properties of the plates, all models with a uniform continental crust share the same kinematic behaviour: the trench starts to advance once the continent arrives at the subduction zone. Hence, the advancing mode in continental collision scenarios is at least partly driven by an intrinsic feature of the system. Moreover, the presence of a weak lower crust within the continental plate can lead to the occurrence of delamination. Indeed, by changing the viscosity of the lower crust, both delamination and slab detachment can occur. Delamination is favoured by a low viscosity value of the lower crust, because this makes the mechanical decoupling easier between crust and lithospheric mantle. These features are observed both in 2D and 3D models, but the numerical results of the 3D models also showed that the rheology of the continental crust has a very strong effect on the dynamics of the whole system, since it influences not only the continental part of plate but also the oceanic sides.

Forward modelling of petrological crust-forming processes on the early Earth

Tonalite-trondhjemite-granodiorite (TTG) gneisses form up to two-thirds of the preserved Archean continental crust and there is considerable debate regarding the primary magmatic processes of the generation of these rocks. The popular theories indicate that these rocks were formed by partial melting of basaltic oceanic crust which was previously metamorphosed to garnet-amphibolite and/or eclogite facies conditions either at the base of thick oceanic crust or by subduction processes.rnThis study investigates a new aspect regarding the source rock for Archean continental crust which is inferred to have had a bulk compostion richer in magnesium (picrite) than present-day basaltic oceanic crust. This difference is supposed to originate from a higher geothermal gradient in the early Archean which may have induced higher degrees of partial melting in the mantle, which resulted in a thicker and more magnesian oceanic crust. rnThe methods used to investigate the role of a more MgO-rich source rock in the formation of TTG-like melts in the context of this new approach are mineral equilibria calculations with the software THERMOCALC and high-pressure experiments conducted from 10–20 kbar and 900–1100 °C, both combined in a forward modelling approach. Initially, P–T pseudosections for natural rock compositions with increasing MgO contents were calculated in the system NCFMASHTO (Na2O–CaO–FeO–MgO–Al2O3–SiO2–H2O–TiO2) to ascertain the metamorphic products from rocks with increasing MgO contents from a MORB up to a komatiite. A small number of previous experiments on komatiites showed the development of pyroxenite instead of eclogite and garnet-amphibolite during metamorphism and established that melts of these pyroxenites are of basaltic composition, thus again building oceanic crust instead of continental crust.rnThe P–T pseudosections calculated represent a continuous development of their metamorphic products from amphibolites and eclogites towards pyroxenites. On the basis of these calculations and the changes within the range of compositions, three picritic Models of Archean Oceanic Crust (MAOC) were established with different MgO contents (11, 13 and 15 wt%) ranging between basalt and komatiite. The thermodynamic modelling for MAOC 11, 13 and 15 at supersolidus conditions is imprecise since no appropriate melt model for metabasic rocks is currently available and the melt model for metapelitic rocks resulted in unsatisfactory calculations. The partially molten region is therfore covered by high-pressure experiments. The results of the experiments show a transition from predominantly tonalitic melts in MAOC 11 to basaltic melts in MAOC 15 and a solidus moving towards higher temperatures with increasing magnesium in the bulk composition. Tonalitic melts were generated in MAOC 11 and 13 at pressures up to 12.5 kbar in the presence of garnet, clinopyroxene, plagioclase plus/minus quartz (plus/minus orthopyroxene in the presence of quartz and at lower pressures) in the absence of amphibole but it could not be explicitly indicated whether the tonalitic melts coexisting with an eclogitic residue and rutile at 20 kbar do belong to the Archean TTG suite. Basaltic melts were generated predominantly in the presence of granulite facies residues such as amphibole plus/minus garnet, plagioclase, orthopyroxene that lack quartz in all MAOC compositions at pressures up to 15 kbar. rnThe tonalitic melts generated in MAOC 11 and 13 indicate that thicker oceanic crust with more magnesium than that of a modern basalt is also a viable source for the generation of TTG-like melts and therefore continental crust in the Archean. The experimental results are related to different geologic settings as a function of pressure. The favoured setting for the generation of early TTG-like melts at 15 kbar is the base of an oceanic crust thicker than existing today or by melting of slabs in shallow subduction zones, both without interaction of tonalic melts with the mantle. Tonalitic melts at 20 kbar may have been generated below the plagioclase stability by slab melting in deeper subduction zones that have developed with time during the progressive cooling of the Earth, but it is unlikely that those melts reached lower pressure levels without further mantle interaction.rn

The Geology of Central Chile's Long-Lived Volcanic Arc

Discussion of subduction zones that are associated with volcanic arcs and major chemical exchanges between the Earth's surface and underlying mantle.

Melting of metasomatized peridotite at 4–6 GPa and up to 1200 °C: an experimental approach

The phase assemblages and compositions in a K-bearing lherzolite + H2O system are determined between 4 and 6 GPa and 850–1200 °C, and the melting reactions occurring at subarc depth in subduction zones are constrained. Experiments were performed on a rocking multi-anvil apparatus. The experiments had around 16 wt% water content, and hydrous melt or aqueous fluid was segregated and trapped in a diamond aggregate layer. The compositions of the aqueous fluid and hydrous melt phases were measured using the cryogenic LA-ICP-MS technique. The residual lherzolite consists of olivine, orthopyroxene, clinopyroxene, and garnet, while diamond (C) is assumed to be inert. Hydrous and alkali-rich minerals were absent from the run products due to preferred dissolution of K2O (and Na2O) to the aqueous fluid/hydrous melt phases. The role of phlogopite in melting relations is, thus, controlled by the water content in the system: at the water content of around 16 wt% used here, phlogopite is unstable and thus does not participate in melting reactions. The water-saturated solidus, i.e., the first appearance of hydrous melt in the K–lherzolite composition, is located between 900 and 1000 °C at 4 GPa and between 1000 and 1100 °C at 5 and 6 GPa. Compositional jumps between hydrous melt and aqueous fluid at the solidus include a significant increase in the total dissolved solids load. All melts/fluids are peralkaline and calcium-rich. The melting reactions at the solidus are peritectic, as olivine, clinopyroxene, garnet, and H2O are consumed to generate hydrous melt plus orthopyroxene. Our fluid/melt compositional data demonstrate that the water-saturated hybrid peridotite solidus lies above 1000 °C at depths greater than 150 km and that the second critical endpoint is not reached at 6 GPa for a K2O–Na2O–CaO–FeO–MgO–Al2O3–SiO2–H2O–Cr2O3(–TiO2) peridotite composition.

Geochemical composition of the Trincheira Complex in the southwestern Amazon craton, Brazil

We document the first-known Mesoproterozoic ophiolite from the southwestern part of the Amazon craton, corresponding to the Trincheira Complex of Calymmian age, and propose a tectonic model that explains many previously enigmatic features of the Precambrian history of this key craton, and discuss its role in the reconstruction of the Columbia supercontinent. The complex comprises extrusive rocks (fine-grained amphibolites derived from massive and pillowed basalts), mafic-ultramafic intrusive rocks, chert, banded iron formation (BIFs), pelites, psammitic and a smaller proportion of calc-silicate rocks. This sequence was deformed, metasomatized and metamorphosed during the development of the Alto Guaporé Belt, a Mesoproterozoic accretionary orogen. The rocks were deformed by a single tectonic event, which included isoclinal folding and metamorphism of the granulite-amphibolite facies. Layered magmatic structures were preserved in areas of low strain, including amygdaloidal and cumulate structures. Metamorphism was pervasive and reached temperatures of 780-853°C in mafic granulites and 680-720°C in amphibolites under an overall pressure of 6.8 kbar. The geochemical composition of the extrusive and intrusive rocks indicates that all noncumulus mafic-ultramafic rocks are tholeiitic basalts. The mafic-ultramafic rocks display moderately to strongly fractionation of light rare earth elements (LREE), near-flat heavy rare earth elements (HREE) patterns and moderate to strong negative high field strength elements (HFSE) anomalies (especially Nb), a geochemical signature typical of subduction zones. The lowest units of mafic granulites and porphyroblastic amphibolites in the Trincheira ophiolite are similar to the modern mid-ocean ridge basalt (MORB), although they locally display small Ta, Ti and Nb negative anomalies, indicating a small subduction influence. This behavior changes to an island arc tholeiites (IAT) signature in the upper units of fine-grained amphibolites and amphibole rich-amphibolites, characterized by progressive depletion in the incompatible elements and more pronounced negative Ta and Nb anomalies, as well as common Ti and Zr negative anomalies. Tectono-magmatic variation diagrams and chondrite-normalized REE and primitive mantle normalized patterns suggest a back-arc to intra-oceanic island arc tectonic regime for the eruption of these rocks. Therefore, the Trincheira ophiolite appears to have originated in an intraoceanic supra-subduction setting composed of an arc-back-arc system. Accordingly, the Trincheira Complex is a record of oceanic crust relics obducted during the collision of the Amazon craton and the Paraguá block during the Middle Mesoproterozoic. Thus, the recognition of the Trincheira ophiolite and suture significantly changes views on the evolution of the southern margin of the Amazon craton, and how it can influence the global tectonics and the reconstruction of the continents.

Trace element concentrations and Li isotope composition in DSDP and ODP sediments from island arcs

We have studied the sedimentary and basaltic inputs of lithium to subduction zones. Various sediments from DSDP and ODP drill cores in front of the Mariana, South Sandwich, Banda, East Sunda and Lesser Antilles island arcs have been analysed and show highly variable Li contents and d7Li values. The sediment piles in front of the Mariana and South Sandwich arcs largely consist of pelagic sediments (clays and oozes). The pelagic clays have high Li contents (up to 57.3 ppm) and Li isotope compositions ranging from +1.3? to +4.1?. The oozes have lower Li contents (7.3-16 ppm) with d7Li values of the diatom oozes from the South Sandwich lower (+2.8? to +3.2?) than those of the radiolarian oozes from the Mariana arc (+8.1? to +14.5?). Mariana sediment also contains a significant portion of volcanogenic material, which is characterised by a moderate Li content (14 ppm) and a relatively heavy isotope composition (+6.4?). Sediments from the Banda and Lesser Antilles contain considerable amounts of continental detritus, and have high Li contents (up to 74.3 ppm) and low d7Li values (around 0?), caused by weathering of continental bedrock. East Sunda sediments largely consist of calcareous oozes. These carbonate sediments display intermediate to high Li contents (2.4-41.9 ppm) and highly variable d7Li values (-1.6? to +12.8?). Basaltic oceanic crust samples from worldwide DSDP and ODP drill cores are characterised by enrichment of Li compared to fresh MORB (6.6-33.1 vs. 3.6-7.5 ppm, respectively), and show a large range in Li isotope compositions (+1.7? to +11.8?). The elemental and isotopic enrichment of Li in altered basalts is due to the uptake of isotopically heavy seawater Li during weathering. However, old oceanic crust samples from Sites 417/418 exhibit lighter Li isotope compositions compared to young basaltic crust samples from Sites 332B and 504B. This lighter Li isotope signature in old crust is unexpected and further research is needed to explore this issue.