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

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

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.

Helium and neon isotopes in oceanic crust of ODP Hole 118-735B

In an attempt to determine the helium and neon isotopic composition of the lower oceanic crust, we report new noble gas measurements on 11 million year old gabbros from Ocean Drilling Program site 735B in the Indian Ocean. The nine whole rock samples analyzed came from 20 to 500 m depth below the seafloor. Helium contents vary from 3.3*10**-10 to 2.5*10**-7 ccSTP/g by crushing and from 5.4*10**-8 to 2.4*10**-7 ccSTP/g by melting. 3He/4He ratios vary between 2.2 and 8.6 Ra by crushing and between 2.9 and 8.2 by melting. The highest R/Ra ratios are similar to the mean mid-ocean ridge basalt (MORB) ratio of 8+/-1. The lower values are attributed to radiogenic helium from in situ alüha-particle production during uranium and thorium decay. Neon isotopic ratios are similar to atmospheric ratios, reflecting a significant seawater circulation in the upper 500 m of exposed crust at this site. MORB-like neon, with elevated 20Ne/22Ne and 21Ne/22Ne ratios, was found in some high temperature steps of heating experiments, but with very small anomalies compared to air. These first results from the lower oceanic crust indicate that subducted lower oceanic crust has an atmospheric 20Ne/22Ne ratio. Most of this neon must be removed during the subduction process, if the ocean crust is to be recirculated in the upper mantle, otherwise this atmospheric neon will overwhelm the upper mantle neon budget. Similarly, the high (U+Th)/3He ratio of these crustal gabbros will generate very radiogenic 4He/3He ratios on a 100 Ma time scale, so lower oceanic crust cannot be recycled into either MORB or oceanic island basalt without some form of processing.

Element contents, sulfur and oxygen isotopes, and molecular biomarkers of phosphatic crust samples from the Chimbote Platform of the shelf off Peru

Authigenic phosphatic laminites enclosed in phosphorite crusts from the shelf off Peru (10°01' S and 10°24' S) consist of carbonate fluorapatite layers, which contain abundant sulfide minerals including pyrite (FeS2) and sphalerite (ZnS). Low d34Spyrite values (average -28.8 per mill) agree with bacterial sulfate reduction and subsequent pyrite formation. Stable sulfur isotopic compositions of sulfate bound in carbonate fluorapatite are lower than that of sulfate from ambient sea water, suggesting bacterial reoxidation of sulfide by sulfide-oxidizing bacteria. The release of phosphorus and subsequent formation of the autochthonous phosphatic laminites are apparently caused by the activity of sulfate-reducing bacteria and associated sulfide-oxidizing bacteria. Following an extraction-phosphorite dissolution-extraction procedure, molecular fossils of sulfate-reducing bacteria (mono-O-alkyl glycerol ethers, di-O-alkyl glycerol ethers, as well as the short-chain branched fatty acids i/ai-C15:0, i/ai-C17:0 and 10MeC16:0) are found to be among the most abundant compounds. The fact that these molecular fossils of sulfate-reducing bacteria are distinctly more abundant after dissolution of the phosphatic laminite reveals that the lipids are tightly bound to the mineral lattice of carbonate fluorapatite. Moreover, compared with the autochthonous laminite, molecular fossils of sulfate-reducing bacteria are: (1) significantly less abundant and (2) not as tightly bound to the mineral lattice in the other, allochthonous facies of the Peruvian crusts consisting of phosphatic coated grains. These observations confirm the importance of sulfate-reducing bacteria in the formation of the phosphatic laminite. Model calculations highlight that organic matter degradation by sulfate-reducing bacteria has the potential to liberate sufficient phosphorus for phosphogenesis.

Vein mineral ages of oceanic crust from the Atlantic Ocean

Vein smectites with large Rb/Sr enrichments from extensively altered basaltic oceanic crust in Deep Sea Drilling Project hole 417A in the western Atlantic define a highly constrained Rb/Sr isochron age of 108 +/- 3 m.y. This age is identical to a less well constrained age of 108 +/- 17 m.y. for vein smectites with lower Rb/Sr enrichments from adjacent hole 418A and to the 108 m.y. age of crust formation derived by paleontological and magnetic anomaly correlation. Reasonable agreement exists between the 87Sr/86Sr ratio of vein calcites from both sites and the seawater 87Sr/86Sr ratio at the time. Pervasive low-temperature alteration in the contrasting environments of sites 417 and 418 appears to be coeval and essentially coincident with basement formation. Alteration may be used to advantage in determining ages of old oceanic crust.

Alkali contents, Sr isotope and U-Th-Pb data for ocean crust vein minerals from DSDP Holes 27-261, 61-462A, and 72-516F

Cation exchange experiments (ammonium acetate and cation resin) on celadonite-smectite vein minerals from three DSDP holes demonstrate selective removal of common Sr relative to Rb and radiogenic Sr. This technique increases the Rb/Sr ratio by factors of 2.3 to 22 without significantly altering the age of the minerals, allowing easier and more precise dating of such vein minerals. The ages determined by this technique (Site 261 - 121.4+/-1.6 m.y.; Site 462A - 105.1+/-2.8 m.y.; Site 516F - 69.9+/-2.4 m.y.) are 34, 54 and 18 m.y. younger, respectively, than the age of crust formation at the site; in the case of site 462A, the young age is clearly related to off-ridge emplacement of a massive sill/flow complex. At the other sites, either the hydrothermal circulation systems persisted longer than for normal crust (10-15 m.y.), or were reactivated by off-ridge igneous activity. Celadonites show U and Pb contents and Pb isotopic compositions little changed from their basalt precursors, while Th contents are significantly lower. Celadonites thus have unusually high alkali/U,Th ratios and low Th/U ratios. If this celadonite alteration signature is significantly imprinted on oceanic crust as a whole, it will lead to very distinctive Pb isotope signatures for any hot spot magmas which contain a component of aged subducted recycled oceanic crust. Initial Sr isotope ratios of ocean crust vein minerals (smectite, celadonite, zeolite, calcite) are intermediate between primary basalt values and contemporary sea water values and indicate formation under seawaterdominated systems with effective water/rock ratios of 20-200.

The petrology of the lower series volcanics of ODP Hole 104-642E

Between 1086.6 and 1229.4 m below seafloor at Site 642 on the Outer Vøring Plateau, a series of intermediate volcanic extrusive flow units and volcaniclastic sediments was sampled. A mixed sequence of dacitic subaerial flows, andesitic basalts, intermediate volcaniclastics, subordinate mid-ocean ridge basalt, (MORB) lithologies, and intrusives was recovered, in sharp contrast to the more uniform tholeiitic T-type MORB units of the overlying upper series. This lower series of volcanics is composed of three chemically distinct groups, (B, A2, A1), rather than the two previously identified. Flows of the dacitic group (B) have trace-element and initial Sr isotope signatures which indicate that their source magma derived from the partial melting of a component of continental material in a magma chamber at a relatively high level in the crust. The relative proportions of crustal components in this complex melt are not known precisely. The most basic group (A2) probably represents a mixture of this material with MORB-type tholeiitic melt. A third group (A1), of which there was only one representative flow recovered, is chemically intermediate between the two groups above, and may suggest a repetition of, or a transition phase in, the mixing processes.

Stable oxygen isotope ratios in zircons from oceanic crust

Lower ocean crust is primarily gabbroic, although 1-2% felsic igneous rocks that are referred to collectively as plagiogranites occur locally. Recent experimental evidence suggests that plagiogranite magmas can form by hydrous partial melting of gabbro triggered by seawater-derived fluids, and thus they may indicate early, high-temperature hydrothermal fluid circulation. To explore seawater-rock interaction prior to and during the genesis of plagiogranite and other late-stage magmas, oxygen-isotope ratios preserved in igneous zircon have been measured by ion microprobe. A total of 197 zircons from 43 plagiogranite, evolved gabbro, and hydrothermally altered fault rock samples have been analyzed. Samples originate primarily from drill core acquired during Ocean Drilling Program and Integrated Ocean Drilling Program operations near the Mid-Atlantic and Southwest Indian Ridges. With the exception of rare, distinctively luminescent rims, all zircons from ocean crust record remarkably uniform d18O with an average value of 5.2 ± 0.5 per mil (2SD). The average d18O(Zrc) would be in magmatic equilibrium with unaltered MORB [d18O(WR) ~5.6-5.7 per mil], and is consistent with the previously determined value for equilibrium with the mantle. The narrow range of measured d18O values is predicted for zircon crystallization from variable parent melt compositions and temperatures in a closed system, and provides no indication of any interactions between altered rocks or seawater and the evolved parent melts. If plagiogranite forms by hydrous partial melting, the uniform mantle-like d18O(Zrc) requires melting and zircon crystallization prior to significant amounts of water-rock interactions that alter the protolith d18O. Zircons from ocean crust have been proposed as a tectonic analog for >3.9 Ga detrital zircons from the earliest (Hadean) Earth by multiple workers. However, zircons from ocean crust are readily distinguished geochemically from zircons formed in continental crustal environments. Many of the >3.9 Ga zircons have mildly elevated d18O (6.0-7.5 per mil), but such values have not been identified in any zircons from the large sample suite examined here. The difference in d18O, in combination with newly acquired lithium concentrations and published trace element data, clearly shows that the >3.9 Ga detrital zircons did not originate by processes analogous to those in modern mid-ocean ridge settings.

Carbon chemistry of oceanic crust

Carbonate mineral precipitation in the upper oceanic crust during low-temperature, off-axis, hydrothermal circulation is investigated using new estimates of the bulk CO2 content of seven DSDP/ODP drill cores. In combination with previously published data these new data show: (i) the CO2 content of the upper ~ 300 m of the crust is substantially higher in Cretaceous than in Cenozoic crust and (ii) for any age of crust, there is substantially more CO2 in Atlantic (slow-spreading) than Pacific (intermediate- to fast-spreading) crust. Modelling the Sr-isotopic composition of the carbonates suggests that > 80% of carbonate mineral formation occurs within < 20 Myr of crust formation. This means that the higher CO2 content of Cretaceous crust reflects a secular change in the rate of CO2 uptake by the crust. Oxygen isotope derived estimates of carbonate mineral precipitation temperatures show that the average and minimum temperature of carbonate precipitation was ~10 °C higher temperatures in the Cretaceous than in the Cenozoic. This difference is consistent with previous estimates of secular change in bottom seawater temperature. Higher fluid temperature within the crust will have increased reaction rates potentially liberating more basaltic Ca and hence enhancing carbonate mineral precipitation. Additionally, if crustal fluid pH is controlled by fluid-rock reaction, the higher Ca content of the Cretaceous ocean will also have enhanced carbonate mineral precipitation. New estimates of the rate of CO2 uptake by the upper ocean crust during the Cenozoic are much lower than previous estimates.

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.