# Particulate-phase mercury emissions from biomass burning and impact on resulting deposition: a modelling assessment

Francesco de Simone, Paulo Artaxo, Mariantonia Bencardino, Sergio Cinnirella, Francesco Carbone, Francesco d'Amore, Aurélien Dommergue, Xin Bin Feng, Christian N. Gencarelli, Ian M. Hedgecock, et al.

► To cite this version:

Francesco de Simone, Paulo Artaxo, Mariantonia Bencardino, Sergio Cinnirella, Francesco Carbone, et al.. Particulate-phase mercury emissions from biomass burning and impact on resulting deposition: a modelling assessment. Atmospheric Chemistry and Physics, 2017, 17, pp.1881-1899. 10.5194/acp-17-1881-2017. insu-03706530

HAL Id: insu-03706530

https://insu.hal.science/insu-03706530

Submitted on 28 Jun 2022

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>

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

# Particulate-phase mercury emissions from biomass burning and impact on resulting deposition: a modelling assessment

Francesco De Simone<fnref n="1" />, Paulo Artaxo<fnref n="2" />, Mariantonia Bencardino<fnref n="1" />, Sergio Cinnirella<fnref n="1" />, Francesco Carbone<fnref n="1" />, Francesco D'Amore<fnref n="1" />, Aurélien Dommergue<fnref n="3" />, Xin Bin Feng<fnref n="4" />, Christian N. Gencarelli<fnref n="1" />, Ian M. Hedgecock<fnref n="1" />, Matthew S. Landis<fnref n="5" />, Francesca Sprovieri<fnref n="1" />, Noriuki Suzuki<fnref n="6" />, Ingvar Wängberg<fnref n="7" />, and Nicola Pirrone<fnref n="8" />

*Correspondence to:* Francesco De Simone (f.desimone@iia.cnr.it)

Received: 29 July 2016 – Discussion started: 1 August 2016

Revised: 9 January 2017 – Accepted: 12 January 2017 – Published: 8 February 2017

**Abstract.** Mercury (Hg) emissions from biomass burning (BB) are an important source of atmospheric Hg and a major factor driving the interannual variation of Hg concentrations in the troposphere. The greatest fraction of Hg from BB is released in the form of elemental Hg ($Hg_{(g)}^0$). However, little is known about the fraction of Hg bound to particulate matter ($Hg^P$) released from BB, and the factors controlling this fraction are also uncertain. In light of the aims of the Minamata Convention to reduce intentional Hg use and emissions from anthropogenic activities, the relative importance of Hg emissions from BB will have an increasing impact on Hg deposition fluxes. Hg speciation is one of the most important factors determining the redistribution of Hg in the atmosphere and the geographical distribution of Hg deposition. Using the latest version of the Global Fire Emissions Database (GFEDv4.1s) and the global Hg chemistry transport model, ECHMERIT, the impact of Hg speciation in BB emissions, and the factors which influence speciation, on Hg deposition have been investigated for the year 2013. The role of other uncertainties related to physical and chemical atmospheric processes involving Hg and the influence of model parametrisations were also investigated, since their

interactions with Hg speciation are complex. The comparison with atmospheric $Hg^P$ concentrations observed at two remote sites, Amsterdam Island (AMD) and Manaus (MAN), in the Amazon showed a significant improvement when considering a fraction of $Hg^P$ from BB. The set of sensitivity runs also showed how the quantity and geographical distribution of $Hg^P$ emitted from BB has a limited impact on a global scale, although the inclusion of increasing fractions $Hg^P$ does limit $Hg_{(g)}^0$ availability to the global atmospheric pool. This reduces the fraction of Hg from BB which deposits to the world's oceans from 71 to 62%. The impact locally is, however, significant on northern boreal and tropical forests, where fires are frequent, uncontrolled and lead to notable Hg inputs to local ecosystems. In the light of ongoing climatic changes this effect could be potentially be exacerbated in the future.

<footnote n="1">¹CNR-Institute of Atmospheric Pollution Research, Division of Rende, UNICAL-Polifunzionale, 87036 Rende, Italy</footnote>

<footnote n="2">²University of Sao Paulo, Sao Paulo, Brazil</footnote>

<footnote n="3">³Univ. Grenoble Alpes, CNRS, IRD, IGE, Grenoble, France</footnote>

<footnote n="4">⁴Institute of Geochemistry, State Key Laboratory of Environmental Geochemistry, Chinese Academy of Sciences, Guiyang, China</footnote>

<footnote n="5">⁵Office of Research and Development, US Environmental Protection Agency, Research Triangle Park, NC, USA</footnote>

<footnote n="6">⁶National Institute for Environmental Studies (NIES), Ministry of Environment, Okinawa, Japan</footnote>

<footnote n="7">⁷IVL, Swedish Environmental Research Inst. Ltd., Göteborg, Sweden</footnote>

<footnote n="8">⁸CNR-Institute of Atmospheric Pollution Research, Area della Ricerca di Roma 1, Via Salaria km 29 300, Monterotondo, 00015 Rome, Italy</footnote>
# 1 Introduction

Emissions from biomass burning (BB) are an important source of mercury (Hg) to the atmosphere (De Simone et al., 2015; Friedli et al., 2009) and a major factor in determining the interannual variations of its tropospheric concentration (Slemr et al., 2016). Although the Hg released by BB varies from year to year, it can amount to up to roughly one third of the anthropogenic emission estimates (AMAP/UNEP, 2013; Friedli et al., 2009; De Simone et al., 2015). With the eventual implementation of the Minamata Convention (http://www.mercuryconvention.org/) and future curbs on industrial emission, as a by-product of industrial emission abatement measures, its relative importance will increase in the coming years. A previous modelling study (De Simone et al., 2015) used the global Hg chemistry model, ECHMERIT, and three BB inventories to assess the distribution of Hg deposition resulting from BB. A large part of the Hg released from BB deposits over oceans, where its re-emission is driven by sea surface temperature, among other factors (Carbone et al., 2016; Andersson et al., 2011), or where it can be converted to toxic methyl mercury (MeHg) compounds, has important implications for the food web and, through fish consumption, also for human health (see Chen et al., 2016, and references therein). The deposition flux of Hg from BB has been shown to be more sensitive to certain factors, in particular the chemical mechanism employed in the model and the choice of emission inventory, than to others such as the vertical profiles of emissions (De Simone et al., 2015). In this previous study all Hg emitted from BB was considered to be $Hg_{(g)}^0$. There is, however, evidence that the fraction of Hg emitted bound to particulates ($Hg^P$) may be sizeable, up to 30%, especially when the fuel moisture content (FMC) is high (Obrist et al., 2007; Finley et al., 2009; Friedli et al., 2009; Wang et al., 2010). These levels, however, remain uncertain since different methodologies have led to different conclusions (Zhang et al., 2013; Obrist et al., 2007). Little is known about the mechanisms that control the speciation of Hg in BB emissions, which leads to uncertainties in the Hg deposition patterns, since the atmospheric lifetime of $Hg^P$ is significantly shorter than $Hg_{(g)}^0$, leading to greater local deposition. Local Hg deposition due to BB could have important repercussions in regions such as the South-East Asia, where there is intensive rice cultivation, which is subject to major BB events, especially during El Niño periods. Hg deposited to rice paddies can be readily converted to toxic MeHg that can accumulate in the grains (Wang et al., 2015; Feng et al., 2008; Meng et al., 2014; Zhang et al., 2010). Moreover, it has been reported that $Hg^P$ from BB deposited to foliage has the ability to enhance MeHg formation (Witt et al., 2009). The aim of this study is to investigate the effects on simulated deposition fluxes of Hg resulting from BB when variations in $Hg^P$ fraction and production processes are considered. The most recent version of the GFED BB emission inventory (van der Werf et al., 2010; Randerson et al., 2012; Mu et al., 2011),

has been included in the global online Hg chemical transport model ECHMERIT to simulate Hg deposition from BB for the year 2013 and to quantify the influence of variations in model inputs, assumptions and parametrisations.

# 2 Methods

## 2.1 The biomass burning inventory

The reference BB inventory in this study, Global Fire Emissions Database version 4 (GFED4.1s), is based on an updated version of the inventory of van der Werf et al. (2010) with burned area from Giglio et al. (2013), and with the addition of small fire-burned area (Randerson et al., 2012). The standard temporal resolution of the emissions files is monthly, but data are provided to distribute these daily, and a diurnal cycle based on Mu et al. (2011) is also available. Daily BB emissions from two other global inventories, GFASv1.2 (Kaiser et al., 2012, 2015) and FINNv1.5 (Wiedinmyer et al., 2011), were also included in the model for sensitivity runs. These three inventories are all compiled using the imagery obtained from the MODIS instruments. However, the way in which the data are filtered or processed yields substantial differences between the final products; see Andela et al. (2013) and references therein for a detailed description of the differences among the inventories.

## 2.2 Experimental set-up

The global Hg chemical transport model ECHMERIT (Jung et al., 2009; De Simone et al., 2014) uses T42 horizontal resolution (roughly 2.8° by 2.8° at the Equator) and 19 vertical levels up to 10 hPa. Hg emissions from BB were included in the model by mapping them to CO emissions using the global averaged enhancement ratio (ER) of $1.96 \times 10^{-7}$, as obtained by Friedli et al. (2009), averaging field measurements from different biomes in various regions around the globe, including in plume measurements from the CARIBIC project (Ebinghaus et al., 2007). Previous modelling studies have used different ERs (De Simone et al., 2015; Holmes et al., 2010), but all these values were well within the range of uncertainty ($0.3–6.0 \times 10^{-7}$; see Wang et al., 2015). ECHMERIT, in the base configuration, includes the oxidation of $Hg_{(g)}^0$ to in $Hg_{(g/aq)}^{II}$ oxidation by $O_3/OH$ in the gas and aqueous phases. $OH$ and $O_3$ concentration fields were imported from MOZART (Model for Ozone and Related chemical Tracers) (Emmons et al., 2010). $Hg^P$ is assumed to be inert, whether it is emitted from anthropogenic activities or BB, and it is subject to transport and deposition processes but is not involved in any chemical reactions. Mechanisms and parametrisations used for calculating the dry and the wet deposition of the different Hg species are the same as described in Jung et al. (2009). Beyond this standard configuration a number of alternative processes and chemical mechanisms have been considered for this study, as explained in
Sect. 2.3. Atmospheric reduction of Hg<sub>(g/aq)</sub><sup>II</sup> to Hg<sub>(g)</sub><sup>0</sup> has been included in many models to regulate the residence time of Hg<sub>(g)</sub><sup>0</sup> in the atmosphere. However, a number of the proposed mechanisms are unlikely to occur under most atmospheric conditions or are based on empirical rates to better match the observations (see Kwon and Selin (2016) for a recent review). Due to this uncertainty, reduction was not included in this study. No further Hg<sup>P</sup> particulate matter (PM) dimension distributions other than the standard log-normal particle size distribution, as described in detail in (Jung et al., 2009), were considered in this study due to large uncertainties regarding the dynamic size range of PM emitted during BB (see Janhäll et al. (2010) and references therein). GFED4.1s provides monthly burned area, fire carbon (C) and dry matter (DM) emissions (http://www.falw.vu/~gwerf/GFED/GFED4/). A script is provided to derive gaseous and PM emissions from DM fields making use of biome-based emission factors based on Akagi et al. (2011) and van der Werf et al. (2010). The resulting emission fields were then interpolated on to the ECHMERIT T42 grid using the mass conserving remapping function included in the Climate Data Operators (https://code.zmaw.de/projects/cdo).

## 2.3 Simulations and their scope

The BASE simulation used as the reference case in this study includes daily BB emissions from GFEDv4.1s, in which a global uniform fraction of Hg<sup>P</sup>, equal to 15% of the total Hg emission is assumed. This value is within the range of observations (Obrist et al., 2007; Finley et al., 2009). However, since there are uncertainties in the proportion of Hg<sup>P</sup> emitted from BB (Zhang et al., 2013), further simulations were carried out with varying fractions of Hg<sup>P</sup> (0, 4 and 30%). Simulations were also conducted mapping the 15% of the total Hg emitted as Hg<sup>P</sup> to the geographical distribution of different proxy chemical species (see Sect. 2.4). The shorter lifetime of Hg<sup>P</sup> with respect to Hg<sub>(g)</sub><sup>0</sup> potentially means that the vertical profile of the emissions could have an impact on the distribution of Hg deposition, as is the case for other specified Hg emission sources (De Simone et al., 2016). Therefore two vertical profile parametrisations, as well as different emission injection time resolutions, were also included in the study. The principal vertical profile used (PBL-Profile) maps the Hg emissions uniformly within the planetary boundary layer (PBL), whereas in the second the vertical profile of the standard version of the ECHAM-HAM model was used (HAM-Profile) (Zhang et al., 2012). The HAM-Profile is equal to PBL-Profile when the PBL height is greater than 4000 m; otherwise 75% of the emissions are placed within the PBL and the remainder in the two layers above the PBL (17 and 8%). This threshold value is arbitrary, but it is the standard configuration of ECHAM6-HAM2 (Zhang et al., 2012; Veira et al., 2015). Biomass burning emissions from GFASv1.2 (Kaiser et al., 2012, 2015) and FINNv1.5 (Wiedinmyer et al., 2011) were also used in the study to assess

uncertainty related to the satellite imagery processing and inventory compilation. Simulations using GFASv1.2 were excluded from subsequent analyses since the low Hg emissions could be due to a technical problem arising from GRIB encoding (see GFAS, 2015). These simulations primarily employ a O<sub>3</sub> / OH Hg<sub>(g)</sub><sup>0</sup> oxidation mechanism. However, since the precise atmospheric Hg oxidation mechanism remains unclear (Hynes et al., 2009; Subir et al., 2011, 2012; Gustin and Jaffe, 2010; Gustin et al., 2015; Ariya et al., 2015), a number of runs were performed using a Br-based oxidation mechanism. Some studies (Steffen et al., 2014; Amos et al., 2012) suggest that the partitioning of reactive Hg species the between gas and particulate phases might be driven by air temperature and on the surface are of the aerosol present in the atmosphere. Therefore, two other simulations were conducted including the temperature-dependent gas-particle partitioning described in Amos et al. (2012), one assuming BB Hg emissions to be only Hg<sub>(g)</sub><sup>0</sup> and another assuming a 15% of BB Hg emissions to be Hg<sup>P</sup>. To estimate the ratio of Hg deposition from BB compared to anthropogenic sources, six further simulations were conducted including only anthropogenic emissions using the EDGAR (Muntean et al., 2014), AMAP2010 (AMAP/UNEP, 2013) and STREETS (Corbitt et al., 2011) inventories, employing the O<sub>3</sub> / OH and Br oxidation mechanisms. This study covers a single year, 2013, chosen due to the availability of measurements from GMOS network (Sprovieri et al., 2016a, b; D'Amore et al., 2015). All simulations were performed for a full year, without the rapid re-emission mechanism (Selin et al., 2008), and were continued without further emissions for another 12 months to allow most of the 2013 Hg emissions to be deposited. Finally, a selection of simulations were re-run including Hg emissions from all sources, BB, anthropogenic emissions from AMAP2010 (AMAP/UNEP, 2013), dynamic ocean emissions, terrestrial emissions and re-emissions as described in De Simone et al. (2014), to evaluate model performance against measurements and to evaluate the assumptions made in this study. A summary of the simulations performed can be found in Table 1.

## 2.4 BB emission speciation

The release of Hg from BB occurs prevalently as Hg<sub>(g)</sub><sup>0</sup>. However, as mentioned previously, a measurable fraction may be emitted as Hg<sup>P</sup> (Obrist et al., 2007; Friedli et al., 2009; Finley et al., 2009; Wang et al., 2010). No significant amounts of gaseous oxidised Hg (Hg<sub>(g)</sub><sup>II</sup>) have so far been detected in BB emissions (Obrist et al., 2007, and references therein). The speciation of Hg emissions is of great importance, since it largely determines the atmospheric lifetime and hence the distance emitted Hg is transported in the atmosphere before deposition, as seen for other speciated Hg sources (Bieser et al., 2014). The fraction of Hg<sup>P</sup> released by BB determined in field and laboratory studies ranges from fractions of a few percent to over 30% (Obrist et al., 2007). The factors deter-
<table><caption>Table 1. Simulations performed.</caption><thead><tr><th>Name</th><th>Inventory (BB emission (Mg))</th><th>Full version</th><th>Emiss. time res.</th><th>Fraction Hg<sup>P</sup> (%)</th><th>Map Hg<sup>P</sup></th><th>Chem. mech.</th><th>Vertical profile</th><th>Scope</th></tr></thead><tbody><tr><td>BASE</td><td>GFED4.1s (390)</td><td rowspan="10">Yes</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Reference</td></tr><tr><td>3-hourly</td><td>GFED4.1s (390)</td><td>3 h</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Emiss. time resol.</td></tr><tr><td>Monthly</td><td>GFED4.1s (390)</td><td>monthly</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Emiss. time resol.</td></tr><tr><td>HAM-Profile</td><td>GFED4.1s (390)</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>HAM</td><td>Vertical profile</td></tr><tr><td>Only first level</td><td>GFED4.1s (390)</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>1st</td><td>Vertical profile</td></tr><tr><td>Only PBL level</td><td>GFED4.1s (390)</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>level of PBL</td><td>Vertical profile</td></tr><tr><td>3 h + HAM-prof</td><td>GFED4.1s (390)</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>HAM</td><td>V. Pr. &amp; E. T. res.</td></tr><tr><td>Hg<sup>P</sup> to PM</td><td>GFED4.1s (390)</td><td>Yes</td><td>daily</td><td>15</td><td>PM</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Hg<sup>P</sup> mapping</td></tr><tr><td>Hg<sup>P</sup> to OC</td><td>GFED4.1s (390)</td><td>Yes</td><td>daily</td><td>15</td><td>OC</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Hg<sup>P</sup> mapping</td></tr><tr><td>Hg<sup>P</sup> to FMC</td><td>GFED4.1s (390)</td><td>Yes</td><td>daily</td><td>variable</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Hg<sup>P</sup> mapping</td></tr><tr><td>NO Hg<sup>P</sup></td><td>GFED4.1s (390)</td><td rowspan="4">Yes</td><td>daily</td><td>0</td><td>NA</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Fraction Hg<sup>P</sup></td></tr><tr><td>4 % Hg<sup>P</sup></td><td>GFED4.1s (390)</td><td>daily</td><td>4</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Fraction Hg<sup>P</sup></td></tr><tr><td>30 % Hg<sup>P</sup></td><td>GFED4.1s (390)</td><td>daily</td><td>30</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Fraction Hg<sup>P</sup></td></tr><tr><td>100 % Hg<sup>P</sup></td><td>GFED4.1s (390)</td><td>daily</td><td>100</td><td>CO</td><td>None</td><td>PBL</td><td>Transport Hg<sup>P</sup></td></tr><tr><td>Partitioning</td><td>GFED4.1s (390)</td><td>Yes</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Partitioning Hg<sup>P</sup>/II</td></tr><tr><td>Partitioning ref.</td><td>GFED4.1s (390)</td><td>Yes</td><td>daily</td><td>0</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Partitioning Hg<sup>P</sup>/II</td></tr><tr><td>Reduction</td><td>GFED4.1s (390)</td><td>Yes</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH + Red.</td><td>PBL</td><td>Chemistry</td></tr><tr><td>Br</td><td>GFED4.1s (390)</td><td rowspan="7">Yes</td><td>daily</td><td>15</td><td>CO</td><td>Br</td><td>PBL</td><td>Chemistry</td></tr><tr><td>Br No Hg<sup>P</sup></td><td>GFED4.1s (390)</td><td>0</td><td>NA</td><td>Br</td><td>PBL</td><td>Chemistry</td></tr><tr><td>Br 30 % Hg<sup>P</sup></td><td>GFED4.1s (390)</td><td>daily</td><td>30</td><td>CO</td><td>Br</td><td>PBL</td><td>Chemistry</td></tr><tr><td>Br Hg<sup>P</sup> to OC</td><td>GFED4.1s (390)</td><td>daily</td><td>15</td><td>OC</td><td>Br</td><td>PBL</td><td>Chemistry</td></tr><tr><td>Br Hg<sup>P</sup> to FMC</td><td>GFED4.1s (390)</td><td>daily</td><td>variable</td><td>CO</td><td>Br</td><td>PBL</td><td>Chemistry</td></tr><tr><td>GFAS</td><td>GFASv1.2 (150; see Sect. 2.3)</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Inventory</td></tr><tr><td>GFAS Br</td><td>GFASv1.2 (150; see Sect. 2.3)</td><td>daily</td><td>15</td><td>CO</td><td>Br</td><td>PBL</td><td>Inventory</td></tr><tr><td>FINN</td><td>FINNv1.5 (550)</td><td rowspan="2">Yes</td><td>daily</td><td>15</td><td>CO</td><td>O<sub>3</sub> / OH</td><td>PBL</td><td>Inventory</td></tr><tr><td>FINN Br</td><td>FINNv1.5 (550)</td><td>daily</td><td>15</td><td>CO</td><td>Br</td><td>PBL</td><td>Chemistry</td></tr><tr><td>AMAPOH</td><td>AMAPOH2010</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>O<sub>3</sub> / OH</td><td>NA</td><td>Ratio to anth. emiss.</td></tr><tr><td>AMAPBr</td><td>AMAPOH2010</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>Br</td><td>NA</td><td>Ratio to anth. emiss.</td></tr><tr><td>EDGAROH</td><td>EDGAR2008</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>O<sub>3</sub> / OH</td><td>NA</td><td>Ratio to anth. emiss.</td></tr><tr><td>EDGARBr</td><td>EDGAR2008</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>Br</td><td>NA</td><td>Ratio to anth. emiss.</td></tr><tr><td>STREETSOH</td><td>STREETS2005</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>O<sub>3</sub> / OH</td><td>NA</td><td>Ratio to anth. emiss.</td></tr><tr><td>STREETSBr</td><td>STREETS2005</td><td></td><td>NA</td><td>NA</td><td>NA</td><td>Br</td><td>NA</td><td>Ratio to anth. emiss.</td></tr></tbody></table>

mining speciation and whether HgP is directly emitted or if it is the product of the oxidation of Hg(g)0 within the plume (Obrist et al., 2007; Webster et al., 2016) are not known. However, foliage, moisture content, fuel type, plant species and combustion proprieties certainly play a role. HgP emissions were found to be well correlated with particulate matter (PM) and organic carbon (OC) emissions (Obrist et al., 2007). Obrist et al. (2007) found that Hg(g)0 is the dominant species in dry fuel combustion, whereas the fraction of HgP becomes appreciable when FMC reaches roughly 30%, above which HgP release appears to increase linearly with FMC. In the inventory used for the BASE case both Hg(g)0 and HgP follow the spatial distribution of CO emissions from BB, and 15% of the emitted Hg is considered to be HgP (see Figs. 1a and 2a). Hg emission fields were also compiled in which the HgP fraction of the total Hg emitted was mapped to OC and PM emissions (see Fig. 2b and c). A further emission field was compiled in which the ratio of Hg(g)0 to HgP is determined by the FMC (Figs. 1b and 2d). A relationship was found to exist between HgP emissions and the fire burn duration and severity as well as combustion conditions (Obrist et al., 2007; Webster et al., 2016). In particular high HgP fractions were observed during smouldering phases, whereas very low or undetectable HgP levels were found during flaming combustion. These potential parametrisations were not investigated here due to the difficulty in finding a suitable

proxy data set. Appendix A contains a more detailed description of the methods used to calculate the different Hg BB emission fields.

# 3 Results

## 3.1 Emissions

The total Hg emitted in 2013 based on the GFED inventory is roughly 400 Mg, which is at the lower end of the initial estimates (675 ± 240 Mg) (Friedli et al., 2009) but is reasonable considering the natural variation of BB activity and the diminishing trend of the CO emission estimates in the latest inventory revisions (up to 50% for some years) (van der Werf et al., 2010). Considering 15% of the emissions to be HgP, in the BASE run this corresponds to approximately 340 Mg Hg(g)0 and 60 Mg HgP. Interestingly the emissions of HgP amount to 58 Mg when relating the HgP fraction to FMC. The exact amount of Hg emitted by BB in the different model runs is detailed in Table 1. The spatial distribution and the vertical profile of the emission injection height, considering the PBL-Profile for Hg(g)0 and HgP in the different cases considered are shown in Figs. 1 and 2. Both the geographical and vertical distributions of the emissions of the Hg species reveal notable differences depending on the methodology used, particularly for HgP. Compared to the cases where HgP emis-
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Geographical distribution (a–b) and PBL-type vertical profiles (c–d) of the Hg<sub>(g)</sub><sup>0</sup> emissions, when mapped to CO (a, c) and when speciation is determined by FMC (b, d). For the emissions mapped to CO, only the speciation (15 : 85 Hg<sup>P</sup> : Hg<sub>(g)</sub><sup>0</sup>) is shown for clarity.</figcaption></figure>

sions are mapped to CO and PM (Fig. 2a–b and e–f), map-
ping Hg<sup>P</sup> to OC and using the FMC to determine the speci-
ation (Fig. 2c–d and g–h) result in enhanced Hg<sup>P</sup> emissions,
above 60° N and over some areas the Amazon, central Africa
and East Asia as evident in Fig. 3. The timing and location
of the enhanced Hg<sup>P</sup> emission at northerly latitudes could be
particularly relevant for Hg deposition to the Arctic. From
Fig. 3 it is evident how the geographical distribution of the
Hg<sup>P</sup> to Hg<sub>(g)</sub><sup>0</sup> emission ratio differs with the assumptions con-
sidered. However, for OC and FMC there is general agree-
ment on the areas where the Hg<sup>P</sup> emissions are relatively
higher, especially in the Northern Hemisphere and particu-
larly for areas above 60° N. The agreement between OC and
FMC is not surprising and is related to the combustion char-
acteristics that enhance OC emissions, i.e. lower combustion
temperatures and the dominance of the smouldering phase
of combustion (Zhang et al., 2013), that are likely to occur
where FMC is greatest.

**3.2 Emission latitudinal profiles**

The latitudinal profiles of Hg<sub>(g)</sub><sup>0</sup> and Hg<sup>P</sup> emissions, using
the different approaches (Sect. 2.4), are shown in Fig. 4a
and b. For those emissions mapped to CO, only the 15 : 85
(Hg<sup>P</sup> : Hg<sub>(g)</sub><sup>0</sup>) speciation is reported for clarity. The dif-
ferences in the latitudinal profiles of the Hg<sub>(g)</sub><sup>0</sup> emissions

(Fig. 4a) are sizeable only for the peaks north of 45° N, where
the FMC-based speciation has an Hg<sub>(g)</sub><sup>0</sup> fraction below 85 %.
The latitudinal profiles of Hg<sup>P</sup> emissions mapped to PM and
CO look very similar over the entire domain (Fig. 4b), apart
from a peak a few degrees north of the Equator. The Hg<sup>P</sup>
emissions mapped to OC and FMC differ from the PM and
CO profiles but are similar to each other between roughly
30° S and 60° N. South of 30° S Hg<sup>P</sup> emissions mapped to
OC are higher, while peak Hg<sup>P</sup> emissions derived from FMC
at 65° N (1.5 g km<sup>-2</sup> yr<sup>-1</sup>) are nearly 30 % greater than those
derived from OC and roughly double those mapped to CO
and PM. Moreover, in the FMC scenario the peak in Hg<sup>P</sup>
emissions at 65° N are greater than the peak seen at 15° S (1.5
vs. 1.4 g km<sup>-2</sup> yr<sup>-1</sup>). As is particularly evident in Fig. 4c, the
most notable differences among the different assumptions
hypothesised are above 60° N, where both the OC and the
FMC cases agree on the location of the greatest Hg<sup>P</sup> emis-
sions probably due to the linkage between OC emissions
and combustion processes favoured by FMC (Zhang et al.,
2013), and between 30 and 45° S, where only OC and PM
are greater than BASE. A previous modelling study focus-
ing on the fate of Hg from BB, where all emissions were
considered as Hg<sub>(g)</sub><sup>0</sup>, showed that the long atmospheric life
of the elemental Hg smoothed the deposition latitudinal pro-
files compared to the emission profiles (De Simone et al.,
2015). The four panels in Fig. 5 compare the normalised lati-
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Geographical distribution (a–d) and PBL-type vertical profiles (e–h) of the Hg$^P$ emissions as injected in the model, when mapped to CO (a, e), PM (b, f) and OC (c, g) and when speciation is determined by FMC (d, h). For the emissions mapped to CO, only the speciation (15: 85 Hg$^P$:Hg$^0$) is shown for clarity.</figcaption></figure>

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

tudinal deposition profiles obtained for the BASE simulation with those obtained from the alternative Hg$^P$ emission scenarios by category. Figure 5a demonstrates the very limited impact of the time resolution used for BB emissions, most likely due to the coarse horizontal resolution of the model. The two vertical emission profiles (Fig. 5b) give deposition fields that are to all effects indistinguishable, even when considering varying temporal resolution of the BB emissions, whereas assuming all emissions to be in the first model level

(with an average height of approximately 35 m) leads to enhanced deposition near emission peaks. In this instance, the maximum deposition coincides with peak emission, at approximately 15° S, whereas in all other cases maximum deposition is shifted towards the Equator.

The similarities in the latitudinal profiles of Hg$^P$ emissions when mapped to CO and PM are reflected in their deposition profiles (Fig. 5c). The relatively greater deposition north of 60° N seen in Fig. 5c, obtained when Hg$^P$ emissions are
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Geographical distribution of the Hg<sup>P</sup>:Hg<sub>(g)</sub><sup>0</sup> emissions ratio, when mapped to PM (a) and OC (b) and when speciation is determined by FMC (c). In the colour bar the levels corresponding to the constant speciations (4:96, 15:85 and 30:70 Hg<sup>P</sup>:Hg<sub>(g)</sub><sup>0</sup>) used for the emissions mapped to CO are indicated.</figcaption></figure>

mapped to OC and when driven by FMC, reflects the peak in HgP emissions at this latitude. The greatest differences in the latitudinal deposition profiles, using the GFED inventory, are seen when varying the percentage of HgP in the emissions (Fig. 5d). Considering emissions to be solely Hg(g)0 yields a relatively smooth profile extending from pole to pole, increasing HgP causes enhanced deposition near BB hotspots. The emission peak at around 50°N remains relatively distinct also in the deposition for all the simulations (although it seen as a shoulder in the 100% Hg(g)0 profile). The peak north of 60°N is more dependent on emission speciation, supporting the previous finding that the location of Hg deposition depends on complex interactions between emission location and the time of year which influences both atmospheric trans-

port patterns and oxidant concentration fields (De Simone et al., 2015).

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Latitudinal profiles of (a) Hg(g)0 emissions when mapped to CO and when speciation is determined by FMC; (b) HgP emissions when mapped to CO, PM and OC and when speciation is driven by FMC; and (c) the relevant ratio HgP: Hg(g)0. For both Hg(g)0 and HgP emissions mapped to CO, only the speciation (15:85 HgP: Hg(g)0) is reported for clarity, whereas in (c) all the speciations are reported.</figcaption></figure>

## 3.3 Geographical distribution of Hg deposition

Due to the uncertainty in the atmospheric oxidation pathway of Hg, simulations were performed using both O<sub>3</sub> / OH and Br oxidation mechanisms to investigate their impact on Hg deposition fields. Figure 6a–d compare the geographical distribution of the modelled Hg deposition field using emission fields with 0% and of 15% HgP, for each of the oxidation mechanisms. The O<sub>3</sub> / OH mechanism leads to enhanced deposition in the tropics, whereas the Br mechanism leads to relatively higher deposition over the South Atlantic and Indian oceans. Assuming a fraction of HgP in the emissions subtracts some Hg(g)0 from the global pool, and this fraction is deposited nearer to emission sources in central Africa, South-East Asia, the Amazon and near the wildfires which
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Latitudinal profiles of the normalised Hg total deposition from the model BASE run, compared with a selection of sensitivity runs, assuming (a–b) different emission time resolution and vertical profile, as well as a combination of both; (c) different Hg<sup>P</sup> emission geographical distributions, as well as different Hg<sub>(g)</sub><sup>0</sup>: Hg<sup>P</sup> ratios. The normalisation was done by maximum.</figcaption></figure>

occur in North America and in North Asia in the northern hemispheric summer. From Fig. 6, it appears that assuming a fraction of the BB emissions to be Hg<sup>P</sup> causes the deposition field simulated using the Br oxidation mechanism to more closely resemble that using the O<sub>3</sub> / OH mechanism. To better understand the combined effect of Hg speciation and oxidation pathway on the deposition distribution, agreement maps were created to highlight the similarities and differences in the distribution of high-deposition (≥μ + 1σ, the average plus 1 standard deviation) model cells in the different simulations as described in De Simone et al. (2014). Figure 7a and b show the agreement maps of the deposition for three different Hg<sup>P</sup> fractions using the two oxidation mechanisms. Using the O<sub>3</sub> / OH mechanism, the number of model cells in which the model predicts high deposition in all three emission speciation scenarios is higher than when using the Br mechanism (631 vs. 248). This is due to the combination of high emissions and high oxidant concentrations in the tropics when using the O<sub>3</sub> / OH mechanism, constraining Hg deposition to a relatively narrow latitude band. Using the Br mechanism, Hg has a greater possibility of being transported to mid- and high latitudes before being oxidised and deposited. In both the oxidation scenarios the higher deposition over the remote areas of North America and North Asia occurs only when the fraction of Hg<sup>P</sup> in the emissions is greater than zero. High local contributions to Hg deposition from BB using the Br mechanism occur more frequently

<table><caption>**Table 2.** Characteristics of ground-based sites measuring Hg<sup>P</sup>.</caption><thead><tr><th>Long name</th><th>Short name</th><th>Lat</th><th>Long</th><th>Elev. (m)</th></tr></thead><tbody><tr><td>Amsterdam Island</td><td>AMD</td><td>-37.8</td><td>77.58</td><td>70</td></tr><tr><td>Cape Hedo</td><td>CHE</td><td>26.86</td><td>128.25</td><td>60</td></tr><tr><td>Longobucco</td><td>LON</td><td>39.39</td><td>16.61</td><td>1379</td></tr><tr><td>Manaus</td><td>MAN</td><td>-2.89</td><td>-59.97</td><td>110</td></tr><tr><td>Mauna Loa</td><td>MAU</td><td>19.54</td><td>-155.58</td><td>3399</td></tr><tr><td>Mt. Changbai</td><td>MCH</td><td>42.4</td><td>128.11</td><td>741</td></tr><tr><td>Mt. Waliguan</td><td>MWA</td><td>36.29</td><td>100.9</td><td>3816</td></tr><tr><td>Rao</td><td>RAO</td><td>57.39</td><td>11.91</td><td>5</td></tr></tbody></table>

when the fraction of Hg<sup>P</sup> is non-zero (purple in Fig. 7b), unlike the O<sub>3</sub> / OH simulations. Figure 8 contrasts the results from the two oxidation mechanisms with varying percentages of Hg<sup>P</sup> and a simulation in which the Hg<sup>P</sup> fraction was assumed to be 100%, so that it behaves as an inert tracer. The agreement maps show clearly that the similarity in the deposition fields increases with increasing Hg<sup>P</sup> fraction, reflected in the number of cells where all three simulations agree (grey in the figure) and the decrease in the number of cells where only one simulation predicts deposition higher than μ + σ (red, blue and yellow).

## 3.4 Constraints from global measurements networks

The output from the simulations including all emissions (as indicated in Table 1) for the year 2013 were compared to measurement data available from GMOS and other monitoring networks. The sites are the same as those used in Travnikov et al. (2016), the measurements from which have been reviewed Sprovieri et al. (2016a, b). Table 6 summarises a selection of metrics from the comparison for total gaseous mercury (TGM; Hg<sub>(g)</sub><sup>0</sup> + Hg<sub>(g)</sub><sup>0</sup>) and for Hg in wet deposition. The results are in line with those obtained from previous studies (De Simone et al., 2015, 2016) focusing on a different time period, and they indicate a generally good agreement between measured and simulated TGM, especially for the run with the Br-driven oxidation mechanism. For the Hg wet deposition fluxes, the results show poorer performance due to the difficulties for coarse-resolution global models to simulate precipitation events correctly (De Simone et al., 2014; Roeckner et al., 2003). Since the different sensitivity runs considering Hg<sup>P</sup> from BB differ by a only a small perturbation in the speciation of total Hg emitted from the BASE (or the relevant reference) case, the results are actually indistinguishable from BASE (or the relevant reference) case. Therefore the table reports the comparison only from runs which yield different results. Also, this means that neither wet deposition nor TGM is the most appropriate variable to assess the validity of any of the assumptions concerning Hg<sup>P</sup> emitted during BB. During 2013, within the GMOS and other Hg monitoring initiatives, a number of measurement sites collected samples of atmospheric Hg<sup>P</sup>. These stations and their precise locations are reported in the Table 2. The result of
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Geographical distribution of the Hg total deposition from model runs including only BB emission sources and assuming two different Hg$^P$ emission fractions, 15% (**a**, **c**) and 0% (**b**, **d**), for the two oxidation mechanisms considered, O$_3$ / OH (**a**–**b**) and Br (**c**–**d**).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Agreement maps of high Hg deposition model cells obtained considering only BB emissions and assuming 0, 15 and 30% to be Hg$^P$ under both the oxidation mechanisms considered, O$_3$ / OH (**a**) and Br (**b**). The maps show the areas where deposition is greater than $\mu + \sigma$.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Agreement maps, under three different speciation scenarios, 0% (a), 15% (b) and 30% (c) Hg<sup>P</sup>, of high Hg deposition model cells obtained considering only BB and using the O<sub>3</sub>/OH, the Br oxidation mechanisms, and a sensitivity run where all Hg BB emissions were considered inert (i.e. all Hg<sup>P</sup>). The deposition field from this “inert” run was retained under the three different speciation scenarios. The maps show the areas where deposition is greater than μ + σ.</figcaption></figure>

the comparison with the measurements from these sites is summarised in Fig. 9. Figure 9a shows the annually averaged surface concentrations of Hg<sup>P</sup> as simulated by the BASE run for 2013. As is evident, surface Hg<sup>P</sup> hotspots are close to the industrial areas of eastern Europe, India, East Asia and South Africa and to areas characterised by significant BB activity, including Indonesia, central Africa and boreal areas of Canada and Asia.

A first analysis to find those areas where the model run, assuming a fraction Hg<sup>P</sup> from BB (i.e. BASE), gives results that are statistically distinguishable from the model run assuming Hg from BB to be only Hg<sub>(g)</sub><sup>0</sup> was performed to identify the measurements sites best suited for further analysis.

The geographical distribution of these differences is reported in panel b of Fig. 9. The areas were the anthropogenic input is the greatest differ little between the simulations (based on a Student t test at 95% level of confidence), as
<figure><img src="image_13.png" /><figcaption>**Figure 9.** (a) Annual averaged surface Hg$^P$ concentrations as simulated by BASE run including all emission sources. (b) Differences in annual averaged surface Hg$^P$ concentrations as simulated by BASE and by NO Hg$^P$ runs, both including emissions from all sources. Black dots indicate that differences are not significant based on a Student *t* test at a 95% confidence interval. Blue bigger points indicate the locations of measurements sites reported in Table 2. Short names are depicted for sites where the differences between BASE and NO Hg$^P$ runs are significant. (c) Scatter plot of annual averaged Hg$^P$ concentrations measured at sites of Table 2 compared with those obtained by different sensitivity runs. The blue circles in the figure indicate values relative to the sites further investigated at an higher temporal resolution (see Fig. 10), whereas the red circles indicate values relative to high-altitude sites affected by processes other than BB.</figcaption></figure>

indicated by dot points in the panel. Most of the stations, depicted by the blue solid points in the same panel, are within these regions and therefore unsuitable for the analysis. Only three stations are in areas where the model results are significantly different. These, the short names of which are reported in the panel, are Amsterdam Island (AMD), Manaus (MAN)

and Mauna Loa (MAU). However, MAU and Mt. Waliguan (MWA) are high-altitude sites and affected by processes other than BB. For both the remaining stations (AMD and MAN), the fraction of Hg$^P$ that is assumed to be emitted by anthropogenic activities, as estimated by AMAP2010 inventory (AMAP/UNEP, 2013), is not sufficient alone to explain the averaged Hg$^P$ concentrations collected over the year, as is evident from Fig. 9c. The inclusion of 30 % Hg$^P$ from BB emissions at MAN and AMD and also the inclusion of 15 % Hg$^P$ from BB as using the FINN inventory at MAN significantly improve the model performances, in terms of the annual average Hg$^P$ concentrations. The result of the comparison between the Hg$^P$ concentrations collected at these two stations with the same modelled at the same points by a selection of sensitivity runs at a finer temporal resolution (daily averages) is reported in the two panels of Fig. 10. The same comparisons for all the stations, among with the box and whisker plot of distributions of the Hg$^P$ concentrations measured and modelled, are reported in Fig. 11. Although the measurement coverage of the year at MAN is sporadic, it is an important station because it is situated in a remote area where the local Hg emissions are due only to ASGM (only Hg$^0_{(g)}$) and BB (Sprovieri et al., 2016b). The consistent reduction of the error between measured and modelled Hg$^P$ concentrations when considering a fraction of particulate bound Hg emitted from BB (NRMSE from 48 to 34 % and 27 for 30 % Hg$^P$ and FINN, respectively) clearly indicates the role of BB on the observed Hg$^P$ values. At AMD (Fig. 10b), the inclusion of the fraction of Hg$^P$ from BB results only in a slightly better agreement with the measurements (NRMSE from 16 to 14 %). However, the Hg$^P$ event matching grows from 25 to 32 %, especially in the last part of the year. These Hg$^P$ events have been associated with BB events in the central Africa in Angot et al. (2014). Peaks was evaluated using the “findpeak” function in MATLAB, available from https://it.mathworks.com/help/signal/ref/findpeaks.html. To summarise, it seems that the emissions of a fraction Hg$^P$ from BB is plausible and supported by the measures of atmospheric Hg$^P$, at least for the period investigated and for the location of the two remote stations AMD and MAN. However, it has to be noted that the uncertainties related to the precise nature of atmospheric Hg$^P$ and to the processes it undergoes in the atmosphere could have an appreciable impact on the model results. For example, the assumption of a temperature-dependent gas-particle Hg$^{II}$ partitioning proposed by Amos et al. (2012) (i.e. the “Partitioning” and “Partitioning ref” runs) yield overall better model agreement with annually average Hg$^P$ concentrations (stars in Fig. 9c). However, comparing the modelled daily average time series with measurements results in clearly poorer performance at both the AMD and MAN stations (see Fig. 12b and c). More importantly, this assumption tends to render statistically indistinguishable (Student *t* test at 95 % level of confidence) the contribution of any eventual Hg$^P$ from BB, as evident from Fig. 12a.
<figure><img src="image_14.png" /><figcaption>**Figure 10.** Temporal evolution of daily averaged surface Hg<sup>P</sup> concentrations measured at Manaus (MAN) and Amsterdam Island (AMD) for the entire 2013, compared with a selection of sensitivity runs.</figcaption></figure>

**3.5 Uncertainty and biomass burning versus anthropogenic impact**

Besides the uncertainty related to the atmospheric Hg oxidation mechanism (Hynes et al., 2009; Subir et al., 2011, 2012; Gustin et al., 2015; Ariya et al., 2015) there are a number of other factors that lead to uncertainty in ascertaining the fate of Hg released by BB. Some of the model assumptions and parametrisations, in particular emission height, made little difference to the eventual deposition fields in the case where emissions from BB were considered to be 100% Hg<sup>0</sup>(g) (De Simone et al., 2015). Other sensitivity studies of the speciation of anthropogenic emissions reveal that varying the fractions of Hg<sup>II</sup>(g) and Hg<sup>P</sup> can result in quite different Hg deposition patterns due to their shorter residence time compared to Hg<sup>0</sup>(g) (De Simone et al., 2016; Bieser et al., 2014).

However, the choice of the two main vertical profile of the BB emissions used in this study, also when combined with the temporal resolution of the emissions, actually has little influence on the final Hg deposition fields. Emitting all of the Hg in a single model layer does have an impact. However, these cases are a little speculative, and therefore not included in the final analysis. The factor which has the greatest influence on the Hg deposition pattern is the choice of emis-

<figure><img src="image_15.png" /><figcaption>**Figure 11.** Left column: temporal evolution of the daily averaged surface Hg<sup>P</sup> concentrations measured at all sites from Table 2 for the entire 2013, compared with the modelled values as simulated by BASE and by NO Hg<sup>P</sup> runs, including emissions from all sources. Right column: box plots of the distribution of the daily averaged surface Hg<sup>P</sup> concentrations, for the entire 2013, as measured and simulated by the different sensitivity runs. Note the logarithmic for both MAU and MWA subplot.</figcaption></figure>

sion inventory, whereas for a given inventory the most important factors are the fraction of Hg<sup>P</sup> and the oxidation mechanism, although as seen in Sect. 3.3 the impact of the oxidation mechanism decreases with increasing Hg<sup>P</sup> fraction. The method of calculating the Hg<sup>P</sup> fraction has a limited impact on deposition on a global scale, with 66% of Hg deposited over the oceans, but the regional impact does change. Using FMC to determine the Hg<sup>P</sup> fraction increases deposition to the Arctic by 16 and 13% (O<sub>3</sub>/OH and Br) and to the Southern Ocean by 30 and 25% (O<sub>3</sub>/OH and Br); see Table 4. Apart from the polar oceans the oceanic basins, most influenced by the fraction of Hg<sup>P</sup> in the BB emissions are the North and South Pacific and the Indian ocean. The to-
<figure><img src="image_16.png" /><figcaption>**Figure 12.** (a) Differences in annual averaged surface HgP concentrations as simulated by Partitioning and by Partitioning ref. runs, both including emissions from all sources and the temperature-dependent HgII gas-particle partitioning as implemented in Amos et al. (2012). Black dots indicate that differences are not significant based on a Student *t* test at a 95 % confidence interval. Bigger blue points indicate the locations of measurements sites reported in Table 2. Temporal evolution of daily averaged surface HgP concentrations are measured at Manaus (MAN) and Amsterdam Island (AMD) for the entire 2013, compared with the modelled values from the same sensitivity runs.</figcaption></figure>

tal deposition to individual basins from the limiting 0 and 30 % HgP cases is included in Table 4. The horizontal pattern correlation method (Santer et al., 1995, 1996) and the non-parametric Kolmogorov–Smirnov two-sample test were used to assess the differences in the deposition fields obtained from the simulations summarised in Table 1, as in De Simone et al. (2015). The results of the comparison of the simulations with the BASE run are presented in Table 3. The results of the Kolmogorov–Smirnov two-sample test were exploited to construct an inspected ensemble, following the approach of Solazzo and Galmarini (2015) and previously employed in De Simone et al. (2015). The ensemble includes

<table><caption>**Table 3.** Horizontal pattern correlation (R) and probabilities that the Hg deposition fields of the different runs belong to the same distribution as the BASE run (P<sub>KS</sub>). The checks in the ensemble column indicate the inclusion of the respective run in the ensemble in Fig. 13.</caption><thead><tr><th></th><th>Sim.</th><th>R</th><th>P<sub>KS</sub></th><th>Ensemble</th></tr></thead><tbody><tr><td rowspan="4">Time resolution and vertical profile</td><td>3-hourly</td><td>1</td><td>1</td><td></td></tr><tr><td>Monthly</td><td>1</td><td>0.99</td><td></td></tr><tr><td>HAM-Profile</td><td>1</td><td>1</td><td></td></tr><tr><td>3 h + HAM-Profile</td><td>1</td><td>1</td><td></td></tr><tr><td rowspan="3">Hg<sup>P</sup> mapping</td><td>Hg<sup>P</sup> to PM</td><td>1</td><td>1</td><td></td></tr><tr><td>Hg<sup>P</sup> to OC</td><td>1</td><td>0.42</td><td>✓</td></tr><tr><td>Hg<sup>P</sup> to FMC</td><td>0.99</td><td>0.45</td><td>✓</td></tr><tr><td rowspan="3">Hg<sup>P</sup> fraction</td><td>NO Hg<sup>P</sup></td><td>0.94</td><td>0.38</td><td>✓</td></tr><tr><td>4 % Hg<sup>P</sup></td><td>0.97</td><td>0.72</td><td>✓</td></tr><tr><td>30 % Hg<sup>P</sup></td><td>0.97</td><td>0.5</td><td>✓</td></tr><tr><td rowspan="2">Inventory</td><td>GFAS</td><td>0.98</td><td>0</td><td>✓</td></tr><tr><td>FINN</td><td>0.96</td><td>0</td><td>✓</td></tr><tr><td rowspan="7">Oxidation mech. and combination</td><td>Br</td><td>0.96</td><td>0</td><td>✓</td></tr><tr><td>Br No Hg<sup>P</sup></td><td>0.81</td><td>0</td><td>✓</td></tr><tr><td>Br 30 % Hg<sup>P</sup></td><td>0.91</td><td>0</td><td>✓</td></tr><tr><td>Br Hg<sup>P</sup> to OC</td><td>0.95</td><td>0</td><td>✓</td></tr><tr><td>Br Hg<sup>P</sup> to FMC</td><td>0.94</td><td>0</td><td>✓</td></tr><tr><td>GFAS Br</td><td>0.94</td><td>0</td><td>✓</td></tr><tr><td>FINN Br</td><td>0.92</td><td>0</td><td>✓</td></tr></tbody></table>

only those simulations with realistic assumptions and deposition fields with little or no probability of belonging to the same distribution. Hg deposition from the resulting ensemble is shown in Fig. 13a. The figure shows how the inclusion of HgP in the BB emissions causes greater deposition near the hotspots of central Africa, Brazil, South-East Asia, North America and North Asia. Nonetheless approximately 70 % of Hg deposition occurs over the oceans, with the Tropical Atlantic, Tropical Pacific and Indian oceans most impacted (see Table 5). Figure 13b compares the BB ensemble results with an ensemble constructed using only anthropogenic emissions, using the EDGAR (Muntean et al., 2014), AMAP2010 (AMAP/UNEP, 2013) and STREETS (Corbitt et al., 2011) inventories (considering both oxidation mechanisms; see Table 1). It can be seen that the contribution of BB to Hg deposition is close to or greater than that from anthropogenic activities in the areas near the locations of wildfires, central Africa, the Amazon, part of the Southern Atlantic and North Asia. The contribution to Hg deposition from BB relative to anthropogenic emissions is greater than 25 % everywhere in the Southern Hemisphere and exceeds 30 % in the South Pacific and South Atlantic (Table 5). As anthropogenic Hg emissions decline the relative impact of BB Hg will rise, as shown in Fig. 14, where the Hg deposition due to BB is compared with Hg deposition from anthropogenic sources in three different emission scenarios for 2035 (see Pacyna et al., 2016, for details of the emission scenarios).
<table><caption>**Table 4.** Hg deposition (Mg) coming from BB to the oceans as obtained by the different runs for the 2013. The last two columns reports the percentage of the total Hg that deposits over sea and land.</caption><thead><tr><th rowspan="2">Run</th><th colspan="8">Total deposition/Mg</th><th colspan="2">%</th></tr><tr><th>N. Atlantic</th><th>S. Atlantic</th><th>N. Pacific</th><th>S. Pacific</th><th>Indian Ocean</th><th>Med. Sea</th><th>Arctic</th><th>S. Ocean</th><th>Sea</th><th>Land</th></tr></thead><tbody><tr><td>BASE</td><td>31.7</td><td>32.5</td><td>75.3</td><td>67.4</td><td>45.9</td><td>1.1</td><td>5.0</td><td>2.3</td><td>66</td><td>34</td></tr><tr><td>NO Hg<sup>P</sup></td><td>32.1</td><td>32.4</td><td>82.0</td><td>74.4</td><td>48.9</td><td>1.2</td><td>4.7</td><td>2.6</td><td>71</td><td>29</td></tr><tr><td>30 % Hg<sup>P</sup></td><td>31.3</td><td>32.5</td><td>69.3</td><td>61.0</td><td>43.2</td><td>1.0</td><td>5.2</td><td>2.0</td><td>62</td><td>38</td></tr><tr><td>Hg<sup>P</sup> to FMC</td><td>31.4</td><td>32.1</td><td>74.3</td><td>66.6</td><td>44.7</td><td>1.1</td><td>5.8</td><td>2.3</td><td>66</td><td>34</td></tr><tr><td>Br No Hg<sup>P</sup></td><td>26.6</td><td>39.4</td><td>75.8</td><td>83.0</td><td>55.3</td><td>1.1</td><td>3.7</td><td>7.6</td><td>74</td><td>26</td></tr><tr><td>Br 30 % Hg<sup>P</sup></td><td>28.0</td><td>36.4</td><td>61.7</td><td>61.1</td><td>44.9</td><td>0.9</td><td>4.8</td><td>4.6</td><td>62</td><td>38</td></tr><tr><td>Br Hg<sup>P</sup> to FMC</td><td>27.3</td><td>36.8</td><td>66.6</td><td>68.8</td><td>47.1</td><td>1.0</td><td>5.6</td><td>5.8</td><td>66</td><td>34</td></tr></tbody></table>

<table><caption>**Table 5.** Mercury deposition (Mg) to the oceans for 2013 from BB and comparison (ratio) with deposition from anthropogenic activities for both oxidation mechanisms.</caption><thead><tr><th>O<sub>3</sub> / OH</th><th>N. Atlantic</th><th>S. Atlantic</th><th>N. Pacific</th><th>S. Pacific</th><th>Indian Ocean</th><th>Med. Sea</th><th>Arctic</th><th>S. Ocean</th></tr></thead><tbody><tr><td>Only BB</td><td>29.8</td><td>29.9</td><td>72.1</td><td>63.0</td><td>43.0</td><td>1.1</td><td>4.7</td><td>2.1</td></tr><tr><td>Only anthropogenic</td><td>144.0</td><td>80.0</td><td>417.7</td><td>206.7</td><td>151.3</td><td>10.0</td><td>34.3</td><td>11.0</td></tr><tr><td>Ratio</td><td>0.21</td><td>0.37</td><td>0.17</td><td>0.31</td><td>0.28</td><td>0.11</td><td>0.14</td><td>0.19</td></tr><tr><td>Br</td><th>N. Atlantic</th><th>S. Atlantic</th><th>N. Pacific</th><th>S. Pacific</th><th>Indian Ocean</th><th>Med. Sea</th><th>Arctic</th><th>S. Ocean</th></tr><tr><td>Only BB</td><td>25.7</td><td>34.7</td><td>65.1</td><td>66.2</td><td>46.2</td><td>0.9</td><td>4.2</td><td>5.1</td></tr><tr><td>Only anthropogenic</td><td>153</td><td>85.33</td><td>457.3</td><td>188.3</td><td>140</td><td>12.33</td><td>34</td><td>27.3</td></tr><tr><td>Ratio</td><td>0.17</td><td>0.41</td><td>0.14</td><td>0.35</td><td>0.33</td><td>0.08</td><td>0.12</td><td>0.19</td></tr></tbody></table>

## 4 Conclusions

That a fraction of HgP is present in BB Hg emissions has been confirmed by several field measurements (Obrist et al., 2007; Finley et al., 2009), and this fact has been suggested as an explanation of high HgP observations at a remote site (Angot et al., 2014), but this is the first time it has been included in a model study to assess its effects on a global scale. A previous modelling study assuming emissions from BB to be 100 % Hg(g) (De Simone et al., 2015) suggested that as much as 75 % of the Hg emitted by BB was deposited to ocean basins, with global implications for food webs and human health. Including a fraction of HgP in the BB Hg emissions has an impact on the geographical distribution of the deposition fluxes for the year analysed, reducing input to the global oceans and some high-latitude regions, while enhancing potentially negative effects on ecosystems close to areas where significant BB occurs. The presence of HgP in the emissions decreases the differences seen in Hg deposition patterns produced by employing different oxidation mechanisms. In the remote areas of North Asia and North America, BB has a strong local impact if the HgP fraction is non-zero. This latter result is independent of the atmospheric oxidation pathway. In simulations with 30 % HgP in the BB emissions, deposition over the Arctic increases by 11 % with respect to 0 % HgP (30 % in the Br simulations) and by 16 % when the HgP fraction is determined by FMC (37 % in the Br simulation).

The fraction of HgP released from BB while having an impact on the land-sea distribution of global Hg deposition, has a more significant impact in particular regions including the polar regions, the South Atlantic and Pacific and Indian oceans. These results apply for the investigated year (2013) and may differ for other years due to the complex interaction of the numerous factors determining the final fate of Hg. However, few alternatives of analysis period exist due to the limited time coverage of global measurement network(s). Indeed the year selected for the analysis allowed for the hypotheses tested in this study to be supported by observations at a number of sites from GMOS, which has extended the observational network in the tropics and the Southern Hemisphere (Sprovieri et al., 2016a, b). The eventual emissions of a fraction of HgP from BB cannot be evaluated by comparison with observed gaseous atmospheric Hg concentrations or Hg in wet precipitation samples due to the very small impact of HgP from BB on both the atmospheric burden and wet deposition relative to all other emissions sources (≈ 1–2 %). Conversely, its contribution to atmospheric HgP is comparable to that of anthropogenic activities and therefore may be investigated. The inclusion in the model run of a fraction of HgP from BB contributes to better model performances at two remote sites, Manaus and Amsterdam Island. However results are not definitive due to the large uncertainty related to HgP emissions and transformation processes. Further modelling and more measurement sites, particularly in remote ar-
<table><caption>**Table 6.** Comparison of the results of BASE and Br simulations including all emissions sources with observations from measurement networks for 2013.</caption><thead><tr><th rowspan="2"></th><th colspan="4">Total gaseous mercury</th><th colspan="4">Wet deposition</th></tr><tr><th colspan="2">Regression</th><th rowspan="2">r</th><th rowspan="2">Stats NRMSE %</th><th colspan="2">Regression</th><th rowspan="2">r</th><th rowspan="2">Stats NRMSE %</th></tr><tr><th></th><th>Intercept</th><th>Slope</th><th>Intercept</th><th>Slope</th></tr></thead><tbody><tr><td>BASE</td><td>0.36</td><td>0.62</td><td>0.72</td><td>10.54</td><td>5.84</td><td>0.04</td><td>0.12</td><td>6.89</td></tr><tr><td>Partitioning</td><td>0.34</td><td>0.7</td><td>0.73</td><td>11.9</td><td>3.71</td><td>0.03</td><td>0.14</td><td>4.76</td></tr><tr><td>Br</td><td>-0.08</td><td>0.96</td><td>0.74</td><td>15.68</td><td>7.1</td><td>0.08</td><td>0.18</td><td>9.12</td></tr></tbody></table>

<figure><img src="image_17.png" /><figcaption>**Figure 13.** Geographical distribution of the total Hg deposition from BB emissions obtained from an ensemble of simulations for the year 2013 (a) in terms of the average (μ) and standard deviation σ of the ensemble. The comparison of the BB simulation with an ensemble of runs including only anthropogenic emissions (De Simone et al., 2016) shows (b) the geographic distribution of the fraction of the BB contribution to the Hg deposition from the anthropogenic sources.</figcaption></figure>

eas, would help reduce some of the uncertainties associated with Hg emissions from BB and constrain these processes. Biomass burning has and will continue to play a significant role in the cycling of legacy Hg, and its relative importance is likely to increase as anthropogenic emissions are reduced and global temperatures rise.

<figure><img src="image_18.png" /><figcaption>**Figure 14.** Ratio of the Hg deposition due to biomass burning with respect to Hg deposition due to anthropogenic emissions for three anthropogenic emissions scenarios for 2035: (a) current policy (CP), (b) new policy (NP) and (c) maximum feasible reduction (MFR).</figcaption></figure>

## 5 Data availability

Mercury data discussed in this paper are reported within the GMOS central database and are available upon request at http://sdi.iia.cnr.it/geoint/publicpage/GMOS/gmos_ historical.zul.
# Appendix A: How Hg emission fields are calculated

## A1 Mapping to CO

When mapped to CO, the emissions of Hg<sub>(g)</sub><sup>0</sup> were calculated from those of CO using a global averaged ER (1.96 × 10<sup>-7</sup>). These were unchanged in the run assuming Hg emissions from BB to be 100% Hg<sub>(g)</sub><sup>0</sup> and divided between Hg<sup>P</sup> and Hg<sub>(g)</sub><sup>0</sup> species, with ratios 4:96, 15:85 and 30:70, in mass, in the runs considering the respective constant fractions of Hg<sup>P</sup>. Consequently, the geographical and temporal distributions of Hg<sub>(g)</sub><sup>0</sup> and Hg<sup>P</sup> BB emissions follow those of CO. For all cases, the GFEDv4 inventory was used, except for those sensitivity runs performed to test the impact of different inventories, FINNv1.5 and GFAS1.4.

## A2 Mapping to OC

When mapped to OC, geographical and temporal distributions of Hg<sub>(g)</sub><sup>0</sup> BB emissions, as well as the total Hg emitted, were calculated as described in Appendix A1. The fractioning of Hg emissions, in mass, between Hg<sup>P</sup> and Hg<sub>(g)</sub><sup>0</sup> species was assumed to be in the ratio 15:85. The Hg<sup>P</sup> emissions so calculated were then geographically and temporally mapped to those of OC from the GFEDv4 inventory.

## A3 Mapping to PM

This mapping method is similar to the one described in Appendix A2, except for the fact the Hg<sup>P</sup> temporal and geographical distributions follow those of PM from the GFEDv4 inventory.

## A4 Emissions speciation determination by FMC

When using this procedure for determining the BB emission speciation between Hg<sub>(g)</sub><sup>0</sup> and Hg<sup>P</sup>, the geographical and temporal distributions of Hg<sub>(g)</sub><sup>0</sup> and Hg<sup>P</sup> BB emissions and the total Hg emitted were calculated in the same way as described in Appendix A1. The main difference is in that the fractioning of Hg emissions, in mass, between Hg<sub>(g)</sub><sup>0</sup> and Hg<sup>P</sup> species were calculated dynamically using the piece wise linear relationship between fuel moisture content empirically determined by relative figure in Obrist et al. (2007). As a proxy for FMC, we used the monthly averaged vegetation water content derived from passive microwave remote sensing data (Advanced Microwave Scanning Radiometer 2 (ASMR2)), and employing the Land Parameter Retrieval Model (LPRM) available at http://gcmd.nasa.gov/search/Metadata.do?Entry=C1235316240-GES_DISC#metadata.
*Competing interests.* The authors declare that they have no conflict of interest.

*Disclaimer.* The US Environmental Protection Agency (EPA) through its Office of Research and Development participated in this research. The views expressed in this paper are those of the authors and do not necessarily reflect the views or policies of EPA. It has been subjected to EPA Agency review and approved for publication. Mention of trade names or commercial products does not constitute endorsement or recommendation for use.

*Acknowledgements.* We are grateful to Sebastian Rast and colleagues at the Max Planck Institute for Meteorology in Hamburg, Germany, for the distribution of their software ECHAM5 and for providing the access to the processed ERA-INTERIM data. We are grateful to Xin Yang for providing the Br/BrO fields from p-Tomcat. The research was performed in the framework of the EU project GMOS (FP7-265113). The authors would also like to thank the referees whose helpful suggestions and comments contributed much to improving the original manuscript.

Edited by: R. Ebinghaus

Reviewed by: M. Cohen, T. Dvonch, and two anonymous referees

## References

Akagi, S. K., Yokelson, R. J., Wiedinmyer, C., Alvarado, M. J., Reid, J. S., Karl, T., Crounse, J. D., and Wennberg, P. O.: Emission factors for open and domestic biomass burning for use in atmospheric models, Atmos. Chem. Phys., 11, 4039–4072, doi:10.5194/acp-11-4039-2011, 2011.

AMAP/UNEP: Technical Background Report for the Global Mercury Assessment 2013, Tech. rep., Arcitic Monitoring and Assessment Programme, Oslo, Norway/UNEP ChemicalsBranch, Geneva, Switzerland, available at: http://www.amap.no/documents/doc/technical-background-report-for-the-global-mercury-assessment-2013/ 848 (last access: 1 February 2017), 2013.

Amos, H. M., Jacob, D. J., Holmes, C. D., Fisher, J. A., Wang, Q., Yantosca, R. M., Corbitt, E. S., Galarneau, E., Rutter, A. P., Gustin, M. S., Steffen, A., Schauer, J. J., Graydon, J. A., Louis, V. L. St., Talbot, R. W., Edgerton, E. S., Zhang, Y., and Sunderland, E. M.: Gas-particle partitioning of atmospheric Hg(II) and its effect on global mercury deposition, Atmos. Chem. Phys., 12, 591–603, doi:10.5194/acp-12-591-2012, 2012.

Andela, N., Kaiser, J., Heil, A., van Leeuwen, T., van der Werf, G., Wooster, M., Remy, S., and Schultz, M.: Assessment of the Global Fire Assimilation System (GFASv1), MACC-II (Monitoring Atmospheric Composition and Climate) project, available at: http://juser.fz-juelich.de/record/186645 (last access: 5 May 2016), 2013.

Andersson, M. E., Sommar, J., Gårdfeldt, K., and Jutterström, S.: Air–sea exchange of volatile mercury in the North Atlantic Ocean, Mar. Chem., 125, 1–7, 2011.

Angot, H., Barret, M., Magand, O., Ramonet, M., and Dommergue, A.: A 2-year record of atmospheric mercury species at a background Southern Hemisphere station on Amsterdam Island, Atmos. Chem. Phys., 14, 11461–11473, doi:10.5194/acp-14-11461-2014, 2014.

Ariya, P. A., Amyot, M., Dastoor, A., Deeds, D., Feinberg, A., Kos, G., Poulain, A., Ryjkov, A., Semeniuk, K., Subir, M., and Toyota, K.: Mercury Physicochemical and Biogeochemical Transformation in the Atmosphere and at Atmospheric Interfaces: A Review and Future Directions, Chem. Rev., 115, 3760–3802, 2015.

Bieser, J., De Simone, F., Gencarelli, C., Geyer, B., Hedgecock, I., Matthias, V., Travnikov, O., and Weigelt, A.: A diagnostic evaluation of modeled mercury wet depositions in Europe using atmospheric speciated high-resolution observations, Environ. Sci. Pollut. R., 21, 9995–10012, 2014.

Carbone, F., Landis, M., Gencarelli, C., Naccarato, A., Sprovieri, F., Simone, F. D., Hedgecock, I., and Pirrone, N.: Sea Surface Temperature variation linked to elemental mercury concentrations measured on Mauna Loa, Geophys. Res. Lett., 43, 7751–7757, doi:10.1002/2016GL069252, 2016.

Chen, C. Y., Driscoll, C. T., Lambert, K. F., Mason, R. P., and Sunderland, E. M.: Connecting mercury science to policy: from sources to seafood, Rev. Environ. Health, 31, 17–20, 2016.

Corbitt, E. S., Jacob, D. J., Holmes, C. D., Streets, D. G., and Sunderland, E. M.: Global source–receptor relationships for mercury deposition under present-day and 2050 emissions scenarios, Environ. Sci. Technol., 45, 10477–10484, 2011.

D'Amore, F., Bencardino, M., Cinnirella, S., Sprovieri, F., and Pirrone, N.: Data quality through a web-based QA/QC system: implementation for atmospheric mercury data from the global mercury observation system, Environmental Science: Processes & Impacts, 17, 1482–1491, 2015.

De Simone, F., Gencarelli, C., Hedgecock, I., and Pirrone, N.: Global atmospheric cycle of mercury: a model study on the impact of oxidation mechanisms, Environ. Sci. Pollut. R., 21, 4110–4123, 2014.

De Simone, F., Cinnirella, S., Gencarelli, C. N., Yang, X., Hedgecock, I. M., and Pirrone, N.: Model study of global mercury deposition from biomass burning, Environ. Sci. Technol., 49, 6712–6721, 2015.

De Simone, F., Gencarelli, C. N., Hedgecock, I. M., and Pirrone, N.: A Modeling Comparison of Mercury Deposition from Current Anthropogenic Mercury Emission Inventories, Environ. Sci. Technol., 50, 5154–5162, 2016.

Ebinghaus, R., Slemr, F., Brenninkmeijer, C. A. M., van Velthoven, P., Zahn, A., Hermann, M., O'Sullivan, D. A., and Oram, D. E.: Emissions of gaseous mercury from biomass burning in South America in 2005 observed during CARIBIC flights, Geophys. Res. Lett., 34, L08813, doi:10.1029/2006GL028866, 2007.

Emmons, L. K., Walters, S., Hess, P. G., Lamarque, J.-F., Pfister, G. G., Fillmore, D., Granier, C., Guenther, A., Kinnison, D., Laepple, T., Orlando, J., Tie, X., Tyndall, G., Wiedinmyer, C., Baughcum, S. L., and Kloster, S.: Description and evaluation of the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4), Geosci. Model Dev., 3, 43–67, doi:10.5194/gmd-3-43-2010, 2010.

Feng, X., Li, P., Qiu, G., Wang, S., Li, G., Shang, L., Meng, B., Jiang, H., Bai, W., Li, Z., and Fu, X.: Human Exposure To Methylmercury through Rice Intake in Mercury Mining Areas,
Guizhou Province, China, Environ. Sci. Technol., 42, 326–332, 2008.

Finley, B., Swartzendruber, P., and Jaffe, D.: Particulate mercury emissions in regional wildfire plumes observed at the Mount Bachelor Observatory, Atmos. Environ., 43, 6074–6083, 2009.

Friedli, H., Arellano, A., Cinnirella, S., and Pirrone, N.: Initial Estimates of Mercury Emissions to the Atmosphere from Global Biomass Burning, Environ. Sci. Technol., 43, 3507–3513, 2009.

GFAS: Global Fire Assimilation System v1.2 documentation, available at: http://www.gmes-atmosphere.eu/about/project_structure/input_data/d_fire/gfas_versions (last access: 1 June 2016), 2015.

Giglio, L., Randerson, J. T., and Werf, G. R.: Analysis of daily, monthly, and annual burned area using the fourth-generation global fire emissions database (GFED4), J. Geophys. Res.-Biogeo., 118, 317–328, 2013.

Gustin, M. and Jaffe, D.: Reducing the uncertainty in measurement and understanding of mercury in the atmosphere, Environ. Sci. Technol., 44, 2222–2227, 2010.

Gustin, M. S., Amos, H. M., Huang, J., Miller, M. B., and Heidecorn, K.: Measuring and modeling mercury in the atmosphere: a critical review, Atmos. Chem. Phys., 15, 5697–5713, doi:10.5194/acp-15-5697-2015, 2015.

Holmes, C. D., Jacob, D. J., Corbitt, E. S., Mao, J., Yang, X., Talbot, R., and Slemr, F.: Global atmospheric model for mercury including oxidation by bromine atoms, Atmos. Chem. Phys., 10, 12037–12057, doi:10.5194/acp-10-12037-2010, 2010.

Hynes, A. J., Donohoue, D. L., Goodsite, M. E., and Hedgecock, I. M.: Our current understanding of major chemical and physical processes affecting mercury dynamics in the atmosphere and at the air-water/terrestrial interfaces, in: Mercury Fate and Transport in the Global Atmosphere: Emissions, Measurements and Models, edited by: Pirrone, N. and Mason, R. P., Chap. 14, 427–457, Springer, 2009.

Janhäll, S., Andreae, M. O., and Pöschl, U.: Biomass burning aerosol emissions from vegetation fires: particle number and mass emission factors and size distributions, Atmos. Chem. Phys., 10, 1427–1439, doi:10.5194/acp-10-1427-2010, 2010.

Jung, G., Hedgecock, I. M., and Pirrone, N.: ECHMERIT V1.0 – a new global fully coupled mercury-chemistry and transport model, Geosci. Model Dev., 2, 175–195, doi:10.5194/gmd-2-175-2009, 2009.

Kaiser, J. W., Heil, A., Andreae, M. O., Benedetti, A., Chubarova, N., Jones, L., Morcrette, J.-J., Razinger, M., Schultz, M. G., Suttie, M., and van der Werf, G. R.: Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power, Biogeosciences, 9, 527–554, doi:10.5194/bg-9-527-2012, 2012.

Kaiser, J. W., Andela, N., Benedetti, A., He, J., Heil, A., Inness, A., Paugam, R., Remy, S., Trigo, I., van der Werf, G. R., and Wooster, M. J.: GFAS fire emissions and smoke in the Copernicus Atmosphere Monitoring Service, Geophys. Res. Abstr., EGU2015-15432, EGU General Assembly 2015, Vienna, Austria, 2015.

Kwon, S. Y. and Selin, N. E.: Uncertainties in Atmospheric Mercury Modeling for Policy Evaluation, Current Pollution Reports, 2, 103–114, 2016.

Meng, B., Feng, X., Qiu, G., Anderson, C. W., Wang, J., and Zhao, L.: Localization and speciation of mercury in brown rice with

implications for Pan-Asian public health, Environ. Sci. Technol., 48, 7974–7981, 2014.

Mu, M., Randerson, J. T., van der Werf, G. R., Giglio, L., Kasibhatla, P., Morton, D., Collatz, G. J., DeFries, R. S., Hyer, E. J., Prins, E. M., Griffith, D. W. T., Wunch, D., Toon, G. C., Sherlock, V., and Wennberg, P. O.: Daily and 3-hourly variability in global fire emissions and consequences for atmospheric model predictions of carbon monoxide, J. Geophys. Res.-Atmos., 116, D24303, doi:10.1029/2011JD016245, 2011.

Muntean, M., Janssens-Maenhout, G., Song, S., Selin, N. E., Olivier, J. G., Guizzardi, D., Maas, R., and Dentener, F.: Trend analysis from 1970 to 2008 and model evaluation of EDGARv4 global gridded anthropogenic mercury emissions, Sci. Total Environment, 494–495, 337–350, 2014.

Obrist, D., Moosmüller, H., Schürmann, R., Chen, L.-W. A., and Kreidenweis, S. M.: Particulate-phase and gaseous elemental mercury emissions during biomass combustion: controlling factors and correlation with particulate matter emissions, Environ. Sci. Technol., 42, 721–727, 2007.

Pacyna, J. M., Travnikov, O., De Simone, F., Hedgecock, I. M., Sundseth, K., Pacyna, E. G., Steenhuisen, F., Pirrone, N., Munthe, J., and Kindbom, K.: Current and future levels of mercury atmospheric pollution on a global scale, Atmos. Chem. Phys., 16, 12495–12511, doi:10.5194/acp-16-12495-2016, 2016.

Randerson, J., Chen, Y., Werf, G., Rogers, B., and Morton, D.: Global burned area and biomass burning emissions from small fires, J. Geophys. Res.-Biogeo., 117, G04012, doi:10.1029/2012JG002128, 2012.

Roeckner, E., Bäuml, G., Bonaventura, L., Brokopf, R., Esch, M., Giorgetta, M., Hagemann, S., Kirchner, I., Kornblueh, L., Manzini, E., Rhodin, A., Schlese, U., Schulzweida, U., and Tompkins, A.: The atmospheric general circulation model ECHAM5. PART I: Model description, MPI-Report No. 349, available at: https://www.mpimet.mpg.de/fileadmin/publikationen/Reports/max_scirep_349.pdf (last access: 5 May 2016), 2003.

Santer, B. D., Taylor, K. E., Wigley, T. M., Penner, J. E., Jones, P. D., and Cubasch, U.: Towards the detection and attribution of an anthropogenic effect on climate, Clim. Dynam., 12, 77–100, 1995.

Santer, B. D., Taylor, K. E., Wigley, T. M. L., Johns, T. C., Jones, P. D., Karoly, D. J., Mitchell, J. F. B., Oort, A. H., Penner, J. E., Ramaswamy, V., Schwarzkopf, M. D., Stouffer, R. J., and Tett, S.: A search for human influences on the thermal structure of the atmosphere, Nature, 382, 39–46, 1996.

Selin, N. E., Jacob, D. J., Yantosca, R. M., Strode, S., Jaeglé, L., and Sunderland, E. M.: Global 3-D land-ocean-atmosphere model for mercury: Present-day versus preindustrial cycles and anthropogenic enrichment factors for deposition, Global Biogeochem. Cy., 22, GB2011, doi:10.1029/2007GB003040, 2008.

Slemr, F., Brenninkmeijer, C. A., Rauthe-Schöch, A., Weigelt, A., Ebinghaus, R., Brunke, E.-G., Martin, L., Spain, T. G., and O’Doherty, S.: El Niño–Southern Oscillation influence on tropospheric mercury concentrations, Geophys. Res. Lett., 43, 1766–1771, 2016.

Solazzo, E. and Galmarini, S.: A science-based use of ensembles of opportunities for assessment and scenario studies, Atmos. Chem. Phys., 15, 2535–2544, doi:10.5194/acp-15-2535-2015, 2015.
Sprovieri, F., Pirrone, N., Bencardino, M., D'Amore, F., Angot, H., Barbante, C., Brunke, E.-G., Arcega-Cabrera, F., Cairns, W., Comero, S., Diéguez, M. D. C., Dommergue, A., Ebinghaus, R., Feng, X. B., Fu, X., Garcia, P. E., Gawlik, B. M., Hageström, U., Hansson, K., Horvat, M., Kotnik, J., Labuschagne, C., Magand, O., Martin, L., Mashyanov, N., Mkololo, T., Munthe, J., Obolkin, V., Islas, M. R., Sena, F., Somerset, V., Spandow, P., Vardè, M., Walters, C., Wängberg, I., Weigelt, A., Yang, X., and Zhang, H.: Five-year records of Total Mercury Deposition flux at GMOS sites in the Northern and Southern Hemispheres, Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-517, in review, 2016a.

Sprovieri, F., Pirrone, N., Bencardino, M., D'Amore, F., Carboni, F., Cinnirella, S., Mannarino, V., Landis, M., Ebinghaus, R., Weigelt, A., Brunke, E.-G., Labuschagne, C., Martin, L., Munthe, J., Wängberg, I., Artaxo, P., Morais, F., Barbosa, H. D. M. J., Brito, J., Cairns, W., Barbante, C., Diéguez, M. D. C., Garcia, P. E., Dommergue, A., Angot, H., Magand, O., Skov, H., Horvat, M., Kotnik, J., Read, K. A., Neves, L. M., Gawlik, B. M., Sena, F., Mashyanov, N., Obolkin, V., Wip, D., Feng, X. B., Zhang, H., Fu, X., Ramachandran, R., Cossa, D., Knoery, J., Marusczak, N., Nerentorp, M., and Norstrom, C.: Atmospheric mercury concentrations observed at ground-based monitoring sites globally distributed in the framework of the GMOS network, Atmos. Chem. Phys., 16, 11915–11935, doi:10.5194/acp-16-11915-2016, 2016b.

Steffen, A., Bottenheim, J., Cole, A., Ebinghaus, R., Lawson, G., and Leaitch, W. R.: Atmospheric mercury speciation and mercury in snow over time at Alert, Canada, Atmos. Chem. Phys., 14, 2219–2231, doi:10.5194/acp-14-2219-2014, 2014.

Subir, M., Ariya, P. A., and Dastoor, A. P.: A review of uncertainties in atmospheric modeling of mercury chemistry I. Uncertainties in existing kinetic parameters – Fundamental limitations and the importance of heterogeneous chemistry, Atmos. Environ., 45, 5664–5676, 2011.

Subir, M., Ariya, P. A., and Dastoor, A. P.: A review of the sources of uncertainties in atmospheric mercury modeling II. Mercury surface and heterogeneous chemistry – A missing link, Atmos. Environ., 46, 1–10, 2012.

Travnikov, O., Angot, H., Artaxo, P., Bencardino, M., Bieser, J., D'Amore, F., Dastoor, A., De Simone, F., Diéguez, M. D. C., Dommergue, A., Ebinghaus, R., Feng, X. B., Gencarelli, C. N., Hedgecock, I. M., Magand, O., Martin, L., Matthias, V., Mashyanov, N., Pirrone, N., Ramachandran, R., Read, K. A., Ryjkov, A., Selin, N. E., Sena, F., Song, S., Sprovieri, F., Wip, D., Wängberg, I., and Yang, X.: Multi-model study of mercury dispersion in the atmosphere: Atmospheric processes and model evaluation, Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-924, in review, 2016.

van der Werf, G. R., Randerson, J. T., Giglio, L., Collatz, G. J., Mu, M., Kasibhatla, P. S., Morton, D. C., DeFries, R. S., Jin, Y., and van Leeuwen, T. T.: Global fire emissions and the contribution of deforestation, savanna, forest, agricultural, and peat fires (1997–2009), Atmos. Chem. Phys., 10, 11707–11735, doi:10.5194/acp-10-11707-2010, 2010.

Veira, A., Kloster, S., Wilkenskjeld, S., and Remy, S.: Fire emission heights in the climate system – Part 1: Global plume height patterns simulated by ECHAM6-HAM2, Atmos. Chem. Phys., 15, 7155–7171, doi:10.5194/acp-15-7155-2015, 2015.

Wang, X., Zhang, H., Lin, C.-J., Fu, X., Zhang, Y., and Feng, X.: Transboundary transport and deposition of Hg emission from springtime biomass burning in the Indo-China Peninsula, J. Geophys. Res.-Atmos., 120, 9758–9771, 2015.

Wang, Y., Huang, J., Zananski, T. J., Hopke, P. K., and Holsen, T. M.: Impacts of the Canadian forest fires on atmospheric mercury and carbonaceous particles in northern New York, Environ. Sci. Technol., 44, 8435–8440, 2010.

Webster, J. P., Kane, T. J., Obrist, D., Ryan, J. N., and Aiken, G. R.: Estimating mercury emissions resulting from wildfire in forests of the Western United States, Sci. Total Environ., 568, 578–586, doi:10.1016/j.scitotenv.2016.01.166, 2016.

Wiedinmyer, C., Akagi, S. K., Yokelson, R. J., Emmons, L. K., Al-Saadi, J. A., Orlando, J. J., and Soja, A. J.: The Fire INventory from NCAR (FINN): a high resolution global model to estimate the emissions from open burning, Geosci. Model Dev., 4, 625–641, doi:10.5194/gmd-4-625-2011, 2011.

Witt, E. L., Kolka, R. K., Nater, E. A., and Wickman, T. R.: Forest fire effects on mercury deposition in the boreal forest, Environ. Sci. Technol., 43, 1776–1782, 2009.

Zhang, H., Feng, X., Larssen, T., Shang, L., and Li, P.: Bioaccumulation of methylmercury versus inorganic mercury in rice (*Oryza sativa* L.) grain, Environ. Sci. Technol., 44, 4499–4504, 2010.

Zhang, K., O'Donnell, D., Kazil, J., Stier, P., Kinne, S., Lohmann, U., Ferrachat, S., Croft, B., Quaas, J., Wan, H., Rast, S., and Feichter, J.: The global aerosol-climate model ECHAM-HAM, version 2: sensitivity to improvements in process representations, Atmos. Chem. Phys., 12, 8911–8949, doi:10.5194/acp-12-8911-2012, 2012.

Zhang, Y., Obrist, D., Zielinska, B., and Gertler, A.: Particulate emissions from different types of biomass burning, Atmos. Environ., 72, 27–35, 2013.