Assessing the spatial and temporal variability of MeHg biogeochemistry and bioaccumulation in the Mediterranean Sea with a coupled 3D model
- National Institute of Oceanography and Applied Geophysics - OGS
- National Institute of Oceanography and Applied Geophysics - OGS
Abstract. Human exposure to mercury (Hg) is a cause of concern, due to the biomagnification of the neurotoxic species monomethylmercury (MMHg) in marine ecosystems. Previous research revealed that commercial fish species in the Mediterranean Sea ecosystems are particularly enriched in Hg, due to a combination of physical and ecological factors. Since the fate of Hg depends on the interactions among several biogeochemical and physical drivers, biogeochemical modelling is crucial to support the integration and interpretation of field data. Here, we develop and apply a coupled transport-biogeochemical-metal bioaccumulation numerical model (OGSTM-BFM-Hg), to simulate the biogeochemical cycling of the main Hg species (HgII, Hg0, MMHg, and DMHg) in seawater, organic detritus, and through the planktonic food web. The model is applied to a 3D domain of the Mediterranean Sea to investigate the spatial and temporal variability of MeHg distribution and bioaccumulation. Model results reproduce the strong vertical and zonal gradients of MeHg concentrations related to primary production consistently with the observations, and highlight the role of winter deep convection and summer water stratification in shaping the MeHg vertical distribution, including sub-surface MeHg maximum. The modelled bioaccumulation dynamics in plankton food webs are characterized by a high spatial and temporal variability that is driven by plankton phenology, and are in agreement with available field data of concentrations in plankton and with other indicators, such as bioconcentration factors (BCFs) and trophic magnification factors (TMFs). Model results pointed out that the increment in water temperature linked to a decline of deep convection can cause an increase in water MeHg concentrations with cascading effects on plankton exposure and bioaccumulation.
- Preprint
(3601 KB) -
Supplement
(4957 KB) - BibTeX
- EndNote
Ginevra Rosati et al.
Status: final response (author comments only)
-
RC1: 'Comment on bg-2022-14', Yanxu Zhang, 28 Feb 2022
A coupled transport-biogeochemical-metal bioaccumulation model is developed in this study and is applied to simulate the biogeochemical cycling of Hg in the Mediterranean Sea. The model results reveal the spatial and temporal variation of methylmercury concentrations in this region and its pattern in the plankton food webs. I find the model scheme and parameters are up to date and has several novel features that have not been considered in previous modeling efforts, which I think are important progress in this field:
- Online coupling with biogeochemical models facilitates the sensitivity analysis of biogeochemical parameters such as POC sinking velocity. The findings that this speed can influence the depth of maximum MeHg concentrations are interesting. This is a parameter that can be constrained by observations, pointing to important future research directions.
- The inclusion of DOC-bound Hg pools, and a 4-pool partitioning scheme: POC, DOC, HgCl, and dissolved phase. Different Kd values are also used for different pools.
In addition, the authors have a detailed analysis of the seasonal cycle of Hg cycling at different locations, well constrained by observed MeHg profiles. This highlights the importance of both hydrodynamic and biogeochemical parameters.
Here are my specific comments:
Line 75: I suggest including a subsection to describe the general biogeoprovinces of the Mediterranean Sea. A very brief introduction to bathymetry and circulation patterns is also helpful. This will lay a basis for the discussion of different sites of the Sea in the later text.
Line 105-115: need a better layout for these equations.
Line 114-115: This sentence is quite misleading. If the KD values reflect a balance between adsorption and remineralization, then why KD is a constant value through the whole water column, given that the remineralization rate varies drastically?
Line 160: It's a first-order reaction approach, where does the 0.118 parameter come from? Any literature or a tunable parameter?
Line 198: Seems that only riverine load in the dissolved phase is considered? How about riverine discharge in the particulate phase, which is the dominant form of riverine Hg discharge?
Section 2.7: It is a very meaningful practice to test the sensitivity of Hg cycling to the POC sinking velocity. Changing this parameter will not only influence Hg but also C. Are there any sediment trap observations that help to constrain the POC sinking flux itself?
-
AC1: 'Comment on bg-2022-14', Ginevra Rosati, 28 Jun 2022
A coupled transport-biogeochemical-metal bioaccumulation model is developed in this study and is applied to simulate the biogeochemical cycling of Hg in the Mediterranean Sea. The model results reveal the spatial and temporal variation of methylmercury concentrations in this region and its pattern in the plankton food webs. I find the model scheme and parameters are up to date and has several novel features that have not been considered in previous modeling efforts, which I think are important progress in this field:
- Online coupling with biogeochemical models facilitates the sensitivity analysis of biogeochemical parameters such as POC sinking velocity. The findings that this speed can influence the depth of maximum MeHg concentrations are interesting. This is a parameter that can be constrained by observations, pointing to important future research directions.
- The inclusion of DOC-bound Hg pools, and a 4-pool partitioning scheme: POC, DOC, HgCl, and dissolved phase. Different Kd values are also used for different pools.
In addition, the authors have a detailed analysis of the seasonal cycle of Hg cycling at different locations, well constrained by observed MeHg profiles. This highlights the importance of both hydrodynamic and biogeochemical parameters.
Here are my specific comments:
Line 75: I suggest including a subsection to describe the general biogeoprovinces of the Mediterranean Sea. A very brief introduction to bathymetry and circulation patterns is also helpful. This will lay a basis for the discussion of different sites of the Sea in the later text.
Line 105-115: need a better layout for these equations.
Line 114-115: This sentence is quite misleading. If the KD values reflect a balance between adsorption and remineralization, then why KD is a constant value through the whole water column, given that the remineralization rate varies drastically?
Line 160: It's a first-order reaction approach, where does the 0.118 parameter come from? Any literature or a tunable parameter?
Line 198: Seems that only riverine load in the dissolved phase is considered? How about riverine discharge in the particulate phase, which is the dominant form of riverine Hg discharge?
Section 2.7: It is a very meaningful practice to test the sensitivity of Hg cycling to the POC sinking velocity. Changing this parameter will not only influence Hg but also C. Are there any sediment trap observations that help to constrain the POC sinking flux itself?
-
RC2: 'Comment on bg-2022-14', Gwenaël Abril, 31 May 2022
Dear Authors
please consider the following comments by an anonymous referee.
Best Regard
G Abril, asssociate editor
________________
Anonymous review
The paper describe a modelling study of Hg cycle in the Mediterranean Sea using a high resolution regional model. It is a very interesting and important modelling effort since It is the first attempt to simulate this contaminant in this region with a coupled dynamical-biogeochemical model. The paper is well written and results are globally well illustrated and discussed. The modelling approach considers the main Hg species in the ocean and the related processes that control their exchanges and redistributions. The analysis focuses on MeHg species that represent the toxic species for ecosystems, and investigates the transfer and bioaccumulation in the low trophic level, the planktonic reservoir (phytoplankton and zooplankton). It provides very new interesting and quantified informations on the spatial and temporal distribution of MeHg species in the different regions of the Mediterranean Sea. I then recommend this paper for publication, but suggest also some major revision on the description of the model. The paper focuses mainly on MeHg species, but the cycling of Hg in the Mediterranean sea is controlled by the distribution and exchange among all the different species that are not well document in this paper, while it is of importance to assess the consistency of the results.
Recently an assessment of Mercury in the Mediterranean Sea has been published (Cossa et al, 2022). It provides constraints on sources of hg in the Mediterranean Sea (atmospheric, riverine, sediments) , exchange fluxes at Gibraltar with the Atlantic Ocean, and budget of THg and MeHg in the western and Eastern basin. I suggest that before publication, , this paper compares in its supplementary material, its modelling results with the budget derived by this Mediterranean Hg assessments, in order to verify that the global modelling approach is coherent with observations.
Cossa et al. Mediterranean Mercury Assessment 2022: An Updated Budget, Health Consequences, and Research Perspectives. Environ. Sci. Technol. 2022, 56, 3840−3862
-
AC2: 'Reply on RC2', Ginevra Rosati, 28 Jun 2022
We thank the anonymous reviewer for his/her positive and constructive comments that will improve the manuscript quality. Please find our replies to the comments attached as supplemental pdf to the present comment.
Kind regards,
Ginevra Rosati on behalf of the coauthors
-
AC2: 'Reply on RC2', Ginevra Rosati, 28 Jun 2022
- AC1: 'Comment on bg-2022-14', Ginevra Rosati, 28 Jun 2022
Ginevra Rosati et al.
Model code and software
OGSTM-BFM-Hg model code Rosati Ginevra, Canu Donata, Lazzari Paolo, Solidoro, Cosimo https://doi.org/10.5281/zenodo.5851442
Ginevra Rosati et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
255 | 99 | 15 | 369 | 40 | 7 | 6 |
- HTML: 255
- PDF: 99
- XML: 15
- Total: 369
- Supplement: 40
- BibTeX: 7
- EndNote: 6
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1