Abstract. Understanding how multiphase processes affect the iron-containing aerosol cycle is key to predicting ocean biogeochemistry changes and hence the feedback effects on climate. For this work, the EC-Earth Earth system model in its climate–chemistry configuration is used to simulate the global atmospheric oxalate (OXL), sulfate (SO42-), and iron (Fe) cycles after incorporating a comprehensive representation of the multiphase chemistry in cloud droplets and aerosol water. The model considers a detailed gas-phase chemistry scheme, all major aerosol components, and the partitioning of gases in aerosol and atmospheric water phases. The dissolution of Fe-containing aerosols accounts kinetically for the solution's acidity, oxalic acid, and irradiation. Aerosol acidity is explicitly calculated in the model, both for accumulation and coarse modes, accounting for thermodynamic processes involving inorganic and crustal species from sea salt and dust. Simulations for present-day conditions (2000–2014) have been carried out with both EC-Earth and the atmospheric composition component of the model in standalone mode driven by meteorological fields from ECMWF's ERA-Interim reanalysis. The calculated global budgets are presented and the links between the (1) aqueous-phase processes, (2) aerosol dissolution, and (3) atmospheric composition are demonstrated and quantified. The model results are supported by comparison to available observations. We obtain an average global OXL net chemical production of 12.615 ± 0.064 Tg yr−1 in EC-Earth, with glyoxal being by far the most important precursor of oxalic acid. In comparison to the ERA-Interim simulation, differences in atmospheric dynamics and the simulated weaker oxidizing capacity in EC-Earth overall result in a ∼ 30 % lower OXL source. On the other hand, the more explicit representation of the aqueous-phase chemistry in EC-Earth compared to the previous versions of the model leads to an overall ∼ 20 % higher sulfate production, but this is still well correlated with atmospheric observations. The total Fe dissolution rate in EC-Earth is calculated at 0.806 ± 0.014 Tg yr−1 and is added to the primary dissolved Fe (DFe) sources from dust and combustion aerosols in the model (0.072 ± 0.001 Tg yr−1). The simulated DFe concentrations show a satisfactory comparison with available observations, indicating an atmospheric burden of ∼0.007 Tg, resulting in an overall atmospheric deposition flux into the global ocean of 0.376 ± 0.005 Tg yr−1, which is well within the range reported in the literature. All in all, this work is a first step towards the development of EC-Earth into an Earth system model with fully interactive bioavailable atmospheric Fe inputs to the marine biogeochemistry component of the model.

Abstract. Understanding how multiphase processes affect the iron-containing aerosol cycle is key to predict ocean biogeochemistry changes and hence the feedback effects on climate. For this work, the EC-Earth Earth system model in its climate-chemistry configuration is used to simulate the global atmospheric oxalate (OXL), sulfate (SO42−), and iron (Fe) cycles, after incorporating a comprehensive representation of the multiphase chemistry in cloud droplets and aerosol water. The model considers a detailed gas-phase chemistry scheme, all major aerosol components, and the partitioning of gases in aerosol and atmospheric water phases. The dissolution of Fe-containing aerosols accounts kinetically for the solution’s acidity, oxalic acid, and irradiation. Aerosol acidity is explicitly calculated in the model, both for accumulation and coarse modes, accounting for thermodynamic processes involving inorganic and crustal species from sea salt and dust. Simulations for present-day conditions (2000–2014) have been carried out with both EC-Earth and the atmospheric composition component of the model in standalone mode driven by meteorological fields from ECMWF’s ERA-Interim reanalysis. The calculated global budgets are presented and the links between the 1) aqueous-phase processes, 2) aerosol dissolution, and 3) atmospheric composition, are demonstrated and quantified. The model results are supported by comparison to available observations. We obtain an average global OXL net chemical production of 12.61 ± 0.06 Tg yr−1 in EC-Earth, with glyoxal being by far the most important precursor of oxalic acid. In comparison to the ERA-Interim simulation, differences in atmospheric dynamics as well as the simulated weaker oxidizing capacity in EC-Earth result overall in a ~30 % lower OXL source. On the other hand, the more explicit representation of the aqueous-phase chemistry in EC-Earth compared to the previous versions of the model leads to an overall ~20 % higher sulfate production, but still well correlated with atmospheric observations. The total Fe dissolution rate in EC-Earth is calculated at 0.806 ± 0.014 Tg Fe yr−1 and is added to the primary dissolved Fe (DFe) sources from dust and combustion aerosols in the model (0.072 ± 0.001 Tg Fe yr−1). The simulated DFe concentrations show a satisfactory comparison with available observations, indicating an atmospheric burden of ∼0.007 Tg Fe, and overall resulting in an atmospheric deposition flux into the global ocean of 0.376 ± 0.005 Tg Fe yr−1, well within the range reported in the literature. All in all, this work is a first step towards the development of EC-Earth into an Earth System Model with fully interactive bioavailable atmospheric Fe inputs to the marine biogeochemistry component of the model.

Changes in atmospheric iron (Fe) deposition to the open ocean affect net primary productivity, nitrogen fixation, and carbon uptake rates. We investigate the changes in soluble Fe (SFe) deposition from the pre-industrial period to the late 21st century using the EC-Earth3-Iron Earth System model, which stands out for its comprehensive representation of the atmospheric oxalate, sulfate, and Fe cycles. We show how anthropogenic activity has modified the magnitude and spatial distribution of SFe deposition by increasing combustion Fe emissions along with atmospheric acidity and oxalate levels. We find that SFe deposition has doubled since the early Industrial Era using the Coupled Model Intercomparison Project Phase 6 (CMIP6) emission inventory, with acidity being the main solubilization pathway for dust Fe, and ligand-promoted (oxalate) processing dominating the solubilization of combustion Fe. We project a global SFe deposition increase of 40% by the late 21st century relative to present day under Shared Socioeconomic Pathway (SSP) 3-7.0, which assumes weak climate change mitigation policies. In contrast, sustainable and business-as-usual SSPs (1-2.6 and 2-4.5) result in 35% and 10% global decreases, respectively. Despite these differences, SFe deposition consistently increases and decreases across SSPs over the (high nutrient low chlorophyl) equatorial Pacific and Southern Ocean (SO), respectively. Future changes in dust and wildfires with climate remains a key challenge for constraining SFe projections. We show that the equatorial Pacific and the SO would be sensitive not only to changes in Australian or South American dust emissions, but also to those in North Africa.