Earth and Planetary Science Letters最新文献

{"title":"Probing Cretaceous-Paleogene crustal thickness in southern Tibet using quartz-zircon chronobarometry","authors":"T. Mark Harrison ,&nbsp;Casey A. Yamamoto ,&nbsp;Ming-Chang Liu","doi":"10.1016/j.epsl.2025.119413","DOIUrl":"10.1016/j.epsl.2025.119413","url":null,"abstract":"<div><div>Knowledge of the crustal thickness history of southern Tibet during the India-Asia collision is key to understanding what deformation mechanisms accommodated northward propagating crustal shortening. Thermoisotopic models suggest that a relatively thin margin of southern Asia persisted from ca. 200 to 45 Ma whereas trace element paleodepth proxies are interpreted to indicate a thinning phase from ca. 100 to 65 Ma when the crust reached ∼30 km-thickness. As neither of these methods directly measures crustal thickness, resolution of this conflict awaits development of a method that can. In this study we place bounds on the thickening history of southern Tibet using coupled quartz-zircon thermobarometry and U-Pb geochronology of granitoid plutons in the Gangdese batholith. We find that crustal thicknesses were as high (or higher) as 60–75 km at 65 Ma, or 10–15 Ma prior to the onset of hard continental collision. Magmatic inflation in the lower crust likely contributed to crustal thickening and thermal weakening of the lower crust, suggesting that pure shear in the lower crust was likely the primary accommodation mechanism for N-S shortening during collision. These data are in partial agreement with results of thermoisotopic models of crustal thickness but contradict thickness histories derived from empirical trace element proxies.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"663 ","pages":"Article 119413"},"PeriodicalIF":4.8,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143929601","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Shallow storage conditions at Krafla IDDP-1 revealed by rhyolite-MELTS geobarometry, and implications for global shallow magmatism","authors":"Lydia J. Harmon ,&nbsp;Guilherme A.R. Gualda ,&nbsp;Blake M. Wallrich ,&nbsp;Calvin F. Miller","doi":"10.1016/j.epsl.2025.119421","DOIUrl":"10.1016/j.epsl.2025.119421","url":null,"abstract":"<div><div>Identifying the storage depths of melt-dominated magma bodies prior to eruption is critical for understanding magma transport, eruption hazards, and magma body longevity. Rhyolite-MELTS has been used effectively to calculate pre-eruptive storage pressures for silicic magma bodies in the upper crust (∼100–350 MPa), but its precision and accuracy in very low-pressure systems (&lt;100 MPa) has not been sufficiently investigated. During the 2009 Krafla IDDP-1 drilling project, magma was surprisingly intersected at 2.1 km depth. Here, we test the use of rhyolite-MELTS geobarometry for this very low-pressure system, using natural Krafla IDDP-1 compositions that were intersected at a known depth. We input the composition of the melt (preserved as glass) and search in pressure and temperature and oxygen fugacity (<em>f_O2</em>) spaces to model the storage conditions of the Krafla magma. For the average composition of the drilled glass, rhyolite-MELTS yields reasonable storage pressures (∼40–50 MPa). After converting calculated pressure to depth, the calculated depths are 1.6–1.9 km. These estimates are only 0.2–0.5 km different from that of the intersected magma, showing that rhyolite-MELTS provides excellent estimates for very shallow magma storage, further strengthened by results from a Monte Carlo analysis. The agreement between rhyolite-MELTS pressures and the drilled depth of the Krafla magma supports the previously calculated very shallow storage pressures in other locations, like the Taupō Volcanic Zone (TVZ), Aotearoa New Zealand. This shallowest storage zone of melt-dominated magmas has significant implications for modeling volcanic unrest and evaluating geothermal and economic resource potential.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"663 ","pages":"Article 119421"},"PeriodicalIF":4.8,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143931816","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Scale-dependent effects of growth rate on temperature reconstruction in a marine biogenic carbonate","authors":"Ronald E. Thresher","doi":"10.1016/j.epsl.2025.119391","DOIUrl":"10.1016/j.epsl.2025.119391","url":null,"abstract":"<div><div>Analyses of a temperature proxy—Mg/Ca ratios—in the calcitic skeletons of bamboo corals (Octocorallia: Keratoisididae and Mopseidae) show that the effects of ambient temperature and growth rate on Mg/Ca (and oxygen isotope ratios) differ depending on the scale of the analysis. At the whole-of-specimen scale, Mg/Ca ratios are canonically positively correlated with temperature and influenced by growth rate-dependent entrapment and diffusion. However, within specimens, Mg/Ca is much more steeply positive with temperature and correlates negatively with growth rate. Within specimens, the slope of the (Mg/Ca)/temperature relationship itself varies with growth rate, which is hypothesized to result from micro-scale shifts in the relative abundance of slow and fast growing crystal types that, in turn, differ in their affinity for Mg++. This variable slope critically affects proxy-based temperature reconstruction, can underlie frequent mis-matches in marine calcifiers between environmental records and proxy-based reconstructions, and suggests cautious application within-specimens of proxies developed and calibrated using geographic “core top” comparisons.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"662 ","pages":"Article 119391"},"PeriodicalIF":4.8,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143923400","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Reconstruction of marine redox landscape during the Cryogenian interglacial oceans using thallium isotopes","authors":"Lulu Wang ,&nbsp;Mengchun Cao ,&nbsp;Yi-bo Lin ,&nbsp;Fei Wu ,&nbsp;Qing Tang ,&nbsp;Feifei Zhang","doi":"10.1016/j.epsl.2025.119419","DOIUrl":"10.1016/j.epsl.2025.119419","url":null,"abstract":"<div><div>The change in marine redox chemistry between the Cryogenian Sturtian and Marinoan glaciations (ca. 663–654 Ma) is crucial for understanding the potential relationship between environmental conditions and the emergence and diversification of early animals. The Nanhua Basin in South China provides a nearly complete sedimentary record of Cryogenian interglacial sedimentation, commonly referred to as the Datangpo Formation. In this study, we present high-resolution thallium isotope (ε²⁰⁵Tl) records from the basal Datangpo Formation using two drill cores in the Nanhua Basin (South China) to constrain changes in marine redox chemistry during the Cryogenian interglacial period. Our ε²⁰⁵Tl values range from −4.9 to −0.7 with a mean of −2.8, which is lower than the upper continental crust value (ε²⁰⁵Tl_UCC = −2.1 ± 0.3) but higher than the modern seawater value (ε²⁰⁵Tl_seawater = −6 ± 0.6). We observed a negative excursion in the lower part of both drill cores and suggest an episode of ocean oxygenation in the basal Datangpo Formation. The ε²⁰⁵Tl values shifted from lower to higher values upsection in both the Gaodi and Changxingpo drill cores, suggesting a return to widespread anoxic conditions. A simple mass-balance model suggests that well-oxygenated, Mn-oxide-rich sediments acted as a significant sink for Tl in the Nanhua Basin (Gaodi: <em>f</em>_oxic = 7 − 14 %, Changxingpo: <em>f</em>_oxic = 13−20 %). Our new ε²⁰⁵Tl data also provide insights for interpreting black shale-hosted Mn-carbonate deposits in the basal Datangpo Formation, with the ε²⁰⁵Tl values from the Mn-carbonate-hosted interval supporting the transformation of Mn-carbonates from Mn-oxide precursors in the wake of Sturtian glaciation. The Tl isotope data indicate that Mn-oxide-rich sediments were a significant sink for Tl in the Cryogenian interglacial oceans; if this reflects oxygenated marine conditions, this transient oxygenation event coincides with the first significant radiation of algae, which likely laid the ecological foundation for the emergence of multicellular organisms. These findings highlight the intricate link between marine redox fluctuations and early biological evolution during the Cryogenian interglacial period.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"662 ","pages":"Article 119419"},"PeriodicalIF":4.8,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143923401","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Earthquake dynamics from pseudotachylyte microstructure","authors":"Silvia Aldrighetti ,&nbsp;Giorgio Pennacchioni ,&nbsp;Giulio Di Toro","doi":"10.1016/j.epsl.2025.119424","DOIUrl":"10.1016/j.epsl.2025.119424","url":null,"abstract":"<div><div>During an earthquake, most of the dissipation of the stored elastic strain energy occurs by fracturing processes and frictional heat along and adjacent to the seismic fault. Information on earthquake energy partitioning and dynamics can be retrieved from the analysis of exhumed faults containing pseudotachylytes (solidified frictional melts produced during seismic faulting). Here, microstructural analysis is carried out on an east-west striking pseudotachylyte-bearing fault of the dextral strike-slip Gole Larghe Fault Zone within the Adamello granitoid pluton (Italian Alps), exhumed from 8–11 km depth. FESEM cathodoluminescence analysis reveals a strong fragmentation of the wall rocks, invisible with other techniques, which decays in the first centimetres from the pseudotachylyte fault-parallel vein. In the northern block the microfracture density is on average low (7448 mm^-2) and the microfractures strike preferentially E-W. In contrast, in the southern block microfracture density is on average high (12,120 mm^-2), and the microfractures strike preferentially N-S. This asymmetric wall rock damage provides evidence that the microfractures developed as result of the dynamic stress field associated with earthquake rupture propagation. The measured surface area associated with wall rock fracturing and the volume of the pseudotachylyte allow the estimate of the energy dissipated in fracturing processes <em>U</em>_S (0.015–1.88 MJ m^-2) and frictional heat <em>Q</em> (32 MJ m^-2), respectively. The comparison between <em>U</em>_S and <em>Q</em> implies that frictional heat is the major energy sink during rupture propagation in these intracontinental earthquakes.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"663 ","pages":"Article 119424"},"PeriodicalIF":4.8,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143929600","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Spatial biases in oxygen-based Phanerozoic seawater temperature reconstructions","authors":"Alexandre POHL ,&nbsp;Thomas W. WONG  HEARING ,&nbsp;Arnaud BRAYARD ,&nbsp;Ethan GROSSMAN ,&nbsp;Michael M. JOACHIMSKI ,&nbsp;Guillaume LE  HIR ,&nbsp;Thomas LETULLE ,&nbsp;Daniel J. LUNT ,&nbsp;Mathieu MARTINEZ ,&nbsp;Emmanuelle PUCEAT ,&nbsp;Guillaume SUAN ,&nbsp;Paul VALDES ,&nbsp;Yannick DONNADIEU","doi":"10.1016/j.epsl.2025.119418","DOIUrl":"10.1016/j.epsl.2025.119418","url":null,"abstract":"<div><div>Stable oxygen isotopes (δ¹⁸O) are routinely used to reconstruct sea-surface temperatures (SSTs) in the geological past, with mineral δ¹⁸O values reflecting a combination of the temperature and oxygen isotope composition of seawater (δ¹⁸O_sw). Temporal variation of mean-ocean δ¹⁸O_sw is usually accounted for following estimates of land-ice volume. Spatial variations in δ¹⁸O_sw, however, are often neglected or corrected using calibrations derived from the present-day or recent past. Geochemical methods for constraining δ¹⁸O_sw and isotope-enabled general circulation model (GCM) simulations are still technically challenging. This lack of constraints on ancient δ¹⁸O_sw is a substantial source of uncertainty for SST reconstructions. Here we use the co-variation of δ¹⁸O_sw and seawater salinity, together with GCM simulations of ocean salinity, to propose estimations of spatial variability in δ¹⁸O_sw over the Phanerozoic. Sensitivity tests of the δ¹⁸O_sw-salinity relationship and climate model, and comparison with results of isotope-enabled GCMs, suggest that our calculations are robust at first order. We show that continental configuration exerts a primary control on δ¹⁸O_sw spatial variability. Complex ocean basin geometries in periods younger than 66 Ma lead to strong inter-basinal contrasts in δ¹⁸O_sw. Latitudinal SST gradients may be steeper than previously suggested during most of the Mesozoic and Cenozoic. This work has limitations, with δ¹⁸O_sw-salinity relationships being less reliable in both low-latitude epicontinental settings and high-latitude regions of deep-water formation. Whilst our calculations are limited use in correcting δ¹⁸O measurements for local δ¹⁸O_sw, they identify the time slices and paleogeographical regions that should be prioritized for future work using isotope-enabled GCMs.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"663 ","pages":"Article 119418"},"PeriodicalIF":4.8,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143922245","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Unexpected megathrust slip evolution revealed by the 2024 Mw 7.1 and the 2025 Mw 6.8 Hyuga-nada earthquakes in southwest Japan","authors":"Xiaolong Zhang ,&nbsp;Shaoyang Li ,&nbsp;Ling Chen","doi":"10.1016/j.epsl.2025.119384","DOIUrl":"10.1016/j.epsl.2025.119384","url":null,"abstract":"<div><div>Subduction faults were thought to slip either seismically as earthquakes or aseismically through creep, representing two end-member frictional behaviors. A dynamic transition from interseismic creep to coseismic failure during rupture had also been proposed. However, the potential for a subduction fault segment to transition between seismic and aseismic mode in a time scale of days to months remains poorly understood, constrained by the remote offshore nature of these zones and limited instrumentation in most subduction regions. Here, we investigate the coseismic and time-dependent postseismic slip associated with the 2024 and 2025 megathrust events in Hyuga-nada, Japan, using dense inland GNSS data, aftershocks and a finite-element model. Our preferred source model of the 2024 event identifies a quasi-circular, thrust-dominated rupture with a geodetic moment magnitude (Mw) of 7.1 and a maximum slip of ∼1.5 m. This event occurred at relatively greater depths (15–30 km), compared to previous M7 events in the region, along the downdip edge of the inferred subducting Kyushu-Palau Ridge. Kinematic analysis of postseismic displacements indicates that the initial 50 days of afterslip were excited concurrently with ample plate-interface aftershocks in the earthquake source region and migrated downdip to the northwest into an area previously characterized by long-term and short-term slow slip events. These overlaps between coseismic and postseismic slip, as well as between postseismic slip and slow slip events, align with earlier observations from the 1996 doublet events in the vicinity of the 2024 rupture, albeit with less dense data. Approximately one month after the mainshock, an additional aseismic slip episode initiated at depths of 60–80 km, coinciding with a burst of M0–2 aftershocks in the same region. Approximately five months later, an Mw 6.8 event occurred within the main afterslip area of the 2024 event, forming a doublet with the 2024 event and largely overlapping with the rupture zone of the December 1996 event. Our findings highlight the complex, depth-dependent interplay between subducting geometrical irregularities and local megathrust rheology, which governs the partitioning of seismic and aseismic slip to accommodate the oceanic plate subduction. This interplay also appears to allow the subduction fault to transition mechanically between fast and slow slip modes on a large scale and over an unexpectedly short timescale. The unexpected dynamic slip behavior underscores the need for a systematic, global reassessment of the effects of subducting topography on megathrust slip behavior and earthquake hazards.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"662 ","pages":"Article 119384"},"PeriodicalIF":4.8,"publicationDate":"2025-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143916759","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Drainage divide formation by uplift gradient onset: The case of the southern sector of the Main Gulf Escarpment, Baja California Peninsula, Mexico","authors":"Arturo Godínez–Tamay,&nbsp;Miguel Castillo,&nbsp;Esperanza Muñoz–Salinas","doi":"10.1016/j.epsl.2025.119411","DOIUrl":"10.1016/j.epsl.2025.119411","url":null,"abstract":"<div><div>Field cases and numerical simulations support that uplift gradients set drainage divide close to the high uplift side. However, the development of a mountain landscape since the onset of one uplift gradient has been poorly documented. We analyze the case of the Main Gulf Escarpment (MGE) on the Baja California Peninsula (Western Mexico). The MGE is a landscape feature produced by El Carrizal Fault (ECF), an active normal fault that runs obliquely to the rifting system of the Gulf of California. By analyzing the rivers of the southern sector of the MGE with different approaches such as linear inversion, <em>k</em>_sn index, and drainage divide positions, we document how the different episodes along ECF tilted a volcanoclastic rock sequence known as Comondú Group, forming a main drainage divide that now corresponds to the MGE. The fault activity, which began in the Late–Miocene, was also advantageous in providing a reference time span to constrict the incision response of rivers to the episodes along ECF.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"662 ","pages":"Article 119411"},"PeriodicalIF":4.8,"publicationDate":"2025-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143923399","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"What is the present-day tectonic uplift rate of the Tibetan Plateau?","authors":"Weilong Rao ,&nbsp;Bin Liu ,&nbsp;Tengfei Feng ,&nbsp;Qiuyu Wang ,&nbsp;He Tang","doi":"10.1016/j.epsl.2025.119389","DOIUrl":"10.1016/j.epsl.2025.119389","url":null,"abstract":"<div><div>The Tibetan Plateau (TP) exhibits millimeter-scale vertical crustal deformation, detected through Global Positioning System (GPS) measurements. Concurrently, significant terrestrial water loss in High Mountain Asia contributes to plateau uplift via hydrological loading. After accounting for loading-induced uplift, it remains unclear whether the TP continues to experience significant tectonic uplift. To resolve this, a global hydrological mass change field was reconstructed using satellite observations and hydrological models, rather than relying on the Gravity Recovery and Climate Experiment (GRACE) data, to quantify vertical loading displacement rates for correcting GPS vertical velocities. Furthermore, a GPS imaging technique that incorporates the spatial structure function (SSF) was used to generate a tectonic deformation image of the TP. The results indicate that the TP is not experiencing uniform and significant tectonic uplift; instead, regional crustal subsidence is observed. The mean tectonic uplift rate is only 0.31 ± 0.32 mm/yr, which is smaller than the 0.46 ± 0.08 mm/yr contributed by loading and glacial isostatic adjustment (GIA).</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"662 ","pages":"Article 119389"},"PeriodicalIF":4.8,"publicationDate":"2025-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143922790","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Behaviour of cadmium isotopes in sulfidic waters and sediments of the Black Sea: Implications for global cadmium cycling and the application of cadmium isotopes as a paleo-oceanographic proxy","authors":"Alexander J. Dickson ,&nbsp;Allison L. Bryan ,&nbsp;Ejin George ,&nbsp;Gideon M. Henderson ,&nbsp;Donald Porcelli ,&nbsp;John Rolison ,&nbsp;Caroline P. Slomp ,&nbsp;Rob Middag ,&nbsp;Claudine H. Stirling","doi":"10.1016/j.epsl.2025.119408","DOIUrl":"10.1016/j.epsl.2025.119408","url":null,"abstract":"<div><div>The cadmium isotope system has found use as a tracer for biological productivity, redox and organic carbon burial in the oceans. There are, however, very few observational constraints on Cd isotopic behaviour in modern sulfidic marine conditions, limiting our understanding of the modern Cd cycle, and our ability to use Cd isotope measurements of ancient sedimentary deposits as a paleoceanographic proxy. Here we study the behaviour of dissolved Cd in the water column of the Black Sea and its incorporation into sediments on the basin floor. The isotopic composition of dissolved Cd in the upper ∼50 m of the Black Sea water column is controlled by a combination of biological uptake and regeneration along with mixing of river water and Mediterranean seawater. Cadmium declines to &lt;2 % of its peak subsurface concentration within the zone of nitrate-reduction but above the sulfide chemocline. The isotopic composition of dissolved Cd below the sulfide chemocline evolves in a manner that is consistent with experimental Cd sulfide fractionation factors, providing field-based support for these earlier studies. In contrast, the stability of the dissolved Cd isotopic signature above the chemocline, despite a large reduction in dissolved Cd concentrations, is explained by diffusion towards the Cd-sulfide sink below the chemocline. Sediments accumulating in the deep Black Sea below the sulfide chemocline have isotopic compositions that are similar to dissolved Cd in its near-surface aqueous concentration maximum. Sediments accumulating at shallower depths have lighter isotopic compositions that are likely affected by non-sulfidic burial pathways mediated by macrofaunal mixing and diagenesis of shallow burial phases. Our study shows that the isotopic composition of bulk sediments accumulating under sulfidic conditions is similar to dissolved Cd in the upper water column, supporting its use as a paleo-chemical proxy for ancient seawater as long as basin-scale controls on open ocean seawater sources are considered. We show that Cd burial in the deep Black Sea, predominantly as CdS, accounts for ∼0.5–15 % of total annual Cd burial in the global ocean.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"662 ","pages":"Article 119408"},"PeriodicalIF":4.8,"publicationDate":"2025-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143923398","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}