Ocean bottom seismic noise : applications for the crust knowledge, interaction ocean-atmosphere and instrumental behaviour

Tese de doutoramento, Ciências Geofísicas e da Geoinformação (Geofisíca), Universidade de Lisboa, Faculdade de Ciências, 2014

Assinatura tectono-sedimentar do auge da compressão bética em Portugal: a descontinuidade sedimentar Valesiano terminal-Turoliano

An Upper Miocene important sedimentary break can be accurately recognised in the Portuguese basins and is reflected by a drastic palaeogeographic change in relation to a large-scale tectonic event of probable uppermost Vallesian-Turolian (9,5 Ma; middle Tortonian) age. The characterisation of the sedimentary record of this tectonic event, as well as its relations with interpreted active faults is made for different situations: Douro (NW border), Mondego, Lower-Tagus and Sado Tertiary basins. The sedimentary record, considered upper Tortonian-Messinian ? (uppermost Vallesian-Turolian ?) is interpreted mainly as endorheic alluvial fans (internal drainage), developed along active NNE-SSW indent-linked strike-slip faults and NE-SW reverse faults. At NE Portugal, proximal fluvial systems of an endorheic hydrographic network drained eastwards to the Spanish Duero interior Basin. The main evidences of the betic compression clímax in Portugal mainland are presented; the interpreted active tectonic structures are in accordance with an intense NNW-SSE crustal shortening, but some regional differences are also documented.

The role of a transcrustal shear zone in orogenic gold mineralization at the Aijanahalli Mine, Dharwar Craton, South India

The Ajjanahalli gold mine is spatially associated with a Late Archean craton-scale shear zone in the eastern Chitradurga greenstone belt of the Dharwar craton, India. Gold mineralization is hosted by an similar to100-m-wide antiform in a banded iron formation. Original magnetite and siderite are replaced by a peak metamorphic alteration assemblage of chlorite, stilpnomelane, minnesotaite, sericite, ankerite, arsenopyrite, pyrite, pyrrhotite, and gold at ca. 300degrees to 350degreesC. Elements enriched in the banded iron formation include Ca, Mg, C, S, An, As, Bi. Cu, Sb, Zn, Pb, Se, Ag, and Te, whereas in the wall rocks As, Cu, Zn, Bi, Ag, and An are only slightly enriched. Strontium correlates with CaO, MgO, CO2, and As, which indicates cogenetic formation of arsenopyrite and Mg-Ca carbonates. The greater extent of alteration in the Fe-rich banded iron formation layers than in the wall rock reflects the greater reactivity of the banded iron formation layers. The ore fluids, as interpreted from their isotopic composition (delta(18)O = 6.5-8.5parts per thousand; initial Sr-87/Sr-86 = 0.7068-0.7078), formed by metamorphic devolatilization of deeper levels of the Chitradurga greenstone belt. Arsenopyrite, chalcopyrite, and pyrrhotite have delta(34)S values within a narrow range between 2.1 and 2.7 per mil, consistent with a sulfur source in Chitradurga greenstone belt lithologies. Based on spatial and temporal relationships between mineralization, local structure development, and sinistral strike-slip deformation in the shear zone at the eastern contact of the Chitradurga greenstone belt, we suggest that the Ajjanahalli gold mineralization formed by fluid infiltration into a low strain area within the first-order structure. The ore fluids were transported along this shear zone into relatively shallow crustal levels during lateral terrane accretion and a change from thrust to transcurrent tectonics. Based on this model of fluid flow, exploration should focus on similar low strain areas or potentially connected higher order splays of the first-order shear zone.

Depleted arc volcanism in the Alboran Sea and shoshonitic volcanism in Morocco: geochemical and isotopic constraints on Neogene tectonic processes

Samples of volcanic rocks from Alboran Island, the Alboran Sea floor and from the Gourougou volcanic centre in northern Morocco have been analyzed for major and trace elements and Sr-Nd isotopes to test current theories on the tectonic geodynamic evolution of the Alboran Sea. The Alboran Island samples are low-K tholeiitic basaltic andesites whose depleted contents of HFS elements (similar to0.5xN-MORB), especially Nb (similar to0.2xN-MORB), show marked geochemical parallels with volcanics from immature intra-oceanic arcs and back-arc basins. Several of the submarine samples have similar compositions, one showing low-Ca boninite affinity. Nd-143/Nd-144 ratios fall in the same range as many island-arc and back-arc basin samples, whereas Sr-87/Sr-86 ratios (on leached samples) are somewhat more radiogenic. Our data point to active subduction taking place beneath the Alboran region in Miocene times, and imply the presence of an associated back-arc spreading centre. Our sea floor suite includes a few more evolved dacite and rhyolite samples with (Sr-87/Sr-86)(0) up to 0.717 that probably represent varying degrees of crustal melting. The shoshonite and high-K basaltic andesite lavas from Gourougou have comparable normalized incompatible-element enrichment diagrams and Ce/Y ratios to shoshonitic volcanics from oceanic island arcs, though they have less pronounced Nb deficits. They are much less LIL- and LREE-enriched than continental arc analogues and post-collisional shoshonites from Tibet. The magmas probably originated by melting in subcontinental lithospheric mantle that had experienced negligible subduction input. Sr-Nd isotope compositions point to significant crustal contamination which appears to account for the small Nb anomalies. The unmistakable supra-subduction zone (SSZ) signature shown by our Alboran basalts and basaltic andesite samples refutes geodynamic models that attribute all Neogene volcanism in the Alboran domain to decompression melting of upwelling asthenosphere arising from convective thinning of over-thickened lithosphere. Our data support recent models in which subsidence is caused by westward rollback of an eastward-dipping subduction zone beneath the westemmost Mediterranean. Moreover, severance of the lithosphere at the edges of the rolling-back slab provides opportunities for locally melting lithospheric mantle, providing a possible explanation for the shoshonitic volcanism seen in northern Morocco and more sporadically in SE Spain. (C) 2004 Elsevier B.V. All rights reserved.

Modeling of in-situ crystallization processes in the Permian mafic layered intrusion of Mont Collon (Dent Blanche nappe, western Alps)

The Mont Collon mafic complex is one of the best preserved examples of the Early Permian magmatism in the Central Alps, related to the intra-continental collapse of the Variscan belt. It mostly consists (> 95 vol.%) of ol+hy-nonnative plagioclase-wehrlites, olivine- and cpx-gabbros with cumulitic structures, crosscut by acid dikes. Pegmatitic gabbros, troctolites and anorthosites outcrop locally. A well-preserved cumulative, sequence is exposed in the Dents de Bertol area (center of intrusion). PT-calculations indicate that this layered magma chamber emplaced at mid-crustal levels at about 0.5 GPa and 1100 degrees C. The Mont Collon cumulitic rocks record little magmatic differentiation, as illustrated by the restricted range of clinopyroxene mg-number (Mg#(cpx)=83-89). Whole-rock incompatible trace-element contents (e.g. Nb, Zr, Ba) vary largely and without correlation with major-element composition. These features are characteristic of an in-situ crystallization process with variable amounts of interstitial liquid L trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distribution in the latter is homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in-situ crystallization equation successfully duplicated the REE concentrations in cumulitic minerals of all rock facies of the intrusion. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 20%, relative to the least evolved facies of the intrusion. L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. However, the in-situ crystallization model reaches its limitations with rock containing high modal content of REE-bearing minerals (i.e. zircon), such as pegmatitic gabbros. Dikes of anorthositic composition, locally crosscutting the layered lithologies, evidence that the Mont Collon rocks evolved in open system with mixing of intercumulus liquids of different origins and possibly contrasting compositions. The proposed model is not able to resolve these complex open systems, but migrating liquids could be partly responsible for the observed dispersion of points in some correlation diagrams. Absence of significant differentiation with recurrent lithologies in the cumulitic pile of Dents de Bertol points to an efficiently convective magma chamber, with possible periodic replenishment, (c) 2005 Elsevier B.V. All rights reserved.

An interpretation of the deep seismic lines from the Penninic Alps of Valais (Switzerland)

The upper part of three deep seismic lines running across the Penninic Swiss Alps of Valais have been studied. Numerous reflectors illustrate the nappe structure of this internal part of the orogen. These reflectors, even at great depths (20-25 km), can be correlated with outcropping geological features and are most likely produced by lithological boundaries rather than by mylonites zones, which are hardly reflective in such an environment. Our interpretations, largely constrained by projections of the outcropping geology, have improved our knowledge of the deep structure of this segment of the Alpine belt, enhancing the importance of the backfolding and the crustal scale deformation phase which produced the Rawil-Valpelline depression and the Aar-Toce culmination. Furthermore we have here the possibility of correlating seismic patterns produced by ductile folds with the outcropping structures.

Petrography, Petrochemistry and petrogenesis of Huronian volcanic rocks of the Elliot Lake region, Ontario

In the Elliot lake region of northern Ontario, Yolcanlc lava piles represent the lowermost units of the Huronian SUpergroup. These rocks general1y trend east-west and belong to the Elliot lake Group. They are s1tuated on the north and south limbs or the QuIrke lake Syncline. The volcanIc rocks of this study contain a secondary minerai assemblage consisting of actinolite, biotite, chlorIte, eptdote/cllnozoislte tttanomagnettte and calcite characteristic of greenschist metamorphism. Compilation of data suggests that metamorphism of the volcanic rocks proceeded between 325- and 425-C and between 2.4 and 4.7 kb. Geochemtcally these lavas represent tholeiitic and calc-alkaline assemblages. The tholeiites are character1sttcally enriched tn Fe and Tt and consist mainly of basalts, basaltic andesites and andesites. These rocks are believed to have formed by the partIal melting of a peridottte source at low P-T. In contrast, the calc-alkaline rocks are depleted in Fe and TI, but show a signIficant enrichment In 51 and Zr; andesIte Is the major rock type for thIs assemblage. I·t Is postUlated that the calc-alkalIne sU1te of rocks was the result of eIther the partIal meltIng of abasaltic·magma at shallow depth, or the melttng of s1al1c crustal materIal due to the added we1ght of tholeiitIc material on an unstable crust and to downwarplng processes Inttlated by convection cells.

The sedimentology of the lower Silurian whirlpool sandstone in subsurface Lake Erie, Ontario

The lower Silurian Whirlpool Sandstone is composed of two main units: a fluvial unit and an estuarine to transitional marine unit. The lowermost unit is made up of sandy braided fluvial deposits, in shallow valleys, that flowed towards the northwest. The fluvial channels are largely filled by cross-bedded, well sorted, quartzose sands, with little ripple crosslaminated or overbank shales. Erosionally overlying this lower unit are brackish water to marine deposits. In the east, this unit consists of estuarine channels and tidal flat deposits. The channels consist of fluvial sands at the base, changing upwards into brackish and tidally influenced channelized sandstones and shales. The estuarine channels flowed to the southwest. Westwards, the unit contains backbarrier facies with extensive washover deposits. Separating the backbarrier facies from shoreface sandstone facies to the west, are barrier island sands represented by barrier-foreshore facies. The barrier islands are dissected by tidal inlets characterized by fining upward abandonment sequences. Inlet deposits are also present west of the barrier island, abandoned by transgression on the shoreface. The sandy marine deposits are replaced to the west by carbonates of the Manitoulin Limestone. During the latest Ordovician, a hiatus in crustal loading during the Taconic Orogeny led to erosional offloading and crustal rebound, the eroded material distributed towards the west, northwest and north as the terrestrial deposits of the fluvial Whirlpool. The "anti-peripheral bulge" of the rebound interfered with the peripheral bulge of the Michigan Basin, nulling the Algonquin Arch, and allowing the detritus of the fluvial Whirlpool to spread onto the Algonquin Arch. The Taconic Orogeny resumed in the earliest Silurian with crustal loading to the south and southeast, and causing tilting of the surface slope in subsurface Lake Erie towards the ii southwest. Lowstand terrestrial deposits were scoured into the new slope. The new crustal loading also reactivated the peripheral bulge of the Appalachian Basin, allowing it to interact with the bulge of the Michigan Basin, raising the Algonquin Arch. The crustal loading depressed the Appalachian basin and allowed transgression to occur. The renewed Algonquin Arch allowed the early Silurian transgression to proceed up two slopes, one to the east and one to the west. The transgression to the east entered the lowstand valleys and created the estuarine Whirlpool. The rising arch caused progradation of the Manitoulin carbonates upon shoreface facies of the Whirlpool Sandstone and upon offshore facies of the Cabot Head Formation. Further crustal loading caused basin subsidence and rapid transgression, abandoning the Whirlpool estuary in an offshore setting.

Geophysical Constraints on Structure and Tectonics of the Eastern Arabian Sea and the Adjoining West Coast of India with Special Reference to the Kerala Basin

The main objective of the present study is to model the gravity fields in terms of lithospheric structure below the western continental margin of India (WCMI) identify zones of crustal mass anomalies and attempt to infer the location of Ocean Continent transition in the Arabian Sea. In this study, the area starting from the western shield margin to the region covering the deep oceanic parts of the Arabian Sea which is bounded by Carlsberg and Cerg and Central Indian ridges in the south, eastern part of the Indus Cone in the west and falling between 630E and 800E longitudes, and 50N - 240N latitudes has been considered. The vast amount of seismic reflection and refraction data in the form of crustal velocities, basement configuration and crustal thicknesses available for the west coast as well as the eastern Arabian Sea has been utilized for this purpose

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.

Geodynamics of the Andaman - Sumatra - Java Trench - Arc System Based on Gravity and Seismotectonic Study

This work aims to study the variation in subduction zone geometry along and across the arc and the fault pattern within the subducting plate. Depth of penetration as well as the dip of the Benioff zone varies considerably along the arc which corresponds to the curvature of the fold- thrust belt which varies from concave to convex in different sectors of the arc. The entire arc is divided into 27 segments and depth sections thus prepared are utilized to investigate the average dip of the Benioff zone in the different parts of the entire arc, penetration depth of the subducting lithosphere, the subduction zone geometry underlying the trench, the arctrench gap, etc.The study also describes how different seismogenic sources are identified in the region, estimation of moment release rate and deformation pattern. The region is divided into broad seismogenic belts. Based on these previous studies and seismicity Pattern, we identified several broad distinct seismogenic belts/sources. These are l) the Outer arc region consisting of Andaman-Nicobar islands 2) the back-arc Andaman Sea 3)The Sumatran fault zone(SFZ)4)Java onshore region termed as Jave Fault Zone(JFZ)5)Sumatran fore arc silver plate consisting of Mentawai fault(MFZ)6) The offshore java fore arc region 7)The Sunda Strait region.As the Seismicity is variable,it is difficult to demarcate individual seismogenic sources.Hence, we employed a moving window method having a window length of 3—4° and with 50% overlapping starting from one end to the other. We succeeded in defining 4 sources each in the Andaman fore arc and Back arc region, 9 such sources (moving windows) in the Sumatran Fault zone (SFZ), 9 sources in the offshore SFZ region and 7 sources in the offshore Java region. Because of the low seismicity along JFZ, it is separated into three seismogenic sources namely West Java, Central Java and East Java. The Sunda strait is considered as a single seismogenic source.The deformation rates for each of the seismogenic zones have been computed. A detailed error analysis of velocity tensors using Monte—Carlo simulation method has been carried out in order to obtain uncertainties. The eigen values and the respective eigen vectors of the velocity tensor are computed to analyze the actual deformation pattem for different zones. The results obtained have been discussed in the light of regional tectonics, and their implications in terms of geodynamics have been enumerated.ln the light of recent major earthquakes (26th December 2004 and 28th March 2005 events) and the ongoing seismic activity, we have recalculated the variation in the crustal deformation rates prior and after these earthquakes in Andaman—Sumatra region including the data up to 2005 and the significant results has been presented.ln this chapter, the down going lithosphere along the subduction zone is modeled using the free air gravity data by taking into consideration the thickness of the crustal layer, the thickness of the subducting slab, sediment thickness, presence of volcanism, the proximity of the continental crust etc. Here a systematic and detailed gravity interpretation constrained by seismicity and seismic data in the Andaman arc and the Andaman Sea region in order to delineate the crustal structure and density heterogeneities a Io nagnd across the arc and its correlation with the seismogenic behaviour is presented.

Modelling of Mercury’s surface composition and remote detection from the orbit with the BepiColombo Mercury Planetary Orbiter

It can be assumed that the composition of Mercury’s thin gas envelope (exosphere) is related to the composition of the planets crustal materials. If this relationship is true, then inferences regarding the bulk chemistry of the planet might be made from a thorough exospheric study. The most vexing of all unsolved problems is the uncertainty in the source of each component. Historically, it has been believed that H and He come primarily from the solar wind, while Na and K originate from volatilized materials partitioned between Mercury’s crust and meteoritic impactors. The processes that eject atoms and molecules into the exosphere of Mercury are generally considered to be thermal vaporization, photonstimulated desorption (PSD), impact vaporization, and ion sputtering. Each of these processes has its own temporal and spatial dependence. The exosphere is strongly influenced by Mercury’s highly elliptical orbit and rapid orbital speed. As a consequence the surface undergoes large fluctuations in temperature and experiences differences of insolation with longitude. We will discuss these processes but focus more on the expected surface composition and solar wind particle sputtering which releases material like Ca and other elements from the surface minerals and discuss the relevance of composition modelling

Ground deformation at Soufriere Hills Volcano, Montserrat during 1998-2000 measured by radar interferometry and GPS

We examine the motion of the ground surface on the Soufriere Hills Volcano, Montserrat between 1998 and 2000 using radar interferometry (InSAR). To minimise the effects of variable atmospheric water vapour on the InSAR measurements we use independently-derived measurements of the radar path delay from six continuous GPS receivers. The surfaces providing a measurable inter-ferometric signal are those on pyroclastic flow deposits, mainly emplaced in 1997. Three types of surface motion can be discriminated. Firstly, the surfaces of thick, valley-filling deposits subsided at rates of 150-120 mm/year in the year after emplacement to 50-30 mm/year two years later. This must be due to contraction and settling effects during cooling. The second type is the near-field motion localised within about one kilometre of the dome. Both subsidence and uplift events are seen and though the former could be due to surface gravitational effects, the latter may reflect shallow (< 1 km) pressurisation effects within the conduit/dome. Far-field motions of the surface away from the deeply buried valleys are interpreted as crustal strains. Because the flux of magma to the surface stopped from March 1998 to November 1999 and then resumed from November 1999 through 2000, we use InSAR data from these two periods to test the crustal strain behaviour of three models of magma supply: open, depleting and unbalanced. The InSAR observations of strain gradients of 75-80 mm/year/krn uplift during the period of quiescence on the western side of the volcano are consistent with an unbalanced model in which magma supply into a crustal magma chamber continues during quiescence, raising chamber pressure that is then released upon resumption of effusion. GPS motion vectors agree qualitatively with the InSAR displacements but are of smaller magnitude. The discrepancy may be due to inaccurate compensation for atmospheric delays in the InSAR data. (c) 2005 Elsevier B.V. All rights reserved.

A stress-controlled mechanism for the intensity of very large magnitude explosive eruptions

Large magnitude explosive eruptions are the result of the rapid and large-scale transport of silicic magma stored in the Earth's crust, but the mechanics of erupting teratonnes of silicic magma remain poorly understood. Here, we demonstrate that the combined effect of local crustal extension and magma chamber overpressure can sustain linear dyke-fed explosive eruptions with mass fluxes in excess of 10^10 kg/s from shallow-seated (4–6 km depth) chambers during moderate extensional stresses. Early eruption column collapse is facilitated with eruption duration of the order of few days with an intensity of at least one order of magnitude greater than the largest eruptions in the 20th century. The conditions explored in this study are one way in which high mass eruption rates can be achieved to feed large explosive eruptions. Our results corroborate geological and volcanological evidences from volcano-tectonic complexes such as the Sierra Madre Occidental (Mexico) and the Taupo Volcanic Zone (New Zealand).

Stress field control of eruption dynamics at a rift volcano: Nyamuragira, D.R.Congo

Using the record of 30 flank eruptions over the last 110 years at Nyamuragira, we have tested the relationship between the eruption dynamics and the local stress field. There are two groups of eruptions based on their duration (< 80days >) that are also clustered in space and time. We find that the eruptions fed by dykes parallel to the East African Rift Valley have longer durations (and larger volumes) than those eruptions fed by dykes with other orientations. This is compatible with a model for compressible magma transported through an elastic-walled dyke in a differential stress field from an over-pressured reservoir (Woods et al., 2006). The observed pattern of eruptive fissures is consistent with a local stress field modified by a northwest-trending, right lateral slip fault that is part of the northern transfer zone of the Kivu Basin rift segment. We have also re-tested with new data the stochastic eruption models for Nyamuragira of Burt et al. (1994). The time-predictable, pressure-threshold model remains the best fit and is consistent with the typically observed declining rate of sulphur dioxide emission during the first few days of eruption with lava emission from a depressurising, closed, crustal reservoir. The 2.4-fold increase in long-term eruption rate that occurred after 1977 is confirmed in the new analysis. Since that change, the record has been dominated by short-duration eruptions fed by dykes perpendicular to the Rift. We suggest that the intrusion of a major dyke during the 1977 volcano-tectonic event at neighbouring Nyiragongo volcano inhibited subsequent dyke formation on the southern flanks of Nyamuragira and this may also have resulted in more dykes reaching the surface elsewhere. Thus that sudden change in output was a result of a changed stress field that forced more of the deep magma supply to the surface. Another volcano-tectonic event in 2002 may also have changed the magma output rate at Nyamuragira.