Atmospheric processing of iron in mineral and combustion aerosols: development of an intermediate-complexity mechanism suitable for Earth system models

Atmos. Chem. Phys., 18, 14175–14196, 2018 https://doi.org/10.5194/acp-18-14175-2018 © Author(s) 2018. This work is distributed under the Creative Commons Attribution 4.0 License. Atmospheric processing of iron in mineral and combustion aerosols: development of an intermediate-complexity mechanism suitable for Earth system models Rachel A. Scanza1,2, Douglas S. Hamilton1, Carlos Perez Garcia-Pando3, Clifton Buck4, Alex Baker5, and Natalie M. Mahowald1 1Department of Earth and Atmospheric Sciences, Cornell University, Ithaca, New York, USA 2Atmospheric Sciences and Global Change Division, Pacific Northwest National Laboratory, Richland, Washington, USA 3Earth Sciences Department, Barcelona Supercomputing Center, Barcelona, Spain 4Department of Marine Sciences, University of Georgia, Athens, Georgia, USA 5School of Environmental Sciences, University of East Anglia, Norwich, UK Correspondence: Rachel A. Scanza ([email protected]) Received: 25 January 2018 – Discussion started: 19 February 2018 Revised: 30 June 2018 – Accepted: 28 August 2018 – Published: 5 October 2018 Abstract. Atmospheric processing of iron in dust and combustion aerosols is simulated using an intermediatecomplexity soluble iron mechanism designed for Earth system models. The solubilization mechanism includes both a dependence on aerosol water pH and in-cloud oxalic acid. The simulations of size-resolved total, soluble and fractional iron solubility indicate that this mechanism captures many but not all of the features seen from cruise observations of labile iron. The primary objective was to determine the extent to which our solubility scheme could adequately match observations of fractional iron solubility. We define a semiquantitative metric as the model mean at points with observations divided by the observational mean (MMO). The model is in reasonable agreement with observations of fractional iron solubility with an MMO of 0.86. Several sensitivity studies are performed to ascertain the degree of complexity needed to match observations; including the oxalic acid enhancement is necessary, while different parameterizations for calculating model oxalate concentrations are less important. The percent change in soluble iron deposition between the reference case (REF) and the simulation with acidic processing alone is 63.8 %, which is consistent with previous studies. Upon deposition to global oceans, global mean combustion iron solubility to total fractional iron solubility is 8.2 %; however, the contribution of fractional iron solubility from combustion sources to ocean basins below 15◦S is approximately 50 %. We conclude that, in many remote ocean regions, sources of iron from combustion and dust aerosols are equally important. Our estimates of changes in deposition of soluble iron to the ocean since preindustrial climate conditions suggest roughly a doubling due to a combination of higher dust and combustion iron emissions along with more efficient atmospheric processing. 1 Introduction Nearly all ocean productivity relies on enzymatic iron for photosynthesis, respiration and nitrogen fixation. The iron biogeochemical cycle is therefore an important modulator of the oceans’ ability to uptake atmospheric CO2. Since approximately 30 % of atmospheric CO2emissions are absorbed by the oceans (Seiler and Crutzen, 1980; Broecker et al., 1979; Siegenthaler and Sarmiento, 1993), changes to the iron cycle may induce a potentially large negative feedback on the climate system (Martin and Fitzwater, 1988; Martin, 1990; De Baar et al., 1995; Jickells and Spokes, 2001; Jickells et al., 2005). Many remote ocean regions are characterized by highnutrient low-chlorophyll (HNLC) conditions; in the equatorial Pacific and Southern Ocean basins, ocean primary productivity (OPP) is iron limited (Martin and Fitzwater, 1988; Boyd et al., 2000). A significant source of new iron (in conPublished by Copernicus Publications on behalf of the European Geosciences Union. 14176 R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols trast to recirculated iron from upwelling, estuary runoff or fluxes from continental margins) to these regions is atmospheric deposition of desert dust (Jickells et al., 2005; Fung et al., 2000); however, iron in dust is typically insoluble and not considered to be bioavailable for OPP (Mahowald et al., 2005; Johnson, 2001; Jickells and Spokes, 2001). Although combustion iron is a much smaller source of new iron, it is considerably more soluble than dust iron (Guieu et al., 2005; Chuang et al., 2005), potentially contributing up to 50 % of bioavailable iron in remote HNLC regions (Luo et al., 2008; Mahowald et al., 2009; Winton et al., 2015; Ito, 2015). Many definitions exist for bioavailable iron although how organisms utilize iron is not well understood (Jickells et al., 2005). Bioavailable iron is typically considered to be dissolved iron in the ferrous (Fe(II)) oxidation state; however, ocean organisms have been observed to utilize iron in different forms (Barbeau et al., 1996). In this study, bioavailable iron is defined as labile iron, or dissolved iron in either the ferrous or ferric oxidation states. Acidic processing of mineral dust and combustion aerosols during atmospheric transport is one potential mechanism for iron dissolution upon deposition, and many studies have observed increased iron liberation from insoluble iron oxides with decreasing pH (Duce and Tindale, 1991; Zhu et al., 1997; Zhuang et al., 1992; Jickells and Spokes, 2001; Desboeufs et al., 2001; Meskhidze et al., 2003). However, acidic processing is likely to work in conjunction with other physical and chemical processes such as photochemical reduction and organic ligand processing. Iron solubility has also been correlated with particle size with increasing fractional iron solubility as mineral dust concentrations decrease (Baker and Croot, 2010; Chen and Siefert, 2004; Baker and Jickells, 2006). This relates to the preferential settling of larger particles during transport along with a higher surface-area-to-volume ratio for smaller particles. Because smaller particles are associated with longer transport times, their probability of undergoing multiple cycles of evaporation and condensation in the atmosphere is increased. Oxalate, the oxidation product of oxalic acid, can act as an organic ligand; this has been observed to enhance iron dissolution when complexation with aerosols occurs in aqueous solutions under moderately acidic conditions (e.g., within clouds) (Cornell and Schindler, 1987; Xu and Gao, 2008; Solmon et al., 2009; Paris et al., 2011). The formation of oxalate in the atmosphere is complicated (Myriokefalitakis et al., 2011) but may increase soluble iron concentrations over 75 % in addition to acidic processing alone (Johnson and Meskhidze, 2013; Myriokefalitakis et al., 2015; Ito, 2015). Previous studies focused on observations have shown a very strong inverse relationship between total iron and solubility (Sholkovitz et al., 2012), where differences in methods of measuring soluble iron appear to be much less important than previously thought (Baker et al., 2014). Observations from Sholkovitz et al. (2012) show a trend of increasing iron solubility with decreasing total iron concentration and conclude that much of this is due to the smaller particle size distribution of combustion aerosols and their ability to enhance dust iron dissolution. Mahowald et al. (2018) show using a simple 1-D plume model that either differential solubility in emissions of combustion aerosols or atmospheric processing of dust iron can match the observed relationship described above. At emission, combustion aerosols are typically smaller and thus remain suspended in the atmosphere longer than the larger dust aerosols. In addition, since combustion aerosols appear to have higher Fe solubility, the plume will increase in solubility as it is transported downwind from continental sources. Alternatively, if only dust aerosols are considered to be sources of iron, but are atmospherically processed downwind from source regions, one will see lower values of iron in the aerosols, as particles fall out, yet the remaining iron will be more soluble (longer atmospheric processing). Quantitative comparison of this simple plume model indicates that either of these two cases can match observations (Box 2, Mahowald et al., 2018). Panel (b) in Box 2 (Mahowald et al., 2018) is reproduced here (Fig. S5 in the Supplement) on a linear scale to facilitate a comparison to Fig. 3 in Sholkovitz et al. (2012). While the combustion case (blue) is the only case able to replicate very high aerosol Fe solubility (>30 %), all three cases match the observed trend of increasing fractional iron solubility with decreasing total aerosol iron concentrations. The type of iron in dust is dependent on mineralogy and surface area. For example, iron associated with hematite and goethite is tightly bound as iron (hydr)oxides and its dissolution even under very acidic conditions is slow. Nanoparticles of hematite and ferrihydrite have been observed to coat clay minerals illite and smectite (Shi et al., 2012; Ito and Xu, 2014); the surface-area-to-volume ratio in nanoparticles (<0.1 µm) is sufficiently high to facilitate rapid dissolution of Fe at low pH. The clay minerals illite, smectite, feldspar and kaolinite also include Fe as substitutions within aluminosilicate mineral lattices and may be appreciable sources of soluble Fe (Journet et al., 2008). Different levels of model complexity have been employed to simulate atmospheric iron processing from very simple models including a first-order rate constant applied to a constant 3.5 % iron in dust (Hand et al., 2004; Luo et al., 2008) to more complex schemes allowing different types of acidic species to interact with mineral dust and combustion aerosols during transport and applying mineral-specific dissolution rates (Meskhidze et al., 2005; Johnson and Meskhidze, 2013; Ito and Xu, 2014; Myriokefalitakis et al., 2015; Ito, 2015). A more recent study developed an iron dissolution scheme with fewer mineral tracers to allow for simulations using Earth system models (Ito and Shi, 2016). In this study, our goal is to develop an iron dissolution scheme of intermediate complexity that can be used in Earth system models. Recent work has emphasized the importance of variable dissolution rates based on pH and chemical composition (Meskhidze et al., 2005; Johnson and Meskhidze, Atmos. Chem. Phys., 18, 14175–14196, 2018 www.atmos-chem-phys.net/18/14175/2018/ R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols 14177 2013; Ito and Xu, 2014; Myriokefalitakis et al., 2015; Ito, 2015) and this requires the advection of multiple chemical species and calculations of chemical equilibrium. We include these complex mechanisms but use simplified parameterizations based on standard aerosol species available in most Earth system models. The type and amount of iron at emission is determined from mineral maps and brittle fragmentation theory (Claquin et al., 1999; Scanza et al., 2015; Kok, 2011). We simulate four types of iron, three for mineral dust aerosols (readily released Fe, moderately soluble Fe, weakly soluble Fe) and one for combustion aerosols (combustion Fe) (Ito and Xu, 2014; Ito, 2015). An acid-processing scheme is employed based on atmospheric concentrations of sulfate along with an in-cloud oxalate mechanism derived from laboratory data from Paris et al. (2011) (Johnson and Meskhidze, 2013; Myriokefalitakis et al., 2015; Ito, 2015). Our method performs as well as more computationally intensive methods in matching available observations of soluble iron. Additionally, the parameterizations and simplifications incorporated in our scheme render it computationally efficient, allowing for multi-decadal to centennial simulations. The study is organized as follows: Sect. 2 describes the climate and dust model we use along with the details of our iron-processing scheme. Section 3 compares the results of the model simulations to available observations of soluble iron using a set of metrics established to evaluate model accuracy. Additionally, results from six sensitivity studies are discussed in order to quantify which model parameterizations are required to best match observations. Section 4 presents the results from four preindustrial simulations. Section 5 includes a discussion and comparison to previous modeling studies. The final section describes where future work should be directed. 2 Methods The Community Earth System Model (CESM) coordinated by the National Center for Atmospheric Research (NCAR) is a global Earth system model used for past, present and future climate simulations. In this study we use the CESM version 1.0.5 with a bulk aerosol model, the Community Atmosphere Model version 4 (CAM4) (Neale et al., 2013). Previous versions of the model have been modified to allow for improved treatment of dust and mineral speciation (Albani et al., 2014; Scanza et al., 2015; Zhang et al., 2015), and here we further modify CAM4 in order to simulate the emission, transport and processing, and deposition of total and soluble iron in desert dust and combustion aerosols. 2.1 Desert dust model CAM4 is used to simulate the emission, transport and deposition of bulk aerosols in the CESM. Major aerosol species include dust, sea salt, black carbon, organic carbon (OC) and sulfate, and they are prescribed as external mixtures. Simulations are performed at a horizontal resolution of 1.9◦×2.5◦ with 56 vertical levels that span the pressure at the surface to around 2 hPa (Computational and Information Systems Laboratory, 2012; Yellowstone: IBM iDataPlex System, Climate Simulation Laboratory; Boulder, CO: National Center for Atmospheric Research, http://n2t.net/ark:/85065/d7wd3xhc, last access: 28 September 2018). Meteorology is driven by NASA’s Goddard Earth Observation System (GEOS-5) (Suarez et al., 2008; Hurrell et al., 2013; Lamarque et al., 2012). The model is run from 2005 through 2011 assuming present-day climate conditions; we use the last five years (2007–2011) for analysis. The desert dust model in CAM4 is a modified version of the Dust Entrainment and Deposition model (DEAD) (Zender et al., 2003) and is exactly as described and used in Albani et al. (2014) and Scanza et al. (2015). Briefly, the emission of dust is computed with an analytic trimodal lognormal probability density function from three source modes to four transport bins (Zender et al., 2003). The three source modes are defined by observed mass median diameter and geometric standard deviation (d’Almeida, 1987; Schulz et al., 1998; Zender et al., 2003). Bin widths are prescribed at 0.1–1.0, 1.0–2.5, 2.5–5.0 and 5.0–10.0 µm (Mahowald et al., 2006; Zender et al., 2003) and have fixed lognormal sub-bin distributions (mass median diameter =3.5 µm, σg=2.0) (Zender et al., 2003). The size distribution across dust bins was modified from the release version of the model to follow the brittle fragmentation theory of vertical dust flux (Kok, 2011) with particle mass fractions of 0.011, 0.087, 0.272 and 0.625 in bins 1–4, respectively (Albani et al., 2014). Dry removal of dust aerosols involves parameterizations for gravitational settling and turbulent mix out and wet removal includes incloud and below-cloud scavenging (Rasch et al., 2000; Zender et al., 2003; Mahowald et al., 2006). The in-cloud aerosol removal rate is proportional to the fraction of cloud water that precipitates; some of the falling precipitation is allowed to re-evaporate (Rasch and Kristjánsson, 1998; Rasch et al., 2000). CAM4 allows for multiple cycles of condensation and evaporation (Gent et al., 2011; Hurrell et al., 2013; Neale et al., 2013) in order to match observational estimates of approximately three in-cloud cycles for aerosols (Lelieveld et al., 1998; Crutzen and Zimmermann, 1991); the model time step is 30 min with 20 sub-time steps for in-cloud chemistry. Modifications from the release version include changes to the scavenging coefficients (from 0.1 to 0.3 for bins 3 (2.5– 5.0 µm) and 4 (5.0–10.0 µm)), to the prescribed particle solubility from 0.15 to 0.3 across all size bins, and to the threshold for dust generation of the leaf area index from 0.15 up to 0.3 (Lancaster and Baas, 1998; Mahowald et al., 2006; Okin, 2008; Albani et al., 2014; Scanza et al., 2015). The dust module was tuned to best match aerosol optical depth (AOD), deposition and surface concentration data (Albani et al., 2014). Similar to other models (Huneeus et al., 2011), our model is unable to simultaneously match the surface concentration www.atmos-chem-phys.net/18/14175/2018/ Atmos. Chem. Phys., 18, 14175–14196, 2018 14178 R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols and deposition data in remote regions (Fig. 3 in Albani et al., 2014), and since here we only show concentration data, it will appear the model over-predicts the surface concentration of dust. 2.2 Mineralogy Mineralogy maps are derived from (Claquin et al., 1999) and are gridded using FAO/UNESCO WGB84 at 50×50arcmin with soil legend 184 from FAO/UNESCO Soil Map of the World (1976; file identifier: f7ccd330-bdce-11db-a0f6-185000d939bc5d8) (Batjes, 1997). The conversion from soil to aerosol mineralogy for each transport bin follows the brittle fragmentation theory of dust (Kok, 2011) and is described in detail in Scanza et al. (2015). The mineralogy maps are based on measurements performed following wet sieving, a process that destroys the mineral aggregates found in unperturbed parent soils. The brittle fragmentation theory reconstructs the mineral aggregation found at emission, allowing for more accurate reproduction of size-resolved dust and iron concentrations (Perlwitz et al., 2015a, b; Pérez García-Pando et al., 2016; Ito and Shi, 2016). We determine iron concentrations at emission from the derived aerosol mineralogy; therefore, we eliminate the need to explicitly advect individual mineral tracers with the exception of calcite in order to parameterize the buffer effect of mineral dust on aerosol pH (Böke et al., 1999). Evaluation of the mineral distributions against observations of dust plume events showed some skill in the simulation of concentrations of minerals (Scanza et al., 2015); in addition, the resulting elemental distributions of iron, aluminium and calcium were improved through the use of the soil map (Zhang et al., 2015). 2.3 Determination of iron Iron in mineral dust is found in several forms and the solubility depends on the type of chemical bonding (or compositional form) (i.e., iron oxides vs. aluminosilicate inclusions) (Lafon et al., 2004; Journet et al., 2008) and on the particle size (Baker and Jickells, 2006). The type of iron is mineral dependent: aluminosilicates contain Fe inclusions as well as surface coatings of nanohematite and ferrihydrite (Shi et al., 2009, 2012; Ito and Xu, 2014). Fe associated with hematite and goethite is tightly bound to oxygen in iron oxides or hydroxides, is typically larger in size and is less soluble. Following previous studies, we define the Fe concentration in mineral dust as follows: 57.5 % Fe in hematite, 11 % in smectite, 4 % in illite, 0.24 % in kaolinite and 0.34 % in feldspar (Journet et al., 2008; Ito and Xu, 2014). Three types of dust Fe are simulated in the model: readily released iron (FeRR), moderately soluble iron (Femed), and weakly soluble or refractory iron (Feslow). The FeRR is assumed here as 2.7 % of the 4 % Fe in illite as ferrihydrite, 5 % of the 11 % Fe in smectite as nanohematite, 4.3 % of the 0.24 % Fe in kaolinite and 3 % of the 0.34 % Fe in feldspar. Femed is defined as the remaining 97.3 % and 95 % Fe in illite and smectite, respectively. Feslow is defined as the remaining 95.7 % Fe in kaolinite, 97 % in feldspar and 100 % of the 57.5 % Fe in hematite (Table S1 in the Supplement). Separate tracers are defined for total moderately soluble iron (Fetmed) and weakly soluble iron (Fetslow) in each of the four CAM4 size bins along with separate tracers for moderately and weakly soluble iron (Fesmed and Fesslow). In the finest particle size bin (0.1–1.0 µm), as an approximation, we assume that FeRR has already been solubilized and is added as Fesmed at emission (Ito and Xu, 2014). In bins 2–4, the FeRR associated with kaolinite and feldspar is also assumed to be solubilized (Ito and Xu, 2014) along with one quarter of the FeRR associated with illite and smectite and is prescribed to Fesmed at emission. The distribution of iron from the silt and clay fractions in the soils follows the brittle fragmentation theory of dust emission (Kok, 2011) and is described in more detail in Scanza et al. (2015) in Table 2a. Total iron in combustion aerosols (Fetcomb) is prescribed and industrial sources of iron are emitted following Luo et al. (2008). Biomass burning is assumed to have a ratio of 0.02 g Fe /g BC in the fine mode and 1.4 g Fe /g BC in the coarse mode (Luo et al., 2008). Fetcomb is partitioned in the model transport bins as follows: fine-mode Fetcomb solely comprises the first bin and the coarse-mode combustion iron is apportioned into bins 2–4 at 20 %, 30 % and 50 %, respectively. Separate tracers account for Fetcomb and Fescomb with 4 % of Fetcomb assumed already soluble and prescribed to Fescomb. 2.4 Iron dissolution Fe dissolution in mineral dust involves complex atmospheric chemical and physical processes. The major processes, i.e., processes that can be observed in laboratory settings, are related to the initial mineral composition, atmospheric temperature and acidity, insolation and concentration of organic acids. In this study, the amount of Fe solubilized in dust is assumed to be irreversible and is calculated using a simplified iron dissolution scheme for three types of iron in dust (FeRR, Femed and Feslow). However, we only explicitly simulate Femed and Feslow both for computational efficiency and because we also assume that a portion of FeRR has already been solubilized at emission. The CAM4 sulfur chemistry model is not configured to interact with dust; thus, the rate of dissolution is only explicitly dependent on temperature and is simplified to a first-order dissolution rate following Eq. (1), RFei=Ki(T ) ×aH+mi×f (∇Gr)×Ai×MWi,(1) d dt[Fesoluble]=RFei×[Feinsoluble],(2) where RFeiis a function of temperature and has units (s−1), i represents either moderately or weakly soluble Fe, Ki(T ) in units of (moles m−2s−1) is the temperature-dependent rate Atmos. Chem. Phys., 18, 14175–14196, 2018 www.atmos-chem-phys.net/18/14175/2018/ R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols 14179 12345678910 10−12 10−11 10−10 10−9 10−8 10−7 10−6 10−5 pH R (1 s ) Iron dissolution reaction rate vs pH (T=298) Illite Hematite Illite [ox=0] Illite [ox=1] Illite [ox=10] Hematite [ox=0] Hematite [ox=1] Hematite [ox=10] −1 Figure 1. Iron dissolution rate for illite and hematite representing moderately soluble iron and weakly soluble iron, respectively. The “+” sign indicates the proton-promoted dissolution rate. Solid, dashed and dotted lines are the addition of the proton-promoted and the oxalatepromoted dissolution rates at [oxalate] =0, [oxalate] =1 and [oxalate] =10 µM, respectively. coefficient (Meskhidze et al., 2005), a(H+)is the proton concentration with an empirical reaction order mi,f (∇Gr) accounts for the change in the dissolution rate with variation from equilibrium (and equals 1 for simplicity; Luo et al., 2008), Aiis the specific surface area of mineral i in units of (m2g−1)and MWiis the molecular weight in units of (g mol−1) for minerali. The pH dependence is parameterized from the concentration of calcite and sulfate at each atmospheric grid box and at each time step. When [calcite]i,j,k >[sulfate]i,j,k, pH is set to 7.5; when model sulfate concentration is greater, pH is set to 2 (Journet et al., 2008; Luo et al., 2008). Femed is comprised of the Fe in illite and smectite; because the mineral abundance of illite is twice that of smectite, we use Killite(T ) and MWillite for the Femed dissolution rate as an additional simplification. The dissolution of Fe in hematite proceeds in three stages (Meskhidze et al., 2005) and we use the fastest of the three stages for Khem(T ) following Ito and Xu (2014) (Fig. 1). The initial concentration of Fescomb is assumed to be 4 % of Fetcomb (Chuang et al., 2005; Luo et al., 2008) in each size bin and we use the Femed dissolution rate for the remaining combustion iron. The impact of including oxalate in an iron dissolution scheme can increase the soluble Fe fraction considerably (Paris et al., 2011; Myriokefalitakis et al., 2011, 2015; Johnson and Meskhidze, 2013; Ito, 2015). Because the formation of atmospheric oxalate is complex, we develop a simple scheme to estimate oxalate concentrations in the model at grid cells within clouds. The concentration of oxalate in µmols L−1is calculated as hC2O2− 4ii,j,k =15 ×OCi,j,k+SOAi,j,k max[OC +SOA].(3) This simple approximation is liable to add a source of uncertainty into our calculations; however, the spatial distribution of oxalate at the surface (Fig. 2a) is comparable to Fig. 2a, c in Myriokefalitakis et al. (2011), which calculates oxalate concentrations using a full-complexity chemical mechanism. In choosing our reference case (REF), multiple simulations were conducted to best match observations of fractional iron solubility and surface distributions of oxalate (Myriokefalitakis et al., 2011); the factor of 15 in Eq. (1) yielded the best results. Previous studies use sulfate as a proxy for calculating model oxalate concentrations (Yu et al., 2005; Johnson and Meskhidze, 2013), and we conduct a sensitivity study to determine the relative importance of this assumption. Because the sources of OC and secondary organic aerosols (SOA) are different than the sources of sulfate, the distribution is different between these two “proxies” for oxalate; the OC +SOA in our model appears to be a better proxy. There tends to be more oxalate in the model simulations in tropical regions (Fig. 2 from Myriokefalitakis et al., 2011), which is better captured in our model simulations using the OC +SOA versus the sulfate proxy for modeled oxalate concentrations. www.atmos-chem-phys.net/18/14175/2018/ Atmos. Chem. Phys., 18, 14175–14196, 2018 14180 R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols − 90 − 45 0 45 90 0.0001 0.0005 0.001 0.005 0.01 0.05 0.1 0.5 1 5 10 −180 −120 −60 0 60 120 180 − 90 − 45 0 45 90 (a) Sfc. oxalate concentration: OC proxy (μmol L ) (b) Sfc. oxalate concentration: SO4 proxy (μmol L ) μM -1 -1 Figure 2. Spatial distribution of annually averaged surface (Sfc.) oxalate concentrations (µM). Oxalate concentrations are derived from modeled organic carbon concentrations (a) and from modeled sulfate concentrations (b). The oxalate-dependent reaction rate is then added to the first-order dissolution rate (Eq. 1) using the model-derived oxalate concentration from Eq. (3): RFei=Ki(T)×aH+mi×f(∇Gr)×Ai×MWi +Ki,oxalate,(4) where Ki,oxalate is Ki,oxalate =ai×hC2O2− 4i+bi.(5) Coefficients aifor illite and hematite are determined from the slope of the best fit of the data in Table 4 in Paris et al. (2011) and bifor illite and hematite correspond to the intercept of the best-fit equations. For illite, a=2.3× 10−7µM−1s−1and b=4.8×10−7s−1and, for hematite, a=9.5×10−9µM−1s−1and b=3.0×10−8s−1. The iron dissolution rate at [C2O2− 4]=0 utilizes Paris et al. (2011) with our method for best fit, which was a result of multiple different simulations where the intercept was examined to find the best match to observations of %Fes. The nonzero intercept for Ki,oxalate accounts for (non-oxalate) in-cloud organic ligand processing. The reference case is simulated under present-day climate conditions and includes the enhancement of oxalate on the rate of soluble Fe formation in cloudy grid boxes only while Eq. (1) is applied at every atmospheric model grid box. Separate tracers are included for dust and for total and soluble iron in dust for Femed and Feslow as well as tracers for total and soluble combustion iron. Calcite is advected to estimate the buffering effect of dust on the pH dependence of the dissolution scheme at each atmospheric model grid box (8 new species ×4 bins +9=41 total tracers added here compared to the original model which has 5 species carried as 13 total tracers to simulate aerosols for climate interactions). 2.5 Sensitivity studies To investigate the impact on soluble Fe formation of various parameters, we include six sensitivity studies (Table 1a). In the first experiment (SS1), the oxalate mechanism is removed and the formation of Fes depends solely on the protonpromoted dissolution rate; the basis for this was to determine if the additional complexity of the oxalate scheme was required to match observations of labile iron. The second sensitivity study (SS2) is identical to the reference case but with model oxalate concentrations estimated from sulfate, was performed because other studies of Fe dissolution use sulfate as a proxy for estimating model oxalate concentrations and follows Eq. (4): hC2O2− 4ii,j,k =15 ×hSO2− 4ii,j,k maxhSO2− 4i,(6) in units of µM. The third sensitivity study (SS3) investigates how important our dissolution mechanism is by comparing to a very simple dissolution scheme from Hand et al. (2004). In this scheme, soluble iron dissolution depends upon cloud presence as parameterized by Eq. (2) in Hand et al. (2004). Total iron is assumed to be 3.5 % of dust and Fes is calculated from kappa ·(Fet −Fes) at each cloudy grid box, where kappa is defined in Eq. (2) in Hand et al. (2004) as P i,j Ci,j Cavg 1 τcld,sol :Ci,j is the cloud fraction at each atmosphere grid box, Cavg is the average cloud fraction at 10◦N and τcld,sol is the soluble decay lifetime (Siefert et al., 1998; Saydam and Senyuva, 2002; Hand et al., 2004). Hand et al. (2004) did not take into account iron from combustion aerosols and there is a wide range of reported dust Fe solubilities at emission 0.1 to 1+% (Luo et al., 2008; Jickells et al., 2005; Shi et al., 2012; Johnson and Meskhidze, 2013; Ito and Xu, 2014). To account for the typically higher solubilities associated with combustion iron, and using the ratio Fes / Fet for dust iron (Table S2d), we assume an initial solubility for combustion and dust Fe of 4 % and 0.89 %, respectively; this corresponds to an enhancement of 4.5 times and hence we chose, somewhat arbitrarily, to increase the dissolution rate by a factor of 5. The total number of tracers for SS3 is 20 (Fetdust, Fetcomb, Fesdust, Fescomb and dust in each of the four size bins). In the fourth sensitivity study (SS4), the spatial dependence of iron on mineralogy is removed and we apply the global average fraction from the reference case for calcite, Fetmed, Fetslow and Fesmed at emission. The fifth and sixth simulations (SS5 and SS6) are identical to REF and Atmos. Chem. Phys., 18, 14175–14196, 2018 www.atmos-chem-phys.net/18/14175/2018/ R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols 14181 Table 1. (a) List of current climate simulation used in this study. [H+] indicates proton-promoted aerosol water iron processing. C2O2− 4OC-proxy is the parameterization for in-cloud iron processing including the oxalate parameterization derived from model concentrations of organic carbon and secondary organic carbon. C2O2− 4SO4-proxy calculates oxalate concentration from model sulfate concentrations. SS4 tracers have no spatial dependence on mineralogy and have prescribed global average fractions from the reference case (REF) at emission. (b) List of preindustrial simulations used in this study. Preindustrial simulations assume pH =2 when model concentrations of sulfate are greater than calcite. (a) [H+] C2O2− 4OC-proxy C2O2− 4SO4-proxy pH REF × × 2 SS1 ×2 SS2 × × 2 SS3 n/a n/a n/a n/a SS4 (no mineralogy) × × 2 SS5 × × 1 SS6 × × 1 (b) [H+] C2O2− 4OC-proxy C2O2− 4SO4-proxy Emissions Dust source PI combustion PI1 × × × PI2 × × × × PI3 × × × × × PI4 × × × × × n/a: not applicable. SS2, respectively; however, pH is set to 1 instead of 2 when sulfate concentrations are greater than calcite concentrations. Because recent studies have identified highly acidic aerosol solutions (Weber et al., 2016; Guo et al., 2016), we wanted to test if the increase in the iron dissolution would better match observations of labile iron solubility. 2.6 Preindustrial simulations Four preindustrial simulations are performed using the iron dissolution mechanism described in Sect. 2.4. Here we refer to “preindustrial chemistry” as historical emissions for CMIP5 (Lamarque et al., 2010), which includes reduced emissions for sulfate, OC and BC. The first simulation (PI1) has preindustrial chemical emissions (1850) but current climate dust and combustion sources. The second simulation (PI2) includes preindustrial chemistry and dust sources with current combustion. PI3 follows PI2 with preindustrial combustion sources and PI4 is identical to PI3 with the sulfate proxy used for calculating model oxalate concentrations instead of the OC +SOA proxy (Table 1). 2.7 Comparison to observations Observations of total and soluble iron are complicated due to various collection methods (e.g., using different collection substrates and sampling with or without particle size segregation), different filter sizes used to define the soluble iron fraction, different solvents used for extraction (which determines the species of iron defined as soluble) and myriad definitions of the form of iron which is in fact soluble. The majority of observations are collected during cruise campaigns (daily means) as surface concentrations or dry deposition of particulate matter. These sampling campaigns inherently fail to capture the ephemeral nature of aerosol transport over ocean basins, varying in both space and time. Aerosol collections are later leached in one or more solvents to extract the soluble iron fraction from the insoluble component. A handful of alternative collection methods have been used, including the sampling of aerosols at 70 m above ground level (Winton et al., 2015) over periods of 1 to 7 weeks, the ACE-Asia campaign that sampled total suspended particles at the surface (Chuang et al., 2005) during dust storms and the collection of rainwater (wet deposition) at 2 m above ground in the remote Southern Ocean (Heimburger et al., 2013). Solvents used for extraction range from acidic solutions and ultra-pure water to alkaline seawater, with obvious differences in soluble iron yield due to the differences in acidity. The definition of soluble iron is rather ambiguous in part due to the field of study (i.e., oceanography, atmospheric science), the leaching solution and the oxidation state of dissolved iron. In this study, we attempt to compare our model only to observations of labile iron, which we define as dissolved iron in either the ferrous or ferric oxidation state. Fine-mode and coarse-mode total and labile iron in the model is compared to available published observations from Baker et al. (2006a, b) and Chen and Siefert (2004) along with previously unpublished fineand coarse-mode total and labile iron observations from Clifton Buck (personal communication, 2013, 2017). Bulk total and labile iron observawww.atmos-chem-phys.net/18/14175/2018/ Atmos. Chem. Phys., 18, 14175–14196, 2018 14182 R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols Table 2. Root mean square error (RMSE) for the reference case (REF) and for sensitivity studies 1–6. REF: proton-promoted and oxalatepromoted iron dissolution; SS1: proton-promoted iron dissolution; SS2: REF with oxalate dependent on sulfate; SS3: simple cloud parameterization for iron dissolution (Hand et al., 2004); SS4: REF with no spatial dependence on mineralogy; SS5: REF, pH =1; SS6: SS2, pH =1. RMSE REF SS1 SS2 SS3 SS4 SS5 SS6 Fesfine 5.9 6.1 5.9 6.6 6.2 5.9 5.9 Fescoarse 26.5 13.8 27.9 5.3 20.7 27.2 28.6 Festotal 25.5 17.1 26.72 9.2 19.4 26.1 27.3 Fetfine 190.6 190.6 190.6 193.2 194.1 190.6 190.6 Fetcoarse 1448.4 1448.6 1448.0 1356.6 1241.8 1448.4 1447.9 Fettotal 2432.1 2432.2 2432.0 1149.2 2362.6 2432.1 2432.0 %Fesfine 19.5 20.3 19.4 20.6 19.9 19.6 19.4 %Fescoarse 7.1 7.6 6.9 7.5 7.1 7.4 7.2 %Festotal 9.7 10.7 9.5 10.5 9.9 10.1 9.9 tions that do not distinguish particle sizes are from supplementary Table 2 provided in Mahowald et al. (2009), from Sholkovitz et al. (2012) and from Baker et al. (2016). A set of metrics is developed to compare model simulations with the available observations to determine which set of parameterizations is necessary. The root mean square error (RMSE) is calculated for model versus observations as sn P i=1 (modi−obsi)2 n for total iron, labile iron and solubility (fine, coarse and fine +coarse). Similarly, the Pearson correlation coefficient is calculated between the model and observations. As a main metric in evaluating model performance, we develop a semiquantitative metric defined as the mean of model values at grid boxes with observations divided by the mean of the observations and will hereafter refer to this metric with the acronym MMO. The MMO does not intend to evaluate the model’s ability to capture observational variability but simply to assess if the model can reasonably estimate the observational mean. 3 Results 3.1 Comparison of observations to the reference case In order to assess the ability of the model to simulate total, soluble iron and fractional iron solubility, model results are compared to available observations for fine- (diameter <1 µm), coarse- (1 µm ≥diameter ≤10 µm) and bulk-mode particles (all sizes) as most available observations do not differentiate particle size. The total iron at emission is calculated based on the distribution of minerals containing iron (Journet et al., 2008; Ito and Xu, 2014; Scanza et al., 2015); Fet is approximately 3.2 % of dust which is comparable to the assumption that 3.5 % of dust is comprised of iron (Hand et al., 2004). Simulated surface concentrations of total iron in the fine and coarse mode are underand over-predicted, respectively, compared to observations. Because 98 % of Fet has prescribed diameters >1 µm, Fet is overestimated in the tropical and northern latitudes (near source regions) and underestimated below 15◦S (Fe <1 µm is more likely to undergo long-range transport). Because there are relatively few observations that distinguish the particle size distribution of iron, the significance of comparing total iron surface concentrations is difficult to assess; however, the model is underestimating the dynamic range in the observations (Fig. S1 in the Supplement). This is likely due in part to the fact that we are comparing annually averaged concentrations to observations that are influenced by daily weather fluctuations. Because the dust model used in this study was optimized to best capture observations of particle size distribution, deposition and aerosol optical depth (Figs. 8, 9 in Albani et al., 2014), surface concentrations of dust are generally over-predicted with the exception of remote Southern Hemisphere (SH) ocean regions. In addition, most dust models exhibit a low bias in SH dust deposition (Huneeus et al., 2010, 2011; Albani et al., 2014). Hence, it is not possible to simultaneously match dust deposition, AOD and surface concentrations. Separating total and soluble iron allows us to identify these signals and to evaluate the fractional iron solubility. Fine-mode labile iron (Fes) surface concentrations are within the range of the observations in the tropical and northern latitudes and are underestimated off the coast of Patagonia; a significant limitation arises from the lack of observations in the regions where iron deposition should have the largest impact (see Sect. 3.3). Coarse-mode Fes is considerably over-predicted by the model, which is consistent with the model bias in dust surface concentrations, resulting in an overestimate in the total concentration of Fes (Fig. S2) (individual fineand coarse-mode Fes distributions not shown). Therefore, the fractional iron solubility at the model surface provides a more meaningful comparison to observations and facilitates the evaluation of the iron-processing mechanism. In addition, examining total iron and percent solubilAtmos. Chem. Phys., 18, 14175–14196, 2018 www.atmos-chem-phys.net/18/14175/2018/ R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols 14183 ity separately allows us to evaluate the distribution and the solubilization of iron as distinct processes, thereby simplifying the identification and resolution of model biases in these two processes. However, additional uncertainty is introduced by considering the distributions of Fet and %Fes separately since, inherently, they are coupled. For example, erroneously high dust emissions would result in fractional solubilities that are biased low because of this coupling. When considering percent labile iron, theory and observations support lower solubilities near dust sources where pH is likely higher (due to the presence of calcite in dust) and where lithogenic iron has not yet undergone significant atmospheric processing; in contrast, because combustion sources are typically associated with emissions of acidic species, we expect the solubility of anthropogenic aerosols to be higher (Li et al., 2017). However, because dust iron emissions are an order of magnitude greater than emissions of combustion iron, lower solubilities near coastal regions should be dominant. In particular, we hypothesize that the simulations will predict higher percentages of Fes over remote ocean regions and, indeed, our simulations are consistent with our hypothesis; this is in line with theory and observations that show the gravitational settling of larger, less soluble iron containing aerosols near their source and long-range transport of smaller aerosols enabling multiple cycles of evaporation and condensation (Jickells et al., 2005; Baker and Croot, 2010). Simulated %Fes <1 µm is small off the coast of North Africa and Patagonia and is an order of magnitude higher in the equatorial Pacific and southern latitude ocean regions (Fig. 3b). The relatively few observations indicate reasonable agreement with the model near and directly downwind of North Africa; however, there are few observations to substantiate the predictions of higher solubilities in the HNLC regions (Fig. 3a). One standard deviation from the mean of daily averaged model concentrations indicates that the model is within range of the observations (Fig. 3c). The spatial distribution of coarse-mode %Fes is similar to the fine mode, and again our comparison is limited by the lack of observations (Fig. 4). Simulations of the sum of the fineand coarsemode %Fes show greater spatial heterogeneity in the Southern Hemisphere and less variation in the Northern Hemisphere (Fig. 5). While the scatterplot comparing annually averaged %Fes versus daily cruise-based observations suggests that the model is not capturing the range of the observations, the daily averaged standard deviations for each model value indicate that our mechanism can generally capture the observations. Thus, with MMO equal to 0.86, we conclude that the soluble-iron-processing mechanism employed is a reasonable representation of the atmospheric processing that dust and combustion aerosols undergo during transport. 3.2 Evaluation of sensitivity studies Atmospheric iron dissolution is very complex and required many simplifications and assumptions in order to simulate the processes involved. In this study, we develop a set of metrics to evaluate the strengths and weaknesses of our proposed iron dissolution scheme and to assess its ability to capture the true processing of atmospheric aerosols. For each simulated variable, we compare the results to available observations and calculate the RMSE (Table 2), the MMO (Table 3) and the correlation coefficient (Table 4) for the reference case and for each sensitivity study. When considering RMSE for fractional iron solubility for each case, REF compared to SS2 (oxalate calculated with a sulfate proxy), SS4 (no dependence on mineral spatial distribution at emission) and SS5– 6 (higher acidity) has similar errors. SS1, the case without oxalate processing, and SS3, which corresponds to the very simple dissolution method from Hand et al. (2004), are both associated with larger errors, indicating that the very simple mechanism (SS3) and mechanism solely based on acidic processing (SS1) are not sufficient representations. The correlation coefficients across all cases are weak, especially for the fractional iron solubilities. In order to counter the ambiguity surrounding which mechanism best matches observations, we believe the semiquantitative MMO metric is more meaningful. A ratio greater than 1 indicates that the model is over-predicting iron while a ratio less than 1 signifies an underestimation compared to the mean of the observations. The reference case was chosen from many different simulations with varying parameters after carefully reviewing the spatial distributions of soluble iron as well as the model metrics. Although SS5 and SS6 most closely match the observational mean for total fractional iron solubility, the pH when the concentration of sulfate is greater than calcite (pH =1) is not typically considered realistic in atmospheric waters. Recent studies that identify highly acidic aerosol solutions warrant more investigation as to whether SS5 and SS6 are in fact better representations of atmospheric iron processing. SS1 and SS3 greatly underestimate the observed mean for fractional iron solubility, leading to the conclusion that the oxalate enhancement within an iron dissolution scheme is a vital inclusion. 3.3 Regional and global iron deposition to ocean basins Soluble iron deposition to global ocean basins from dust and combustion aerosols is a significant source of new iron, a limiting nutrient in many regions. This represents an important biogeochemical process that may indirectly impact climate via a cooling effect by promoting the uptake of atmospheric CO2. The global ocean is divided into 12 regions (Table S4) defined by Gregg et al. (2003) in order to quantify deposition to ocean regions that are most likely to be iron limited. Annually averaged global deposition of total and soluble iron from dust to the oceans is 17.4 % and 0.54 Tg yr−1, respectively, and corresponds to the dissolution of 3.1 % of dust iron during atmospheric transport. Including the contribution from combustion iron increases the soluble iron deposition to 3.3 % of total iron. This is within the ranges reported (model: www.atmos-chem-phys.net/18/14175/2018/ Atmos. Chem. Phys., 18, 14175–14196, 2018 14190 R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols −90 −45 0 45 90 −1.5 −1 −0.5 0 0.5 1 1.5 −1.5 −1 −0.5 0 0.5 1 1.5 −90 −45 0 45 90 −1.5 −1 −0.5 0 0.5 1 1.5 −1.5 −1 −0.5 0 0.5 1 1.5 −180 −90 0 90 180 −90 −45 0 45 90 −1.5 −1 −0.5 0 0.5 1 1.5 −180 −90 0 90 180 −1.5 −1 −0.5 0 0.5 1 1.5 (a) REF/SS1 (b) REF/SS2 (c) REF/SS3 (d) REF/SS4 (e) REF/SS5 (f) REF/SS6 2.0+ 1.75 1.5 1.25 1.0 0.75 0.5 0.25 0.0 Figure 8. The spatial distribution of the Fes deposition ratio of the reference case over the six sensitivity studies. SS1 corresponds to the case with proton-promoted iron dissolution and no processing via organic ligand dissolution enhancement. The processing mechanism for SS2 is identical to the reference case; however, it uses a sulfate proxy for calculating model oxalate concentrations instead of the OC proxy used in the reference case. SS3 corresponds to a simplified dissolution mechanism defined in Hand et al. (2004); a decay rate is applied when iron comes into contact with a cloud. The SS4 simulation has the same processing mechanism as the reference case; however, the spatial dependence of total and soluble iron on mineralogy at emission is replaced with the average values for these variables from the reference case. SS5 and SS6 are identical to the reference case and SS2, respectively, with the exception of prescribing the pH at 1 when sulfate concentrations exceed calcite concentrations. White indicates no change, red indicates that the deposition in the reference case is higher than the sensitivity study and blue indicates higher deposition compared to the reference case. A comparison of our results with the results obtained from previous modeling studies enables some determination of the robustness of the iron dissolution scheme used in this study. SS1 is compared to a similar study (Ito and Xu, 2014) since this study did not include an oxalate scheme. The deposition of total and soluble iron from dust and combustion sources is calculated for ocean regions in the western and eastern North Pacific in Table S2e and compared to Table 6 in Ito and Xu (2014). Because we used similar definitions for Fet and a similar dissolution scheme for acidic processing of iron, our total and soluble iron reasonably matches their data. In both the western and eastern North Pacific basins, Ito and Xu (2014) report dust iron solubility of around 2 %; we calculate dust iron solubility of around 1.4 % and 1.6 % for the two basins, respectively. This disparity may be a result of not including separate tracers and dissolution rates for FeRR along with some simplifications we made in the dissolution rate (i.e., not accounting for the deviation from equilibrium for Fet →Fes). The sources of combustion in this study are different and we did not include a combustion-specific processing scheme. For lack of a better understanding of combustion iron dissolution, combustion iron in this study is mobilized with the rate used for Fesmed. Johnson and Meskhidze (2013) include an oxalateprocessing scheme based on oxalate-promoted dissolution data reported in Paris et al. (2011). It is difficult to compare our results since we interpreted the Paris et al. (2011) data differently to calculate our coefficients for the oxalate dissolution rates for illite and hematite. We predict higher dissolved iron deposition from dust to the global ocean basins compared to their prediction (0.54 versus 0.21 Tg yr−1) and for the range in Fes deposition reported here (0.21– 0.41 Tg yr−1) from Luo et al. (2008), Luo and Gao (2010), and Okin et al. (2011). Johnson and Meskhidze (2013) report a 75 % increase in Fes deposition when the oxalate mechanism is included, which is similar to our prediction of a 63 % increase (Table 6b). A recent study that included both dust and combustion iron along with acidic processing, photochemical reduction and organic ligand processing (Myriokefalitakis et al., 2015) reports Fet emissions from mineral dust and combustion at ∼35 and 2 Tg yr−1, respectively. In fact, the method by which this study determined the iron content at emission is nearly identical to the method used in our study; we calculate ∼57 and 2 Tg yr−1of total iron, which is very similar to their emissions if we account for the 38 % difference in Atmos. Chem. Phys., 18, 14175–14196, 2018 www.atmos-chem-phys.net/18/14175/2018/ R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols 14191 −90 −45 0 45 90 −1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1 −1 80 − 90 090 1 80 −90 −45 0 45 90 −1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1 −90 −45 0 45 90 5e−05 0.0005 0.005 0.05 0.5 5 50 (a) PI3 Fes deposition (b) REF/PI3 (c) (REF/SS2)/(PI3/PI4) 2.0 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 5e1 5e0 5e-1 5e-2 5e-3 5e-4 5e-5 Tg yr-1 Figure 9. Spatial distribution of soluble iron deposition for the preindustrial simulation with the year 1850 chemical emissions, dust source strengths and combustion emissions (PI3) (Tg yr−1) in panel (a). Spatial distribution of the Fes deposition ratio of the reference case over PI3 (b), and the difference between the sulfate and OC proxies for current over preindustrial (c). dust emission between the two studies. However, the global annual flux of Fet from dust and combustion in our study is nearly 3 times higher than the 0.496 Tg yr−1reported here. A significant difference between our mechanism involves the oxalate processing; we assume that moderately soluble iron (iron from illite and smectite) undergoes organic ligand processing with a rate determined from Paris et al. (2011) for illite and weakly soluble iron (iron from kaolinite, feldspar and hematite) is processed with a rate derived from the hematite data. This study only processes illite and hematite and neglects the iron from smectite, feldspar and kaolinite, the first of which has significant iron concentrations as reported in Journet et al. (2008). In addition, because we used a simple parameterization for pH, we do not simulate the decrease in the oxalate mechanism at low pH. Regardless, our ability to match observations, despite the simplified parameterizations we employed compared with the full chemistry scheme used here, is at least as good if not better than the metric values reported in the Supplement. Ito (2015) includes an oxalate-processing mechanism for the dissolution of iron in combustion aerosols. The global annual deposition for labile iron reported is 1.07 Tg yr−1, which is closer to the 1.4 Tg yr−1estimated in this study. Theory and observations of iron processing show an inverse trend in total and soluble iron concentrations with fractional iron solubility (Sholkovitz et al., 2012). The results from our work suggest that this relationship cannot be uniquely constrained by considering either dust or combustion sources of iron alone. Figure 10 shows the total fractional iron solubility versus total iron for dust iron and combustion iron. In order to illustrate the difficulty in determining which process is responsible for the inverse trend, for each observation, we define a box spanning the standard deviation of modeled fractional iron solubility and of modeled total iron. When the observation lies inside or outside of this box, we can calculate the percentage of model points that are within the range of the observations. For dust processing only, for combustion only and for the total (dust +combustion), 75.4 %, 73.8 % and 73.3 % of model values are within the observational range, respectively. If only combustion sources were included and assumed to be more soluble, the modeled values for soluble iron would all be shifted upwards in Fig. 10 and could match available observations quite well. 6 Conclusion A medium-complexity soluble-iron-processing mechanism is developed for use in an Earth system model and includes a proton-promoted aerosol-water-processing scheme along with an in-cloud organic-ligand-promoted dissolution scheme. We define a semi-quantitative metric as the model mean at points with observations divided by the observational mean (MMO); fractional iron solubility MMO is 0.86, indicating that, while the model is not capturing all of the observational variability, it is in reasonable agreement with the observational mean. Several sensitivity studies are performed to ascertain the degree of complexity needed to match observations. A simple first-order decay rate from Hand et al. (2004) is too slow compared to observations than our intermediate-complexity mechanism. Our results indicate that the acidic processing alone is insufficient in matching observations and that the additional ligand-promoted mechanism was required; the mechanism employed is reasonably able to capture some but not all of the observational features of soluble iron and iron solubility, similar to the more complicated mechanisms employed in previous studies (Meskhidze et al., 2005; Johnson and Meskhidze, 2013; Ito and Xu, 2014; Myriokefalitakis et al., 2015; Ito, 2015). The in-cloud organic-ligand-processing mechanism is likely www.atmos-chem-phys.net/18/14175/2018/ Atmos. Chem. Phys., 18, 14175–14196, 2018 14192 R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols 10−1 100101102103104 10−1 100 101 102 10−1 100101102103104 10−1 100 101 102 10−1 100101102103104 10−3 10−2 10−1 100 101 102 FETatOBS (ng/m3) Fractional iron solubility (%) 220 200 2000 0.05 0.5 5 50 %Fes at OBS %Fes at OBS−dust %Fes at OBS−comb %Fes(OBS) ∗ %Fes(OBS) O O %Fes at OBS ∗ %Fes at OBS (outside 1 stdd) -- -- 1 SD (daily avg.) ∗ %Fes(OBS) O O %Fes at OBS-dust ∗ %Fes at OBS (outside 1 stdd) -- -- 1 SD (daily avg.) ∗ %Fes(OBS) O O %Fes at OBS-comb ∗ %Fes at OBS (outside 1 SD) -- -- 1 SD (daily avg.) Fractional iron solubility (%) Fractional iron solubility (%) FetAtOBS (ng m ) FetAtOBS (ng m ) (a) (b) (d) (c) -3 -3 Figure 10. Scatterplots of total iron vs. fractional iron solubility. For panels (a)–(c), black points correspond to observations, blue circles correspond to annually averaged model values at grid boxes where we have observational data. Red stars indicate model values that fall outside the range (predicted from the average of 1 standard deviation from daily averaged values for both total iron and fractional iron solubility) in relation to observed values. Panel (a) corresponds to fractional iron solubility from both dust and combustion, panel (b) is for fractional iron solubility from dust only and panel (c) is for combustion iron only. Panel (d) corresponds to the zoomed-in version of Fig. 1 in Sholkovitz et al. (2012) shown logarithmically. FetAtOBS: the model value for Fet at the geographic location where Fet was observed. to be responsible for the majority of the atmospheric processing (>60 %). An additional objective was characterizing the relative importance of the contribution of iron from combustion aerosols versus iron from dust to different ocean basins. Upon deposition to global oceans, global mean combustion iron solubility to total fractional iron solubility is 8.2 %; however, the contribution of fractional iron solubility from combustion sources to ocean basins below 15◦S is approximately 50 %. We estimate approximately a doubling of soluble iron deposition to the oceans in the current climate relative to the preindustrial, due to increased dust, combustion iron and acidifying compounds, and a ∼1.7-fold increase in the iron-limited regions of the ocean (Tables 5–7). Even in the preindustrial, when industrial sources of iron are much lower, combustion sources of iron are likely to be important in some regions. Building the framework to model iron dissolution in a full-complexity Earth system model was the main goal of this study; however, more observations are needed in remote ocean regions to evaluate model performance. Additionally, more laboratory and field experiments are necessary to better Atmos. Chem. Phys., 18, 14175–14196, 2018 www.atmos-chem-phys.net/18/14175/2018/ R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols 14193 constrain the importance of combustion iron and atmospheric processing of mineral aerosols with both inorganic and organic acids. The sensitivity studies indicate that including an organic-ligand-processing mechanism is the most important parameter in order to match observations. Further investigation into atmospheric aerosol solution acidity is necessary in order to distinguish whether the sensitivity studies with higher acidity are better representations of soluble iron processing in the atmosphere. Data availability. The observational data and the model output for the base case and the PI simulations are available at http://www.geo.cornell.edu/eas/PeoplePlaces/Faculty/mahowald/ dust/scanzaetal2018 (Mahowald, 2018). The Supplement related to this article is available online at https://doi.org/10.5194/acp-18-14175-2018supplement. Author contributions. RAS and NMM designed the experiments; RAS performed model simulations; and RAS, NMM, and DSH analyzed results. AB and CB provided observations, and all authors revised the paper. Competing interests. The authors declare that they have no conflict of interest. Special issue statement. This article is part of the special issue “Anthropogenic dust and its climate impact”. It is not associated with a conference. Acknowledgements. We would like to acknowledge the support of DOE DE-SC0006735 and NSF 1049033. Carlos Pérez GarcíaPando acknowledges long-term support from the AXA Research Fund through the AXA Chair on Sand and Dust Storms, as well as the support received through the Ramón y Cajal program (grant RYC-2015-18690) of the Spanish Ministry of Economy and Competitiveness. Edited by: Jianping Huang Reviewed by: two anonymous referees References Albani, S., Mahowald, N., Perry, A., Scanza, R., Zender, C., Heavens, N., Maggi, V., Kok, J., and Otto-Bliesner, B.: Improved dust representation in the Community Atmosphere Model, J. Adv. Model. Earth Sy., 6, 541–570, 2014. Baker, A. and Croot, P.: Atmospheric and marine controls on aerosol iron solubility in seawater, Mar. Chem., 120, 4–13, 2010. Baker, A. and Jickells, T.: Mineral particle size as a control on aerosol iron solubility, Geophys. Res. Lett., 33, L17608, https://doi.org/10.1029/2006GL026557, 2006. Baker, A., French, M., and Linge, K.: Trends in aerosol nutrient solubility along a west–east transect of the Saharan dust plume, Geophys. Res. Lett., 33, L07805,https://doi.org/10.1029/2005GL024764, 2006a. Baker, A., Jickells, T., Witt, M., and Linge, K.: Trends in the solubility of iron, aluminium, manganese and phosphorus in aerosol collected over the Atlantic Ocean, Mar. Chem., 98, 43–58, 2006b. Baker, A. R., Laskina, O., and Grassian, V. H.: Processing and Ageing in the Atmosphere, in: Mineral Dust, Springer, 75–92, 2014. Baker, A. R., Thomas, M., Bange, H. W., and Plasencia Sánchez, E.: Soluble trace metals in aerosols over the tropical southeast Pacific offshore of Peru, Biogeosciences, 13, 817–825, https://doi.org/10.5194/bg-13-817-2016, 2016. Barbeau, K., Moffett, J. W., Caron, D. A., Croot, P. L., and Erdner, D. L.: Role of protozoan grazing in relieving iron limitation of phytoplankton, 1996. Batjes, N.: A world dataset of derived soil properties by FAO– UNESCO soil unit for global modelling, Soil Use Manage., 13, 9–16, 1997. Böke, H., Göktürk, E. H., Caner-Saltık, E. N., and Demirci, ¸S.: Effect of airborne particle on SO 2–calcite reaction, Appl. Surf. Sci., 140, 70–82, 1999. Boyd, P. W., Watson, A. J., Law, C. S., Abraham, E. R., Trull, T., Murdoch, R., Bakker, D. C., Bowie, A. R., Buesseler, K., and Chang, H.: A mesoscale phytoplankton bloom in the polar Southern Ocean stimulated by iron fertilization, Nature, 407, 695–702, 2000. Broecker, W. S., Takahashi, T., Simpson, H., and Peng, T.-H.: Fate of fossil fuel carbon dioxide and the global carbon budget, Science, 206, 409–418, 1979. Chen, Y. and Siefert, R. L.: Seasonal and spatial distributions and dry deposition fluxes of atmospheric total and labile iron over the tropical and subtropical North Atlantic Ocean, J. Geophys. Res.- Atmos., 109, D09305, https://doi.org/0.1029/2003JD003958, 2004. Chuang, P., Duvall, R., Shafer, M., and Schauer, J.: The origin of water soluble particulate iron in the Asian atmospheric outflow, Geophys. Res. Lett., 32, L07813, https://doi.org/10.1029/2004GL021946, 2005. Claquin, T., Schulz, M., and Balkanski, Y.: Modeling the mineralogy of atmospheric dust sources, J. Geophys. Res.-Atmos., 104, 22243–22256, 1999. Cornell, R. and Schindler, P.: Photochemical dissolution of goethite in acid/oxalate solution, Clay. Clay Miner., 35, 347–352, 1987. Crutzen, P. J. and Zimmermann, P. H.: The changing photochemistry of the troposphere, Tellus B, 43, 136–151, 1991. d’Almeida, G. A.: On the variability of desert aerosol radiative characteristics, J. Geophys. Res.-Atmos., 92, 3017–3026, 1987. De Baar, H. J., De Jong, J. T., Bakker, D. C., Löscher, B. M., Veth, C., Bathmann, U., and Smetacek, V.: Importance of iron for plankton blooms and carbon dioxide drawdown in the Southern Ocean, Nature, 373, 412–415, 1995. Desboeufs, K., Losno, R., and Colin, J.-L.: Factors influencing aerosol solubility during cloud processes, Atmos. Environ., 35, 3529–3537, 2001. www.atmos-chem-phys.net/18/14175/2018/ Atmos. Chem. Phys., 18, 14175–14196, 2018 14194 R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols Duce, R. A. and Tindale, N. W.: Atmospheric transport of iron and its deposition in the ocean, Limnol. Oceanogr., 36, 1715–1726, 1991. Fung, I. Y., Meyn, S. K., Tegen, I., Doney, S. C., John, J. G., and Bishop, J. K.: Iron supply and demand in the upper ocean, Global Biogeochem. Cy., 14, 281–295, 2000. Gent, P. R., Danabasoglu, G., Donner, L. J., Holland, M. M., Hunke, E. C., Jayne, S. R., Lawrence, D. M., Neale, R. B., Rasch, P. J., and Vertenstein, M.: The community climate system model version 4, J. Climate, 24, 4973–4991, 2011. Gregg, W. W., Conkright, M. E., Ginoux, P., O’Reilly, J. E., and Casey, N. W.: Ocean primary production and climate: Global decadal changes, Geophys. Res. Lett., 30, 1809, https://doi.org/10.1029/2003GL016889, 2003. Guieu, C., Bonnet, S., Wagener, T., and Loÿe-Pilot, M. D.: Biomass burning as a source of dissolved iron to the open ocean?, Geophys. Res. Lett., 32, L19608, https://doi.org/10.1029/2005GL022962, 2005. Guo, H., Sullivan, A. P., Campuzano-Jost, P., Schroder, J. C., Lopez-Hilfiker, F. D., Dibb, J. E., Jimenez, J. L., Thornton, J. A., Brown, S. S., and Nenes, A.: Fine particle pH and the partitioning of nitric acid during winter in the northeastern United States, J. Geophys. Res.-Atmos., 121, 10355–10376, https://doi.org/10.1002/2016JD025311, 2016. Hand, J., Mahowald, N., Chen, Y., Siefert, R., Luo, C., Subramaniam, A., and Fung, I.: Estimates of atmospheric-processed soluble iron from observations and a global mineral aerosol model: Biogeochemical implications, J. Geophys. Res.-Atmos., 109, D17205, https://doi.org/10.1029/2004JD004574, 2004. Heimburger, A., Losno, R., and Triquet, S.: Solubility of iron and other trace elements in rainwater collected on the Kerguelen Islands (South Indian Ocean), Biogeosciences, 10, 6617–6628, https://doi.org/10.5194/bg-10-6617-2013, 2013. Huneeus, N., Schulz, M., Balkanski, Y., Griesfeller, J., Kinne, S., Prospero, J., Bauer, S., Boucher, O., Chin, M., and Dentener, F.: Global dust model intercomparison in AeroCom phase I, Atmos. Chem. Physics Discuss., 10, 23781–23864, https://doi.org/10.5194/acpd-10-23781-2010, 2010. Huneeus, N., Schulz, M., Balkanski, Y., Griesfeller, J., Prospero, J., Kinne, S., Bauer, S., Boucher, O., Chin, M., Dentener, F., Diehl, T., Easter, R., Fillmore, D., Ghan, S., Ginoux, P., Grini, A., Horowitz, L., Koch, D., Krol, M. C., Landing, W., Liu, X., Mahowald, N., Miller, R., Morcrette, J.-J., Myhre, G., Penner, J., Perlwitz, J., Stier, P., Takemura, T., and Zender, C. S.: Global dust model intercomparison in AeroCom phase I, Atmos. Chem. Phys., 11, 7781–7816, https://doi.org/10.5194/acp11-7781-2011, 2011. Hurrell, J. W., Holland, M. M., Gent, P. R., Ghan, S., Kay, J. E., Kushner, P., Lamarque, J.-F., Large, W. G., Lawrence, D., and Lindsay, K.: The community earth system model: a framework for collaborative research, B. Am. Meteorol. Soc., 94, 1339– 1360, 2013. Ito, A.: Atmospheric processing of combustion aerosols as a source of bioavailable iron, Environ. Sci. Tech. Let., 2, 70–75, 2015. Ito, A. and Shi, Z.: Delivery of anthropogenic bioavailable iron from mineral dust and combustion aerosols to the ocean, Atmos. Chem. Phys., 16, 85–99, https://doi.org/10.5194/acp-16-852016, 2016. Ito, A. and Xu, L.: Response of acid mobilization of ironcontaining mineral dust to improvement of air quality projected in the future, Atmos. Chem. Phys., 14, 3441–3459, https://doi.org/10.5194/acp-14-3441-2014, 2014. Jickells, T., An, Z., Andersen, K. K., Baker, A., Bergametti, G., Brooks, N., Cao, J., Boyd, P., Duce, R., and Hunter, K.: Global iron connections between desert dust, ocean biogeochemistry, and climate, Science, 308, 67–71, 2005. Jickells, T. D. and Spokes, L. J.: Atmospheric iron inputs to the oceans, IUPAC series on analytical and physical chemistry of environmental systems, 7, 85–122, 2001. Johnson, K.: Iron supply and demand in the upper ocean: Is extraterrestrial dust a significant source of bioavailable iron?, Global Biogeochem. Cy., 15, 61–63, 2001. Johnson, M. S. and Meskhidze, N.: Atmospheric dissolved iron deposition to the global oceans: effects of oxalate-promoted Fe dissolution, photochemical redox cycling, and dust mineralogy, Geosci. Model Dev., 6, 1137–1155, https://doi.org/10.5194/gmd6-1137-2013, 2013. Journet, E., Desboeufs, K. V., Caquineau, S., and Colin, J. L.: Mineralogy as a critical factor of dust iron solubility, Geophys. Res. Lett., 35, L07805, https://doi.org/10.1029/2007GL031589, 2008. Kok, J. F.: A scaling theory for the size distribution of emitted dust aerosols suggests climate models underestimate the size of the global dust cycle, P. Natl. Acad. Sci. USA, 108, 1016–1021, 2011. Lafon, S., Rajot, J.-L., Alfaro, S. C., and Gaudichet, A.: Quantification of iron oxides in desert aerosol, Atmos. Environ., 38, 1211– 1218, 2004. Lamarque, J.-F., Bond, T. C., Eyring, V., Granier, C., Heil, A., Klimont, Z., Lee, D., Liousse, C., Mieville, A., Owen, B., Schultz, M. G., Shindell, D., Smith, S. J., Stehfest, E., Van Aardenne, J., Cooper, O. R., Kainuma, M., Mahowald, N., McConnell, J. R., Naik, V., Riahi, K., and van Vuuren, D. P.: Historical (1850–2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols: methodology and application, Atmos. Chem. Phys., 10, 7017–7039, https://doi.org/10.5194/acp-10-7017-2010, 2010. Lamarque, J.-F., Emmons, L. K., Hess, P. G., Kinnison, D. E., Tilmes, S., Vitt, F., Heald, C. L., Holland, E. A., Lauritzen, P. H., Neu, J., Orlando, J. J., Rasch, P. J., and Tyndall, G. K.: CAM-chem: description and evaluation of interactive atmospheric chemistry in the Community Earth System Model, Geosci. Model Dev., 5, 369–411, https://doi.org/10.5194/gmd-5369-2012, 2012. Lancaster, N. and Baas, A.: Influence of vegetation cover on sand transport by wind: field studies at Owens Lake, California, Earth Surf. Proc. Land., 23, 69–82, 1998. Lelieveld, J., Crutzen, P. J., and Dentener, F. J.: Changing concentration, lifetime and climate forcing of atmospheric methane, Tellus B, 50, 128–150, 1998. Li, W., Xu, L., Liu, X., Zhang, J., Lin, Y., Yao, X., Gao, H., Zhang, D., Chen, J., and Wang, W.: Air pollution–aerosol interactions produce more bioavailable iron for ocean ecosystems, Sci. Adv., 3, e1601749, https://doi.org/10.1126/sciadv.1601749, 2017. Luo, C. and Gao, Y.: Aeolian iron mobilisation by dust–acid interactions and their implications for soluble iron deposition to the ocean: a test involving potential anthropogenic organic acidic species, Environ. Chem., 7, 153–161, 2010. Atmos. Chem. Phys., 18, 14175–14196, 2018 www.atmos-chem-phys.net/18/14175/2018/ R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols 14195 Luo, C., Mahowald, N., Bond, T., Chuang, P., Artaxo, P., Siefert, R., Chen, Y., and Schauer, J.: Combustion iron distribution and deposition, Global Biogeochem. Cy., 22, GB1012, https://doi.org/10.1029/2007GB002964, 2008. Mahowald, N. M.: Observation data and model ouput, available at: http://www.geo.cornell.edu/eas/PeoplePlaces/Faculty/ mahowald/dust/scanzaetal2018, last access: 3 October 2018. Mahowald, N., Hamilton, D., Mackey, K., Moore, J., Baker, A., Scanza, R., and Zhang, Y.: Aerosol trace metal leaching and impacts on marine microorganisms, Nat. Commun., 9, 2614, https://doi.org/10.1038/s41467-018-04970-7, 2018. Mahowald, N. M., Baker, A. R., Bergametti, G., Brooks, N., Duce, R. A., Jickells, T. D., Kubilay, N., Prospero, J. M., and Tegen, I.: Atmospheric global dust cycle and iron inputs to the ocean, Global Biogeochem. Cy., 19, GB4025, https://doi.org/10.1029/2004GB002402, 2005. Mahowald, N. M., Muhs, D. R., Levis, S., Rasch, P. J., Yoshioka, M., Zender, C. S., and Luo, C.: Change in atmospheric mineral aerosols in response to climate: Last glacial period, preindustrial, modern, and doubled carbon dioxide climates, J. Geophys. Res.- Atmos., 111, D10202, https://doi.org/10.1029/2005JD006653, 2006. Mahowald, N. M., Engelstaedter, S., Luo, C., Sealy, A., Artaxo, P., Benitez-Nelson, C., Bonnet, S., Chen, Y., Chuang, P. Y., and Cohen, D. D.: Atmospheric iron deposition: Global distribution, variability, and human perturbations*, Annu. Rev. Mar. Sci., 1, 245–278, 2009. Mahowald, N. M., Kloster, S., Engelstaedter, S., Moore, J. K., Mukhopadhyay, S., McConnell, J. R., Albani, S., Doney, S. C., Bhattacharya, A., Curran, M. A. J., Flanner, M. G., Hoffman, F. M., Lawrence, D. M., Lindsay, K., Mayewski, P. A., Neff, J., Rothenberg, D., Thomas, E., Thornton, P. E., and Zender, C. S.: Observed 20th century desert dust variability: impact on climate and biogeochemistry, Atmos. Chem. Phys., 10, 10875–10893, https://doi.org/10.5194/acp-10-10875-2010, 2010. Martin, J. H.: Glacial-interglacial CO2change: The iron hypothesis, Paleoceanography, 5, 1–13, 1990. Martin, J. H. and Fitzwater, S.: Iron deficiency limits phytoplankton growth in the north-east Pacific subarctic, Nature, 331, 947–975, 1988. Meskhidze, N., Chameides, W., Nenes, A., and Chen, G.: Iron mobilization in mineral dust: can anthropogenic SO2emissions affect ocean productivity?, Geophys. Res. Lett., 30, 2085, https://doi.org/10.1029/2003GL018035, 2003. Meskhidze, N., Chameides, W., and Nenes, A.: Dust and pollution: a recipe for enhanced ocean fertilization?, J. Geophys. Res.- Atmos., 110, D03301,https://doi.org/10.1029/2004JD005082, 2005. Mulitza, S., Heslop, D., Pittauerova, D., Fischer, H. W., Meyer, I., Stuut, J.-B., Zabel, M., Mollenhauer, G., Collins, J. A., and Kuhnert, H.: Increase in African dust flux at the onset of commercial agriculture in the Sahel region, Nature, 466, 226–228, 2010. Myriokefalitakis, S., Tsigaridis, K., Mihalopoulos, N., Sciare, J., Nenes, A., Kawamura, K., Segers, A., and Kanakidou, M.: In-cloud oxalate formation in the global troposphere: a 3-D modeling study, Atmos. Chem. Phys., 11, 5761–5782, https://doi.org/10.5194/acp-11-5761-2011, 2011. Myriokefalitakis, S., Daskalakis, N., Mihalopoulos, N., Baker, A. R., Nenes, A., and Kanakidou, M.: Changes in dissolved iron deposition to the oceans driven by human activity: a 3-D global modelling study, Biogeosciences, 12, 3973–3992, https://doi.org/10.5194/bg-12-3973-2015, 2015. Neale, R. B., Richter, J., Park, S., Lauritzen, P. H., Vavrus, S. J., Rasch, P. J., and Zhang, M.: The mean climate of the Community Atmosphere Model (CAM4) in forced SST and fully coupled experiments, J. Climate, 26, 5150–5168, 2013. Okin, G. S.: A new model of wind erosion in the presence of vegetation, J. Geophys. Res.-Earth, 113, F02S10, https://doi.org/10.1029/2007JF000758, 2008. Okin, G. S., Baker, A. R., Tegen, I., Mahowald, N. M., Dentener, F. J., Duce, R. A., Galloway, J. N., Hunter, K., Kanakidou, M., and Kubilay, N.: Impacts of atmospheric nutrient deposition on marine productivity: Roles of nitrogen, phosphorus, and iron, Global Biogeochem. Cy., 25, GB2022, https://doi.org/10.1029/2010GB003858, 2011. Paris, R., Desboeufs, K., and Journet, E.: Variability of dust iron solubility in atmospheric waters: Investigation of the role of oxalate organic complexation, Atmos. Environ., 45, 6510–6517, 2011. Pérez García-Pando, C., Miller, R. L., Perlwitz, J. P., Rodríguez, S., and Prospero, J. M.: Predicting the mineral composition of dust aerosols: Insights from elemental composition measured at the Izaña Observatory, Geophys. Res. Lett., 43, 10520–10529, https://doi.org/10.1002/2016GL069873, 2016. Perlwitz, J. P., Pérez García-Pando, C., and Miller, R. L.: Predicting the mineral composition of dust aerosols – Part 1: Representing key processes, Atmos. Chem. Phys., 15, 11593–11627, https://doi.org/10.5194/acp-15-11593-2015, 2015a. Perlwitz, J. P., Pérez García-Pando, C., and Miller, R. L.: Predicting the mineral composition of dust aerosols – Part 2: Model evaluation and identification of key processes with observations, Atmos. Chem. Phys., 15, 11629–11652, https://doi.org/10.5194/acp-15-11629-2015, 2015b. Rasch, P. and Kristjánsson, J.: A comparison of the CCM3 model climate using diagnosed and predicted condensate parameterizations, J. Climate, 11, 1587–1614, 1998. Rasch, P., Feichter, J., Law, K., Mahowald, N., Penner, J., Benkovitz, C., Genthon, C., Giannakopoulos, C., Kasibhatla, P., and Koch, D.: A comparison of scavenging and deposition processes in global models: results from the WCRP Cambridge Workshop of 1995, Tellus B, 52, 1025–1056, 2000. Saydam, A. and Senyuva, H.: Deserts: Can they be the potential suppliers of bioavailable iron?, Geophys. Res. Lett., 29, https://doi.org/10.1029/2001GL013562, 2002. Scanza, R. A., Mahowald, N., Ghan, S., Zender, C. S., Kok, J. F., Liu, X., Zhang, Y., and Albani, S.: Modeling dust as component minerals in the Community Atmosphere Model: development of framework and impact on radiative forcing, Atmos. Chem. Phys., 15, 537–561, https://doi.org/10.5194/acp-15-537-2015, 2015. Schulz, M., Balkanski, Y. J., Guelle, W., and Dulac, F.: Role of aerosol size distribution and source location in a threedimensional simulation of a Saharan dust episode tested against satellite-derived optical thickness, J. Geophys. Res.-Atmos., 103, 10579–10592, 1998. Seiler, W. and Crutzen, P. J.: Estimates of gross and net fluxes of carbon between the biosphere and the atmosphere from biomass burning, Climatic Change, 2, 207–247, 1980. Shi, Z., Krom, M. D., Bonneville, S., Baker, A. R., Jickells, T. D., and Benning, L. G.: Formation of iron nanoparticles and increase www.atmos-chem-phys.net/18/14175/2018/ Atmos. Chem. Phys., 18, 14175–14196, 2018 14196 R. A. Scanza et al.: Atmospheric processing of iron in mineral and combustion aerosols in iron reactivity in mineral dust during simulated cloud processing, Environ. Sci. Technol., 43, 6592–6596, 2009. Shi, Z., Krom, M. D., Jickells, T. D., Bonneville, S., Carslaw, K. S., Mihalopoulos, N., Baker, A. R., and Benning, L. G.: Impacts on iron solubility in the mineral dust by processes in the source region and the atmosphere: A review, Aeolian Res., 5, 21–42, 2012. Sholkovitz, E. R., Sedwick, P. N., Church, T. M., Baker, A. R., and Powell, C. F.: Fractional solubility of aerosol iron: Synthesis of a global-scale data set, Geochim. Cosmochim. Acta, 89, 173–189, 2012. Siefert, R. L., Johansen, A. M., Hoffmann, M. R., and Pehkonen, S. O.: Measurements of trace metal (Fe, Cu, Mn, Cr) oxidation states in fog and stratus clouds, J. Air Waste Manage., 48, 128– 143, 1998. Siegenthaler, U. and Sarmiento, J.: Atmospheric carbon dioxide and the ocean, Nature, 365, 119–125, 1993. Solmon, F., Chuang, P., Meskhidze, N., and Chen, Y.: Acidic processing of mineral dust iron by anthropogenic compounds over the north Pacific Ocean, J. Geophys. Res.-Atmos., 114, D02305, https://doi.org/10.1029/2008JD010417, 2009. Suarez, M. J., Rienecker, M., Todling, R., Bacmeister, J., Takacs, L., Liu, H., Gu, W., Sienkiewicz, M., Koster, R., and Gelaro, R.: The GEOS-5 Data Assimilation System-Documentation of Versions 5.0.1, 5.1.0, and 5.2.0, 2008. Weber, R. J., Guo, H., Russell, A. G., and Nenes, A.: High aerosol acidity despite declining atmospheric sulfate concentrations over the past 15 years, Nat. Geosci., 9, 282–285, https://doi.org/10.1038/ngeo2665, 2016. Winton, V., Bowie, A., Edwards, R., Keywood, M., Townsend, A., van der Merwe, P., and Bollhöfer, A.: Fractional iron solubility of atmospheric iron inputs to the Southern Ocean, Mar. Chem., 177, 20–32, 2015. Xu, N. and Gao, Y.: Characterization of hematite dissolution affected by oxalate coating, kinetics and pH, Appl. Geochem., 23, 783–793, 2008. Yu, J. Z., Huang, X.-F., Xu, J., and Hu, M.: When aerosol sulfate goes up, so does oxalate: Implication for the formation mechanisms of oxalate, Environ. Sci. Technol., 39, 128–133, 2005. Zender, C. S., Bian, H., and Newman, D.: Mineral Dust Entrainment and Deposition (DEAD) model: Description and 1990s dust climatology, J. Geophys. Res.-Atmos., 108, 4416, https://doi.org/10.1029/2002JD002775, 2003. Zhang, Y., Mahowald, N., Scanza, R. A., Journet, E., Desboeufs, K., Albani, S., Kok, J. F., Zhuang, G., Chen, Y., Cohen, D. D., Paytan, A., Patey, M. D., Achterberg, E. P., Engelbrecht, J. P., and Fomba, K. W.: Modeling the global emission, transport and deposition of trace elements associated with mineral dust, Biogeosciences, 12, 5771–5792, https://doi.org/10.5194/bg-125771-2015, 2015. Zhu, X., Prospero, J., and Millero, F.: Diel variability of soluble Fe (II) and soluble total Fe in North African dust in the trade winds at Barbados, J. Geophys. Res., 102, 21297–21305, 1997. Zhuang, G., Yi, Z., Duce, R. A., and Brown, P. R.: Link between iron and sulphur cycles suggested by detection of Fe (n) in remote marine aerosols, 1992. Atmos. Chem. Phys., 18, 14175–14196, 2018 www.atmos-chem-phys.net/18/14175/2018/