Subduction zones represent the interface between Earth’s interior (crust and mantle) and exterior (atmosphere and oceans), where carbon and other volatile elements are actively cycled between Earth reservoirs by plate tectonics. Helium is a sensitive tracer of volatile sources and can be used to deconvolute mantle and crustal sources in arcs; however it is not thought to be recycled into the mantle by subduction processes. In contrast, carbon is readily recycled, mostly in the form of carbon-rich sediments, and can thus be used to understand volatile delivery via subduction. Further, carbon is chemically-reactive and isotope fractionation can be used to determine the main processes controlling volatile movements within arc systems. Here, we report helium isotope and abundance data for 42 deeply-sourced fluid and gas samples from the Central Volcanic Zone (CVZ) and Southern Volcanic Zone (SVZ) of the Andean Convergent Margin (ACM). Data are used to assess the influence of subduction parameters (e.g., crustal thickness, subduction inputs, and convergence rate) on the composition of volatiles in surface volcanic fluid and gas emissions. He isotopes from the CVZ backarc range from 0.1 to 2.6 R A ( n = 23), with the highest values in the Puna and the lowest in the Sub-Andean foreland fold-and-thrust belt. Atmosphere-corrected He isotopes from the SVZ range from 0.7 to 5.0 R A ( n = 19). Taken together, these data reveal a clear southeastward increase in 3 He/ 4 He, with the highest values (in the SVZ) falling below the nominal range associated with pure upper mantle helium (8 ± 1 R A ), approaching the mean He isotope value for arc gases of (5.4 ± 1.9 R A ). Notably, the lowest values are found in the CVZ, suggesting more significant crustal inputs (i.e., assimilation of 4 He) to the helium budget. The crustal thickness in the CVZ (up to 70 km) is significantly larger than in the SVZ, where it is just ∼40 km. We suggest that crustal thickness exerts a primary control on the extent of fluid-crust interaction, as helium and other volatiles rise through the upper plate in the ACM. We also report carbon isotopes from ( n = 11) sites in the CVZ, where δ 13 C varies between −15.3‰ and −1.2‰ [vs. Vienna Pee Dee Belemnite (VPDB)] and CO 2 / 3 He values that vary by over two orders of magnitude (6.9 × 10 8 –1.7 × 10 11 ). In the SVZ, carbon isotope ratios are also reported from ( n = 13) sites and vary between −17.2‰ and −4.1‰. CO 2 / 3 He values vary by over four orders of magnitude (4.7 × 10 7 –1.7 × 10 12 ). Low δ 13 C and CO 2 / 3 He values are consistent with CO 2 removal (e.g., calcite precipitation and gas dissolution) in shallow hydrothermal systems. Carbon isotope fractionation modeling suggests that calcite precipitation occurs at temperatures coincident with the upper temperature limit for life (122°C), suggesting that biology may play a role in C-He systematics of arc-related volcanic fluid and gas emissions.
more »
« less
Reversal of carbonate-silicate cation exchange in cold slabs in Earth’s lower mantle
Abstract The stable forms of carbon in Earth’s deep interior control storage and fluxes of carbon through the planet over geologic time, impacting the surface climate as well as carrying records of geologic processes in the form of diamond inclusions. However, current estimates of the distribution of carbon in Earth’s mantle are uncertain, due in part to limited understanding of the fate of carbonates through subduction, the main mechanism that transports carbon from Earth’s surface to its interior. Oxidized carbon carried by subduction has been found to reside in MgCO 3 throughout much of the mantle. Experiments in this study demonstrate that at deep mantle conditions MgCO 3 reacts with silicates to form CaCO 3 . In combination with previous work indicating that CaCO 3 is more stable than MgCO 3 under reducing conditions of Earth’s lowermost mantle, these observations allow us to predict that the signature of surface carbon reaching Earth’s lowermost mantle may include CaCO 3 .
more »
« less
- Award ID(s):
- 1751664
- NSF-PAR ID:
- 10315710
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
SUMMARY Determinations of seismic anisotropy, or the dependence of seismic wave velocities on the polarization or propagation direction of the wave, can allow for inferences on the style of deformation and the patterns of flow in the Earth’s interior. While it is relatively straightforward to resolve seismic anisotropy in the uppermost mantle directly beneath a seismic station, measurements of deep mantle anisotropy are more challenging. This is due in large part to the fact that measurements of anisotropy in the deep mantle are typically blurred by the potential influence of upper mantle and/or crustal anisotropy beneath a seismic station. Several shear wave splitting techniques are commonly used that attempt resolve seismic anisotropy in deep mantle by considering the presence of multiple anisotropic layers along a raypath. Examples include source-side S-wave splitting, which is used to characterize anisotropy in the deep upper mantle and mantle transition zone beneath subduction zones, and differential S-ScS and differential SKS-SKKS splitting, which are used to study anisotropy in the D″ layer at the base of the mantle. Each of these methods has a series of assumptions built into them that allow for the consideration of multiple regions of anisotropy. In this work, we systematically assess the accuracy of these assumptions. To do this, we conduct global wavefield modelling using the spectral element solver AxiSEM3D. We compute synthetic seismograms for earth models that include seismic anisotropy at the periods relevant for shear wave splitting measurements (down to 5 s). We apply shear wave splitting algorithms to our synthetic seismograms and analyse whether the assumptions that underpin common measurement techniques are adequate, and whether these techniques can correctly resolve the anisotropy incorporated in our models. Our simulations reveal some inaccuracies and limitations of reliability in various methods. Specifically, explicit corrections for upper mantle anisotropy, which are often used in source-side direct S splitting and S-ScS differential splitting, are typically reliable for the fast polarization direction ϕ but not always for the time lag δt, and their accuracy depends on the details of the upper mantle elastic tensor. We find that several of the assumptions that underpin the S-ScS differential splitting technique are inaccurate under certain conditions, and we suggest modifications to traditional S-ScS differential splitting approaches that lead to improved reliability. We investigate the reliability of differential SKS-SKKS splitting intensity measurements as an indicator for lowermost mantle anisotropy and find that the assumptions built into the splitting intensity formula can break down for strong splitting cases. We suggest some guidelines to ensure the accuracy of SKS-SKKS splitting intensity comparisons that are often used to infer lowermost mantle anisotropy. Finally, we suggest a new strategy to detect lowermost mantle anisotropy which does not rely on explicit upper mantle corrections and use this method to analyse the lowermost mantle beneath east Asia.
-
more » « less
Rock recycling within the forearcs of subduction zones involves subduction of sediments and hydrated lithosphere into the upper mantle, exhumation of rocks to the surface, and erosion to form new sediment. The compositions of, and inclusions within detrital minerals revealed by electron microprobe analysis and Raman spectroscopy preserve petrogenetic clues that can be related to transit through the rock cycle. We report the discovery of the ultrahigh-pressure (UHP) indicator mineral coesite as inclusions in detrital garnet from a modern placer deposit in the actively exhuming Late Miocene–Recent high- and ultrahigh-pressure ((U)HP) metamorphic terrane of eastern Papua New Guinea. Garnet compositions indicate the coesite-bearing detrital garnets are sourced from felsic protoliths. Carbonate, graphite, and CO2inclusions also provide observational constraints for geochemical cycling of carbon and volatiles during subduction. Additional discoveries include polyphase inclusions of metastable polymorphs of SiO2(cristobalite) and K-feldspar (kokchetavite) that we interpret as rapidly cooled former melt inclusions. Application of elastic thermobarometry on coexisting quartz and zircon inclusions in six detrital garnets indicates elastic equilibration during exhumation at granulite and amphibolite facies conditions. The garnet placer deposit preserves a record of the complete rock cycle, operative on <10-My geologic timescales, including subduction of sedimentary protoliths to UHP conditions, rapid exhumation, surface uplift, and erosion. Detrital garnet geochemistry and inclusion suites from both modern sediments and stratigraphic sections can be used to decipher the petrologic evolution of plate boundary zones and reveal recycling processes throughout Earth’s history.
-
more » « less
Summary The Earth’s magnetic field is generated by a dynamo in the outer core and is crucial for shielding our planet from harmful radiation. Despite the established importance of the core-mantle boundary heat flux as driver for the dynamo, open questions remain about how heat flux heterogeneities affect the magnetic field. Here, we explore the distribution of core-mantle boundary heat flux on Earth and its changes over time using compressible global 3-D mantle convection models in the geodynamic modeling software ASPECT. We discuss the use of the consistent boundary flux method as a tool to more accurately compute boundary heat fluxes in finite element simulations and the workflow to provide the computed heat flux patterns as boundary conditions in geodynamo simulations. Our models use a plate reconstruction throughout the last 1 billion years—encompassing the complete supercontinent cycle—to determine the location and sinking speed of subducted plates. The results show how mantle upwellings and downwellings create localized heat flux anomalies at the core-mantle boundary that can vary drastically over Earth’s history and depend on the properties and evolution of the lowermost mantle as well as the surface subduction zone configuration. The distribution of hot and cold structures at the core-mantle boundary changes throughout the supercontinent cycle in terms of location, shape and number, indicating that these structures fluctuate and might have looked very differently in Earth’s past. We estimate the resulting amplitude of spatial heat flux variations, expressed by the ratio of peak-to-peak amplitude to average heat flux, q*, to be at least 2. However, depending on the material properties and the adiabatic heat flux out of the core, q* can easily reach values >30. For a given set of material properties, q* generally varies by 30-50% over time. Our results have implications for understanding the Earth’s thermal evolution and the stability of its magnetic field over geological timescales. They provide insights into the potential effects of the mantle on the magnetic field and pave the way for further exploring questions about the nucleation of the inner core and the past state of the lowermost mantle.
-
more » « less
Abstract Nitrogen is considered to be transported from Earth′s surface to the top of the lower mantle through subduction. However, little is known on the transportation and fate of subducted nitrogen to the Earth′s interior during slab‐mantle interactions. In this study, the stability of subducted sedimentary nitrogen in the reduced mantle was investigated to 35 GPa and 1600 K by laser‐heated diamond anvil cell experiments and first‐principles calculations. Our results showed that subducted nitrogen‐bearing silicates and fluids could not coexist with the metallic iron or iron‐rich alloys, and reacted with them to form different products at high pressure‐temperature conditions. Combining our results with previous data, we re‐determined the relative stability of iron‐light element binary compounds to 35 GPa and 1600 K to be Fe‐O > Fe‐N > Fe‐S > Fe‐C. This stability sequence contributes to explaining the observation that iron nitrides are trapped as inclusions in sulfur‐depleted lower‐mantle diamonds and are absent in sulfur‐rich ones. The recycling efficiency of subducted sedimentary nitrogen is strongly related to the availability of the metallic iron of the reduced mantle. Hydration of the metallic iron limits the storage of nitrogen in it and contributes to recycling nitrogen to Earth′s surface. Therefore, unlike subducted continental sediments, subducted marine sediments are unlikely to transport a large amount of surficial nitrogen to the metallic iron of the reduced mantle in which nitrogen could reside over long geologic periods.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41467-021-21761-9
