Evaluation of aerosol number concentrations in NorESM with improved nucleation parameterization

The Norwegian Earth System Model (NorESM) is evaluated against atmospheric observations of aerosol number concentrations. The model is extended to include an explicit mechanism for new particle formation and secondary organic aerosol (SOA) formation from biogenic precursors. Three nucleation mechanisms are included in NorESM: binary sulfuric acid, activation type, and organic nucleation. Ten model experiments are conducted to study the sensitivity of the simulated aerosol number concentrations to nucleation, SOA formation, black carbon size distribution and model meteorology. Simulated vertical profiles are evaluated against 12 flight campaigns. Comparison of monthly averaged aerosol number concentrations against 60 measurement sites reveals that the model with explicit nucleation and SOA scheme performs well in terms of correlation coefficient, R2 = 0.41 and a bias of −6%. NorESM generally overestimates the amplitude of the seasonal cycle, possibly due to underestimated sinks or exaggerated sensitivity to biogenic precursors

The effect of harmonized emissions on aerosol properties in global models – an AeroCom experiment

The effects of unified aerosol sources on global aerosol fields simulated by different models are examined in this paper. We compare results from two AeroCom experiments, one with different (ExpA) and one with unified emissions, injection heights, and particle sizes at the source (ExpB). Surprisingly, harmonization of aerosol sources has only a small impact on the simulated diversity for aerosol burden, and consequently optical properties, as the results are largely controlled by model-specific transport, removal, chemistry (leading to the formation of secondary aerosols) and parameterizations of aerosol microphysics (e.g. the split between deposition pathways) and to a lesser extent on the spatial and temporal distributions of the (precursor) emissions. The burdens of black carbon and especially sea salt become more coherent in ExpB only, because the large ExpA diversity for these two species was caused by few outliers. The experiment also indicated that despite prescribing emission fluxes and size distributions, ambiguities in the implementation in individual models can lead to substantial differences. These results indicate the need for a better understanding of aerosol life cycles at process level (including spatial dispersal and interaction with meteorological parameters) in order to obtain more reliable results from global aerosol simulations. This is particularly important as such model results are used to assess the consequences of specific air pollution abatement strategies

Analysis and quantification of the diversities of aerosol life cycles within AeroCom

Simulation results of global aerosol models have been assembled in the framework of the AeroCom intercomparison exercise. In this paper, we analyze the life cycles of dust, sea salt, sulfate, black carbon and particulate organic matter as simulated by sixteen global aerosol models. The diversities among the models for the sources and sinks, burdens, particle sizes, water uptakes, and spatial dispersals have been established. These diversities have large consequences for the calculated radiative forcing and the aerosol concentrations at the surface. The AeroCom all-models-average emissions are dominated by the mass of sea salt (SS), followed by dust (DU), sulfate (SO_4), particulate organic matter (POM), and finally black carbon (BC). Interactive parameterizations of the emissions and contrasting particles sizes of SS and DU lead generally to higher diversities of these species, and for total aerosol. The lower diversity of the emissions of the fine aerosols, BC, POM, and SO_4, is due to the use of similar emission inventories, and does therefore not necessarily indicate a better understanding of their sources. The diversity of SO_4-sources is mainly caused by the disagreement on depositional loss of precursor gases and on chemical production. The diversities of the emissions are passed on to the burdens, but the latter are also strongly affected by the model-specific treatments of transport and aerosol processes. The burdens of dry masses decrease from largest to smallest: DU, SS, SO_4, POM, and BC. The all-models-average residence time is shortest for SS with about half a day, followed by S_O4 and DU with four days, and POM and BC with six and seven days, respectively. The wet deposition rate is controlled by the solubility and increases from DU, BC, POM to SO_4 and SS. It is the dominant sink for SO_4, BC, and POM, and contributes about one third to the total removal rate coefficients of SS and DU species. For SS and DU we find high diversities for the removal rate coefficients and deposition pathways. Models do neither agree on the split between wet and dry deposition, nor on that between sedimentation and turbulent dry Deposition. We diagnose an extremely high diversity for the uptake of ambient water vapor that influences the particle size and thus the sink rate coefficients. Furthermore, we find little agreement among the model results for the partitioning of wet removal into scavenging by convective and stratiform rain. Large differences exist for aerosol dispersal both in the vertical and in the horizontal direction. In some models, a minimum of total aerosol concentration is simulated at the surface. Aerosol dispersal is most pronounced for SO4 and BC and lowest for SS. Diversities are higher for meridional than for vertical dispersal, they are similar for a given species and highest for SS and DU. For these two components we do not find a correlation between vertical and meridional aerosol dispersal. In addition the degree of dispersals of SS and DU is not related to their residence times. SO_4, BC, and POM, however, show increased meridional dispersal in models with larger vertical dispersal, and dispersal is larger for longer simulated residence times

An AeroCom initial assessment – optical properties in aerosol component modules of global models

The AeroCom exercise diagnoses multi-component aerosol modules in global modeling. In an initial assessment simulated global distributions for mass and mid-visible aerosol optical thickness (aot) were compared among 20 different modules. Model diversity was also explored in the context of previous comparisons. For the component combined aot general agreement has improved for the annual global mean. At 0.11 to 0.14, simulated aot values are at the lower end of global averages suggested by remote sensing from ground (AERONET ca. 0.135) and space (satellite composite ca. 0.15). More detailed comparisons, however, reveal that larger differences in regional distribution and significant differences in compositional mixture remain. Of particular concern are large model diversities for contributions by dust and carbonaceous aerosol, because they lead to significant uncertainty in aerosol absorption (aab). Since aot and aab, both, influence the aerosol impact on the radiative energy-balance, the aerosol (direct) forcing uncertainty in modeling is larger than differences in aot might suggest. New diagnostic approaches are proposed to trace model differences in terms of aerosol processing and transport: These include the prescription of common input (e.g. amount, size and injection of aerosol component emissions) and the use of observational capabilities from ground (e.g. measurements networks) or space (e.g. correlations between aerosol and clouds)

Fast and slow precipitation responses to individual climate forcers: a PDRMIP multi-model study

Precipitation is expected to respond differently to various drivers of anthropogenic climate change. We present the first results from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP), where nine global climate models have perturbed CO2, CH4, black carbon, sulfate, and solar insolation. We divide the resulting changes to global mean and regional precipitation into fast responses that scale with changes in atmospheric absorption and slow responses scaling with surface temperature change. While the overall features are broadly similar between models, we find significant regional intermodel variability, especially over land. Black carbon stands out as a component that may cause significant model diversity in predicted precipitation change. Processes linked to atmospheric absorption are less consistently modeled than those linked to top-of-atmosphere radiative forcing. We identify a number of land regions where the model ensemble consistently predicts that fast precipitation responses to climate perturbations dominate over the slow, temperature-driven responses

Rapid Adjustments Cause Weak Surface Temperature Response to Increased Black Carbon Concentrations

We investigate the climate response to increased concentrations of black carbon (BC), as part of the Precipitation Driver Response Model Intercomparison Project (PDRMIP). A tenfold increase in BC is simulated by nine global coupled‐climate models, producing a model median effective radiative forcing of 0.82 (ranging from 0.41 to 2.91) W m⁻², and a warming of 0.67 (0.16 to 1.66) K globally and 1.24 (0.26 to 4.31) K in the Arctic. A strong positive instantaneous radiative forcing (median of 2.10 W m⁻² based on five of the models) is countered by negative rapid adjustments (−0.64 W m⁻² for the same five models), which dampen the total surface temperature signal. Unlike other drivers of climate change, the response of temperature and cloud profiles to the BC forcing is dominated by rapid adjustments. Low‐level cloud amounts increase for all models, while higher‐level clouds are diminished. The rapid temperature response is particularly strong above 400 hPa, where increased atmospheric stabilization and reduced cloud cover contrast the response pattern of the other drivers. In conclusion, we find that this substantial increase in BC concentrations does have considerable impacts on important aspects of the climate system. However, some of these effects tend to offset one another, leaving a relatively small median global warming of 0.47 K per W m⁻²—about 20% lower than the response to a doubling of CO₂. Translating the tenfold increase in BC to the present‐day impact of anthropogenic BC (given the emissions used in this work) would leave a warming of merely 0.07 K

Microphysical explanation of the RH-dependent water affinity of biogenic organic aerosol and its importance for climate

This is the final version of the article. Available from American Geophysical Union via the DOI in this record.A large fraction of atmospheric organic aerosol (OA) originates from natural emissions that are oxidized in the atmosphere to form secondary organic aerosol (SOA). Isoprene (IP) and monoterpenes (MT) are the most important precursors of SOA originating from forests. The climate impacts from OA are currently estimated through parameterizations of water uptake that drastically simplify the complexity of OA. We combine laboratory experiments, thermodynamic modeling, field observations, and climate modeling to (1) explain the molecular mechanisms behind RH-dependent SOA water-uptake with solubility and phase separation; (2) show that laboratory data on IP- and MT-SOA hygroscopicity are representative of ambient data with corresponding OA source profiles; and (3) demonstrate the sensitivity of the modeled aerosol climate effect to assumed OA water affinity. We conclude that the commonly used single-parameter hygroscopicity framework can introduce significant error when quantifying the climate effects of organic aerosol. The results highlight the need for better constraints on the overall global OA mass loadings and its molecular composition, including currently underexplored anthropogenic and marine OA sources.The data presented in the paper will be available through the Bolin Centre database (http://bolin.su.se/data/). The EC H2020 European Research Council ERC (ERC-StGATMOGAIN-278277 and ERC-StG-QAPPA-335478) and integrated project 641816 CRESCENDO Svenska Forskningsrådet Formas (Swedish Research Council Formas) (2015-749), Knut och Alice Wallenbergs Stiftelse (Knut and Alice Wallenberg Foundation Wallenberg Fellowship AtmoRemove), Academy of Finland (grants 272041 and 259005), Natural Environment Research Council (NERC grants NE/M003531/1 and NE/J02175X/1), Norwegian Research Council (EVA grant 229771), Natural Sciences and Engineering Research Council of Canada (NSERC, grant RGPIN/04315-2014), National Science Foundation (NSF, grants ATM-1242258, AGS-1242932, and AGS-1360834), U.S. Environmental Protection Agency (EPA, STAR grant R835410), National Oceanic and Atmospheric Administration (NOAA, CPO award 538NA10OAR4310102), Electric Power Research Institute (EPRI, grant 10004734), U.S. Department of Energy (DOE, grants BER/ASR DE-SC0016559 and DE-SC0012792), Georgia Institute of Technology, and NordForsk (Nordic Centre of Excellence eSTICC) are gratefully acknowledged for funding. The climate model simulations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) at the National Supercomputing Centre. Benjamin Murphy is acknowledged for useful discussions

The Norwegian Earth System Model, NorESM1-M – Part 1: Description and basic evaluation of the physical climate

The core version of the Norwegian Climate Center's Earth System Model, named NorESM1-M, is presented. The NorESM family of models are based on the Community Climate System Model version 4 (CCSM4) of the University Corporation for Atmospheric Research, but differs from the latter by, in particular, an isopycnic coordinate ocean model and advanced chemistry–aerosol–cloud–radiation interaction schemes. NorESM1-M has a horizontal resolution of approximately 2° for the atmosphere and land components and 1° for the ocean and ice components. NorESM is also available in a lower resolution version (NorESM1-L) and a version that includes prognostic biogeochemical cycling (NorESM1-ME). The latter two model configurations are not part of this paper. Here, a first-order assessment of the model stability, the mean model state and the internal variability based on the model experiments made available to CMIP5 are presented. Further analysis of the model performance is provided in an accompanying paper (Iversen et al., 2013), presenting the corresponding climate response and scenario projections made with NorESM1-M

Historical Changes and Reasons for Model Differences in Anthropogenic Aerosol Forcing in CMIP6

The Radiative Forcing Model Intercomparison Project (RFMIP) allows estimates of effective radiative forcing (ERF) in the Coupled Model Intercomparison Project phase six (CMIP6). We analyze the RFMIP output, including the new experiments from models that use the same parameterization for anthropogenic aerosols (RFMIP-SpAer), to characterize and better understand model differences in aerosol ERF. We find little changes in the aerosol ERF for 1970–2014 in the CMIP6 multi-model mean, which implies greenhouse gases primarily explain the positive trend in the total anthropogenic ERF. Cloud-mediated effects dominate the present-day aerosol ERF in most models. The results highlight a regional increase in marine cloudiness due to aerosols, despite suppressed cloud lifetime effects in that RFMIP-SpAer experiment. Negative cloud-mediated effects mask positive direct effects in many models, which arise from strong anthropogenic aerosol absorption. The findings suggest opportunities to better constrain simulated ERF by revisiting the optical properties and long-range transport of aerosols