Metamorphic growth and recrystallization of zircons in extremely O-depleted rocks during eclogite-facies metamorphism: Evidence from U-Pb ages, trace elements, and O-Hf isotopes

A combined in situ SIMS and LA-(MC)-ICPMS study of U-Pb ages, trace elements, O and Lu-Hf isotopes was conducted for zircon from eclogite-facies metamorphic rocks in the Sulu orogen. The two microbeam techniques sampled various depths of zircon domains, revealing different element and isotope relationships between residual magmatic cores and new metamorphic rims and thus the geochemical architecture of metamorphic zircons which otherwise cannot be recognized by the single microbeam technique. This enables discrimination of metamorphic growth from different subtypes of metamorphic recrystallization. Magmatic cores with U-Pb ages of 769 ± 9 Ma have positive δ¹⁸O values of 0.1-10.1‰, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high REE contents, and steep MREE-HREE patterns with negative Eu anomalies. They are interpreted as crystallizing from positive δ¹⁸O magmas during protolith emplacement. In contrast, newly grown domains have concordant U-Pb ages of 204 ± 4 to 252 ± 7 Ma and show negative δ¹⁸O values of −10.0‰ to −2.2‰, low Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, low REE contents, and flat HREE patterns with weak to no Eu anomalies. They are interpreted as growing from negative δ¹⁸O fluids that were produced by metamorphic dehydration of high-T glacial-hydrothermally altered rocks during continental subduction-zone metamorphism. Differences in δ¹⁸O between different domains within single grains vary from 0.8‰ to 12.5‰, suggesting different degrees of O isotope exchange between the positive δ¹⁸O magmatic core and the negative δ¹⁸O metamorphic fluid during the metamorphism. The magmatic zircons underwent three subtypes of metamorphic recrystallization, depending on their accessibility to negative δ¹⁸O fluids. The zircons recrystallized in solid-state maintained positive δ¹⁸O values, and REE and Lu-Hf isotopes of protolith zircon, but their U-Pb ages are lowered. The zircons recrystallized through dissolution exhibit negative δ¹⁸O values similar to the metamorphic growths, almost completely reset U-Pb ages, and partially reset REE systems. The zircons recrystallized through replacement show variably negative δ¹⁸O values, and partially reset REE, and U-Pb and Lu-Hf isotopic systems. Therefore, this study places robust constraints on the origin of metamorphic zircons in eclogite-facies rocks and provides a methodological framework for linking the different types of metamorphic zircons to petrological processes during continental collision.     

Zircon U-Pb and oxygen isotope evidence for a large-scale O depletion event in igneous rocks during the Neoproterozoic

Oxygen isotope analysis and U-Pb dating were carried out on zircons from granite, granitic gneiss and eclogite in the Dabie-Sulu orogenic belt of east-central China. The results show a wide variation in zircon δ¹⁸O values from −10.9 to 8.5‰. Most of the values are lower than normal mantle zircon and almost half have prominently negative values. The low δ¹⁸O zircons have protolith ages of 700 to 800 Ma and metamorphic ages of 205 to 250 Ma, respectively. Fluid availability within the metaigneous rocks dictates the extent of metamorphic recrystallization and overgrowth. The igneous zircons have preserved their magmatic zoning and middle Neoproterozoic U-Pb age during the Triassic metamorphism, indicating low fluid availability. Widespread low δ¹⁸O values are identified in the magmatic zircons of middle Neoproterozoic age, within an outcrop area of over 20,000 km² along the northern edge of the South China Block. The low δ¹⁸O zircons record the presence of large volumes of low δ¹⁸O igneous rocks that were derived from remelting of meteoric-hydrothermally altered rocks at some time between 700 and 800 Ma.The U-Pb ages for metaigneous protoliths and granites are correlated not only with the timing of rifting accompanying the breakup of Rodinia, but also with contemporary glacial deposits in the South China Block at paleolatitudes of 30 to 40°N. Melting of glacial ice and snow is suggested, instead of the direct involvement of meteoric water, to produce the low δ¹⁸O fluid with oxygen isotopic signatures like the cold-climate meteoric water. The rifting created conditions favorable to anatexis of meteoric-hydrothermally altered rocks. Glaciated regions supplied copious water for the water-rock interaction during magma emplacement along rifting zones. Both rifting and glaciation favored the generation of the low δ¹⁸O magmas in the region. The low δ¹⁸O zircons are thus interpreted to have crystallized from the low δ¹⁸O magmas of middle Neoproterozoic age, and onset of the generation of low δ¹⁸O magmas in the northern margin of the South China Block is estimated to occur at about 758 ± 15 Ma. The large-scale remelting of hydrothermally altered crust not only results from repetitive emplacement of mafic magmas along the rifting zones with protracted episodes of water-rock interaction, but also involves rift systems that rapidly introduce large volumes of fluid through confined pathways and traps in a short space of time in response to tectonic triggers. Occurrence of the large-scale ¹⁸O depletion during the middle Neoproterozoic may be a manifestation of the cold paleoclimate related to the snowball Earth event.     

U-Pb, Hf and O isotope evidence for two episodes of fluid-assisted zircon growth in marble-hosted eclogites from the Dabie orogen

A combined study of internal structure, U-Pb age, and Hf and O isotopes was carried out for metamorphic zircons from ultrahigh-pressure eclogite boudins enclosed in marbles from the Dabie orogen in China. CL imaging identifies two types of zircon that are metamorphically new growth and recrystallized domain, respectively. The metamorphic zircons have low Th and U contents with low Th/U ratios, yielding two groups of ²⁰⁶Pb/²³⁸U age at 245 ± 3 to 240 ± 2 Ma and 226 ± 4 to 223 ± 2 Ma, respectively. Anomalously high δ¹⁸O values were obtained for refractory minerals, with 9.9 to 21.4‰ for garnet and 16.9‰ for zircon. This indicates that eclogite protolith is sedimentary rocks capable of liberating aqueous fluid for zircon growth during continental subduction-zone metamorphism. Most of the zircons are characterized by very low ¹⁷⁶Lu/¹⁷⁷Hf ratios of 0.000001-0.000028, indicating their growth in association with garnet recrystallization. A few of them falling within the older age group have comparatively high ¹⁷⁶Lu/¹⁷⁷Hf ratios of 0.000192-0.000383, suggesting their growth prior to the formation of garnet in the late stage of subduction. The variations in the Lu/Hf ratios for zircons can thus be used to correlate with garnet growth during eclogite-facies metamorphism. In either case, the zircons have variable ε_Hf (t) values for individual samples, suggesting that their protolith is heterogeneous in Hf isotope composition with localized fluid availability in the bulk processes of orogenic cycle. Nevertheless, a positive correlation exists between ²⁰⁶Pb/²³⁸U ages and Lu-Hf isotope ratios for the metamorphically recrystallized zircons, suggesting that eclogite-facies metamorphism in the presence of fluid has the identical effect on zircon Lu-Hf and U-Th-Pb isotopic systems. We conclude that the zircons of the older group grew in the presence of fluid during the subduction prior to the onset of peak ultrahigh-pressure metamorphism, whereas the younger zircons grew in the presence of fluid released during the initial exhumation toward high-pressure eclogite-facies regime.     

威海地区超高压变质花岗片麻岩锆石U-Pb定年和氧同位素研究

在苏鲁超高压变质带东北端山东威海地区皂埠镇发现锆石δ18O值低至-7.8‰左右的花岗片麻岩, 与前人在苏鲁超高压变质带西南端江苏东海青龙山地区发现的锆石δ18O值为-7‰~-9‰左右的花岗片麻岩一致.对这些低δ18O值花岗片麻岩进行了锆石SHRIMP法UPb定年和系统的激光氟化法矿物氧同位素分析, 结果对低δ18O值锆石成因和花岗片麻岩的原岩性质提供了制约.研究得到: (1) 这些低δ18O值锆石以新元古代岩浆锆石为主, 但部分岩浆锆石在印支期超高压变质作用过程中发生了不同程度的重结晶作用.δ18O值为-7.08‰的岩浆核锆石UPb定年得到的花岗片麻岩原岩谐和年龄和不一致线上交点年龄分别为(760±49) Ma和(751±27) Ma, 变质谐和年龄和不一致线下交点年龄分别为(232±4) Ma和(241±33) Ma, 指示其原岩为新元古代花岗岩并经历了印支期变质作用; (2) 锆石δ18O值在局部范围内变化于-7.76‰~5.40‰之间, 低δ18O值岩浆锆石表明它们是从新元古代低δ18O值岩浆中直接结晶形成, 锆石δ18O值的局部变化表明其原岩岩浆的氧同位素组成具有不均一性, 指示低δ18O值岩浆源区物质曾经在地表与极度亏损18O的大气降水发生过不同程度的高温水岩反应; (3) 低δ18O值花岗片麻岩在印支期板块俯冲和折返过程中基本没有与外部发生显著的氧同位素交换, 在退变质作用过程中花岗片麻岩内部缓冲流体对原岩岩浆锆石的δ18O值影响不大.威海皂埠镇地区和东海青龙山地区的花岗片麻岩在原岩时代、变质时代和氧同位素组成等方面基本相同, 指示它们应具有相同的原岩性质, 并经历了相同的变质作用和水岩相互作用过程.因此, 极度亏损18O的新元古代双峰式基性-酸性岩浆岩可能分布于整个大别-苏鲁造山带.      

Rate of growth of the preserved North American continental crust: Evidence from Hf and O isotopes in Mississippi detrital zircons

Detrital zircons from the Mississippi River have been analyzed for U-Th-Pb, Lu-Hf and O isotopes to constrain the rate of growth of the preserved North American continental crust. One hundred and forty two concordant zircon U/Pb dates on grains mounted in epoxy, obtained by Excimer laser ablation ICP-MS method, resolved six major periods of zircon crystallization: 0-0.25, 0.3-0.6, 0.95-1.25, 1.3-1.5, 1.65-1.95 and 2.5-3.0 Ga. These age ranges match the ages of the recognized lithotectonic units of the North American continent in the hinterland of the Mississippi River. Ninety-six zircons mounted on tape, which show no age zonation and were within 7.5% of concordance, were selected to represent the six U/Pb age time intervals and analyzed for Lu-Hf and O isotope by laser ablation MC-ICP-MS and SHRIMP II, respectively. The δ¹⁸O values of the zircons show a small step increase in the maximum δ¹⁸O values at the Archean-Proterozoic boundary from 7.5‰ in the Archean to 9.5‰, and rarely 13‰, in the Proterozoic and Phanerozoic. However, the average value of δ¹⁸O in zircons changes little with time, showing that the increase in the maximum δ¹⁸O values between 2.5 and 2.0 Ga, which can be attributed to an increase in the sediment content of the source regions of younger granitoids, is largely balanced by an increase in zircons with anomalously low δ¹⁸O, which can be attributed to hydrothermally altered crust in the granitoid source region.εHf_i values for the zircons range from 13.1 to −26.9. Zircons derived from juvenile crust, which we define as having mantle δ¹⁸O (4.5-6.5‰) and lying within error of the Hf depleted mantle growth curve, are rare or absent in the Mississippi basin. The overwhelming majority of zircons crystallized from melted pre-existing continental crust, or mantle-derived magmas that were contaminated by continental crust. The average time difference between primitive crust formation and remelting for each of the recognized lithotectonic time intervals, which is defined as crustal incubation time in this study, is 890 ± 460 Myr. There is also a suggestion that the crustal incubation time increases with decreasing age in the Mississippi basin, which is consistent with the declining role of radioactive heat production in the lower crust with time.The average Hf model age (1.94 Ga), weighted by fraction of zircons in the river load is in reasonable agreement with the Nd model age (1.7 Ga) for the Mississippi River. However, if the zircons are weighted by the area of North America covered by the six recognized periods of zircon crystallization the average model age is 2.35 Ga, which compares favorably with an area weighted Nd model age of 2.36 Ga. Our preferred approach is to use the measured O isotope values to constrain variations in the ¹⁷⁶Lu/¹⁷⁷Hf ratio of the granitic source region from which the zircons crystallized, making the assumption that zircons with mantle-like O isotopic ratios have higher ¹⁷⁶Lu/¹⁷⁷Hf than zircons with higher O isotope values. This method gives an average Hf model age of 2.53 Ga, which is 180 Myr older than the constant ¹⁷⁶Lu/¹⁷⁷Hf calculation.The area weighted zircon Hf model ages show two distinct periods of crust formation for the North American continent, 1.6-2.2 and 2.9-3.4 Ga. At least 50% of the preserved North American continental crust was extracted from the mantle by 2.9 Ga and 90% by 1.6 Ga. Two similar periods of crustal growth are also recognized in Gondwana (Hawkesworth C. J. and Kemp A. I. S. (2006) Using hafnium and oxygen isotopes in zircons to unravel the record of crustal evolution. Chem. Geol.226, 144-162.), suggesting that these may be periods of global continental crustal growth. However, we stress that more data from other continents are required before the hypothesis of episodic global continental growth can be accepted with confidence.     

Unusual Uranium‐rich Fluid Flow During Exhumation of Deeply Subducted Continental Crust

In‐situ SIMS analyses of O and U‐Pb isotopes were carried out for zircons from a quartz vein hosted by ultrahigh‐pressure metagranite (UHP) in the Dabie orogen. The results are integrated to decipher the property of unusual U‐rich aqueous fluids and their effects on both metamorphic and magmatic zircons during exhumation of the UHP metagranite. In CL images, most zircon grains show distinct core‐rim structures. Relict cores are bright and exhibit oscillatory or patchy zonation, giving Neoproterozoic upper‐intercept ages of 795 ± 26 Ma. Newly grown rims are dark and exhibit no zoning, yielding Triassic concordant ages of 215 ± 5 Ma. The cores give Th contents of 59 to 463 ppm and U contents of 98 to 558 ppm, with Th/U ratios of 0.263 to 1.423. The rims yield reduced Th contents of 11 to 124 ppm but significantly elevated U contents of 1051 to 3531 ppm, with Th/U ratios of 0.010 to 0.035. Comparison with the cores of magmatic origin, the unusual enrichment in U but depletion in Th in the rims of metamorphic origin are interpreted as zircon growth from Cl‐rich oxidized vein‐forming aqueous fluids that were produced by dehydration reactions of the wallrock during continental exhumation. The cores have variably positive δ¹⁸O values with concordant or discordant Neoproterozoic U‐Pb ages, suggesting their solid‐state modification of both O and U‐Pb isotopes through interaction with the fluids. The rims yield negative δ¹⁸O values, indicating their growth from the negative δ¹⁸O fluids. Taken together, the proposed Cl‐rich oxidized negative‐δ¹⁸O vein‐forming aqueous fluids have such an ability to not only cause variable metamorphic recrystallization in the relict magmatic zircons but also produce dramatic fractionation of U over Th in the metamorphic zircons during quartz veining, and potentially impact on the overlain metasomatite in the mantle wedge.     

Iron and magnesium isotopic compositions of peridotite xenoliths from Eastern China

We present high-precision measurements of Mg and Fe isotopic compositions of olivine, orthopyroxene (opx), and clinopyroxene (cpx) for 18 lherzolite xenoliths from east central China and provide the first combined Fe and Mg isotopic study of the upper mantle. δ⁵⁶Fe in olivines varies from 0.18‰ to −0.22‰ with an average of −0.01 ± 0.18‰ (2SD, n = 18), opx from 0.24‰ to −0.22‰ with an average of 0.04 ± 0.20‰, and cpx from 0.24‰ to −0.16‰ with an average of 0.10 ± 0.19‰. δ²⁶Mg of olivines varies from −0.25‰ to −0.42‰ with an average of −0.34 ± 0.10‰ (2SD, n = 18), opx from −0.19‰ to −0.34‰ with an average of −0.25 ± 0.10‰, and cpx from −0.09‰ to −0.43‰ with an average of −0.24 ± 0.18‰. Although current precision (∼±0.06‰ for δ⁵⁶Fe; ±0.10‰ for δ²⁶Mg, 2SD) limits the ability to analytically distinguish inter-mineral isotopic fractionations, systematic behavior of inter-mineral fractionation for both Fe and Mg is statistically observed: Δ⁵⁶Fe_ol-cpx = −0.10 ± 0.12‰ (2SD, n = 18); Δ⁵⁶Fe_ol-opx = −0.05 ± 0.11‰; Δ²⁶Mg_ol-opx = −0.09 ± 0.12‰; Δ²⁶Mg_ol-cpx = −0.10 ± 0.15‰. Fe and Mg isotopic composition of bulk rocks were calculated based on the modes of olivine, opx, and cpx. The average δ⁵⁶Fe of peridotites in this study is 0.01 ± 0.17‰ (2SD, n = 18), similar to the values of chondrites but slightly lower than mid-ocean ridge basalts (MORB) and oceanic island basalts (OIB). The average δ²⁶Mg is −0.30 ± 0.09‰, indistinguishable from chondrites, MORB, and OIB. Our data support the conclusion that the bulk silicate Earth (BSE) has chondritic δ⁵⁶Fe and δ²⁶Mg.The origin of inter-mineral fractionations of Fe and Mg isotopic ratios remains debated. δ⁵⁶Fe between the main peridotite minerals shows positive linear correlations with slopes within error of unity, strongly suggesting intra-sample mineral-mineral Fe and Mg isotopic equilibrium. Because inter-mineral isotopic equilibrium should be reached earlier than major element equilibrium via chemical diffusion at mantle temperatures, Fe and Mg isotope ratios of coexisting minerals could be useful tools for justifying mineral thermometry and barometry on the basis of chemical equilibrium between minerals. Although most peridotites in this study exhibit a narrow range in δ⁵⁶Fe, the larger deviations from average δ⁵⁶Fe for three samples likely indicate changes due to metasomatic processes. Two samples show heavy δ⁵⁶Fe relative to the average and they also have high La/Yb and total Fe content, consistent with metasomatic reaction between peridotite and Fe-rich and isotopically heavy melt. The other sample has light δ⁵⁶Fe and slightly heavy δ²⁶Mg, which may reflect Fe-Mg inter-diffusion between peridotite and percolating melt.     

Zircons from kimberlite: New insights from oxygen isotopes, trace elements, and Ti in zircon thermometry

Zircons found in mantle-sourced kimberlite provide probes into the isotopic chemistry of the asthenosphere and subcontinental lithospheric mantle. However, little is known about the conditions of formation of these zircons. A suite of 88 zircons found in kimberlites from Africa, Siberia, Brazil, and the United States have been analyzed for their Ti concentration and selected zircons were analyzed for their Rare Earth Element (REE) concentrations by ion microprobe. In addition, precise and accurate laser-fluorination oxygen isotope data were obtained for zircons from Brazil (5.1 ± 0.3‰, 1SD) and the Midwest United States (5.3 ± 0.3‰), yielding mantle-like δ¹⁸O values similar to published data for Africa (5.2 ± 0.3‰) and Siberia (5.3 ± 0.2‰). Most megacrysts in this study preserve fine-scale, oscillatory zoning in CL and are generally homogenous in oxygen isotopic composition, consistent with preservation of primary compositions. A few zircons from Brazil show some evidence of chemical zoning due to recrystallization. The Ti content of mantle zircon is in general low with average compositions from each locality of 13 ± 8.4 ppm (1SD, Kaapvaal craton), 12 ± 8.7 ppm (Siberian platform), 18 ± 11 ppm (Brazil), and 4.8 ±4.3 ppm (United States). The recently calibrated Ti in zircon thermometer yields an average temperature of 744 ±62 °C (1SD) for the average of 13 ± 9 ppm Ti, with no correction for pressure, a_TiO2, or a_SiO2. The Ti content of zircons found within rutile nodules from the Orapa kimberlite (Kaapvaal craton) is almost indistinguishable from those with no constraint on a_TiO2, suggesting that reduced a_TiO2 is not responsible for lower than expected mantle temperatures. The average temperature in this study corresponds to ∼3 GPa on a 40 mW/m² cratonic geotherm. If correct, this would suggest that zircon megacrysts from all four cratons formed in the shallow lithospheric mantle. However, there are several possibly confounding effects to this thermometer, including: a pressure correction and disequilibrium zircon growth. Zircons from rutile nodules have REE contents that span the range of mantle zircon REE and are similar to both zircon megacrysts and zircons from metasomatic assemblages.     

Environmental memory and a possible seasonal bias in the stable isotope composition of (U-Th)/He-dated goethite from the Canadian Arctic

Goethite (Ax-2) from Axel Heiberg Island (∼80°N) on the margin of the Arctic Ocean is the dominant mineral in a sample of “petrified” Eocene wood, but U, Th, and He measurements suggest that the goethite (α-FeOOH) crystallized in the latest Miocene/Pliocene (ca. 5.5 to 2.8 Ma). Measured δD and δ¹⁸O values of Ax-2 are −221 (±6)‰ and −9.6 (±0.5)‰, respectively. The inferred δD and δ¹⁸O values of the ancient water were about −139‰ and −18.6‰, respectively, with a calculated temperature of crystallization of 3 (±5)°C, which compares with the modern summer (J-J-A) temperature of 3 °C and contrasts with a modern MAT of −19 °C. Published results from various biological proxies on nearby Ellesmere Island indicate a Pliocene (∼4 Ma) MAT of either −6 or −0.4 °C and corresponding seasonal amplitudes of about 18 or 13 °C. A conductive heat flow model suggests that a temperature of 3 °C could represent goethite crystallization at depths of ∼100-200 cm (for MAT = −6 °C) or ∼250-450 cm (for MAT = −0.4 °C) over seasonally restricted intervals of time.The δ¹⁸O value of the Ax-2 water (−18.6‰) is more positive than the modern J-J-A precipitation (−22‰). In combination, the paleotemperatures and δ¹⁸O values of ancient waters (from Ax-2 and published results from three Eocene or Pliocene proxy sites on Axel Heiberg and Ellesmere Islands) are consistent with a warm season bias in those isotopic proxies. The results are also consistent with higher proportions of J-J-A precipitation in the annual total. If so, this emphasizes the importance of seasonality at high latitudes even in times of warmer global climates, and suggests that the Arctic hydrologic cycle, as expressed in the seasonal distribution and isotopic composition of precipitation (perhaps modified by a warmer Arctic Ocean), differed from modern.The δ¹³C value of the Fe(CO₃)OH component in the Ax-2 goethite is +6.6‰, which is much more positive than expected if crystallizing goethite incorporated CO₂ derived primarily from oxidation of relict Eocene wood with δ¹³C values of about −24‰. This apparent paradox may be resolved if the goethite is a product of oxidation of ¹³C-rich siderite, which had previously replaced wood in an Eocene methanogenic burial environment. Thus, the goethite retains a carbon isotope “memory” of a diagenetic Eocene event, but a δD and δ¹⁸O record of the latest Miocene/Pliocene Arctic climate.     

Zircon U–Pb ages, Hf and O isotopes constrain the crustal architecture of the ultrahigh-pressure Dabie orogen in China

The crustal structure of the Dabie orogen was reconstructed by a combined study of U–Pb ages, Hf and O isotope compositions of zircons from granitic gneiss from North Dabie, the largest lithotectonic unit in the orogen. The results were deciphered from metamorphic history to protolith origin with respect to continental subduction and exhumation. Zircon U–Pb dating provides consistent ages of 751 ± 7 Ma for protolith crystallization, and two group ages of 213 ± 4 to 245 ± 17 Ma and 126 ± 4 to 131 ± 36 Ma for regional metamorphism. Majority of zircon Hf isotope analyses displays negative ε_Hf(t) values of − 5.1 to − 2.9 with crust Hf model ages of 1.84 to 1.99 Ga, indicating protolith origin from reworking of middle Paleoproterozoic crust. The remaining analyses exhibit positive ε_Hf(t) values of 5.3 to 14.5 with mantle Hf model ages of 0.74 to 1.11 Ga, suggesting prompt reworking of Late Mesoproterozoic to Early Neoproterozoic juvenile crust. Zircon O isotope analyses yield δ¹⁸O values of − 3.26 to 2.79‰, indicating differential involvement of meteoric water in protolith magma by remelting of hydrothermally altered low δ¹⁸O rocks. North Dabie shares the same age of Neoproterozoic low δ¹⁸O protolith with Central Dabie experiencing the Triassic UHP metamorphism, but it was significantly reworked at Early Cretaceous in association with contemporaneous magma emplacement. The Rodinia breakup at about 750 Ma would lead to not only the reworking of juvenile crust in an active rift zone for bimodal protolith of Central Dabie, but also reworking of ancient crust in an arc-continent collision zone for the North Dabie protolith. The spatial difference in the metamorphic age (Triassic vs. Cretaceous) between the northern and southern parts of North Dabie suggests intra-crustal detachment during the continental subduction. Furthermore, the Dabie orogen would have a three-layer structure prior to the Early Cretaceous magmatism: Central Dabie in the upper, North Dabie in the middle, and the source region of Cretaceous magmas in the lower.     

Effects of aridity and vegetation on plant-wax δD in modern lake sediments

We analyzed the deuterium composition of individual plant-waxes in lake sediments from 28 watersheds that span a range of precipitation D/H, vegetation types and climates. The apparent isotopic fractionation (ε_a) between plant-wax n-alkanes and precipitation differs with watershed ecosystem type and structure, and decreases with increasing regional aridity as measured by enrichment of ²H and ¹⁸O associated with evaporation of lake waters. The most negative ε_a values represent signatures least affected by aridity; these values were −125 ± 5‰ for tropical evergreen and dry forests, −130‰ for a temperate broadleaf forest, −120 ± 9‰ for the high-altitude tropical páramo (herbs, shrubs and grasses), and −98 ± 6‰ for North American montane gymnosperm forests. Minimum ε_a values reflect ecosystem-dependent differences in leaf water enrichment and soil evaporation. Slopes of lipid/lake water isotopic enrichments differ slightly with ecosystem structure (i.e. open shrublands versus forests) and overall are quite small (slopes = 0-2), indicating low sensitivity of lipid δD variations to aridity compared with coexisting lake waters. This finding provides an approach for reconstructing ancient precipitation signatures based on plant-wax δD measurements and independent proxies for lake water changes with regional aridity. To illustrate this approach, we employ paired plant-wax δD and carbonate-δ¹⁸O measurements on lake sediments to estimate the isotopic composition of Miocene precipitation on the Tibetan plateau.     

In situ U-Pb, O and Hf isotopic compositions of zircon and olivine from Eoarchaean rocks, West Greenland: New insights to making old crust

The sources and petrogenetic processes that generated some of the Earth’s oldest continental crust have been more tightly constrained via an integrated, in situ (U-Pb, O and Hf) isotopic approach. The minerals analysed were representative zircon from four Eoarchaean TTG tonalites and two felsic volcanic rocks, and olivine from one harzburgite/dunite of the Itsaq Gneiss Complex (IGC), southern West Greenland. The samples were carefully chosen from localities with least migmatisation, metasomatism and strain. Zircon was thoroughly characterized prior to analysis using cathodoluminescence, scanning electron, reflected and transmitted light imaging. The zircon from all but one sample showed only minor post-magmatic recrystallisation. ²⁰⁷Pb/²⁰⁶Pb dating of oscillatory-zoned zircon using SHRIMP RG (n = 142) indicates derivation of the felsic igneous rocks from different batches of magma at 3.88, 3.85, 3.81, 3.80 and 3.69 Ga.Analyses of ¹⁸O/¹⁶O compositions of olivine from a harzburgite/dunite (n = 8) using SHRIMP II in multi-collector mode, indicate that the oxygen isotopic composition of this sample of Eoarchaean mantle (δ¹⁸O_Ol = 6.0 ± 0.4‰) was slightly enriched in ¹⁸O, but not significantly different from that of the modern mantle. Zircon δ¹⁸O measurements from the six felsic rocks (n = 93) record mean or weighted mean compositions ranging from 4.9 ± 0.7‰ to 5.1 ± 0.4‰, with recrystallised domains showing no indication of oxygen isotopic exchange during younger tectonothermal events. δ¹⁸O_Zr compositions indicate that the primary magmas were largely in equilibrium with the mantle or mantle-derived melts generated at similar high temperatures, while calculated tonalite δ¹⁸O_WR compositions (6.7-6.9‰) resemble those of modern adakites.LA-MC-ICPMS zircon ¹⁷⁶Hf/¹⁷⁷Hf analyses were obtained from six samples (n = 122). Five samples record weighted mean initial ε_Hf compositions ranging from to 0.5 ± 0.6 to −0.1 ± 0.7 (calculated using λ¹⁷⁶Lu = 1.867 × 10⁻¹¹ yr⁻¹), while one sample records a composition of 1.3 ± 0.7, indicating the magmas were generated from a reservoir with a time averaged, near chondritic Lu/Hf. The derivation of TTG magmas from a chondritic Lu/Hf source implies either that there was not voluminous continental crustal growth nor major mantle differentiation leading to Lu/Hf fractionation during the Hadean or Eoarchaean, or alternatively that rapid recycling of an early formed crust allowed the early mantle to maintain a chondritic Lu/Hf.Previous studies have demonstrated that ancient TTG rocks were mostly produced by dehydration melting of mafic rocks within the stability field of garnet, probably in flatly-subducted or buried oceanic crust. The oxygen isotopic signatures measured here at high spatial resolution allow the source materials to be better defined. Melting of a mixed mafic source consisting of ∼80% unaltered gabbro (δ¹⁸O_WR = 5.5‰) with ∼20% hydrothermally altered gabbro/basalt (δ¹⁸O_WR = 4.0‰) would produce tonalite magmas within the average compositional range observed. ¹⁸O-enriched components such as altered shallow basaltic oceanic crust and pelagic or continental sediments were not present in the sources of these TTG melts. The absence of high ¹⁸O signatures may indicate either the rarity of low temperature altered sediments, or their effective removal from the down-going slab.     

Growth rate of the preserved continental crust: II. Constraints from Hf and O isotopes in detrital zircons from Greater Russian Rivers

Detrital zircons from the Ob, Yenisey, Lena, Amur, Volga, Dnieper, Don and Pechora rivers have been analyzed for U-Th-Pb, O and Lu-Hf isotopes to constrain the growth rate of the preserved continental crust in Greater Russia. Four major periods of zircon crystallization, 0.1-0.55, 0.95-1.3, 1.45-2.0 and 2.5-2.9 Ga, were resolved from a compilation of 1366 zircon U/Pb ages. The Archean zircons have δ¹⁸O values lying between 4.53‰ and 7.33‰, whereas Proterozoic and Phanerozoic zircons have a larger range of δ¹⁸O values in each of the recognized U/Pb time intervals with maximum δ¹⁸O values up to 12‰. We interpret the zircons with δ¹⁸O between 4.5‰ and 6.5‰ to have been derived from a magmatic precursor that contains little or no sedimentary component. The variable δ¹⁸O values of the zircons were used to constrain the ¹⁷⁶Lu/¹⁷⁷Hf ratios of the crustal source region of the zircons, which, in turn, were used to calculate Hf model ages (T_DM^V). The crustal incubation time, the time difference between primitive crust formation (dated by T_DM^V) and crustal melting (dated by zircon U/Pb age), varies between 300 to 1000 Myr for the majority of analyzed zircon grains, but can be up to 2500 Myr. The average T_DM^V Hf model age weighted by the fraction of zircons in the river load is 2.12 Ga, which is in reasonable agreement with the area-weighted average of 2.25 Ga. The T_DM^V Hf model age crustal growth curve for zircons with mantle-like δ¹⁸O values (4.5-6.5‰), weighted by area, shows that growth of the Great Russian continental crust started at 4.2 Ga, and that there are two principal periods of crustal growth, 3.6-3.3 Ga and 0.8-0.6 Ga, which are separated by an interval of low but more or less continuous growth. An alternative interpretation, in which the average ¹⁷⁶Lu/¹⁷⁷Hf ratio (0.0115) of the continental crust is used for the Paleoproterozoic zircons from the Lena River, lowers the average T_DM^V age of these grains by about 500 Myr and delays the onset of significant crustal growth to 3.5 Ga.The two principal growth periods recognized in Greater Russia differ from those identified from the Gondwana and the Mississippi river basin, which show peaks at 1.7-1.9 and 2.9-3.1 Ga (Hawkesworth and Kemp, 2006a) and 1.6-2.2 and 2.9-3.4 Ga (Wang et al., 2009), respectively. The older 3.6-3.3 Ga or 3.5-3.3 Ga peak for Greater Russia is slightly older than the older Gondwana-Mississippi peaks, whereas the younger 0.8-0.6 Ga peak is distinctly younger than the youngest peak in either Gondwana or the Mississippi river basin. This suggests that the two major peaks of crustal growth identified in Gondwana and the Mississippi river basin may not be global periods of enhanced continental growth and that the major periods of crustal growth may differ from continent to continent.     

Nitrogen isotope fractionations during progressive metamorphism: A case study from the Paleozoic Cooma metasedimentary complex, southeastern Australia

The well-studied Paleozoic Cooma metamorphic complex in southeastern Australia is characterized by a uniform siliciclastic protolith, of uniform age, with a continuous range of metamorphic grade from subgreenschist- to upper amphibolite-facies, and migmatite-grade in an annular pattern around the Cooma granodiorite. Those conditions are optimal for investigating variations of N concentrations and δ¹⁵N values during progressive metamorphism. Nitrogen concentrations decrease and δ¹⁵N increases with increasing metamorphic grade (sub-chlorite zone: 120 ppm N, δ¹⁵N = 2.3‰; chlorite zone: 110 ppm N, δ¹⁵N = 3.0‰; biotite and andalusite zone: 85 ppm N, δ¹⁵N = 3.8 ‰; sillimanite and migmatite zones: 40 ppm N, δ¹⁵N = 10.7‰). Covariation of K and N contents is consistent with N substituting for K as NH₄⁺ in micas. Observed trends of increasing δ¹⁵N values with decreasing nitrogen concentrations can be explained by a continuous release of nitrogen depleted in ¹⁵N with progressive metamorphism, which causes an enrichment of ¹⁵N in the residual nitrogen of the rock. Equilibrium models for Rayleigh distillation and batch volatilisation for data of the greenschist and amphibolite facies metasedimentary rocks can be explained by N₂-NH₄⁺ exchange at temperatures of 300-600 °C, whereas observed large fractionations for the upper amphibolite-facies and melt products in the migmatite-grade samples may be interpreted as NH₃-NH₄⁺ exchanges at temperature of 650-730 °C. Lower values in the highest grade zones may also stem in part from input of ¹⁵N-depleted fluids from the granodiorite.The magnitude of isotope fractionation of nitrogen is about 1-2‰ during progressive metamorphism of metasedimentary rocks from sub-chlorite zone to biotite-andalusite zone, which is consistent with previous studies. Consequently, the large spread of δ¹⁵N values in Archean greenschist-facies metasedimentary rocks of −6‰ to 30‰ can be accounted for by variable mixtures of mantle plume-dominated volatiles with a δ¹⁵N of −5‰, and a ¹⁵N-enriched marine sedimentary kerogen component inherited from a CI chondrite veneer having δ¹⁵N of 30‰ to 42‰.     

北苏鲁威海地区超基性岩的原岩形成时代和超高压变质时代

超基性岩是苏鲁超高压变质地体中一类特殊且十分重要的岩石类型,它们通常呈规模不一的块状、条带状或不规则透镜状 (体) 赋存于区域大面积出露的花岗质片麻岩中。锆石中矿物包体激光拉曼测试、阴极发光图像分析和不同性质锆石LA-ICP-MS U-Pb定年等综合研究结果表明,北苏鲁威海地区含橄榄石辉石岩 (样品W1和W2) 中锆石的成因十分复杂,可进一步划分3种不同类型锆石。其中第一类锆石呈自形-半自形晶,阴极发光图像显示清晰的岩浆结晶环带,矿物包体主要为Ol+Cpx+Ap, 记录的²⁰⁷Pb/²⁰⁶Pb年龄为1835~1845Ma,应代表含橄榄石辉石岩的原岩形成时代;第二类为变质重结晶锆石,呈半自形-他形晶,阴极发光图像显示模糊的岩浆结晶环带,矿物包体与第一类完全一致,记录的²⁰⁶Pb/²³⁸U年龄变化范围大,为250~784Ma之间,表明部分继承性岩浆结晶锆石明显受到后期岩浆-变质热事件的影响而发生不完全重结晶和Pb丢失,进而使其记录的年龄相对偏新;第三类锆石呈他形晶,为典型的变质锆石,阴极发光图像十分均匀,矿物包体相对少见,主要为Grt+Cpx,记录的²⁰⁶Pb/²³⁸U年龄为230~234Ma, 且与苏鲁地体榴辉岩及其围岩中含柯石英锆石微区记录的超高压变质年龄 (225~235Ma) 十分一致,应代表含橄榄石辉石岩的峰期超高压变质时代。超基性岩中超高压变质锆石的准确识别表明苏鲁地体在峰期超高压变质阶段的确存在流体,流体的存在对超高压变质锆石的形成起着至关重要的作用。该项研究不仅准确厘定北苏鲁威海地区超基性岩的原岩形成时代和超高压变质时代,而且对于深入探讨苏鲁-大别超高压地体流体行为、演化规律及其水-岩相互作用机理具有重要的科学意义。     

O/O and D/H in goethite from a North American Oxisol of the Early Eocene climatic optimum

An Early Eocene Oxisol in the Ione Formation of California formed in a coastal continental weathering environment at a paleolatitude of ∼38°N. The dominant minerals in the Oxisol are goethite, quartz, and kaolinite. Material balance calculations were applied to new measurements of chemical composition, D/H, and ¹⁸O/¹⁶O ratios of Oxisol samples to determine the δD (−150 ± 3‰) and δ¹⁸O (−2.4 ± 0.3‰) values of the goethite (α-FeOOH). These data, in combination with the global meteoric water line (MWL), yielded an isotopic temperature of 21(±4) °C. The nominal value of 21 °C contrasts with the modern mean annual temperature (MAT) of 16 °C in that area. The warmer temperature is consistent with formation of the goethite during the Early Eocene climatic optimum. The isotopic composition of the goethite and a temperature of 21 °C imply ancient water with a δD value of −61(±4)‰ and a δ¹⁸O value of −8.9(±0.5)‰. This Early Eocene δ¹⁸O (or δD) value is more negative than values in the range of isotopic scatter observed for modern global precipitation at sites with a MAT of 21 °C.At times of warm global climates, the location of a near-surface atmospheric isotherm would generally shift relative to its location under modern climatic conditions. A simple Rayleigh-type condensation model indicates that, if one “follows the isotherm”, the associated scatter in δD and δ¹⁸O of precipitation in very warm global climates should shift (for a given isotherm) to more negative values that may be detectable in proxy records. The isotopic results from the goethite of the Early Eocene Oxisol appear to add to evidence in support of this idea.     

Metamorphic growth and recrystallization of zircon: Distinction by simultaneous in-situ analyses of trace elements,U–Th–Pb and Lu–Hf isotopes in zircons from eclogite-facies rocks in the Sulu orogen

Simultaneous in-situ analyses of trace elements, U–Th–Pb and Lu–Hf isotopes were carried out on distinct domains of zircons in ultrahigh-pressure (UHP) eclogite-facies metamorphic rocks from the main hole of the Chinese Continental Scientific Drilling (CCSD) in the Sulu orogen. For the first time, trace elements are directly linked to Lu–Hf isotopes in metamorphic zircons with reference to their U–Pb dates. This enables methodological integration to distinguish four types of metamorphic zircon: solid-state, replacement and dissolution recrystallizations of protolith zircons, and new growth from the aqueous fluid. Metamorphically grown zircons are characterized by concordant U–Pb ages for the metamorphism, flat HREE patterns typical of the garnet effect, low contents of REE (especially HREE), Y, Nb + Ta and Th + U, high contents of Hf, low (Lu/Gd)_N, Lu/Hf and Th/U (< 0.1) ratios, and elevated ¹⁷⁶Hf/¹⁷⁷Hf ratios relative to solid-state recrystallized zircons. This suggests the effects of both garnet and fluid on the growth of metamorphic zircons. In contrast, metamorphic recrystallization has reset the U–Th–Pb isotope system of protolith zircons to different extents, depending on the extents of fluid action during metamorphism. Solid-state recrystallized zircons exhibit the lowest degrees of resetting and thus almost inherit all geochemical features from the protolith zircons, which are characterized by discordant U–Pb ages close to or below the protolith age, steep MREE–HREE patterns typical of magmatic origin, high contents of trace elements and their ratios, and low ¹⁷⁶Hf/¹⁷⁷Hf ratios. On the other hand, dissolution recrystallized zircons show the highest degrees of reworking and thus have concordant or nearly concordant U–Pb ages for the metamorphism, steep MREE–HREE patterns, lowered contents of trace elements such as REE, Th, U, Y, Nb, Ta and Ti relative to the protolith zircons, and almost unchanged Hf isotope ratios. Replacement recrystallized zircons display intermediate degrees of reworking and thus have their many features of elements and isotopes in between. While the metamorphic growth in the presence of both garnet and fluid is characterized by both depletion of HREE with flat pattern and the low contents of trace elements, the metamorphic recrystallization in the presence of aqueous fluid is indicated by gradual decreases of MREE to HREE without the flat HREE pattern. Therefore, the simultaneous in-situ analyses of metamorphic zircons have the advantage over single-term analyses in making distinction between the new growth and the different types of recrystallization.     

Variations of Li and Mg isotope ratios in bulk chondrites and mantle xenoliths

We present whole rock Li and Mg isotope analyses of 33 ultramafic xenoliths from the terrestrial mantle, which we compare with analyses of 30 (mostly chondritic) meteorites. The accuracy of our new Mg isotope ratio measurement protocol is substantiated by a combination of standard addition experiments, the absence of mass independent effects in terrestrial samples and our obtaining identical values for rock standards using two different separation chemistries and three different mass-spectrometric introduction systems. Carbonaceous, ordinary and enstatite chondrites have irresolvable mean stable Mg isotopic compositions (δ²⁵Mg = −0.14 ± 0.06; δ²⁶Mg = −0.27 ± 0.12‰, 2SD), but our enstatite chondrite samples have lighter δ⁷Li (by up to ∼3‰) than our mean carbonaceous and ordinary chondrites (3.0 ± 1.5‰, 2SD), possibly as a result of spallation in the early solar system. Measurements of equilibrated, fertile peridotites give mean values of δ⁷Li = 3.5 ± 0.5‰, δ²⁵Mg = −0.10 ± 0.03‰ and δ²⁶Mg = −0.21 ± 0.07‰. We believe these values provide a useful estimate of the primitive mantle and they are within error of our average of bulk carbonaceous and ordinary chondrites. A fuller range of fresh, terrestrial, ultramafic samples, covering a variety of geological histories, show a broad positive correlation between bulk δ⁷Li and δ²⁶Mg, which vary from −3.7‰ to +14.5‰, and −0.36‰ to + 0.06‰, respectively. Values of δ⁷Li and δ²⁶Mg lower than our estimate of primitive mantle are strongly linked to kinetic isotope fractionation, occurring during transport of the mantle xenoliths. We suggest Mg and Li diffusion into the xenoliths is coupled to H loss from nominally anhydrous minerals following degassing. Diffusion models suggest that the co-variation of Mg and Li isotopes requires comparable diffusivities of Li and Mg in olivine. The isotopically lightest samples require ∼5-10 years of diffusive ingress, which we interpret as a time since volatile loss in the host magma. Xenoliths erupted in pyroclastic flows appear to have retained their mantle isotope ratios, likely as a result of little prior degassing in these explosive events. High δ⁷Li, coupled with high [Li], in rapidly cooled arc peridotites may indicate that these samples represent fragments of mantle wedge that has been metasomatised by heavy, slab-derived fluids. If such material is typically stirred back into the convecting mantle, it may account for the heavy δ⁷Li seen in some oceanic basalts.     

Correlated microanalysis of zircon: Trace element, δO, and U-Th-Pb isotopic constraints on the igneous origin of complex >3900 Ma detrital grains

The origins of >3900 Ma detrital zircons from Western Australia are controversial, in part due to their complexity and long geologic histories. Conflicting interpretations for the genesis of these zircons propose magmatic, hydrothermal, or metamorphic origins. To test the hypothesis that these zircons preserve magmatic compositions, trace elements [rare earth elements (REE), Y, P, Th, U] were analyzed by ion microprobe from a suite of >3900 Ma zircons from Jack Hills, Western Australia, and include some of the oldest detrital zircons known (4400-4300 Ma). The same ∼20 μm domains previously characterized for U/Pb age, oxygen isotope composition (δ¹⁸O), and cathodoluminescence (CL) zoning were specifically targeted for analysis. The zircons are classified into two types based on the light-REE (LREE) composition of the domain analyzed. Zircons with Type 1 domains form the largest group (37 of 42), consisting of grains that preserve evolved REE compositions typical of igneous zircon from crustal rocks. Grains with Type 1 domains display a wide range of CL zoning patterns, yield nearly concordant U/Pb ages from 4400 to 3900 Ma, and preserve a narrow range of δ¹⁸O values from 4.7‰ to 7.3‰ that overlap or are slightly elevated relative to mantle oxygen isotope composition. Type 1 domains are interpreted to preserve magmatic compositions. Type 2 domains occur in six zircons that contain spots with enriched light-REE (LREE) compositions, here defined as having chondrite normalized values of La_N > 1 and Pr_N > 10. A subset of analyses in Type 2 domains appear to result from incorporation of sub-surface mineral inclusions in the analysis volume, as evidenced by positively correlated secondary ion beam intensities for LREE, P, and Y, which are anti-correlated to Si, although not all Type 2 analyses show these features. The LREE enrichment also occurs in areas with discordant U/Pb ages and/or high Th/U ratios, and is apparently associated with past or present radiation damage. The enrichment is not attributed to hydrothermal alteration, however, as oxygen isotope ratios in Type 2 domains overlap with magmatic values of Type 1 domains, and do not appear re-set as might be expected from dissolution or ion-exchange processes operating at variable temperatures. Thus, REE compositions in Type 2 domains where mineral inclusions are not suspected are best interpreted to result from localized enrichment of LREE in areas with past or present radiation damage, and with a very low fluid/rock ratio. Correlated in situ analyses allow magmatic compositions in these complex zircons to be distinguished from the effects of secondary processes. These results are additional evidence for preservation of magmatic compositions in Jack Hills zircons, and demonstrate the benefits of detailed imaging in studies of complicated detrital zircons of unknown origin. The data reported here support previous interpretations that the majority of >3900 Ma zircons from the Jack Hills have an origin in evolved granitic melts, and are evidence for the existence of continental crust very early in Earth’s history.     

Stable isotope and petrologic evidence for open-system degassing during the climactic and pre-climactic eruptions of Mt. Mazama, Crater Lake, Oregon

Evaluation of the extent of volatile element recycling in convergent margin volcanism requires delineating likely source(s) of magmatic volatiles through stable isotopic characterization of sulfur, hydrogen and oxygen in erupted tephra with appropriate assessment of modification by degassing. The climactic eruption of Mt. Mazama ejected approximately 50 km³ of rhyodacitic magma into the atmosphere and resulted in formation of a 10-km diameter caldera now occupied by Crater Lake, Oregon (lat. 43°N, long. 122°W). Isotopic compositions of whole-rocks, matrix glasses and minerals from Mt. Mazama climactic, pre-climactic and postcaldera tephra were determined to identify the likely source(s) of H₂O and S. Integration of stable isotopic data with petrologic data from melt inclusions has allowed for estimation of pre-eruptive dissolved volatile concentrations and placed constraints on the extent, conditions and style of degassing.Sulfur isotope analyses of climactic rhyodacitic whole rocks yield δ³⁴S values of 2.8-14.8‰ with corresponding matrix glass values of 2.4-13.2‰. δ³⁴S tends to increase with stratigraphic height through climactic eruptive units, consistent with open-system degassing. Dissolved sulfur concentrations in melt inclusions (MIs) from pre-climactic and climactic rhyodacitic pumices varies from 80 to 330 ppm, with highest concentrations in inclusions with 4.8-5.2 wt% H₂O (by FTIR). Up to 50% of the initial S may have been lost through pre-eruptive degassing at depths of 4-5 km. Ion microprobe analyses of pyrrhotite in climactic rhyodacitic tephra and andesitic scoria indicate a range in δ³⁴S from −0.4‰ to 5.8‰ and from −0.1‰ to 3.5‰, respectively. Initial δ³⁴S values of rhyodacitic and andesitic magmas were likely near the mantle value of 0‰. Hydrogen isotope (δD) and total H₂O analyses of rhyodacitic obsidian (and vitrophyre) from the climactic fall deposit yielded values οf −103 to −53‰ and 0.23-1.74 wt%, respectively. Values of δD and wt% H₂O of obsidian decrease towards the top of the fall deposit. Samples with depleted δD, and mantle δ¹⁸O values, have elevated δ³⁴S values consistent with open-system degassing. These results imply that more mantle-derived sulfur is degassed to the Earth’s atmosphere/hydrosphere through convergent margin volcanism than previously attributed. Magmatic degassing can modify initial isotopic compositions of sulfur by >14‰ (to δ³⁴S values of 14‰ or more here) and hydrogen isotopic compositions by 90‰ (to δD values of −127‰ in this case).