# Methylmercury produced in upper oceans accumulates in deep Mariana Trench fauna

Ruoyu Sun, Jingjing Yuan, Jeroen E Sonke, Yanxu Zhang, Tong Zhang, Wang Zheng, Shun Chen, Mei Meng, Jiubin Chen, Yi Liu, et al.

► To cite this version:

Ruoyu Sun, Jingjing Yuan, Jeroen E Sonke, Yanxu Zhang, Tong Zhang, et al. Methylmercury produced in upper oceans accumulates in deep Mariana Trench fauna. Nature Communications, 2020, 11, 10.1038/s41467-020-17045-3. hal-03008002

HAL Id: hal-03008002

https://hal.science/hal-03008002

Submitted on 16 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-020-17045-3

OPEN

<figure><img src="image_3.png" /></figure>

Check for updates

# Methylmercury produced in upper oceans accumulates in deep Mariana Trench fauna

Ruoyu Sun<sup>1</sup>, Jingjing Yuan<sup>1</sup>, Jeroen E. Sonke<sup>2</sup>, Yanxu Zhang<sup>3</sup>, Tong Zhang<sup>4</sup>, Wang Zheng<sup>1</sup>, Shun Chen<sup>5</sup>, Mei Meng<sup>1</sup>, Jiubin Chen<sup>1</sup>, Yi Liu<sup>1✉</sup>, Xiaotong Peng<sup>5✉</sup> & Congqiang Liu<sup>1</sup>

Monomethylmercury (MMHg) is a potent toxin that bioaccumulates and magnifies in marine food webs. Recent studies show abundant methylated Hg in deep oceans (>1000 m), yet its origin remains uncertain. Here we measured Hg isotope compositions in fauna and surface sediments from the Mariana Trench. The trench fauna at 7000-11000 m depth all have substantially positive mass-independent fractionation of odd Hg isotopes (odd-MIF), which can be generated only in the photic zone via MMHg photo-degradation. Given the identical odd-MIF in trench fauna and North Pacific upper ocean (<1000 m) biota MMHg, we suggest that the accumulated Hg in trench fauna originates exclusively from MMHg produced in upper oceans, which penetrates to depth by sorption to sinking particles. Our findings reveal little in-situ MMHg production in deep oceans and imply that anthropogenic Hg released at the Earth's surface is much more pervasive across deep oceans than was previously thought.

<footnote n="1">¹Institute of Surface-Earth System Science, School of Earth System Science, Tianjin University, 300072 Tianjin, China. ²Laboratoire Géosciences Environnement Toulouse, CNRS/Institute for Research and Development/Université Paul Sabatier-Toulouse III, 31400 Toulouse, France. ³School of Atmospheric Sciences, Nanjing University, 210023 Jiangsu, China. ⁴College of Environmental Science and Engineering, Ministry of Education Key Laboratory of Pollution Processes and Environmental Criteria, Tianjin Key Laboratory of Environmental Remediation and Pollution Control, Nankai University, 38 Tongyan Rd., 300350 Tianjin, China. ⁵Deep Sea Science Division, Institute of Deep Sea Science and Engineering, Chinese Academy of Sciences, 572000 Sanya, Hainan, China. ✉email: liuyigeo@tju.edu.cn; xtpeng@idsse.ac.cn</footnote>
Humans are exposed to toxic monomethylmercury (MMHg) by the consumption of marine fish¹. Ocean circulation, biogeochemistry, and marine ecology play a crucial role in MMHg exposure to humans. Natural and anthropogenic releases of Hg to air, land, and water eventually enter the oceans by atmospheric deposition and terrestrial discharge²,³. Within the oceans, inorganic Hg(II) (IHg(II)) is reduced to gaseous Hg(0) through biotic and abiotic processes, but also transformed to its methylated forms, MMHg and dimethyl-Hg (DMHg), potentially by anaerobic bacteria⁴,⁵. Particulate organic matter (POM) scavenges IHg(II) in surface oceans and releases it at depth⁶,⁷. This so-called biological pump results in a macronutrient-like vertical distribution of total Hg concentrations, with low values in epipelagic waters of the surface ocean (0–100 m), increasing through oxygen-depleted mesopelagic waters of the intermediate ocean (100–1000 m), and remaining high in deep ocean waters (>1000 m)²,⁸⁻¹¹.

MMHg makes up a large fraction of methylated Hg in pelagic marine waters, and its levels and chemical speciation determine the burden available to marine food webs¹⁰⁻¹³. Investigating where and how MMHg is produced in oceans is critically important to understand global Hg cycling and human exposure. At present, MMHg is thought to be produced in situ in the oxygen-depleted zone of mesopelagic waters during microbial remineralization of sinking POM derived from phytoplankton primary production in the photic zone⁹,¹²,¹⁴⁻¹⁶. A fraction of mesopelagic MMHg can be transported to surface oceans via diffusion and upwelling, and in situ MMHg production in oxygenated, epipelagic waters has also been demonstrated⁸,¹⁰,¹⁷⁻²⁰. Because MMHg in the epipelagic waters is readily photo-degraded, a depletion of MMHg concentrations is observed in surface oceans²¹,²². Mesopelagic MMHg can penetrate into deep oceans via downwelling or via complexation to or incorporation in sinking POM (i.e., biological pump). However, field data and models indicate that the downward transport of mesopelagic MMHg via downwelling and sinking POM is possibly limited, if not negligible¹²,²³,²⁴. Yet, available data, show abundant MMHg and DMHg in deep oceans⁸,¹²,¹³, implying that in situ Hg methylation and/or di-methylation in deep oceans is plausible. Due to sparse incubation and measurement of deep waters, the origin (surface or deep), lifetime and fate of deep ocean MMHg and DMHg are therefore largely unknown.

Hg stable isotopes exhibit mass-dependent fractionation (MDF, represented by δ²⁰²Hg) in all investigated abiotic and biotic Hg transformations, and large mass-independent fractionation of odd-mass number isotopes (odd-MIF, represented by Δ¹⁹⁹Hg and Δ²⁰¹Hg) in abiotic, photochemical Hg transformations²⁵,²⁶. Significant even-mass number MIF (even-MIF, represented by Δ²⁰⁰Hg and Δ²⁰⁴Hg) has also been observed primarily in atmospheric precipitations, which is proposed to be related to photochemical oxidation of gaseous Hg(0) in the upper atmosphere²⁷,²⁸. Due to lack of significant Hg isotope fractionation during accumulation and trophic transfer of MMHg in aquatic food chains²⁹,³⁰, Hg isotope compositions of marine biota have demonstrated great potential in tracing the sources, production and degradation of MMHg in oceans¹⁶,³¹,³².

The hadal zone of oceans hosts a great diversity of endemic fauna, potentially serving as useful vectors to constrain MMHg dynamics and sources in deep oceans. Here, we measured concentrations and isotope compositions of Hg in endemic amphipods (mostly *Hirondellea gigas*) at water depths of 7000–11,000 m below the surface (mbs) and in benthic sediments at water depths of 5500–9200 mbs from the oligotrophic Mariana Trench in the North Pacific Ocean (NPO, 11.5°N, 142.5°E). Two amphipod samples (*Hirondellea gigas*) and one snailfish sample at water depths of ~8000 mbs from the nearby Yap Trench (9.5°N, 138.5°E)

were also collected and measured (Fig. 1 and Supplementary Fig. 1). We find that both MDF and MIF of trench fauna are comparable to those of upper NPO biota MMHg, suggesting that the accumulated Hg in trench fauna originates from MMHg produced in upper oceans.

## Results

**Total Hg and MMHg concentrations.** Total Hg concentrations in trench amphipods have an average of 547 ± 230 ng g⁻¹ (1 SD, 235 to 1070 ng g⁻¹, dry weight, Supplementary Table 1), which is comparable to their counterparts from the abyssal Arctic Ocean (250–1000 ng g⁻¹)³³, but considerably higher than benthic organisms in freshwater and coastal ecosystems (mostly <400 ng g⁻¹ on average)³⁴,³⁵. The relatively high Hg concentrations in trench amphipods likely relate to their low tissue turnover rates and long life spans³⁶,³⁷. Total Hg concentrations of amphipods show no correlation with their depths of occurrence (R² = 0.00, *P* = 0.85), and an insignificantly positive correlation with their body length, a surrogate of life span (R² = 0.12, *P* = 0.08, Supplementary Fig. 2a). MMHg concentrations were measured in nearly half of amphipods, and vary from 8 to 403 ng g⁻¹. The fraction of MMHg (MMHg%) is highly variable, from 2 to 59% of total Hg. A significantly negative correlation is observed between MMHg% and body length of amphipods (R² = 0.48, *P* = 0.02, Supplementary Fig. 2b). The snailfish that feeds on amphipods has a total Hg concentration of 970 ng g⁻¹ and MMHg concentration of 809 ng g⁻¹, which indicates that its MMHg% (83%) is higher than the corresponding values of amphipods (Supplementary Tables 1 and 2).

**Hg isotope compositions.** The most striking observation of Hg isotopes is that both MDF (δ²⁰²Hg: −0.05 to 0.54‰, *n* = 28) and odd-MIF (Δ¹⁹⁹Hg: 1.26 to 1.70‰, Δ²⁰¹Hg: 1.01 to 1.37‰, *n* = 28) of trench amphipods are comparable to those of NPO fishes at 300–600 mbs¹⁶,³², despite their dramatic difference in depths (Fig. 2). The variability of MDF and odd-MIF for amphipods is only ~0.5‰, and their variations are not correlated with MMHg % (2 to 59%, Fig. 3), depths of occurrence (7000–11,000 mbs), sampling locations (Yap and Mariana trenches) or sampling time (July 2016 to March 2017) (Supplementary Fig. 3 and Supplementary Table 2). However, positive correlations are seen between δ²⁰²Hg and Hg concentrations (235–1070 ng g⁻¹) (R² = 0.27, *P* = 0.01), and between Δ¹⁹⁹Hg and body length (1–4.5 cm) (R² = 0.18, *P* = 0.04) (Supplementary Fig. 3). Although these positive correlations suggest that physiological processes in amphipods may affect their Hg isotope compositions, they can only explain 20–30% of Hg isotope variations. Hg isotope compositions of the snailfish are within those of amphipods (Fig. 2). In contrast, the trench sediments have negative δ²⁰²Hg (−0.96 ± 0.27‰, 1 SD, *n* = 5) and very small, positive odd-MIF (Δ¹⁹⁹Hg = 0.20 ± 0.07‰, Δ²⁰¹Hg = 0.18 ± 0.04‰, 1 SD, *n* = 5) (Supplementary Table 3). The even-MIF of trench fauna (amphipod and snailfish) is characterized by positive Δ²⁰⁰Hg (0.02–0.13‰) and negative Δ²⁰⁴Hg (−0.18 to −0.03‰), which is very similar to those of trench sediments, and upper marine particles and fishes from the upper NPO¹⁶,³² (Supplementary Fig. 4).

## Discussion

The similar MIF signatures in deep fauna and surface fishes reveal important insight into hadal Hg provenance. First, Δ²⁰⁰Hg has thus far been used as a conservative tracer at the Earth's surface, and is thought to be generated during net Hg(0) oxidation in the upper atmosphere²⁸,³⁸. The identical, positive Δ²⁰⁰Hg in deep and surface fauna (Supplementary Fig. 4) suggests that the IHg (II) precursor to MMHg in the deep and surface fauna is similar
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Location map showing the samples from Mariana/Yap trenches. a Overview of sampling locations, with the inset **b** for Yap Trench and inset **c** for Mariana Trench. Sediment (green squares), amphipod (orange circles), and snailfish (blue triangle) samples are marked with filled symbols.</figcaption></figure>

and likely of atmospheric Hg origin. Second, odd-MIF mainly accompanies photoreduction of IHg(II) and photodegradation of MMHg in aqueous environments³⁹,⁴⁰. The Δ¹⁹⁹Hg/Δ²⁰¹Hg ratio is a diagnostic of photochemical reaction types, with an experimental value of ~1.3 for MMHg photodegradation and ~1.0 for IHg(II) photoreduction³⁹. The deep fauna all have positive odd-MIF values, and Δ¹⁹⁹Hg and Δ²⁰¹Hg are highly correlated along a linear slope of 1.15 ± 0.08 (1SE, Fig. 4a). This Δ¹⁹⁹Hg/Δ²⁰¹Hg slope is comparable to slopes (~1.2) reported for fishes from upper NPO and other coastal and open marine waters¹⁶,³¹,³²,⁴¹, suggesting that the odd-MIF in the trench fauna is produced by similar surface ocean MMHg photodegradation processes. Experimental photodegradation of aqueous MMHg has also been suggested to induce a characteristic Δ¹⁹⁹Hg/δ²⁰²Hg ratio of 2.43 ± 0.10³⁹. A similar slope is observed in NPO fishes (2.13 ± 0.19, 1SE, Fig. 4b)¹⁶,³². Most of our trench fauna samples lie on the NPO fish slope, but do not show sufficient variation in both

Δ¹⁹⁹Hg and δ²⁰²Hg to confirm the Δ¹⁹⁹Hg/δ²⁰²Hg slope as a diagnostic. Minor δ²⁰²Hg deviations in some fauna samples (Fig. 4b) are likely caused by small amounts of non-chemical, microbial degradation of MMHg, which results in no MIF but increases the δ²⁰²Hg of residual MMHg accumulated in fauna⁴².

The isotope compositions of total Hg in trench amphipods are representative of the isotope signatures of MMHg in the same amphipods, despite the relatively low MMHg% of 2–59%. The reason for this is that most of the IHg(II) in the amphipods is likely produced by in vivo degradation of MMHg rather than assimilation of IHg(II) from the ambient environments. This is supported by multiple lines of evidence. (1) The significantly negative correlation between MMHg% and body length of amphipods (Supplementary Fig. 2b) indicates that in vivo demethylation of MMHg likely occurs in the amphipods. MMHg demethylation is commonly observed in organisms with a
<figure><img src="image_5.png" /></figure>

**Fig. 2** Depth profiles of δ<sup>202</sup>Hg and Δ<sup>199</sup>Hg in trench samples. Samples measured in this study are marked with filled symbols. Δ<sup>201</sup>Hg (not shown) shows similar vertical trend as Δ<sup>199</sup>Hg. Also shown are the marine particles and fishes (open symbols) from North Pacific Subtropical Gyre<sup>16,32</sup>, and the best-fitted regression curves of water depths with δ<sup>202</sup>Hg (R<sup>2</sup> = 0.812, *p* = 1.11E-16, *n* = 37 samples) and Δ<sup>199</sup>Hg (R<sup>2</sup> = 0.749, *p* = 1.32E-14, *n* = 37 samples) for fishes. The dashed lines show the water depths predicted by the regression curves using the δ<sup>202</sup>Hg and Δ<sup>199</sup>Hg values of trench fauna. The typical analytic uncertainty of our fauna (sediment) samples was 0.08‰ (0.08‰, two standard deviation) for δ<sup>202</sup>Hg and 0.10‰ (0.04‰, two standard deviation) for Δ<sup>199</sup>Hg.

gastro-intestinal tract including isopods, fish and mammals<sup>43-45</sup>. (2) The trench fauna MIF values exhibit a Δ<sup>199</sup>Hg/Δ<sup>201</sup>Hg regression slope of 1.15 that is diagnostic of MMHg photo-degradation rather than IHg(II) photoreduction. This suggests that most of the IHg(II) fraction in amphipods inherits its Δ<sup>199</sup>Hg and Δ<sup>201</sup>Hg signatures from ingested MMHg, following in vivo demethylation. (3) Trench amphipods have very consistent, high δ<sup>202</sup>Hg (0.27 ± 0.14‰, 1 SD) and Δ<sup>199</sup>Hg (1.47 ± 0.13‰, 1 SD), and exhibit little variation in both signatures in spite of large variations in MMHg% (Fig. 3). The snailfish with 83% of MMHg % has nearly the same Hg isotope composition as trench amphipods which have varying and low MMHg%. If significant fraction of IHg(II) was to be assimilated from the surrounding environments, we would expect a more scattered Hg isotope distribution due to the mixing of IHg(II) and MMHg. We would expect that the trench fauna with high IHg(II)% to have low δ<sup>202</sup>Hg and Δ<sup>199</sup>Hg, because IHg(II) in the oceans has relatively low δ<sup>202</sup>Hg and Δ<sup>199</sup>Hg, as seen from our measured seafloor sediments and previously reported open ocean particles<sup>32</sup> and sediments<sup>46,47</sup>, and coastal seawater<sup>48,49</sup>. (4) Insignificant differences (<0.1‰) in Hg isotope compositions are observed between the whole tissues and the separated tissues (muscle, lipid, and gut contents) of amphipods, although their Hg concentrations and speciation vary significantly (Supplementary Table 2). This suggests that the internal distribution or in vivo transformation of Hg species do not significantly fractionate Hg isotopes between MMHg and IHg(II) in amphipods. In the following we discuss the origin of fauna MMHg, based on all available open ocean Hg isotope data, and make a number of deductions that constrain the deep ocean MMHg cycling.

IHg(II) in sinking marine particles of different sizes and depths in the NPO has a constant Δ¹⁹⁹Hg of 0.16 ± 0.09‰ (1 SD)³², which is found back at depth in Mariana sediments (0.20 ± 0.07‰, 1 SD) (Fig. 2). Potential MMHg produced in situ in trench sediments or deep water column should therefore inherit

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 3** Trench fauna δ<sup>202</sup>Hg, Δ<sup>199</sup>Hg and monomethylmercury fraction.</figcaption></figure>

Linear correlations of a δ<sup>202</sup>Hg versus monomethylmercury fraction (MMHg%), and b of Δ<sup>199</sup>Hg versus monomethylmercury fraction (MMHg%). The typical analytic uncertainty of our fauna samples was 0.08‰ (two standard deviation) for δ<sup>202</sup>Hg and 0.10‰ (two standard deviation) for Δ<sup>199</sup>Hg. *n* is the number of samples.

the low Δ¹⁹⁹Hg of ~0.2‰. MMHg produced from hydrothermal vent fluids is expected to have near-zero Δ¹⁹⁹Hg, because Hg of geological origin generally has near-zero MIF⁵⁰. The trench amphipods are benthic scavengers, which mainly reside in the top of the sediments surrounding the hydrothermal vent fields, and feed on sinking particles. If MMHg is produced in the sediments, diffusing to the water column, or in hydrothermal fluids or in the deep water column, it would transfer the near-zero to very small Δ¹⁹⁹Hg to the amphipods. The occurrence of significant, 1.47 ± 0.13‰ (1 SD), Δ¹⁹⁹Hg in the deepest marine fauna on Earth suggests that this Δ¹⁹⁹Hg is not produced in situ in the deep ocean environments. Instead, elevated MMHg Δ¹⁹⁹Hg can only be acquired during photodegradation of MMHg in the euphotic zone (Fig. 2), since potential microbial methylation of IHg(II) and non-photodegradation of MMHg in deep waters would not induce MIF⁴²,⁵¹⁻⁵³. We therefore suggest that the trench MMHg is produced in the surface ocean and subsequently transport to depth before entering the deep ocean food web.

IHg(II) and MMHg can be transported from surface to inter-
mediate and even deep oceans by sinking particles, and by sub-
ducting water masses during deep water formation⁷,⁵⁴. The latter
transport pathway is unlikely in the Mariana/Yap trench system
because there is little deep water formation in the NPO⁵⁵. The
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 4 Correlations of mercury isotope values in trench samples. a $\Delta^{201}\text{Hg}$ versus $\Delta^{199}\text{Hg}$, and b $\delta^{202}\text{Hg}$ versus $\Delta^{199}\text{Hg}$. Also shown are the marine particles and fishes from North Pacific Subtropical Gyre$^{16,32}$, and the linear regression lines of Hg isotope compositions for fishes$^{16,32}$. The inset in a shows the linear regression line of $\Delta^{201}\text{Hg}$ versus $\Delta^{199}\text{Hg}$ in trench fauna. Symbols for samples are the same as in Fig. 2. SE is the standard error, and *n* is the number of samples. The arrows show the experimental isotope fractionation trajectories during microbial methylation of IHg(II)$^{51–53}$, and microbial and photodegradation of MMHg$^{39,42}$. The typical 2 SD ana- lytic uncertainty of our fauna (sediment) samples was 0.08% (*c*$_{e}$, two standard deviation) for $\delta^{202}\text{Hg}$, and 0.10% (*c*$_{e}$, two standard deviation) for $\Delta^{199}\text{Hg}$ and $\Delta^{201}\text{Hg}$.</figcaption></figure>

lifetime of MMHg with a surface ocean $\Delta^{199}\text{Hg}$ of ~1.5‰ would have to be thousands of years for the MMHg to arrive in the Mariana Trench by diffusion and advection, which is also unli- kely⁴. The mixing of MMHg produced in epipelagic and meso- pelagic waters has been previously proposed to account for the observed declines of $\delta^{202}\text{Hg}$ and $\Delta^{199}\text{Hg}$ values in NPO fishes with their feeding depths (Fig. 2)¹⁶. We observe that $\Delta^{199}\text{Hg}$ values in trench fauna at 7000–11000 mbs are comparable to those of NPO fishes at 300–600 mbs (Fig. 2). We use mean NPO surface and intermediate dissolved MMHg concentrations to estimate a weighted average $\Delta^{199}\text{Hg}$ of 1.44 ± 0.75‰ (1 SD) in surface/intermediate food web MMHg (see Methods). The close correspondence of the weighted mean $\Delta^{199}\text{Hg}$ of NPO surface/ intermediate MMHg to trench fauna $\Delta^{199}\text{Hg}$ (1.47 ± 0.13‰, 1 SD) suggests that MMHg is transported to depth on sinking

particles. The importance of particles to move IHg(II) from the upper ocean to greater depths have been recently highlighted⁶,⁷. A study on amphipods from the Mariana Trench and two other Pacific trenches revealed that the bomb $^{14}$C signals of the surface ocean could rapidly penetrate into the deepest ocean via fast sinking of surface POM⁵⁶. It is thus likely that, similar to IHg(II), the MMHg produced in the surface and intermediate waters above the trenches transfers to the deep ocean through sinking of surface POM as well. A conceptual model of Hg cycling in the Mariana/Yap trench system is shown in Fig. 5, illustrating the methylation and degradation of Hg, and how the surface MMHg odd-MIF propagates to deep food webs.

The $\delta^{202}\text{Hg}$ values of trench fauna at 7000–11000 mbs (0.27 ± 0.14‰, 1 SD) are ~0.3‰ higher than the weighted mean $\delta^{202}\text{Hg}$ (−0.05 ± 0.30‰, 1 SD) of NPO surface/intermediate MMHg (see Methods). This magnitude of the $\delta^{202}\text{Hg}$ shift is comparable to the ~0.4‰ shift observed during microbial degradation of MMHg in culture, which increases the $\delta^{202}\text{Hg}$ value of residual MMHg⁴². Thus, microbial degradation of MMHg likely occur during MMHg transport from the upper ocean to depth. The ensuing question is why MMHg is degraded, but not produced in the deep ocean that occupies more than 70% of the whole Ocean’s volume. Both non-photochemical degradation and production of MMHg are predominantly mediated by microbes. The microbes that produce MMHg are thought to be mostly anaerobic bacteria with *hgcAB* genes, while the microbes that degrade MMHg could be either aerobic or anaerobic bacteria⁵⁷,⁵⁸. The microbes that are capable of MMHg degradation are thus more widespread in nature, including oceans. The deep marine waters are cold, nutrient-limited, and oxygen-rich⁵⁹ (Fig. 5), which might inhibit the methylation ability of anaerobic bacteria⁶⁰. The bacteria in deep oceans that are capable to demethylate MMHg also likely exist in the guts of trench amphipods⁴⁴,⁴⁵, and cause the in vivo demethylation of MMHg we infer above.

The Hg isotope compositions of trench fauna clearly indicate that bioavailable MMHg in the planet’s deepest ocean is produced in the upper ocean, and that consequently little MMHg is pro- duced in the deep ocean. We use a binary mixing model based on $\Delta^{199}\text{Hg}$ to estimate that 37–48% and 52–63% of MMHg in the trench fauna is derived from the surface and the intermediate ocean, respectively (see Methods). Several centuries of human activities have severely perturbed natural Hg cycling, elevating atmospheric Hg levels at the Earth’s surface by a factor of 3–7$^{2,61}$. Currently, ~300 Mmol of anthropogenic Hg is estimated to reside in oceans, with two thirds in the upper ocean and the remainder subducting in North Atlantic Deep Water and Antarctic Bottom Water². The tight MMHg linkage between trench fauna and the upper ocean suggest that anthropogenic Hg is much more per- vasive across the world’s oceans than was previously thought, and may reach the deepest ocean by fast particle sinking. Ongoing surface ocean warming was recently suggested to increase MMHg levels in the upper ocean⁶², which would likely rapidly propagate to the hadal ecosystems.

Our findings not only inform on the source of MMHg in the deep ocean, but also allow us to speculate about the production of DMHg in the deep ocean. Mean oceanic dissolved MMHg and DMHg profiles show similar concentration trends at depth throughout the world’s oceans$^{10–13,24}$. The limited global varia- tion in MMHg and DMHg concentrations indicates slow trans- formations, supported by our finding that little IHg(II) methylation takes place in the deep NPO. The exact link between MMHg and DMHg in marine waters is presently unclear. State of the art marine Hg models exclude direct DMHg formation from IHg(II)$^{24,63}$, supported by slow observed rates¹⁷, but include microbial methylation of IHg(II) to form MMHg, and subsequent microbial methylation of MMHg to form DMHg in deep waters.
<figure><img src="image_10.png" /><figcaption>**Fig. 5** Conceptual illustration of mercury cycling in trench systems. Atmospheric inorganic Hg<sup>II</sup> (IHg(II)) is deposited to the epipelagic waters where it is microbially (1) or photochemically (2) reduced to Hg(0), microbially methylated to monomethylmercury (MMHg) (3), and scavenged as Hg<sup>P</sup> by sinking particles (4). Photo-reduction of IHg(II) results in residual IHg(II) with small odd-MIF (2), while photodegradation of MMHg results in residual MMHg with large odd-MIF (5). The sinking particles are decomposed and release IHg(II) in the mesopelagic waters (6) where microbial methylation produces MMHg that inherits the small odd-MIF of IHg(II) (7). The mixing of MMHg produced in epipelagic and mesopelagic waters in the intermediate ocean (8) results in MMHg with moderate odd-MIF<sup>16</sup>, which is then transported to the deep ocean by sinking particles (9). Since no/little MMHg is produced in the deep ocean, the MMHg with moderate odd-MIF is incorporated into deep food webs (10). Also shown are typical depth profiles of NPO dissolved MMHg concentrations<sup>10,12</sup>, Mariana dissolved oxygen [O<sub>2</sub>]<sup>59</sup>, and Mariana seawater salinity and temperature measured in this study. The hadal food web is based on Gerringer<sup>37</sup>.</figcaption></figure>

Here, we find that deep ocean MMHg is supplied from the upper ocean by particle sinking rather than in situ microbial methylation. DMHg is a dissolved gas and is unlikely to be transported on particles. If in situ microbial MMHg production in the deep waters is limited, then microbial methylation of MMHg to form DMHg is likely also limited. This leads us to speculate that deep ocean DMHg is abiotically produced, possibly on particle surfaces as supported by recent laboratory experiments<sup>64</sup>.

**Methods**

**Sampling.** Endemic amphipods were captured by bait traps installed on the deep-sea lander vehicles (@Tianya, @Yuanwei, @Wanquan) that were operated by the Chinese Academy of Sciences during several cruises between July 2016 and March 2017. The lander vehicles were deployed to the seafloor of Mariana (11.5°N, 142.5°E) and Yap trenches (9.5°N, 138.5°E) in the NPO at depths of 7000–11,000 mbs for 12–24 h and monitored by cameras (Fig. 1 and Supplementary Table 1). The bait was enclosed in an isolated mesh bag of a funnel trap, which allowed for odor plume diffusion and prevented consumption of the bait. The captured amphipods were frozen immediately at –80 °C upon loading on deck. The amphipods mainly comprised of *Hirondellea gigas* with body lengths of 1–4.5 cm except one amphipod *Alicella gigantea* with a body length of 12 cm (Supplementary Fig. 1). After returning to laboratory, amphipods were rinsed with Milli-Q water before tissue separation, freeze-drying and homogenization. The amphipods mostly had empty gut contents, and were processed as whole tissues after removing the heads to avoid potential contamination by bait. To increase the sampling representativeness, several amphipod individuals of similar sizes were usually combined to one sample. Three large-size amphipod samples were separated for muscle, lipid and gut contents using pre-cleaned surgical scissors and tweezers to assess whether internal distribution or in vivo transformation of Hg species induce Hg isotope fractionation (Supplementary Table 2). We also trapped one snailfish that primarily fed on amphipods in the Yap Trench, for which we separated its muscle and processed it in the same way as amphipods. Ocean floor sediments at depths of 5500–9200 mbs were collected with a box corer that was launched by a geological winch during one cruise operated by the Chinese Academy of Sciences in 2016. Immediately after recovery, short sediment cores were collected by pushing PVC tubes (L = 80 cm, ID = 9 cm) into the sediment blocks. The top 0–6 cm layers of sediments from all the short cores were then subsampled and stored at –20 °C before freeze-drying and homogenization.

**Total Hg and MMHg concentration measurement.** Total Hg concentrations in the homogenized freeze-dried samples were measured by a Lumex RA-915F Mercury Analyzer. The samples were combusted in the PYRO-915+ unit, and the evaporated Hg(0) was then purified and transferred to an analytical cell before measurement by differential Zeeman atomic absorption spectrometry (AAS). Total Hg concentrations in the digested sample solutions used for Hg isotope ratio analysis were measured by a Tekran 2600 cold vapor atomic fluorescence spectrometer (CV-AFS) according to US-EPA method 1631E. MMHg was determined in amphipod and snailfish samples by a Tekran 2700 gas chromatograph-CV-AFS after solvent extraction using KOH/CH<sub>3</sub>OH solution, and ethylation by NaBEt<sub>4</sub> in closed purge vessels<sup>35</sup>. The combined 2 SD analytic uncertainties of total Hg concentrations, as evaluated by periodically measured certified reference materials (CRMs, GBW07405, soil; DORM-4, dogfish muscle), were <5% for both Lumex RA-915F and Tekran 2600. The combined 2 SD analytic uncertainties for MMHg concentrations of DORM-4, were <15%. The relative standard deviations of sample/CRM duplicates were <5% for THg and <8% for MMHg concentrations.

**Hg isotope ratio measurement.** According to the measured total Hg concentrations in freeze-dried samples, ~0.2–0.6 g of each fauna sample was digested using 4 mL double-distilled HNO<sub>3</sub> and 4 mL super-pure H<sub>2</sub>O<sub>2</sub> in a 30 mL Teflon lined vessel heated by a programmed microwave digestion system. Sediments were pre-concentrated for Hg using the combustion-trapping method<sup>65</sup>, and 10 mL of 40% double-distilled acid (2HNO<sub>3</sub>/1HCl, v/v) was used to trap the volatilized Hg (0). Procedural blanks and CRMs were processed with the samples in a same manner. The procedural blanks accounted for <1% of Hg mass in the samples, and the Hg digestion and preconcentration recoveries were in the range of 88–110% for both the samples and procedural CRMs. The processed sample solutions were diluted with Milli-Q water to Hg concentrations of 0.5–2 ng g<sup>-1</sup>, and were measured for Hg isotope ratios by coupling a customized cold vapor generation system to multi-collector inductively coupled plasma mass spectrometry (MC-ICPMS, Nu Plasma 3D at Tianjin University, China). The typical sensitivity for <sup>202</sup>Hg was ~2 V per ng g<sup>-1</sup> Hg at a solution uptake rate of 0.8 mL min<sup>-1</sup>. Instrumental mass bias was corrected by both an internal NIST 997 Tl standard solution (supplied via Aridus II desolvation nebulizer system) using the exponential fractionation law and NIST 3133 Hg standard-sample bracketing method. The bracketed NIST 3133 solutions were matched to the sample solutions within 5% in both acid matrix and Hg concentrations. The Faraday cups were positioned to simultaneously collect all seven Hg isotopes and two Tl isotopes. Acquisition time was 7 min (5 blocks, 20 cycles, 4.2 s of integration time) with 3 min of initial uptake time. Between samples,
the system was washed with the sample matrix solution for 7 min to ensure that the blank signals were <1% of the preceding sample or standard signals.

Hg isotope ratio is expressed as δ$^{xxx}$Hg (‰, xxx = 199, 200, 201, 202, 204) by normalizing to a common NIST 3133 Hg standard:

$$ \delta^{xxx}\text{Hg}(\text{‰}) = \left[ \frac{(^{xxx}\text{Hg}/^{198}\text{Hg})_{\text{sample}}}{(^{xxx}\text{Hg}/^{198}\text{Hg})_{\text{NIST 3133}}} - 1 \right] \times 1000 \quad (1) $$

MIF value is denoted as A$^{xxx}$Hg (‰, xxx = 199, 200, 201, 204), representing the difference between the measured δ$^{xxx}$Hg value and that predicted from δ$^{202}$Hg using a kinetic MDF law$^{66}$:

$$ \Delta^{xxx}\text{Hg}(\text{‰}) = \delta^{xxx}\text{Hg} - \Delta^{xxx}\beta \times \delta^{202}\text{Hg} \quad (2) $$

The mass-dependent scaling factor xxxβ is 0.2520 for $^{199}$Hg, 0.5024 for $^{200}$Hg, 0.7520 for $^{201}$Hg, and 1.4930 for $^{204}$Hg.

Hg isotope ratios of secondary standard NIST 3177 solution and procedural CRMs (GBW07310, DORM-4) analyzed during different analytic sessions (Supplementary Table 4) were in agreement with those reported in previous studies$^{35,65,66}$. The typical 2 SD analytic uncertainties of samples were estimated as the larger 2 SD uncertainties of Hg isotope ratios in NIST 3177 and GBW07310 (sediment) or DORM-4 (amphipod and snaifish). The 2 SD uncertainties of Hg isotope ratios in samples with replicate analyses were applied as the analytic uncertainties only when they were larger than the typical 2 SD analytic uncertainties.

**Estimate of weighted MMHg isotope values of upper NPO.** The weighted MMHg isotope values (δ$^{202}$Hgupper, Δ$^{199}$Hgupper, Δ$^{201}$Hgupper) of the upper NPO (0–1000 mbs) are estimated by the following binary mixing equations using a Monte Carlo simulation approach (n = 10,000 times) through the pseudorandom number generation function of the MatLab software (R2016b, MathWorks):

$$ \delta^{202}\text{Hgupper} = (C_{\text{surf}} \times \delta^{202}\text{Hg}_{\text{surf}} + C_{\text{int}} \times \delta^{202}\text{Hg}_{\text{int}})/(C_{\text{surf}} + C_{\text{int}}) \quad (3) $$

$$ \Delta^{199}\text{Hgupper} = (C_{\text{surf}} \times \Delta^{199}\text{Hg}_{\text{surf}} + C_{\text{int}} \times \Delta^{199}\text{Hg}_{\text{int}})/(C_{\text{surf}} + C_{\text{int}}) \quad (4) $$

$$ \Delta^{201}\text{Hgupper} = (C_{\text{surf}} \times \Delta^{201}\text{Hg}_{\text{surf}} + C_{\text{int}} \times \Delta^{201}\text{Hg}_{\text{int}})/(C_{\text{surf}} + C_{\text{int}}) \quad (5) $$

In which, δ$^{202}$Hg$_{surf}$, Δ$^{199}$Hg$_{surf}$ and Δ$^{201}$Hg$_{surf}$ of MMHg in the surface ocean (after photo-demethylation) are represented by those of NPO fishes at 0–100 mbs (δ$^{202}$Hg = 1.06 ± 0.35‰, Δ$^{199}$Hg = 3.84 ± 1.15‰, Δ$^{201}$Hg = 3.13 ± 0.97‰, 1 SD)$^{16,32}$, and δ$^{202}$Hg$_{int}$, Δ$^{199}$Hg$_{int}$ and Δ$^{201}$Hg$_{int}$ of MMHg in the intermediate ocean (without photo-demethylation) are represented by those of NPO marine particles (mainly IHg(II)) after correcting for −0.5‰ shift of δ$^{202}$Hg during microbial methylation of IHg(II)$^{16,32,51}$ (δ$^{202}$Hg = −0.62 ± 0.13‰, Δ$^{199}$Hg = 0.16‰ ± 0.09‰, Δ$^{201}$Hg = 0.13‰ ± 0.10‰, 1 SD). C$_{surf}$ and C$_{int}$ are represented by mean NPO surface (18 ± 7 pM, 1 SD) and intermediate (38 ± 19 pM, 1 SD) MMHg concentrations, respectively (see Supplementary Note 1). The calculated upper NPO MMHg has a weighted mean of −0.05 ± 0.30‰ (1 SD) for δ$^{202}$Hg, 1.44 ± 0.75‰ (1 SD) for Δ$^{199}$Hg and 1.16 ± 0.60‰ (1 SD) for Δ$^{201}$Hg.

**Estimate of MMHg contributions in trench fauna.** Since total Hg in the trench fauna is derived from MMHg in the upper ocean, MMHg contributions (f$_{surf}$: fraction from the surface ocean; f$_{int}$: fraction from the intermediate ocean) to the trench fauna can be estimated by the following binary mixing model using a Monte Carlo simulation approach (n = 10,000 times) through the pseudorandom number generation function of the MatLab software (R2016b, MathWorks):

$$ \Delta^{199}\text{Hg}_{\text{sam}} = f_{\text{surf}} \times \Delta^{199}\text{Hg}_{\text{surf}} + f_{\text{int}} \times \Delta^{199}\text{Hg}_{\text{int}} \quad (6) $$

$$ f_{\text{surf}} + f_{\text{int}} = 1 \quad (7) $$

where Δ$^{199}$Hg$_{sam}$, Δ$^{199}$Hg$_{surf}$, and Δ$^{199}$Hg$_{int}$ represent Δ$^{199}$Hg values of trench fauna (measured), MMHg produced in the surface ocean (3.84 ± 1.15‰, 1 SD) and MMHg produced in the intermediate ocean (0.16‰ ± 0.09‰), respectively$^{16,32}$.

**Description of statistical analysis.** All the statistical analyses were performed using OriginPro 9 except for mixing models by MatLab R2016b. Nonlinear exponential fitting was applied for Fig. 2 and linear fitting for other figures. The fitting curves are bounded by 95% confidence bands, and Pearson's R-Square and P-values are calculated by algorithms of the software. Two-side analysis of variance (ANOVA) was used to assess if the regression slope is significantly different from zero at the 0.05 level.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

The authors declare that the main data supporting the findings of this study are available within the article and its Supplementary Information file. Extra data are available from the corresponding author upon request.

## Code availability

Binary mixing model codes are available at the following link: [https://github.com/ruoyusun/Mixing-models].

Received: 7 February 2020; Accepted: 8 June 2020; Published online: 07 July 2020

## References

1. Sunderland, E. M., Li, M. & Bullard, K. Decadal changes in the edible supply of seafood and methylmercury exposure in the United States. *Environ. Health Perspect.* **126**, 017006 (2018).

2. Lamborg, C. H. et al. A global ocean inventory of anthropogenic mercury based on water column measurements. *Nature* **512**, 65–68 (2014).

3. Streets, D. G. et al. Total mercury released to the environment by human activities. *Environ. Sci. Technol.* **51**, 5969–5977 (2017).

4. Fitzgerald, W. F., Lamborg, C. H. & Hammerschmidt, C. R. Marine biogeochemical cycling of mercury. *Chem. Rev.* **107**, 641–662 (2007).

5. Podar, M. et al. Global prevalence and distribution of genes and microorganisms involved in mercury methylation. *Sci. Adv.* **1**, e1500675 (2015).

6. Lamborg, C. H., Hammerschmidt Chad, R. & Bowman Katlin, L. An examination of the role of particles in oceanic mercury cycling. *Philos. T. R. Soc. A* **374**, 20150297 (2016).

7. Zhang, Y., Jaeglé, L. & Thompson, L. Natural biogeochemical cycle of mercury in a global three-dimensional ocean tracer model. *Global Biogeochem. Cycles* **28**, 2014GB004814 (2014).

8. Cossa, D. et al. Mercury in the Southern Ocean. *Geochim. Cosmochim. Acta* **75**, 4037–4052 (2011).

9. Sunderland, E. M., Krabbenhoft, D. P., Moreau, J. W., Strode, S. A. & Landing, W. M. Mercury sources, distribution, and bioavailability in the North Pacific Ocean: insights from data and models. *Global Biogeochem. Cycles* **23**, GB2010 (2009).

10. Hammerschmidt, C. R. & Bowman, K. L. Vertical methylmercury distribution in the subtropical North Pacific Ocean. *Mar. Chem.* **132–133**, 77–82 (2012).

11. Bowman, K. L., Hammerschmidt, C. R., Lamborg, C. H. & Swarr, G. Mercury in the North Atlantic Ocean: The U.S. GEOTRACES zonal and meridional sections. *Deep Sea Res. PT. II* **116**, 251–261 (2015).

12. Munson, K. M., Lamborg, C. H., Swarr, G. J. & Saito, M. A. Mercury species concentrations and fluxes in the Central Tropical Pacific Ocean. *Global Biogeochem. Cycles* **29**, 656–676 (2015).

13. Bowman, K. L., Hammerschmidt, C. R., Lamborg, C. H., Swarr, G. J. & Agather, A. M. Distribution of mercury species across a zonal section of the eastern tropical South Pacific Ocean (U.S. GEOTRACES GP16). *Mar. Chem.* **186**, 156–166 (2016).

14. Cossa, D., Averty, B. & Pirrone, N. The origin of methylmercury in open Mediterranean waters. *Limnol. Oceanogr.* **54**, 837–844 (2009).

15. Choy, C. A., Popp, B. N., Kaneko, J. J. & Drazen, J. C. The influence of depth on mercury levels in pelagic fishes and their prey. *Proc. Natl Acad. Sci. USA* **106**, 13865–13869 (2009).

16. Blum, J. D., Popp, B. N., Drazen, J. C., Anela Choy, C. & Johnson, M. W. Methylmercury production below the mixed layer in the North Pacific Ocean. *Nat. Geosci.* **6**, 879–884 (2013).

17. Lehherr, I., Louis, V. L. S., Hintelmann, H. & Kirk, J. L. Methylation of inorganic mercury in polar marine waters. *Nat. Geosci.* **4**, 298–302 (2011).

18. Heimbürger, L.-E. et al. Methyl mercury distributions in relation to the presence of nano- and picophytoplankton in an oceanic water column (Ligurian Sea, North-western Mediterranean). *Geochim. Cosmochim. Acta* **74**, 5549–5559 (2010).

19. Monperrus, M. et al. Mercury methylation, demethylation and reduction rates in coastal and marine surface waters of the Mediterranean Sea. *Mar. Chem.* **107**, 49–63 (2007).

20. Munson, K. M., Lamborg, C. H., Boiteau, R. M. & Saito, M. A. Dynamic mercury methylation and demethylation in oligotrophic marine water. *Biogeosciences* **15**, 6451–6460 (2018).

21. Mason, R. P. & Sullivan, K. A. The distribution and speciation of mercury in the South and equatorial Atlantic. *Deep Sea Res. PT. II* **46**, 937–956 (1999).

22. Black, F. J., Poulin, B. A. & Flegal, A. R. Factors controlling the abiotic photo-degradation of monomethylmercury in surface waters. *Geochim. Cosmochim. Acta* **84**, 492–507 (2012).

23. Kim, H. et al. Methylmercury mass budgets and distribution characteristics in the western Pacific Ocean. *Environ. Sci. Technol.* **51**, 1186–1194 (2017).

24. Zhang, Y., Soerensen, A. L., Schartup, A. T. & Sunderland, E. M. A global model for methylmercury formation and uptake at the base of marine food webs. *Global Biogeochem. Cycles* **34**, e2019GB006348 (2020).
25. Sun, R. et al. Modelling the mercury stable isotope distribution of Earth surface reservoirs: implications for global Hg cycling. *Geochim. Cosmochim. Acta* **246**, 156–173 (2019).

26. Blum, J. D., Sherman, L. S. & Johnson, M. W. Mercury isotopes in Earth and environmental sciences. *Annu. Rev. Earth Planet. Sci.* **42**, 249–269 (2014).

27. Sun, G. et al. Mass-dependent and -independent fractionation of mercury isotope during gas-phase oxidation of elemental mercury vapor by atomic Cl and Br. *Environ. Sci. Technol.* **50**, 9232–9241 (2016).

28. Chen, J., Hintelmann, H., Feng, X. & Dimock, B. Unusual fractionation of both odd and even mercury isotopes in precipitation from Peterborough, ON, Canada. *Geochim. Cosmochim. Acta* **90**, 33–46 (2012).

29. Kwon, S. Y. et al. Absence of fractionation of mercury isotopes during trophic transfer of methylmercury to freshwater fish in captivity. *Environ. Sci. Technol.* **46**, 7527–7534 (2012).

30. Kwon, S. Y., Blum, J. D., Chirby, M. A. & Chesney, E. J. Application of mercury isotopes for tracing trophic transfer and internal distribution of mercury in marine fish feeding experiments. *Environ. Toxicol. Chem.* **32**, 2322–2330 (2013).

31. Madigan, D. J. et al. Mercury stable isotopes reveal influence of foraging depth on mercury concentrations and growth in Pacific bluefin tuna. *Environ. Sci. Technol.* **52**, 6256–6264 (2018).

32. Motta, L. C. et al. Mercury cycling in the North Pacific Subtropical Gyre as revealed by mercury stable isotope ratios. *Global Biogeochem. Cycles* **33**, 777–794 (2019).

33. Bidleman, T. F. et al. Scavenging amphipods: sentinels for penetration of mercury and persistent organic chemicals into food webs of the Deep Arctic Ocean. *Environ. Sci. Technol.* **47**, 5553–5561 (2013).

34. Yan, H., Li, Q., Yuan, Z., Jin, S. & Jing, M. Research progress of mercury bioaccumulation in the aquatic food chain, china: a review. *Bull. Environ. Contam. Toxicol.* **102**, 612–620 (2019).

35. Meng, M. et al. Mercury isotope variations within the marine food web of Chinese Bohai Sea: implications for mercury sources and biogeochemical cycling. *J. Hazard. Mater.* **384**, 121379 (2020).

36. Jamieson, A. J., Fujii, T., Mayor, D. J., Solan, M. & Priede, I. G. Hadal trenches: the ecology of the deepest places on Earth. *Trends Ecol. Evol.* **25**, 190–197 (2010).

37. Gerringer, M. E. On the success of the hadal snailfishes. *Integ. Organ. Biol.* **1**, 1–18 (2019).

38. Obrist, D. et al. Tundra uptake of atmospheric elemental mercury drives Arctic mercury pollution. *Nature* **547**, 201–204 (2017).

39. Bergquist, B. A. & Blum, J. D. Mass-dependent and -independent fractionation of hg isotopes by photoreduction in aquatic systems. *Science* **318**, 417–420 (2007).

40. Zheng, W. & Hintelmann, H. Mercury isotope fractionation during photoreduction in natural water is controlled by its Hg/DOC ratio. *Geochim. Cosmochim. Acta* **73**, 6704–6715 (2009).

41. Senn, D. B. et al. Stable isotope (N, C, Hg) study of methylmercury sources and trophic transfer in the northern Gulf of Mexico. *Environ. Sci. Technol.* **44**, 1630–1637 (2010).

42. Kritee, K., Barkay, T. & Blum, J. D. Mass dependent stable isotope fractionation of mercury during mer mediated microbial degradation of monomethylmercury. *Geochim. Cosmochim. Acta* **73**, 1285–1296 (2009).

43. Masbou, J. et al. Hg-stable isotope variations in marine top predators of the Western Arctic Ocean. *ACS Earth Space Chem.* **2**, 479–490 (2018).

44. Nolde, N., Drobne, D., Horvat, M. & Jereb, V. Reduction and methylation of mercury in the terrestrial isopod Porcellio scaber (Crustacea) and its environment. *Environ. Toxicol. Chem.* **24**, 1697–1704 (2005).

45. Li, H. et al. Intestinal methylation and demethylation of mercury. *Bull. Environ. Contam. Toxicol.* **102**, 597–604 (2019).

46. Gehrke, G. E., Blum, J. D. & Meyers, P. A. The geochemical behavior and isotopic composition of Hg in a mid-Pleistocene western Mediterranean sapropel. *Geochim. Cosmochim. Acta* **73**, 1651–1665 (2009).

47. Ogrinc, N., Hintelmann, H., Kotnik, J., Horvat, M. & Pirrone, N. Sources of mercury in deep-sea sediments of the Mediterranean Sea as revealed by mercury stable isotopes. *Sci. Rep.* **9**, 11626 (2019).

48. Štrok, M., Baya, P. A. & Hintelmann, H. The mercury isotope composition of Arctic coastal seawater. *C. R. Geosci.* **347**, 368–376 (2015).

49. Štrok, M., Hintelmann, H. & Dimock, B. Development of pre-concentration procedure for the determination of Hg isotope ratios in seawater samples. *Anal. Chim. Acta* **851**, 57–63 (2014).

50. Sherman, L. S. et al. Mercury isotopic composition of hydrothermal systems in the Yellowstone Plateau volcanic field and Guaymas Basin sea-floor rift. *Earth Planet. Sci. Lett.* **279**, 86–96 (2009).

51. Rodriguez-González, P. et al. Species-specific stable isotope fractionation of mercury during hg(ii) methylation by an anaerobic bacteria (desulfobulbus propionicus) under dark conditions. *Environ. Sci. Technol.* **43**, 9183–9188 (2009).

52. Janssen, S. E., Schaefer, J. K., Barkay, T. & Reinfelder, J. R. Fractionation of mercury stable isotopes during microbial methylmercury production by iron- and sulfate-reducing bacteria. *Environ. Sci. Technol.* **50**, 8077–8083 (2016).

53. Perrot, V. et al. Identical Hg isotope mass dependent fractionation signature during methylation by sulfate-reducing bacteria in sulfate and sulfate-free environment. *Environ. Sci. Technol.* **49**, 1365–1373 (2015).

54. Archer, D. E. & Blum, J. D. A model of mercury cycling and isotopic fractionation in the ocean. *Biogeosciences* **15**, 6297–6313 (2018).

55. Lumpkin, R. & Speer, K. Global ocean meridional overturning. *J. Phys. Oceanogr.* **37**, 2550–2562 (2007).

56. Wang, N. et al. Penetration of bomb 14C into the deepest ocean trench. *Geophys. Res. Lett.* **46**, 5413–5419 (2019).

57. Parks, J. M. et al. The genetic basis for bacterial mercury methylation. *Science* **339**, 1332–1335 (2013).

58. Hsu-Kim, H., Kucharzyk, K. H., Zhang, T. & Deshusses, M. A. Mechanisms regulating mercury bioavailability for methylating microorganisms in the aquatic environment: a critical review. *Environ. Sci. Technol.* **47**, 2441–2456 (2013).

59. Nunoura, T. et al. Hadal biosphere: insight into the microbial ecosystem in the deepest ocean on Earth. *Proc. Natl Acad. Sci. USA* **112**, E1230–E1236 (2015).

60. Bowman, K. L. et al. Distribution of mercury-cycling genes in the Arctic and equatorial Pacific Oceans and their relationship to mercury speciation. *Limnol. Oceanogr.* **65**, S310–S320 (2020).

61. Amos, H. M. et al. Observational and modeling constraints on global anthropogenic enrichment of mercury. *Environ. Sci. Technol.* **49**, 4036–4047 (2015).

62. Schartup, A. T. et al. Climate change and overfishing increase neurotoxicant in marine predators. *Nature* **572**, 648–650 (2019).

63. Soerensen, A. L. et al. A mass budget for mercury and methylmercury in the Arctic Ocean. *Global Biogeochem. Cycles* **30**, 560–575 (2016).

64. Jonsson, S., Mazrui, N. M. & Mason, R. P. Dimethylmercury formation mediated by inorganic and organic reduced sulfur surfaces. *Sci. Rep.* **6**, 27958 (2016).

65. Sun, R., Enrico, M., Heimbürger, L.-E., Scott, C. & Sonke, J. E. A double-stage tube furnace-acid-trapping protocol for the pre-concentration of mercury from solid samples for isotopic analysis. *Anal. Bioanal. Chem.* **405**, 6771–6781 (2013).

66. Blum, J. D. & Bergquist, B. A. Reporting of variations in the natural isotopic composition of mercury. *Anal. Bioanal. Chem.* **388**, 353–359 (2007).

## Acknowledgements

We thank P. Li for assistance with MC-ICPMS operation and Hg stable isotope analysis; J. Martin and Y. Zhang for helping constructing the Matlab codes of mixing models; and Y. Tong, Y. Liu and Tianne for discussion of this work. This study is supported by the National Key Research and Development Plan (2016YFC0300503; 2017YFC0212700), the National Natural Science Foundation of China (41625012; 41602345; 41922019; U1612442; 41961144028), and the Tianjin Natural Science Foundation (17JCYBJC23500).

## Author contributions

R.S., Y.L., and X.P. conceived and designed this project, S.C. and X.P. conducted the field sampling, and R.S., J.Y., M.M., and T.Z. carried out all the measurements. R.S. and J.E.S. led manuscript writing with supports from Y.L., W.Z., S.C.,Y.Z., J.C., and C.L. All the authors contributed to the data interpretation and manuscript revision.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41467-020-17045-3.

Correspondence and requests for materials should be addressed to Y.L. or X.P.

Peer review information *Nature Communications* thanks Kathleen Munson and other, anonymous, reviewer(s) for their contributions to the peer review of this work. Peer review reports are available.

Reprints and permission 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.
<figure><img src="image_11.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020