The global monitoring plan of the Minamata Convention on Mercury was established to generate long-term data necessary for evaluating the effectiveness of regulatory measures at a global scale. After 25 years of monitoring (since 1995), Mace Head is one of the atmospheric monitoring stations with the longest mercury record and has produced sufficient data for the analysis of temporal trends of total gaseous mercury (TGM) in Europe and the North Atlantic. Using concentration-weighted trajectories for atmospheric mercury measured at Mace Head as well as another five locations in Europe, Amderma, Andøya, Villum, Waldhof and Zeppelin, we identify the regional probabilistic source contribution factor and its changes for the period of 1996 to 2019.
Temporal trends indicate that concentrations of mercury in the atmosphere in Europe and the North Atlantic have declined significantly over the past 25 years at a non-monotonic rate averaging 0.03
Mercury is a toxic pollutant of crucial concern to public health globally. Due to its neurotoxicity, bioaccumulation, and long-range atmospheric
transport, mercury was added to the priority list of several international agreements and conventions dealing with environmental protection, including
the Minamata Convention on Mercury (e.g. Driscoll et al., 2013). Following the entry into force of the Stockholm Convention (SC) in 2004 accompanied by the Minamata Convention in 2013 to restrict releases of mercury and its compounds into the environment, a global monitoring plan was devised to evaluate the effectiveness of regulatory measures at regional and global scales. At this time, regions such as western Europe and North America have already established monitoring networks for mercury in air and precipitation, some of which have been in operation since the 1990s (Schmeltz et al., 2011; Gay et al., 2013; EMEP, 2020;
During the past decades, atmospheric mercury concentrations in the Northern Hemisphere decreased substantially (Slemr et al., 2003; Cole et al., 2014;
Steffen et al., 2015; Weigelt et al., 2015; Weiss-Penzias et al., 2016; Marumoto et al., 2019; Custodio et al., 2020). This downward trend has been attributed to decreasing emissions from the North Atlantic Ocean due to decreasing mercury concentrations in subsurface water (Soerensen et al., 2012)
and more recently to decreasing global anthropogenic emissions, mainly due to the decline in mercury release from commercial products (Horowitz et al., 2014) and the changes in
As reported by Lyman et al. (2020), the mercury emission into the atmosphere is continuously changing. Its monitoring is needed to track the trends, identify persistent and new sources, and assess the efficacy of mercury pollution control policies.
In a 5-year source apportionment study, Custodio et al. (2020) show that a factor with a high load of long-lived anthropogenic atmospheric species could explain the decrease in total gaseous mercury (TGM) at Mace Head. This decrease is consistent with a decrease in the anthropogenic mercury emissions inventory in Europe and North America (Horowitz et al., 2014). Wu et al. (2016) estimated that China's emissions also decreased since 2012, which could have a hemispheric effect. However, the downward trend of global anthropogenic mercury emissions needs to be confirmed by atmospheric observations, and a long-term evaluation of the time series of still not unknown sources and its implication should be assessed.
This study reports continuous long-term temporal trends of TGM in the North Atlantic, Arctic, and Europe, reporting mercury atmospheric concentrations at Mace Head (1995–2019), Amderma (2001–2017), Andøya (2010–2019), Villum (1999–2019), Waldhof (2005–2019), and Zeppelin
(2000–2019). Here, we combine a long-time series of atmospheric mercury observed at these sites with calculated 120
This paper aims to evaluate the TGM trend on a continental scale and the contribution of the baseline factor as a driver of the downward trend in mercury for the North Atlantic and Europe.
Based on long-range Lagrangian reanalysis backward trajectories and receptor modelling, we investigate the trends and sources of mercury in the atmosphere, assessing the inter-annual variability in the light of atmospheric transport patterns and changes in the regional emissions. In addition, we exploit atmospheric mercury temporal variability, which can be used as additional constraints to improve the ability of models to predict the cycling of mercury in the atmosphere.
Data from six sites in Europe and Greenland with the longest records of atmospheric mercury concentrations were selected for this study: Mace Head (data available for 1995–2019), Zeppelin (2000–2019), Waldhof (2006–2019), Villum (2008–2019), Andøya (2010–2019), and Amderma (2001–2013). Mace Head and Waldhof are mid-latitude stations, and Zeppelin, Amderma, and Villum can be classified as Arctic ones. Andøya, though at a latitude comparable to that of Amderma, behaves more like a mid-latitude station because the ocean around it is ice-free for most of the year. At all sites mercury had been measured by a Tekran instrument (Tekran Inc, Toronto, Canada). More details will be given at the end of the section.
The Mace Head Global Atmosphere Watch (GAW) station (53
The Zeppelin GAW station is located on the ridge of the Zeppelin mountain (78
Waldhof (52
Villum Research Station is located at the military outpost Station Nord. It is located in the furthermost north-eastern corner of Greenland on the north–south-oriented peninsula of Princess Ingeborg Halvø (81
Andøya Observatory (69.3
Amderma Polar Station is located near the Amderma settlement of the Arkhangelsk Arctic region of Russia near the coast of the Kara Sea
(69
At all sites mercury was measured using a Tekran 2537 A and/or B instrument (Tekran Inc, Toronto, Canada, mostly Model A, at Mace Head and Villum also Model B), an automated dual-channel, single-amalgamation, cold vapour atomic fluorescence (CVAFS) analyser. The instrument has two gold cartridges. While mercury is collected on one of them during the sampling period, the other is being analysed by thermodesorption and CVAFS detection. The functions of the cartridges are then alternated, allowing for quasi-continuous measurement. The instruments are usually protected by an upstream polytetrafluorethylene (PTFE) filter against dust and aerosols.
As discussed by Slemr et al. (2016), gaseous oxidized mercury (GOM) compounds are collected on the gold cartridges and were found to be converted to
elemental mercury (GEM), probably during the thermodesorption. The instrument is thus able to measure TGM provided that GOM compounds reach the cartridges. This is frequently not the case because the GOM compounds are sticky and can thus be removed on the way from the inlet
to the cartridges (Lyman et al., 2020). The instruments are usually protected by an upstream PTFE filter (mostly 0.2
Speciated mercury measurements made at Waldhof between 2009 and 2011 provided median concentrations of 6.3
To evaluate the spatial coverage and sources of air sampled at the six stations, three-dimensional reanalysis air mass back trajectories at arrival heights of 50 and 500
Concentration-weighted trajectories (CWTs) are an approach which can be used to indicate the probability of a grid cell contribution to pollution events (Cheng et al., 2013). It is based on a statistical model and can incorporate meteorological information into its analysis scheme to identify the average concentration in areas for pollutants based on a conditional probability that an air parcel that passed through a cell with a gradient
concentration displays a high concentration at the trajectory end point (Ashbaugh et al., 1985; Byčenkienė et al., 2014). The CWTs obtained at this study are a function of average mercury concentrations that were obtained every 12
The source apportionment for Mace Head was performed based on the mass conservation principle with the inclusion of potential rotated infinity
matrix transformation producing factors with a chemical profile signed by tracer species linked to its source. The full description of positive matrix factorization (PMF) and its reconstruction consideration, chemical species considered, uncertainties, and constraining of factors are presented in Custodio et al. (2020). In this
study, the PMF was applied to the Mace Head daily data. The species considered in the factorization and their mass loaded in each factor are displayed
in Fig. S3 in the Supplement in the article supplement section. In addition, the reconstructed gaseous mercury and the observation displayed
an
In this study, the assessment was performed on annual bases, and the concentrations in grid cells were calculated by counting the average concentration of trajectory segment end points that terminate within each cell as described by Byčenkienė et al. (2014) and Tang et al. (2018).
Summary of time series of GEM (
In this section, we present the time series and trends of GEM concentrations from a data set covering the periods from February 1996 to December 2019 (Mace Head), July 2001 to March 2017 (Amderma), February 2000 to December 2019 (Zeppelin), January 2006 to December 2019 (Waldhof), January 2004 to December 2019 (Andøya), and June 1999 to December 2019 (Villum). At Villum the measurements covered only 6 months (spring, summer, and early autumn) in 1999–2002, and no measurements are available for the years 2003–2008 (Skov et al., 2020). The data are summarized in Fig. 1.
GEM concentrations and their frequency distributions shown in Fig. 1 display distinct differences between the stations. GEM concentrations at Villum,
Amderma, and Zeppelin decrease frequently to values near 0 (minima of 0.0, 0.0, and 0.1
The distribution of GEM concentrations at Waldhof, a mid-latitude station in central Europe, is by contrast skewed to higher values because of frequent events with local and regional pollution (Weigelt et al., 2013). The average and median GEM concentrations at Waldhof are the highest of all the investigated stations, and the average GEM concentration is substantially higher than the median one.
The frequency distribution at Andøya is nearly symmetric, neither skewed to low nor to high GEM concentrations, although a pronounced seasonal variation can be observed. At a latitude comparable to that of Amderma, there are no pronounced depletion events at Andøya because it is exposed to the gulf stream and as such is free of ice for most of the year. Events with local and regional pollution are also largely missing at Andøya (95th
percentile of 1.79
Left panels: the density map of atmospheric mercury concentrations in different seasons. Right panels: normalized annual variation of the mercury concentrations at Arctic stations (Amderma, Villum, Zeppelin) and at the mid-latitude ones (Mace Head, Waldhof, Andøya). The shaded areas are the 95 % confidence intervals for the monthly mean.
Figure 2 shows similar seasonal variations at Mace Head, Waldhof, and Andøya, with the maximum GEM concentrations in late winter and early spring and the minimum ones in late summer and early autumn. Similar seasonal variation has been observed at most of the mid-latitude sites in the Northern Hemisphere (e.g. Cole et al., 2014; Weigelt et al., 2015; Sprovieri et al., 2016, Angot et al., 2016). It is usually accompanied by a summer maximum in
wet deposition (Gratz et al., 2009; Prestbo and Gay, 2009; Zhang and Jaeglé, 2013; Sprovieri et al., 2017), which is caused by faster oxidation of
Seasonal variations in mercury at Amderma, Villum, and Zeppelin are influenced by polar depletion events in spring and the subsequent re-emission of the deposited mercury from snow in summer, which result in pronounced GEM minima in April and May and maxima in July (Steffen et al., 2008, 2015; Dommergue
et al., 2010; Cole and Steffen, 2010; Cole et al., 2014; Angot et al., 2016; Skov et al., 2020). A similar pattern is also observed at Alert (Cole
et al., 2014). Note the larger amplitude of seasonal variation at Arctic stations (0.8–1.2
Figure 2 shows density maps which are based on the seasonal mean mercury concentration associated with respective trajectories which arrived synchronously at all six stations. The northern parts of the spring and summer panels show over the Arctic Ocean the lowest and highest mercury concentrations, respectively, which is consistent with the spring polar mercury depletion and summer emission of the mercury deposited during the depletion events. The highest GEM concentrations over the middle of the North Atlantic occur in winter and the lowest ones in summer and autumn, which is consistent with the seasonal variations at Mace Head and Andøya. High GEM levels over large parts of Europe occur in all seasons. The highest concentrations by level and extension occur in winter and spring and are somewhat lower in summer and autumn.
Kernel regression of GEM at Amderma, Andøya, Mace Head, Villum, Waldhof, and Zeppelin for the periods of 2001–2013, 2010–2019, 1995–2019, 2008–2019, 2006–2019, and 2000–2019, respectively. The smooth lines and shaded areas represent the kernel regression at the 95 % significance level. The thin lines show the monthly time series of GEM after removing annual cycles with amplitudes of 0.49, 0.23, 0.17, 0.30, 22, and 0.25
Comparison of GEM trends and average concentrations at Zeppelin, Waldhof, Andøya, Amderma, and Villum with those at Mace Head. The trends (
Figure 3 and Fig. S1 in the Supplement show the kernel regression of mercury concentrations at Mace Head, Amderma, Andøya, Villum, Waldhof, and Zeppelin. Both figures show a non-monotonic concentration change with temporary increases to
intermediate maxima at Waldhof and Zeppelin and the most pronounced ones at Villum with a maximum in 2013. The overall trend of GEM concentrations at all sites points in the downward direction. Table 1 summarizes the overall trends calculated by least square fit (LSQF) from monthly medians and compares them with those at Mace Head over the same periods of available measurements. Averages of monthly medians over the same periods are also listed. Monthly medians
were chosen to reduce the influence of depletion events at polar stations and of pollution events at mid-latitude stations. Mace Head was taken as a benchmark because of the longest and most complete data record. In addition, the trend at Mace Head represents the baseline trend (Weigelt et al.,
2015). All trends in the table are significant at the
GEM concentration at Mace Head decreased with an annual rate of
Of the Arctic stations, GEM concentration at Zeppelin decreased by only
Based on LSQF, the GEM at the Arctic stations Amderma and Villum behaves differently. The downward trends of
Concentration level (concentration-weighted trajectory) of GEM (
The higher level of atmospheric mercury at Villum in 2013 is consistent with an elevated mercury level over Greenland in that year, as deduced from back-trajectory analyses shown in Fig. 4. A large subglacial source of mercury at Greenland was recently reported by Hawkings et al. (2021). The increase in GEM at Villum in 2010 and 2013, which drives the trend up during this period, corresponds to two periods of negative extremes at Arctic Oscillation (AO). The extreme on AO and North Atlantic Oscillation (NAO) can enhance the mercury discharge from ice to the atmosphere. Bevis et al. (2019) report an anomalous ice mass loss at Greenland in the 2010–2014 epoch. The abrupt ice melting was driven mainly by changes in air temperature and solar radiation caused by atmospheric circulation anomalies.
In addition, the negative phase of the summertime NAO index increases the prevalence of high-pressure, clear-sky conditions, enhancing surface absorption of solar radiation and decreasing snowfall, and it causes the advection of warm air from southern latitudes into Greenland. These changes promote higher air temperatures, a more extended ablation season, and enhanced melt ice (Fettweis et al., 2013). In 2014/2015, when the AO indexes again turned positive and NAO negative, significant ice loss was reestablished (Bevis et al., 2019).
The back trajectories of air masses calculated for each site were combined with the measured concentration at a 12
The most revealing detail in the observed trend of GEM is displayed in Fig. 4, where it is noticeable that the downward trend is ongoing on a regional scale. This decrease could represent a change in the balance between sources and sinks of mercury in the atmosphere.
The downward trend seems to be driven by decreasing concentrations in continental Europe. This phenomenon is observed mainly after 2005, when data from Waldhof are considered. The downward trend in mercury concentration is observed in all trajectories, even in remote areas, indicated by the yellow fades to green. This phenomenon can be explained only by reductions in global atmospheric mercury sources. In addition, Fig. 4 also shows that the decrease is more pronounced in the hotspot areas identified as anthropogenic sources, where the colour shifts from dark to light red in plots from 2005 to 2019.
The later downward trend at Zeppelin and Villum (Figs. 3 and S1) suggests that these remote, high-latitude stations are less affected by direct European continental emission.
The seemingly non-monotonic downward trends with inter-annual ups and downs observed in this study are not well explained. However, an inspection of the Mace Head data (e.g. in Figs. 3 and 4) reveals that this trend is composed of two segments: one starting in 1999 and ending approximately in 2010 and a second one in 2014 after a biennial upward tendency. It could be premature to assume that the atmospheric mercury trend can be driven simply by a political decision. However, it can be seen that the two important GEM trend deflections in 1999 and 2014 coincide with COUNCIL DIRECTIVE 1999/31/EC, a European Union (EU) directive that regulates waste management of landfills in the EU, and the mercury international treaty (Minamata Convention on Mercury) designed to protect human health and the environment from anthropogenic emissions and releases of mercury approved on 10 October 2013. Continental and international environmental treaties are the result of long political and societal debate, and commitment to such a deal could reflect an already established control policy at the national level.
For example, in 1990 the United States Clean Air Act put mercury on a list of toxic pollutants that needed to be controlled to the greatest possible extent, forcing industries that release high concentrations of mercury into the environment to install maximum achievable control technologies (MACT). In 2005, the EPA promulgated a regulation that added power plants to the list of sources that should be controlled and instituted in the nation, and in 2011 new rules for coal-fired power plants were announced by the EPA (State of New Jersey et al., 2008; Castro and Sherwell, 2015).
Time series of Europe and North Atlantic mercury emission. Emission inventory provided by the Emissions Database for Global Atmospheric Research (EDGARv4.tox2, 2018). The inventory data are available at
Additionally, in 2007 the European Union implemented new mercury control measures, banning mercury in new non-electrical measuring devices, such as thermometers and barometers (Jones, 2007).
We note that Waldhof, a continental station close to anthropogenic sources in Europe, corroborates the interpretation of an anthropogenic emission-driven mercury trend. This station shows a more pronounced GEM decrease between 2005 and 2010 compared to the years since then.
An accurate emissions inventory is essential for interpreting trends in atmospheric concentrations and assessing the effectiveness of mercury
pollution control policies (Lyman et al., 2020). However, the trends in GEM concentrations have not always been consistent with those of global
anthropogenic emissions inventories, whose uncertainties ranged from
Zhang et al. (2016) presented a revised inventory of Hg emissions for the estimation of artisanal and small-scale gold mining emissions, and,
accounting for the change in
Natural sources can contribute up to 40 % of the atmospheric mercury budget (Pirrone et al., 2010); however, a trend in such a source has not been observed or reported in the literature so far.
Based on the GEM associated with each air mass trajectory, we investigated the impact of atmospheric circulation on continental Europe and the North Atlantic Ocean and observe distinct concentration patterns for the ocean and continental regions. We observed for example that air masses arriving at Mace Head from central Europe show distinct trends. We compared the regional patterns of GEM with other pollutants (
Figure 4 shows the concentration-weighted trajectory maps for GEM measurements over Mace Head, Amderma, Andøya, Villum, Waldhof, and Zeppelin. It can be seen that the highest concentrations are almost exclusively from air masses over central Europe. Exceptions are 1997 to 2000, which indicate high levels of GEM in air masses coming from the north-west. However, it should be mentioned that CWTs for this period computed only Mace Head and Villum (1999–2000) data.
The results also show a lower level of GEM in air mass segments over the North Atlantic region. This region is constantly associated with a sink of anthropogenic pollutants.
Time series (thin lines) and percentile average contribution (pie) of factors solved by PMF in the GEM reconstruction for Mace Head from 1996 to 2019, baseline (green) combustion (red), ocean (blue) and unknown factor (grey). The smooth lines and shaded areas represent the kernel regression at the 95 % significance level. The thin lines show the monthly time series with annual cycles removed.
Based on our analysis so far, our hypothesis is that the mercury concentration in North Atlantic air masses is affected by the intensity of transport from important regional and global sources and also by temporal changes in these sources. For example, the high mercury concentrations observed in the late 1990s coincide with higher contributions from continental air masses. During 2001, a noticeable reduction in the Mace Head GEM concentration was observed, corresponding to a lesser influence of continental European air masses. This was due not only to a lower frequency of air masses from continental Europe, but also lower concentrations of GEM in those air masses compared to previous years. A similar phenomenon was observed in the trend during 2005/2006 and 2008 to 2010, when an increase and decrease in the inter-annual trend corresponded to higher and lower CWTs in air masses coming from continental Europe (Fig. S2 in the Supplement).
In a 5-year source apportionment of mercury at Mace Head, Custodio et al. (2020) show that a factor with a high load of anthropogenic species could explain downward trends of atmospheric mercury. The downward trend of that factor was associated with a reduction in emissions due to cleaner manufacturing processes involving mercury and regulations limiting the emissions from coal-fired power plants since the 1980s as well as a reduction in the release of mercury from commercial products since the 1990s (Streets et al., 2011; Zhang et al., 2016).
Downward trend of GEM (blue) and baseline factor (red) at Mace Head. The smooth lines and shaded areas represent the kernel regression and 95 % significance level. The thin lines show the monthly time series with annual cycles removed. On the bottom right the correlation regression is shown with the distribution of each variable and the value of the correlation plus the significance level as stars.
Here we extend the source apportionment analysis back to 1996. The extended reconstruction of the main sources of mercury back to 1996, shown in
Fig. 6, displays a similar apportionment pattern to that reported by Custodio et al. (2020). The source apportionment indicates a baseline factor
characterized by a high load of anthropogenic species accounting for 65 % of GEM mass. The baseline factor has already been proposed as the driving factor for mercury trends at Mace Head by Custodio et al. (2020). In this study, this factor displays a downward trend of 2.7
One important consideration to take into account is that the baseline factor is interpreted as a global mercury budget from several sources which were not solved by PMF; such a factor could also take into account the strength of non-modulated extreme events or periodic oscillations such as El Niño–Southern Oscillation (ENSO) as speculated by Slemr et al. (2020) and references therein. Those events can be a reason for increased rotation in the mercury trend, imposing significance and raising the correlation.
The Global Mercury Assessment inventory (AMAP/UNEP, 2019) estimates a contribution of combustion sources to atmospheric mercury at 24 %. In this
study the combustion factor, which was indicated by a high load of
However, as reported by Custodio et al. (2020), this trend should be taken with caution since the combustion factor was fingerprinted by
The ocean factors account for 12 % of total GEM mass at Mace Head and were identified by a high load of
As reported by Custodio et al. (2020) and references therein, the residence time of mercury in the ocean is substantially longer than in the atmosphere, ranging from years to decades or millennia. Human activity has substantially increased the oceanic mercury reservoir and consequently is affecting the fluxes of mercury between the sea and atmosphere (Strode et al., 2007).
The acidification of oceans, climate change, excess nutrient inputs, and pollution are fundamentally changing the ocean's biogeochemistry (Doney, 2010) and will certainly also influence mercury ocean–air fluxes in a still unknown direction.
This study shows an upward trend in the oceanic factor after 2010, as can be seen in Fig. 6; however, its significance, implication, and causes remain to be determined.
Due to a lack of source markers that could allow the propagation of the eigenvector from axis rotation to reconstruct more realistically the complexity of mercury sources, only four factors solved our factorization. However, such an approach provided a valuable method to evaluate mercury fluxes.
A conundrum in the observed negative trend in mercury in Europe and North Atlantic over the past 2 decades is explained in this study by a decrease in anthropogenic emissions. The significant decline in concentrations of GEM over the past 2 decades demonstrates that regulatory measures across Europe have been successful in reducing the atmospheric concentration of this species, although an extensive fossil fuel use and a legacy of stockpiles in the environment continue to pose a challenge.
These results show the transport pattern of atmospheric mercury, reveal that a baseline factor with a high load of long-lived anthropogenic species dominates the source of mercury in the North Atlantic, and highlights the need for continued monitoring of the GEM and its sources. This study brings a monitoring concept for mercury on a continental scale which can be extended to a global monitoring plan by integration of the mercury monitoring network, potentially identifying hotspot concentration areas and their change over time.
This large-scale, long-term trend data evaluation can be used for assessing the effectiveness of the Minamata Convention.
More specific conclusions include the following.
Enhancement of mercury in the air masses over Greenland in summer during epochs of atmospheric circulation anomalies Mercury downward trends of 2 The observed GEM downward trend in the North Atlantic and Arctic seems to be driven by decreasing concentration in continental Europe. A baseline factor with a high load of anthropogenic species drives the mercury trend down by a strength of 97 % at the level of 0.001 ( Combustion sources could account for 20 % of GEM with a decreasing trend, and ocean sources account for 12 % with a slightly increasing
trend.
Codes used in this study are available by request to the corresponding author.
Data from this study are publicly available, as described in the experimental section.
The supplement related to this article is available online at:
DC proposed the article, processed data, and wrote the article. KAP provided data and evaluated the findings. TGS provided data and supported the writing and discussions. FFP provided data and participated in the discussion. IS supported the calculation in scripts and data assimilation besides providing meteorological and Lagrangian analysis. KAP supported the trajectory calculation and discussion. HS provided data and discussion in its interpretation. JB and RE endorsed and supported the article preparation, respectively.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This article is part of the special issue “Research results from the 14th International Conference on Mercury as a Global Pollutant (ICMGP 2019), MercOx project, and iGOSP and iCUPE projects of ERA-PLANET in support of the Minamata Convention on Mercury (ACP/AMT inter-journal SI)”. It is not associated with a conference.
The authors are grateful to the Amderma, Andøya, Mace Head, Villum, Waldhof, and Zeppelin Observatory for all data provided. The authors are grateful to the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model and the READY website (
This research has been supported by the European Commission's Horizon 2020 Framework Programme (iGOSP ERA-PLANET and e-shape “EUROGEOSS” Showcase projects).The article processing charges for this open-access publication were covered by the Helmholtz-Zentrum Hereon.
This paper was edited by Aurélien Dommergue and reviewed by two anonymous referees.