EC-Earth3-AerChem: A global climate model with interactive aerosols and atmospheric chemistry participating in CMIP6

This paper documents the global climate model EC-Earth3-AerChem, one of the members of the EC-Earth3 family of models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6). EC-Earth3-AerChem has interactive aerosols and atmospheric chemistry and contributes to the Aerosols and Chemistry Model Intercomparison Project (AerChemMIP). In this paper, we give an overview of the model, describe in detail how it differs from the other EC-Earth3 configurations, and outline the new features compared with the previously documented version of the model (EC-Earth 2.4). We explain how the model was tuned and spun up under preindustrial conditions and characterize the model's general performance on the basis of a selection of coupled simulations conducted for CMIP6. The net energy imbalance at the top of the atmosphere in the preindustrial control simulation is on average -0.09 W m-2 with a standard deviation due to interannual variability of 0.25 W m-2, showing no significant drift. The global surface air temperature in the simulation is on average 14.08 ∼ C with an interannual standard deviation of 0.17 ∼ C, exhibiting a small drift of 0.015 ± 0.005 ∼ C per century. The model's effective equilibrium climate sensitivity is estimated at 3.9 ∼ C, and its transient climate response is estimated at 2.1 ∼ C. The CMIP6 historical simulation displays spurious interdecadal variability in Northern Hemisphere temperatures, resulting in a large spread across ensemble members and a tendency to underestimate observed annual surface temperature anomalies from the early 20th century onwards. The observed warming of the Southern Hemisphere is well reproduced by the model. Compared with the ECMWF (European Centre for Medium-Range Weather Forecasts) Reanalysis version 5 (ERA5), the surface air temperature climatology for 1995-2014 has an average bias of -0.86 ± 0.05 ∼ C with a standard deviation across ensemble members of 0.35 ∼ C in the Northern Hemisphere and 1.29 ± 0.02 ∼ C with a corresponding standard deviation of 0.05 ∼ C in the Southern Hemisphere. The Southern Hemisphere warm bias is largely caused by errors in shortwave cloud radiative effects over the Southern Ocean, a deficiency of many climate models. Changes in the emissions of near-term climate forcers (NTCFs) have significant effects on the global climate from the second half of the 20th century onwards. For the SSP3-7.0 Shared Socioeconomic Pathway, the model gives a global warming at the end of the 21st century (2091-2100) of 4.9 ∼ C above the preindustrial mean. A 0.5 ∼ C stronger warming is obtained for the AerChemMIP scenario with reduced emissions of NTCFs. With concurrent reductions of future methane concentrations, the warming is projected to be reduced by 0.5 ∼ C

Evaluation of global simulations of aerosol particle and cloud condensation nuclei number, with implications for cloud droplet formation

A total of 16 global chemistry transport models and general circulation models have participated in this study; 14 models have been evaluated with regard to their ability to reproduce the near-surface observed number concentration of aerosol particles and cloud condensation nuclei (CCN), as well as derived cloud droplet number concentration (CDNC). Model results for the period 2011-2015 are compared with aerosol measurements (aerosol particle number, CCN and aerosol particle composition in the submicron fraction) from nine surface stations located in Europe and Japan. The evaluation focuses on the ability of models to simulate the average across time state in diverse environments and on the seasonal and short-term variability in the aerosol properties. There is no single model that systematically performs best across all environments represented by the observations. Models tend to underestimate the observed aerosol particle and CCN number concentrations, with average normalized mean bias (NMB) of all models and for all stations, where data are available, of -24 % and -35 % for particles with dry diameters > 50 and > 120 nm, as well as -36 % and -34 % for CCN at supersaturations of 0.2 % and 1.0 %, respectively. However, they seem to behave differently for particles activating at very low supersaturations (<0.1 %) than at higher ones. A total of 15 models have been used to produce ensemble annual median distributions of relevant parameters. The model diversity (defined as the ratio of standard deviation to mean) is up to about 3 for simulated N-3 (number concentration of particles with dry diameters larger than 3 nm) and up to about 1 for simulated CCN in the extra-polar regions. A global mean reduction of a factor of about 2 is found in the model diversity for CCN at a supersaturation of 0.2 % (CCN0.2) compared to that for N-3, maximizing over regions where new particle formation is important. An additional model has been used to investigate potential causes of model diversity in CCN and bias compared to the observations by performing a perturbed parameter ensemble (PPE) accounting for uncertainties in 26 aerosol-related model input parameters. This PPE suggests that biogenic secondary organic aerosol formation and the hygroscopic properties of the organic material are likely to be the major sources of CCN uncertainty in summer, with dry deposition and cloud processing being dominant in winter. Models capture the relative amplitude of the seasonal variability of the aerosol particle number concentration for all studied particle sizes with available observations (dry diameters larger than 50, 80 and 120 nm). The short-term persistence time (on the order of a few days) of CCN concentrations, which is a measure of aerosol dynamic behavior in the models, is underestimated on average by the models by 40 % during winter and 20 % in summer. In contrast to the large spread in simulated aerosol particle and CCN number concentrations, the CDNC derived from simulated CCN spectra is less diverse and in better agreement with CDNC estimates consistently derived from the observations (average NMB -13 % and -22 % for updraft velocities 0.3 and 0.6 m s(-1), respectively). In addition, simulated CDNC is in slightly better agreement with observationally derived values at lower than at higher updraft velocities (index of agreement 0.64 vs. 0.65). The reduced spread of CDNC compared to that of CCN is attributed to the sublinear response of CDNC to aerosol particle number variations and the negative correlation between the sensitivities of CDNC to aerosol particle number concentration (partial derivative N-d/partial derivative N-a) and to updraft velocity (partial derivative N-d/partial derivative w). Overall, we find that while CCN is controlled by both aerosol particle number and composition, CDNC is sensitive to CCN at low and moderate CCN concentrations and to the updraft velocity when CCN levels are high. Discrepancies are found in sensitivities partial derivative N-d/partial derivative N-a and partial derivative N-d/partial derivative w; models may be predisposed to be too "aerosol sensitive" or "aerosol insensitive" in aerosol-cloud-climate interaction studies, even if they may capture average droplet numbers well. This is a subtle but profound finding that only the sensitivities can clearly reveal and may explain intermodel biases on the aerosol indirect effect.Peer reviewe

Evaluation of Global Simulations of Aerosol Particle and Cloud Condensation Nuclei Number, with Implications for Cloud Droplet Formation

A total of 16 global chemistry transport models and general circulation models have participated in this study; 14 models have been evaluated with regard to their ability to reproduce the near-surface observed number concentration of aerosol particles and cloud condensation nuclei (CCN), as well as derived cloud droplet number concentration (CDNC). Model results for the period 2011-2015 are compared with aerosol measurements (aerosol particle number, CCN and aerosol particle composition in the submicron fraction) from nine surface stations located in Europe and Japan. The evaluation focuses on the ability of models to simulate the average across time state in diverse environments and on the seasonal and short-term variability in the aerosol properties. There is no single model that systematically performs best across all environments represented by the observations. Models tend to underestimate the observed aerosol particle and CCN number concentrations, with average normalized mean bias (NMB) of all models and for all stations, where data are available, of -24% and -35% for particles with dry diameters > 50 and > 120nm, as well as -36% and -34% for CCN at supersaturations of 0.2% and 1.0%, respectively. However, they seem to behave differently for particles activating at very low supersaturations (< 0.1%) than at higher ones. A total of 15 models have been used to produce ensemble annual median distributions of relevant parameters. The model diversity (defined as the ratio of standard deviation to mean) is up to about 3 for simulated N3 (number concentration of particles with dry diameters larger than 3 nm) and up to about 1 for simulated CCN in the extra-polar regions. A global mean reduction of a factor of about 2 is found in the model diversity for CCN at a supersaturation of 0.2% (CCN(0.2)) compared to that for N3, maximizing over regions where new particle formation is important. An additional model has been used to investigate potential causes of model diversity in CCN and bias compared to the observations by performing a perturbed parameter ensemble (PPE) accounting for uncertainties in 26 aerosol-related model input parameters. This PPE suggests that biogenic secondary organic aerosol formation and the hygroscopic properties of the organic material are likely to be the major sources of CCN uncertainty in summer, with dry deposition and cloud processing being dominant in winter. Models capture the relative amplitude of the seasonal variability of the aerosol particle number concentration for all studied particle sizes with available observations (dry diameters larger than 50, 80 and 120nm). The short-term persistence time (on the order of a few days) of CCN concentrations, which is a measure of aerosol dynamic behavior in the models, is underestimated on average by the models by 40% during winter and 20% in summer