Arquitetura crustal da bacia de Almada no contexto das bacias da margem lesste da América do Sul

A Bacia de Almada, localizada no estado da Bahia, compartilha características similares com as outras bacias da margem leste do Brasil, quando é analisada segundo aspectos como os processos sedimentares e o regime de esforço dominante durante a sua formação. Observa-se uma diferença marcante em relação as outras bacias quando é analisada sob a ótica da composição da crosta transicional, uma vez que não se registra atividade vulcânica durante a fase rifte. A aquisição de um extenso levantamento sísmico 3D, com cabos de 6 km de comprimento e 9.2 segundos de tempo de registro (tempo sísmico duplo), resultaram em imagens sísmicas de boa qualidade das estruturas profundas do rifte. Adicionalmente, estudos de modelagem gravimétrica foram integrados com a análise sísmica para corroborar o modelo geológico. A Bacia de Almada é parte dos sistemas de rifte continentais, desenvolvidos durante o Berriasiano até o Aptiano, que antecederam a quebra do continente do Gondwana, evoluindo posteriormente para uma margem passiva divergente. O processo do rifteamento desenvolveu cinco sub-bacias de orientação NNE-SSO, desde posições terrestres até marinhas profundas, produzindo um arcabouço estrutural complexo. Os perfis da sísmica profunda mostram o afinamento progressivo da crosta continental até espessuras da ordem de 5 km, abaixo da sub-bacia mais oriental, com fatores de estiramento crustal próximo a 7 antes do desenvolvimento de crosta oceânica propriamente dita. As imagens sísmicas de boa qualidade permitem também o reconhecimento de sistemas de falhas lístricas que se iniciam na crosta superior, evoluem atravessando a crosta e conectando as sub-bacias para finalizar em um descolamento horizontal na crosta inferior estratificada. Adicionalmente, a bacia apresenta um perfil assimétrico, compatível com mecanismos de cisalhamento simples. As margens vulcânicas (VM) e não vulcânicas (NVM), são os extremos da análise composicional das margens divergentes continentais. Na Bacia de Almada não se reconhecem os elementos arquiteturais típicos das VM, tais como são as grandes províncias ígneas, caracterizadas por cunhas de refletores que mergulham em direção ao mar e por intenso vulcanismo pré- e sin-rifte nas bacias. Embora a margem divergente do Atlântico Sul seja interpretada tradicionalmente como vulcânica, o segmento do rifte ao sul do Estado da Bahia apresenta características não-vulcânicas, devido à ausência destes elementos arquiteturais e aos resultados obtidos nas perfurações geológicas que eventualmente alcançam a seqüência rifte e embasamento. Regionalmente a margem divergente sul-americana é majoritariamente vulcânica, embora a abundância e a influência do magmatísmo contemporâneo ao rifte seja muito variável. Ao longo da margem continental, desde a Bacia Austral no sul da Argentina, até a Bacia de Pernambuco no nordeste do Brasil, podem ser reconhecidos segmentos de caráter vulcânico forte, médio e não vulcânico. Nos exemplos clássicos de margens não vulcânicas, como a margem da Ibéria, a crosta transicional é altamente afinada podendo apresentar evidências de exumação de manto. Na Bacia de Almada, a crosta transicional apresenta importante estiramento embora não haja evidências concretas de exumação de manto. Os mecanismos responsáveis pela geração e intrusão dos grandes volumes de magma registrados nas margens divergentes são ainda sujeitos a intenso debate. Ao longo da margem divergente sul-americana há evidências da presença dos mecanismos genéticos de estiramento litosférico e impacto de plumas. Alternativamente estes dois mecanismos parecem ter tido um papel importante na evolução tectônica da margem sudeste e sul, diferenciando-as da margem continental onde foi implantada a Bacia de Almada.

Integração de métodos geofísicos na modelagem crustal da Bacia de Santos

A Bacia de Santos é uma bacia de margem passiva gerada durante o Neocomiano, associado com o evento de ruptura do Supercontinente de Gondwana. A espessura da seção sedimentar pode alcançar 15 km, incluindo a seqüência rifte, e foi depositada sobre uma crosta afinada cuja espessura original era de cerca de 35 km. Foi realizada uma modelagem crustal em 8 linhas, aproximadamente perpendiculares ao mergulho deposicional da bacia, objetivando o reconhecimento da geometria da fase rifte e a estimativa do estiramento crustal. Foram utilizados durante a modelagem quatro horizontes sísmicos em profundidade, mapas magnético e gravimétrico além de dados da profundidade do embasamento magnético. Foi reconhecido um pronunciado estiramento crustal em toda a bacia (fatores de estiramento crustal variando entre 1.2 e 3.1), evidenciando ampla acomodação da deformação. Os resultados da modelagem necessitaram, ainda, de espessas e contínuas camadas de rochas vulcânicas para ajuste do modelo. Estas cunhas de rochas vulcânicas, limitadas lateralmente por falhas normais, são feições comuns em bacias rifte. Os resultados da modelagem foram confrontados com dados da profundidade do embasamento magnético, obtidos anteriormente, e evidenciaram forte correlação, sobretudo nas áreas rasas da bacia.

Modelagem crustal da Bacia de Santos pela integração de métodos geofísicos

Pós-graduação em Geologia Regional - IGCE

O graben de Palestina: contribuição à estratigrafia e estrutura do estágio rifte na bacia do Araripe, nordeste do Brasil

The Palestina Graben is one of the NE-trending asymmetric grabens of the Araripe Basin. This basin rests on the precambrian terrains of the Transversal Zone, Borborema Province, immediately to the south of the Patos Lineament. It is part of the Interior Basins province of Northeastern Brazil, being related to the fragmentation of the Gondwana supercontinent and the opening of the South Atlantic ocean. The Palestina Graben trends NE-SW and presents an asymmetric geometry, controled by the NW extensional eocretaceous strain. The graben borders display distinct geometries. The SE border is a flexural margin, characterized by the non conformity of the eopaleozoic Mauriti Formation (the oldest unit of the basin) overlying the crystalline basement, but also affected by normal faults with small displacements. On the opposite, the NW border is continuous and rectilinear, being marked by normal faults with major displacements, that control the general tilting of the layers to the NW. In this sense, the Mauriti Formation is overlain by the Brejo Santo, Missão Velha (which also occurs in the Brejo Santo-Mauriti horst, to the NW of the fault border) and Abaiara formations, the latter restricted to the graben. The interpretation of available gravity data and a seismic line indicates that the main fault has a variable dip slip component, defining two deeper portions within the graben, in which the sedimentary column can reach thicknesses of up to 2 km. Regarding to the stratigraphy of Araripe Basin in the study area, the sedimentary package includes three distinct tectonosequences. The Paleozoic Syneclisis Tectonosequence is composed by the Mauriti Formation, deposited by a braided fluvial system. The Jurassic Tectonosequence, whose tectonic setting is still debatable (initial stage of the Neocomian rift, or a pre-rift syneclisis ?), is represented by the Brejo Santo Formation, originated in a distal floodplain related to ephemeral drainages. The Rift Tectonosequence, of neocomian age, includes the Missão Velha Formation, whose lower section is related to a braided to meandering fluvial system, outlining the Rift Initiation Tectonic Systems Tract. The upper section of the Missão Velha Formation is separated from the latter by a major unconformity. This interval was originated by a braided fluvial system, overlain by the Abaiara Formation, a deltaic system fed by a meandering fluvial system. Both sections correspond to the Rift Climax Tectonic Systems Tract. In the area, NE-trending normal to oblique faults are associated with NW transfer faults, while ENE to E-W faults display dominant strike slip kinematics. Both NE and E-W fault sets exhibit clear heritage from the basement structures (in particular, shear zones), which must have been reactivated during the eocretaceous rifting. Faults with EW trends display a dominant sinistral shear sense, commonly found along reactivated segments of the Patos Lineament and satellyte structures. Usually subordinate, dextral directional movements, occur in faults striking NNW to NE. Within this framework bearing to the Palestina Graben, classical models with orthogonal extension or pull-apart style deserve some caution in their application. The Palestina Graben is not limited, in its extremeties, by E-W transcurrent zones (as it should be in the case of the pull-apart geometry), suggesting a model close to the classic style of orthogonal opening. At the same time, others, adjacent depocenters (like the Abaiara-Jenipapeiro semi-graben) display a transtensional style. The control by the basement structures explains such differences

New insights into Crustal Contributions to Large-volume Rhyolite Generation in the Mid-Tertiary Sierra Madre Occidental Province, Mexico, revealed by U Pb Geochronology

Voluminous (≥3·9 × 105 km3), prolonged (∼18 Myr) explosive silicic volcanism makes the mid-Tertiary Sierra Madre Occidental province of Mexico one of the largest intact silicic volcanic provinces known. Previous models have proposed an assimilation–fractional crystallization origin for the rhyolites involving closed-system fractional crystallization from crustally contaminated andesitic parental magmas, with <20% crustal contributions. The lack of isotopic variation among the lower crustal xenoliths inferred to represent the crustal contaminants and coeval Sierra Madre Occidental rhyolite and basaltic andesite to andesite volcanic rocks has constrained interpretations for larger crustal contributions. Here, we use zircon age populations as probes to assess crustal involvement in Sierra Madre Occidental silicic magmatism. Laser ablation-inductively coupled plasma-mass spectrometry analyses of zircons from rhyolitic ignimbrites from the northeastern and southwestern sectors of the province yield U–Pb ages that show significant age discrepancies of 1–4 Myr compared with previously determined K/Ar and 40Ar/39Ar ages from the same ignimbrites; the age differences are greater than the errors attributable to analytical uncertainty. Zircon xenocrysts with new overgrowths in the Late Eocene to earliest Oligocene rhyolite ignimbrites from the northeastern sector provide direct evidence for some involvement of Proterozoic crustal materials, and, potentially more importantly, the derivation of zircon from Mesozoic and Eocene age, isotopically primitive, subduction-related igneous basement. The youngest rhyolitic ignimbrites from the southwestern sector show even stronger evidence for inheritance in the age spectra, but lack old inherited zircon (i.e. Eocene or older). Instead, these Early Miocene ignimbrites are dominated by antecrystic zircons, representing >33 to ∼100% of the dated population; most antecrysts range in age between ∼20 and 32 Ma. A sub-population of the antecrystic zircons is chemically distinct in terms of their high U (>1000 ppm to 1·3 wt %) and heavy REE contents; these are not present in the Oligocene ignimbrites in the northeastern sector of the Sierra Madre Occidental. The combination of antecryst zircon U–Pb ages and chemistry suggests that much of the zircon in the youngest rhyolites was derived by remelting of partially molten to solidified igneous rocks formed during preceding phases of Sierra Madre Occidental volcanism. Strong Zr undersaturation, and estimations for very rapid dissolution rates of entrained zircons, preclude coeval mafic magmas being parental to the rhyolite magmas by a process of lower crustal assimilation followed by closed-system crystal fractionation as interpreted in previous studies of the Sierra Madre Occidental rhyolites. Mafic magmas were more probably important in providing a long-lived heat and material flux into the crust, resulting in the remelting and recycling of older crust and newly formed igneous materials related to Sierra Madre Occidental magmatism.

Role of kilometer-scale weak circular heterogeneities on upper crustal deformation patterns: Evidence from scaled analogue modeling and the Sudbury Basin, Canada

The Sudbury Basin is a non-cylindrical fold basin occupying the central portion of the Sudbury Impact Structure. The impact structure lends itself excellently to explore the structural evolution of continental crust containing a circular region of long-term weakness. In a series of scaled analogue experiments various model crustal configurations were shortened horizontally at a constant rate. In mechanically weakened crust, model basins formed that mimic several first-order structural characteristics of the Sudbury Basin: (1) asymmetric, non-cylindrical folding of the Basin, (2) structures indicating concentric shortening around lateral basin termini and (3) the presence of a zone of strain concentration near the hinge zones of model basins. Geometrically and kinematically this zone corresponds to the South Range Shear Zone of the Sudbury Basin. According to our experiments, this shear zone is a direct mechanical consequence of basin formation, rather than the result of thrusting following folding. Overall, the models highlight the structurally anomalous character of the Sudbury Basin within the Paleoproterozoic Eastern Penokean Orogen. In particular, our models suggest that the Basin formed by pure shear thickening of crust, whereas transpressive deformation prevailed elsewhere in the orogen. The model basin is deformed by thickening and non-cylindrical synformal buckling, while conjugate transpressive shear zones propagated away from its lateral tips. This is consistent with pure shear deformation of a weak circular inclusion in a strong matrix. The models suggest that the Sudbury Basin formed as a consequence of long-term weakening of the upper crust by meteorite impact.

The causes of sinuous crustal-scale deformation patterns in hot orogens: Evidence from scaled analogue experiments and the southern Central Andes

The cause of upper-crustal segmentation into rhomb-shaped, shear zone-bound domains associated with contractional sedimentary basins in hot, wide orogens is not well understood. Here we use scaled multilayered analogue experiments to investigate the role of an orogen-parallel crustal-strength gradient on the formation of such structures. We show that the aspect ratio and size of domains, the sinuous character and abundance of transpressional shear zones vary with the integrated mechanical strength of crust. Upper-crustal deformation patterns and the degree of strain localization in the experiments are controlled by the ratio between the brittle and ductile strength in the model crust as well as gradients in tectonic and buoyancy forces. The experimental results match the first-order kinematic and structural characteristics of the southern Central Andes and provide insight on the dynamics of underlying deformation patterns in hot, wide orogens.

Heat-producing element enrichment in granitic rocks, the role of crustal composition and evolution

This study greatly enhanced our knowledge of the potential for geothermal energy development in Queensland as a viable clean energy source in the coming decades. Key outcomes of the project were understanding the first-order controls on the concentration of the heat-producing elements: uranium, thorium and potassium in granitic rocks, and constraining where rocks with the greatest heat-producing potential lie at depth in Queensland. Importantly, new temperature and heat flow maps for southwest Queensland were developed that will greatly assist future exploration efforts.

Low δ18O zircons from the Bruneau-­Jarbidge eruptive center: A key to crustal anatexis along the track of the Yellowstone hotspot

The Bruneau-Jarbidge eruptive center (BJEC) in the central Snake River Plain, Idaho, USA consists of the Cougar Point Tuff (CPT), a series of ten, high-temperature (900-1000°C) voluminous ignimbrites produced over the explosive phase of volcanism (12.8-10.5 Ma) and more than a dozen equally high-temperature rhyolite lava flows produced during the effusive phase (10.5-8 Ma). Spot analyses by ion microprobe of oxygen isotope ratios in 210 zircons demonstrate that all of the eruptive units of the BJEC are characterized by zircon δ¹⁸O values ≤ 2.5‰, thus documenting the largest low δ¹⁸O silicic volcanic province known on Earth (>10⁴ km³). There is no evidence for voluminous normal δ¹⁸O magmatism at the BJEC that precedes generation of low δ¹⁸O magmas as there is at other volcanic centers that generate low δ¹⁸O magmas such as Heise and Yellowstone. At these younger volcanic centers of the hotspot track, such low δ¹⁸O magmas represent ~45 % and ~20% respectively of total eruptive volumes. Zircons in all BJEC tuffs and lavas studied (23 units) document strong δ¹⁸O depletion (median CPT δ¹⁸OZrc = 1.0‰, post-CPT lavas = 1.5‰) with the third member of the CPT recording an excursion to minimum δ¹⁸O values (δ¹⁸OZrc= -1.8‰) in a supereruption > 2‰ lower than other voluminous low δ¹⁸O rhyolites known worldwide (δ¹⁸OWR ≤0.9 vs. 3.4‰). Subsequent units of the CPT and lavas record a progressive recovery in δ¹⁸OZrc to ~2.5‰ over a ~ 4 m.y. interval (12 to 8 Ma). We present detailed evidence of unit-to-unit systematic patterns in O isotopic zoning in zircons (i.e. direction and magnitude of Δcore-rim), spectrum of δ¹⁸O in individual units, and zircon inheritance patterns established by re-analysis of spots for U-Th-Pb isotopes by LA-ICPMS and SHRIMP. In conjunction with mineral thermometry and magma compositions, these patterns are difficult to reconcile with the well-established model for "cannibalistic" low δ¹⁸O magma genesis at Heise and Yellowstone. We present an alternative model for the central Snake River Plain using the modeling results of Leeman et al. (2008) for ¹⁸O depletion as a function of depth in a mid-upper crustal protolith that was hydrothermally altered by infiltrating meteoric waters prior to the onset of silicic magmatism. The model proposes that BJEC silicic magmas were generated in response to the propagation of a melting front, driven by the incremental growth of a vast underlying mafic sill complex, over a ~5 m.y. interval through a crustal volume in which a vertically asymmetric δ¹⁸OWR gradient had previously developed that was sharply inflected from ~ -1 to 10‰ at mid-upper crustal depths. Within the context of the model, data from BJEC zircons are consistent with incremental melting and mixing events in roof zones of magma reservoirs that accompany surfaceward advance of the coupled mafic-silicic magmatic system.

Archaean crustal evolution in eastern Finland : New geochronological and geochemical constraints from the Kuhmo terrain and the Nurmes belt

This thesis summarises the results of four original papers concerning U-Pb geochronology and geochemical evolution of Archaean rocks from the Kuhmo terrain and the Nurmes belt, eastern Finland. The study area belongs to a typical Archaean granite-greenstone terrain, composed of metavolcanic and metasedimentary rocks in generally N-S trending greenstone belts as well as a granitoid-gneiss complex with intervening gneissic and migmatised supracrustal and plutonic rocks. U-Pb data on migmatite mesosomes indicate that the crust surrounding the Tipasjärvi-Kuhmo-Suomussalmi greenstone belt is of varying age. The oldest protolith detected for a migmatite mesosome from the granitoid-gneiss complex is 2.94 Ga, whereas the other dated migmatites protoliths have ages of 2.84 2.79 Ga. The latter protoliths are syngenetic with the majority of volcanic rocks in the adjacent Tipasjärvi-Kuhmo-Suomussalmi greenstone belt. This suggests that the genesis of some of the volcanic rocks within the greenstone belt and surrounding migmatite protoliths could be linked. Metamorphic zircon overgrowths with ages of 2.84 2.81 Ga were also obtained. The non-migmatised plutonic rocks in the Kuhmo terrain and in the Nurmes belt record secular geochemical evolution, typical of Archaean cratons. The studied tonalitic rocks have ages of 2.83 2.75 Ga and they have geochemical characteristics similar to low-Al and high-Al TTD (tonalite-trondhjemite-dacite). The granodiorites, diorites, and gabbros with high Mg/Fe and LILE-enriched characteristics were mostly emplaced between 2.74 2.70 Ga and they exhibit geochemical characteristics typical of Archaean sanukitoid suites. The latest identified plutonic episode took place at 2.70 2.68 Ga, when compositionally heterogeneous leucocratic granitoid rocks, with a variable crustal component, were emplaced. U-Pb data on migmatite leucosomes suggest that leucosome generation may have been coeval with this latest plutonic event. On the basis of available U-Pb and Sm-Nd isotopic data it appears that the plutonic rocks of the Kuhmo terrain and the Nurmes belt do not contain any significant input from Palaeoarchaean sources. A characteristic feature of the Nurmes belt is the presence of migmatised paragneisses, locally preserving primary edimentary structures, with sporadic amphibolite intercalations. U-Pb studies on zircons indicate that the precursors of the Nurmes paragneisses were graywackes that were deposited between 2.71 Ga and 2.69 Ga and that they had a prominent 2.75 2.70 Ga source. Nd isotopic and whole-rock geochemical data for the intercalated amphibolites imply MORB sources. U-Pb data on zircons from the plutonic rocks and paragneisses reveal that metamorphic zircon growth took place at 2.72 2.63 Ga. This was the last tectonothermal event related to cratonisation of the Archaean crust of eastern Finland.

Physical Properties of Deep Drill Cores : Implications for Meteorite Impact Effects and Crustal Structures

Physical properties provide valuable information about the nature and behavior of rocks and minerals. The changes in rock physical properties generate petrophysical contrasts between various lithologies, for example, between shocked and unshocked rocks in meteorite impact structures or between various lithologies in the crust. These contrasts may cause distinct geophysical anomalies, which are often diagnostic to their primary cause (impact, tectonism, etc). This information is vital to understand the fundamental Earth processes, such as impact cratering and associated crustal deformations. However, most of the present day knowledge of changes in rock physical properties is limited due to a lack of petrophysical data of subsurface samples, especially for meteorite impact structures, since they are often buried under post-impact lithologies or eroded. In order to explore the uppermost crust, deep drillings are required. This dissertation is based on the deep drill core data from three impact structures: (i) the Bosumtwi impact structure (diameter 10.5 km, 1.07 Ma age; Ghana), (ii) the Chesapeake Bay impact structure (85 km, 35 Ma; Virginia, U.S.A.), and (iii) the Chicxulub impact structure (180 km, 65 Ma; Mexico). These drill cores have yielded all basic lithologies associated with impact craters such as post-impact lithologies, impact rocks including suevites and breccias, as well as fractured and unfractured target rocks. The fourth study case of this dissertation deals with the data of the Paleoproterozoic Outokumpu area (Finland), as a non-impact crustal case, where a deep drilling through an economically important ophiolite complex was carried out. The focus in all four cases was to combine results of basic petrophysical studies of relevant rocks of these crustal structures in order to identify and characterize various lithologies by their physical properties and, in this way, to provide new input data for geophysical modellings. Furthermore, the rock magnetic and paleomagnetic properties of three impact structures, combined with basic petrophysics, were used to acquire insight into the impact generated changes in rocks and their magnetic minerals, in order to better understand the influence of impact. The obtained petrophysical data outline the various lithologies and divide rocks into four domains. Based on target lithology the physical properties of the unshocked target rocks are controlled by mineral composition or fabric, particularly porosity in sedimentary rocks, while sediments result from diverse sedimentation and diagenesis processes. The impact rocks, such as breccias and suevites, strongly reflect the impact formation mechanism and are distinguishable from the other lithologies by their density, porosity and magnetic properties. The numerous shock features resulting from melting, brecciation and fracturing of the target rocks, can be seen in the changes of physical properties. These features include an increase in porosity and subsequent decrease in density in impact derived units, either an increase or a decrease in magnetic properties (depending on a specific case), as well as large heterogeneity in physical properties. In few cases a slight gradual downward decrease in porosity, as a shock-induced fracturing, was observed. Coupled with rock magnetic studies, the impact generated changes in magnetic fraction the shock-induced magnetic grain size reduction, hydrothermal- or melting-related magnetic mineral alteration, shock demagnetization and shock- or temperature-related remagnetization can be seen. The Outokumpu drill core shows varying velocities throughout the drill core depending on the microcracking and sample conditions. This is similar to observations by Kern et al., (2009), who also reported the velocity dependence on anisotropy. The physical properties are also used to explain the distinct crustal reflectors as observed in seismic reflection studies in the Outokumpu area. According to the seismic velocity data, the interfaces between the diopside-tremolite skarn layer and either serpentinite, mica schist or black schist are causing the strong seismic reflectivities.

Crustal Deformation and Seismic History Associated with the 2004 Indian Ocean Earthquake: A Perspective from the Andaman–Nicobar Islands

The Indian Ocean earthquake of 26 December 2004 led to significant ground deformation in the Andaman and Nicobar region, accounting for ~800 km of the rupture. Part of this article deals with coseismic changes along these islands, observable from coastal morphology, biological indicators, and Global Positioning System (GPS) data. Our studies indicate that the islands south of 10° N latitude coseismically subsided by 1–1.5 m, both on their eastern and western margins, whereas those to the north showed a mixed response. The western margin of the Middle Andaman emerged by >1 m, and the eastern margin submerged by the same amount. In the North Andaman, both western and eastern margins emerged by >1 m. We also assess the pattern of long-term deformation (uplift/subsidence) and attempt to reconstruct earthquake/tsunami history, with the available data. Geological evidence for past submergence includes dead mangrove vegetation dating to 740 ± 100 yr B.P., near Port Blair and peat layers at 2–4 m and 10–15 m depths observed in core samples from nearby locations. Preliminary paleoseismological/tsunami evidence from the Andaman and Nicobar region and from the east coast of India, suggest at least one predecessor for the 2004 earthquake 900–1000 years ago. The history of earthquakes, although incomplete at this stage, seems to imply that the 2004-type earthquakes are infrequent and follow variable intervals

Peraluminous sapphirine-cordierite pods in Mg-rich orthopyroxene granulite from southern India: Implications for lower crustal processes

Sapphirine-cordierite intergrowths occur as pods within garnet-absent, high-Mg orthopyroxene-granulite xenoliths in the Kambam valley, Madurai Block, southern India. Whereas the cores of the pods are composed of sapphirine (X-Mg = 0.871-0.897) - cordierite (X-Mg = 0.892-0.931) intergrowth along with rutile, zircon and monazite, the rims are characterized by cordierite, apatite, plagioclase, K-feldspar, quartz and minor calcite. The surrounding matrix comprises orthopyroxene (maximum Al2O3 4.1 wt.%, X-Mg 0.848-0.850), plagioclase, biotite and quartz, similar to the assemblage in the surrounding charnockites. Sapphirine in the Kambam rocks is characterized by high Al contents with an end-member composition in the range of 7:9:3 and 3:5:1. The occurrence of peraluminous sapphirine in association with cordierite and in the absence of phases such as sillimanite and garnet is distinct from ultrahigh-temperature assemblages in other localities within the Madurai Block. The peraluminous composition of the pods suggests that these domains could represent cryptic pathways through which aluminous melts migrated. The reaction of such peraluminous melts with Mg-rich orthopyroxene in the host granulite at temperatures of 1025 degrees C and pressures around 8 kbar as computed from phase equilibria modeling followed by an isobaric cooling is inferred to have generated the sapphirine-cordierite pods. The unusual high-Mg orthopyroxene granulite suggests interaction of supracrustal rocks with mafic magmas, which probably acted as the heat source for the partial melting of lower crust and UHT metamorphism.

High-temperature, high-pressure granulites (retrogressed eclogites) in the central region of the Lewisian, NW Scotland: Crustal-scale subduction in the Neoarchaean

Eclogites and associated high-pressure (HP) rocks in collisional and accretionary orogenic belts preserve a record of subduction and exhumation, and provide a key constraint on the tectonic evolution of the continents. Most eclogites that formed at high pressures but low temperatures at > 10-11 kbar and 450-650 degrees C can be interpreted as a result of subduction of cold oceanic lithosphere. A new class of high-temperature (HT) eclogites that formed above 900 degrees C and at 14 to 30 kbar occurs in the deep continental crust, but their geodynamic significance and processes of formation are poorly understood. Here we show that Neoarchaean mafic-ultramafic complexes in the central granulite facies region of the Lewisian in NW Scotland contain HP/HT garnet-bearing granulites (retrogressed eclogites), gabbros, Iherzolites, and websterites, and that the HP granulites have garnets that contain inclusions of omphacite. From thermodynamic modeling and compositional isopleths we calculate that peak eclogite-facies metamorphism took place at 24-22 kbar and 1060-1040 degrees C. The geochemical signature of one (G-21) of the samples shows a strong depletion of Eu indicating magma fractionation at a crustal level. The Sm-Nd isochron ages of HP phases record different cooling ages of ca. 2480 and 2330 Ma. We suggest that the layered mafic-ultramafic complexes, which may have formed in an oceanic environment, were subducted to eclogite depths, and exhumed as HP garnet-bearing orogenic peridotites. The layered complexes were engulfed by widespread orthogneisses of tonalite-trondhjemite-granodiorite (TTG) composition with granulite facies assemblages. We propose two possible tectonic models: (1) the fact that the relicts of eclogitic complexes are so widespread in the Scourian can be taken as evidence that a >90 km x 40 km-size slab of continental crust containing mafic-ultramafic complexes was subducted to at least 70 km depth in the late Archaean. During exhumation the gneiss protoliths were retrogressed to granulite facies assemblages, but the mafic-ultramafic rocks resisted retrogression. (2) The layered complexes of mafic and ultramafic rocks were subducted to eclogite-facies depths and during exhumation under crustal conditions they were intruded by the orthogneiss protoliths (TTG) that were metamorphosed in the granulite facies. Apart from poorly defined UHP metamorphic rocks in Norway, the retrogressed eclogites in the central granulite/retrogressed eclogite facies Lewisian region, NW Scotland have the highest crustal pressures so far reported for Archaean rocks, and demonstrate that lithospheric subduction was transporting crustal rocks to HP depths in the Neoarchaean. (C) 2012 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.