Pyrogenic iron: The missing link to high iron solubility in aerosols

See allHide authors and affiliations

Science Advances  01 May 2019:
Vol. 5, no. 5, eaau7671
DOI: 10.1126/sciadv.aau7671


Atmospheric deposition is a source of potentially bioavailable iron (Fe) and thus can partially control biological productivity in large parts of the ocean. However, the explanation of observed high aerosol Fe solubility compared to that in soil particles is still controversial, as several hypotheses have been proposed to explain this observation. Here, a statistical analysis of aerosol Fe solubility estimated from four models and observations compiled from multiple field campaigns suggests that pyrogenic aerosols are the main sources of aerosols with high Fe solubility at low concentration. Additionally, we find that field data over the Southern Ocean display a much wider range in aerosol Fe solubility compared to the models, which indicate an underestimation of labile Fe concentrations by a factor of 15. These findings suggest that pyrogenic Fe-containing aerosols are important sources of atmospheric bioavailable Fe to the open ocean and crucial for predicting anthropogenic perturbations to marine productivity.


Marine primary productivity, nitrogen fixation, and particulate organic matter fluxes in the global open ocean are often controlled by the availability of the micronutrient iron (Fe) (1). Primary productivity is most sensitive to the amount of bioavailable Fe input in high-nutrient, low-chlorophyll regions, such as the sub-Arctic and equatorial Pacific, and the Southern Ocean. Atmospheric Fe-containing particles are important sources of bioavailable Fe to the ocean and can thus affect global climate via their influence on primary productivity (2). The solubility of Fe in ambient aerosols (i.e., the ratio of labile to total Fe) largely controls the bioavailability of this micronutrient after its deposition to the ocean (35). The term “labile” Fe in total Fe is operationally defined as the fraction of aerosol Fe that can be leached from the particulate phase into solution, in line with the widely held assumption that this form of Fe is readily bioavailable and therefore affects marine primary productivity (6).

Many factors can potentially affect aerosol Fe solubility near particle source regions and the subsequent enhancement of Fe solubility during atmospheric transport (7). Among these, atmospheric acidity is considered to be a key factor contributing to transformations of aerosol Fe from relatively insoluble to labile forms in the atmosphere (6). Notable progress has been made in our understanding of atmospheric inputs of labile Fe from natural and anthropogenic sources to the surface oceans (8). However, there are still large uncertainties regarding the relative importance of different sources of aerosol Fe and the effects of atmospheric processing on the solubility of the Fe delivered to the surface ocean. Furthermore, the experimental method for the determination of labile Fe in aerosols differs among research groups with regard to the leaching techniques used (e.g., batch versus flow-through), the leaching solutions and extraction times, and the analytical techniques (9). However, these differences do not appear to be the sole cause of the wide range of reported Fe solubilities (10). Generally, an increase in Fe solubility with a decrease in total aerosol Fe concentration appears to be a robust feature in global aerosol datasets (10, 11).

The natural emission of dimethylsulfide from the ocean surface and its oxidation to sulfuric acid was initially hypothesized to be linked to the photochemical transformation of Fe from mineral aerosols in the marine atmosphere (3). Further, the acid mobilization of Fe from mineral dust aerosols was also proposed via sulfuric acid formation from anthropogenic sources during long-range transport (7). The alternative hypothesis points to aerosol particle size as the driver for the enhancement of Fe solubility in dust-dominated aerosols (11), although this effect could not be reproduced in a subsequent study (12).

Atmospheric chemical transport models have considered the dissolution rates for mineral dust from laboratory experiments as a kinetic process, which is dependent on the pH, the mineral composition of Fe-containing particles, organic species (as Fe-reducing agents), and solar radiation (4, 1315). The Fe dissolution scheme was improved on the basis of laboratory studies considering the types of Fe species associated with mineral dust source materials (1620). Following these studies, three Fe pools were characterized as follows: labile Fe (e.g., ferrihydrite and nanosized Fe oxides), less labile Fe (e.g., heterogeneous inclusion of nano-Fe grains in aluminosilicates), and refractory Fe (e.g., crystalline and micrometer-sized iron oxides). The labile Fe fraction in aerosols is regarded as the most readily bioavailable form of Fe and often used as the input from atmospheric models to ocean biogeochemistry models. The fraction of less labile Fe in aerosols may be regarded as potentially bioavailable forms that can be photochemically transformed to labile forms in the atmosphere or after deposition into the ocean. “Less labile Fe” may be operationally distinguished from “labile Fe” by using specific leaching agents (2123), which do not access the refractory Fe fraction (22, 23). Therefore, the labile and less labile content in aerosols may set the upper and lower limits for aerosol Fe solubility, respectively, as it applies to the marine environment.

Elevated levels of Fe solubility (>10%) have been observed in aerosols dominated by combustion sources (2426). The physical and chemical properties of Fe oxides in combustion aerosols (e.g., coal fly ashes) are different from those of naturally occurring minerals due to combustion (27, 28). For instance, ferric sulfate and aggregated nanocrystals of magnetite (Fe3O4) are dominant components of Fe in oil fly ash (29, 30), formed via high-temperature combustion (>800°C) followed by sulfuric acid condensation and potentially modified by filtration of particles from the stack gas, depending on what emission controls are in place (31). Laboratory experiments suggested that Fe in coal fly ash could be mainly labile because of chemical transformations during acidic processing (27). The enhanced Fe solubility of combustion-derived aerosols has been incorporated into the atmospheric Fe transport models in two different ways (32, 33). One of the atmospheric chemical transport models discussed here has adopted a higher initial solubility for Fe sulfate derived from laboratory experiments on oil fly ashes only, and faster dissolution rates for all combustion aerosols, compared to mineral dust (33). Other models only assume higher initial solubilities in combustion aerosols including the oil fly ashes (18, 32).

Given the key role of Fe in ocean biogeochemistry (1) and the importance of atmospheric transport in delivering Fe to the remote oceans, it is critical to develop atmospheric Fe transport models that can realistically simulate total and labile Fe deposition fluxes to the oceans. These models can then provide a tool to consider how Fe deposition fluxes may have changed in the past and/or could change in the future and how anthropogenic perturbations might affect ocean biogeochemistry. Our companion article described four atmospheric Fe models, which participated in the Joint Group of Experts on the Scientific Aspects of Marine Environmental Protection (GESAMP) intercomparison study, and analyzed the differences in atmospheric concentrations and deposition fluxes between models, as well as between models and observations compiled from multiple field campaigns (34). The four models parameterize natural and anthropogenic Fe sources, atmospheric processing, and deposition differently (15, 18, 19, 28). We pointed out that the comparison of monthly mean model results with the shorter-term (e.g., daily) observations during different sampling periods introduced inaccuracies because of the variable nature of high aerosol Fe solubility.

Here, we performed a statistical analysis using the maximum likelihood method to compare Fe solubility between field data and model information (see Materials and Methods). We show better agreement in Fe solubility between model estimates and field data when higher solubility is simulated in pyrogenic Fe than lithogenic Fe for the global ocean and four separate oceanic regions from more to less air polluted conditions. We conclude that high Fe solubility in aerosols is mainly attributed to labile Fe released from pyrogenic Fe oxides. The description of the models and details of the statistical approach are provided in Materials and Methods.


Comparison of simulated Fe solubility to field data over the global ocean

We compiled data for total and labile aerosol Fe concentrations simulated by four models (see Materials and Methods) and from digestions and leaches of aerosol samples collected during multiple field campaigns to improve our understanding of the factors controlling solubility (fig. S1 and tables S1 and S2). We used these data to assess the variability in the solubility of aerosol Fe using a diverse range of methodologies. To reduce the uncertainties associated with source attribution, some simulated data, for which maximum likelihood estimates (MLEs) of total Fe concentrations do not fall within ±2σo of the measurements, are not used for the comparison with field data, as is described in Materials and Methods (fig. S2 and table S2). Here, σo represents the SD of Fe concentration for measurements. This procedure can eliminate the cases that result in a good agreement in labile Fe concentration by compensating the underestimates in Fe solubility with the overestimates in Fe concentration and vice versa, which would bias the relative contribution of combustion and mineral dust sources of Fe to labile Fe. Because direct association of model estimates to a specific cruise track and time period can introduce biases when used as a reference for larger regions and different time periods, the maximum likelihood method was implemented to estimate the state of the atmosphere from observations and models.

The ordinary arithmetic averages of simulated Fe solubility and weighted arithmetic averages of Fe solubility by simulated Fe concentration from the models are compared with the field data in table S3. The weighted mean values are calculated from the cumulated values of labile Fe and total Fe for all data points. The minimum (Min), maximum (Max), number of data points, mean error (Error), root mean square errors (RMSE), mean bias (Bias), and correlation coefficient (r) of Fe solubility between the model results and field data are also listed in table S3. The field data show a wider variability in Fe solubility compared to those in simulated Fe solubility under different conditions, ranging from the minimum (0.02%) to the maximum (98%) with an SD of 8.2 for the arithmetic averages. For each model, the arithmetic average of Fe solubility weighted by the total Fe concentration for all grid points with available observations, and its SD are lower than the corresponding ordinary arithmetic average of the Fe solubilities for the same grid points. This reflects the fact that high Fe concentrations with low solubilities near mineral dust source regions contribute more to the Fe concentrations than others such as combustion aerosols and aged mineral dust aerosols that are sampled over remote regions. The ratios of arithmetic averages to weighted averages for the field data (5.0) and the ratios of their SDs (104) are more consistent with those of the Integrated Massively Parallel Atmospheric Chemical Transport (IMPACT) model when considering both combustion aerosols and mineral dust (5.1 and 115) than those with mineral dust only (1.9 and 40) (table S3). Moreover, better agreement of the model results with the field data of labile Fe concentrations can be seen for IMPACT when considering combustion aerosols plus mineral dust than with mineral dust only (fig. S3 versus fig. S4). At the same time, IMPACT output fits the field data best (correlation coefficient of 0.38 summarized in table S3) and also has the highest solubility estimates for combustion aerosols. This suggests that enhanced and more variable aerosol Fe solubility is associated with mixing of combustion aerosols and its atmospheric processing during transport.

The high Fe solubility (>10%) at low Fe concentration in the data was successfully simulated only by the IMPACT model (Fig. 1). The IMPACT model results suggest that combustion aerosols substantially contribute to labile Fe concentrations (>90%, red color in Fig. 1 and fig. S5). These results support the argument that the source-dependent composition of aerosols is an important factor in controlling the sporadically high Fe solubility (>10%) (25, 29, 35) under polluted conditions in the Northern Hemisphere. On the other hand, the Tracer Model 4 of the Environmental Chemical Processes Laboratory (TM4-ECPL) model assigned high Fe solubility to combustion aerosols but does not simulate high solubility at low concentrations in Fig. 1, mainly because of faster deposition rates for labile Fe than those for total Fe from combustion sources and the absence of the Fe dissolution process for combustion aerosols in TM4-ECPL. Moreover, the IMPACT model prescribed lower initial Fe solubility but simulated high solubility at low concentrations due to Fe released from Fe-containing combustion aerosols during transport. These results suggest that agreements with field data depend not only on the source chemical composition but also on the atmospheric processing during transport. The apparent deviation of Fe solubility in IMPACT from the 1-to-1 ratio is possibly due to the prescribed Fe solubility at emission by the arithmetic mean of field data for oil fly ashes (58 ± 22%) (fig. S5) (29, 30, 36). In addition to the wide range in the initial Fe solubility for oil fly ashes, the Fe dissolution rates for combustion aerosols significantly vary between different samples used for laboratory experiments (27, 30). This suggests that IMPACT may need to resolve the actual variability in the degree of aerosol Fe lability at emission and in the subsequent atmospheric processes to better represent the aerosol source characteristics.

Fig. 1 Atmospheric concentration of total aerosol Fe (ng m−3) versus Fe solubility (%) for simulated estimates (colored circles) and field data (black circles).

(A to E) The color represents the (combustion)/(combustion + dust) ratio for labile Fe concentration in aerosols.

Iron solubility over the northern Indian Ocean

Observations from the northern Indian Ocean provide a valuable dataset to evaluate model performance, as Fe solubilities in aerosols sampled over the Bay of Bengal are considerably higher than those sampled over the neighboring Arabian Sea (Fig. 2 and fig. S1) (26, 37). The Fe solubility correlates linearly with the concentration of acidic species over the Indo-Gangetic Plain, suggesting that acid processing of Fe-containing aerosols influences Fe solubility during long-range transport (38). Moreover, enhanced levels of Fe solubility have been observed in aerosols dominated by biomass burning and fossil fuel combustion sources over the Bay of Bengal (38). IMPACT, TM4-ECPL, and the Community Atmosphere Model version 4 (CAM4), all of which include combustion aerosols, show to some extent the contrast between high Fe solubility in the Bay of Bengal and low solubility in the Arabian Sea, while the Goddard Earth Observing System with Chemistry (GEOS-Chem) model, which does not include combustion aerosols, shows lower Fe solubility over the Bay of Bengal. However, there are large differences regarding the relative importance of different sources of labile Fe (Fig. 2), due to differences in the relative source strength of dust and combustion Fe, and in solubility estimates for dust and combustion aerosols in the models. All models overestimate the low values of the Fe solubility in field data over the Arabian Sea by more than a factor of 10. This is primarily due to the prescribed solubility of labile Fe for dust aerosols in the models at their primary emissions and the subsequent enhancement of Fe solubility estimated from the simulated amount of atmospheric processing during transport in the models. We note that the low solubilities of the field data are outliers (triangles in fig. S5) in the modeled Fe solubility versus the field data. The Fe solubility from the field data used in this study (0.08 ± 0.09%) is significantly lower than those reported by other research groups over the Arabian Sea [1.4 ± 0.5% (39) and 0.4 ± 1.0% (40)]. These differences in the Fe solubility between investigators presumably reflect methodological and analytical differences in the measurements; these uncertainties could potentially be resolved if standard methods and reference material analysis were used by different investigators.

Fig. 2 Comparison of simulated atmospheric Fe solubility with field data (black circles) over the Bay of Bengal and Arabian Sea.

(A to D) The color represents the simulated (combustion)/(combustion + dust) ratio (percent) for labile Fe concentration in aerosols.

Iron solubility over the North Atlantic Ocean

The large body of available field data on aerosol Fe over the North Atlantic offers another useful area in which to compare the models. Air mass back trajectories (AMBTs) were used to classify the source regions and transport pathways for the aerosol samples collected over the North Atlantic (23, 41). However, AMBTs were not sufficiently discriminating to identify the aerosol sources or the potential effects of atmospheric processing on Fe solubility (23). In Fig. 3, the aerosol samples are classified as sourced by air masses from North Africa or other air mass regimes (23). The “degree of atmospheric processing” from the field data can be derived from the ratio of Fe solubility derived from two different leaching methods (see Materials and Methods) (2123). All models calculate that mineral dust was the major source of labile Fe at high total Fe concentrations for samples largely influenced by North African air masses and reproduced the low solubility of 0.4 ± 0.1% of the field data. The internal mixing of alkaline components in mineral dust with Fe-containing minerals can lead to higher pH and thus suppression of Fe dissolution (4, 14, 18). In North African mineral dust aerosols, the ratio of labile Fe to the sum of labile and less labile Fe concentrations indicates a weak degree of Fe processing in aerosols (5.6 ± 2.2% on average from Eq. 3 in Materials and Methods). By contrast, non–North African air mass has a strong degree of Fe processing in aerosols (36 ± 15% on average). At the same time, the three models (IMPACT, TM4-ECPL, and CAM4) estimate a large contribution of combustion aerosols to the labile aerosol Fe at higher Fe solubility and lower aerosol Fe concentrations (Fig. 3). Moreover, IMPACT reproduced the sporadically high Fe solubility observed in aerosols from the other air mass regimes, specifically marine and high-latitude air masses.

Fig. 3 Fe solubility versus Fe concentration (ng m−3) for field data and models over the North Atlantic.

The colors for field data (A) represent percentage contributions of labile Fe to the sum of less labile and labile Fe concentrations in aerosols. The colors for the models (B to D) represent the simulated contribution of combustion aerosols to the sum of combustion and dust aerosols to labile Fe. The field data in light gray color represent no data for the 25% acetic acid leach. The aerosol samples are classified as sourced by two different air masses from North Africa (triangles) or other air mass regimes (circles) including North America, Europe, marine (no interaction with major land masses within 5 days before sample collection), or high latitudes (air masses originating north of 50°N) (23).

Iron solubility over the equatorial Atlantic Ocean

The simulated Fe solubility is further compared with the field data (black symbols) in the equatorial Atlantic (Fig. 4). The atmospheric residence time, and thus the time available for atmospheric processing, of North African mineral dust aerosols increases as the dust is transported to the west, but there is no significant trend in Fe solubility for the field data with longitude, around 10°N across the tropical Atlantic (42). This low solubility of 2.0 ± 1.5% with no longitudinal trend for aerosol Fe is captured by IMPACT, TM4-ECPL, and GEOS-Chem.

Fig. 4 Comparison of simulated atmospheric Fe solubility (colored points) with field data (black points) over the equatorial Atlantic.

(A to D) The color represents the simulated (combustion)/(combustion + dust) ratio (percent) for labile Fe concentration in aerosols.

Iron solubility over the Southern Ocean

Observations of Fe solubility near coastal areas may be affected by regional anthropogenic emissions, whereas those over the remote ocean may be influenced more by atmospheric processing during long-range transport. Over the Southern Ocean (Fig. 5), none of the models was able to reproduce the high Fe solubility (>10%) occasionally reported for aerosols (12 of 42 samples), because there is no aerosol Fe dissolution mechanism under high-pH and low-oxalate conditions in IMPACT, TM4-ECPL, and GEOS-Chem. Additionally, the high Fe solubility of the field data in the Southern Ocean exceeds the sum of labile and less labile solubility (6.0 ± 1.0%) for North African dust (23) and the Fe solubility in glacial weathering products (2 to 3%) (29). Thus, such a wide range of solubility (from 0.2 to 48%) in this region as derived from observations cannot be explained by considering only chemical aging of mineral dust aerosol as measured in the laboratory. Different types of Fe-containing minerals such as pyrogenic Fe oxides must be considered in conjunction with additional aerosol processing such as plume chemistry to achieve these high solubilities. Moreover, averaged labile Fe concentration from all modeled estimates (0.065 ± 0.078 ng m−3) underestimates the observed mean over the Southern Ocean (0.967 ± 1.049 ng m−3) by a factor of 15.

Fig. 5 Fe solubility versus Fe concentration (ng m−3) for field data (black squares), IMPACT (red squares), TM4-ECPL (orange squares), CAM4 (purple squares), and GEOS-Chem (blue squares) over the Southern Ocean (>45°S).

The vertical lines in symbols correspond to ±1 SD for Fe solubility in aerosols.

We therefore suggest some mechanisms to explain these high solubilities. Pyrogenic Fe, which can be transformed to labile Fe during atmospheric processing, may play an important role in providing labile Fe to the Southern Ocean. However, further studies are required to identify and quantify such aerosol sources in the Southern Hemisphere (4345). The higher solubility may originate from aerosol interactions with emissions from open biomass burning, which is common in the Southern Hemisphere (43, 44). Additionally, continental Fe-containing aerosols have longer lifetimes over remote oceanic regions in the Southern Hemisphere, such as in the outflow of South America, southern Africa, and Australia (34). Thus, aerosols originating from land can experience atmospheric processing over longer periods than in the Northern Hemisphere. An additional possibility is that some Fe dissolution processes may be missing from the models. For instance, dissolution of nanosized colloidal Fe oxides in ice is not taken into account in our models (46). Mineral dust in Antarctica snow samples was shown to have low Fe solubility (~0.7%) (47), which was consistent with the IMPACT model estimates, while the Fe-containing aerosols affected by fires may be associated with sporadic high Fe solubility (44). Thus, we recommend ice-enhanced dissolution experiments of representative combustion aerosols (e.g., biomass burning aerosols and volcanic ashes) under conditions relevant to ice particles in the atmosphere over the Southern Ocean. Studies of stable Fe isotopes may offer additional constraints on the contribution from combustion sources to aerosol Fe (48).

Uncertainty in aerosol Fe solubility stems from the scarcity of observations in the Southern Hemisphere. There are fewer field data in the Southern Hemisphere, which are sparsely distributed and collected over longer sampling periods during research cruises, partly due to low concentration of aerosols. Thus, these results should be interpreted with a degree of caution. We suggest that more observations are needed before drawing any firm conclusions on the Fe solubility and the models’ performance in the Pacific, Indian, and Southern oceans. The high solubility in regions of low Fe concentration (<100 ng m−3) was thought to be a robust feature of all the available data regardless of the leaching methods. However, low Fe solubilities were also found at low Fe concentrations over the Southern Ocean (Fig. 5) in some recent field data (10, 43, 49, 50). Given the key role of the Southern Ocean in the Earth system, it is important to gain a better understanding of pyrogenic and lithogenic Fe sources in the region alongside improved understanding of atmospheric processing of Fe-containing aerosols to resolve these discrepancies between model estimates and field data.

Simulated contribution of combustion aerosols to labile Fe

The contribution of combustion aerosols to total Fe deposition flux to the ocean, as calculated by the Fe deposition fluxes originating from dust and combustion sources (figs. S6 and S7), is shown in Fig. 6. Three models (IMPACT, TM4-ECPL, and CAM4) can calculate the contribution from combustion aerosols to total Fe and to labile Fe deposition fluxes. The models agree that mineral dust is the major source of Fe deposited to the North Atlantic, Arabian Sea, and South Atlantic downwind of the arid and semiarid regions of North Africa, the Middle East, and Patagonia (fig. S6). However, the models estimate a large range in the contribution of combustion aerosols to the deposition fluxes of total (from 2.6 to 7.1% on a global mean) and labile Fe (from 9 to 44% on a global mean) (Fig. 6). For total Fe, this range in estimates could be due to differences in the relative emission strengths of dust versus combustion-derived Fe (fig. S6), the mineralogical composition at the source region, the size distribution of airborne dust particles, and the meteorological data in the models (34). Because minerals in soils differ in their Fe content, spanning from 0.2 to 70%, Fe content in aerosols depends on the database of mineral composition in clay- and silt-sized soils, and the global mean Fe content in mineral dust emissions ranges from 2.6 to 3.5% (34).

Fig. 6 Proportion of pyrogenic Fe in total aerosol Fe and labile Fe.

Percentage of total aerosol Fe from combustion aerosols (A, D, and G) and the percentage of labile Fe from combustion aerosols to the total atmospheric Fe deposition flux (combustion and dust) calculated from three models with atmospheric processing (B, E, and H) and without atmospheric processing (C, F, and I).

The differences in the contribution of combustion aerosols to the total deposition fluxes among model results are significantly larger for labile Fe (22% for IMPACT, 44% for TM4-ECPL, and 9.3% for CAM4) than for total Fe (2.6% for IMPACT, 7.1% for TM4-ECPL, and 3.1% for CAM4) and increase with distance from the major dust source regions (Fig. 6). This is largely due to the differences in the assumed solubility of dust and combustion aerosols upon emission (table S1) and the subsequent enhancement of Fe solubility dynamically calculated from the simulated amount of atmospheric processing during transport. To show the sensitivity of simulated Fe solubility to atmospheric processing, the initial solubilities of 0.45% for dust and 4% for combustion aerosols (32) are prescribed with no atmospheric processing in Fig. 6 (C, F, and I). When atmospheric processing is not considered to enhance Fe solubility, the IMPACT model shows lower contribution of combustion aerosols than that with atmospheric processing over the remote ocean, because the faster dissolution schemes at lower pH than those at higher pH primarily work for wet combustion aerosols during atmospheric transport. The IMPACT model results are reasonably consistent with observations, which indicate much more increases in aerosol Fe solubilities over the open ocean from those near sources of biomass burning and coal combustion (24, 26, 43), compared to that for mineral dust during atmospheric transport (23, 42). When atmospheric processing is neglected, TM4-ECPL model shows higher contribution of combustion aerosols over the Southern Ocean, because simulated labile Fe over the ocean is due only to primary emissions (initial solubility) and labile Fe is removed faster than total Fe by wet deposition during atmospheric transport. Therefore, reducing uncertainties in model estimates of labile Fe deposition fluxes primarily requires improved parameterization for the factors that control aerosol Fe solubility (table S1) and the relative contribution of combustion aerosols to the labile fraction of aerosol Fe (Fig. 6).


The models need to consider the emissions of pyrogenic Fe-containing aerosols and their atmospheric transformation to simulate the high Fe solubility at low total aerosol Fe concentrations in the Northern Hemisphere. We found that the models did not reproduce the high Fe solubility over the Southern Ocean, suggesting potentially missing sources of pyrogenic Fe-containing aerosols and their dissolution processes. Further investigations of processes enhancing Fe solubility are needed in aerosols, rain (liquid cloud), and snow (ice cloud) over the oceans. A standard analytical method for assessing aerosol Fe solubility for polluted and clean environments is a high priority for future studies. Furthermore, our study suggests that assessments of atmospheric fluxes of labile Fe to the open oceans should consider Fe-containing aerosols emanating from both pyrogenic and lithogenic sources for accurately predicting anthropogenic perturbations to oceanic nutrient biogeochemistry.


Global atmospheric models

The global atmospheric models used in this study are summarized in table S1 (15, 18, 19, 28). The models participating in the intercomparison studies, within the framework of the GESAMP Working Group 38, differ in how the Fe emission and dissolution schemes are coupled and in the complexity of those schemes such as the inclusion or exclusion of combustion aerosols (34). Neither meteorological conditions nor emission inventories in these simulations were prescribed to a specific year. We bias-corrected the size-resolved mineral dust loading and deposition fluxes calculated by each model using recently obtained constraints on size-resolved dust loading (34, 51).

Fe solubility from mineral dust in the four models (and for combustion aerosols in IMPACT and CAM4) increases with atmospheric processing. The amount of dissolved Fe is primarily determined by two factors: (i) content of labile Fe and (ii) dissolution rates of less labile Fe. Only a fraction of the less labile Fe can be dissolved during aerosol lifetimes in the atmosphere. Laboratory studies have demonstrated that the dissolution rate of Fe from minerals is strongly dependent on the proton concentration in bulk solution, the mineral surface concentration of oxalate, solar radiation, and the mineralogy (16, 20, 27, 30). IMPACT, TM4-ECPL, and GEOS-Chem involved thermodynamic equilibrium modules to estimate the acidity in the aqueous phase of hygroscopic particles (7, 17, 18). On the other hand, CAM4 used two pH values (i.e., pH = 2 or pH = 7.5), which depended on the simulated ratio of calcite to sulfate (13). Note that coating of aerosols with inorganic and organic matter, and restructuring of the aggregates during the dispersion of the plume released by the combustion sources, is part of the so-called plume chemistry and is a sub–grid scale phenomenon in coarse-scale models. Thus, models that prescribed the Fe solubility for combustion aerosols at the emission level (table S1) implicitly assumed that atmospheric processing of aerosols also produced labile Fe within each grid box at the emission level. Therefore, the emissions treated in the models as primary sources also included secondary sources occurring within the large grid box, which could not be separated into primary emission and rapid secondary formation in the atmosphere. The Fe solubilities were derived from the sum of mineral dust and combustion aerosols in IMPACT, TM4-ECPL, and CAM4, but only for mineral dust in GEOS-Chem.

Constraining estimates of Fe concentration

The daily averages of Fe concentration and solubility at the surface from four models were used for a comparison with the ambient measurements over the oceans (fig. S1 and table S2) (10, 21, 23, 24, 26, 37, 4143, 49, 50, 5259). When the field data uncertainties were not available, regional averages of relative SDs (SD divided by the average number of replicates) were used to estimate the SDs of the measurements. Mineral dust and biomass burning emissions are highly episodic and exhibit seasonal and interannual variability (5). It is, therefore, problematic to compare the monthly mean values of model estimates with the field data from different times of year, or different years. Moreover, the field observations have shown relatively high solubility at low concentration in aerosols influenced by combustion sources (2426). Thus, aerosols at high atmospheric concentration can be more representative of mineral dust aerosols, while mixed aerosols from different sources, i.e., dust and combustion aerosols, may be more representative at low atmospheric concentration. This suggests that modeled Fe concentrations should be consistent with the measurements to assess the relative contribution of each source. However, direct association of model estimates to a specific cruise track and time period can introduce biases when used as a reference for larger regions and different time periods. The maximum likelihood method has been operationally implemented to atmospheric data assimilation systems to estimate the state of the atmosphere from observations and models. Therefore, the MLEs of Fe concentration, sa, with error variance, σa2, were derived from the measurements, so, and a priori daily estimates modeled in the same month as the measurements, sb, with corresponding error variances, σo2 and σb2,σb2σo2+σb2 × (sosb)(1)σa2=σo2σb2σo2+σb2(2)

Using the probability function of the MLEs of Fe concentration assuming a normal distribution (MATLAB function), the weighted averages and SDs of Fe concentration and solubility were calculated. The a posteriori data were used as simulated estimates for the comparison with the measurements when MLEs of Fe concentration fell within ±2σo of the measurements (fig. S1 and table S2).

Degree of atmospheric processing of Fe

Fe-containing aerosols from combustion sources are characterized by high solubility in solution (i.e., leaching media), compared to mineral dust, because the chemical form of Fe in the combustion aerosol may be inherently more labile (e.g., ferric sulfate and nanocrystals of Fe oxides). Because the aerosol water associated with combustion aerosols can be very acidic (pH generally below 3) (60) in the presence of enough organic ligands, Fe dissolution is significantly facilitated via proton-promoted, ligand-promoted, and photoreductive mechanisms, and labile Fe is maintained in solution, resulting in relatively high solubility for combustion aerosols (27, 33). Therefore, the contribution of labile Fe to the sum of labile and less labile Fe may be viewed as the degree of atmospheric processing before deposition into the ocean. To assess the degree of atmospheric photochemical processing of Fe from the less labile forms, the ratio of Fe solubility estimates from two leaching methods (ultrahigh-purity water leach and 25% acetic acid leach) (23) was used according to the following equation(Degree of atmospheric processing of Fe)=[Fe]water leach[Fe]water leach+[Fe]acid leach(3)

The pH value of the 25% acetic acid leach (pH 2.1) is consistent with the pH estimates of water surrounding deliquescent aerosols using a thermodynamic equilibrium model at a rural location in the southeastern United States during 1998–2013 (60).


Supplementary material for this article is available at

Fig. S1. Fe solubility (%) of aerosols collected during multiple field campaigns for the field data and the simulated estimates from models used in this study.

Fig. S2. Comparison of simulated estimates and observed values for total Fe concentration (ng m−3).

Fig. S3. Comparison of simulated estimates with combustion aerosols plus mineral dust and observed values for labile Fe concentration (ng m−3).

Fig. S4. Comparison of simulated estimates with mineral dust only and observed values for labile Fe concentration (ng m−3).

Fig. S5. Model estimates of Fe solubility in aerosols versus field data collected during multiple field campaigns (triangles for Arabian Sea and circles for others).

Fig. S6. Deposition fluxes of total Fe (ng Fe m−2 s−1) from dust and combustion sources to present-day oceans.

Fig. S7. Deposition fluxes of labile Fe (ng Fe m−2 s−1) from dust and combustion sources to present-day oceans.

Table S1. Summary of averaged Fe solubility in atmospheric models used in this study.

Table S2. Latitude, longitude, month, and Fe solubility (%) in aerosols for the intercomparison study.

Table S3. Comparison of Fe solubility (%) between the simulated estimates from models and field data.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We are grateful to T. Wagener for providing his field dataset. W.M.L. would like to acknowledge the GEOTRACES program for the North Atlantic aerosol samples. Funding: Support for this research was provided to A.I. by JSPS KAKENHI grant number JP16K00530 and Integrated Research Program for Advancing Climate Models (MEXT). S.M. acknowledges financial support by the Marie-Curie H2020-MSCA-IF-2015 grant ODEON (ID 705652). Part of model data post-processing was carried out on the Dutch national e-infrastructure with the support of the SURF Cooperative. N.M.M., R.A.S., and D.S.H. acknowledge support from DOE DE-SC0006791. This work was supported by grants to J.F.K. (NSF-CLD-1552519), W.M.L. (NSF-OCE-0929919 and NSF-OCE-1034764), and Y.G. (NSF-OPP-0944589). Support to Y.F. is provided by Argonne National Laboratory under the U.S. Department of Energy contract DE-AC02-06CH11357. This paper resulted from the deliberations of United Nations GESAMP Working Group 38, “The Atmospheric Input of Chemicals to the Ocean.” We thank the ICSU Scientific Committee on Oceanic Research, the U.S. NSF, the Global Atmosphere Watch of the World Meteorological Organization, the International Maritime Organization, and the University of East Anglia for their support. Author contributions: This study was initiated by participants in a GESAMP Working Group 38 workshop in 2017 (A.I., S.M., M.K., N.M.M., A.R.B., T.J., M.S., Y.G., R.U.S., M.M.G.P., C.G., A.N., and R.A.D.). A.I. conceived the statistical analysis and analyzed the data. S.M. manipulated model fields to a common format. A.I., S.M., M.K., N.M.M., R.A.S., D.S.H., N.M., and M.S.J. contributed model products. A.R.B., T.J., S.M., S.B., Y.G., R.U.S., C.S.B., W.M.L., A.R.B., C.G., and Y.F. contributed measured data. J.F.K. provided the analysis for dust correction factors for model data. All authors contributed to writing the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article