Article

communications earth & environment
A Nature Portfolio journal
https://doi.org/10.1038/S43247-026-03211-w
Anthropogenic emissions of volatile Cd detected in western tropical North Atlantic surface seawater
* Check for updates
Hui Xu®1 , Mark Rehkàmper®1E, Yihang Huang ®1, Zijian Jia1, Katharina Kreissig1, Barry J. Coles1, Rebekah E. T. Moore1, Arianna Olivelli1,2,3, Rob Middag ®45, Alex Baker®6, Rachel Shelley®e &
Tina van de Flierdt®1
The Cd concentrations and isotope compositions of open ocean surface waters are generally thought to be governed by internal cycling, and particularly by the balance between upwelling of more Cd-rich deeper water masses and Cd uptake by phytoplankton. Here we present a new dataset of coupled Cd isotope compositions and concentrations for seawater depth profiles sampled in the western tropical Atlantic Ocean during Leg 2 of the GEOTRACES GA02 section. A box model for the Cd source and sink fluxes of the oligotrophic surface waters of the study area shows that the observed light Cd isotope compositions and low Cd concentrations are a consequence of biological Cd uptake and atmospheric deposition of isotopically light anthropogenic Cd. Aerosols enriched in anthropogenic Cd thereby contributed at least 19%, and possibly more than 45%, to the dissolved surface water Cd inventory during the sampling period. This reveals that anthropogenic emissions of volatile Cd can have a key impact on the distribution of Cd in open ocean surface waters.

'Department of Earth Science and Engineering, Imperial College London, London, UK. 2Grantham Institute-Climate Change and the Environment, Imperial College London, London, UK. 3Flanders Marine Institute (VLIZ), Ostend, Belgium. "Department of Ocean Systems, NIOZ Royal Netherlands Institute for Sea Research, Den Burg, The Netherlands. 5Centre for Isotope Research, University of Groningen, Groningen, The Netherlands. “Centre for Ocean and Atmospheric Science, School of Environmental Sciences, University of East Anglia, Norwich, UK. e-mail: h.xu21@imperial.ac.uk; markrehk@imperial.ac.uk

Cadmium is a trace element with a marine distribution that is closely correlated with that of the nutrient PO43, as both have low dissolved concentrations in surface waters and higher concentrations at depth1-3. The similar distributions suggest that Cd may be a marine micronutrient, such that the low surface water Cd concentrations are a consequence of biological uptake. The mechanism of Cd uptake remains unclear, but evidence suggests that Cd may substitute for Zn or Co in the enzyme carbonic anhydrase under Zn-limited conditions4. Some marine diatoms have a Cd-cofactored carbonic anhydrase, allowing them to substitute Zn or Co with Cd5,6. Alternatively, or in addition, non-specific uptake of Cd can occur if phytoplankton are unable to differentiate the element from other essential divalent metals5-7. Regardless of the exact mechanism, it is evident that in the surface ocean, there is net removal of isotopically light Cd by phytoplankton89, in agreement with culture experiments10611. Closed-system Rayleigh fractionation has been used to describe the Cd isotope fractionation recorded in high nutrient low chlorophyll (HNLC) surface waters, where Cd is rapidly removed from the surface waters by biological uptake12-15. In contrast, open-system steady state conditions were employed to account for Cd isotope fractionations in regions where the removal of dissolved Cd from the surface ocean is balanced by inputs via the upwelling of deeper water masses1316-19. Following uptake and incorporation into biomass, Cd is subsequently released into the dissolved phase through remineralization in oxyclines20. The regeneration of Cd and Ρ is not always fully linked, leading to subsurface Cd:P maxima and distributions of particulate Cd that are not strictly coupled to PO43 or O221· Despite these complexities, biological uptake and remineralization remain key in shaping Cd and POf distributions in the ocean3,20.
In the deeper ocean, recent Cd isotope studies of seawater have provided new insights into the interplay between biological and physical processes in the biogeochemical cycling of Cd19,22,23. Cadmium has a long, but uncertain, oceanic residence time, estimated at between 2 x 104 years24 and 3 x 105 years25. This results in relatively uniform Cd isotope compositions for waters from depths >2000 m, with a mean 8114Cd value of 0.25 ± 0.04% (1 SD, η = 538; 8114Cd is the relative deviation of the 114Cd/110Cd ratio determined for a sample from the NIST SRM 3108 Cd isotope reference material in parts per 1000) for a global deep water dataset that excludes the North Atlantic13. The deep North Atlantic, which is dominated by North Atlantic Deep Water (NADW), has a slightly but significantly elevated 6114Cd of 0.43 ± 0.09%o (1 SD, H = 215)822,2637. This higher 8114Cd value most likely stems primarily from the shallow NADW source waters of Labrador
Sea Water, Iceland-Scotland Outflow Water, and Denmark Strait Overflow Water, which are expected to be isotopically heavy from biological uptake of Cd near the surface8'19'28'29. In addition, the 8114Cd signature of NADW may reflect inputs from Antarctic Intermediate Water and Subtropical Mode Water, as the latter have relatively high 8114Cd values, which they can imprint on NADW as they move northward and mix8,919.
While processes within the ocean’s interior influence the marine Cd distribution, external inputs, such as atmospheric deposition and riverine inputs, are unlikely to have a clear impact on basin-wide or global scale deep waters due to the long residence time of Cd in the ocean interior26'29,30. However, such inputs could affect Cd distributions in surface waters, where Cd concentrations are generally low and Cd has a much shorter residence time of months to years3132. Sieber et al. 202311 recently inferred that the relatively light Cd isotope compositions (with Sll4Cd = 0%ο) observed in low-Cd surface seawater samples from the equatorial Pacific Ocean along the GEOTRACES GP15 section reflect the combined effects of biological uptake, which lowers Cd concentrations, and atmospheric deposition, which supplies isotopically light Cd. This conclusion is supported by Cd data for leachates of 10 aerosol samples collected along the section, which are strongly enriched in Cd and characterized by δ11’Cd values of -1.91 ± 0.39% to —0.07 ± 0.25%. Such characteristics do not appear to be unusual, as aerosol samples collected along the GA03 section in the North Atlantic and the GA06 section in the Eastern Tropical Atlantic also featured high Cd contents, most likely from anthropogenic emissions, with 8114Cd values of between —0.54 ± O.13%o and 0.19 ± 0.21%3334. Notably, biological uptake of Cd in surface waters typically leaves the remaining dissolved Cd pool isotopically heavier, while the input and solubilization of anthropogenic Cd, which is often isotopically light, could balance this effect. This interplay between biological fractionation and external inputs may therefore provide an explanation for low Cd concentrations coupled with light Cd isotope signatures of surface waters.
As for riverine Cd, major rivers have variable Cd concentrations and, based on limited data, Cd isotope compositions that are influenced by local geological and biogeochemical settings, and anthropogenic inputs35. For instance, relatively pristine Siberian rivers exhibit Cd concentrations ranging from 0.02 to 0.06 nmol kg 1, with a mean 8114Cd of 0.2 ± 0.1%o, the boreal Kalix River in Sweden has a higher [Cd] of 0.24 nmol kg ’and a lower 8114Cd of —0.38 ± 0.10%o, whilst the acid mine drainage-affected Odiel and Tinto Rivers in southern Spain revealed variable Cd concentrations of between 230 and <1 nmol kg 1, depending on proximity to the estuary, but uniform 8114Cd values of about 0.02%38.36. These findings highlight the heterogeneity of riverine Cd inputs, with variations in riverine Cd concentrations significantly exceeding those observed in surface seawater

Fig. 1 | Sampling locations for seawater depth profiles analyzed in this study along the central part of the GEOTRACES section GA02 (Leg 2; cruise ΡΕ321). The map was produced using Ocean Data View'4. The arrows indicate the predominant wind directions in the area.
(1 x 10-4 to 0.8 nmol kg 1)35. Such variable riverine Cd could locally influence Cd distributions and isotope compositions in coastal and adjacent surface waters, even if their impacts on the open ocean remain limited35.
Here, new coupled Cd isotope and concentration data are presented for seawater from four stations in the western tropical North Atlantic Ocean that were sampled during Leg 2 of GEOTRACES section GA02 (cruise ΡΕ321, 13 June to 4 July 2010) to investigate the processes that impact their distributions. The new results show that horizontal mixing plays a dominant role in shaping the distribution of Cd in deep waters within this hydro- logically complex region, where water masses from the Northern and Southern Hemispheres converge. At thermocline and intermediate depths, the relative contributions of remineralization and mixing processes are evaluated to explain the variations in Cd concentrations and isotope compositions. Whilst the upwelling of Cd-rich deeper waters is generally thought to be the main source of Cd to the surface waters919, the light Cd isotope compositions found in the study area suggest that atmospheric sources, highly influenced by anthropogenic activities, strongly impacted the Cd inventory of the surface waters in this oligotrophic open ocean region.
Results
Hydrography
The hydrography of the GA02 Leg 2 section in the western North Atlantic Ocean is extensively described in previous publications337-39. In this study, we applied an adapted version of the extended optimum multiparameter (eOMP) model originally developed by Middag et al.3, which builds on the framework of Tomczak et al. (1981)40. Our adaptation incorporates five tracers—potential temperature, salinity, NO, , Si and O2 (see Supplementary Note 1), and the effects of remineralization. The model results of the eOMP analysis indicate that the main water masses present in the region are North Atlantic Central Water (NACW), South Atlantic Central Water (SACW), AAIW, Upper Circumpolar Deep Water (UCDW), NADW, and Antarctic Bottom Water (AABW) (Figs.l and 2, Supplementary Fig. 1 and Supplementary Note 1).
Cadmium concentrations and isotope compositions
The Cd concentration and isotope depth profiles for Stations 21, 30, 36, and 40 along the GA02 Leg 2 section display consistent patterns for both Cd concentrations and 8114Cd (Figs. 2 and 3; Table 1). In the top -100 m, all profiles are characterized by low Cd levels and particularly low 0l4Cd, with a subsurface 8114Cd maximum observed at between 101 and 200 m depth. In detail, the Cd concentrations of the surface waters in this study (with depths of about 25 to 100 m) range from 0.4 χ 10-3 to 3.0 x 10-3 nmol kg 1, with isotopically lighter Cd (6"'Cd between -0.42 ± 0.37% and 0.37 ± 0.17%o, n = 8) compared to deeper waters (8114Cd between 0.25 ± 0.02 and 0.63 ± 0.05, η = 38). Subsurface waters, particularly within the main thermocline and intermediate waters of NACW, SACW and AAIW/UCDW at depths of about 100 to 1000 m, exhibit Cd concentrations that generally increase with depth from 0.012 to 0.63 nmol kg 1 and show higher 8114Cd values of between 0.45 and 0.6 l%o. At depths exceeding ~1000 m, both Cd concentrations and 8114Cd align closely with water masses distribution. Northern-sourced water of NADW is characterized by Cd concentrations of 0.30 ± 0.05 nmol kg 1 and 8114Cd values of 0.41 ± Ο.Ο8%ο (2 SD, η = 14). In contrast, higher Cd concentrations of 0.61 ± 0.02 nmol kg 1 and lower 8114Cd of 0.27 ± 0.04%o (2 SD, η = 2) are found at depth at the southernmost Station 40 within the AABW layer. The impact of the AABW layer is seen to a lesser extent at the northernmost stations (Figs. 2 and 3; Table 1).
Discussion
As shown above, the GA02 seawater samples from Stations 21,30,36, and 40 along the western boundary of the North Atlantic encompass water masses from both the Northern and Southern Hemispheres (Figs. 1 and 2). The control of water mass mixing on deep water Cd concentrations and isotope compositions was described previously for the South Atlantic Ocean. For this ocean basin, the observed Cd isotope variations were shown to reflect

Fig. 2 I Distributions of a Cd concentrations, b S114Cd, c phosphate concentrations, and d oxygen concentrations, along the GA02 Leg 2 transect. Distributions of a Cd concentrations, b 8114Cd, c phosphate concentrations, and d oxygen concentrations, along the GA02 Leg 2 transect. The white contours represent different values of the variables. AABW Antarctic Bottom Water, AAIW Antarctic

Intermediate Water, NACW North Atlantic Central Water, NADW North Atlantic Deep Water, UCDW Upper Circumpolar Deep Water. The map was produced using Ocean Data View-4. The [Cd] and 8114Cd data are from this study, whilst the phosphate and oxygen concentrations are from the GEOTRACES Intermediate Data Product (IDP2021v2)48.

the mixing proportions of NADW, AABW, and AAIW, highlighting the dominant role of large-scale ocean circulation in shaping Cd isotope distributions in the deep ocean19. The current study provides a particularly clear illustration of the systematic Cd isotope variations with latitude for the deepwater flow along the western boundary of the tropical North Atlantic, where water from the Northern Hemisphere moves southward and water from the Southern Hemisphere moves northward (Fig. 2). In addition to horizontal mixing of water masses, remineralization acts to release Cd from biological particles back into the dissolved form in the oxygen minimum layers of the western Atlantic Ocean3.
To quantify the impact of horizontal mixing and remineralization in the study area, the eOMP model was coupled with Monte Carlo simulations (Supplementary Note 1)341. The relative proportion of each water mass present was obtained from the eOMP model (Supplementary Figs. 1 and 2)3. Each water mass was thereby assigned an endmember Cd concentration and 8114Cd value, based on the best available data (Supplementary Table 1). For water masses formed in the North Atlantic, datasets from GA02 Leg 2 are used. For NACW, formed through winter subduction at the southern flank of the Gulf Stream42, the closest available Cd values are from GA02 Leg 2 Station 21 (31.7°N, 64.2°W). Formation of NADW is initiated by subduction of surface waters at high latitudes in the North Atlantic Ocean43. However, due to the lack of direct 8114Cd data for the formation region, the closest available Cd data are also from GA02 Leg 2 Station 21, which are used as endmember NADW Cd values. Datasets from the western South Atlantic cruise GA02 Leg 3 are employed to define the Cd endmember values of the South Atlantic-sourced water masses SACW, AAIW, UCDW, and AABW19. South Atlantic Central Water is formed in the region near the South American coast between 30°S and 45°S44, and its endmember composition is based on data for samples from GA02 Leg 3 Station 6 (40.0°S, 42.4°W)19. The formation of AAIW and AABW involves the subduction of HNLC surface waters in the Polar Frontal Zone and near the Antarctic coast, respectively. In contrast, UCDW originates from the subduction and northward advection of Circumpolar Deep Water, which upwells around Antarctica. For AAIW, UCDW, and AABW, the closest available endmember Cd values are from GA02 Leg 3 Station 2 (48.9°S, 48.7°W)B.
The Cd endmember values and the proportions of each water mass present at each location, as derived with the refined eOMP—Monte Carlo approach41 (Supplementary Note 1), were used to predict the Cd concentrations and isotope compositions based on water mass mixing alone for all samples below the surface water. The calculated results from the eOMP-Monte Carlo model, denoted as [Cd]caI(mix) and 8114Cdcal(mix) are compared with the measured data (denoted by the subscript “meas”) in Fig. 3 (see details in Supplementary Table 2). Overall, the model reproduces the observed Cd concentrations and 8114Cd values very well at depths exceeding 1002 m, confirming that conservative mixing dominates the Cd distribution at these depths as shown in Fig. 3. Similar to the results obtained using the original eOMP approach for the full GA02 section3, the modeled [Cd]cal(mix) values are significantly lower than the measured concentrations [Cd]meas for samples collected at depths that record a strong deficiency in dissolved oxygen concentrations, O2-def, with O2- def = [Ozlcal(mix) - [O2]meas. Notably, at these depths the 8114Cdmeas and 8114Cdcal(mix) values are either identical within error or very similar (Fig. 3). The clear discrepancy between calculated Cd concentrations from water mass mixing and measured Cd concentrations only at depths with low dissolved oxygen concentrations, is indicative of remineralization of organic particles, which release dissolved Cd back into the ocean320.
The Cd concentration that is contributed by remineralization, [Cd]remin, can be determined from [Cd]meas and [Cd]eal(mix):
[Cd]remin = [Cd]meas - [Cd]cal(mix) (1)
Calculations that apply Eq. (1) indicate that the concentration of remineralized Cd ranges from 0.03 to 0.31 nmol kg1, contributing 22-86% ([Cd]remin/[Cd]meas) of the total measured Cd inventories for samples where mixing calculations indicate an oxygen deficiency (i.e., high O2-def; Supplementary Table 2). This demonstrates that remineralization is an important source of dissolved Cd in the study area at depths with O2-def values higher than 100 umol kg 1 (Fig. 3). The observation that conservative water mass mixing alone appears sufficient to account for the observed 8114Cdmeas values suggests that the remineralized Cd had a similar isotope composition to that contributed by conservative mixing. This inference is explored below.
As 0";Cd data for biological particles are not available for the GA02 section, the remineralized 8114Cd, denoted as S114Cdremin, is best estimated based on an open-system steady state model for the upper Atlantic Ocean. The data available for North Atlantic thermocline and

Fig. 3 I Depth profiles of Cd concentrations, 8114Cd and dissolved oxygen deficiency O2-def for samples of this study. Top row: Cd concentrations; middle row: 8114Cd, and bottom row: O2-def (see text for definition). Diamonds represent Station 21,squares represent Station 30, triangles represent Station 36, and inverted triangles represent Station 40. The black and blue symbols in the top and middle panels
represent the measured data and the modeled [Cd]cal(mix) and 8114Cdcal(mix) values from eOMP-Monte Carlo model, respectively (see text for details). The gray shaded regions highlight the high oxygen deficiency zone associated with remineralization. The error bars for [Cd] and 8114Cd denote the lsd and 2sd precision, respectively.
Table 1 | Cadmium concentrations and isotope compositions, as well as temperatures, dissolved O2 concentrations and water mass assignments, for the GA02 Leg 2 seawater samples analyzed in this study
Station	Longitude (Έ)	Latitude (°N)	Depth (m)	Temperature (°C)a	Oxygen (umol kg 1)a	[Cd]
(nmol kg 1)	611,Cd (%)	2sdb	nc	Water massd	
21	-64.2	31.7	25	21.7	225	0.003	0.09	0.3	1	SSWe	
21	-64.2	31.7	150	18.3	214	0.012	0.61	0.2	3	NACW	
21	-64.2	31.7	253	18	210	0.02	0.58	0.09	2	NACW	
21	-64.2	31.7	398	17.5	208	0.03	0.57	0.09	2	NACW	
21	-64.2	31.7	790	10	148	0.32	0.46	0.04	8	NACW	
21	-64.2	31.7	1249	4.7	244	0.27	0.45	0.04	4	NADW	
21	-64.2	31.7	1749	3.9	258	0.27	0.44	0.04	4	NADW	
21	-64.2	31.7	2501	3.2	255	0.29	0.42	0.04	4	NADW	
21	-64.2	31.7	3500	2.3	263	0.3	0.39	0.04	4	NADW	
21	-64.2	31.7	4249	2.2	259	0.34	0.36	0.04	5	NADW	
21	-64.2	31.7	4525	2.2	256	0.36	0.34	0.03	5	NADW	
30	-57.6	18.6	20	28.8	197	0.0005	-0.16	0.41	2	SSW	
30	-57.6	18.6	50	27.3	206	0.0007	-0.42	0.37	1	SSW	
30	-57.6	18.6	101	25.5	197	0.0004	0.03	0.3	1	SSW	
30	-57.6	18.6	249	20.9	158	0.01	0.6	0.31	2	NACW	
30	-57.6	18.6	401	14.5	131	0.16	0.56	0.05	8	NACW	
30	-57.6	18.6	749	7.7	129	0.48	0.46	0.05	7	AAIW/UCDW	
30	-57.6	18.6	998	5.8	160	0.46	0.44	0.02	2	NACW	
30	-57.6	18.6	2000	3.5	251	0.28	0.44	0.03	3	NADW	
30	-57.6	18.6	3001	2.7	253	0.31	0.43	0.02	3	NADW	
30	-57.6	18.6	4501	2.1	251	0.38	0.34	0.03	4	AABW	
30	-57.6	18.6	5003	1.9	244	0.46	0.32	0.02	4	AABW	
36	-48.9	7.8	25	28.6	194	0.0015	0.37	0.17	1	SSW	
36	-48.9	7.8	50	28.4	196	0.0017	0.34	0.17	1	SSW	
36	-48.9	7.8	101	20.8	137	0.06	0.63	0.05	3	NACW	
36	-48.9	7.8	201	10.3	121	0.36	0.58	0.1	2	SACW	
36	-48.9	7.8	300	8.7	129	0.44	0.51	0.04	7	SACW	
36	-48.9	7.8	501	6.9	122	0.56	0.5	0.04	6	SACW	
36	-48.9	7.8	1002	4.9	149	0.56	0.44	0.04	7	AAIW/UCDW	
36	-48.9	7.8	1500	4.4	218	0.34	0.44	0.04	4	NADW	
36	-48.9	7.8	2001	3.6	252	0.28	0.42	0.04	4	NADW	
36	-48.9	7.8	3000	2.7	253	0.31	0.38	0.04	5	NADW	
36	-48.9	7.8	3749	2.3	258	0.32	0.37	0.04	4	NADW	
36	-48.9	7.8	4250	2.2	254	0.35	0.35	0.04	5	NADW	
40	-39.7	1.1	27	28	195	0.0013	0.25	0.38	1	SSW	
40	-39.7	1.1	77	27.1	189	0.0019	0.26	0.12	3	SSW	
40	-39.7	1.1	151	15.1	156	0.14	0.58	0.06	4	NACW	
40	-39.7	1.1	251	12.6	146	0.22	0.58	0.06	4	SACW	
40	-39.7	1.1	403	8.5	138	0.44	0.51	0.04	7	AAIW	
40	-39.7	1.1	1002	4.6	153	0.63	0.45	0.07	7	AAIW/UCDW	
40	-39.7	1.1	1501	4.6	189	0.31	0.44	0.07	6	NADW	
40	-39.7	1.1	1999	4	244	0.28	0.42	0.04	5	NADW	
40	-39.7	1.1	3002	3	251	0.31	0.33	0.06	5	NADW	
40	-39.7	1.1	4001	2.5	255	0.35	0.34	0.05	6	NADW	
40	-39.7	1.1	4365	1.4	233	0.6	0.28	0.09	7	AABW	
40	-39.7	1.1	4417	1.2	226	0.62	0.26	0.05	7	AABW	
aThe temperature and oxygen concentration data are from the GEOTRACES Intermediate Data Product (IDP2021v2) 37.
b2 x standard deviation was calculated from η individual sample measurements. If η was <2, the 2sd precision was calculated from repeated bracketing analyses of NIST SRM 3108 Cd solutions in the same measurement session.
cNumber of individual Cd isotope runs; the quoted 6114Cd values are mean results.
The definition of water masses is from the revised eOMP model (Supplementary Note 1), and the assignment represents the water mass with the highest calculated proportion.
eSSW denotes the surface waters.

Fig. 4 I Measured 8114Cdmeas> modeled 8114Cdcal(mix), and calculated
6114Cdcal(mix+remin)plotted against measured [Cd]. Measured 8114Cdmeas data are represented by full colored circles with black borders, modeled 8114Cdcal(mix) values based on water mass mixing from the eOMP-Monte Carlo model are shown as open circles with colored borders, whilst the calculated 8114Cdcal(mix + remin) values from water mass mixing + remineralization (see Eq. 2 in main text) are denoted by crosses with colored borders. The errorbars for 8114Cd are 2sd. The symbol colors represent the dissolved oxygen deficiency, O2-def (see text for definition). Only data for samples from depths between 101 and 1002 m are plotted in the figure. The data sources are provided in Supplementary Table 2. The errorbars for 8114Cd denote the 2sd precision.

Fig. 5 | Plot of Cd isotope compositions versus Cd concentrations with data available for the Atlantic Ocean (2010-2015). Four results for North Atlantic samples with 8114Cd > 2.0%o at low [Cd] are cut off but shown in Supplementary Fig. 4 with an expanded y-axis. Orange, blue and purple symbols represent samples of this study from surface waters (<101 m), thermocline and intermediate water masses (102 - 1002 m), and deep water masses (>1002 m), respectively; these samples were collected at Stations 21 (diamonds), 30 (squares), 36 (triangles), and 40 (inverted triangles). The error bars for 8114Cd denote the 2sd precision. The gray symbols denote literature data for different regions of the Atlantic Ocean: Tropical Atlantic22, North Atlantic8, western South Atlantic19, and eastern South Atlantic17.

isotope compositions. For deep waters, water mass mixing alone is sufficient to explain both Cd concentrations and isotope compositions, in accord with the results of previous studies of the Atlantic and other ocean basins1516”8,45.
Along the GA02 Leg 2 section, the surface waters from depths of about 20 to 100 m have low Cd concentrations of 0.4 χ 10-3 to 3.0 x 10-3 nmol kg1, coupled with relatively low 8114Cd values of —0.42 to 0.37% (n = 8) compared to the deeper waters, which have a 8114Cd of between 0.40 and 0.60% (Fig. 5 and Table 1). Notably, both the Cd concentrations and 8114Cd values of the surface seawater samples are at the lower end of the ranges reported for Atlantic Ocean surface waters (Fig. 5). The surface water 8114Cd values of our samples are, furthermore, at the lower end of the global range, and the lowest values coincide with the lowest Cd concentrations, although the relationship is not statistically significant given the analytical uncertainties8,9,14. Literature observations, in contrast, suggest that preferential uptake of lighter Cd isotopes by phytoplankton at opensystem steady state or closed-system Rayleigh conditions should generate low dissolved Cd concentrations coupled with 8114Cd values that are either nearly unfractionated (0.50-0.60%0) or clearly higher (>l%o), respectively, compared to upwelling deeper waters with typical 8114Cd of 0.40-0.50% (Supplementary Note 2; Supplementary Figs. 3 and 4)91214. In the following, possible origins for the unusual surface water trend of low 8114Cd at low Cd concentrations (Fig. 5) are explored and discussed.
First, it is conceivable that analytical artefacts are responsible for the low observed 8114Cd of the surface waters, given the very low Cd masses (of <1 ng per sample) that were available for the isotopic analyses of these samples. Previous multiple-collector inductively coupled plasma mass spectrometry (MC-ICP-MS) studies that reported high 8114Cd values of up to 5%o for surface waters have been questioned due to potential spectral interferences18,46. In this study, a new chemical procedure effectively removed elements that can cause spectral interference (Ge, Ga, Zr, Mo, In and Sn)47, ensuring unbiased 6114Cd results, even for low Cd concentrations (0.4 x 10“ -3.0 x 10-3 nmolkg”). Our results also align well with

Fig. 6 I Box model for the cycling of Cd in the surface waters of the western tropical North Atlantic. The black arrows represent Cd input fluxes to the surface waters, from atmospheric deposition (Fatm) and upwelling (Fup). The gray arrows represent Cd output fluxes from the surface waters, encompassing downwelling (Fdown) and particle sinking (Fsbp). The values of fatm and fsbp denote the Cd mass fractions provided by atmospheric deposition to the total Cd input flux and by sinking biological particles to the total Cd output flux, respectively. The [Cd] and δ "Cd data denote the Cd concentrations and isotope compositions of the fluxes or reservoirs.
As such, atmospheric emissions enriched in anthropogenic Cd are evaluated in the following as a possible cause of the unusually low 8114Cd values observed in the surface waters of the GEOTRACES GA02 Leg 2 section in the western tropical North Atlantic.
As atmospheric fluxes to the surface ocean will increase dissolved Cd concentrations, the observed surface water trend of low 6114Cd with low Cd concentrations (Fig. 5), cannot be explained by aerosol deposition alone. Instead, it is more likely a consequence of the combined effects of atmospheric deposition of isotopically light Cd, which dominates the 8114Cd fingerprint, and biological Cd uptake, which leads to an overall reduction in dissolved Cd concentrations. A box model is used to quantitatively evaluate the surface water Cd systematics of the western tropical North Atlantic Ocean (Fig. 6), to show that atmospheric Cd deposition is indeed a key factor, as discussed below.
For the surface water box model, the Cd input fluxes are from atmospheric deposition and water upwelling, the Cd output fluxes are ascribed to sinking biological particles and water downwelling (Fig. 6), and it is assumed that horizontal advection is not the major process that controls Cd cycling in the surface waters’2. Considering a steady-state condition for Cd cycling within the box, the mass balance of Cd fluxes and their isotope compositions can be described by the following equations:
Farm + Fw = Fdown + Fsbp (3)
^χδ114^ +Fup x6114Cdul - Fd0„ x6114Cd„+Fibpx6114Cdsbp (4)
where Fatm and Fup represent the Cd input fluxes from atmospheric deposition and upwelling; Fdown and Fsbp are the Cd output fluxes from downwelling and sinking biological particles; 8114Cdatm and 8114Cdtai refer to the Cd isotope compositions of atmospheric deposition and upwelled thermocline waters; and d"Cdssw and 8114Cdsbp denote the Cd isotope compositions of the surface waters and sinking biological particles. The parameters Jatm and / bp are employed below to represent the mass fractions of Cd delivered by atmospheric deposition relative to the total Cd input flux and removed by sinking biological particles relative to the total Cd output flux, respectively:
fatm = Faun/(Faum + Fup)	(5)
sbp Fsbp/(Fdown + Fsbp)	(6)
As the results of the current study suggest that biological uptake of Cd in the Atlantic mainly reflects open-system steady-state conditions (Supplementary Note 2 and Supplementary Fig. 3), the isotope fractionation factor 8114Cdpio-sw can be used to describe the biological fractionation of 8114Cd in the surface waters:
e,,4Cdblo_sw = Abo/sw - 1 = i114Cdsbp - 8114Cd,	(7)
Equations (3) to (7) can furthermore be combined to arrive at the following equation (see detailed derivation in Supplementary Note 3):
811*Cd,e =/am x (5114Cdaun - 811Cd,g) -/Sbp x e114Cdbio_sw + /"'cd,
(8)
The box model, based on Eq. (8), evaluates the impact of atmospheric Cd inputs on 8l4Cdssw. The aim is to determine whether reasonable parameter values, and especially for 6114Cdatrn, can account for the observed surface water Cd isotope compositions, with 8114Cdssw values of —0.42% to 0.37%o (Table 1). Focusing first on the output fluxes, the biological isotope fractionation factor 8 Cdbio-sw has a simple relationship to the fractionation factor Cbio/sw (Eq. 7) . With an Abio/sw value of 0.9998, which is appropriate for many regions of the Atlantic Ocean (Supplementary Note 2 and Supplementary Fig. 3)8,17,19,22, this is equivalent to e114Cdhio-sw = -0.2%. The

Fig. 7 I Modeled relationships between the Cd isotope compositions of the surface waters (6114Cdssw) and atmospheric deposition (6114Cdatm). The panels show results assuming an isotope fractionation factor dbio/sw = 0.9998 (equilavent to s11 Cdbio-sw = —0.2%) for atmospheric deposition with Cd enrichement factors of EF = 1 (a), 50 (b), 100 (c) and 200 (d), corresponding to Cd input flux mass fractions from atmospheric deposition, Jatm, of between 0.06% and 48%. Results are shown for

Cd output flux mass fractions to sinking biological particles, & bp of 0.1 to 0.9 (blue, orange, purple lines). The gray shaded boxes show the range of 8114Cdssw values measured for the surface waters in this study (between -0.42% and 0.37%o). The dashed lines denote the maximum (and for EF = 200 also the minimum) 61HCdatm values of atmospheric deposition that are in accord with the δ' 'Cd,,,,, data obtained for the surface water samples.

value of /sbp is poorly constrained and a wide range of values of between 0.1 and 0.9 are therefore considered.
Considering the Cd input fluxes, a value of 0.60%0 is applied for S114Cdtcl based on the results presented in this study (Fig. 3, Table 1). The fraction of atmospheric Cd deposition fatm, is constrained by the relative Cd fluxes to the surface waters from the atmospheric deposition and water mass upwelling (Eq. 5). As the western tropical North Atlantic is not a major upwelling zone ’4”, an average global upwelling rate of 1 m year 1 is employed, based on a mean oceanic depth of 4000 m and an oceanic overturning time of 4000 years19. Along the GA02 Leg 2 section, the base of the thermocline has a mean [Cd] of ~0.4 nmol kg1; using a seawater density of 1025 kg m\ Fup is hence 46 ug m 2 year 1 (see details in Supplementary Table 3).
The atmospheric Cd flux Fatm was estimated based on the yearly average dust deposition of 2.97 g m2 year 1 in the study area during the GA02 Leg 2 cruise, as derived from the dissolved Al concentration of surface seawater56. For further calculations, the Cd concentration of the aerosol matter is varied from the natural geogenie “background” abundance, with up to 200-fold enrichments of the element, to explore the impact of aerosols that are significantly enriched in anthropogenic Cd. If the Cd in the dust originates only from natural sources, the average upper continental crust [Cd]ucc of 0.09 ugg is a reasonable estimate for the dust, whilst anthropogenic Cd enrichments in the aerosols are quantified by an enrichment factor defined as EF = (Cd/Al)sample/(Cd/Al)ucc57. Whether atmospheric deposition can be responsible for the low observed 8ll4Cdssw values is explored using figures that summarize the modeling results in plots of 8114Cdssw versus 8114Cdatm for different EF values (Fig. 7). In detail, results are shown for EF values of 1, 50, 100, and 200, with Cd solubilities following aerosol deposition of 10% and 80% for Cd of natural and anthropogenic origin, respectively, in accord with published results33-58'59. This yields overall Cd solubilities of 10.0%, 78.6%, 79.3%, and 79.7% for EF values of 1, 50, 100, and 200, respectively. The corresponding Cd deposition rates range from 0.03 to 43 ug m 2 year 1 (Supplementary Table 3), with 8114Cdatm values of between —2.0 and 1.0%13,5160.61.
For EF = 1, corresponding to a Cd deposition rate of 0.03 ug m 2 year 1 and Jatm = 0.06%, atmospheric Cd deposition is unable to account for the low observed surface water 6"4Cdss„ values of the tropical western Atlantic Ocean, for all considered 8114Cdatm values (Fig. 7a). When the aerosol Cd EF is increased to 50, equivalent to a Cd deposition rate of 11 ug m 2 year 1 and Jatm = 19%, 8114Cdatm values of less than —0.76%o can account for 6114Cdssw data that exceed 0.10%0 (Fig. 7b). With an aerosol Cd EF of 100, equivalent to a Cd deposition rate of 21 ug m 2 year 1 and /atm = 32%, 8114Cdatm values of less than—0.18%o are sufficient to explain most of the δ11^^, results (> -0.20%) obtained in this study (Fig. 7c). Finally, for EF = 200, corresponding to a Cd deposition rate of 43 ug m 2 year 1 and Jatm = 48%, 8114Cdatm values that range from 0.09% to -1.58% can account for all 8114Cdssw results of this study (Fig. 7d).
The modeling hence shows that the low observed 6114Cdssw can only be explained if the aerosols are characterized by high EF values for Cd and low 6114Cdatm. The model was also tested with dbio/sw values of 0.9995-0.9997, in accord with previous Cd isotope studies of the Atlantic Ocean81719, and such calculations produced data consistent with those obtained with Abio/ „„-0.9998 (Supplementary Fig. 6). The results of the modeling are discussed in the context of relevant available Cd EF and 8114Cdatm data for marine aerosols, to assess whether atmospheric deposition can be responsible for the low observed d1I4Cdssw values in the western tropical North Atlantic.
Notably, Cd EF data determined for total digested and ammonium acetate leachates of aerosols collected during the GA02 Leg 2 section are fully in accord with the modeled results. In detail, the Cd EFs for the total aerosol digests, also calculated using [Cd]ucc = 0.09 ugg, range from 4 to 370, suggesting strong anthropogenic signatures (Supplementary Table 4)57. The observed variability, with some lower Cd EF values for aerosols collected in our study area, thereby does not question the modeling results. This reflects that surface seawater integrates over the Cd deposition received for several months as a consequence of the Cd residence time in the marine surface layer (0.18-0.35 years; Supplementary Note 4). As such, the 8114Cd of surface seawater is not expected to mirror the EF values of individual aerosol sampling events, but will reflect the cumulative effect of anthropogenically impacted air masses from the same broad source region (given the stable wind pattern of the area; Fig. 1) but with Cd EFs that will inevitably show some temporal variability. Furthermore, the aerosol leachates have high [Cd]/[A1] ratios of 187 to 10,513 ugg1 (Supplementary Table 4) and which are thus approximately 150- to 10,000-fold higher than the upper continental crust [Cd]/[A1] ratio of 1.1 ugg157· The analyses also support the conclusion that anthropogenic Cd is more soluble than mineral-bound Cd, as the aerosol leachates released at least 50%, and in many cases more than 80%, of the total Cd (Supplementary Table 4). The measured solubility is hence very similar to the maximum overall Cd solubility of ~80% assumed for aerosols in the modeling (Supplementary Tables 3 and 4). Considering samples from other Atlantic Ocean cruises, total digested aerosols from the zonal GA03 North Atlantic section showed Cd EFs of approximately 14 within the North African dust plume, but values exceeding 10,000 outside its reach34. Aerosol samples from near Bermuda in the western North Atlantic exhibited Cd EF values of 130-230 due to anthropogenic North American emissions50, whilst analyses of aerosols from the eastern tropical Atlantic (GEOTRACES cruise GA06) showed Cd EFs ranging from 4 to 107, with 61I4Cd values of between -0.54 ± 0.13% and 0.19 ± 0.20%o33. Meanwhile, Pacific Ocean (GEOTRACES cruise GP15) aerosol samples exhibited higher EF values of approximately 200 to 10,000 for Cd, and displayed 8114Cd of -1.91 ± 0.39 to -0.07 ± 0.25%o13. The Cd EF values of 1-200 that are used in our box model calculations are therefore within the range of previously reported marine aerosol compositions.
The Cd present in aerosols is either of natural or anthropogenic origin. As Cd has a low boiling point (of 1040 Κ), it is readily emitted to the atmosphere by high-temperature anthropogenic activities, akin to Pb. Processes such as combustion of fossil fuels, metal smelting, and waste incineration hence contribute to the large anthropogenic emissions of Cd to the atmosphere61·6'. Such anthropogenic emissions, and particularly emissions from waste incineration and lead smelters, exhibit a relatively large range of 8114Cd values, with published data ranging from -0.70 ± 0.13%o to 0.40 ± 0.19%o66,68·69. In contrast, natural Cd in aerosols is expected to be less variable in ó'’Cd, as the upper continental crust has a well-defined mean 8114Cd value of 0.03 ± 0.10%o60. The above findings on the 8114Cd values of natural and anthropogenic aerosols, hence, suggest that the latter commonly feature lower 8114Cd compared to natural dust, but further work is needed to ascertain whether this is indeed a general characteristic of Atlantic Ocean aerosols. Whilst the GA02 aerosol data of this study show no clear correlation between Cd EF and Cd solubility, further studies should also investigate previous findings, which suggest that anthropogenic Cd in atmospheric deposition may be up to 10 times more soluble than natural Cd, implying preferential dissolution of the former in marine surface waters’71"'".
Taken together, the presented evidence hence suggests that atmospheric deposition of Cd from anthropogenic sources likely plays an important role in shaping the Cd concentrations and isotope compositions of the highly Cd-depleted North Atlantic surface waters. The box model developed in this study indicates that atmospheric inputs contributed at least 19% and possibly up to 48% of the surface water Cd inventory, for atmospheric deposition with Cd EF values of 50 to 200 (Fig. 7). Such atmospheric inputs have a clear impact on the Cd cycling in this open ocean setting, with a surface water Cd residence time of less than 0.35 years determined from the results of this study (Supplementary Note 4). However, further analyses of marine aerosols from this and other regions for both Cd concentrations and isotope compositions are needed to better constrain the impact of anthropogenic emissions on the surface water Cd systematics of the western tropical North Atlantic and to investigate whether such processes are also observable in other ocean basins.
In summary, this study provides a comprehensive analysis of variations in Cd concentrations and 6114Cd values determined for seawater samples from four depth profiles in the western tropical North Atlantic. For the deep waters of the study area, the mixing of water masses derived from the Northern and Southern Hemispheres determines the observed Cd concentrations and 8114Cd values. At intermediate depths, the latter parameters are controlled by both internal biogeochemical cycling and ocean circulation, whereby in-situ remineralization accounts for about 20-85% of the dissolved seawater Cd inventory in the low oxygen thermocline and intermediate waters. Viewed in the context of surface water data for the Atlantic and other ocean basins, surface water of the study area is characterized by low Cd concentrations of 0.4 x 10-3 to 3.0 x 10 3nmol kg-1 (n = 8), as typically found in oligotrophic settings, as a consequence of dissolved Cd uptake by phytoplankton. However, these surface waters also have unusually low 8114Cd values of —0.42%o to 0.37%o. A box model approach shows that the low observed 8114Cd values are most likely a consequence of atmospheric Cd inputs that are enriched in isotopically light Cd from anthropogenic emissions, and which overprint any isotope fractionations from biological uptake. Such emissions appear to contribute more than -20% of the surface water Cd inventory in the tropical western Atlantic Ocean. Our results thus highlight that anthropogenic Cd emissions are sufficient to clearly impact Cd cycling in the surface Atlantic Ocean, particularly in oligotrophic open ocean settings with low Cd.
Methods
Sample collection and preparation
The seawater samples were collected from the western tropical North Atlantic at four stations between 13 June and 4 July 2010, during RV Pelagia cruise ΡΕ321 for GEOTRACES section GA02 Leg 2, from Bermuda to Fortaleza (Brazil). Seawater collection utilized an ultraclean all-titanium frame CTD trace metal sampling system with novel PVDF samplers 0. Following collection, the samples were filtered through 0.2 um membranes and acidified on board with distilled HCl to a pH of about 23.
Further sample processing took place in ISO Class 4 laminar flow hoods within an ISO Class 6 clean room in the MAGIC Laboratories at Imperial College London, UK. All water used for processing was of 18.2 M Ω cm grade from a Millipore purification system, whilst mineral acids were purified in-house by sub-boiling distillation in Teflon or quartz glass stills.
The separation of Cd from the seawater matrix followed previously described methods for seawater volumes of >1 L47. Briefly, 1 to 11 L seawater samples were spiked with a nlCd-113Cd double spike solution, with detailed information on seawater sample volumes and the respective ratios of spike- derived to natural Cd provided in Supplementary Table 5. After the samples were left to stand for at least 7 days with daily shaking to ensure full equilibration of the spike-derived and seawater Cd, 1 ml of 0.05 M H2O2 was added to per liter of seawater to break down any remaining organic 
compounds. Following addition of 1 ml of 5 M NH4OAc per liter of seawater, a variable volume of 5 M aqueous NH3 solution was added to adjust the pH to 4.8 ± 0.2, depending on the initial seawater volume and pH. The samples were then equilibrated with Nobias Chelate PA-1 resin which adsorbs Cd and other trace metals from seawater at high yields at these conditions71. The resin was separated from the seawater in a filtration rig, whereupon the adsorbed trace metals, including Cd, were eluted from the resin using the same procedure that is used in the laboratory for the separation of seawater Pb72. The metal fraction was evaporated, refluxed in a mixture of 16 M HNO, and ~10 M H2O2 to break down organic compounds released from the resin, dried again, and then redissolved in 2 M HC1 for further purification with a 2-stage anion exchange column chromatography procedure (Supplementary Table 6). The purified Cd fractions were evaporated to dryness, dried down twice with a few drops of 16 M HNO, and then redissolved in 0.1 Μ HNO3 for the mass spectrometric analyses.
Determination of Cd isotope compositions and concentrations The Cd isotope measurements were performed following established procedures27,47,73 with a Nu Instruments Nu Plasma II MC-ICP-MS instrument, using a Cetac autosampler and either a Nu Instruments DSN- 100 or a Cetac Aridus II desolvation system for sample introduction. Four cadmium isotopes, luCd, 112Cd, 113Cd, and 114Cd, were measured to determine the Cd isotope compositions. To correct for potential isobaric interferences from 112Sn, 113In, and 114Sn, the 115In and 117Sn ion beams were measured during the Cd isotope analysis. All ion beams were monitored with Faraday collectors equipped with 10110 resistors. The sensitivities for Cd in the measurements were 400-900 V ug 1 ml 1 using the DSN-100 and 600-1000 V ug 1 ml-1 with the Aridus II system.
The isotope measurements employed solutions with Cd concentrations of 10-60 ng ml-1, with lower concentrations of 2-5 ng ml1 employed for the surface water samples with low [Cd] of 0.4 x 10-3 to 3.0 χ 10-3 nmolkg-1 (Supplementary Table 5). The sample analyses were carried out interspersed between and relative to bracketing standard solutions that were prepared as mixtures of the NIST SRM 3108 Cd isotope reference material with the Cd double spike, with closely matched ratios of spike-derived to natural Cd (deviations less than 5%) and Cd concentrations (deviations less than 10%). The secondary isotope reference material BAM 1012 Cd, as well as the in-house seawater quality control material SwQC, were analyzed using the same procedures to monitor data quality. Following the measurements, the reduction of the double spike data was carried out offline using an iterative approach, which corrects for both the monitored isobaric interferences and any laboratory- or instrument-induced mass discrimination27,47,73. The obtained Cd isotope compositions are expressed as 8114Cd values relative to NIST SRM 3108 Cd:

8114Cd =

/ (114Cd/110Cd)sample A l(114Cd/110Cd)	1 x10 \ CO/ CC/NIST SRM 3108/
The 8114Cd values obtained in repeated analyses of BAM I012 Cd and the SwQC seawater sample are provided in Supplementary Table 7. The Cd concentrations were determined using the isotope dilution method, from the measured and mass bias-corrected 114Cd/111Cd ratios of the spiked sample solutions. The average Cd yield from the seawater samples was 80%. Total Cd blanks for the complete analytical procedure were 10 to 20 pg. As a result, the blanks typically increased the natural Cd concentrations of the samples by only 0.02-0.07%; the maximum blank contribution, for a surface water sample, was 0.42%. The blank contributions were subtracted from the measured Cd concentrations but not corrected for the Cd isotope data, as the Cd isotope composition of the blank is not well-defined, and estimated blank corrections were smaller than the analytical uncertainties of the data.
Given the challenging nature of the Cd isotope analyses, especially for the low Cd concentration samples from the surface waters, substantial analytical efforts were expended on quality control. These efforts,
documented in detail in Supplementary Note 5 and Supplementary Figs. 7-9, encompass (i) a detailed investigation and evaluation of spectral interferences, (ii) multiple analyses of Cd standard solutions and the SwQC seawater sample using run solutions with Cd concentrations of between 2 ng ml-1 and 60 ng ml-1, and (iii) intercalibration of our [Cd] and 8114Cd depth profiles with published results for samples from the same locations.
Data availability
All source data supporting the findings of this study are provided in the tables in the main manuscript and Supplementary Information. In addition, these data are publicly available via figshare (https://doi.org/10.6084/m9. figshare.30958382).
Code availability
Code and results of the eOMP analysis are currently available upon request.
Received: 1 May 2025; Accepted: 9 January 2026;
Published online: 22 January 2026
References
1.		Boyle, E. A., Sclater, F. & Edmond, J. M. On the marine geochemistry of cadmium. Nature 263, 42-44 (1976).
2.		Bruland, K. W., Knauer, G. A. & Martin, J. H. Cadmium in northeast Pacific waters 1. Limnol. Oceanogr. 23, 618-625 (2003).
3.		Middag, R., van Heuven,S.M.A. C., Bruland, K.W. & de Baar, H. J.W. The relationship between cadmium and phosphate in the Atlantic Ocean unravelled. Earth Planet. Sci. Lett. 492, 79-88 (2018).
4.		Price, N. M. & Morel, F. M. M. Cadmium and cobalt substitution for zinc in a marine diatom. Nature 344, 658-660 (1990).
5.		Cullen, J. T., Lane, T. W., Morel, F. M. M. & Sherrell, R. M. Modulation of cadmium uptake in phytoplankton by seawater CO2 concentration. Nature 402, 165-167 (1999).
6.		Lane, T. W. & Morel, F. M. A biological function for cadmium in marine diatoms. Proc. Natl. Acad. Sci. USA 97, 4627-4631 (2000).
7.		Horner, T. J., Lee, R. B., Henderson, G. M. & Rickaby, R. E. Nonspecific uptake and homeostasis drive the oceanic cadmium cycle. Proc. Natl. Acad. Sci. USA 110, 2500-2505 (2013).
8.		Conway, T. M. & John, S. G. Biogeochemical cycling of cadmium isotopes along a high-resolution section through the North Atlantic Ocean. Geochim.Cosmochim. Acta 148, 269-283 (2015).
9.		Ripperger, S., Rehkàmper, M., Porcelli, D. & Halliday, A. N. Cadmium isotope fractionation in seawater — A signature of biological activity. Earth Planet. Sci. Lett. 261, 670-684 (2007).
10.		John, S. G. & Conway, T. M. A role for scavenging in the marine biogeochemical cycling of zinc and zinc isotopes. Earth Planet. Sci. Lett. 394, 159-167 (2014).
11.		Lacan, F., Francois, R., Ji, Υ. & Sherrell, R. M. Cadmium isotopic composition in the ocean. Geochim. Cosmochim. Acta 70, 5104-5118 (2006).
12.		Abouchami, W. et al. Modulation of the Southern Ocean cadmium isotope signature by ocean circulation and primary productivity. Earth Planet. Sci. Lett. 305, 83-91 (2011).
13.		Sieber, M. et al. Biological, physical, and atmospheric controls on the distribution of cadmium and its isotopes in the Pacific Ocean. Glob. Biogeochem. Cycles 37, e2022GB007441 (2023).
14.		Xue, Z. et al. Cadmium isotope variations in the Southern Ocean. Earth Planet. Sci. Lett. 382, 161-172 (2013).
15.		Zhang, R., Jensen, L. T., Fitzsimmons, J. N., Sherrell, R. M. & John, S. Dissolved cadmium and cadmium stable isotopes in the western Arctic Ocean. Geochim. Cosmochim. Acta 258, 258-273 (2019).
16.		George, E., Stirling, C. H., Gault-Ringold, M., Ellwood, M. J. & Middag, R. Marine biogeochemical cycling of cadmium and cadmium isotopes in the extreme nutrient-depleted subtropical gyre of the South West Pacific Ocean. Earth Planet. Sci. Lett. 514, 84-95 (2019).
17.		Guinoiseau, D. et al. Importance of cadmium sulfides for biogeochemical cycling of Cd and its isotopes in oxygen deficient zones—A case study of the Angola Basin. Glob. Biogeochem. Cycles 33, 1746-1763 (2019).
18.		Janssen, D. J., Abouchami, W., Galer, S. J. G. & Cullen, J. T. Fine-scale spatial and interannual cadmium isotope variability in the subarctic northeast Pacific. Earth Planet. Sci. Lett. 472, 241-252 (2017).
19.		Xie, R. C. et al. Non-Rayleigh control of upper-ocean Cd isotope fractionation in the western South Atlantic. Earth Planet. Scl. Lett. 471, 94-103 (2017).
20.		de Souza, G. F., Vance, D., Sieber, M., Conway, Τ. M.& Little, S. H. Reassessing the influence of particle-hosted sulphide precipitation on the marine cadmium cycle. Geochim. Cosmochim. Acta 322, 274-296 (2022).
21.		Laubach, A. et al. Particulate cadmium accumulation in the Mesopelagic Ocean. Glob. Biogeochem. Cycles 39, e2024GB008281 (2025).
22.		Xie, R. C., Galer, S. J. G., Abouchami, W. & Frank, M. Limited impact of eolian and riverine sources on the biogeochemical cycling of Cd in the tropical Atlantic. Chem. Geol. 511, 371-379 (2019).
23.		Zheng, L., Minami, T., Takano, S., Ho, T. Y. & Sohrin, Y. Sectional distribution patterns of Cd, Ni, Zn, and Cu in the North Pacific Ocean: relationships to nutrients and importance of scavenging. Glob. Biogeochem. Cycles 35, e2020GB006558 (2021).
24.		Martin J-M, Whitfield M. The significance of the river input of chemical elements to the ocean, in: (eds Wong CS, Boyle Ε, Bruland KW, Burton JD, Goldberg ED) Trace Metals in Sea Water (Springer, 1983).
25.		Li, Y.-H. A brief discussion on the mean oceanic residence time of elements. Geochim. Cosmochim. Acta 46, 2671-2675 (1982).
26.		Boyle, E. A. et al. GEOTRACES IC1 (BATS) contamination-prone trace element isotopes Cd, Fe, Pb, Zn, Cu, and Mo intercalibration. Limnol. Oceanogr. Methods 10, 653-665 (2012).
27.		Xue, Z., Rehkamper, M., Schonbachler, M., Statham, P. J. & Coles, B. J. A new methodology for precise cadmium isotope analyses of seawater. Anal. Bioanal. Chem. 402, 883-893 (2012).
28.		Abouchami, W. et al. Biogeochemical cycling of cadmium isotopes in the Southern Ocean along the Zero Meridian. Geochim. Cosmochim. Acta 127, 348-367 (2014).
29.		Bruland, K. W. Oceanographic distributions of cadmium, zinc, nickel, and copper in the North Pacific. Earth Planet. Sci. Lett. 47, 176-198 (1980).
30.		Duce, R. A. et al. The atmospheric input of trace species to the world ocean. Glob. Biogeochem. Cycles 5, 193-259 (2012).
31.		Yang, S.-C., Lee, D.-C. & Ho, Τ.-Υ. The isotopic composition of Cadmium in the water column of the South China Sea. Geochim. Cosmochim. Acta 98, 66-77 (2012).
32.		Hayes, C.T. et al. Replacement times of a spectrum of elements in the North Atlantic based on thorium supply. Glob. Biogeochem. Cycles 32, 1294-1311 (2018).
33.		Bridgestock, L. et al. The Cd isotope composition of atmospheric aerosols from the Tropical Atlantic Ocean. Geophys. Res. Lett. 44, 2932-2940 (2017).
34.		Shelley, R. U., Morton, P. L. & Landing, W. M. Elemental ratios and enrichment factors in aerosols from the US-GEOTRACES North Atlantic transects. Deep Sea Res. Part II Top. Stud. Oceanogr. 116, 262-272 (2015).
35.		Packman, H. et al. Tracing acid mine drainage and estuarine Zn attenuation using Cd and Zn isotopes. Geochim. Cosmochim. Acta 360, 36-56 (2023).
36.		Lambelet, M. et al. Isotopic analysis of Cd in the mixing zone of Siberian rivers with the Arctic Ocean—New constraints on marine Cd cycling and the isotope composition of riverine Cd. Earth Planet. Sci. Lett. 361, 64-73 (2013).
37.		Lambelet, M. et al. Neodymium isotopic composition and concentration in the western North Atlantic Ocean: Results from the GEOTRACES GA02 section. Geochim. Cosmochim. Acta 177, 1-29 (2016).
38.		Middag, R. et al. Dissolved aluminium in the ocean conveyor of the West Atlantic Ocean: Effects of the biological cycle, scavenging, sediment resuspension and hydrography. Mar. Chem. 177, 69-86 (2015).
39.		Rijkenberg, M. J. et al. The distribution of dissolved iron in the West Atlantic Ocean. PLoS ONE 9, e101323 (2014).
40.		Tomczak, M. A multi-parameter extension of temperature/salinity diagram techniques for the analysis of non-isopycnal mixing. Prog. Oceanogr. 10, 147-171 (1981).
41.		Olivelli, A. et al. Vertical transport of anthropogenic lead by reversible scavenging in the South Atlantic Ocean. Earth Planet. Sci. Lett. 646, 118980 (2024).
42.		Klein, B. & Hogg, Ν. On the variability of 18 Degree Water formation as observed from moored instruments at 55°W. Deep Sea Res. Part I Oceanogr. Res. Pap. 43, 1777-1806 (1996).
43.		Stramma, L. et al. Deep water changes at the western boundary of the subpolar North Atlantic during 1996 to 2001. Deep Sea Res. Parti Oceanogr. Res. Pap. 51, 1033-1056 (2004).
44.		Kuhlbrodt, T. et al. On the driving processes of the Atlantic meridional overturning circulation. Rev. Geophy. 45, RG2001 (2007).
45.		Gault-Ringold, M., Adu, T., Stirling, C. H., Frew, R. D. & Hunter, K. A. Anomalous biogeochemical behavior of cadmium in subantarctic surface waters: mechanistic constraints from cadmium isotopes. Earth Planet. Sci. Lett. 341-344, 94-103 (2012).
46.		Sieber, M. et al. Physical and biogeochemical controls on the distribution of dissolved cadmium and its isotopes in the Southwest Pacific Ocean. Chem. Geol. 511, 494-509 (2019).
47.		Huang, Y. et al. Efficient Cd separation protocols for high-precision cadmium isotope analyses of diverse samples by double spike MC- ICP-MS. Talanta 285, 127296 (2025).
48.		Group GIDP. The GEOTRACES Intermediate Data Product 2021v2 (IDP2021v2) (NERC EDS British Oceanographic Data Centre NOC, 2023).
49.		Boyle, E. A., Huested, S. S. & Grant, B. The chemical mass balance of the amazon plume—II. Copper, nickel, and cadmium. Deep Sea Res. Part A Oceanogr. Res. Pap. 29, 1355-1364 (1982).
50.		Fishwick, M. P. et al. Impact of surface ocean conditions and aerosol provenance on the dissolution of aerosol manganese, cobalt, nickel and lead in seawater. Mar. Chem. 198, 28-43 (2018).
51.		Rehkamper, M., Wombacher, F., Horner, T.J. & Xue, Z. Natural and anthropogenic Cd isotope variations. In Handbook of Environmental Isotope Geochemistry 125-154 (Springer, 2012).
52.		Yang, S.-C., Zhang, J., Sohrin, Y. & Ho, Τ.-Υ. Cadmium cycling in the water column of the Kuroshio-Oyashio Extension region: insights from dissolved and particulate isotopic composition. Geochim. Cosmochim. Acta 233, 66-80 (2018).
53.		Hayes, J.M. An introduction to isotopic calculations, 2543 (Woods Hole Oceanographic Institution, 2004).
54.		Liao, F., Liang, X., Li, Y. & Spall, M. Hidden upwelling systems associated with major western boundary currents. J. Geophys. Res. Oceans 127, e2021JC017649 (2022).
55.		McClain, C. R. & Firestone, J. An investigation of Ekman upwelling in the North Atlantic. J. Geophys. Res. Oceans 98, 12327-12339 (2012).
56.		Dulaquais, G. et al. Contrasting biogeochemical cycles of cobalt in the surface western Atlantic Ocean. Glob. Biogeochem. Cycles 28, 1387-1412 (2014).
57.		Rudnick RL, Gao S. Composition of the continental crust, in: Treatise on Geochemistry (Elsevier, 2014).
58.		Hsieh C-C, You C-F, Ho Τ-Υ. The solubility and deposition flux of East Asian aerosol metals in the East China Sea: the effects of aeolian transport processes. Mar. Chem. 253, 104268 (2023).
59.		Annibaldi, A., Truzzi, C., Illuminati, S., Bassotti, E. & Scarponi, G. Determination of water-soluble and insoluble (dilute-HCl-extractable) fractions of Cd, Pb and Cu in Antarctic aerosol by square wave anodic stripping voltammetry: distribution and summer seasonal evolution at Terra Nova Bay (Victoria Land). Anal. Bioanal. Chem. 387, 977-998 (2007).
60.		Pickard, H. et al. The cadmium and zinc isotope compositions of the silicate Earth - Implications for terrestrial volatile accretion. Geochim. Cosmochim. Acta 338, 165-180 (2022).
61.		Shiel, A. E., Weis, D. & Orians, Κ. J. Evaluation of zinc, cadmium and lead isotope fractionation during smelting and refining. Sci. Total Environ. 408, 2357-2368 (2010).
62.		Cheng, K. et al. Atmospheric emission inventory of cadmium from anthropogenic sources. Int. J. Environ. Sci. Technol. 11, 605-616 (2013).
63.		Nriagu, J. 0. & Pacyna, J. M. Quantitative assessment of worldwide contamination of air, water and soils by trace metals. Nature 333, 134-139 (1988).
64.		Pacyna, E. G. et al. Current and future emissions of selected heavy metals to the atmosphere from anthropogenic sources in Europe. Atmos. Environ. 41, 8557-8566 (2007).
65.		Pacyna, J. M. & Pacyna, E. G. An assessment of global and regional emissions of trace metals to the atmosphere from anthropogenic sources worldwide. Environ. Rev. 9, 269-298 (2001).
66.		Cloquet, C., Carignan, J., Libourel, G., Sterckeman, Τ. & Perdrix, E. Tracing source pollution in soils using cadmium and lead isotopes. Environ. Sci. Technol. 40, 2525-2530 (2006).
67.		Gao, B. et al. Precise determination of cadmium and lead isotopic compositions in river sediments. Anal. Chim. Acta 612, 114-120 (2008).
68.		Wang, P. et al. Apportionment of sources of heavy metals to agricultural soils using isotope fingerprints and multivariate statistical analyses. Environ. Pollut. 249, 208-216 (2019).
69.		Wen, H. et al. Tracing sources of pollution in soils from the Jinding Pb-Zn mining district in China using cadmium and lead isotopes. Appl. Geochem. 52, 147-154 (2015).
70.		Rijkenberg, M. J. et al. PRISTINE”, a new high volume sampler for ultraclean sampling of trace metals and isotopes. Mar. Chem. 177, 501-509 (2015).
71.		Biller, D.V. & Bruland, K.W. Analysis of Mn,Fe, Co, Ni, Cu,Zn, Cd, and Pb in seawater using the Nobias-chelate ΡΑ1 resin and magnetic sector inductively coupled plasma mass spectrometry (ICP-MS). Mar. Chem. 130-131, 12-20 (2012).
72.		Griffiths, A. et al. Evaluation of optimized procedures for high- precision lead isotope analyses of seawater by multiple collector inductively coupled plasma mass spectrometry. Anal. Chem. 92, 11232-11241 (2020).
73.		Ripperger, S. & Rehkamper, M. Precise determination of cadmium isotope fractionation in seawater by double spike MC-ICPMS. Geochim. Cosmochim. Acta 71, 631-642 (2007).
74.		Schlitzer, R. Ocean data view, https://odv.awi.de/ (2022).
Acknowledgements
The first author is supported by an Imperial-CSC PhD scholarship. We are grateful to the captain, crew and cruise participants of the ΡΕ321 cruise on RV Pelagia (GEOTRACES section GA02, Leg 2) for the collection of the seawater samples. Steven van Heuven is thanked for sharing his MATLAB code for the eOMP analysis. Greg de Souza and Dylan Rood are thanked for the informal reviews of the manuscript, which led to improvements. The GEOTRACES 2021 Intermediate Data Product version 2 (IDP2021v2) represents an international collaboration and is endorsed by the Scientific Committee on Oceanic Research (SCOR). The many researchers and funding agencies responsible for the collection of data and quality control are thanked for their contributions to the IDP2021v2.
Author contributions
H.X. and M.R. conceptualized and designed the study. H.X. conducted the experiments, performed the data analyses, and prepared the initial manuscript draft. Η.Χ., M.R., Y.H., and Z.J. contributed to the interpretation of the results. Κ.Κ. and B.J. supported laboratory and instrument maintenance. R.E.T.M. facilitated data collection, Α.Ο. contributed to the eOMP-Monte Carlo analysis, and R.M. oversaw sample collection. Α.Β. and R.S. analyzed the aerosols from the GEOTRACES GA02 Leg 2 cruise. H.X., M.R., and T.v.d.F contributed most to later versions of the manuscript, with all authors supplying comments and suggestions.
Competing interests
The authors declare no competing interests.
Additional information
Supplementary information The online version contains supplementary material available at
https://doi.Org/10.1038/s43247-026-03211-w.
Correspondence and requests for materials should be addressed to Hui Xu or Mark Rehkamper.
Peer review information Communications Earth & Environment thanks Ruifang Xie and the other anonymous reviewer(s) for their contribution to the peer review of this work. Primary handling editors: Sophia Johannessen and Joe Aslin. [A peer review file is available].
Reprints and permissions information is available at http://www.nature.com/reprints
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.Org/licenses/by/4.0/.
© The Author(s) 2026
2 SD, n = 5), the difference is minor considering the measurement
Following the addition of remineralized Cd with 8114Cdremin = 0.45%o,
the Cd isotope composition that is expected for a water mass mixture that
was also affected by subsequent remineralization, 8114Cdcal(mix + remind can
be calculated as:
8 Cdeal(ris-tremin) = ([Cdlaitmin) x 814Cdawlkmio + [Cd]emin x 5"4Cdrem„,)/[Cd]meas
(2)
After accounting for remineralization, the calculated 8114Cdcal(mix + remin) Values deviate only slightly from the measured 8114Cdmeas results, typically by less than -0.1l%o (Fig. 4, Supplementary Table 2). This shows that due to the limited extent of biological isotope fractionation, remineralization had only a minor impact on the Cd isotope signatures of the sampled subsurface waters in the western tropical Atlantic. This conclusion is in accord with the results of Guinoiseau et al.1', who used a similar approach to estimate a 8114Cdremin of 0.36%o to 0.56%o for the eastern South Atlantic Ocean, and this altered the 8114Cd signatures of remineralization-affected waters by less than 0.08%0 when remineralization contributed 30% to the Cd inventory.
Overall, the new results confirm earlier findings that remineralization had a major impact on Cd concentrations of intermediate depth waters that record a strong deficiency in dissolved oxygen concentrations3. However, the new 8114Cd data reveal that this process has only a minor influence on 8114Cd, as the dissolved and particulate Cd inventories had similar Cd
previous 8114Cd values determined for surface seawater samples from the Tropical Atlantic22 and eastern South Atlantic Ocean17 (range: -0.61 ± 0.11 to 0.06 ± 0.22%o) by thermal ionization mass spectrometry (TIMS), and from the equatorial Pacific Ocean by MC-ICP-MS13. This indicates that natural processes are likely responsible for 8114Cd values observed in the surface seawater samples of this study.
Internal marine processes, such as Cd uptake by phytoplankton in Zn- depleted environments and Cd complexation by organic ligands, can only explain δ11'Cd values as low as about 0.40%ο for eastern and western South Atlantic Ocean17,19, but not the much lower values observed for the western tropical Atlantic Ocean in this study and for the equatorial Pacific in a previous investigation13. Since internal processes are unlikely to be responsible, external inputs from either rivers or atmospheric deposition are considered in the following.
Despite the proximity of the study area to the Amazon River, our samples were collected outside of the main Amazon plume, where salinity exceeds 35.5, indicating minimal riverine influence (<2%). Moreover, Cd concentrations within the plume are not systematically elevated (Supplementary Fig. 5j", suggesting no clear riverine impact on Cd distribution, consistent with previous studies22,49.
Given the limited riverine contribution, aerosol deposition emerges as the key external source that could affect the surface water 8114Cd values. Although previous studies have suggested that atmospheric deposition of natural Cd is negligible compared to the Cd upwelling fluxes in the Tropical Atlantic Ocean, where strong equatorial upwelling dominates Cd supply22, aerosol-derived inputs may still play a role in regions outside the main upwelling zone. Recent evidence furthermore indicates that anthropogenic Cd emissions can readily overwhelm natural atmospheric Cd inputs, with urban and marine aerosols showing strong enrichments of Cd coupled with δ"'Cd < 0.0%633.34,50,51. Furthermore, Sieber et al.13 proposed that the low 8114Cd (~0.0%o) of oligotrophic equatorial Pacific Ocean surface waters may originate from atmospheric deposition.