Thermal-elastic stresses and the criticality of the continental crust

Heating or cooling can lead to high stresses in rocks due to the different thermal-elastic properties of minerals. In the upper 4 km of the crust, such internal stresses might cause fracturing. Yet it is unclear if thermal elasticity contributes significantly to critical stresses and failure deeper in Earth's continental crust, where ductile creep causes stress relaxation. We combined a heating experiment conducted in a Synchrotron microtomograph (Advanced Photon Source, USA) with numerical simulations to calculate the grain-scale stress field in granite generated by slow burial. We find that deviatoric stresses >100 MPa can be stored during burial, with relaxation times from 100's to 1000's ka, even in the ductile crust. Hence, grain-scale thermal-elastic stresses may serve as nuclei for instabilities, thus rendering the continental crust close to criticality.

The origin of cyclical high- and low-Cr basalts with continental crust trace element signatures in the Kalkarindji Large Igneous Province

The 510 million year old Kalkarindji Large Igneous Province correlates in time with the first major extinction event after the Cambrian explosion of life. Large igneous provinces correlate with all major mass extinction events in the last 500 million years. The genetic link between large igneous provinces and mass extinction remains unclear. My work is a contribution towards understanding magmatic processes involved in the generation of Large Igneous Provinces. I concentrate on the origin of variation in Cr in magmas and have developed a model in which high temperature melts intrude into and assimilate large amounts of upper continental crust.

Low delta18O zircon grains in the Neoarchean Rum Jungle Complex, northern Australia: An indicator of emergent continental crust

The timing of widespread continental emergence is generally considered to have had a dramatic effect on the hydrological cycle, atmospheric conditions, and climate. New secondary ion mass spectrometry (SIMS) oxygen and laser-ablation–multicollector–inductively coupled plasma–mass spectrometry (LA-MC-ICP-MS) Lu-Hf isotopic results from dated zircon grains in the granitic Neoarchean Rum Jungle Complex provide a minimum time constraint on the emergence of continental crust above sea level for the North Australian craton. A 2535 ± 7 Ma monzogranite is characterized by magmatic zircon with slightly elevated δ18O (6.0‰–7.5‰ relative to Vienna standard mean ocean water [VSMOW]), consistent with some contribution to the magma from reworked supracrustal material. A supracrustal contribution to magma genesis is supported by the presence of metasedimentary rock enclaves, a large population of inherited zircon grains, and subchondritic zircon Hf (εHf = −6.6 to −4.1). A separate, distinct crustal source to the same magma is indicated by inherited zircon grains that are dominated by low δ18O values (2.5‰–4.8‰, n = 9 of 15) across a range of ages (3536–2598 Ma; εHf = −18.2 to +0.4). The low δ18O grains may be the product of one of two processes: (1) grain-scale diffusion of oxygen in zircon by exchange with a low δ18O magma or (2) several episodes of magmatic reworking of a Mesoarchean or older low δ18O source. Both scenarios require shallow crustal magmatism in emergent crust, to allow interaction with rocks altered by hydrothermal meteoric water in order to generate the low δ18O zircon. In the first scenario, assimilation of these altered rocks during Neoarchean magmatism generated low δ18O magma with which residual detrital zircons were able to exchange oxygen, while preserving their U-Pb systematics. In the second scenario, wholesale melting of the altered rocks occurred in several distinct events through the Mesoarchean, generating low δ18O magma from which zircon crystallized. Ultimately, in either scenario, the low δ18O zircons were entrained as inherited grains in a Neoarchean granite. The data suggest operation of a modern hydrological cycle by the Neoarchean and add to evidence for the increased emergence of continents by this time

Temperature-Time History of Subducted Continental-Crust, Mount Olympos Region, Greece

In the Mt. Olympos region of northeastern Greece, continental margin strata and basement rocks were subducted and metamorphosed under blueschist facies conditions, and thrust over carbonate platform strata during Alpine orogenesis. Subsequent exposure of the subducted basement rocks by normal faulting has allowed an integrated study of the timing of metamorphism, its relationship to deformation, and the thermal history of the subducted terrane. Alpine low-grade metamorphic assemblages occur at four structural levels. Three thrust sheets composed of Paleozoic granitic basement and Mesozoic metasedimentary cover were thrust over Mesozoic carbonate rocks and Eocene flysch; thrusting and metamorphism occurred first in the highest thrust sheets and progressed downward as units were imbricated from NE to SW. 40Ar/39Ar spectra from hornblende, white mica, and biotite samples indicate that the upper two units preserve evidence of four distinct thermal events: (1) 293–302 Ma crystallization of granites, with cooling from >550°C to <325°C by 284 Ma; (2) 98–100 Ma greenschist to blueschist-greenschist transition facies metamorphism (T∼350–500°C) and imbrication of continental thrust sheets; (3) 53–61 Ma blueschist facies metamorphism and deformation of the basement and continental margin units at T<350–400°C; (4) 36–40 Ma thrusting of blueschists over the carbonate platform, and metamorphism at T∼200–350°C. Only the Eocene and younger events affected the lower two structural packages. A fifth event, indicated by diffusive loss profiles in microcline spectra, reflects the beginning of uplift and cooling to T<100–150°C at 16–23 Ma, associated with normal faulting which continued until Quaternary time. Incomplete resetting of mica ages in all units constrains the temperature of metamorphism during continental subduction to T≤350°C, the closure temperature for Ar in muscovite. The diffusive loss profiles in micas and K-feldspars enable us to “see through” the younger events to older events in the high-T parts of the release spectra. Micas grown during earlier metamorphic events lost relatively small amounts of Ar during subsequent high pressure-low temperature metamorphism. Release spectra from phengites grown during Eocene metamorphism and deformation record the ages of the Ar-loss events. Alpine deformation in northern Greece occurred over a long time span (∼90 Ma), and involved subduction and episodic imbrication of continental basement before, during, and after the collision of the Apulian and Eurasian plates. Syn-subduction uplift and cooling probably combined with intermittently higher cooling rates during extensional events to preserve the blueschist facies mineral assemblages as they were exhumed from depths of >20 km. Extension in the Olympos region was synchronous with extension in the Mesohellenic trough and the Aegean back-arc, and concurrent with westward-progressing shortening in the external Hellenides.

Contrasting PT paths within the Barchi-Kol UHP terrain (Kokchetav Complex): Implications for subduction and exhumation of continental crust

The Barchi-Kol terrain is a classic locality of ultrahigh-pressure (UHP) metamorphism within the Kokchetav metamorphic belt. We provide a detailed and systematic characterization of four metasedimentary samples using dominant mineral assemblages, mineral inclusions in zircon and monazite, garnet zonation with respect to major and trace elements, and Zr-in-rutile and Ti-in-zircon temperatures. A typical diamond-bearing gneiss records peak conditions of 49 ± 4 kbar and 950–1000 °C. Near isothermal decompression of this rock resulted in the breakdown of phengite associated with a pervasive recrystallization of the rock. The same terrain also contains mica schists that experienced peak conditions close to those of the diamond-bearing rocks, but they were exhumed along a cooler path where phengite remained stable. In these rocks, major and trace element zoning in garnet has been completely equilibrated. A layered gneiss was metamorphosed at UHP conditions in the coesite field, but did not reach diamond-facies conditions (peak conditions: 30 kbar and 800–900 °C). In this sample, garnet records retrograde zonation in major elements and also retains prograde zoning in trace elements. A garnet-kyanite-micaschist that reached significantly lower pressures (24 ± 2 kbar, 710 ± 20 °C) contains garnet with major and trace element zoning. The diverse garnet zoning in samples that experienced different metamorphic conditions allows to establish that diffusional equilibration of rare earth element in garnet likely occurs at ~900–950 °C. Different metamorphic conditions in the four investigated samples are also documented in zircon trace element zonation and mineral inclusions in zircon and monazite. U-Pb geochronology of metamorphic zircon and monazite domains demonstrates that prograde (528–521 Ma), peak (528–522 Ma), and peak to retrograde metamorphism (503–532 Ma) occurred over a relatively short time interval that is indistinguishable from metamorphism of other UHP rocks within the Kokchetav metamorphic belt. Therefore, the assembly of rocks with contrasting P-T trajectories must have occurred in a single subduction-exhumation cycle, providing a snapshot of the thermal structure of a subducted continental margin prior to collision. The rocks were initially buried along a low geothermal gradient. At 20–25 kbar they underwent near isobaric heating of 200 °C, which was followed by continued burial along a low geothermal gradient. Such a step-wise geotherm is in good agreement with predictions from subduction zone thermal models.

Buoyancy-driven, rapid exhumation of ultrahigh-pressure metamorphosed continental crust

Preservation of ultrahigh-pressure (UHP) minerals formed at depths of 90–125 km require unusual conditions. Our subduction model involves underflow of a salient (250 ± 150 km wide, 90–125 km long) of continental crust embedded in cold, largely oceanic crust-capped lithosphere; loss of leading portions of the high-density oceanic lithosphere by slab break-off, as increasing volumes of microcontinental material enter the subduction zone; buoyancy-driven return toward midcrustal levels of a thin (2–15 km thick), low-density slice; finally, uplift, backfolding, normal faulting, and exposure of the UHP terrane. Sustained over ≈20 million years, rapid (≈5 mm/year) exhumation of the thin-aspect ratio UHP sialic sheet caught between cooler hanging-wall plate and refrigerating, downgoing lithosphere allows withdrawal of heat along both its upper and lower surfaces. The intracratonal position of most UHP complexes reflects consumption of an intervening ocean basin and introduction of a sialic promontory into the subduction zone. UHP metamorphic terranes consist chiefly of transformed, yet relatively low-density continental crust compared with displaced mantle material—otherwise such complexes could not return to shallow depths. Relatively rare metabasaltic, metagabbroic, and metacherty lithologies retain traces of phases characteristic of UHP conditions because they are massive, virtually impervious to fluids, and nearly anhydrous. In contrast, H2O-rich quartzofeldspathic, gneissose/schistose, more permeable metasedimentary and metagranitic units have backreacted thoroughly, so coesite and other UHP silicates are exceedingly rare. Because of the initial presence of biogenic carbon, and its especially sluggish transformation rate, UHP paragneisses contain the most abundantly preserved crustal diamonds.

A refined solution to Earth's hidden niobium: implications for evolution of continental crust and mode of core formation

New high-precision niobium (Nb) and tantalum (Ta) concentration data are presented for early Archaean metabasalts, metabasaltic komatiites and their erosion products (mafic metapelites) from SW Greenland and the Acasta gneiss complex, Canada. Individual datasets consistently show sub-chondritic Nb/Ta ratios averaging 15.1+/-11.6. This finding is discussed with regard to two competing models for the solution of the Nb-deficit that characterises the accessible Earth. Firstly, we test whether Nb could have sequestered into the core due to its slightly siderophile (or chalcophile) character under very reducing conditions, as recently proposed from experimental evidence. We demonstrate that troilite inclusions of the Canyon Diablo iron meteorite have Nb and V concentrations in excess of typical chondrites but that the metal phase of the Grant, Toluca and Canyon Diablo iron meteorites do not have significant concentrations of these lithophile elements. We find that if the entire accessible Earth Nb-deficit were explained by Nb in the core, only ca. 17% of the mantle could be depleted and that by 3.7 Ga, continental crust would have already achieved ca. 50% of its present mass. Nb/Ta systematics of late Archaean metabasalts compiled from the literature would further require that by 2.5 Ga, 90% of the present mass of continental crust was already in existence. As an alternative to this explanation, we propose that the average Nb/Ta ratio (15.1+/-11.6) of Earth's oldest mafic rocks is a valid approximation for bulk silicate Earth. This would require that ca. 13% of the terrestrial Nb resided in the Ta-free core. Since the partitioning of Nb between silicate and metal melts depends largely on oxygen fugacity and pressure, this finding could mean that metal/silicate segregation did not occur at the base of a deep magma ocean or that the early mantle was slightly less reducing than generally assumed. A bulk silicate Earth Nb/Ta ratio of 15.1 allows for depletion of up to 40% of the total mantle. This could indicate that in addition to the upper mantle, a portion of the lower mantle is depleted also, or if only the upper mantle were depleted, an additional hidden high Nb/Ta reservoir must exist. Comparison of Nb/Ta systematics between early and late Archaean metabasalts supports the latter idea and indicates deeply subducted high Nb/Ta eclogite slabs could reside in the mantle transition zone or the lower mantle. Accumulation of such slabs appears to have commenced between 2.5 and 2.0 Ga. Regardless of these complexities of terrestrial Nb/Ta systematics, it is shown that the depleted mantle Nb/Th ratio is a very robust proxy for the amount of extracted continental crust, because the temporal evolution of this ratio is dominated by Th-loss to the continents and not Nb-retention in the mantle. We present a new parameterisation of the continental crust volume versus age curve that specifically explores the possibility of lithophile element loss to the core and storage of eclogite slabs in the transition zone. (C) 2003 Elsevier Science B.V. All rights reserved.

A new estimate for the composition of weathered young upper continental crust from alluvial sediments, Queensland, Australia

We present new major element, trace element and Nd-isotope data for 30 alluvial sediments collected from 25 rivers in Queensland, E Australia. Samples were chosen to represent drainage from the region's most important lithologies, including Tertiary intraplate volcanic rocks, a Cretaceous igneous province (and sedimentary rocks derived thereof) as well as Proterozoic blocks. In most chemical and isotopic aspects, the alluvial sediments represent binary or ternary mixing relationships, with absolute abundances implied to reflect the proportion of lithologies in the catchments. When averaged, the studied sediments differ from other proxies of upper continental crust (UCC) mainly in their relative middle rare earth element enrichment (including an elevated Sm/Nd ratio), higher relative Eu abundance and higher Nb/Ta ratio. These features are inherited from eroded Tertiary intraplate basalts, which commonly form topographic highs in the studied region. Despite the high degree of weathering strong to excellent coherence between similarly incompatible elements is found for all samples. From this coherence, we suggest revisions of the following upper crustal element ratios: Y/Ho = 26.2, Yb/Tm = 6.37, Th/W = 7.14, Th/Tl = 24 and Zr/Hf = 36.9. Lithium, Rb, Cs and Be contents do not seem depleted relative to UCC, which may reflect paucity of K-feldspar in the eroded catchments. Nickel, Cr, Pb, Cu and Zn concentrations are elevated in polluted rivers surrounding the state capital. River sediments in the Proterozoic Georgetown Inlier are elevated in Pb, Cu and Zn but this could be a natural phenomenon reflecting abundant sulphide mineralisation in the area. Except for relative Sr concentrations, which broadly anticorrelate with mean annual rainfall in catchments, there is no obvious relationship between the extent of weathering and climate types, which range from and to tropical. The most likely explanation for this observation is that the weathering profiles in many catchments are several Myr old, established during the much wetter Miocene period. The studied sediment compositions (excluding those from the Proterozoic catchments) are used to propose a new trace element normalisation termed MUQ (MUd from Queensland), which serves as an alternative to UCC proxies derived from sedimentary rocks. Copyright (C) 2005 Elsevier Ltd

Thinning of the Goban Spur continental margin and formation of early oceanic crust: Constraints from forward modelling and inversion of marine magnetic anomalies

The deep seismic reflection profile Western Approaches Margin (WAM) cuts across the Goban Spur continental margin, located southwest of Ireland. This non-volcanic margin is characterized by a few tilted blocks parallel to the margin. A volcanic sill has been emplaced on the westernmost tilted block. The shape of the eastern part of this sill is known from seismic data, but neither seismic nor gravity data allow a precise determination of the extent and shape of the volcanic body at depth. Forward modelling and inversion of magnetic data constrain the shape of this volcanic sill and the location of the ocean-continent transition. The volcanic body thickens towards the ocean, and seems to be in direct contact with the oceanic crust. In the contact zone, the volcanic body and the oceanic magnetic layer display approximately the same thickness. The oceanic magnetic layer is anomalously thick immediately west of the volcanic body, and gradually thins to reach more typical values 40 km further to the west. The volcanic sill would therefore represent the very first formation of oceanic crust, just before or at the continental break-up. The ocean-continent transition is limited to a zone 15 km wide. The continental magnetic layer seems to thin gradually oceanwards, as does the continental crust, but no simple relation is observed between their respective thinnings.