Multi-Disciplinary Applications of Oceanographic Geophysical Data Collection

Thesis (Ph.D.)--University of Washington, 2016-08

U-Pb SHRIMP zircon dating of high-grade rocks from the Upper Allochthonous Terrane of Bragança and Morais Massifs (NE Portugal); geodynamic consequences

Bragança and Morais Massifs are part of the mega-klippen ensemble of NW Iberia, comprising a tectonic pile of four allochthonous units stacked above the Central-Iberian Zone autochthon. On top of this pile, the Upper Allochthonous Terrane (UAT) includes different high-grade metamorphic series whose age and geodynamic meaning are controversial. Mafic granulites provided U–Pb zircon ages at 399±7 Ma, dating the Variscan emplacement of UAT. In contrast,U–Pb zircon ages of ky- and hb-eclogites, felsic/intermediate HP/HT-granulites and orthogneisses (ca. 500–480 Ma) are identical to those of gabbros (488 ± 10 Ma) and Grt-pyroxenites (495 ± 8 Ma) belonging to a mafic/ultramafic igneous suite that records upper mantle melting and mafic magma crustal underplating at these times. Gabbros intrude the high-grade units of UAT and did not underwent the HP metamorphic event experienced by eclogites and granulites. These features and the zircon dates resemblance among different lithologies, suggest that extensive age resetting of older events may have been correlative with the igneous suite emplacement/crystallisation. Accordingly, reconciliation of structural, petrological and geochronological evidence implies that the development and early deformation of UAT high-grade rocks should be ascribed to an orogenic cycle prior to ≈500 Ma. Undisputable dating of this cycle is impossible, but the sporadic vestiges of Cadomian ages cannot be disregarded. The ca. 500–480 Ma time-window harmonises well with the Lower Palaeozoic continental rifting that trace the VariscanWilson Cycle onset and the Rheic Ocean opening. Subsequent preservation of the high heat-flowregime, possibly related to the Palaeotethys back-arc basin development (ca. 450–420 Ma), would explain the 461 ± 10 Ma age yielded by some zircon domains in felsic granulites, conceivably reflecting zircon dissolution/ recrystallisation till Ordovician times, long before the Variscan paroxysm (ca. 400–390 Ma). This geodynamic scenario suggests also that UAT should have been part of Armorica before its emplacement on top of Iberia after Palaeotethys closure.

Geoquímica dos filões máficos de Lagoa (Maciço de Morais, Trás-os-Montes, NE Portugal)

No Maciço de Morais, os gnaisses de Lagoa são intruídos por filões máficos, de idade precâmbrica-câmbrica, que cortam um "fabric" metamórfico precoce naqueles gnaisses. Os filões máficos de Lagoa têm características geoquímicas de toleítos continentais e testemunham processos ACF. A actividade magmática em que esses filões se inserem reflecte ascenção astenosférica subsequente ao sobre-espessamento crustal precâmbrico, correspondendo a uma manifestação precoce do regime tectónico/magmático que irá caracterizar, durante o Paleozóico inferior, o ciclo varisco ibérico.

Thermal regime of the Southeast Indian Ridge between 88°E and 140°E: Remarks on the subsidence of the ridge flanks

The flanks of the Southeast Indian Ridge are characterized by anomalously low subsidence rates for the 0-25 Ma period: less than 300 m Ma(-1/2) between 101 degrees E and 120 degrees E and less than 260 m Ma(-1/2) within the Australian-Antarctic Discordance (AAD), between 120 degrees E and 128 degrees E. The expected along-axis variation in mantle temperature (similar to 50 degrees C) is too small to explain this observation, even when the temperature dependence of the mantle physical properties is accounted for. We successively analyze the effect on subsidence of different factors, such as variations in crustal thickness; the dynamic contribution of an old, detached slab supposedly present within the mantle below the AAD; and depletion in phi(m), a parameter here defined as the "ubiquitously distributed melt fraction" within the asthenosphere. These effects may all contribute to the observed, anomalously low subsidence rate of the ridge flanks, with the most significant contribution being probably related to the depletion in phi(m). However, these effects have a deep-seated origin that cannot explain the abruptness of the transition across the fracture zones that delineate the boundaries of the AAD, near 120 degrees E and near 128 degrees E, respectively.

The multistage crystallization of zircon in calc-alkaline granitoids : U-Pb age constraints on the timing of Variscan tectonic activity in SW Iberia

CL imaging and U–Th–Pb data for a population of zircons from two of the Évora Massif granitoids (Ossa-Morena Zone, SW Iberia) show that both calc-alkaline granitoids have zircon populations dominated by grains with cores and rims either showing or not showing differences in Th/U ratio, and having ages in the range ca. 350–335 Ma (Early Carboniferous). Multistage crystallization of zircon is revealed in two main growth stages (ca. 344–342 Ma and ca. 336–335 Ma), well represented by morphologically complex zircons with cores and rims with different ages and different Th/U ratios that can be explained by: (1) crystallization from melts with different compositions (felsic peraluminous to felsic-intermediate metaluminous; 0.001 Th/U ratio < 0.5) and (2) transient temperature fluctuations in a system where anatectic felsic melts periodically underwent injection of more mafic magmas at higher temperatures. The two studied calc-alkaline granitoids do not include inherited zircons (pre-Carboniferous), probably because they were formed at the highest grade of metamorphism (T 837 °C; granulite facies) and/or because they were derived from inheritance-poor felsic and mafic rocks from a previous cycle, as suggested by the internal structures of zircon cores. These Variscan magmatic rocks with crystallization ages estimated at ca. 336–335 Ma are spatially and temporally related to high-temperature metamorphism, anatexis, processes of interaction between crustal- and mantle-derived magmas and intra-orogenic extension that acted in SW Iberia during the Early Carboniferous.

Post-rift evolution of the Gulf of Lion margin tested by stratigraphic modelling

The sedimentary architecture of basins and passive margins is determined by a complex interaction of parameters, including subsidence, eustasy, and sediment supply. A quantification of the post-rift (20 Ma-0 Ma) vertical movements of the Gulf of Lion (West Mediterranean) is proposed here based on the stratigraphic study of sedimentary paleomarkers using a large 3D grid of reflection seismic data, correlations with existing drillings, and refraction data. Post-rift subsidence was measured by the direct use of sedimentary geometries analysed in 3D and validated by numerical stratigraphic modelling. Three domains of subsidence were found: on the continental shelf and slope, subsidence corresponds to a seaward tilting with different amplitudes, whereas the deep basin subsides purely vertically. We show that these domains fit with the deeper crustal domains highlighted by previous geophysical data, and that post-break-up subsidence follows the initial hinge lines of the rifting phase. Subsidence rates are quantified on each domain for each stratigraphic interval. At a constant distance from the rotational hinge line, the Plio-Quaternary subsidence rate is constant on the shelf overall. Conversely, Miocene subsidence rates are very different on the eastern and western shelves. Stratigraphic simulations focused on the Messinian salinity crisis (MSC) were also performed. Their results are discussed together with our post-rift subsidence estimates in order to provide ideas and hypotheses for future detailed quantifications of Miocene subsidence, including isostatic readjustments linked to the MSC.

Quantifying uplift and denudation of a thermally heterogeneous crust: a detailed multi-thermochronometeric study of central west Britain

Topography is often thought as exclusively linked to mountain ranges formed by plates collision. It is now, however, known that apart from compression, uplift and denudation of rocks may be triggered by rifting, like it happens at elevated passive margins, and away from plate boundaries by both intra-plate stress causing reactivation of older structures, and by epeirogenic movements driven by mantle dynamics and initiating long-wavelength uplift. In the Cenozoic, central west Britain and other parts of the North Atlantic margins experienced multiple episodes of rock uplift and denudation that have been variable both at spatial and temporal scales. The origin of topography in central west Britain is enigmatic, and because of its location, it may be related to any of the processes mentioned above. In this study, three low temperature thermochronometers, the apatite fission track (AFT) and apatite and zircon (U-Th-Sm)/He (AHe and ZHe, respectively) methods were used to establish the rock cooling history from 200◦C to 30◦C. The samples were collected from the intrusive rocks in the high elevation, high relief regions of the Lake District (NW England), southern Scotland and northern Wales. AFT ages from the region are youngest (55–70 Ma) in the Lake District and increase northwards into southern Scotland and southwards in north Wales (>200 Ma). AHe and ZHe ages show no systematic pattern; the former range from 50 to 80 Ma and the latter tend to record the post-emplacement cooling of the intrusions (200–400 Ma). The complex, multi-thermochronometric inverse modelling suggests a ubiquitous, rapid Late Cretaceous/early Palaeogene cooling event that is particularly marked in Lake District and Criffell. The timing and rate of cooling in southern Scotland and in northern Wales is poorly resolved as the amount of cooling was less than 60◦C. The Lake District plutons were at >110◦C prior to the early Palaeogene; cooling due to a combined effect of high heat flow, from the heat producing granite batholith, and the blanketing effect of the overlying low conductivity Late Mesozoic limestones and mudstones. Modelling of the heat transfer suggests that this combination produced an elevated geothermal gradient within the sedimentary rocks (50–70◦C/km) that was about two times higher than at the present day. Inverse modelling of the AFT and AHe data taking the crustal structure into consideration suggests that denudation was the highest, 2.0–2.5 km, in the coastal areas of the Lake District and southern Scotland, gradually decreasing to less than 1 km in the northern Southern Uplands and northern Wales. Both the rift-related uplift and the intra-plate compression poorly correlate with the timing, location and spatial distribution of the early Palaeogene denudation. The pattern of early Palaeogene denudation correlates with the thickness of magmatic underplating, if the changes of mean topography, Late Cretaceous water depth and eroded rock density are taken into consideration. However, the uplift due to underplating alone cannot fully justify the total early Palaeogene denudation. The amount that is not ex- plained by underplating is, however, roughly spatially constant across the study area and can be referred to the transient thermal uplift induced by the mantle plume arrival. No other mechanisms are required to explain the observed pattern of denudation. The onset of denudation across the region is not uniform. Denudation started at 70–75 Ma in the central part of the Lake District whereas the coastal areas the rapid erosion appears to have initiated later (65–60 Ma). This is ~10 Ma earlier than the first vol- canic manifestation of the proto-Iceland plume and favours the hypothesis of the short period of plume incubation below the lithosphere before the volcanism. In most of the localities, the rocks had cooled to temperatures lower than 30◦C by the end of the Palaeogene, suggesting that the total Neogene denudation was, at a maximum, several hundreds of metres. Rapid cooling in the last 3 million years is resolved in some places in southern Scotland, where it could be explained by glacial erosion and post-glacial isostatic uplift.

Mafic and felsic plutonism in the Elvas region (SW Iberia, Portugal): two distinct Nd isotope sources and geodynamics

The occurrence of mafic (mainly gabbros and diorites) and felsic (syenites and granites) rocks, in close spatial association, in the Elvas region, at the northern part of the Ossa-Morena Zone, could be interpreted as a single bimodal (alkaline) plutonic complex. However, in spite of scarce isotopic (Sm-Nd) data, the co-magmatic origin of both rock groups (mafic and felsic) has already been questioned [1]. Based on the mineral chemistry of primary clinopyroxenes (Di–Hd, %En: 45.5 – 27.2) and representative whole-rock analyses, gabbros and diorites of the Elvas massif show a transitional character between alkaline and non-alkaline fields and wide compositions: SiO2 (42.47 – 58.00 wt%); TiO2 (0.24 – 1.68 wt%); Y/Nb (4.0 – 10.7); Th (0.1 – 6.8 ppm); Zr (18.6 – 576.9 ppm). The felsic group is composed by highly differentiated rocks which correspond to distinct levels of silica saturation and alkalinity. Peralkaline syenites usually present sodic (riebeckite) and sodic-calcic (aegirine-augite, ferrowinchite) inosilicates and reveal quite variable compositions: SiO2 (57.50 – 72.07 wt%); TiO2 (0.10 – 1.45 wt%); Th (1.7 – 67.0 ppm); Zr (133.0 – 4800.0 ppm). The alkaline granites show hedenbergite as the characteristic inosilicate, presenting relatively common compositions: SiO2 (61.85 – 78.06 wt%); TiO2 (0.21 – 0.58 wt%); Th (11.8 – 38.4 ppm); Zr (317.3 – 1234.6 ppm) [2]. Recent Sm-Nd isotopic results, on a total of 18 whole-rock samples (6 mafites and 12 felsites), allow new and more consistent interpretation concerning the petrogenesis of these plutonic rocks. Assuming an age of 490 Ma [3], the felsic rocks provide (0.6 < Nd490 < 4.3), similar to other contemporary (per)alkaline rocks of this region [4], reflecting magmatic extractions from time-integrated depleted mantle sources followed by variable and incomplete mixing (and/or AFC-type) processes with enriched, probably crustal sources. This alkaline/peralkaline magmatism is thought to represent the main regional record of the rifting event which presumably led to the opening of the Rheic Ocean. On the other hand, the mafic plutonic rocks of the Elvas massif cannot represent the magmatic precursors of these syenites and granites as they show completely distinct Nd isotopic ratios (3.7 < Nd490 < 1.2) indicating important contribution of long-term enriched (crustal) sources. Instead, considering the age and the Nd isotopic signature of other mafic plutonic unit emplaced nearby (the Campo Maior massif: ca. 370 Ma; 6.0 < Nd370 < 5.2) [5], and recalculating the isotopic ratios of the Elvas massif for the same age (4.3 < Nd370 < 1.6), it is plausible to consider that these plutons (Campo Maior and Elvas) can be coeval and representative of the Variscan magmatism in this region. In such hypothesis, the differences between these isotopic values could be explained, on a time-integrated basis, either by magmatic sources for the Elvas massif less enriched in LREE than the sources involved in the Campo Maior massif, or, if both plutonites share similar depleted mantle sources, by magmatic differentiation paths considerably affected by crustal contamination processes, which reached higher degrees in the Campo Maior massif.

Savani: A variable resolution whole-mantle model of anisotropic shear velocity variations based on multiple data sets

International audience

CONSTRAINTS ON THE LITHOSPHERIC STRUCTURE OF MID OCEAN RIDGES FROM OCEANIC CORE COMPLEX MORPHOLOGY

The Mid-oceanic ridge system is a feature unique to Earth. It is one of the fundamental components of plate tectonics and reflects interior processes of mantle convection within the Earth. The thermal structure beneath the mid-ocean ridges has been the subject of several modeling studies. It is expected that the elastic thickness of the lithosphere is larger near the transform faults that bound mid-ocean ridge segments. Oceanic core complexes (OCCs), which are generally thought to result from long-lived fault slip and elastic flexure, have a shape that is sensitive to elastic thickness. By modeling the shape of OCCs emplaced along a ridge segment, it is possible to constraint elastic thickness and therefore the thermal structure of the plate and how it varies along-axis. This thesis builds upon previous studies that utilize thin plate flexure to reproduce the shape of OCCs. I compare OCC shape to a suite of models in which elastic thickness, fault dip, fault heave, crustal thickness, and axial infill are systematically varied. Using a grid search, I constrain the parameters that best reproduce the bathymetry and/or the slope of ten candidate OCCs identified along the 12°—15°N segment of the Mid-Atlantic Ridge. The lithospheric elastic thicknesses that explains these OCCs is thinner than previous investigators suggested and the fault planes dip more shallowly in the subsurface, although at an angle compatible with Anderson’s theory of faulting. No relationships between model parameters and an oceanic core complexes location within a segment are identified with the exception that the OCCs located less than 20km from a transform fault have slightly larger elastic thickness than OCCs in the middle of the ridge segment.

Geologia, petrologia, geocronologia e mineralizações sulfetadas do complexo ézio, província mineral de carajás, Brasil

Dissertação (mestrado)—Universidade de Brasília, Instituto de Geociências, Pós-Graduação em Geologia, 2015.

The TOMO-ETNA experiment: an imaging active campaign at Mt. Etna volcano. Context, main objectives, working-plans and involved research projects

The TOMO-ETNA experiment was devised to image of the crust underlying the volcanic edifice and, possibly, its plumbing system by using passive and active refraction/reflection seismic methods. This experiment included activities both on-land and offshore with the main objective of obtaining a new high-resolution seismic tomography to improve the knowledge of the crustal structures existing beneath the Etna volcano and northeast Sicily up to Aeolian Islands. The TOMO ETNA experiment was divided in two phases. The first phase started on June 15, 2014 and finalized on July 24, 2014, with the withdrawal of two removable seismic networks (a Short Period Network and a Broadband network composed by 80 and 20 stations respectively) deployed at Etna volcano and surrounding areas. During this first phase the oceanographic research vessel “Sarmiento de Gamboa” and the hydro-oceanographic vessel “Galatea” performed the offshore activities, which includes the deployment of ocean bottom seismometers (OBS), air-gun shooting for Wide Angle Seismic refraction (WAS), Multi-Channel Seismic (MCS) reflection surveys, magnetic surveys and ROV (Remotely Operated Vehicle) dives. This phase finished with the recovery of the short period seismic network. In the second phase the Broadband seismic network remained operative until October 28, 2014, and the R/V “Aegaeo” performed additional MCS surveys during November 19-27, 2014. Overall, the information deriving from TOMO-ETNA experiment could provide the answer to many uncertainties that have arisen while exploiting the large amount of data provided by the cutting-edge monitoring systems of Etna volcano and seismogenic area of eastern Sicily.

Avaliação da contaminação por elementos metálicos dos sedimentos do Estuário Santos - São Vicente

Three approaches were applied to evaluate metal contamination in 41 sediment samples from the Santos - São Vicente Estuarine System: normalization to Al, statistical analysis and sediment quality guidelines (SQGs). The results showed increases in the concentrations of Zn, Ni, Pb, Cd, Cr and Hg, which seemed to be associated with human activities. The levels of Al, Fe and Co probably were associated with crustal material or natural weathering processes. About 45% of the samples presented concentrations exceeding TEL-ERL, levels occasionally associated with adverse biological effects. Four of these samples presented concentrations above PEL-ERM, levels frequently associated with adverse biological effects.

Contribuição de uma seção gravimétrica transversal ao estudo da estruturação litosférica na porção setentrional da Província Borborema, NE do Brasil

Este trabalho apresenta um modelo gravimétrico 2,5D gerado a partir de 315 novas estações gravimétricas levantadas ao longo de uma seção transversal NW-SE com 750 km de extensão na porção setentrional da Província Borborema, NE do Brasil. A modelagem gravimétrica foi aplicada separadamente nas componentes regional e residual do campo gravitacional. A modelagem 2,5D das anomalias regionais revelou que a profundidade da interface crosta-manto varia de 28 a 32 km, considerando uma densidade média de 2,8 g/cm3 para a crosta continental e de 3,3 g/cm3 para o manto litosférico. As anomalias residuais de alta frequência foram interpretadas a partir do contraste de densidade da crosta superior, com uma espessura não superior a 10 km, e uma ampla associação litológica, com densidades variando de 2,55 a 2,9 g/cm3. A configuração geotectônica litosférica atual da Província Borborema é claramente resultado da ruptura dos continentes Sul-Americano e Africano no Mesozoico, na qual boa parte dos vestígios das estruturas tectônicas de grandes profundidades formadas durante a orogênese Brasiliana/Pan-Africana foi mascarada pelo último episódio tectônico responsável pela fragmentação do Gondwana Ocidental. _______________________________________________________________________________________ ABSTRACT

A proveniência dos sedimentos e cinzas vulcânicas dos sedimentos permianos da bacia do Paraná : implicações para a história geológica do sul-sudoeste de Gondwana

Tese (doutorado)—Universidade de Brasília, Instituto de Geociências, Programa de Pós-Graduação em Geologia, 2016.