Iron from coal combustion particles dissolves much faster than mineral dust under simulated atmospheric acidic conditions

. Mineral dust is the largest source of aerosol iron (Fe) to the offshore global ocean, but acidic processing of coal ﬂy ash (CFA) in the atmosphere could be an important source of soluble aerosol Fe. Here, we determined the Fe speciation and dissolution kinetics of CFA from Aberthaw (United Kingdom), Krakow (Poland), and Shandong (China) in solutions which simulate atmospheric acidic processing. In CFA PM 10 fractions, 8 %–21.5 % of the total Fe was found to be hematite and goethite kinetics Fe in CFA up to 7 times a Saharan dust sample Based on these laboratory data, we developed a new scheme for the proton- and oxalate-promoted Fe dissolution of CFA, which was implemented into the global atmospheric chemical transport model IMPACT (Integrated Massively Parallel Atmospheric Chemical Transport). The revised model showed a better agreement with observations of Fe solubility in aerosol particles over the Bay of due to the initial rapid release of Fe and the suppression of the oxalate-promoted dissolution at low pH. The improved model enabled us to predict sensitivity to a more dynamic range of pH changes, between anthropogenic combustion and aerosols.


Introduction
The availability of iron (Fe) limits primary productivity in high-nutrient, low-chlorophyll (HNLC) regions of the global ocean, including the subarctic North Pacific, the eastern Equatorial Pacific, and the Southern Ocean (Boyd et al., 2007;Martin, 1990). In other regions of the global ocean, such as the subtropical North Atlantic, the Fe input may affect primary productivity by stimulating nitrogen fixation (Mills et al., 2004;Moore et al., 2006). These areas are particularly sensitive to changes in the supply of bioavailable Fe. Atmospheric aerosols are an important source of soluble (and, thus, potentially bio-accessible) Fe to the offshore global ocean. The deposition of bio-accessible Fe to the ocean can alter biogeochemical cycles and increase the carbon uptake, consequently affecting the climate (e.g. Jickells and Moore, 2015;Jickells et al., 2005;Kanakidou et al., 2018;Mahowald et al., 2010;Shi et al., 2012). In general, bio-accessible Fe consists of aerosol-dissolved Fe, and Fe nanoparticles which can be present in the original particulate matter and/or formed during atmospheric transport as a result of cycling into and out of clouds (Shi et al., 2009). It is, in addition, possible that other more refractory forms of Fe could be solubilized in the surface waters by zooplankton (Schlosser et al., 2018) or the microbial community (Rubin et al., 2011).
The Fe transported in the atmosphere is largely derived from lithogenic sources, which contribute around 95 % of the total Fe in suspended particles (e.g. Shelley et al., 2018), and most studies so far have concentrated on the atmospheric processing of mineral dust (e.g. Cwiertny et al., 2008;Fu et al., 2010;Ito and Shi, 2016;Shi et al., 2011aShi et al., , 2015. Mineral dust has low Fe solubility (dissolved Fe / total Fe × 100) near the source regions, generally below 1 % (e.g. Shi et al., 2011;Sholkovitz et al., 2009Sholkovitz et al., , 2012, increasing somewhat as a result of processes occurring during atmospheric transport (e.g. Baker et al., 2021Baker et al., , 2020. Other sources of bio-accessible Fe to the ocean are from combustion sources such as biomass burning, coal combustion, oil combustion, and metal smelting (e.g. Ito et al., 2018;Rathod et al., 2020). Although these sources are only a small fraction of the total Fe in atmospheric particulates, the Fe solubility of pyrogenic sources can be 1-2 orders of magnitude higher than in mineral dust (Ito et al., 2021b and references therein) and, thus, can be important in promoting carbon uptake. However, the Fe solubility of pyrogenic sources varies considerably, depending on the particular sources, with higher values observed for oil combustion and biomass burning than coal combustion sources (Ito et al., 2021b and references therein). Wang et al. (2015) estimated that coal combustion emitted around ∼ 0.9 Tg yr −1 of Fe into the atmosphere (on average for , contributing up to ∼ 86 % of the total anthropogenic Fe emissions. A more recent study, which has included metal smelting as an atmospheric Fe source, estimated that coal combustion emitted ∼ 0.7 Tg yr −1 of Fe for the year 2010, contributing around 34 % of the total anthropogenic Fe atmospheric loading (Rathod et al., 2020). Although the use of coal as a principal energy source has been recently reduced as a result of concerns about air quality and global warming, coal is still an important energy source in a number of countries, in particular in the Asia-Pacific region (BP, 2020). In China, most of the total energy is supplied by coal, contributing over 50 % of the global coal consumption in 2019, followed by India (12 %), and the USA (8 %). Germany and Poland are the largest coal consumers in Europe, accounting together for around 40 % of the European usage (BP, 2020). South Africa is also among the principal countries for coal consumption (BP, 2020) and is a source of Fe-bearing particles to the anaemic Southern Ocean (e.g. Ito et al., 2019).
Coal fly ash (CFA) is a by-product of coal combustion. This generally consists of glassy spherical particles (e.g. Brown et al., 2011) which are formed through different transformations (decomposition, fusion, agglomeration, and volatilization) of mineral matter in coal during combustion (e.g. Jones, 1995) and are transported with the flue gases undergoing rapid solidification. CFA is co-emitted with acidic gases such as sulfur dioxide (SO 2 ), nitrogen oxides (NO x ), and carbon dioxide (CO 2 ; e.g. Munawer, 2018).
During long-range transport, CFA particles undergo atmospheric processing with the CFA surface coated by acidic species, such as sulfuric acid (H 2 SO 4 ) and oxalic acid (H 2 C 2 O 4 ), in atmospheric aerosols. Aged CFA particles are hygroscopic and absorb water at typical relative humidity in the marine atmosphere. As a result, a thin layer of water with high acidity, low pH, and high ionic strength is formed around the particles (Meskhidze et al., 2003;Spokes and Jickells, 1995;Zhu et al., 1992). In addition, ammonia (NH 3 ), which is a highly hydrophilic gas, can also partition into the aerosol phase, react with H 2 SO 4 , and form ammonium sulfate ((NH 4 ) 2 SO 4 ), which is an important inorganic salt contributing to the high ionic strength in aged atmospheric aerosols (Seinfeld and Pandis, 2016).
At low pH conditions, Fe solubility in aerosols increases, as the high concentration of protons (H + ) weakens the Fe-O bonds facilitating the detachment of Fe from the surface lattice (Furrer and Stumm, 1986). Li et al. (2017) provided the first observational evidence that acidification leads to the release of Fe from anthropogenic particles.
In addition to these inorganic processes, organic ligands can also enhance atmospheric Fe dissolution by forming soluble complexes with Fe (e.g. Cornell and Schwertmann, 2003). For example, H 2 C 2 O 4 is an important organic species in aerosols (e.g. Kawamura and Bikkina, 2016). Laboratory studies have demonstrated that H 2 C 2 O 4 increases Fe solubility of aerosol sources (Chen and Grassian, 2013;Johnson and Meskhidze, 2013;Paris and Desboeufs, 2013;Paris et al., 2011;Xu and Gao, 2008). Recently, observations over the Bay of Bengal indicate that H 2 C 2 O 4 contributes to the increase in dissolved Fe in atmospheric water (Bikkina et al., 2020).
To simulate the Fe dissolution in CFA, it is necessary to determine the dissolution kinetics under realistic conditions. Previous studies have investigated the Fe dissolution kinetics of CFA under acidic conditions. Chen et al. (2012) simulated acidic and cloud processing of certified CFA. Fu et al. (2012) determined the dissolution kinetics of CFA samples at pH 2, while Chen and Grassian (2013) investigated the effect of organic species (e.g. oxalate and acetate) at pH 2-3. These studies showed that high acidity and the presence of oxalate enhanced Fe dissolution at the surface of CFA particles, similar to those reported in mineral dust Chen and Grassian, 2013;Fu et al., 2012;Ito and Shi, 2016;Shi et al., 2011a). They also demonstrated that there are large differences in dissolution rates in different types of CFA, likely related to the Fe speciation.
Furthermore, high ionic strength, commonly seen in aerosol water, affects the activity of molecular species present in solution; consequently, it can significantly impact the Fe dissolution behaviour. Recent studies have considered the effect of the high ionic strength on the Fe dissolution kinetics of CFA under acidic conditions. For example, the Fe solubility of CFA samples was measured at pH 1-2 with high sodium chloride (NaCl) concentrations (Borgatta et al., 2016) and with high sodium nitrate (NaNO 3 ) concentrations (Kim et al., 2020). In real atmospheric conditions, NaCl or NaNO 3 are unlikely to be the main drivers of high ionic strength in aged CFA. Although NaCl can coagulate with dust particles in the marine boundary layer (Zhang et al., 2003), the ageing of CFA is primarily by the uptake of secondary species, particularly sulfate and ammonia (Li et al., 2003). Ito and Shi (2016) found that, at low pH and a high concentration of (NH 4 ) 2 SO 4 , the Fe solubility of mineral dust is likely to be enhanced by the adsorption of sulfate ions on the particle surface. However, to date the effect of high (NH 4 ) 2 SO 4 concentrations on the Fe dissolution behaviour in combustion sources in the presence or absence of oxalate remains unknown.
In this study, laboratory experiments were conducted to determine the dissolution kinetics of coal combustion emission products (i.e. CFA) during simulated atmospheric acidic processing in the presence of (NH 4 ) 2 SO 4 and oxalate, which are commonly found in atmospheric aerosols. In particular, we investigated the effect of high (NH 4 ) 2 SO 4 concentrations on the proton-promoted and oxalate-promoted Fe dissolution at low pH conditions. Our study also determined the Fe-bearing phases present in the CFA and compared them to those present in mineral dust. The experimental results enabled us to develop a new Fe release scheme for CFA sources, which was then implemented into the global atmospheric chemical transport model IMPACT (Integrated Massively Parallel Atmospheric Chemical Transport). The model results were compared with observations of Fe solubility in aerosol particles over the Bay of Bengal from Bikkina et al. (2020).

Sample collection and subsequent size fractionation
CFA samples were collected from the electrostatic precipitators at three coal-fired power stations at different locations in the United Kingdom (Aberthaw ash), Poland (Krakow ash), and China (Shandong ash). The bulk samples were resuspended to obtain aerosol fractions representative of particles emitted into the atmosphere. A custom-made resuspension system was used to collect the PM 10 fraction (particles with an aerodynamic diameter smaller than 10 µm), which is shown in Fig. S1 in the Supplement. Around 20 g of sample was placed into a glass bottle and injected at regular intervals (2-5 s) into a glass reactor (∼ 70 L) by flushing the bottle with pure nitrogen. The air in the reactor was pumped at a flow rate of 30 L min −1 into a PM 10 sampling head. Particles were collected on 0.6 µm polycarbonate filters and transferred into centrifuge tubes. The system was cleaned manually and flushed for 10 min with pure nitrogen before loading a new sample. A soil sample from Libya (Soil 5;32.29237 • N, 22.30437 • E) was dry sieved to 63 µm and used as an analogue for a Saharan mineral dust precursor to make a comparison between CFA and mineral dust.

Fe dissolution kinetics
The Fe dissolution kinetics of the CFA samples were determined by time-dependent leaching experiments. We followed a similar methodology as in Ito and Shi (2016). PM 10 fractions were exposed to H 2 SO 4 solutions at pH 1, 2, or 3 in the presence of H 2 C 2 O 4 and/or (NH 4 ) 2 SO 4 to simulate acidic processing in aerosol conditions. The concentration of H 2 C 2 O 4 in the experiment solutions was chosen based on the molar ratio of oxalate and sulfate in PM 2.5 (particles with an aerodynamic diameter smaller than 2.5 µm) from observations over the East Asia region (Yu et al., 2005). Around 50 mg of CFA was leached in 50 mL of acidic solution to obtain a particles / liquid ratio of 1 g L −1 . The sample solution was mixed continuously in a rotary mixer, in the dark, at room temperature. A volume of 0.5 mL was sampled at fixed time intervals (2.5, 15, and 60 min and 2, 6, 24, 48, 72, and 168 h after the CFA sample was added to the experiment solution) and filtered through 0.2 µm pore size syringe filters. The dissolved Fe concentration in the filtrate was determined using the ferrozine method (Viollier et al., 2000). Leaching experiments were also conducted on the Libyan dust precursor sample. The relative standard deviation (RSD) at each sampling time varied from 4 % to 15 % (n = 7).
The pH of all the experiment solutions was calculated using the E-AIM III model for aqueous solutions (Wexler and Clegg, 2002). This was partly because the high ionic strength generated by the elevated concentration of (NH 4 ) 2 SO 4 prevents electrochemical sensors from making accurate pH measurements. For the experiment solutions with no (NH 4 ) 2 SO 4 , the pH was measured by a pH meter before adding the ash and at the end of the experiments. The solution pH increased after adding the ash, and the change in pH was used to estimate the buffer capacity of alkaline minerals in the samples, including for example calcium carbonates (CaCO 3 ), lime (CaO), and portlandite (Ca(OH) 2 ). The estimated concentration of the H + buffered was used to input the concentration of H + into the E-AIM model. For each experiment, the pH was calculated before adding the CFA samples and at the end of the experiments. The pH of the original solution before adding the samples was estimated from the molar concentrations (mol L −1 ) of H 2 SO 4 , H 2 C 2 O 4 , and (NH 4 ) 2 SO 4 used to prepare the solution. The model inputs included the total concentrations of H + (without H 2 C 2 O 4 contribution), NH + 4 , SO 2− 4 , and H 2 C 2 O 4 . For the experiment solutions with no (NH 4 ) 2 SO 4 , we calculated the final pH by reducing the total H + concentration input into the model to match the pH measured at the end of the experiments. The buffered H + was then estimated from the difference between the original and final H + concentration input into the model. To determine the final pH of the solutions with high ionic strength, the H + concentration input in the model was calculated as the difference between the H + concentration in the original solution and the buffered H + estimated at low ionic strength.
For the solution with no (NH 4 ) 2 SO 4 , the difference between calculated and measured pH is < 7 %. Table S1 in the Supplement reports the concentrations of H 2 SO 4 , H 2 C 2 O 4 , and (NH 4 ) 2 SO 4 in the experiment solutions, the original and final pH from model estimates (including H + concentrations and activities), and the pH measurements for the solution with low ionic strength.

Sequential extractions
The content of Fe oxide species in the samples was determined by Fe sequential extraction (Baldo et al., 2020;Poulton and Canfield, 2005;Raiswell et al., 2008;Shi et al., 2011b). The Fe oxide species included highly reactive amorphous Fe oxide-hydroxide (FeA), crystalline Fe oxidehydroxide, mainly goethite and hematite (FeD), and Fe associated with magnetite (FeM).
To extract FeA, samples were leached in an ascorbate solution buffered at pH 7.5 (Raiswell et al., 2008;Shi et al., 2011b). The ascorbate solution contained a deoxygenated solution of 50 g L −1 sodium citrate, 50 g L −1 sodium bicarbonate, and 10 g L −1 of ascorbic acid. Around 30 mg of CFA was leached for 24 h in 10 mL of ascorbate extractant and mixed continuously in a rotary mixer. The extraction solution was then filtered through a 0.2 µm membrane filter. In order to extract FeD, the residue was leached for 2 more hours in a dithionite solution buffered at pH 4.8 (50 g L −1 sodium dithionite in 0.35 M acetic acid and 0.2 M sodium citrate; Raiswell et al., 2008;Shi et al., 2011b).
For the extraction of FeM, the CFA samples were first leached for 2 h, using a citrate-buffered dithionite solution which, not sequential to the ascorbate extraction, will remove FeD and FeA. The residue collected after filtration was then leached for 6 h in a solution of 0.2 M ammonium oxalate ((NH 4 ) 2 C 2 O 4 ) and 0.17 M H 2 C 2 O 4 at pH 3.2 (Poulton and Canfield, 2005). The Fe extractions were all carried out in the dark at room temperature. The Fe concentration in the filtered extraction solutions was measured using the ferrozine method (Viollier et al., 2000) or by inductively coupled plasma optical emission spectrometry (ICP-OES) analysis for the solutions containing a high concentration of oxalate.
The total Fe content in the samples was determined by microwave digestion in concentrated nitric acid (HNO 3 ), followed by inductively coupled plasma mass spectrometry (ICP-MS) analysis. A detailed description of the digestion method is provided in the Supplement (Text S1). The total Fe content obtained for the Arizona Test Dust (ATD; ISO 12103-1, A1 Ultrafine Test Dust; Powder Technology Inc.) was comparable with the latest consensus value for the total Fe in ATD, which indicates a good recovery (94.0 % ± 1.5 %). The recovery of Fe assessed using a standard reference material for urban particulate matter (National Institute of Standards and Technology Standard Reference Material -NIST SRM 1648A) was 89.0 % ± 0.4 %. It is possible that some of the Fe in aluminosilicate minerals are not fully digested, but the uncertainty associated with this analytical method is very small, particularly when we compare this with the large uncertainty in simulated Fe solubility in models.
The sequential extraction techniques were tested using the ATD. The weight percentage (wt %) of Fe obtained for each extract using the ATD was 0.057 ± 0.002 for FeA, 0.394 ± 0.045 for FeD, 0.047 ± 0.006 for FeM (n = 7), and 3.501 ± 0.056 for the total Fe (n = 3). A summary of the results for the ATD is reported in Table S2.

X-ray absorption near-edge structure (XANES) analysis
We collected XANES spectra to qualitatively examine the Fe speciation in the CFA samples. The XANES spectra at the Fe K edge were collected at the Diamond Light Source beamline I18. A Si(111) double-crystal monochromator was used in the experiments. The beam size was 400 µm × 400 µm. The XANES spectra were collected from 7000 to 7300 eV at a resolution varying from 0.2 eV for 3 s in proximity to the Fe K edge (7100-7125 eV) to 5 eV for 1 s from 7100 to 7300 eV. Powder samples were suspended in methanol and deposited on Kapton ® tape. The analysis was repeated three times. We measured the XANES spectra of the CFA PM 10 fractions and mineral standards, including hematite, magnetite, and illite. Data were processed using the Athena programme, which part of the software package Demeter (version 0.9.26; Ravel and Newville, 2005).

Model description
This study used the IMPACT model (Ito et al., 2021a, and references therein). The model simulates the emission, chemistry, transport, and deposition of Fe-containing aerosols and the precursor gases of inorganic and organic acids. The coating of acidic species on the surface of Fe-containing aerosols promotes the release of soluble Fe in the aerosol deliquescent layer and enhances the aerosol Fe solubility (Li et al., 2017). On the other hand, the external mixing of oxalaterich aerosols with Fe-rich aerosols can suppress the oxalatepromoted Fe dissolution at a low concentration of oxalate near the source regions (Ito, 2015). However, the internal mixing of alkaline minerals, such as calcium carbonate with Fe-containing dust aerosols, can suppress the Fe dissolution (Ito and Feng, 2010). Since CFA particles are co-emitted with acidic species, the transformation of relatively insoluble Fe in coal combustion aerosols into dissolved Fe is generally much faster than that for mineral dust aerosols during their atmospheric lifetime (Ito, 2015;Ito and Shi, 2016). Additionally, the size of CFA particles is substantially smaller than that of mineral dust. Thus, we adopted an observationally constrained parameter for the dry deposition scheme (Emerson et al., 2020) to improve the simulation of dry deposition velocity of fine particles.
To improve the accuracy of our simulations of Fecontaining aerosols, we revised the online Fe dissolution schemes in the original model (Ito et al., 2021a) in conjunction with a more dynamic range of pH estimates. To apply the Fe dissolution schemes for high ionic strength in aerosols, we used the mean activity coefficient for pH estimate (Pye et al., 2020). Moreover, the dissolution rate was assumed to be de-pendent of pH for highly acidic solutions (pH < 2), unlike in the former dissolution scheme (Ito, 2015), which allowed us to predict the sensitivity of Fe dissolution to pH lower than 2.
To validate the new dissolution scheme, we compared our model results with observations of Fe solubility in PM 2.5 aerosol particles over the Bay of Bengal (Bikkina et al., 2020).

Fe dissolution kinetics
We determined that Krakow ash had the largest buffer capacity, around 0.008 mol of buffered H + per litre, which was related to the content of alkaline minerals in the sample. The buffer capacity of Aberthaw and Shandong ash was ∼ 10 times smaller than that of Krakow ash, at around 0.0007 mol of buffered H + per litre. Leaching Krakow ash in 0.005 M H 2 SO 4 , the initial concentration of H + was similar to the concentration of the buffered H + . As a result, the solution pH raised from approximatively 2.1 to 2.7 corresponding to a pH change of around 20 % (Table S1). For all the other experimental conditions, the pH change was below 12 % (Table S1). At the pH conditions used in this study (pH 1-3), acid buffering was fast and likely occurred within the first 1-2 h. We assumed that the calculated final pH was representative of the solution pH over the duration of the experiments. The leaching experiments were conducted up to 168 h to better capture the dissolution curve in the kinetic model but also to consider the tropospheric lifetime of aerosol particles.
Dissolved Fe at different time intervals is reported as the percent of Fe, which is the fraction of Fe dissolved to the total Fe content (FeT) in the CFA samples. For all samples, a fast dissolution rate was observed at the beginning of the experiment. In the case of Krakow ash, the dissolution plateau was reached after 2 h of leaching in 0.005 M H 2 SO 4 , as sufficient Fe may be dissolved from the highly reactive Fe species to suppress the dissolution of less reactive Fe. For that sample/initial condition, the pH increased to 2.7, and no more Fe was dissolved, leading to a total Fe solubility of ∼ 9 % over the duration of the experiment (7 d; Fig. 1a). Dissolving Krakow ash in 0.01 M H 2 SO 4 (Fig. 1a), the experiment solution had a final calculated pH of 2.1. The total Fe solubility was 34 % at pH 2.1, which is almost 4 times higher than that at pH 2.7 (in 0.005 M H 2 SO 4 ). The dissolution of Aberthaw and Shandong ash was slower compared to Krakow ash (Figs. 1b and 2c, respectively). The leaching of Aberthaw and Shandong ash in 0.005 M H 2 SO 4 resulted in solutions with a pH of around 2.2. At this pH, the total Fe solubility was 18 % for Aberthaw ash and 21 % for Shandong ash, which is 9-10 times higher than the total Fe solubility at pH 2.9 (in 0.001 M H 2 SO 4 ), which is around 2 % for both samples. The final pH of the experiment solutions is also reported, which was calculated using the E-AIM III model for aqueous solutions (Wexler and Clegg, 2002), accounting for the buffer capacity of the CFA samples (Experiments 1-2 in Table S1). The experiments conducted at around pH 2 are in red, while the experiments at around pH 3 are in black. The data uncertainty was estimated using the error propagation formula.
The experimental treatment of dissolved Fe from Krakow ash in 0.05 H 2 SO 4 solution with 1 M (NH 4 ) 2 SO 4 ( Fig. 1a) resulted in a final predicted pH of 2.1. At that pH, the total Fe solubility of Krakow ash increased from 34 %, with no (NH 4 ) 2 SO 4 , to 48 %, with a high (NH 4 ) 2 SO 4 concentra- shown. The final pH of the experiment solutions is also reported, which was calculated using the E-AIM III model for aqueous solutions (Wexler and Clegg, 2002) and accounting for the buffer capacity of the CFA samples (Experiments 1 and 3-4 at around pH 2). The data uncertainty was estimated using the error propagation formula.
tion. The total Fe solubility of Krakow ash was around 28 % at pH 3.0 with 1 M (NH 4 ) 2 SO 4 (Fig. 1a), which is 3 times higher than that at pH 2.7 with no (NH 4 ) 2 SO 4 . At around pH 2, the total Fe solubility of Aberthaw (Fig. 1b) and Shandong ash (Fig. 1c) increased by around 20 % and 30 % in the presence of (NH 4 ) 2 SO 4 . By contrast, the total Fe solubility shown. The final pH of the experiment solutions is also reported, which was calculated using the E-AIM III model for aqueous solutions (Wexler and Clegg, 2002) and accounting for the buffer capacity of the CFA samples (Experiment 7 at pH 1.0, Experiment 3 at pH 2.0, and Experiment 3 at pH 2.9 in Table S1). The data uncertainty was estimated using the error propagation formula. at pH 3.1 with 1 M (NH 4 ) 2 SO 4 was 7.5 % for Aberthaw ash (Fig. 1b) and 14 % for Shandong ash (Fig. 1c), respectively, which was between 4 and 7 times higher than in the experiments carried out at pH 2.9 without (NH 4 ) 2 SO 4 .
The Fe dissolution of the CFA samples in H 2 SO 4 solutions with 0.01 M H 2 C 2 O 4 (at around pH 2) is shown in Fig. 2. The total Fe solubility of Krakow ash at pH 1.9 with 0.01 M H 2 C 2 O 4 was 61 % (Fig. 2a), which was almost 2 times higher than that at pH 2.1 but without H 2 C 2 O 4 (Fig. 2a). For Aberthaw ash, the oxalate contribution to the dissolution process led to a total Fe solubility of 30 % at pH 2.0 ( Fig. 2b), which was 70 % higher than in the experiment carried out in 0.005 M H 2 SO 4 (∼ pH 2.2; Fig. 2b). Shandong ash dissolution behaviour was not affected by the presence of oxalate (Fig. 2c).
We also investigated the effect of a high (NH 4 ) 2 SO 4 concentration on oxalate-promoted dissolution. In Fig. 2a, the total Fe solubility of Krakow ash decreased from 61 % at pH 1.9 in the presence of oxalate to 54 % at pH 2.0 with oxalate and (NH 4 ) 2 SO 4 . For Aberthaw ash, the total Fe solubility at pH 2.0 decreased from 30 % in the presence of oxalate to 19 % after the addition of (NH 4 ) 2 SO 4 (Fig. 2b). Figure 3 shows the Fe dissolution behaviour of Krakow ash at different pH conditions in the presence of 1 M (NH 4 ) 2 SO 4 and H 2 C 2 O 4 (0.01-0.03 M, depending on the solution pH). The total concentration of oxalate ions was calculated using the E-AIM model and was similar at different pH conditions, i.e. 0.015 at pH 1.0 (Experiment 7 in  (Wexler and Clegg, 2002) and accounting for the buffer capacity of the CFA samples (Experiments 5-8 in Table S1). The data uncertainty was estimated using the error propagation formula. Table S3), 0.009 at pH 2.0, and 0.01 at pH 2.9 (Experiment 3 in Table S3). The highest total Fe solubility was observed at pH 1.0 (∼ 67 %). At pH 2.0, the total Fe solubility decreased to 54 %, and no substantial variations were observed between pH 2.0 and 2.9 (54 %-51 %). At pH 1.0, the concentration of H + was considerably higher compared to pH 2.0-2.9, leading to a faster dissolution rate. The total concentration of oxalate ions was 1.5-1.6 times higher in the solution at pH 1.0 than at pH 2.0-2.9, which may also contribute to the faster dissolution rate. C 2 O −2 4 concentration increased with rising pH. Although the concentration of H + was lower at pH 2.9 than at pH 2.0, the E-AIM model estimated that C 2 O −2 4 contributed around 35 % of the total oxalate concentration at pH 2.9, which was 4.5 times higher than at pH 2.0 (Experiment 3 in Table S3). The similar dissolution behaviour at pH 2.0 and 2.9 conditions may reflect the combination of these two opposite factors, with a higher concentration of C 2 O −2 4 but lower concentration of H + at pH 2.9 compared to 2.0.
We determined the Fe dissolution behaviour of Krakow ash at pH 1.0 in the presence of oxalate and increasing concentrations of (NH 4 ) 2 SO 4 . The ash was leached in H 2 SO 4 solutions, with 0.03 M H 2 C 2 O 4 at pH 1.0, while the concentration of (NH 4 ) 2 SO 4 varied from 0 to 1.5 M. In Fig. 4, the total Fe solubility of Krakow ash in the presence of oxalate was 75 % at pH 1.0 and decreased to 68 % after the addition of 0.5 M (NH 4 ) 2 SO 4 . Higher (NH 4 ) 2 SO 4 concentrations did not affect the Fe dissolution behaviour in the presence of oxalate at pH 1.0.

Fe speciation
The Fe-bearing phases in the CFA samples determined through sequential extractions are shown in Fig. 5c. The Fe speciation in the Libyan dust precursor is added for comparison. Krakow ash had a total Fe (FeT) content of 5.2 %, while FeT in Aberthaw and Shandong ash was 3.1 % and 1.6 %, respectively. Amorphous Fe (FeA/FeT) was 6.5 % in Krakow ash, 2 % in Aberthaw ash, and 4.6 % in Shandong ash. The CFA samples showed very different dithionite Fe (FeD/FeT) content, 21.5 % in Krakow ash, 8 % in Aberthaw ash, and 14.8 % in Shandong ash. The content of magnetite (FeM/FeT) was considerably higher in Krakow ash (22.4 %) compared to Aberthaw (2.9 %) and Shandong (4.5 %) ash. About 50 %-87 % of Fe was contained in other phases, most likely in aluminosilicates. Overall, CFA had more magnetite and highly reactive amorphous Fe and less dithionite Fe than the Libyan dust precursor sample.
In Fig. 5a-b, the Fe K-edge XANES spectra of Krakow and Aberthaw ash showed a single peak in the pre-edge region at around 7114.3 and 7114.6 eV, respectively. In the edge region, Aberthaw ash showed a broad peak at around 7132.2 eV, while the peak of Krakow ash was slightly shifted to 7132.9 eV and narrower. The pre-edge peak at around 7115.4 eV suggests that the Fe was mainly found as Fe(III). The spectral features of Aberthaw and Krakow ash are different from those of the hematite, magnetite, and illite standards, suggesting that the glass fraction was dominant and controlled their spectral characteristics, which is consistent with the results of the Fe sequential extractions. The XANES Fe K-edge spectra of the CFA samples have some common features with those of Icelandic dust but tend to differ from the mineral dust sourced in the Saharan dust source region. In the pre-edge region of the spectrum, Icelandic dust (sample D3 in Fig. 5a-b) showed a main peak at around 7114.4 eV and a second, less intense, peak at around 7112.7 eV, while a broad peak was observed at around 7131.9 eV in the edge region (Baldo et al., 2020). A mineral dust sample from the western Sahara (WS dust in Fig. 5a-b) showed a distinct double peak in the pre-edge region at around 7113.9 and 7115.2 eV and a main peak in the edge region at around 7133.3 eV (Baldo et al., 2020). The similarities between Icelandic ash and CFA could be because aluminium silicate glass is dominant in these samples (e.g. Baldo et al., 2020;Brown et al., 2011), while Fe-bearing phases in mineral dust from the Saharan region are primarily iron oxide minerals such as hematite and goethite, clay minerals, and feldspars (e.g. Shi et al., 2011b).

Fe dissolution scheme
Based on the laboratory experiments carried out on the CFA samples, we implemented a three-step dissolution scheme for proton-promoted and oxalate-promoted Fe dissolution (Table 1). The Fe dissolution kinetics were described as follows (Ito, 2015): where RFe i is the dissolution rate of individual mineral i, k i is the rate constant (mol Fe g −1 s −1 ), a(H + ) is the H + activity in solution, and m i represents the empirical reaction order for protons. The function f i (0 ≤ f i ≤ 1) accounts for the suppression of mineral dissolution by competition for oxalate between surface Fe and dissolved Fe, as follows (Ito, 2015): where [Fe] is the molar concentration (mol L −1 ) of Fe 3+ dissolved in a solution, and [lig] is the molar concentration of a ligand (e.g. oxalate). f i was set to 1 for the proton-promoted dissolution.
The scheme assumes three rate constants, namely fast, intermediate, and slow, for the proton-promoted and the proton + oxalate-promoted dissolution (Table 1). These were obtained by fitting the parameters to our measurements for Krakow ash in H 2 SO 4 and (NH 4 ) 2 SO 4 at pH 2-3, with and without oxalate (Experiments 2 and 3 in Table S1), which are shown in Fig. 6. The fast rate constant represents a highly reactive Fe species such as amorphous Fe oxyhydroxides, Fe carbonates, and Fe sulfates. The intermediate rate constant can be applied to nanoparticulate Fe oxides, while more stable phases including, for example, Fe aluminosilicate and crystalline Fe oxides have generally slower rates (Ito and Shi, 2016;Shi et al., 2011aShi et al., , b, 2015. Similarly, we predicted the dissolution kinetics of Aberthaw ash and Shandong ash (Fig. 7). The dissolution kinetics of Krakow ash were calculated based also on the experimental results at pH 1.0, which are shown in Fig. S2 in comparison with kinetics predicted at pH 2.0 and pH 2.9 conditions. The contribution of the oxalate-promoted dissolution to dissolved Fe was derived as the difference between the estimated dissolution rates for the proton + oxalate-promoted dissolution and the proton-promoted dissolution as follows: The Fe dissolution rates were predicted at a wider range of pH, using Eqs.
Since RFe i(oxalate) is less than 0 at low pH (< 2), this equation applies to highly acidic conditions. As a result, the predicted amount of dissolved Fe was smaller when using the dissolution rate for the proton + oxalate-promoted dissolution, RFe i(proton + oxalate) , rather than the rate for the protonpromoted dissolution, RFe i(proton) , at pH < 2. Accordingly, The FeT (as wt %) is given below each sample column. The data uncertainty was estimated using the error propagation formula, with 4 % for FeA/FeT, 11 % for FeD/FeT, 12 % for FeM/FeT, and 2 % for FeT. the dissolution rate, RFe i , was less dependent on the pH compared to RFe i(proton) at highly acidic conditions, possibly due to the competition for the formation of surface complexes. At pH > 2 when oxalate does promote Fe dissolution, the following equation applies: when RFe i(oxalate) > 0. (5)

Aerosol Fe solubility over the Bay of Bengal
The new dissolution scheme was applied in the IMPACT atmospheric chemistry transport model to predict the Fe solubility in atmospheric particles collected over the Bay of Ben-gal, which is an area for which there are detailed field measurements available (Bikkina et al., 2020;Kumar et al., 2010;Srinivas and Sarin, 2013;Srinivas et al., 2012) and multimodelling analyses have been done (Ito et al., 2019). It, thus, represents a test for our experimental results in actual field conditions. A total of three sensitivity simulations were performed to explore the effects of the uncertainties associated with the dissolution schemes and mineralogical component of Fe. In addition, the former setting (Ito et al., 2021a) was used in the IMPACT model for comparison. For all simulations, the total Fe emissions from anthropogenic combustion sources and biomass burning were estimated using the Fe emission inventory of Ito et al. (2018), a k(pH, T ) is the pH and temperature-dependent far-from-equilibrium rate constant (mol Fe g −1 s −1 ). The Fe dissolution scheme assumes three rate constants of fast, intermediate, and slow for the proton-and oxalate-promoted dissolution. The parameters were fitted to our measurements for Krakow ash. b E(pH) = −1.56 × 10 3 × pH + 1.08 × 10 4 . The parameters were fitted to the measurements for soils (Bibi et al., 2014). c m is the reaction order with respect to aqueous-phase protons, which was determined by a linear regression from our experimental data in the pH range between 2 and 3 for proton-and oxalate-promoted dissolution schemes. is also reported, which was calculated using the E-AIM III model for aqueous solutions (Wexler and Clegg, 2002) and accounting for the buffer capacity of the CFA samples (Experiments 2-3 in Table S1). The experiments conducted at around pH 2 are in red, while the experiments at around pH 3 are in black. The data uncertainty was estimated using the error propagation formula.
including also emissions from the iron and steel industry, whereas Fe emissions from mineral dust sources were dynamically simulated (Ito et al., 2021a). In Test 0, we ran the model without the upgrades of the dissolution scheme discussed in Sect. 2.4 and additionally applied the photo-induced dissolution scheme for both combustion and dust aerosols (Ito, 2015;Ito and Shi, 2016), which was turned off in Test 1, Test 2, and Test 3 due to the lack of laboratory measurements under high ionic strength. To estimate the aerosol pH, we applied a H + activity coefficient of 1 for Test 0, while Figure 7. Comparison between the Fe dissolution kinetics of Krakow, Aberthaw, and Shandong ashes predicted using Eq. (1) and measured in (a-c) H 2 SO 4 solutions at around pH 2 with 1 M (NH 4 ) 2 SO 4 (Experiment 2 at around pH 2 in Table S1), (d-f) H 2 SO 4 solutions at around pH 3 with 1 M (NH 4 ) 2 SO 4 (Experiment 2 at around pH 3 in Table S1), and (g-i) H 2 SO 4 solutions at pH 2.0 with 0.01 M H 2 C 2 O 4 and 1 M (NH 4 ) 2 SO 4 (Experiment 3 at pH 2.0 in Table S1). The molar concentrations of H 2 SO 4 , H 2 C 2 O 4 , and (NH 4 ) 2 SO 4 in the experiment solutions are shown. The final pH of the experiment solutions is also reported, which was calculated using the E-AIM III model for aqueous solutions (Wexler and Clegg, 2002) and accounting for the buffer capacity of the CFA samples. the mean activity coefficient from Pye et al. (2020) was used for the other tests. The dissolution rate was assumed to be pH independent for highly acidic solutions (pH < 2; Ito, 2015) in Test 0, based on the laboratory measurements in , while no pH threshold was considered in Test 1, Test 2, and Test 3, as the total dissolution (proton + oxalate) was suppressed at pH < 2 from the predicted dissolution rate.
In Test 1, we used the new dissolution scheme accounting for the proton-and oxalate-promoted dissolution of Krakow ash for all combustion aerosols in the model (Table 1). The dissolution kinetics were calculated using the base mineralogy for anthropogenic Fe emissions reported in Table S11 of Rathod et al. (2020). The Fe composition of wood was used for open biomass burning (Matsuo et al., 1992). In this simulation, three Fe pools were considered. Sulfate Fe in Rathod et al. (2020) was assumed to be fast pool, magnetite Fe as intermediate pool, and hematite, goethite, and clay as slow pool. In Test 2, we calculated the dissolution kinetics only considering the proton-promoted dissolution. In Test 3, the Fe pools were as determined here for Krakow ash, with ascorbate Fe (FeA) as fast pool, magnetite Fe (FeM) as intermediate pool, and hematite plus goethite Fe (FeD) and other Fe as slow pool (Fig. 5). FeA contains highly reactive Fe species with fast dissolution rates (Raiswell et al., 2008;Shi et al., 2011b). FeM appeared to work well for the different fly ash samples in the dissolution scheme as intermediate Fe pool. FeD is associated with crystalline Fe oxides which are mostly highly insoluble (Raiswell et al., 2008;Shi et al., 2011b), thus it was considered to be a slow pool in the dissolution scheme. We assumed other Fe to be mostly Fe-bearing aluminosilicates and considered this to be a slow Fe pool.
Observations of total Fe concentration and Fe solubility in PM 2.5 along the cruise tracks over the Bay of Bengal for the period extending from 27 December 2008 to 26 January 2009 (Bikkina et al., 2020) were compared with temporally and regionally averaged data from model estimates. The daily averages of model results were calculated from hourly mass concentrations in the air over the surface ocean along the cruise tracks. The concentration of total Fe observed over the Bay of Bengal varies from 145 ± 144 ng m −3 over the northern Bay of Bengal (27 December 2008-10 January 2009) to 55 ± 23 ng m −3 over the southern Bay of Bengal (11-26 January 2009;Bikkina et al., 2020). In Fig. 8, the modelled concentrations of total Fe exhibit a similar variability to that of measurements with relatively higher values over the northern Bay of Bengal (59 ± 29 ng m −3 in different sensitivity simulations) compared to the southern Bay of Bengal (20 ± 12 ng m −3 in different sensitivity simulations). However, the modelled concentrations of total Fe were underestimated by a factor of 2.9 ± 1.5. The model reproduced the source apportion of Fe ( Fig. 8; Table S4), which is qualitatively derived from previous observational studies, indicating that the concentrations of total Fe in aerosols over the northern Bay of Bengal are influenced by the emissions of dust and combustion sources from the Indo-Gangetic Plain (Kumar et al., 2010), whereas combustion sources (e.g. biomass burning and fossil-fuel) from southeastern Asia are dominant over the southern Bay of Bengal (Kumar et al., 2010;Srinivas and Sarin, 2013). On the other hand, the model could not reproduce the peak in the total Fe concentration (1.8 % of the Fe content in PM 2.5 sample) reported around 29 December 2008. The total Fe observed in PM 10 (430 ng m −3 ) on 29 December 2008 is lower than that measured on the day before (667 ng m −3 ) and the day after (773 ng m −3 ), whereas that in PM 2.5 peaked on 29 December 2008 (Srinivas et al., 2012). Thus, the extreme value recorded only for PM 2.5 on this date may be an outlier. But we do not have sufficient data to confirm this. One of the possibilities is that the sample collected aerosol particles from a mixture of different aerosols sources (e.g. dust and anthropogenic aerosols). This reflects one of the challenges of modelling such a dynamic parameter.
The comparison of Fe solubility using the same total Fe emissions directly represents the effect of the new dissolution scheme on PM 2.5 . The aerosol Fe solubility measured over the southern Bay of Bengal is higher than that over the northern Bay of Bengal, respectively, with 32 % ± 11 % and 15 % ± 7 % (Bikkina et al., 2020), and model estimates showed a similar trend (Fig. 9). In Fig. 9 and Table S5, the calculated Fe solubilities over the northern Bay of Bengal in Test 1 (11 % ± 4 %), Test 2 (17 % ± 5 %), and Test 3 (17 % ± 6 %) were in good agreement with observations. The aerosol Fe solubility over the southern Bay of Bengal was better captured in Test 1 (30 % ± 5 %) and Test 3 (37 % ± 7 %), whereas Test 0 showed higher variability (37 % ± 22 %). The proton-promoted dissolution scheme in Test 2 significantly overestimated the Fe solubility over the Bay of Bengal ( Fig. 9 and Table S5). The aerosol Fe solubility was largely overestimated in all scenarios after 22 January 2009, as open biomass burning sources become dominant ( Fig. 8 and Table S4).
The comparison between observations and model predictions of aerosol Fe solubility over the Bay of Bengal is shown in Fig. S3. The agreement between the measurements and model predictions was the best in Test 1 and Test 3. These exhibited a good correlation with observations (R = 0.49 in Test 1 and R = 0.54 in Test 3), and the lowest root mean squared difference between the simulated and observed Fe solubilities (root mean square error (RMSE) equal to 11 in Test 1 and RMSE = 12 in Test 3). In Test 0, the model estimates showed a greater difference from observations (RMSE = 21) and poor correlation (R = 0.26).  Bikkina et al. (2020). Model estimates of Test 0, Test 1, Test 2, and Test 3 were calculated along the cruise tracks using the IMPACT model. In Test 0, we run the model without upgrades (Ito et al., 2021a) and apply the proton-promoted, oxalate-promoted, and photoinduced dissolution schemes for combustion aerosols in Table S6 (Ito, 2015). The proton + oxalate dissolution scheme (Table 1) was applied in Tests 1 and 3, while proton-promoted dissolution is used for Test 2. We adopted the base mineralogy for anthropogenic Fe emissions (Rathod et al., 2020) in Tests 1 and 2. In Test 3, the Fe speciation of Krakow ash was used for all combustion sources. The small white square within the box shows the mean. The solid line within the box indicates the median. The lower and upper hinges correspond to the 25th and 75th percentiles. The whiskers above and below the box indicate the 1.5× interquartile range, and the data outside this range are plotted individually.

Dissolution behaviour of Fe in CFA
In this study, the Fe dissolution kinetics of CFA samples from the UK, Poland, and China were investigated under simulated atmospheric acidic conditions. A key parameter in both the atmosphere and the simulation experiments is the pH of the water interacting with the CFA particles. The lower the pH of the experimental solution, the faster the dissolution and eventually the higher the amount of Fe dissolved. Our results showed a strong pH dependence in low ionic strength conditions, with higher dissolution rates at lower pH. For example, reducing the solution pH from 2.7 to 2.1, the Fe solubility of Krakow ash in H 2 SO 4 only increased by a factor of 4 (Fig. 1a) over the duration of the experiments, while the Fe solubility of Aberthaw and Shandong ash increased by 9-10 times from pH 2.9 to 2.2 (Fig. 1b-c). This enhancement is higher than that observed in studies conducted on mineral dust samples, which showed that one pH unit can lead to 3-4 times difference in dissolution rates (Ito and Shi, 2016;Shi et al., 2011a). Furthermore, Chen et al. (2012) reported that the Fe solubility of the certified CFA 2689 only increased by 10 % from pH 2 to 1, after 50 h of dissolution in acidic media. The Fe solubility of CFA (PM 10 fractions) after 6 h at pH 2 was 6 %-10 % for Aberthaw and Shandong ash, respectively, and 28 % for Krakow ash (Fig. 1). The Fe in our CFA samples initially dissolved faster than the samples used by Fu et al. (2012), who reported 2.9 %-4.2 % Fe solubility in bulk CFA from three coal-fired power plants in China after 12 h leaching at pH 2. These results suggest that there are considerable variabilities in the pH-dependent dissolution of Fe in CFA. This could be due to differences in the Fe speciation between CFA samples and/or the different leaching media used.
Our results showed that high ionic strength has a major impact on the dissolution rates of CFA at low pH (i.e. pH 2-3). The Fe solubility of CFA increased by approximatively 20 %-40 % in the presence of 1 M (NH 4 ) 2 SO 4 at around pH 2 over the duration of the experiments and by a factor from 3 to 7 at around pH 3 conditions (Fig. 1). At high ionic strength, the activity of ions in solution is reduced; thus, in order to maintain similar pH conditions, the H + concentration has to be increased (Table S1). Although Fe dissolution was primarily controlled by the concentration of H + , the high concentration of sulfate ions could also be an important factor contributing to Fe dissolution, particularly when the concentration of H + in the system was low (e.g. pH 3). Previous research found that the high ability of anions to form soluble complexes with metals can enhance Fe dissolution (Cornell et al., 1976;Cornell and Schwertmann, 2003;Furrer and Stumm, 1986;Hamer et al., 2003;Rubasinghege et al., 2010;Sidhu et al., 1981;Surana and Warren, 1969). Sulfate ions adsorbed on the particles surface form complexes with Fe (e.g. Rubasinghege et al., 2010). This may increase the surface negative charge favouring the absorption of H + and thereby increase Fe dissolution at the particle surface. In addition, the formation of surface complexes may weaken the bonds between Fe and the neighbouring ions (Cornell et al., 1976;Furrer and Stumm, 1986;Sidhu et al., 1981). Cwiertny et al. (2008) reported that, at pH 1-2, the high ionic strength generated by NaCl up to 1 M did not influence Fe dissolution of mineral dust particles. However, Ito and Shi (2016) showed that the high ionic strength resulting from the addition of 1 M (NH 4 ) 2 SO 4 in leaching solu-6058 C. Baldo et al.: Fe from coal combustion particles dissolves much faster than mineral dust tions at pH 2-3 enhanced the Fe dissolution of dust particles, which was also observed here for the CFA samples. Borgatta et al. (2016) compared the Fe solubility of CFA from the USA Midwest, northeastern India, and Europe in an acidic solution (pH 1-2) containing 1 M NaCl. The Fe solubility measured after 24 h varied from 15 % to 70 % in different CFA (bulk samples) at pH 2 with 1 M NaCl, which was considerably higher than that observed at pH 2 with 1 M NaNO 3 (< 20 %; Kim et al., 2020). Both studies did not investigate the impact of ionic strength on the dissolution behaviour, i.e. by comparing the dissolution at low and high ionic strength. Note that both studies did not specify how the pH conditions were maintained at pH 2. Here, we considered the most important sources of high ionic strength in aerosol water and simulated Fe dissolution in the presence of (NH 4 ) 2 SO 4 and H 2 C 2 O 4 under acidic conditions. We emphasize that the pH under high ionic strength here is estimated from a thermodynamic model, similar to those implemented in the IMPACT model.
The presence of oxalate enhanced the Fe dissolution in Krakow and Aberthaw ash at around pH 2 but not in Shandong ash (Fig. 2). The effect of oxalate on the Fe dissolution kinetics has also been studied by Chen and Grassian (2013) at pH 2 (11.6 mM H 2 C 2 O 4 ). After 45 h of leaching, the Fe solubility of the certified CFA 2689 increased from 16 % in H 2 SO 4 at pH 2 to 44 % in H 2 C 2 O 4 at the same pH (Chen and Grassian, 2013). Therefore, the enhancement in the Fe solubility of CFA in the presence of the oxalate observed in this study (from no impact in Shandong ash to doubled dissolution in Krakow ash) is lower than the 2.8 times increase in Fe solubility reported for the certified CFA 2689 (Chen and Grassian, 2013). Since no data are available in Chen and Grassian (2013), we are unable to make a comparison with the other two certified CFA samples. The Fe solubility of Krakow ash after 48 h of leaching at pH 1.9 with 0.01 M H 2 C 2 O 4 (Fig. 2a) was 53 %, which is within the range of Fe solubilities observed in Chen and Grassian (2013) for the certified CFA samples at similar pH and H 2 C 2 O 4 concentrations (from 44 % to 78 %), whereas the Fe solubility of Aberthaw and Shandong ash (Fig. 2b-c; 18 %-17 % after 48 h of leaching at pH 2.0 with 0.01 M H 2 C 2 O 4 ) was considerably lower than that of certified CFA (Chen and Grassian, 2013). These results suggest a large variability in the effects of oxalate on the Fe dissolution rates in different types of CFA.
Our results also indicated that high (NH 4 ) 2 SO 4 concentrations suppress oxalate-promoted Fe dissolution of CFA (Fig. 2), which was not considered in previous research. At pH 1.9 in the presence of oxalate, the Fe solubility of Krakow ash decreased by around 10 % after the addition of (NH 4 ) 2 SO 4 , while the Fe solubility of Aberthaw ash decreased by 35 % (Fig. 2). We used the E-AIM model to estimate the concentration of oxalate ions and their activity (Table S3). The pH influences the speciation of H 2 C 2 O 4 in the solution (e.g. Lee et al., 2007). H 2 C 2 O 4 is the main species below pH 2, whereas HC 2 O − 4 is dominant between pH 2-4. Above pH 4, C 2 O −2 4 is the principal species. In our experiments, H 2 C 2 O 4 is mainly present as HC 2 O − 4 at around pH 2 (Experiments 3-4 in Table S3). In the presence of (NH 4 ) 2 SO 4 , the activity coefficient of HC 2 O − 4 was reduced by approximatively 35 %-38 % (Experiments 3 in Table S3). Increasing the ionic strength lowers the activity of the oxalate ions, but at the same time favours the dissociation of the acid. At around pH 2 conditions, the E-AIM model estimated that the activity of C 2 O −2 4 was reduced by around 1 order of magnitude in the presence of (NH 4 ) 2 SO 4 , while its concentration increased 12-15 times (Experiments 3 in Table S3). The adsorption of anions can reduce the oxalate adsorption on the particle surface due to electrostatic repulsion, which results in slower release of Fe (Eick et al., 1999). The precipitation of ammonium hydrogen oxalate (NH 4 HC 2 O 4 ) can also occur in the system, but this is very soluble and easily re-dissolves forming soluble oxalate species (Lee et al., 2007). We speculate that the high concentration of sulfate ions is likely to be responsible for inhibiting the oxalate-promoted dissolution by reducing the oxalate adsorption on the particle surface. At pH 1 in the presence of oxalate, increasing the concentration of (NH 4 ) 2 SO 4 from 0.5 to 1.5 M did not affect the Fe dissolution behaviour of the CFA samples (Fig. 4). As previously discussed, the adsorption of sulfate ions on the particle surface may inhibit oxalate-promoted dissolution. However, once the saturation coverage is reached, increasing the concentration of anions has no further effect on the dissolution rate (Cornell et al., 1976).
Fe speciation is an important factor affecting the Fe dissolution behaviour. CFA particles have very different chemical and physical properties depending, for example, on the nature of the coal burnt, the combustion conditions, the cooling process, and the particle control devices implemented at the power stations (e.g. Blissett and Rowson, 2012;Yao et al., 2015). This is likely the reason why the Fe speciation observed in the CFA samples analysed in this study from different locations varied considerably (Fig. 5). In the CFA samples, the Fe dissolution curves for different pH and ionic strengths generally showed the greatest rate of Fe release within the first 2 h, followed by a slower dissolution, reaching almost a plateau at the end of the experimental run. This indicates the presence of multiple Fe-bearing phases in CFA particles with a wide range of reactivity. Initially, highly reactive phases were the main contribution to dissolved Fe. As the dissolution continued, more refractory phases became the dominant source of dissolved Fe (Shi et al., 2011a). Scanning electron microscopy (SEM) analysis conducted on CFA samples showed that CFA particles are mostly spherical (e.g. Chen et al., 2012;Dudas and Warren, 1987;Valeev et al., 2018;Warren and Dudas, 1989), with Fe oxide aggregates on the surface Valeev et al., 2018). The analysis of the CFA samples processed in aqueous solutions at low pH suggests that, initially, Fe dissolved from the reactive external glass coating (Dudas and Warren, 1987;Warren and Dudas, 1989) and from the Fe oxide aggregates on the particle surface Valeev et al., 2018). Subsequently, Fe is likely realized from the structure of the aluminium silicate glass Dudas and Warren, 1987;Valeev et al., 2018;Warren and Dudas, 1989) and crystalline Fe oxide phases (Warren and Dudas, 1989). Overall, Krakow ash showed the fastest dissolution rates, but the dissolution of highly reactive Fe species as FeA is insufficient to account for the high Fe solubility observed at low pH. Our results showed that once the FeA dissolved, additional Fe was dissolved from more refractory Fe-bearing phases. The modelled dissolution kinetics obtained using FeM as intermediate pool were in good agreement with the measurements (Figs. 7 and S2). FeM is likely to be primary magnetite but may contain a fraction of the more reactive aluminosilicate glass. Our model results suggest that magnetite in CFA particles may be more soluble than shown in Marcotte et al. (2020). It is possible that, in real CFA samples, the physicochemical properties of minerals, including, for example, crystal size, degree of crystallinity, cationic, and anionic substitution in the lattice which influence the Fe dissolution behaviour (e.g. Schwertmann, 1991), are likely to be different from those of the reference minerals analysed in Marcotte et al. (2020). In order to investigate the links between Fe solubility and Fe speciation/mineralogy, more work is needed to determine the Fe mineralogy in CFA samples at emission and after atmospheric processing in combination with solubility experiments.
Finally, the modelled dissolution kinetics obtained using the new dissolution scheme for CFA (Table 1) showed better agreement with laboratory measurements than when using the original scheme (Ito, 2015;Fig. 10). In Fig. 10a-b, we compared the Fe dissolution kinetics of Krakow ash at around pH 2 and 3 with 1 M (NH 4 ) 2 SO 4 , which was calculated using the proton-promoted dissolution scheme in Table 1 with the dissolution kinetics calculated at similar pH but using the proton-promoted dissolution scheme for combustion aerosols in Ito (2015; Table S6). The dissolution scheme in Ito (2015) was based on laboratory measurements conducted at low ionic strength  and assumed a single Fe-bearing phase in combustion aerosol particles, while the new dissolution scheme considered the high ionic strength of aerosol water and assumed three rate constants for the fast, intermediate, and slow kinetics of the different Fe-bearing phases present in CFA particles. The Fe dissolution kinetics obtained using the new dissolution scheme showed a better agreement with measurements and was enhanced compared to the model estimates obtained using the original dissolution scheme (Ito, 2015) for low ionic strength conditions (Fig. 10a-b). Figure 10c-d show the Fe dissolution kinetics of Krakow ash at pH 2.0 and 2.9, with 0.01 M H 2 C 2 O 4 and 1 M (NH 4 ) 2 SO 4 , calculated using the proton-and oxalate-promoted dissolution scheme in Table 1 and the dissolution kinetics calculated at a similar pH and H 2 C 2 O 4 concentration but using the scheme in Ito (2015; i.e. single-phase dissolution; see Table S6). The Fe dissolution kinetics predicted using the new dissolution scheme had a much better agreement with measurements. Figure 10e shows the suppression of the oxalate-promoted dissolution at pH 2.0 and high (NH 4 ) 2 SO 4 concentrations. At pH 2, the proton-promoted dissolution was comparable to the proton + oxalate-promoted dissolution (Fig. 10e), with RFe (oxalate) close to zero (see Eq. 3). At pH 2.9, the proton + oxalate-promoted dissolution was higher than the proton + oxalate-promoted dissolution (Fig. 10f), with RFe (oxalate) > 0 (Eq. 5).
Moreover, the new three-step dissolution scheme better captured the initial fast dissolution of CFA (Fig. 10), which was also observed in previous research (Borgatta et al., 2016;Chen et al., 2012;Chen and Grassian, 2013;Fu et al., 2012;Kim et al., 2020; except for the certified CFA 2689 in , which showed increasing dissolution rates over the duration of the experiment). Furthermore, the new scheme enabled us to account for the different Fe speciation determined in the CFA samples, which could be a key factor contributing to the different Fe dissolution behaviour observed in the present study and in the literature (Borgatta et al., 2016;Chen et al., 2012;Chen and Grassian, 2013;Fu et al., 2012;Kim et al., 2020). In Fig. 7, the dissolution kinetics of Aberthaw and Shandong ash calculated using the dissolution rates in Table 1 and the Fe-bearing phases determined in the samples showed a good agreement with the measurements.

Comparison with mineral dust
High ionic strength also impacted the dissolution rates of the Libyan dust precursor sample at low pH (Fig. S4). At around pH 2 conditions, the proton-promoted Fe dissolution of Libyan dust was enhanced by ∼ 40 % after the addition of (NH 4 ) 2 SO 4 . At around pH 2 and with 0.01 M H 2 C 2 O 4 , the Fe solubility of Libyan dust decreased by ∼ 30 % in the presence of (NH 4 ) 2 SO 4 . Overall, the Fe solubility of Libyan dust was lower compared to that observed in the CFA samples. After 168 h of leaching at pH 2.1, with 1 M (NH 4 ) 2 SO 4 , the Fe solubility of Libyan dust was 7.2 % (Fig. S4), which was from around 3 to 7 times lower compared to that of the CFA samples ( Fig. 1). At around pH 2 conditions in the presence of oxalate and high (NH 4 ) 2 SO 4 concentration, the Fe solubility of Libyan dust rose to ∼ 13.6 % (Fig. S4), which is still 4 times lower than that of Krakow ash and around 1.5 lower than Aberthaw and Shandong ash (Fig. 2). The Fe solubilities of the Libyan dust observed in this study are comparable with those of the Tibesti dust (Tibesti Mountains, Libya; 25.583333 • N, 16.516667 • E) in Ito and Shi (2016) at similar experimental conditions.
The enhanced Fe solubility in CFA compared to mineral dust could be primarily related to the different Fe speciation (Fig. 5). CFA contained more highly reactive Fe and magnetite but less hematite and goethite than mineral dust.
Although mineral dust is the largest contribution to aerosol Fe while CFA accounts for only a few percent, the at- Figure 10. Comparison between the Fe dissolution kinetics of Krakow ash calculated using the original (Ito, 2015) and the new dissolution scheme (Tables 1 and S6). (a-b) Proton-promoted Fe dissolution in H 2 SO 4 solutions with 1 M (NH 4 ) 2 SO 4 at pH 2.1 (a) and at pH 3.0 (b) (Experiment 2 at pH 2.1 and Experiment 2 at pH 3.0 in Table S1). (c-d) Proton + oxalate-promoted Fe dissolution in H 2 SO 4 solutions, with 0.01 M H 2 C 2 O 4 and 1 M (NH 4 ) 2 SO 4 at pH 2.0 (c), and at pH 2.9 (d) (Experiment 3 at pH 2.0 and Experiment 3 at pH 2.9 in Table S1). The Fe dissolution kinetics were predicted using the rate constants in Table 1 calculated in this study (open circles) and the dissolution scheme for combustion aerosols in Ito (2015;cross marks). Note that the dissolution scheme in Ito (2015) was calculated based on laboratory measurements conducted at low ionic strength. (e-f) Contribution of the oxalate-promoted dissolution to dissolved Fe estimated using Eq. (3). The molar concentrations of H 2 SO 4 , H 2 C 2 O 4 , and (NH 4 ) 2 SO 4 in the experiment solutions are shown. The final pH of the experiment solutions is also reported, which was calculated using the E-AIM III model for aqueous solutions (Wexler and Clegg, 2002) and accounting for the buffer capacity of the CFA samples. Figure 11. Percentage contribution of anthropogenic combustion (ANTHRO) aerosol to the atmospheric dissolved Fe concentration near the ground surface from (a) Test 0 and (b) Test 1 for December 2008 and January 2009. In Test 0, we ran the model without upgrades in the Fe dissolution scheme (Ito et al., 2021a) and applied the proton-promoted, oxalate-promoted, and photoinduced dissolution schemes for combustion aerosols in Table S6 (Ito, 2015). The proton + oxalate dissolution scheme (Table 1) was applied in Test 1, and we adopted the base mineralogy for anthropogenic Fe emissions (Rathod et al., 2020). mospheric processing of CFA may result in a larger-thanexpected contribution of bio-accessible Fe deposited to the surface ocean. It is, thus, important to quantify the amount and nature of CFA in atmospheric particles.

Comparison of modelled Fe solubility with field measurements
The model results obtained using the new dissolution scheme for the proton + oxalate-promoted dissolution (Table 1) in Test 1 and Test 3 provided a better estimate of aerosol Fe solubility over the Bay of Bengal than the other tests (Figs. 9 and S3). At the same time, the new model improved the agreement of aerosol Fe solubility from Test 0 (68 % ± 5 %) to Test 1 (35 % ± 2 %) and Test 3 (47 % ± 1 %) with the field data (25 % ± 3 %) but still overestimated it after 22 January 2009, when open biomass burning sources become dominant (Bikkina et al., 2020), as also shown in Fig. 8 and and January 2009. In Test 0, we ran the model without upgrades in the Fe dissolution scheme (Ito et al., 2021a) and applied the proton-promoted, oxalate-promoted, and photoinduced dissolution schemes for combustion aerosols in Table S6 (Ito, 2015). The proton + oxalate dissolution scheme (Table 1) was applied in Test 1, and we adopted the base mineralogy for anthropogenic Fe emissions (Rathod et al., 2020). Table S4. This could be due to the unrepresentative Fe speciation used in Test 1 and Test 3 for biomass burning over the Bay of Bengal. To reduce the uncertainty in model predictions, emission inventories could be improved through a comprehensive characterization of Fe species in combustion aerosol particles. The revised model also enabled us to predict sensitivity to a more dynamic range of pH changes, particularly between anthropogenic combustion and biomass burning by the suppression of the oxalate-promoted dissolution at pH lower than 2. In Test 0, the dissolution rate was assumed to be independent from the pH for extremely acidic solutions (pH < 2). The results show that the proton-promoted dissolution scheme in Test 2 significantly overestimated aerosol Fe solubility (Figs. 9 and S3), which indicates the suppression of the proton + oxalate-promoted dissolution at pH < 2. In Fig. S5, the model estimates of aerosol Fe solubility over the Bay of Bengal considerably improved in Test 1 (RMSE 11) compared to Test 0 (RMSE 21), but more work is needed to improve size-resolved Fe emission, transport, and deposition. The model results in Test 1 indicate a larger contribution of anthropogenic combustion sources to the atmospheric Fe loading over East Asia (Fig. 11) but a smaller contribution of biomass burning sources downwind from tropical regions (Fig. 12). We demonstrated that the implementation of the new Fe dissolution scheme, including a rapid Fe release at the initial stage and highly acidic conditions, enhanced the model estimates. However, in Test 1, we turned off the photoreductive dissolution scheme (Ito, 2015) which was based on the laboratory measurements in Chen and Grassian (2013).
To determine the photoinduced dissolution kinetics of CFA particles, it is necessary to account for the effect of high concentration of (NH 4 ) 2 SO 4 on photo-reductive dissolution rate, which should be considered in future research. Data availability. The new dissolution schemes for the protonpromoted and oxalate-promoted dissolution are reported in Table 1. Table S1 reports the concentrations of H 2 SO 4 , H 2 C 2 O 4 , and (NH 4 ) 2 SO 4 in the experiment solutions, the original and final pH from model estimates (including H + concentrations and activities), and the pH measurements for the solution with low ionic strength. Table S3 contains the summary of the concentration and activity of total oxalate ions, C 2 O 2− 4 , and HC 2 O − 4 in the experiment solutions calculated using the E-AIM III model. The observations of the mass concentration of total Fe, dissolved Fe, and Fe solubility for the fine mode (PM 2.5 ) over the Bay of Bengal are from Bikkina et al. (2020) and are available at https://doi.org/10.1021/acsearthspacechem.0c00063. The modelled mass concentrations of total Fe in aerosol particles and the aerosol Fe solubilities over the Bay of Bengal are reported in Tables S4 and S5, respectively. The Fe speciation, the measurements of the Fe dissolution kinetics, and the results of the IMPACT model for each sensitivity simulation (Test 0-3) can be downloaded at https://doi.org/10.25500/edata.bham.00000702 (Baldo, 2022).
Author contributions. CB, ZS, and AI designed the experiments and discussed the results. ZS supervised the experimental and data analyses. CB conducted the experiments and the data analysis with contributions from ZS, AI, MDK, and NDa. NDa, ZS, and KI performed the XANES measurements. AI developed the model of the dissolution kinetics and performed the model simulations. Krakow and Aberthaw ash were provided by TJ, while Shandong ash was provided by WL. Soil 5 from Libya was collected by NDr. CB prepared the article, with contributions from MDK and all the other co-authors.

Competing interests.
The contact author has declared that neither they nor their co-authors have any competing interests.

Disclaimer.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements. We acknowledge Diamond Light Source for the time on Beamline/Lab I18 under various proposals (SP22244-1, SP12760-1, and SP10327-1). We are grateful to Srinivas Bikkina, Manmohan Sarin, and their colleagues, for kindly providing the observational data set during the SK-254 cruise, supported by the Geosphere-Biosphere programme funded by the Indian Space Research Organization (Bengaluru, India).
Financial support. Clarissa Baldo has been funded by the Natural Environment Research Council (NERC) CENTA studentship (grant no. NE/L002493/1). Support for this research was provided to Akinori Ito by JSPS KAKENHI (grant no. 20H04329) and the Integrated Research Program for Advancing Climate Models (TOUGOU; grant no. JPMXD0717935715) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan.
Review statement. This paper was edited by Ari Laaksonen and reviewed by Rachel Shelley and Morgane Perron.