Model evaluation of short-lived climate forcers for the Arctic Monitoring and Assessment Programme: a multi-species, multi-model study

Abstract

While carbon dioxide is the main cause for global warming, modeling short-lived climate forcers (SLCFs) such as methane, ozone, and particles in the Arctic allows us to simulate near-term climate and health impacts for a sensitive, pristine region that is warming at 3 times the global rate. Atmospheric modeling is critical for understanding the long-range transport of pollutants to the Arctic, as well as the abundance and distribution of SLCFs throughout the Arctic atmosphere. Modeling is also used as a tool to determine SLCF impacts on climate and health in the present and in future emissions scenarios.

In this study, we evaluate 18 state-of-the-art atmospheric and Earth system models by assessing their representation of Arctic and Northern Hemisphere atmospheric SLCF distributions, considering a wide range of different chemical species (methane, tropospheric ozone and its precursors, black carbon, sulfate, organic aerosol, and particulate matter) and multiple observational datasets. Model simulations over 4 years (2008–2009 and 2014–2015) conducted for the 2022 Arctic Monitoring and Assessment Programme (AMAP) SLCF assessment report are thoroughly evaluated against satellite, ground, ship, and aircraft-based observations. The annual means, seasonal cycles, and 3-D distributions of SLCFs were evaluated using several metrics, such as absolute and percent model biases and correlation coefficients. The results show a large range in model performance, with no one particular model or model type performing well for all regions and all SLCF species. The multi-model mean (mmm) was able to represent the general features of SLCFs in the Arctic and had the best overall performance. For the SLCFs with the greatest radiative impact (CH${}_{\mathrm{4}}$, O${}_{\mathrm{3}}$, BC, and SO${}_{\mathrm{4}}^{\mathrm{2}-}$), the mmm was within $\pm$25 % of the measurements across the Northern Hemisphere. Therefore, we recommend a multi-model ensemble be used for simulating climate and health impacts of SLCFs.

Of the SLCFs in our study, model biases were smallest for CH${}_{\mathrm{4}}$ and greatest for OA. For most SLCFs, model biases skewed from positive to negative with increasing latitude. Our analysis suggests that vertical mixing, long-range transport, deposition, and wildfires remain highly uncertain processes. These processes need better representation within atmospheric models to improve their simulation of SLCFs in the Arctic environment. As model development proceeds in these areas, we highly recommend that the vertical and 3-D distribution of SLCFs be evaluated, as that information is critical to improving the uncertain processes in models.

1 Introduction

The Arctic atmosphere is warming 3 times more quickly than the global average . Arctic warming is a manifestation of global warming, and the main driver for this is the increasing carbon dioxide (CO${}_{\mathrm{2}}$) radiative forcing . Arctic warming is amplified by sea ice and snow feedbacks and affected by local radiative forcings in the Arctic, including radiative forcings by short-lived climate forcers (SLCFs), such as methane, black carbon, and tropospheric ozone . The remote pristine Arctic environment is sensitive to the long-range transport of atmospheric pollutants and deposition . At the same time, it is difficult to carry out in situ measurements and satellite observations over the Arctic. The majority of the Arctic surface is ocean covered with sea ice that is usually adrift for most of the year. The Arctic environment is also harsh. These aspects have historically kept surface-based measurements sparse. The overwhelming majority of the satellite observations either depend on the visible spectrum, are limited by the presence of clouds, or have very low sensitivity in the lower troposphere where the atmospheric processes mainly determine the fate of the pollutants. Many satellite measurements also do not have good coverage in the Arctic, given their orbital parameters or problems measuring areas with high albedo .

Modeling the Arctic atmosphere comes with its own challenges due to extreme meteorological conditions, its great distance from major global pollution sources, poorly known local emissions, high gradients in physical and chemical fields, and a singularity in some model grids at the pole. Models have been improving in the last 2 decades, but many models still have inaccurate results in the Arctic . That said, there has recently been a number of improvements in numerous models that have allowed for better representation of certain processes . In this study, model simulations for the 2021 Arctic Monitoring and Assessment Programme (AMAP) SLCF assessment report have been thoroughly evaluated by comparison to several freely available observational datasets in the Northern Hemisphere and assessed in more detail in the Arctic. In order to support the integrated assessment of climate and human health for AMAP, 6 SLCF species (methane – CH${}_{\mathrm{4}}$, ozone – O${}_{\mathrm{3}}$, black carbon – BC, sulfate – SO${}_{\mathrm{4}}^{\mathrm{2}-}$, organic aerosol – OA, and fine particulate matter – PM${}_{\mathrm{2.5}}$) and 2 O${}_{\mathrm{3}}$ precursors (carbon monoxide – CO; nitrogen dioxide – NO${}_{\mathrm{2}}$) from 18 atmospheric or Earth system models are compared to numerous observational datasets (from three satellite instruments, seven monitoring networks, and nine measurement campaigns) for 4 years (2008–2009 and 2014–2015), with the goal of answering the following questions.

• 1.

  How well do the AMAP SLCF models perform in the context of measurements and their associated uncertainty?

• What do the best-performing models have in common?

• Are there regional patterns in the model biases?

• Are there patterns in the model biases between SLCF species?

• 2.

  How does the model performance impact model applications, such as simulated climate and health impacts?

• 3.

  What processes should be improved or studied further for better model performance?

Out of scope of this study are any sensitivity tests by the models to assess different components of model errors. Also out of scope are the models' simulations of aerosol optical properties and cloud properties (e.g., cloud fraction, cloud droplet number concentration, cloud scavenging), though those parameters do have a large impact on climate and a tight relationship with some SLCFs. Their initial evaluation can be found in (chap. 7). Estimates of effective radiative forcings of SLCFs in the Arctic by the AMAP participating models are also provided elsewhere .

The next section summarizes the models used in this study, with more information in the Appendix. Section  summarizes the measurements used for model evaluation. Section  presents our model evaluation for each SLCF species, followed by a summary of all SLCFs. Finally, Sect.  is the conclusion where the questions posed above are answered.

2 Models

In this section we briefly describe the models used for the AMAP SLCF study and refer the reader to Appendix  for individual model descriptions and further information. All models were run globally with the same anthropogenic emissions dataset (see Sect. ), and most were run for the years 2008–2009 (as was done for the 2015 AMAP assessment report) and 2014–2015 (to evaluate more recent model results) inclusive for this evaluation, as these were years with numerous Arctic measurements. Unless otherwise indicated, all model output was monthly-averaged.

The models used for this study are summarized in Table . As is shown in the table, not all models provided all SLCF species, and not all models provided all 4 years. There were eight chemical transport models (CTMs), two chemistry–climate models (CCMs), three global climate models (GCMs), and five Earth system models (ESMs). Many models used specified or nudged meteorology, which allows the day-to-day variability of the model meteorology to be more closely aligned with the historical evolution of the atmosphere than occurs in a free-running model. The ERA-Interim reanalysis was the most commonly used meteorology (in 7 out of 18 models), but some were free-running (simulating their own meteorology) and some used other reanalysis products (Table ).

Table 1

Summary of models used in this study. GCM: global climate model, CCM: chemistry–climate model, ESM: Earth system model, CTM: chemical transport model.

Name	Type	Meteorology	Simulation period	SLCF output	Primary reference(s)
CanAM5-PAM	GCM	nudged to ERA-	1990–2015	BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$,
		Interim reanalysis		AOD, AAOD, AE
CESM2.0	ESM	free-running	2008–2009, 2014–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, BC,
				SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$, AOD,
				AAOD, AE
CIESM-MAM7	GCM	nudged to ERA-	1990–2015	BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$,
		Interim reanalysis		AOD
CMAM	CCM	nudged to ERA-	1990–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$
		Interim reanalysis
DEHM	CTM	nudged to ERA-	1990–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$,
		Interim reanalysis		BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$,
				AOD, AAOD, AE
ECHAM6-SALSA	GCM	nudged to ERA-	2008–2009, 2014–2015	BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$,
		Interim reanalysis		AOD, AAOD, AE
EMEP MSC-W	CTM	driven by 3-hourly	1990–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$,
		ECMWF met		BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$, AOD
FLEXPART	Lagrangian CTM	driven by 3-hourly	2014–2015	BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$
		ECMWF met
GEM-MACH	online CTM	driven by GEM	2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, BC,
		numerical forecast		SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$
GEOS-Chem	CTM	Driven by GEOS	2008–2009, 2014–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$,
		meteorology		BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$,
				AOD, AAOD, AE
GISS-E2.1	ESM	nudged to NCEP	1990–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$,
		reanalysis		BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$,
				AOD, AAOD, AE
MATCH	CTM	ERA-Interim	2008–2009, 2014–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, BC,
		reanalysis		SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$, AOD, AAOD, AE
MATCH-SALSA-RCA4	CCM	RCA4	2008–2009, 2014–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, BC,
				SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$, AOD, AAOD, AE
MRI-ESM2	ESM	nudged to JRA55	1990–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$,
		reanalysis		BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$,
				AOD, AAOD
NorESM1-happi	ESM	free-running	2008–2009, 2014–2015	BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, AOD,
				AAOD, AE
Oslo-CTM	CTM	driven by 3-hourly	2008–2009, 2014–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$,
		ECMWF meteorology		BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$
UKESM1	CCM $\mathit{\&}$ ESM	nudged to ERA-	1990–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, CH${}_{\mathrm{4}}$,
		Interim reanalysis		BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, PM${}_{\mathrm{2.5}}$,
				AOD, AAOD, AE
WRF-Chem	CCM $\mathit{\&}$ CTM	nudged to NCEP FNL	2014–2015	O${}_{\mathrm{3}}$, CO, NO${}_{\mathrm{2}}$, BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$,
		reanalysis		OA, PM${}_{\mathrm{2.5}}$, AOD, AAOD, AE

2.1 Emissions

All models used the same anthropogenic emissions dataset, which is called ECLIPSE (Evaluating the Climate and Air Quality Impacts of Short-Lived Pollutants) v6B. These emissions were created using the IIASA-GAINS (International Institute for Applied Systems Analysis – Greenhouse gas – Air pollution Interactions and Synergies) model , which provides emissions of long-lived greenhouse gases and shorter-lived species in a consistent framework. These historical emissions were provided for the years 1990 to 2015 at 5-year intervals, as well as the years 2008–2009 and 2014. Those models that simulated the 1990–2015 time period linearly interpolated the emissions for the years in between. The ECLIPSEv6b emissions include many pollutants, such as CH${}_{\mathrm{4}}$, CO, NO${}_x$, BC, and SO${}_{\mathrm{2}}$. They include the significant sulfur emission reductions that have taken place since the 1980s . Global anthropogenic BC emissions are estimated to be 6.5 Tg in 2010 and 5.9 Tg in 2020, and global anthropogenic SO${}_{\mathrm{2}}$ emissions are estimated to be 90 Tg in 2010 but declined significantly over the subsequent decade to 50 Tg . The reductions are mainly due to stringent emissions standards in the energy and industrial sectors, as well as reduced coal use in the residential sector . Global anthropogenic methane emissions were 340 Tg in 2015 and 350 Tg in 2020, and they are expected to continue to increase, unlike BC and SO${}_{\mathrm{2}}$. The largest methane sources in 2015 were agriculture (42 % of total emissions), oil and gas (extraction and distribution) (18 %), waste (18 %), and energy production (including coal mining) (16 %) . CO and NO${}_x$ emissions have been declining steadily and are expected to continue declining in the future.

In comparison to the CMIP6 emissions , ECLIPSEv6b emissions have additionally taken into account the recent declines in emissions from Asia of SO${}_{\mathrm{2}}$, BC, and NO${}_x$ due to recent control measures, whereas those declines in the CMIP6 emissions were unrealistically small . The inclusion of emissions from the flaring sector in Russia was a significant improvement, which was not present in the previous version of ECLIPSE emissions that was used in the report.

For non-agricultural fire emissions, many models utilized the CMIP6 fire emissions, which are based on monthly GFED (Global Fire Emissions Database) v4.1 . About half of the models included volcanic emissions or stratospheric aerosol concentrations from the CMIP6 dataset or other sources, and the other half did not include volcanic emissions, which mainly impact SO${}_{\mathrm{2}}$ and thus modeled SO${}_{\mathrm{4}}^{\mathrm{2}-}$. The emissions from the October to December 2014 Honoluraun volcano eruption were included by six models in a separate set of simulations. Similar differences in biogenic and agricultural waste emissions appear in these model simulations, and all are summarized in Table .

Table 2

Summary of emissions used in the models

Model	Biogenic	Volcanic	Forest fire	Agricultural waste
				burning
CanAM5-PAM	none	specified climatological emissions	CMIP6	ECLIPSEv6b
		and CMIP6 stratospheric aerosol
CESM2.0	MEGANv2.1	CMIP6	CMIP6	ECLIPSEv6b
CIESM-MAM7	none	CMIP6	CMIP6	ECLIPSEv6b
CMAM	none	none	CMIP6	ECLIPSEv6b
DEHM	MEGANv2	none	GFAS	ECLIPSEv6b
ECHAM6-SALSA	GEIA inventory	3-D emissions based	CMIP6	ECLIPSEv6b
	(PM only)	on AeroCom III
EMEP MSC-W	EMEP scheme	degassing from Ethna, Stromboli,	FINN	ECLIPSEv6b
		Eyjafjallajökull (2010), Grimsvotn (2011),
		Holuhraun (2014, 2015)
FLEXPART	none	none	CMIP6	ECLIPSEv6b
GEM-MACH	BEIS v3.09	none	CFFEPS	ECLIPSEv6b outside
				NA, US NEI and Canadian APEI
GEOS-Chem	MEGANv2.1 with updates	NASA/GMAO	GFEDv4.1	ECLIPSEv6b
GISS-E2.1	isoprene,	AeroCom	CMIP6	ECLIPSEv6b
	ORCHIDEE terpenes, online DMS,
	SS and dust
MATCH	MEGANv2	climatological and Honoluraun	CMIP6	ECLIPSEv6b
MATCH-SALSA-RCA4	MEGANv	climatological and Honoluraun	CMIP6	ECLIPSEv6b
MRI-ESM2		CMIP6 stratospheric aerosol and Honoluraun	CMIP6	ECLIPSEv6b
NorESM1-happi		CMIP6	CMIP6	ECLIPSEv6b
Oslo-CTM	MEGAN-MACC at 2010	AeroCom ()	GFEDv4	ECLIPSEv6b
UKESM1	isoprene and monoterpenes interactive	climatology, CMIP6	CMIP6	CMIP6
	with land surface vegetation scheme
WRF-Chem	MEGANv2.1	none	GFED	ECLIPSEv6b

2.2 Chemistry

This section contains a summary of models' chemistry schemes, and we refer the reader to Appendix and references therein for more details.

2.2.1 Methane

All participating models that provided CH${}_{\mathrm{4}}$ output prescribed CH${}_{\mathrm{4}}$ concentrations based on box model results from for 2015 and from for years prior to 2015. The former utilized the ECLIPSE v6B anthropogenic CH${}_{\mathrm{4}}$ emissions (Sect. ), along with assumptions for the natural emissions , to provide as input to models' surface or boundary layer CH${}_{\mathrm{4}}$ concentrations. Models then allow CH${}_{\mathrm{4}}$ to take part in photochemical processes, such as the production of tropospheric O${}_{\mathrm{3}}$.

2.2.2 Tropospheric chemistry

There is a wide range of tropospheric gas-phase chemistry implemented in the models. Air-quality-focused models, such as DEHM, EMEP MSC-W, GEM-MACH, GEOS-Chem, MATCH, and WRF-Chem, have detailed HO${}_x$–NO${}_x$–hydrocarbon O${}_{\mathrm{3}}$ chemistry, with speciated volatile organic compounds (VOCs) and secondary aerosol formation. The GISS-E2.1, MRI-ESM2, and UKESM1 ESMs also use this level of tropospheric chemistry. In contrast, climate-focused models like CanAM5-PAM, CIESM-MAM7, ECHAM-SALSA, and NorESM1 contain bare minimum gas-phase chemistry and use prescribed O${}_{\mathrm{3}}$ fields (e.g., CanAM5-PAM uses CMAM climatological O${}_{\mathrm{3}}$ fields). The CCMs are somewhere in between, with simplified tropospheric and stratospheric chemistry so that they could be run for longer time periods. For example, CMAM's tropospheric chemistry consists only of CH${}_{\mathrm{4}}$–NO${}_x$–O${}_{\mathrm{3}}$ chemistry, with no VOCs.

2.2.3 Stratospheric chemistry

Only a subset of the participating models have a fully simulated stratosphere. CMAM, MRI-ESM2, GISS-E2.1, OsloCTM, and UKESM1 contain a relatively complete description of the HO${}_x$, NO${}_x$, Cl${}_x$, and Br${}_x$ chemistry that controls stratospheric ozone along with the longer-lived source gases such as CH${}_{\mathrm{4}}$, N${}_{\mathrm{2}}$O, and chlorofluorocarbons (CFCs). Other models have a simplified stratosphere, such as GEOS-Chem, which has a linearized stratospheric chemistry scheme (Linoz, ), and WRF-Chem, which specifies stratospheric concentrations from climatologies – both of which do not simulate stratospheric chemistry. Finally, several models have no stratosphere or stratospheric chemistry at all (e.g., CIESM-MAM7, GEM-MACH, DEHM, and EMEP MSC-W).

2.2.4 Aerosols

Most models contain speciated aerosols: mineral dust (also known as crustal material), sea salt, BC, OA (sometimes separated into primary and secondary), SO${}_{\mathrm{4}}^{\mathrm{2}-}$, nitrate (NO${}_{\mathrm{3}}^{-}$), and ammonium (NH${}_{\mathrm{4}}^{+}$). However, some, like CanAM5-PAM and UKESM1, do not simulate NO${}_{\mathrm{3}}^{-}$ and NH${}_{\mathrm{4}}^{+}$, but assume all is in the form (NH${}_{\mathrm{4}}$)${}_{\mathrm{2}}$SO${}_{\mathrm{4}}^{\mathrm{2}-}$. OA, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, NO${}_{\mathrm{3}}^{-}$, and NH${}_{\mathrm{4}}^{+}$ are involved in chemical reactions interacting with the gas-phase chemistry. Aerosol size distributions are either prescribed or discretized into lognormal modes or size sections. How the aerosol size distribution varies in space and time depends on many different processes, including emission, aerosol microphysics, aerosol–cloud interactions, and removal. How these processes are parameterized depends on the model, and we refer the reader to the Appendix and the references therein for more detail.

3 Measurements

We have utilized many freely available observational datasets of SLCFs to evaluate the models. General descriptions are given below under the broad headings of surface monitoring, satellite, and campaign datasets, and there is some additional information in Appendix .

3.1 Surface monitoring datasets

3.1.1

CH${}_{\mathrm{4}}$ and O${}_{\mathrm{3}}$

Global surface CH${}_{\mathrm{4}}$ measurements were obtained from the World Data Centre for Greenhouse Gases (WDCGG). These measurements were made via gas chromatography, which has a $<$ 1 % uncertainty range. Surface in situ O${}_{\mathrm{3}}$ measurements are typically made via various types of UV absorption monitors, employing the Beer–Lambert law to relate UV absorption of O${}_{\mathrm{3}}$ at 254 nm directly to the concentration of O${}_{\mathrm{3}}$ in the sample air (e.g., ), which have approximately a 3 % or 1–2 ppbv uncertainty range. We obtained surface O${}_{\mathrm{3}}$ measurements from various networks: the National Air Pollutant Surveillance Program (NAPS) and the Canadian Pollutant Monitoring Network (CAPMON) for Canada, the Chemical Speciation Network (CSN) for the US, the Beijing Air Quality and Hong Kong Environmental Protection Agency for China, the Climate Monitoring and Diagnostics Laboratory (CMDL) for some global sites, the European Monitoring and Evaluation Programme (EMEP), and some individual Arctic monitoring stations like Villum Research Station and Zeppelin Mountain. Many of these measurements were downloaded from the EBAS database. The Arctic O${}_{\mathrm{3}}$ measurement locations are shown in Fig. .

Figure 1

Locations of Arctic surface in situ measurements, including (a) O${}_{\mathrm{3}}$, (b) BC in brown, and ice cores in black, as well as (c) SO${}_{\mathrm{4}}^{\mathrm{2}-}$.

[Figure omitted. See PDF]

3.1.2

CO, NO, and NO${}_{\mathrm{2}}$

CO and NO${}_x$ measurements were obtained from the same monitoring networks as O${}_{\mathrm{3}}$. CO instrumentation is similar to that for O${}_{\mathrm{3}}$; however, it uses gas filter correlation to relate infrared absorption of CO at 4.6 $\mathrm{\mu }$m to the concentration of CO in the sample air . For NO${}_x$, the instrument deploys the characteristic chemiluminescence produced by the reaction between NO and O${}_{\mathrm{3}}$, the intensity of which is proportional to the NO concentration. NO${}_{\mathrm{2}}$ measurements are approximated using its thermal reduction to NO by a heated (350${}^{\circ }$C) molybdenum converter . Note that this method has an estimated bias of about 5 %–20 % because of sensitivity to other oxidized nitrogen species, and this has not been corrected for. The bias is on the lower end for high-NO${}_x$ conditions and in the low-NO${}_x$ Arctic can be up to 100 % uncertainty.

3.1.3 BC and OA

There are various BC measurement methods exploiting different properties of BC and thus measuring different quantities : elemental carbon (EC) determined by thermal and/or thermal–optical methods, equivalent BC (eBC) by optical absorption methods, and refractory BC (rBC) by incandescence methods. Table  in Appendix  lists the different measurement techniques and instruments that the different monitoring networks and individual Arctic monitoring stations use. As BC emission inventories, including ECLIPSEv6b, are mainly based on emission factors derived from thermal and/or thermal–optical methods, modeled BC is consequently representative of EC.

The different types of BC measurements (EC, eBC, and rBC) usually agree with each other within a factor of 2 . However, it has been shown that, as the aerosol ages, the complex state of mixing of BC particles causes eBC to increase relative to EC . The absorption and scattering cross sections of coated BC particles vary by more than a factor of 2 due to different coating structures. found an increase of 20 %–250 % in absorption during aging, significantly depending on coating morphology and aging stages. Thus, this complexity impacts model–measurement comparisons at remote Arctic locations where one would expect eBC to have a high, positive uncertainty.

We obtained BC from the Canadian Aerosol Baseline Measurement (CABM) network for Canada, Interagency Monitoring of Protected Visual Environments (IMPROVE) network for the US, the EMEP network for Europe, and individual Arctic locations. To our knowledge, there were no other freely accessible BC measurements. The major observing networks EMEP, CABM, and IMPROVE measure EC with approximately 10 % uncertainty . However, given the complexities in different BC measurement types, as mentioned above, the overall uncertainty is about 200 %.

Another complexity with model evaluation of BC is that some of the eBC measurements that models are compared to were made from collected particulate matter with different maximum diameters (e.g., PM${}_{\mathrm{1}}$, PM${}_{\mathrm{2.5}}$, and PM${}_{\mathrm{10}}$). These are included in Table  for each of the measurement locations. From the models we use BC from PM${}_{\mathrm{2.5}}$, as most of the BC is expected to be in the submicron mode.

Organic carbon (OC) is also measured via thermal and/or thermal–optical methods using the same instrumentation as for EC detection in IMPROVE, CABM, NAPS, and EMEP measurement networks. These OC measurements have approximately 20 % or less uncertainty . Models output organic aerosol (OA), which includes OC and organic matter and is related to OC via a factor of 1.4 , though this factor has been reported as a range from 1.4 to 2.1 in the literature, depending on the source of OC and OA . Nevertheless, we applied a conversion factor of 1.4 to the OC measurements before comparing the modeled OA.

Arctic BC measurement locations are shown in Fig. , and many of these Arctic aerosol measurements were discussed in . We also evaluated modeled BC deposition by comparing it to BC deposition derived from ice core measurements (D4, ACT2: ; Humboldt: ; Summit: ; NGT_B19, ACT11D: ). All of the ice core locations are also shown in Fig. . Deposition fluxes are not a measured value but are derived from the EC concentrations in ice and precipitation estimates.

3.1.4

SO${}_{\mathrm{4}}^{\mathrm{2}-}$

Surface in situ SO${}_{\mathrm{4}}^{\mathrm{2}-}$ measurements in the major observing networks typically use ion chromatography methods, which have an approximately 3 % uncertainty range . However, SO${}_{\mathrm{4}}^{\mathrm{2}-}$ measurements have been shown to have up to 20 % analytical uncertainty . SO${}_{\mathrm{4}}^{\mathrm{2}-}$ datasets were obtained from IMPROVE, EMEP, and CABM networks, often via the EBAS database.

SO${}_{\mathrm{4}}^{\mathrm{2}-}$ deposition was also derived from the same ice core measurements mentioned above for BC deposition (D4, ACT2: ; Humboldt: ; Summit: ; NGT_B19, ACT11D: ). The Arctic SO${}_{\mathrm{4}}^{\mathrm{2}-}$ measurement locations are shown in Fig. .

3.1.5

PM${}_{\mathrm{2.5}}$

Surface in situ PM${}_{\mathrm{2.5}}$ measurements are usually made via gravimetric analysis of particulate matter collected on a filter (e.g., Teflon substrate), which has around a 1 %–6 % uncertainty range . These data were obtained from Beijing Air Quality and US Embassy data for China, NAPS for Canada , IMRPROVE for the US, and EMEP and/or EBAS for Europe.

3.2 Satellite datasets

Satellite observations are useful for evaluating models on larger horizontal spatial scales and for evaluating the three-dimensional atmosphere – not the surface concentrations. Observations from three satellite instruments were used to evaluate model trace gas distributions in the free troposphere and, when appropriate, the lower stratosphere. These were the Tropospheric Emission Spectrometer version 7 (TES; ), the Atmospheric Chemistry Experiment–Fourier Transform Spectrometer version 4.1 (ACE-FTS; ), and the Measurements of Pollution in the Troposphere version 8 (MOPITT; ). The vertical profiles of trace gas volume mixing ratios are derived or retrieved from the satellite-measured emission or absorption spectra, with varying degrees of vertical sensitivity. These remote techniques typically have about a 15 % uncertainty in the measurements (e.g., ), though this depends on the specific instrument and the species retrieved (e.g., ).

Note that while TES and MOPITT have global spatial coverage, their coverage does not extend up into the high Arctic. The TES instrument on NASA's Aura satellite measures vertical profiles of trace gases such as O${}_{\mathrm{3}}$, CH${}_{\mathrm{4}}$, NO${}_{\mathrm{2}}$, CO, and HNO${}_{\mathrm{3}}^{-}$ from 2004–present. After interpolating all models and TES results to a 1${}^{\circ }$ $\times$ 1${}^{\circ }$ horizontal grid, the monthly mean CH${}_{\mathrm{4}}$ and O${}_{\mathrm{3}}$ from the TES lite products were matched in space and time with models. Models were smoothed with the TES monthly mean averaging kernels prior to comparisons with satellite data. TES measurements started in 2004, stopped in late 2015, and had poorer coverage in the last few years.

A similar comparison method was used for MOPITT data. The MOPITT instrument on NASA's Terra satellite measures CO from 2000 to the present.

The ACE-FTS instrument on CSA's SCISAT satellite has measured the trace gases O${}_{\mathrm{3}}$, CO, NO, NO${}_{\mathrm{2}}$, and CH${}_{\mathrm{4}}$, among over 30 others from 2004–present. SCISAT has a high-inclination orbit, giving its instruments better coverage in the Arctic. ACE-FTS is a limb-sounding instrument measuring the solar absorption spectra of dozens of trace gas concentrations from the upper troposphere to the thermosphere. This gives us the opportunity to evaluate the 3-D model output in a region of the atmosphere where the radiative forcing of ozone is at its highest. Evaluating models with ACE-FTS measurements also provides insight into models' transport and upper-tropospheric chemistry. As was shown in , 3-hourly model output (rather than monthly mean output) is required for accurate comparisons to ACE-FTS data; thus, only models that provided output at this time frequency were compared to ACE-FTS measurements. The model output was sampled to match the times and locations of ACE-FTS measurements. We used an updated version of the advanced method in . Instead of taking the model output at the closest time to the ACE-FTS measurement time, the model output was linearly interpolated onto the ACE-FTS time. This reduces the bias introduced by diurnal cycles, which can cause certain volume mixing ratios (VMRs; e.g., that of NO and NO${}_{\mathrm{2}}$) to vary significantly between model output times. As in , the model output is also interpolated vertically in log pressure space and bilinearly in latitude and longitude to account for spatial variation between model grid points.

3.3 Measurement campaigns

Finally, there were air- and ship-based measurement campaigns of black carbon that were used for model evaluation. Aircraft campaigns allow vertical profiles of chemical species to be evaluated, and ship campaigns allow for in situ measurements in the remote Arctic seas.

3.3.1 Aircraft campaigns

The flight paths of the aircraft used for model evaluation of BC are shown in Fig. . The aircraft campaigns include A-FORCE , ARCPAC , ARCTAS , EUCAARI , HIPPO , NETCARE , and PAMARCMIP . Most of these aircraft campaigns occurred during boreal spring and summer months (April to July) except for one (HIPPO) occurring in January and November, and most occurred during the 2008–2009 time period, with only one (NETCARE) occurring during 2014–2015. All of these aircraft campaigns measured rBC from single-particle soot photometers (SP2) .

Figure 2

Flight tracks of BC aircraft campaigns used in this study.

[Figure omitted. See PDF]

The AMAP models that submitted 3-hourly BC output were linearly interpolated onto the aircraft locations in space and time using the Community Intercomparison Suite (CIS; ) in order to provide representative comparisons and robust evaluation.

3.3.2 Ship campaigns

There were three ship-based measurement campaigns in 2014–2015. These were the two Japanese campaigns (MR14-05 and MR15-03 cruises of R/V Mirai) in September of 2014 and 2015 (track from Japan to north of Alaska; ) and the Russian campaign in October 2015 (track north of Russia, from Arkhangelsk to Severnaya Zemlya and back; ) – both are shown in Sect.  (Fig. ). Models that provided 3-hourly BC output were compared to these observations. The Russian measurements of aerosol eBC concentrations were determined continuously using an Aethalometer purposely designed by MSU/CAO . Light attenuation caused by the particles depositing on a quartz fiber filter was measured, and the light attenuation coefficient of the collected aerosol was calculated. eBC concentrations were determined continuously by converting the time-resolved light attenuation to the eBC mass corresponding to the same attenuation and characterized by a specific mean mass attenuation coefficient, in calibration with AE33 (Magee Scientific).

The Japanese measurements provide rBC (refractory BC). showed that rBC and eBC are linearly related; thus, in order to compare the observations to models, we converted rBC to eBC via a factor of 1.8 (eBC $=$ 1.8 $\times$ rBC; ).

4 Model–measurement comparisons

In this section, we evaluate modeled SLCFs from the 18 participating models with a focus on performance in the Northern Hemisphere midlatitudes (defined for our purposes as 30–60${}^{\circ }$ N) and the Arctic (defined here as $>$ 60${}^{\circ }$ N for simplicity). Unless otherwise noted, the observations are compared to the model grid box that they are located in, and when more than one observation location occurs in the same model grid box, those observations are averaged first before the comparison. We look at spatial patterns in the model biases, as well as the vertical distribution and the seasonal cycles for each species, but first we start by providing a multi-species summary of the annual mean model biases in the surface air.

4.1 Multi-species summary

The 2014–2015 average modeled percent biases for surface concentrations of SLCFs are shown in Fig.  for each model and the multi-model mean (mmm). This figure is based on the model comparisons at the surface observation locations that will be shown in subsequent sections (Figs. , , , , , , , , and and additional American observations from the IMPROVE network for BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, and OA).

Figure 3

Mean 2014–2015 model percent biases for each model and the mmm for surface SLCF concentrations as well as BC and SO${}_{\mathrm{4}}^{\mathrm{2}-}$ deposition at (a) midlatitudes and (b) the Arctic. Note that the color scale is not linear.

[Figure omitted. See PDF]

Figure b shows that, for surface Arctic concentrations, no one model performs best for all species but that the mmm performs particularly well. It also shows that the model biases vary quite a bit among SLCF species for both the midlatitudes and the Arctic. It is important to note that there are many more measurement locations at midlatitudes compared to in the Arctic. BC, CH${}_{\mathrm{4}}$, O${}_{\mathrm{3}}$, and PM${}_{\mathrm{2.5}}$ have the smallest model biases out of the SLCFs of this study, whereas OA, CO, and NO${}_{\mathrm{2}}$ have larger model biases.

We find for half of the SLCF species that the mmm percent bias decreases with latitude (Fig. ). O${}_{\mathrm{3}}$, NO${}_{\mathrm{2}}$, BC, and SO${}_{\mathrm{4}}$ have a negative slope in the bias vs. latitude figure. So if the mmm bias was high at the midlatitudes, it is close to zero in the Arctic (O${}_{\mathrm{3}}$), and if the mmm bias was near zero at midlatitudes, it is negative in the Arctic (NO${}_{\mathrm{2}}$, BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$). This implies that there is not enough long-range transport from the midlatitude source regions to the Arctic. That said, the mmm CH${}_{\mathrm{4}}$ bias stays relatively constant with latitude, and we will see in Sect.  that this result is model-dependent. The CO, PM${}_{\mathrm{2.5}}$, and OA mmm biases have an increasing trend with latitude. However, both CO and PM${}_{\mathrm{2.5}}$ have no observation locations in the high Arctic, so those results cannot represent long-range transport. OA only has one observation location in the high Arctic, and its bias is very large overall, so issues other than long-range transport are at play, as we will see in the following discussion (Sect. ).

Figure 4

Multi-model mean percent biases for surface SLCF concentrations versus latitude (in 10${}^{\circ }$ latitude bins).

[Figure omitted. See PDF]

Of course, there are spatial, temporal, and model differences in the results, so we will now explore model performance for each SLCF in more detail in the next subsections.

4.2 Methane

Measured annual mean surface methane is shown in the top left panel of Fig. , along with model biases in the rest of the panels. Recall that unlike the rest of the SLCF species in this study, CH${}_{\mathrm{4}}$ concentrations were prescribed in these models from the same CH${}_{\mathrm{4}}$ dataset . That said, the different decisions by modelers on how those CH${}_{\mathrm{4}}$ global concentrations are distributed make differences in how these models compare to measurements. The mean model biases are small and mainly positive; in the midlatitudes, the multi-model mean bias is $+$145 ppbv (or $+$8.5 %), and in the Arctic, the mmm bias is 24 ppbv (or 1.3 %), which means that the models simulate the magnitude of surface CH${}_{\mathrm{4}}$ well – though still outside the $<$ 1 % measurement uncertainty range. There is a gradient in CH${}_{\mathrm{4}}$ VMRs (higher in the Northern Hemisphere and lower in the Southern Hemisphere) that is seen in the measurements (Fig. , top left) and reported in the literature (e.g., ), though it is not well captured by CMAM, MRI-ESM2, and UKESM1 models, which are all biased low in the Northern Hemisphere and biased higher towards the south. That is because of the simplifications made in these models' distributions of CH${}_{\mathrm{4}}$. For example, CMAM used a single global average CH${}_{\mathrm{4}}$ concentration that is interpolated linearly in time from once-yearly values.

Figure 5

Measured surface-level methane (top left) (ppbv, left color bar) and (remaining panels) model biases (model minus measurement, also ppbv, right color bar) for 2014–2015.

[Figure omitted. See PDF]

Figure also shows that observed annual mean surface CH${}_{\mathrm{4}}$ ranges geographically from about 1500 to 2100 ppbv depending on location; however, the models have a much smaller range due to their prescribing CH${}_{\mathrm{4}}$ concentrations as a lower boundary input. For example CMAM CH${}_{\mathrm{4}}$ volume mixing ratios only span about $\pm$3 ppbv around 1836 ppbv. The span of MRI-ESM2 surface CH${}_{\mathrm{4}}$ is even smaller. GEOS-Chem, GISS-E2.1, and OsloCTM have a more realistic range of 1700–2000 ppbv, though they still do not get the full variability that is seen in surface CH${}_{\mathrm{4}}$ mixing ratios close to CH${}_{\mathrm{4}}$ sources. However, in the free troposphere (above the boundary layer), we have TES satellite measurements of CH${}_{\mathrm{4}}$ that show that CH${}_{\mathrm{4}}$ is much more smoothly distributed aloft. Thus, the simplification of prescribing CH${}_{\mathrm{4}}$ concentrations in the models is more realistic there (Fig. , showing the 600 hPa level in the mid-troposphere). Additionally, Fig.  better illustrates the latitudinal gradient in CH${}_{\mathrm{4}}$ over the globe and its lack in some models, which have more negative biases in the Northern Hemisphere and more positive biases in the Southern Hemisphere. Other models, such as GISS-E2.1, do a good job of capturing the global distribution of CH${}_{\mathrm{4}}$.

Figure 6

TES measurements (top left) of CH${}_{\mathrm{4}}$ in the mid-troposphere (at 600 hPa, ppbv, left color bar) and (remaining panels) model biases for 2008–2009 (model minus measurement, also ppbv, right color bar). Results for 2014–2015 are similar but had less spatial coverage by the satellite. Gray areas have no data (either from the model, TES, or both).

[Figure omitted. See PDF]

In the Arctic, the vertical cross section of CH${}_{\mathrm{4}}$ VMRs over time as measured by the ACE-FTS in the middle to upper troposphere and in the stratosphere is shown in Fig. S1. There is a large decrease in CH${}_{\mathrm{4}}$ above the tropopause at around 300–100 hPa. The models are all biased low around 300 hPa and high around 100 hPa. This pattern is true for midlatitudes as well as in the Arctic and may imply that the altitude of the modeled tropopause is too low. This same conclusion was also found in via comparisons of these models' simulations to ozonesonde measurements and in our satellite O${}_{\mathrm{3}}$ comparison in the next section. The CH${}_{\mathrm{4}}$ model–measurement correlation coefficients for ACE-FTS are relatively high (e.g., $R=\mathrm{0.48}$ to 0.86 depending on the model).

Therefore, the general model evaluation for CH${}_{\mathrm{4}}$ indicates that because models do not explicitly model CH${}_{\mathrm{4}}$ emissions, they do not simulate the surface-level variability of CH${}_{\mathrm{4}}$ VMRs. Models differ in their global distribution of CH${}_{\mathrm{4}}$; thus, only some contain the north–south CH${}_{\mathrm{4}}$ gradient. Those that do not have the largest underestimations of Arctic tropospheric CH${}_{\mathrm{4}}$. The CH${}_{\mathrm{4}}$ evaluation also implies that the modeled tropopause height may be too low.

4.3 Ozone

Tropospheric O${}_{\mathrm{3}}$ is the third most important greenhouse gas , and it is a regional pollutant that causes damage to human health and ecosystems. O${}_{\mathrm{3}}$ is a secondary pollutant formed in the troposphere via photochemical oxidation of volatile organic compounds in the presence of nitrogen oxides (NO${}_x=$ NO $+$ NO${}_{\mathrm{2}}$). As such, models must simultaneously simulate the meteorological conditions, precursor species distributions, and photochemistry correctly in order to accurately simulate O${}_{\mathrm{3}}$. That said, since surface O${}_{\mathrm{3}}$ is an important contributor to poor air quality, there is significant pressure for models to simulate it accurately, particularly in the heavily populated midlatitudes (e.g., for air quality forecasting). Only models with prognostic O${}_{\mathrm{3}}$ are included in this section.

Figure shows the in situ summertime mean O${}_{\mathrm{3}}$ measurements (top left panel), and the model biases (remaining panels) and the same for the annual mean is shown in the Supplement (Fig. S2). These include averaging O${}_{\mathrm{3}}$ from hourly observations (day and night) and 3-hourly or monthly modeled O${}_{\mathrm{3}}$ depending on which were available for each model. Surface O${}_{\mathrm{3}}$ is overpredicted by most models, which has been documented previously . It has been shown that models can have problems producing low O${}_{\mathrm{3}}$ overnight . In the Arctic, simulated surface O${}_{\mathrm{3}}$ has more mixed results. Annual mean concentrations are of the order of 40 ppbv, and individual model biases range from $-\mathrm{20}$ % to $+$52 % globally on average for 2014–2015. The multi-model mean has a bias of $+$11 % for the Arctic, but this is not uniformly spatially distributed. All models overestimated surface O${}_{\mathrm{3}}$ in Alaska (mainly due to the overestimation of summertime concentrations, discussed below), and most models have too little O${}_{\mathrm{3}}$ at the Greenland location and in northern Europe.

Figure 7

Summertime (JJA) (top left) mean in situ surface O${}_{\mathrm{3}}$ measurements (ppbv, left color bar) and (remaining panels) model biases for 2014–2015 (model minus measurement, also ppbv, right color bar). Results for 2008–2009 are similar but are not available for China.

[Figure omitted. See PDF]

The models were all able to represent the summertime peak in the midlatitude seasonal cycle (not shown). In contrast to the more polluted midlatitudes, where surface O${}_{\mathrm{3}}$ peaks in the summertime due to photochemical production being at a maximum, Arctic O${}_{\mathrm{3}}$ is more influenced by the Brewer–Dobson circulation, bringing a maximum of tropospheric O${}_{\mathrm{3}}$ in the springtime due to photochemical production , descent from the stratosphere, and more long-range transport of O${}_{\mathrm{3}}$ to the Arctic. Figure , shows this springtime peak in both the western (a, longitude $<$ 0${}^{\circ }$) and eastern (b, longitude $>$ 0${}^{\circ }$) Arctic in the measurements. However, the models only capture that seasonal cycle in the eastern Arctic (Fig. b), implying that the models represent large-scale circulation and possibly stratosphere to troposphere exchange well. But it is interesting to note that the models that have sophisticated representation of stratosphere–troposphere exchange (such as CMAM, MRI-ESM2, UKESM1) do not particularly stand out as better performers in Fig.  compared to models that do not simulate the stratosphere (such as DEHM, MATCH, MATCH-SALSA). Thus, its impact on surface O${}_{\mathrm{3}}$ may be very small.

Figure 8

Surface O${}_{\mathrm{3}}$ monthly range that occurs at the locations in Fig.  above 60${}^{\circ }$ N. The measurements are the black and white boxes and whiskers, and the models are the colored box and whiskers. (a) The western Arctic and (b) the eastern Arctic for 2014–2015. Thick horizontal lines indicate the median O${}_{\mathrm{3}}$ VMR in each month, and the box extends to the interquartile range. The whiskers extend to the minimum and maximum monthly mean O${}_{\mathrm{3}}$ VMR.

[Figure omitted. See PDF]

In the western Arctic (Alaska mainly, Fig. a), models overestimate summertime Arctic O${}_{\mathrm{3}}$, likely due to overpredicting the impact of wildfire emissions on tropospheric O${}_{\mathrm{3}}$ concentrations, which is a research topic with high uncertainty . Another possibility is that modeled O${}_{\mathrm{3}}$ dry deposition over boreal vegetation is underestimated .

Some Arctic locations are more inclined to get springtime surface O${}_{\mathrm{3}}$ depletion due to bromine explosions and halogen chemistry . None of the model simulations in this study contain the necessary tropospheric halogen chemistry to simulate those events, which partly explains why some models in Fig.  (bottom) overestimate springtime O${}_{\mathrm{3}}$ concentrations. That particular feature is explored further on a site-to-site basis in .

The next subsection shows that both the O${}_{\mathrm{3}}$ precursors CO and NO${}_{\mathrm{2}}$ are underestimated compared to measurements at all global locations. This has implications for simulated tropospheric O${}_{\mathrm{3}}$ chemistry.

Free-tropospheric O${}_{\mathrm{3}}$ – satellite comparisons. Aircraft-based measurements and ozonesondes can provide insight into the vertical distribution of O${}_{\mathrm{3}}$, and these have been well documented (e.g., ). However, model grid boxes may not be representative of those fine-spatial-scale measurements. In this study, we examine how the model biases change in the vertical when compared to satellite measurements, which have a larger, “smoothed out” spatial sensitivity due to their viewing geometry and retrieval methods. Specifically, we compare modeled O${}_{\mathrm{3}}$ to TES and ACE-FTS satellite-based retrievals. These satellite instruments also have better global coverage than aircraft and sonde-based measurements.

The model fractional biases compared to TES measurements from near the surface up to 100 hPa are shown in Fig.  for the Arctic (left) and midlatitudes (right). All models' simulated fractional biases have similar vertical profiles for both the Arctic and midlatitudes, with greater negative values at lower levels and a more positive “bulge” of about 10 % around 300 hPa in the Arctic and about 5 % around 200 hPa at midlatitudes. That bulge in model biases at 300 hPa was also seen to a greater degree (50 %–70 %) when comparing these model simulations to Arctic ozonesonde measurements in . Compared to TES, which has much lower vertical resolution, the results are not as striking but are consistent with ACE-FTS measurements. On average, the models have a negative bias at all vertical levels in the Arctic region and in lower troposphere in the midlatitude region, whereas a positive bias is seen in the upper troposphere below 60${}^{\circ }$ N. This is consistent for the two time periods (2008–2009 and 2014–2015).

Figure 9

Vertical distribution of models' O${}_{\mathrm{3}}$ percent biases (model minus measurement over measurement) for 2008–2009 compared to the TES measurements; (a) average for midlatitudes (30–60${}^{\circ }$ N) and (b) average for Arctic latitudes ($>$ 60${}^{\circ }$ N).

[Figure omitted. See PDF]

Given that Fig.  shows positive biases near the midlatitudes, while Fig.  shows lower O${}_{\mathrm{3}}$ in the free troposphere, these results imply that there is not enough vertical lifting and/or mixing of tropospheric O${}_{\mathrm{3}}$ in most of the models. However, the TES measurements have been shown to be biased high by approximately 13 % throughout the troposphere , which is the same amount that the mmm is low. Similarly, ACE-FTS O${}_{\mathrm{3}}$ has an uncertainty range of $+$5 %–10 % when compared to O${}_{\mathrm{3}}$ from other satellite limb-view observations . The ACE-FTS comparison for O${}_{\mathrm{3}}$ can be found in the Supplement (Fig. S3), showing higher model biases around 300–100 hPa (except for GEOS-Chem) and good agreement below that. Therefore, overall, participating models simulate the free-tropospheric O${}_{\mathrm{3}}$ reasonably well and within the uncertainly limits of the observations.

Therefore, the general model evaluation for O${}_{\mathrm{3}}$ indicates that all models overestimate surface O${}_{\mathrm{3}}$ at midlatitudes, and that, combined with a lack of O${}_{\mathrm{3}}$ transport to the Arctic, results in modeled Arctic O${}_{\mathrm{3}}$ VMRs having relatively little bias (the right answer for the wrong reason). The summertime evaluation implies that models overestimate the O${}_{\mathrm{3}}$ produced and transported by wildfires in the western Arctic. The O${}_{\mathrm{3}}$ evaluation also implies that the modeled tropopause height may be too low.

4.4

O${}_{\mathrm{3}}$ precursors: carbon monoxide and nitrogen oxides

Figures  and show the comparisons of the multi-model medians (MMMs) to the surface in situ measurements. The figures for each model appear in the Supplement (Figs. S4 and S5), but only the MMMs are shown here since the spatial patterns were very similar for all models. The multi-model annual mean underpredicts both CO and NO${}_{\mathrm{2}}$ by approximately $-$55 % in the Northern Hemisphere for 2014–2015. The 2015 AMAP report showed similar findings for simulated surface CO, as have other studies , pointing to a likely underestimation of CO emissions and possibly shorter modeled lifetimes of CO due to an overestimation in OH . The annual mean surface CO underestimation is mainly dominated by the wintertime (e.g., the mmm bias in DJF is $-$92 %), when it has been reported that CO emissions from combustion are too low (e.g., ). All the models exhibit a large negative bias over China, which is consistent with the study by and is attributed to the enhanced destruction of CO by OH radicals, but it was also found in and that bottom-up CO emission inventories in Asia are greatly underestimated.

Figure 10

Mean CO volume mixing ratios (ppbv, left color bar) at surface measurement sites and MMM bias (MMM minus measurement ppbv, right color bar) for 2014–2015. Results from 2008–2009 are similar and not shown.

[Figure omitted. See PDF]

Figure 11

Mean NO${}_{\mathrm{2}}$ volume mixing ratios (ppbv, left color bar) at surface measurement sites and MMM bias (MMM minus measurement in ppbv, right color bar) for 2014–2015. Results from 2008–2009 are similar and not shown.

[Figure omitted. See PDF]

In the free troposphere, we compare modeled CO to that measured by MOPITT. Figure  shows these comparisons for the summertime (JJA) mean at the 600 hPa level, and Fig. S6 in the Supplement shows the same for the springtime (MAM) mean. Unlike in the winter and spring when all models are biased low at midlatitudes, there is more variability in the summertime biases, with about half the models overestimating free-tropospheric CO and the other half underestimating it. In the Arctic region, all models are biased high in the summer but low in the springtime. In the winter, MOPITT does not have coverage in the Arctic. discussed the fact that models had high biases in the tropospheric column of CO compared to MOPITT measurements in the outflow from Asia and low biases north of there due to lack of transport. The study also suggested that summertime CO transport out of Asia is too zonal when comparing models to IASI CO columns. At 600 hPa, where CO concentrations are lower and the atmosphere is more mixed, these features do not appear.

Figure 12

Summertime (JJA) (top left) mean MOPITT CO at 600 hPa (ppbv, left color bar) and (remaining panels) model biases (model minus measurement ppbv, right color bar) for 2014–2015.

[Figure omitted. See PDF]

In the upper troposphere and stratosphere, modeled CO and NO${}_x$ monthly time series are compared to measurements from the ACE-FTS satellite instrument (where NO${}_x=$ NO $+$ NO${}_{\mathrm{2}}$, which are measured separately), and those results are shown as Taylor diagrams in Fig. S7, along with O${}_{\mathrm{3}}$ and CH${}_{\mathrm{4}}$ at 150 hPa, which is in the upper troposphere–lower stratosphere (UTLS) region. The contours show the model's overall skill as defined in . Only the models that simulate the stratosphere were included, and the results show that there is a range in model performance by SLCF species, with no one model performing best for all. Comparison statistics for CO were poorer than those for O${}_{\mathrm{3}}$, CH${}_{\mathrm{4}}$, and NO${}_x$.

4.5 Black carbon

In this section, we examine the spatial and seasonal distributions of BC using ground-based measurements, which are primarily available in North America, Europe, and several locations in the Arctic, but are also available from two ship-based campaigns and several aircraft campaigns. Given the limited global data available for both BC and SO${}_{\mathrm{4}}^{\mathrm{2}-}$ (e.g., we could find none freely available for Asia), we focus the plots on the Arctic region here, and given that the magnitude of BC and SO${}_{\mathrm{4}}^{\mathrm{2}-}$ does not span a wide range throughout the Arctic, we show model biases as percent rather than absolute differences as was done in previous sections for trace gas species shown globally. We also analyze the BC model–measurement comparisons, keeping in mind that because there are various definitions and measurement types for BC, we consider agreement within a factor of 2 to be within the uncertainty range (Sect. ).

Figure  (top left panel) shows annual mean surface-level concentrations of black carbon (BC) at nine Arctic observation stations and (remaining panels) the model percent biases there. The annual mean BC concentrations are of the order of less than 1 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$, and most models tend to underestimate BC in the high Arctic while overestimating in Alaska and Scandinavia. This result could be partially explained by the discrepancy caused by the use of EC and eBC data, which are not the same (Sect. ). As aerosols age during transport to high-Arctic locations, their eBC (based on absorption converted to mass) gets more and more of a positive bias compared to EC. As models are more representative of EC, they will not be able to agree with eBC measurements in aged air at high-Arctic remote stations, such as Gruvebadet, Zeppelin, Alert, and Utqiaġvik. This is in contrast to the Alaskan and European stations, which are closer to sources where BC is more fresh; thus, the eBC measurements have lower uncertainty.

That said, a few models (CanAM5-PAM, DEHM, and FLEXPART) overestimate BC concentrations in the high Arctic. Overall individual model biases range $\pm$100 % at individual sites.

Figure 13

Mean BC concentrations ($\mathrm{\mu }$g m${}^{-\mathrm{3}}$, top color bar) at surface Arctic measurement sites and model bias (as (model–measurement) $/$ measurement in percent, bottom color bar) for 2014–2015. Results from 2008–2009 are similar and not shown.

[Figure omitted. See PDF]

The underestimation of high-Arctic atmospheric BC concentrations may be related to excessive BC deposition further south; however, there are very few BC deposition measurements. In the Arctic, we can evaluate total (wet $+$ dry) modeled deposition via derived ice core measurements. There were six ice cores on Greenland and one in the European Arctic in Spitsbergen (Lomosovfonna). Figure  shows that all models overestimate BC deposition fluxes at the ice core locations. While the ice cores contain BC data starting in the year 1750, only data after 1990 have been used to match the modeled time period (1990–2015, 1995–2015, 2008–2009, and 2014–2015, depending on the model).

Figure 14

Annual average BC deposition flux values for the seven ice core locations (Fig. b) for each model based on values from 2008–2009 and 2014–2015. The observed fluxes are plotted in black, and a black line indicates the level of the average observed flux; the black dashed line is the model mean for each location. The period used for plotting is based on all available years after 1990; the title indicates the last available year form the ice core record.

[Figure omitted. See PDF]

The measured BC deposition flux values on Greenland vary with elevation (lower fluxes at higher elevation). Summit (3177 m a.s.l.) has an average of 285 $\mathrm{\mu }$g m${}^{-\mathrm{2}}$ yr${}^{-\mathrm{1}}$ in contrast to ACT2 (2461 m a.s.l.) with 676 $\mathrm{\mu }$g m${}^{-\mathrm{2}}$ yr${}^{-\mathrm{1}}$. BC deposition is highest in the European Arctic at Spitsbergen with 856 $\mathrm{\mu }$g m${}^{-\mathrm{2}}$ yr${}^{-\mathrm{1}}$. For all seven ice cores used in this comparison the averaged model mean is 3 times as high as the observations. At D4 (2728 m a.s.l.) the modeled mean corresponds best to the observation, with a mean bias of $+$83 %. At ACT11 (2296 m a.s.l.) the models have 4 times the deposition flux compared to the measurements. Generally though, the model mean is skewed higher by FLEXPART and DEHM (Fig. ), which also had higher atmospheric BC concentrations. A few models simulated less BC deposition than observed at these sites, and these models also underestimated BC atmospheric concentrations. Thus, it is difficult to conclude that deposition biases are a cause for atmospheric biases when the two are inter-related parameters.

The seasonal cycles of surface-level atmospheric BC concentrations at several Arctic locations are shown in Fig. . As seen in , some models still underestimate wintertime BC, but many models now show similar seasonality as the observations. The multi-model mean also captures the monthly variations well including the summertime peak at some Alaskan sites caused by fire emissions. The multi-model mean Arctic BC is underestimated in winter ($-$24 %) and overestimated in the summer ($+$32 %), though overall, this is an improvement in model performance in simulating Arctic BC since the 2015 AMAP assessment of black carbon and ozone as climate forcers – the latter of which had $-$59 % winter bias and $+$88 % summertime bias . However, it is difficult to make direct comparisons to that report as those values were for a smaller set of Arctic locations, different observation periods, and with a different set of models (though many overlapping). The model improvement may be due to the improved anthropogenic emissions of BC, particularly from northern Russia, where flaring emission factors were increased in ECLIPSEv6B compared to those used for the 2015 AMAP report.

Figure 15

Modeled (thin colored lines) and measured (thick black line) monthly mean BC concentrations (in $\mathrm{\mu }$g m${}^{-\mathrm{3}}$) at surface Arctic measurement sites in 2014–2015. The multi-model mean is shown by the black dashed line.

[Figure omitted. See PDF]

Most models have reasonable spatial correlation with the measurements across the Arctic in that they correctly simulate the range of BC concentrations that appear across the Arctic (e.g., higher concentrations at Hurdal, lower concentrations at Zeppelin), as shown in Fig. S8. However, there are still large differences and low $R$ values in the statistics shown in Fig.  S8.

There are positive model biases at midlatitudes (in North America and Europe; not shown) for surface-level BC. The vertical analysis of BC from the aircraft campaigns (below) provides further insight and support for the suggestion that models do not have adequate long-range transport of the pollutants from their sources in the midlatitude; thus, they do not simulate enough pollution in the Arctic.

Vertical profiles of BC – aircraft campaigns. Gridded BC output at 3-hourly intervals was provided by 11 of the participating models and was compared to aircraft campaign measurements of BC. The interpolation of model output to flight track coordinates was carried out by tools from the Community Intercomparison Suite (CIS; ), which co-located the extracted model tracks with their corresponding observational values.

Figure shows the median vertical profiles of BC concentrations from the aircraft measurements and from the models. At midlatitudes, from 0–2 km, all of the models agree well with the measurements. However, BC concentrations decline steeply in a few models (e.g., MATCH-SALSA, EMEP-MSC-W, and GEOS-Chem) above 2 km. It would appear that they do not have enough vertical lifting of BC and/or perhaps too short of a BC lifetime. Indeed, one of these is EMEP MSC-W, in which the short BC lifetime was previously reported in . That said, in , the Oslo-CTM and ECHAM models were shown to overestimate the BC lifetime. In our case, OsloCTM is not shown in Fig.  because it did no provide BC at 3-hourly timescales. But ECHAM-SALSA results are consistent with the study in that they particularly overestimate BC in the upper altitudes in both the midlatitudes and the Arctic, implying too long a lifetime and too much long-range transport into the upper Arctic atmosphere. The measured BC profile at midlatitudes drops off more quickly around the tropopause at 11–12 km, and, except for CanAM5-PAM, the models do not reproduce this drop.

Figure 16

Median vertical profiles of observed (heavy black line) and modeled (colored lines) BC concentrations for all aircraft campaigns combined, separated into (left) midlatitudes and (right) the Arctic. The multi-model median is shown by the dashed black line.

[Figure omitted. See PDF]

In the Arctic profiles (Fig. ), the modeled and observed profiles do not decline with altitude throughout the troposphere, but the observed median BC concentration does drop sharply around 9 km – again near the Arctic tropopause, and again, the only model to capture that change is CanAM5-PAM. In the Arctic comparisons, the models that did not simulate enough BC aloft at the midlatitudes stand out as having larger underestimates of BC in the Arctic. For example, MATCH-SALSA and EMEP MSC-W have very low BC throughout the Arctic vertical profile. These results are consistent with the surface BC underestimation in Fig. . Therefore, the underestimation seen in the Arctic for those two models is likely due to a lack of long-range transport from the midlatitudes, as well as errors in BC deposition mentioned above. In addition, the study showed that different parts of the Arctic BC vertical profile are sensitive to BC transported from different areas of the world. For example, the lower-tropospheric BC is influenced by emissions transported from North America, Russia–Belarus–Ukraine, Europe, and East Asia, whereas upper-tropospheric Arctic BC is mainly influenced by transport from South Asia. Thus, the differences in the model results could be related to differences in how they simulate these transport pathways.

found that, overall, considerable differences in wet deposition efficiencies in the models exist and are a leading cause of differences in simulated BC burdens. Results from their model sensitivity experiments indicated that convective scavenging outside the Arctic reduces the mean altitude of BC residing in the Arctic, making it more susceptible to scavenging by stratiform (layer) clouds in the Arctic. Consequently, scavenging of BC in convective clouds outside the Arctic acts to substantially increase the overall efficiency of BC wet deposition in the Arctic, which leads to low BC burdens compared to simulations without convective BC scavenging. also found that convective scavenging at middle and subtropical latitudes removes a significant fraction of BC. In contrast, BC concentrations in the upper troposphere are only weakly influenced by wet deposition in stratiform clouds, whereas lower-tropospheric concentrations are highly sensitive – these are consistent with the results we find in this study, wherein the multi-model median is too high above about 9 km and too low from 0–9 km. Indeed, the MATCH and MATCH-SALSA models, for example, assume reduced scavenging of aerosol in mixed-phase clouds following , which increases long-range transport to the Arctic. It is odd then that MATCH is one of the better=performing models in Fig.  and MATCH-SALSA is not. Despite the large range in modeled vertical BC concentrations, the multi-model median is close to the observed throughout the troposphere at both midlatitudes and in the Arctic.

Arctic seas analysis – ship campaigns. From the ship-based measurements, we see that there is a consistent model overestimate of BC in the Pacific region (Japan cruise) where measured concentrations are very low (Fig. ). Indeed, report that BC concentrations were in the range 0–66 ng m${}^{-\mathrm{3}}$, with an overall mean value of just 1.0 $\pm$ 1.2 ng m${}^{-\mathrm{3}}$. The models, possibly due to their coarse resolutions, were not able to simulate such low background BC concentrations. However, even the model with the highest resolution (GEM-MACH at 15 km resolution) overestimated BC in the Pacific – though that limited-area model in that region near the boundary would have been heavily influenced by the upwind, coarser-resolution boundary conditions that were assumed (1${}^{\circ }$ $\times$ 1${}^{\circ }$ MOZART4 chemical boundary conditions). The high bias in the Pacific may be due to all models overestimating the amount of BC that gets transported off the Asian continent. That model bias is consistent with our low-altitude comparisons of the models to the HIPPO aircraft campaign measurements, which were taken over the northwestern Pacific (Fig. ). The BC overestimate over the Pacific was also found in the study looking at simulated BC from the AeroCom global model intercomparison initiative compared to HIPPO measurements.

Conversely, the modeled BC concentrations generally agree with measurements in the Russian Arctic Ocean, though they are biased slightly low for the most part. attribute the higher BC concentrations measured near the Kara Straight (north of 70${}^{\circ }$ N) to gas-flaring emissions, and near Arkhangelsk (White Sea), important sources were midlatitude biomass burning, transportation, and combustion (residential and commercial). Since models were able to simulate this well, their improvement is likely due to improved Russian anthropogenic emissions in ECLIPSE v6B (Sect. , ) compared to previous emissions datasets, which did not include enough Russian flaring emissions. The best model results were from ECHAM-SALSA and MATCH when compared to all of the ship campaign data.

Figure 17

Observed BC concentrations (top) ($\mathrm{\mu }$g m${}^{-\mathrm{3}}$, top color bar) along the ship paths and (remaining panels) the modeled vs. measured 3 h average BC concentrations along the ship paths. Note the logarithmic scale.

[Figure omitted. See PDF]

Therefore, the general model evaluation for BC indicates that while there is a large variability in models results, they tend to overestimate surface BC at midlatitudes (including over the Pacific Ocean) and underestimate surface BC in the Arctic. Again, these results point to a lack of transport to the Arctic and, in this case, too much BC deposition along the way. While we were only able to evaluate BC deposition in the Arctic in this study, those results support the hypothesis of some models having too much BC deposition. The BC vertical profile evaluation also implies that the modeled tropopause height may be too low.

4.6 Sulfate

We used monthly mean surface level observations of SO${}_{\mathrm{4}}^{\mathrm{2}-}$ from 18 Arctic sites to evaluate the models. Figure  shows that, similar to BC, the SO${}_{\mathrm{4}}^{\mathrm{2}-}$ concentrations in the high Arctic are underestimated by most of the models. A few models overestimate SO${}_{\mathrm{4}}^{\mathrm{2}-}$ in Scandinavia and Alaska.

Figure 18

Mean measured SO${}_{\mathrm{4}}^{\mathrm{2}-}$ concentrations ($\mathrm{\mu }$g m${}^{-\mathrm{3}}$, top color bar) at surface Arctic measurement sites and model bias (as model–measurement $/$ measurement in percent, bottom color bar) for 2014–2015. Results from 2008–2009 are similar and not shown.

[Figure omitted. See PDF]

The model underestimations of SO${}_{\mathrm{4}}^{\mathrm{2}-}$ could be mainly due to higher efficiencies of models in removing aerosol during the long-range transport to the high Arctic. This is consistent with a previous study based on AMAP 2015 model simulations that found that the convective wet deposition outside the Arctic region may have led to different seasonal cycles of aerosol concentrations in the Arctic . Dimethylsulfide (DMS), a naturally occurring source of sulfur from marine algae emissions, could also be misrepresented in models. However, this source would be more pronounced in the summer when there is less sea ice in the high Arctic, and it does not appear as though models are underestimating only in the summertime (Fig. ). Rather, some models appear too low in the winter and spring – which points towards underestimating local Arctic sources and to a lack of transport from midlatitudes as being the key issues. Despite the individual model differences in representing the seasonal cycle, the multi-model mean compares well with observations at most locations. However, the multi-model mean SO${}_{\mathrm{4}}^{\mathrm{2}-}$ is significantly underestimated at Alert and Irafoss sites. Mean model biases for all Arctic sites range from $-$65 % to $+$80 % among different models, and correlation coefficients are typically around 0.5 (Fig. S9).

Figure 19

Modeled (thin colored lines) and measured (thick black line) monthly mean SO${}_{\mathrm{4}}^{\mathrm{2}-}$ concentrations ($\mathrm{\mu }$g m${}^{-\mathrm{3}}$) at surface Arctic measurement sites in 2014–2015. The multi-model mean is shown by the black dashed line.

[Figure omitted. See PDF]

The high GEOS-Chem bias in the summertime seen in Fig.  was first reported in and found to be due to problems with cloud pH and cloud liquid water in the summer over the Arctic. In , the summertime Arctic surface SO${}_{\mathrm{4}}^{\mathrm{2}-}$ concentrations are reduced by a factor of 2 by reducing the cloud liquid water content to a uniform value of 1 $\times$ 10${}^{-\mathrm{7}}$ g m${}^{-\mathrm{3}}$ north of 65${}^{\circ }$ N in the model. The version of GEOS-Chem in our study is more recent and uses the offline cloud liquid water content from both GEOS-FP and MERRA2. We did not scale this variable down, which may be a reason for the high GEOS-Chem sulfate bias in Fig. .

The seasonal cycle for observations grouped together is shown in Fig. S10, showing a consistent seasonal cycle for 2008–2009 as seen in the observations. Most models (e.g., CanAM5-PAM, DEHM, MATCH, OsloCTM) are able to capture the seasonal cycle well. However, several models (e.g., CESM, CIESM-MAM7, ECHAM-SALSA, and EMEP-MSCW) strongly underestimate observed springtime peak values. Conversely, the models and the measurements showed a weaker seasonal cycle during the 2014–2015 time period (Fig. S11). It may be partly due to the local pollution sources in the Arctic during wintertime (e.g., Fairbanks). Those highly localized pollution events, caused by local emissions getting trapped in a stable boundary layer, occur on scales that are smaller than the model resolutions employed here can represent. Many models are also missing chemical formation processes for SO${}_{\mathrm{4}}^{\mathrm{2}-}$ in the absence of sunlight, which may explain underestimations seen in winter (e.g., ). An evaluation of SO${}_{\mathrm{2}}$ could help with our understanding but was beyond the scope of this study. A lack of dark chemistry may be true for organic aerosol, as discussed in the next section as well.

From October to December 2014, the Honoluraun volcanic eruption may have elevated SO${}_{\mathrm{4}}^{\mathrm{2}-}$ concentrations at some locations in the Arctic. However, in our model–measurement comparisons, there does not appear to be a large underestimate during those months, which implies that model performance was not impeded by not including those volcanic emissions.

As mentioned above, the uncertainty in wet deposition could be a significant factor in atmospheric SO${}_{\mathrm{4}}^{\mathrm{2}-}$ model biases. Previous studies have shown that models have too much washout in winter and not enough wet deposition in summer, leading to a “flatter” seasonal cycle than observed (e.g., Fig. ; ). As with BC in the previous section, the SO${}_{\mathrm{4}}^{\mathrm{2}-}$ deposition was evaluated here in the same manner. The average measured SO${}_{\mathrm{4}}^{\mathrm{2}-}$ deposition fluxes from ice cores for all locations (only Greenland was available here) are 18 mg m${}^{-\mathrm{2}}$ yr${}^{-\mathrm{1}}$. The lowest observed fluxes are found at D4 (12 mg m${}^{-\mathrm{2}}$ yr${}^{-\mathrm{1}}$) and the highest at ACT11D (30 mg m${}^{-\mathrm{2}}$ yr${}^{-\mathrm{1}}$). The model average for all locations is overestimated by around 20 % compared to measured fluxes. This is similar to the model biases in atmospheric SO${}_{\mathrm{4}}^{\mathrm{2}-}$ concentrations.

Figure 20

SO${}_{\mathrm{4}}^{\mathrm{2}-}$ deposition fluxes for the Greenland ice core locations shown in Fig. b. The observed fluxes are plotted in black, and a black line indicates the level of the average observed flux; the black dashed line is the multi-model mean at that location. The period used for plotting is based on all available years after 1990; the title indicates the last available year from the ice core record.

[Figure omitted. See PDF]

Therefore, the general model evaluation for SO${}_{\mathrm{4}}^{\mathrm{2}-}$ indicates that while there is a large variability in models results, as with BC, models underestimate surface SO${}_{\mathrm{4}}^{\mathrm{2}-}$ in the Arctic. The evaluation of SO${}_{\mathrm{4}}^{\mathrm{2}-}$ deposition in the Arctic is similar to BC, with both overestimating deposition.

4.7 Organic aerosol

Unfortunately, there is only one high-Arctic station with OA measurements (Alert, NV, Canada); however, there are a few additional stations measuring OA in the sub-Arctic (still $>$ 60${}^{\circ }$ N). These are all shown in Fig. . The seasonal cycles are shown in Fig. , and the model vs. measurement scatter plot along with some comparison statistics are presented in Fig. S12.

Figure 21

Mean OA concentrations ($\mathrm{\mu }$g m${}^{-\mathrm{3}}$, top color bar) at surface Arctic measurement sites and model bias (as model–measurement $/$ measurement in percent, bottom color bar) for 2014–2015. Results from 2008–2009 are similar and not shown.

[Figure omitted. See PDF]

Model biases have a large range of $\pm$200 % at the different locations, but the multi-model mean for the region is $+$65 %. At midlatitudes (30–60${}^{\circ }$ N), measurements are conducted mainly in the US, where the multi-model mean bias is $+$83 % for the 2014–2015 average (not shown).

Figure 22

Monthly Arctic OA from models (red) and measurements (black) for 2014–2015.

[Figure omitted. See PDF]

Several models (CanAM5-PAM, DEHM, CIESM-MAM7, ECHAM-SALSA, GEOS-Chem, MRI-ESM2, NorESM, OsloCTM, and UKESM1) are able to simulate the summertime peak in Arctic OA concentrations; however, the other seven models in Fig.  simulate a seasonal cycle that is too flat or peaks at the wrong time (e.g., the CESM seasonal cycle peaks too early in the year).

Figures and S11 both show that most models consistently overestimate Alaskan OA and underestimate European OA, consistent with our assessment of other species showing that models are likely overestimating wildfire influence in summertime Alaska. All models, except EMEP-MSC-W, used the GFED fire emissions inventory. EMEP-MSC-W used the FINN fire emissions inventory, which for BC$+$OC has been shown to be significantly lower than the GFED emissions . As a result, EMEP-MSC-W model biases in Alaska are lower. However, they are not the lowest. MATCH-SALSA has the lowest OA model biases in Alaska, despite that model using GFED emissions. The comparison statistics in Fig. S12 show highly varying results.

4.8 Fine particulate matter

PM${}_{\mathrm{2.5}}$ is partly connected to direct and indirect climate effects via its interactions with clouds. It is mainly composed of BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, NO${}_{\mathrm{3}}^{-}$, and NH${}_{\mathrm{4}}^{+}$, as well as crustal material (dust), sea salt, and water, though the water component is often dried off during measurements. Model biases of those species will contribute to the total PM${}_{\mathrm{2.5}}$ biases.

While BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, and OA were discussed above, it is beyond the scope of this project to evaluate the other major PM species, which, aside from water, have a smaller radiative impact. Note that the analysis in this section is focused on sub-Arctic and midlatitude sites closer to human populations rather than remote high-Arctic sites due to a lack of data. PM${}_{\mathrm{2.5}}$ is not a typical parameter included in the longer-term remote Arctic observations. Since PM${}_{\mathrm{2.5}}$ has important health impacts, it is well measured at air quality monitoring networks.

Figure 23

Measured ground-level PM${}_{\mathrm{2.5}}$ concentrations ($\mathrm{\mu }$g m${}^{-\mathrm{3}}$) and model biases (as model–measurement $/$ measurement in percent, bottom color bar) for 2014–2015. The upper color bar represents observations, and the lower bar represents model biases.

[Figure omitted. See PDF]

The model PM${}_{\mathrm{2.5}}$ biases at several locations in the United States, Europe, and China are within the 60 %–80 % range. However, some models (CanAM5-PAM, CIESM-MAM7, GEOS-CHEM, GEM-MACH, and Oslo-CTM) show biases larger than 200 %, especially in the western US and Alaska. The large inter-model differences in Fig.  are likely due to uncertainty in mineral dust. CanAM5-PAM has a particularly large contribution of dust to PM${}_{\mathrm{2.5}}$. The study showed that dust regions were globally one of the largest areas of diversity in PM${}_{\mathrm{2.5}}$ concentrations between different CMIP6 models. The EMEP MSC-W results are consistent with EMEP annual evaluations for Europe, where the model underestimates PM${}_{\mathrm{2.5}}$ by 10 %–25 %, including a few Arctic sites in Norway and Finland (, and annual model evaluation reports at https://www.emep.int/publ/common_publications.html, last access: 14 April 2022).

The simulated surface-level SLCFs were quite sensitive to the different meteorological conditions such as boundary layer stability and levels of photochemistry, which differed between the two time periods chosen in this study. 2014–2015 was also a time period with more wildfires compared to 2008–2009. For example, according to the CMIP6 emission data that were used in most of the models, emissions of BC from Canadian wildfires in 2014–2015 were 340 % higher than in 2008–2009, whereas the emissions from the USA and Russia were similar for these years. Given the very intense wildfire emissions and low anthropogenic emissions in northern Canada in 2014–2015, differences in simulated PM${}_{\mathrm{2.5}}$ concentrations over Canada and Alaska can be partly attributed to differences in simulations of wildfire aerosols in the models.

Some models (CanAM5-PAM, CESM, CIESM-MAM7, GEOS-Chem, and WRF-Chem) simulate higher PM${}_{\mathrm{2.5}}$ and more variable PM${}_{\mathrm{2.5}}$ in the summertime (e.g., Fig. S13 in the Supplement). While this is seen to some extent in the observations, this may be due to the way fire emissions and sea salt emissions are treated in these models. Fire emissions, fire plume injection height, plume rise, and wet deposition of fire pollutants are all highly uncertain model processes and a subject of ongoing research (e.g., ). Indeed, the individual model PM${}_{\mathrm{2.5}}$ Arctic biases are more tightly clustered for 2008–2009 when there were fewer fires. showed that Arctic PM${}_{\mathrm{2.5}}$ seasonal pollution is mainly due to local air pollution in the winter and due to fires in the summer.

Figure 24

Annual mean PM${}_{\mathrm{2.5}}$ comparisons between station observations and model simulations for the year 2015.

[Figure omitted. See PDF]

Figure  shows that the annual mean simulated PM${}_{\mathrm{2.5}}$ concentrations compare well with observations and the correlation coefficients are relatively high ($R=\mathrm{0.8}$ or higher for all models). The high concentrations in China and low concentrations in the US and Europe are captured by the models, providing confidence for health impact assessments that utilize these model results.

5 Conclusions

In this study, we evaluated the SLCF simulation capabilities of 18 models that were used in the 2022 AMAP SLCF assessment report. Our conclusions are grouped into the questions we aimed to answer in the Introduction.

5.1 How well do the AMAP SLCF models perform in the context of measurements and their associated uncertainty?

Recall that the in situ SLCF measurements had the following reported uncertainties: CH${}_{\mathrm{4}}$ 1 %, O${}_{\mathrm{3}}$ 3 %, CO 5 %, NO${}_{\mathrm{2}}$ 5 %–100 %, BC 200 %, SO${}_{\mathrm{4}}^{\mathrm{2}-}$ 20 %, OA 20 %, and PM${}_{\mathrm{2.5}}$ 1 %–6 %. However, since the variability in measurements from different techniques was only really taken into account for the BC uncertainty, and since we are comparing annual mean results to each other, it is not a fair comparison to say that models and measurements agree with each other if model biases are within the reported measurement uncertainty range. However, we do use those numbers as a rough guideline for “good” model performance in the absence of other quantitative criteria.

Some model annual mean biases were within those uncertainty ranges. For example, CMAM, MRI-ESM2, and UKESM1 simulate Arctic CH${}_{\mathrm{4}}$ to within 1 %, thus agreeing with the CH${}_{\mathrm{4}}$ measurements. However, at midlatitudes, they are all out of range at around $+$6 %–10 %. MATCH and WRF-Chem simulated midlatitude O${}_{\mathrm{3}}$ to within 2 %, but only MATCH-SALSA was within 3 % in the Arctic region. The Arctic NO${}_{\mathrm{2}}$ measurements are highly uncertain at around 100 %, so all of the models agreed with Arctic NO${}_{\mathrm{2}}$ measurements. However, in the higher-NO${}_{\mathrm{2}}$ midlatitude environment, NO${}_{\mathrm{2}}$ measurement uncertainty is at the smaller end of the range. OsloCTM and WRF-Chem midlatitude NO${}_{\mathrm{2}}$ biases were within 10 %. All models agree with BC measurements in both midlatitudes and the Arctic, as all biases are less than 80 %. CESM, CIESM-MAM7, DEHM, MATCH, UKESM1, and WRF-Chem all simulate midlatitude SO${}_{\mathrm{4}}^{\mathrm{2}-}$ to within 20 %. But only CanAM5-PAM and MRI-ESM2 do the same for the Arctic region. OA had some of the largest model biases (Figs.  and ), though ECHAM-SALSA, EMEP-MSC-W, GISS-E2.1, MATCH, and OsloCTM are all within 20 % in the Arctic but none at midlatitudes. Finally, with the small uncertainty in PM${}_{\mathrm{2.5}}$, only CIESM-MAM7 in the midlatitudes and GISS-E2.1 in the Arctic agree within 2 %.

To summarize the mmm annual mean performance, it “matches” surface observations in the Arctic for CH${}_{\mathrm{4}}$, NO${}_{\mathrm{2}}$, BC, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, OA, and PM${}_{\mathrm{2.5}}$ – as such, the mmm has the best overall performance for the Arctic. In the midlatitudes the mmm matches surface observations for BC and SO${}_{\mathrm{4}}^{\mathrm{2}-}$ only.

Regarding the comparisons of trace gases to the TES, MOPITT, and ACE-FTS satellite measurements (which have roughly 5 %–20 % uncertainty), models agree well. Free-tropospheric distributions of trace gases are somewhat easier to simulate, as common problems like a boundary layer that is too stable or too much deposition do not negatively impact the free-tropospheric SLCF distributions. The variability in the free troposphere is smaller compared to at the surface as well. It is also because of the previously noted difference between the spatial range that remote measurements cover being more akin to the spatial scale of model grid boxes compared to the point measurements from in situ observations.

5.1.1 What do the best-performing models have in common?

There were no models that performed best for all SLCF species and for all regions, highlighting the fact that it is difficult for any one model to bring together numerous complex processes and get results comparable to observations for all SLCFs. This would involve simulating aerosols and chemistry together with the right transport processes, meteorology, and clouds, which is difficult, especially for a remote region like the Arctic where parameterizations might have been built on datasets that are not always applicable there. In addition, studies like that of have shown that there was no clear change in model skill (in that case for OA) with increasing model complexity.

However, several models such as CanAM5-PAM, DEHM, NorESM, and MATCH have better representation of the vertical distribution of BC. DEHM and MATCH also had relatively small biases throughout the O${}_{\mathrm{3}}$ tropospheric profile (CanAM5-PAM and NorESM did not simulate gas-phase SLCFs). MATCH in particular has the smallest surface O${}_{\mathrm{3}}$ bias at midlatitudes, which may be related to its high vertical resolution in the boundary layer (the lowest two layers are 20 m thick and the four lowest layers are below 150 m). These models are a mix of air quality and climate-focused models; thus, it is important to note that there is no obvious difference between climate and air quality model biases for annual mean SLCFs. Despite the lack of complex tropospheric chemistry, CMAM had some of the lowest O${}_{\mathrm{3}}$ biases at both midlatitudes and the Arctic. This may imply that the more complex chemistry is not needed in the context of climatological tropospheric O${}_{\mathrm{3}}$ for climate studies (though of course, O${}_{\mathrm{3}}$ on shorter timescales would need more complete tropospheric O${}_{\mathrm{3}}$ chemistry). In the lower stratosphere, however, models with simplified or climatological O${}_{\mathrm{3}}$ schemes did not perform as well as models that had full stratospheric chemistry included.

5.1.2 Are there regional patterns in the model biases?

Generally speaking, when comparing the midlatitude model biases to those of the Arctic, they skew more negative (Fig. ), implying a lack of long-range transport to the Arctic. The best Arctic results for BC throughout the tropospheric profile were from models (CanAM5-PAM, DEHM, MATCH, WRF-Chem) that simulated the vertical mixing of BC well at midlatitudes. The first three of the four models were nudged to the ERA-Interim analysis, and WRF-Chem was nudged to the NCEP FNL analysis. The study showed that a key determinant in model differences for peroxyacetyl nitrate (PAN) export relative to CO was the meteorology used in the models. Their results implied that the ERA-Interim models had more efficient vertical transport and mixing in midlatitude source regions compared to GEOS-driven models. In the current study, Arctic BC was greatly underestimated throughout the Arctic profile by MATCH-SALSA, EMEP MSC-W, and GEM-MACH, which had BC concentrations at midlatitudes dropping much too low in the free troposphere. Those models had different sources of meteorology (Table ), and some may have insufficient convection schemes. For example, GEM-MACH is missing subgrid-scale deep and shallow convection, which is important for the exchange between the planetary boundary layer and free troposphere and consequently for transport at middle and high latitudes. This subject could be studied further via sensitivity tests with and without nudged meteorology, while keeping the aerosol physics the same.

The summertime evaluations of surface O${}_{\mathrm{3}}$, BC, and OA all imply that models overestimate the amount of these pollutants coming from wildfires in the western Arctic. This could be due to uncertainties in wildfire emissions, fire plume transport, or, in the case of O${}_{\mathrm{3}}$ and secondary OA, the plume chemistry. It is also likely that wet deposition of fire pollutants is underestimated if fire aerosol size is too small and due to climate models lacking pyrocumulus clouds and precipitation. The model overestimations of SLCFs in the summer could be due to a combination of all of these uncertainties.

Our analysis from ACE-FTS, TES, and the aircraft datasets show that CH${}_{\mathrm{4}}$, O${}_{\mathrm{3}}$, and BC model biases all imply that the modeled tropopause height is likely too low. Tropospheric species like CH${}_{\mathrm{4}}$ and BC should drop rapidly above the tropopause, but model biases increase sharply at that point. Stratospheric species like O${}_{\mathrm{3}}$ should increase rapidly above the tropopause, but model biases decline sharply there.

5.1.3 Are there patterns in the model biases between SLCF species?

Some patterns one might expect between SLCF species were not demonstrated in this study's results. For example, both O${}_{\mathrm{3}}$ precursors, CO and NO${}_{\mathrm{2}}$, are too low in models, though that underestimation is worse in the winter. Despite that, surface O${}_{\mathrm{3}}$ tends to be overestimated, though that overestimation is mainly in the summer. Also, the west–east pattern in Arctic surface O${}_{\mathrm{3}}$ was overestimated in Alaska and underestimated in Scandinavia. However, for CO and NO${}_{\mathrm{2}}$ those skewed biases are reversed.

At midlatitudes SO${}_{\mathrm{4}}^{\mathrm{2}-}$, BC, and OA are biased high in most models, yet despite that, PM${}_{\mathrm{2.5}}$ is underestimated. Thus, the PM${}_{\mathrm{2.5}}$ biases must be significantly influenced by the other, mainly natural, aerosol components.

In addition to expected patterns, there were no other discernible patterns in model biases between SLCF species.

5.2 How does the model performance impact model applications, such as simulated climate and health impacts?

In the report, these models go on to be used to simulate future emission scenarios, and from those results, the future temperature changes due to these SLCFs are predicted. They are further used to determine future changes to human health due to the changes in SLCFs. Given the model evaluation of this study, we have determined that using the mmm to predict SLCF climate impacts is generally robust. Considering the SLCFs with the greatest radiative impact (CH${}_{\mathrm{4}}$, O${}_{\mathrm{3}}$, BC, and SO${}_{\mathrm{4}}$), the mmm was within $\pm$25 % of the measurements across the Northern Hemisphere. The mmm also performed very well for PM${}_{\mathrm{2.5}}$, which is a main component considered for human health impacts. Thus, for climate studies, wherein relatively large regions and time periods are considered, the model performance is sufficient.

5.3 What processes should be improved or studied further for better model performance?

The model evaluation in this study brought about results that have been reported in previous publications, and several notable issues remain. Here we recommend some future work that may help improve model performance.

• Models simulate too much surface O${}_{\mathrm{3}}$ at midlatitudes, and this may be due inadequate treatment of dry deposition and/or not including parameterizations for the shade provided by vegetation that reduces photochemistry, as reported in . That said, MATCH had the smallest midlatitude surface O${}_{\mathrm{3}}$ bias without accounting for canopy shading; hence, precursor emissions, vertical mixing, deposition, and O${}_{\mathrm{3}}$ chemistry all have a role in model O${}_{\mathrm{3}}$ results, and errors in these may sometimes cancel out.

• There are a number of indications that simulated boundary layers are too stable (not enough vertical lifting of SLCFs, too much O${}_{\mathrm{3}}$ titration, too much BC and SO${}_{\mathrm{4}}$ deposition). Therefore, increased convection at midlatitudes may be needed. However, this hypothesis is opposite to that found in , which found that excessive tropical convection caused CMIP5 models to overestimate BC aloft. It is therefore important to evaluate models specifically for export and long-range transport events driven by different mechanisms (e.g., frontal export, convective lifting), which is a focus within the PACES initiative .

• The O${}_{\mathrm{3}}$, BC, and PM${}_{\mathrm{2.5}}$ model biases were all high in the Alaskan summertime, implying that many models may simulate too much pollution from wildfires there. Models need improved wildfire parameters for emissions, plume height, plume chemistry, and aerosol–cloud processes. For example, the fire emissions inventories GFED4, GFASv1.2, and FINNv1.5 vary by up to a factor of 3 for BC emissions ; light and temperature attenuation under smoke plumes means less O${}_{\mathrm{3}}$ is produced than precursor concentrations may imply, and plume rise and injection height need to be accurate for long-range transport.

• Modeled deposition is highly uncertain, and there is evidence here that some models have too much deposition of BC at midlatitudes. However, deposition measurements are scarce, even at midlatitudes, and more of those measurements are needed to constrain models.

• Additional, preferably long-term, OA and PM${}_{\mathrm{2.5}}$ measurements are needed in the high Arctic. Both are expected to be important to Arctic conditions in the future, with increasing wildfires and shipping influencing the Arctic atmosphere, and the lack of those measurements is problematic for constraining models.

• An evaluation of SO${}_{\mathrm{2}}$ would help to determine if the model biases in SO${}_{\mathrm{4}}^{\mathrm{2}-}$ are due to transport, emission uncertainty, or if it can be explained by the uncertainty in chemistry. The removal of particles represents a large uncertainty, but without SO${}_{\mathrm{2}}$ (and DMS) it cannot be concluded that the removal is too fast.

Therefore, we conclude that sensitivity tests for the abovementioned model processes will be important for further understanding and improving model performance for SLCFs. But just as important is having additional Arctic measurements and the continuation of existing Arctic measurements in order to assess the model improvements.

Appendix A Model descriptions

The 18 models used in this study are described in each subsection below. Table  contains further information summarizing the models' setup for the AMAP SLCF simulations.

Table A1

Information about models' spatial setups.

Model	Horizontal resolution	Scale (global or regional)
CanAM5-PAM	128$\times$64, Gaussian grid, T63	global
CESM	1.9${}^{\circ }$ $\times$ 2.5${}^{\circ }$ lat–long grid	global
CIESM-MAM7	0.9${}^{\circ }$ $\times$ 1.25${}^{\circ }$ lat–long grid	global
CMAM	96 $\times$ 48 Gaussian grid, T47	global
DEHM	50 km, $>$ 150 $\times$ 150 grid points	polar stereographic
ECHAM-SALSA	T63	global
EMEP-MSC-W	0.5${}^{\circ }$ $\times$ 0.5${}^{\circ }$ regular long–lat	global
FLEXPART	met. input data: 1${}^{\circ }$ $\times$ 1${}^{\circ }$	global
GEM-MACH	0.1375${}^{\circ }$ (or 15 km)	rotated Arctic LAM
GEOS-Chem	2${}^{\circ }$ $\times$ 2.5${}^{\circ }$	global
GISS-E2.1	2${}^{\circ }$ $\times$ 2.5${}^{\circ }$	global
MATCH	186 $\times$ 186, 0.75${}^{\circ }$	rotated lat–long regional
MATCH-SALSA	188 $\times$ 198, 0.75${}^{\circ }$	rotated lat–long regional
MRI-ESM2	TL159 (AGCM), TL95 (aerosol), T42 (ozone)	global
NorESM	0.9${}^{\circ }$ $\times$ 1.25${}^{\circ }$	global
Oslo-CTM	2.25${}^{\circ }$ $\times$ 2.25${}^{\circ }$	global
UKESM1	145x192 (1.875${}^{\circ }$ $\times$ 1.25${}^{\circ }$)	global
WRF-Chem	100 km	regional–Arctic

A1 CanAM5-PAM

The Canadian Atmospheric Model version 5 (CanAM5), with the Piecewise lognormal approximation Aerosol Model (PAM), was used. CanAM5-PAM is an improved version of CanAM4 . The improvements include a higher vertical resolution, improved parameterizations for land surface and snow processes, DMS emissions, and clear-sky radiative transfer. CanAM5-PAM has 49 vertical levels extending up to 1 hPa with a resolution of approximately 100 m near the surface. Model simulations are performed using a spectral resolution of T63, which is equivalent to the horizontal resolution of approximately 2.8${}^{\circ }$ $\times$ 2.8${}^{\circ }$. The model uses separate parameterizations for layer and convective clouds. Aerosol microphysical processes are based on the piecewise lognormal approximation . The model simulates binary homogeneous nucleation of sulfuric acid and water vapor. Newly formed particles grow by condensation and coagulation.

A detailed description of parameterizations of ocean DMS flux to atmosphere, oxidation, and removal processes is provided in . In-cloud production of sulfate requires O${}_{\mathrm{3}}$ and hydrogen peroxide (H${}_{\mathrm{2}}$O${}_{\mathrm{2}}$) as oxidants , with oxidant (OH, NO${}_{\mathrm{3}}^{-}$, H${}_{\mathrm{2}}$O${}_{\mathrm{2}}$, O${}_{\mathrm{3}}$) concentrations specified as climatological results from CMAM. Dry deposition of aerosol depends on concentrations of aerosols in the near-surface model layer . Wet deposition includes in-cloud scavenging in both convective clouds and layer clouds, as well as below-cloud scavenging.

Cloud droplet number concentrations are calculated based on the assumption of a parcel of air which ascends from the subcloud layer into the cloud layer with a characteristic vertical velocity ; the standard deviation of the subgrid-scale cloud vertical velocity probability distribution is parameterized using the approach by . Aerosol particles that are suspended in the parcel of air may activate and grow into cloud droplets by condensation of water vapor. A numerically efficient solution of the condensational droplet growth equation (e.g., ) is employed for this purpose. In grid cells that are affected by clouds, CanAM5-PAM accounts for cloud albedo and lifetime effects (first and second aerosol indirect effects) as well as semi-direct effects.

A2 CESM

The Community Earth System Model version 2 is an ESM that can be configured in many different ways. The configuration applied for this assessment utilized the Community Atmosphere Model (CAM) version 6 and Modal Aerosol Model (MAM4) with four mixed-species aerosol modes . CAM6 employs a spectral element dynamical core . Type 0 and Type 1 CESM runs were conducted at 1.9${}^{\circ }$ $\times$ 2.5${}^{\circ }$ horizontal resolution, while Type 3 runs were conducted at 0.9${}^{\circ }$ $\times$ 1.25${}^{\circ }$, all with 32 vertical layers. For Type 0 and Type 1 simulations, CESM version 2.0 was used with CAM6-chem representations of chemical reactions , enabling prognostic simulation of tropospheric ozone concentrations, along with a volatility basis set (VBS) parameterization for the formation of secondary organic aerosol (SOA) and stratospheric chemistry. CAM6-chem is coupled to the interactive Community Land Model (CLM5), which provides biogenic emissions, calculated online using the MEGANv2.1 algorithm , and handles dry deposition. Tracked aerosol species simulated by MAM4 include sulfate, primary and aged black carbon and organic matter, dust, sea salt, and secondary organic aerosols. Both sea salt and dust emissions are calculated on-line and are highly sensitive to the surface wind speed . These runs were also forced with prescribed sea surface temperatures (SSTs) and sea ice concentrations, created from merged Reynolds–HADISST products as in . Type 3 transient runs utilized CESM version 2.1.1 without atmospheric chemistry and with fully coupled atmosphere, ocean, land, and sea ice components (component set BSSP245cmip6), as applied to simulate future scenarios for CMIP6. All CESM runs specified global mean mixing ratios of methane and carbon dioxide.

A3 CIESM-MAM7

CIESM-MAM7 is the Community Integrated Earth System Model (CIESM) using the Modal Aerosol Model (MAM7) with seven mixed-species aerosol modes . The current CIESM version 1.1 (see Table 1 of ) is based on the NCAR Community Earth System Model (CESM version 1.2.1) with several novel developments and modifications aiming to overcome some persistent systematic biases, such as the double Intertropical Convergence Zone (ITCZ) problem and underestimated marine boundary layer clouds. CIESM-MAM7 employs a finite-volume dynamical core with 0.9${}^{\circ }$ $\times$ 1.25${}^{\circ }$ for horizontal resolution and 31 layers for vertical resolution. The large-scale meteorology (horizontal wind field) is nudged towards ERA-Interim reanalysis data and the relaxation time is set to 6 h. In CIESM-MAM7, the primary emissions of black carbon (BC), organic carbon (OC), ammonia (NH${}_{\mathrm{3}}$), volatile organic compounds (VOCs), sulfur dioxide (SO${}_{\mathrm{2}}$), and oxidizing gases (H${}_{\mathrm{2}}$O${}_{\mathrm{2}}$, O${}_{\mathrm{3}}$, OH) are prescribed by the input data uniformly provided by the AMAP-SLCF group. The emission amounts of dust (DU) and sea salt (SS) are calculated online. Aerosol size distributions in CIESM-MAM7 are described by the seven overlapping lognormal distributions, including Aitken, accumulation, primary carbon, fine dust and sea salt, coarse dust, and sea salt modes. The geometric standard deviation of each mode is prescribed (see Table 1 of ). A simplified gas- and liquid-phase chemistry is included in CIESM-MAM7. SO${}_{\mathrm{2}}$ and dimethyl sulfate (DMS) can be oxidized to sulfuric acid gas (H${}_{\mathrm{2}}$SO${}_{\mathrm{4}}^{\mathrm{2}-}$) and then condenses to form the sulfate aerosols, while the evolution of oxidizing gases is not considered. Primary organic matter (POM) and BC are emitted to the primary carbon mode, then aged and transferred to the accumulation mode by condensation of H${}_{\mathrm{2}}$SO${}_{\mathrm{4}}^{\mathrm{2}-}$, NH${}_{\mathrm{3}}$, and semi-volatile organics and by coagulation with Aitken and accumulation modes. The effect of stratospheric sulfate aerosol from volcanic emission on radiative forcing is considered by following the CMIP6 procedure . No specific stratospheric chemistry is included in CIESM-MAM7.

A4 CMAM

The Canadian Middle Atmosphere Model (CMAM) is based on the third-generation CanAM model, with the model lid raised to approximately 95 km and the necessary radiative processes for the mesosphere included . A representation of gas-phase chemistry has also been included that contains a relatively complete description of the HO${}_x$, NO${}_x$, Cl${}_x$, and Br${}_x$ chemistry that controls stratospheric ozone along with the longer-lived source gases such as CH${}_{\mathrm{4}}$, N${}_{\mathrm{2}}$O, and CFCs . For the troposphere the chemical mechanism can considered as methane–NO${}_x$ chemistry as it does not include the chemistry of larger volatile organic compounds. The model does, however, include a description of associated tropospheric chemical processes such as wet and dry deposition, interactive NO${}_x$ emissions from lightning, corrections of clear-sky photolysis rates for clouds, and N${}_{\mathrm{2}}$O${}_{\mathrm{5}}$ hydrolysis on prescribed sulfate aerosol distribution using the reaction probabilities of . The simulation analyzed here used a “specified dynamics” setup, wherein the model horizontal winds and temperature are nudged towards a meteorological reanalysis dataset that represents the observed historical evolution of the atmosphere. In this way the day-to-day variability of the model meteorology is much more closely aligned with the historical evolution of the atmosphere than would be possible in a free-running model. Here CMAM was nudged to 6-hourly fields from the ERA-Interim reanalysis on all model levels below 1 hPa and with a relaxation time constant of 24 h.

A5 DEHM

The Danish Eulerian Hemispheric Model (DEHM) is a 3-D Eulerian atmospheric chemistry transport model developed at the Department of Environmental Science at Aarhus University in Denmark. The model domain covers the Northern Hemisphere using a polar stereographic projection with a grid resolution of 150 km $\times$ 150 km. It includes nesting capabilities to make simulations with a higher grid resolution in a limited area of the domain, and in this work an Arctic sub-domain with 50 km $\times$ 50 km has been applied covering the Arctic area down to about 40–54${}^{\circ }$ N. The model has 29 vertical levels in sigma coordinates; the lowest 15 levels are below 2000 m above the surface. The lowest model levels are 22 m thick, and the top of the model domain is at 100 hPa, i.e., the whole troposphere and very lowest part of the stratosphere. DEHM includes a SO${}_x$–NO${}_x$–VOC–ozone chemistry, with 71 components including secondary organic aerosol (SOA), the use of the VBS mechanism, and nine particulates including hydrophobic and hydroscopic BC, primary organic aerosols, primary anthropogenic dust, PM${}_{\mathrm{2.5}}$ fraction, and coarse fraction of PM${}_{\mathrm{10}}$ of sea salt and Pb. CH${}_{\mathrm{4}}$ is a prognostic species for which the boundary conditions have a large influence. The model is driven by meteorological data from a numerical weather prediction model from the WRF model version 3.9, with 1 h resolution. The WRF model system is driven by reanalysis data from the ERA-Interim made by ECMWF by nudging.

A6 ECHAM-SALSA

ECHAM-SALSA is the general aerosol–climate model ECHAM-HAMMOZ (ECHAM6.3-HAM2.3-MOZ1.0) using the Sectional Aerosol module for Large Scale Applications (SALSA) to solve the aerosol microphysics. ECHAM6 computes the atmospheric circulation and fluxes using a semi-Lagrangian transport scheme. In the setup used here, the large-scale meteorology (vorticity, divergence, and surface pressure; relaxation times of 24, 6, and 48 h, respectively) was nudged towards ERA-Interim reanalysis data . In SALSA the aerosol size distribution is modeled using 10 size sections (or bins), which span particle sizes between 3 nm and 10 $\mathrm{\mu }$m. The size distribution is further divided into a soluble and an insoluble sub-population, which are treated as externally mixed. Within one size bin of one sub-population, all aerosol particles are considered internally mixed. In its standard setup, SALSA describes the aerosol compounds, BC, organic carbon (OC), SO${}_{\mathrm{4}}^{\mathrm{2}-}$, SS, and mineral dust (DU). In the model, BC, OC, SS, and DU are emitted as primary particles, while SO${}_{\mathrm{4}}^{\mathrm{2}-}$ is emitted as either SO${}_{\mathrm{2}}$ or as DMS, which are oxidized using a simplified chemistry to form H${}_{\mathrm{2}}$SO${}_{\mathrm{4}}$, which then either nucleates or condenses onto existing particles. BC, OC, and SO${}_{\mathrm{2}}$ emissions are prescribed using input files, while SS and DU emissions are computed online. All greenhouse gas concentrations are fixed to pre-defined concentrations. The model resolution for the simulations performed here was T63 (roughly 2${}^{\circ }$ by 2${}^{\circ }$), further using 47 hybrid sigma-pressure levels.

A7 EMEP MSC-W

The EMEP MSC-W model is a 3-D Eulerian chemistry transport model developed at the Norwegian Meteorological Institute within the framework of the UN Convention on Long-range Transboundary Air Pollution. It is described in detail in . Although the model has traditionally been aimed at simulations of acidification, eutrophication, and air quality over Europe, global modeling has been performed and evaluated against observations for many years . The model uses 20 vertical levels defined as eta-hybrid coordinates. The 10 lowest levels are within the planetary boundary layer (with the bottom layer being 92 m thick), and the top of the model domain is at 100 hPa. Model updates since , resulting in EMEP model version rv4.33 as used here, have been described in and references cited therein. The main revisions were made to the parameterizations of coarse NO${}_{\mathrm{3}}^{-}$ formation on sea salt and dust aerosols, N${}_{\mathrm{2}}$O${}_{\mathrm{5}}$ hydrolysis on aerosols, and additional gas–aerosol loss processes for O${}_{\mathrm{3}}$, HNO${}_{\mathrm{3}}^{-}$, and HO${}_{\mathrm{2}}$. The EMEP model, including a user guide, is publicly available as open-source code at https://github.com/metno/emep-ctm (last access: 14 April 2022). EMEP-modeled PM${}_{\mathrm{2.5}}$ and PM${}_{\mathrm{10}}$ include primary and secondary aerosols, both anthropogenic and natural. Secondary aerosol consists of inorganic sulfate, nitrate, ammonium, and SOA; the latter is formed from both anthropogenic and biogenic emissions using the VBS scheme detailed in and . The model also calculates sea salt aerosols and windblown dust particles from soil erosion. Aerosol optical depth (AOD) is calculated based on the mass concentrations of individual aerosols multiplied by corresponding mass extinction coefficients. In these simulations, we did not use the BC and OC emissions from EclipseV6b directly, but applied EclipseV6b PM${}_{\mathrm{2.5}}$ and coarse PM emissions instead, which were split into elementary carbon (EC), organic matter (OM) (here assumed inert), and the remaining inorganic dust. The EC and OM emissions in the fine and coarse fractions were further divided into fossil fuel and wood-burning compounds for each country and source sector. The split applied to the PM emissions is the same as used in EMEP operational runs (IIASA, personal communication). A total of 80 % of emitted EC is assumed to be hydrophobic, aging to become hydrophilic within 1 to 1.5 d. As in , the organic matter to organic carbon ratio of emissions by mass is assumed to be 1.3 for fossil fuel sources and 1.7 for wood-burning sources. Note that different wildfire emissions were used here, i.e., from FINN (the Fire INventory from NCAR version 15). The EMEP model runs were driven by 3-hourly meteorological data from the ECMWF IFS model at 0.5${}^{\circ }$ $\times$ 0.5${}^{\circ }$ resolution.

A8 FLEXPART

The Lagrangian particle dispersion model FLEXPART version 10.4 releases computational particles that are simulated forward in time following 3-hourly ECMWF meteorological fields with 137 vertical layers and a spatial resolution of 1${}^{\circ }$ $\times$ 1${}^{\circ }$. For each year around 330 million particles were released to calculate turbulent diffusion , unresolved mesoscale motions , and convection . A recently updated wet deposition scheme taking into account in-cloud and below-cloud removal was used . Gravitational settling for spherical BC particles with an aerosol mean diameter of 0.25 $\mathrm{\mu }$m, a normalized standard deviation of 3.3, and a particle density of 1500 kg m${}^{-\mathrm{3}}$ is used in the calculation of dry deposition. The surface concentration and deposition fields were retrieved on a monthly basis on a resolution of 0.5${}^{\circ }$ $\times$ 0.5${}^{\circ }$.

A9 GEM-MACH

GEM-MACH (Global Environmental Multiscale model–Modelling Air quality and CHemistry) is the Environment and Climate Change Canada (ECCC) air quality prediction model. It consists of an online tropospheric chemistry module embedded within ECCC's GEM numerical weather forecast model . The chemistry module includes a comprehensive representation of air quality processes, such as gas-phase, aqueous-phase, and heterogeneous chemistry and aerosol processes (e.g., ). Specifically, gas-phase chemistry is represented by a modified ADOM-II mechanism with 47 species and 114 reactions ; inorganic heterogeneous chemistry is parameterized by a modified version of the ISORROPIA algorithm of , as described in detail in ; SOA formation is parameterized using a two-product, overall, or instantaneous aerosol yield formation ; aerosol microphysical processes, including nucleation and condensation (sulfate and SOA), hygroscopic growth, coagulation, and dry deposition and sedimentation, are parameterized as in ; the representation of cloud processing of gases and aerosols includes uptake and activation, aqueous-phase chemistry, and wet removal .

Aerosol chemical composition is represented by eight components: sulfate, nitrate, ammonium, elemental carbon (EC), primary organic aerosol (POA), secondary organic aerosol (SOA), crustal material (CM), and sea salt; aerosol particles are assumed to be internally mixed. A sectional approach is used for representing aerosol size distribution. For the 2015 Arctic simulation, a 12-bin (between 0.01 and 40.96 $\mathrm{\mu }$m in diameter, logarithmically spaced: 0.01–0.02, 0.02–0.04, 0.04–0.08, 0.08–0.16, 0.16–0.32, 0.32–0.64, 0.64–1.28, 1.28–2.56, 2.56–5.12, 5.12–10.24, 10.24–20.48, and 20.48–40.96 $\mathrm{\mu }$m) configuration is used.

The Type 0 simulation was conducted for the year 2015 over a limited-area model (LAM) domain on a rotated lat–long grid at 0.1375${}^{\circ }$ $\times$ 0.1375${}^{\circ }$ (or 15 km) resolution covering the Arctic ($>$ 60${}^{\circ }$ N) and extending to the southern US–Canada border. Some of the model upgrades for the Arctic simulation are described in . Anthropogenic emissions used are based on a combination of North American emission inventories (specifically, the 2016 US National Emission Inventories and 2015 Canadian national Air Pollution Emission Inventories) and global ECLIPSE v6b 2015 baseline emissions. North American wildfire emissions are processed using the Canadian Forest Fire Emission Prediction System (CFFEPS) from satellite-detected fire hotspot data (MODIS, AVHRR, and VIIRS). CFFEPS consists of a fire growth model, a fire emissions model, and a thermodynamic-based model to predict the vertical penetration height of a smoke plume from fire energy (see , for details). Biogenic emissions are calculated online in GEM-MACH based on the algorithm from BEIS version 3.09 with BELD3-format vegetation land cover. Sea salt emissions are computed based on .

The chemical lateral boundary conditions were from MOZART-4/GEOS-5 (). The meteorology was initialized daily (at 00:00 UTC) using the Canadian Meteorological Centre's global objective analyses.

A10 GEOS-Chem

GEOS-Chem is a global three-dimensional chemical transport model driven by assimilated meteorological observation from the Goddard Earth Observing System (GEOS) of the NASA Data Assimilation Office (DAO), which was first introduced in 2001 . GEOS-Chem is a grid-independent model which operates on a 1-D column with default or user-specified horizontal grid points, vertical grid points, and time step. GEOS-Chem Classic can use archived GEOS meteorological data on a rectilinear latitude–longitude grid to compute horizontal and vertical transport and use Open-MP in parallelization. Two types of assimilated meteorological data from the NASA Global Modeling and Assimilation Office (GMAO) can be used to drive the offline mode of GEOS-Chem. The first type are operational data starting from 2012, the GEOS Forward Processing (GEOS-FP, the native resolution of which was 0.25${}^{\circ }$ $\times$ 0.3125${}^{\circ }$; ). The second type is the consistent Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA-2, starting from 1979–present; ), with the native resolution 0.5${}^{\circ }$ $\times$ 0.625${}^{\circ }$. Both types of meteorological data have 72 hybrid sigma-pressure levels with the top at 0.01 hPa, 3-hourly temporal resolution for 3-D fields, and 1 h resolution for 2-D fields. The advection scheme of GEOS-Chem uses the TPCORE advection scheme on the latitude–longitude grid, while the convective transport uses the convective mass flux described by . The wet deposition scheme in GEOS-Chem is based on for water-soluble aerosols and for gases. The dry deposition is based on the resistance-in-series scheme of . Aerosol deposition is from . Emissions of dust aerosol, lightning NO${}_x$, biogenic VOCs, soil NO${}_x$, and sea salt aerosol are dependent on the local meteorological conditions. The CEDS global inventory provides the default anthropogenic emissions, while EDGAR v4.3.2 is also available as an alternate option to CEDS. Future anthropogenic emissions following the RCP scenarios (), aircraft emissions , ships emission (from CEDS), and lighting NO${}_x$ emissions are also included and configured at run time using the HEMCO module described . Biogenic VOC emissions in GEOS-Chem are from the MEGAN v2.1 inventory . The chemical solver in the standard GEOS-Chem simulation uses the Kinetic PreProcessor (KPP) as implemented in GEOS-Chem. The gas phase in the troposphere in GEOS-Chem includes a detailed HO${}_x$–NO${}_x$–VOC–ozone–halogen–aerosol tropospheric chemistry mechanism, which generally follows JPL/IUPAC recommendations including PAN , isoprene , halogens , and Criegees . A linearized stratospheric chemistry scheme has been implemented since GEOS-Chem v9.0. The model will read from archived 3-D monthly mean production rates and losing frequency for each species at the beginning of each month. The Linoz chemistry is also applied as a recommended option for the stratospheric ozone layer. The original sulfate–nitrate–ammonium aerosol simulation in GEOS-Chem coupled to gas-phase chemistry . The black carbon simulation , organic aerosol , complex SOA , the aqueous-phase isoprene SOA scheme , and the dust simulation are also implemented into GEOS-Chem. The dust size distributions are from . GEOS-Chem v12.3.2 with uniform 2${}^{\circ }$ $\times$ 2.5${}^{\circ }$ MERRA-2 meteorological data for 2008–2009, GEOS-FP meteorological data for 2014–2015, and ECLIPSEv6b emissions was used in this study.

A11 GISS-E2.1

The NASA Goddard Institute of Space Studies (GISS) Earth system model (ESM), GISS-E2.1, is a fully coupled ESM. A full description of GISS-E2.1 as well as evaluation of its coupled climatology during the satellite era (1979–2014) and the historical ensemble simulation of the atmosphere and ocean component models (1850–2014) are described in . GISS-E2.1 has a horizontal resolution of 2${}^{\circ }$ in latitude by 2.5${}^{\circ }$ in longitude and 40 vertical layers extending from the surface to 0.1 hPa in the lower mesosphere. The tropospheric chemistry scheme used in GISS-E2.1 includes inorganic chemistry of O${}_x$, NO${}_x$, HO${}_x$, CO, and organic chemistry of CH${}_{\mathrm{4}}$ and higher hydrocarbons using the CBM4 scheme , as well as the stratospheric chemistry scheme , which includes chlorine and bromine chemistry together with polar stratospheric clouds. The meteorology was nudged to the NCEP reanalysis.

In the present work, we used the OMA (the one-moment aerosol scheme) . OMA is a mass-based scheme in which aerosols are assumed to remain externally mixed and have a prescribed and constant size distribution, with the exception of sea salt that has two distinct size classes and dust that is described by a sectional model with an option from four to six bins. The OMA scheme treats sulfate, nitrate, ammonium, carbonaceous aerosols (black carbon and organic carbon, including the NO${}_x$-dependent formation of SOA and methanesulfonic acid formation), dust, and sea salt. The model includes secondary organic aerosol production, as described by . The default dust configuration that is used in this work includes five bins, including a clay and four silt ones, from submicron to 16 $\mathrm{\mu }$m in size. The first three dust size bins can be coated by sulfate and nitrate aerosols . OMA only includes the first aerosol indirect effect. The aerosol number concentration that impacts clouds is obtained from the aerosol mass as described in .

The natural emissions of sea salt, DMS, isoprene, and dust are calculated interactively. Anthropogenic dust sources are not represented in ModelE2.1. Dust emissions vary spatially and temporally only with the evolution of climate variables like wind speed and soil moisture . The version of the model we use in this work uses prescribed sea surface temperature (SST) as well as sea ice thickness and extent during the historical period .

A12 MATCH

MATCH – Multiscale Atmospheric Transport and Chemistry – is an offline, Eulerian, 3-D chemistry transport model developed at the Swedish Meteorological and Hydrological Institute. MATCH can be run on global to urban domains to study a range of atmospheric chemistry and air quality problems, but for this study model runs were performed for the 20–90${}^{\circ }$ N region focusing on long-transport to the Arctic. ERA-Interim reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF) were used as meteorological input to the model. The 6-hourly data (3-hourly for precipitation) were extracted from the ECMWF archives on a 0.75${}^{\circ }$ $\times$ 0.75${}^{\circ }$ rotated latitude–longitude grid. The original data had 60 levels, but the 38 lowest levels, reaching about 16 km in the Arctic, were used in the model.

The scheme for gas-phase tropospheric chemistry and bulk aerosols as described in was used. Methane concentrations were prescribed. Boundary conditions at the top of the model and at the lateral boundaries for a range of species including ozone were based on monthly mean values from the Copernicus Atmospheric Monitoring Service. The aerosol scheme was extended with BC and OC simulated as three fractions: fresh, hydrophobic and aged, and hydrophilic. A total of 80 % of anthropogenic emissions from all sectors were emitted into the hydrophobic and 20 % into the hydrophilic fraction, except for fire–biomass combustion for which 100 % was emitted into the hydrophilic component following . Scavenging and aging were parameterized following ; i.e., aging is proportional to OH and scavenging in mixed-phase clouds is reduced. The hydrophobic fraction is assumed to be 5 % activated in the scavenging scheme, while the hydrophilic fraction is 100 % activated. If the clouds are mixed-phase, then the scavenging efficiency is scaled by the ratio of cloud ice water content to total cloud water content assuming zero scavenging for 100 % ice clouds.

A13 MATCH-SALSA-RCA4

The chemistry transport model, MATCH described above, is online-coupled to the aerosol dynamics model, SALSA . SALSA describes the whole chain from nucleation to the growth and deposition of particles and computes the size distribution, number concentration, and chemical composition of the aerosol species. A sectional representation of the aerosol size distribution is considered with three main size ranges (a: 3–50 nm, b: 50–700 nm, and c: $>$ 700 nm), and each range is again subdivided into smaller bins and into soluble and insoluble bins, adding up to a total of 20 bins. The seasonally varying emissions are based on the sector-wise ECLIPSE inventory. Isoprene emissions are modeled online depending on the meteorology based on the methodology by . The terpene emissions ($\mathit{\alpha }$-pinene) are taken from the modeled fields by the EMEP model. Sea salt is parameterized following the scheme of but modified for varying particle sizes, wherein the scheme is used if the particle diameter is 1 $\mathrm{\mu }$m and the scheme is used otherwise. The coupling of MATCH with SALSA and the evaluation of this model setup is described in detail in . A cloud activation model that computes 3-D CDNCs (cloud droplet number concentrations) based on the prognostic parameterization scheme of specifically designed for aerosol representation with sectional bins is embedded in the MATCH-SALSA model. This scheme simulates the efficiency of an aerosol particle to be converted to a cloud droplet depending on the number concentration and chemical composition of the particles given the updraft velocity and supersaturation of the air parcel. The updraft velocity is computed as the sum of the grid mean vertical velocity and turbulent kinetic energy (TKE) for stratiform clouds . These CDNCs are then offline-coupled to a regional climate model, RCA4 , that provides us with information on cloud properties such as cloud cover, cloud droplet radii, cloud liquid water path, and radiative fluxes. RCA4 is run with 6-hourly ERA-Interim meteorology, and the 3-hourly RCA4 meteorological fields along with the fields needed to calculate updraft velocity are used to drive the MATCH-SALSA cloud activation model. The CDNCs are then used to rerun the RCA4 model to obtain the cloud properties and radiative effects. The validation and more details of this model set up are described in .

A14 MRI-ESM2

MRI-ESM2 (Meteorological Research Institute (MRI) Earth System Model version 2.0, developed by the MRI of the Japan Meteorological Agency) consists of four major component models: an atmospheric general circulation model (MRI-AGCM3.5) with land processes, an ocean–sea ice general circulation model (MRI.COMv4), and aerosol and atmospheric chemistry models . However, we do not couple OGCM in this study's simulations. MRI-ESM2 uses different horizontal resolutions but employs the same vertical resolution in each atmospheric component model as follows: TL159 (approximately 120 km), TL95 (approximately 180 km), and T42 (approximately 280 km) in the MRI-AGCM3.5, the aerosol model, and the atmospheric chemistry model, respectively, all with 80 vertical layers (from the surface to a model top of 0.01 hPa) in a hybrid sigma-pressure coordinate system. Each component model is interactively coupled by a coupler, which enables an explicit representation of the effects of gases and aerosols on the climate system. The atmospheric chemistry component model in MRI-ESM2 is the MRI Chemistry–Climate Model version 2.1 (MRI-CCM2.1), which calculates the evolution and distribution of ozone and other trace gases in the troposphere and in the middle atmosphere. The model calculates a total of 90 gas-phase chemical species and 259 chemical reactions in the atmosphere. The aerosol component model in MRI-ESM2 is the Model of Aerosol Species in the Global Atmosphere mark-2 revision 4-climate (MASINGAR mk-2r4c) that calculates atmospheric aerosol physical and chemical processes and treats the following species: non-sea-salt sulfate, BC, OC, sea salt, mineral dust, and aerosol precursor gases (e.g., sulfur dioxide and dimethyl sulfide). The size distributions of sea salt and mineral dust are divided into 10 discrete bins, and the sizes of the other aerosols are represented by lognormal size distributions. The model assumes external mixing for all aerosol species; however, in the radiation process in MRI-AGCM3.5, hydrophilic BC is assumed to be internally mixed with sulfate with a shell-to-core volume ratio of 2; the optical properties of hydrophilic BC are calculated based on Mie theory with a core–shell aerosol treatment, in which a concentric BC core is surrounded by a uniform coating shell composed of other aerosol compounds . MRI-ESM2 employs a BC aging parameterization that calculates the variable conversion rate of BC from hydrophobic BC to hydrophilic BC, which generally depends on the production rate of condensable materials such as sulfate. In the radiation and cloud processes in MRI-ESM2, sulfate is assumed to be (NH${}_{\mathrm{4}}^{+}$)${}_{\mathrm{2}}$SO${}_{\mathrm{4}}^{\mathrm{2}-}$ and OC is assumed to be organic matter (OM) by lumping OC species using an OM-to-OC factor of 1.4. MRI-ESM2 represents the activation of aerosols into cloud droplets based on the parameterizations, and detailed descriptions and evaluations of the cloud processes and cloud representations in MRI-ESM2 are given by . Evaluations of the effective radiative forcing (ERF) of anthropogenic gases and aerosols in present-day conditions relative to preindustrial conditions globally and in the Arctic using MRI-ESM2 are given by . The simulations in this study were performed from January 2008 (or January 1990) to December 2015 after a 1-year spin-up run using the prescribed SST and sea ice data (provided by the AMIP experiment in CMIP6, https://www.wcrp-climate.org/modelling-wgcm-mip-catalogue/modelling-wgcm-mips-2/240-modelling-wgcm-catalogue-amip, last access: 14 April 2022). The horizontal wind fields were nudged toward the 6-hourly Japanese 55-year Reanalysis (JRA55) data (https://jra.kishou.go.jp/JRA-55/index_en.html, last access: 14 April 2022) in the simulation. We used the monthly anthropogenic emissions from the ECLIPSE V6B emission dataset and the monthly biomass burning emissions from the CMIP6 in the simulations. Major volcanic aerosols are given by the stratospheric aerosol dataset used in the CMIP6 experiments . A second simulation with volcanic SO${}_{\mathrm{2}}$ emission including Holuhraun eruption was also performed for 2014–2015.

A15 NorESM1

NorESM1 is based on the fourth version of the Community Climate System Model (CCSM4) , with coupled models for the atmosphere, ocean, land, and sea ice. Here, we have used a 1${}^{\circ }$ horizontal resolution in the atmosphere (0.95${}^{\circ }$ latitude by 1.25${}^{\circ }$ longitude, version NorESM1-Happi). The model has 26 vertical levels on a hybrid sigma-pressure coordinate up to the model top at 2.194 hPa. The model calculates the life cycles of a range of natural and anthropogenic aerosol components from emissions and physicochemical processing in air and cloud droplets. The only prescribed aerosol concentrations are stratospheric sulfate from explosive volcanoes. The direct and indirect aerosol effects on climate are calculated by parameterization of aerosol interactions with schemes for radiation and warm cloud microphysics . The model uses a prognostic calculation of cloud droplet numbers, allowing for competition effects between aerosols of different hygroscopic property and size.

A16 OsloCTM

The Oslo CTM3 is an offline global three-dimensional chemistry transport model driven by 3-hourly meteorological forecast data from the Integrated Forecast System (IFS) model at the European Centre for Medium-Range Weather Forecasts (ECMWF). The Oslo CTM3 consists of a tropospheric and stratospheric chemistry scheme as well as aerosol modules for sulfate, nitrate, black carbon, primary organic carbon, secondary organic aerosols, mineral dust, and sea salt .

A17 UKESM1

UKESM1 (United Kingdom Earth System Model) is a fully coupled Earth system model with a coupled atmosphere–ocean physical climate model (HadGEM3-GC3.1) at its core . For UKESM1 various Earth system components are incorporated with the physical climate model including ocean biogeochemistry, an interactive stratosphere–troposphere chemistry and aerosol scheme, and terrestrial carbon and nitrogen cycles coupled to interactive vegetation. The model has a horizontal resolution of 135 km at the midlatitudes (1.875${}^{\circ }$ $\times$ 1.25${}^{\circ }$), with 85 levels on a terrain-following hybrid height coordinate system, ranging in height from the surface to a model top of 85 km. The combined stratosphere–troposphere United Kingdom Chemistry and Aerosol (UKCA) scheme is used within UKESM1 and is fully described and evaluated in .

The chemical scheme in UKCA is built upon the scheme described for the stratosphere in and that for the troposphere described in . Chemical reactions are included within UKCA for odd oxygen (O${}_x$), nitrogen (NOy), hydrogen (HO${}_x=$ OH $+$ HO${}_{\mathrm{2}}$), CO, CH${}_{\mathrm{4}}$, and short-chain non-methane volatile organic compounds (NMVOCs), including isoprene. Reactions involving NMVOCs are simulated as discrete species. UKCA includes an interactive photolysis scheme, as well as representations of both wet and dry deposition, for gas and aerosol species. Additional chemical reactions for DMS, SO${}_{\mathrm{2}}$, and monoterpenes (C${}_{\mathrm{10}}$H${}_{\mathrm{16}}$) are included to enable coupling to the aerosol scheme within UKCA. A two-moment aerosol microphysical scheme, GLOMAP (Global Model of Aerosol Processes; ), is used to simulate four aerosol components (SO${}_{\mathrm{3}}$, BC, organic matter, sea salt) across five lognormal modes, ranging from submicron to supermicron sizes. Mineral dust is simulated separately using a six-bin mass-only scheme, ranging in size from 0.6 to 60 $\mathrm{\mu }$m in diameter . Ammonium nitrate is not currently included within the UKCA aerosol scheme. The formation of secondary organic aerosol (SOA) is included based on a fixed yield rate of 26 % from the products of monoterpene oxidation. The higher fixed yield value accounts for the underlying uncertainty in SOA formation and the absence of anthropogenic, marine, and isoprene sources.

Precursor emission fluxes are either prescribed using specified input files or calculated interactively using online meteorological variables within UKESM1. Methane is represented by using prescribed global concentrations. Interactive emission fluxes are calculated online for sea salt, DMS, dust, lightning NO${}_x$, and biogenic volatile organic compounds (BVOCs). Emissions of isoprene and monoterpenes from the natural environment are calculated online by coupling to the land surface scheme within UKESM1. Simulations provided by UKESM1 and used in the AMAP assessment have been undertaken using different configurations. For this study, the experiment with UKESM1 has been set up using an atmosphere-only configuration that is nudged to ECMWF reanalysis (ERA-Interim) of temperature and wind fields above the boundary layer. Prescribed values of sea surface temperatures and sea ice are used for each year of simulation based on historical simulations conducted as part of CMIP6 using the fully coupled atmosphere–ocean configuration of UKESM1. For other ancillary inputs a multi-year climatology was used, equivalent to an AMIP-type simulation.

A18 WRF-Chem

WRF-Chem (Weather Research and Forecasting model with online coupled chemistry) is used to simulate the transport and chemical transformation of trace gases and aerosols simultaneously with the meteorology. The model dynamics (WRF) are non-hydrostatic. The model version used for AMAP is WRF-Chem version 3.8.1 also including updates reported in and . The simulation was performed on a polar stereographic projection with a horizontal resolution of 100 km and 50 vertical hybrid terrain-following vertical pressure levels using hydrostatic pressure. The center of the domain is placed at the North Pole, and the latitude at the domain's outside boundary varies from 7${}^{\circ }$ S to 7${}^{\circ }$ N. The WRF-Chem chemical lateral boundary conditions are from MOZART-4/GEOS-5 (). Pressure at the model top is set to 50 hPa with stratospheric concentrations (e.g., ozone) taken from climatologies. The model was run with Morrison double-moment scheme microphysics, longwave and shortwave radiative effects treated by the RRTMG scheme, and the Kain–Fritsch cumulus potential (KF-CuP) cumulus parameterization scheme. The model was run with the SAPRC-99 chemical scheme providing gas-phase tropospheric reactions including VOCs and NO${}_x$, coupled with the MOSAIC eight-bin sectional scheme including VBS treatments for SOA. Methane concentrations are prescribed. Stratospheric or tropospheric halogen chemistry is not included. It was run using anthropogenic emissions from ECLIPSE v6b and the GFED fire emissions. Boundary and initial meteorological conditions were given by the global NCEP Final Analysis (FNL) and used to nudge the temperature, relative humidity, and winds at every dynamical time step above the planetary boundary layer.

Appendix B Observational datasets

The following datasets were used to evaluate the models in this study. As BC measurements vary by instrument, Table  summarizes the different Arctic BC datasets used in this study.

Table B1

Information about Arctic BC measurements used for model evaluation.

Location or network	Method	Comments/references
IMPROVE	EC via thermo-optical
EMEP	EC via thermo-optical from PM${}_{\mathrm{2.5}}$ and PM${}_{\mathrm{10}}$
CABM	EC via thermal evolution method from
	total suspended particle (2005–2011) and
	PM${}_{\mathrm{1}}$ (2011 to present). At Alert, also eBC via
	Aethalometer for PM1
Gruvebadet Lab	eBC via PSAP from PM1
Zeppelin Mountain	eBC via Aethalometer
Utqiaġvik (formerly Barrow)	eBC via Aethalometer and via PSAP from PM1
Japanese Arctic cruise	rBC via SP2 from PM10
Russian Arctic cruise	eBC via Aethalometer
Aircraft campaigns	rBC from SP2

For SLCFs other than Arctic BC, Table  summarizes some information about the observation networks.

Table B2

Information about SLCF measurements from monitoring networks used for model evaluation in the Northern Hemisphere.

Network	Long name	Species measured	Time period	Comments/references
acronym
CABM	Canadian Baseline Monitoring	SO${}_{\mathrm{4}}^{\mathrm{2}-}$, BC, and OC	2000 to present	six sites in Canada
	network
CSN	Chemical Speciation Network	O${}_{\mathrm{3}}$, NO${}_{\mathrm{2}}$, SO${}_{\mathrm{2}}$, CO, PM${}_{\mathrm{2.5}}$, PM${}_{\mathrm{10}}$,	2001 to present	data from the Air Quality
		SO${}_{\mathrm{4}}^{\mathrm{2}-}$, NO${}_{\mathrm{3}}^{-}$, NH${}_{\mathrm{4}}^{+}$, EC, and OC		System, which centralizes access to
				numerous US datasets
CAWNET	China Atmospheric Watch Network	PM${}_{\mathrm{2.5}}$, PM${}_{\mathrm{10}}$, SO${}_{\mathrm{4}}^{\mathrm{2}-}$, NO${}_{\mathrm{3}}^{-}$, NH${}_{\mathrm{4}}^{+}$,	2000 to present
		EC $/$ BC, OC, O${}_{\mathrm{3}}$, CO, NO${}_x$, CH${}_{\mathrm{4}}$, SO${}_{\mathrm{2}}$
EMEP	European Monitoring and	PM${}_{\mathrm{2.5}}$, PM${}_{\mathrm{10}}$, OC, BC, and SO${}_{\mathrm{4}}^{\mathrm{2}-}$	1993 to present	park of GAW
	Evaluation Programme
GAW	Global Atmosphere Watch	SO${}_{\mathrm{4}}^{\mathrm{2}-}$, BC, OC, O${}_{\mathrm{3}}$, NO${}_x$, CO,	1993 to present	a portal of measurements
		SO${}_{\mathrm{2}}$, AOD, CN
IMPROVE	Interagency Monitoring of Protected	SO${}_{\mathrm{4}}^{\mathrm{2}-}$, EC $/$ BC, OC, PM${}_{\mathrm{2.5}}$,	1987 to present	data obtained at
	Visual Environments	PM${}_{\mathrm{10}}$, and NO${}_{\mathrm{3}}^{-}$		http://views.cira.colostate.edu/fed/Auth/Login.aspx?ReturnUrl=/fed/Express/ImproveData.aspx (last access: 14 April 2022)
NAPS	National Air Pollution	PM${}_{\mathrm{2.5}}$, PM${}_{\mathrm{10}}$, SO${}_{\mathrm{2}}$, HNO${}_{\mathrm{3}}^{-}$, NH${}_{\mathrm{3}}$, SO${}_{\mathrm{4}}^{\mathrm{2}-}$,	1974 to present
	Surveillance Network	EC $/$ BC, OC, O${}_{\mathrm{3}}$, NO, NO${}_{\mathrm{2}}$, SO${}_{\mathrm{2}}$, CO
WHO	World Health Organization (WHO)	PM${}_{\mathrm{2.5}}$ and PM${}_{\mathrm{10}}$	1985 to 2011

Code and data availability

The models' output files in NetCDF format from the simulations used in this project can be found here: http://crd-data-donnees-rdc.ec.gc.ca/CCCMA/products/AMAP/ .

Some of the models' code are available online at the following locations. CanAM5-PAM: https://gitlab.com/cccma . CESM2: https://www.cesm.ucar.edu/models/cesm2/ . ECHAM-SALSA: The codes used for the ECHAM-SALSA simulations are available from the ECHAM-HAMMOZ repository under https://redmine.hammoz.ethz.ch/projects/hammoz/repository/1/show/echam6-hammoz/branches/fmi/AMAP/AMAP_evaluation , after obtaining the HAMMOZ license. FLEXPART: https://www.flexpart.eu . GEOS-Chem: http://wiki.seas.harvard.edu/geos-chem/index.php/GEOS-Chem_12#12.3.2 . GISS-E2.1: https://www.giss.nasa.gov/tools/modelE/ . NorESM: https://github.com/NorESMhub/NorESM . Oslo CTM: https://github.com/NordicESMhub/OsloCTM3 . The other models' code may be available upon request.

The model evaluation programs can be found on GitLab here: https://gitlab.com/cynwhaley/amap-slcf-model-evaluation .

The surface monitoring datasets are available online. WDCGG for CH${}_{\mathrm{4}}$: https://gaw.kishou.go.jp/login/user . EBAS for European (EMEP) and several Arctic locations: http://ebas.nilu.no/ . NAPS: https://open.canada.ca/data/en/dataset/1b36a356-defd-4813-acea-47bc3abd859b . IMPROVE: https://views.cira.colostate.edu/fed/Express/ImproveData.aspx . Beijing Air Quality for China: https://quotsoft.net/air/ . PM${}_{\mathrm{2.5}}$ from the US embassy in China from the data portal: http://www.stateair.net .

The satellite measurement data used in this study are available online. ACE-FTS v4.1 measurements are available, following registration, from http://www.ace.uwaterloo.ca . TES: https://tes.jpl.nasa.gov/tes/data/products/lite . MOPITT: https://www2.acom.ucar.edu/mopitt/products .

The supplement related to this article is available online at: https://doi.org/10.5194/acp-22-5775-2022-supplement.

Author contributions

CHW organized the model–measurement comparisons in the AMAP SLCF expert group, led the trace gas model evaluations, and wrote the paper with contributions from co-authors. RM developed and ran the aerosol model–measurement comparison scripts, including the ship-based BC measurements, and did the aerosol analysis. KvS led the AMAP SLCF modeling strategy and developed and ran the CanAM5-PAM model. SE and NE provided FLEXPART model output, and SE did the model–measurement comparisons for BC and SO${}_{\mathrm{4}}^{\mathrm{2}-}$ deposition. DWP developed the CIS tool for BC aircraft model–measurement comparisons, and BW used the tool and did the model–aircraft comparisons and analysis. LNS did the ACE-FTS trace gas model comparisons, and KAW provided guidance on ACE-FTS data quality, usage, and model–measurement comparisons. MF co-led the AMAP model strategy and ran CESM. LP provided additional CESM runs. DAP developed and ran CMAM. YP and MW provided CIESM-MAM7 model output; JC provided DEHM model output; TK provided ECHAM-SALSA model output; ST and MG provided EMEP-MSC-W model output and some analysis. UI, GF, and KT provided the GISS-E2.1 model output; WG and SB provided the GEM-MACH model output; JSF, RYC, and XD provided the GEOS-Chem model output; JL provided the MATCH model output; MAT provided the MATCH-SALSA model output; NO provided the MRI-ESM2 model output; MS and SK provided the NorESM model output; RS provided the OsloCTM model output; SA and ST provided UKESM1 model output; and KSL, JCR, TO, and LM provided WRF-Chem model output. SS and LH provided Alert datasets; VV and SB provided Gruvebadet datasets; OP provided Russian cruise ship measurements; FT and YK provided the Japanese ship measurements; HS and AM provided Villum Research Station datasets; JS provided aerosol datasets and expertise. TWG and SMD wrote and developed the trace gas model–measurement comparison scripts.

Competing interests

At least one of the (co-)authors is a member of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.

Disclaimer

PM${}_{\mathrm{2.5}}$ observation from the US embassy in China belongs to the US Department of State, is not fully verified or validated, and could be subject to changes, corrections, or errors. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “Arctic climate, air quality, and health impacts from short-lived climate forcers (SLCFs): contributions from the AMAP Expert Group (ACP/BG inter-journal SI)”. It is not associated with a conference.

Acknowledgements

The work reflected in this publication was produced with the financial support of the Arctic Monitoring and Assessment Programme (AMAP). The authors would like to thank all operators and technicians at the Arctic stations for the collection of observational data. Thanks are also due to the following people for providing us with data: Mauro Mazzola, Stefania Gilardoni, and Angelo Lupi from the Institute of Polar Sciences for eBC measurements at Gruvebadet lab; and Mirko Severi from the University of Florence for SO${}_{\mathrm{4}}^{\mathrm{2}-}$ measurements at Gruvebadet. Fairbanks aerosol measurements came from William Simpson. For Villum data, we acknowledge the Aarhus University Department of Environmental Science (ENVS). NOAA/ESRL/GMD, EMEP (http://ebas.nilu.no, last access: 14 April 2022), and the WMO GAW network are acknowledged for Barrow and Zeppelin observational datasets. The ACE-FTS is a Canadian-led mission mainly supported by the Canadian Space Agency (CSA). IMPROVE is a collaborative association of state, tribal, and federal agencies, as well as international partners. The US Environmental Protection Agency is the primary funding source, with contracting and research support from the National Park Service. The Air Quality Group at the University of California, Davis, is the central analytical laboratory, with ion analysis provided by the Research Triangle Institute and carbon analysis provided by the Desert Research Institute. Julia Schmale holds the Ingvar Kamprad Chair for Extreme Environments Research sponsored by Ferring Pharmaceuticals. The ECHAM-HAMMOZ model is developed by a consortium composed of ETH Zurich, the Max-Planck-Institut für Meteorologie, Forschungszentrum Jülich, the University of Oxford, the Finnish Meteorological Institute, and the Leibniz Institute for Tropospheric Research; it is managed by the Center for Climate Systems Modeling (C2SM) at ETH Zurich.

Financial support

Assessments from the Russian ship-based campaign were performed with the support of RFBR project no. 20-55-12001 and according to the development program of the Interdisciplinary Scientific and Educational School of M.V. Lomonosov Moscow State University “Future Planet and Global Environmental Change”. Development of the methodology for aethalometric data treatment was supported by RSF project no. 19-77-30004. The BC observations on R/V Mirai were supported by the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan (Arctic Challenge for Sustainability (ArCS) project). Contributions by SMHI were funded by the Swedish Environmental Protection Agency under contract NV-03174-20 and the Swedish Climate and Clean Air Research program (SCAC) as well as partly by the Swedish National Space Board (NORD-SLCP, grant agreement ID: 94/16) and the EU Horizon 2020 project Integrated Arctic Observing System (INTAROS, grant agreement ID: 727890). Work on ACE-FTS analysis was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC). Julia Schmale received funding from the Swiss National Science Foundation (project no. 200021_188478). Duncan Watson-Parris received funding from NERC projects NE/P013406/1 (A-CURE) and NE/S005390/1 (ACRUISE) as well as funding from the European Union's Horizon 2020 research and innovation program iMIRACLI under Marie Skłodowska-Curie grant agreement no. 860100. LATMOS has been supported by the EU iCUPE (Integrating and Comprehensive Understanding on Polar Environments) project (grant agreement no. 689443) under the European Network for Observing our Changing Planet (ERA-Planet), as well as access to IDRIS HPC resources (GENCI allocation A009017141) and the IPSL mesoscale computing center (CICLAD: Calcul Intensif pour le CLimat, l’Atmosphère et la Dynamique) for model simulations. Naga Oshima was supported by the Japan Society for the Promotion of Science KAKENHI (grant nos. JP18H03363, JP18H05292, and JP21H03582), the Environment Research and Technology Development Fund (grant nos. JPMEERF20202003 and JPMEERF20205001) of the Environmental Restoration and Conservation Agency of Japan, the Arctic Challenge for Sustainability II (ArCS II) under program grant no. JPMXD1420318865, and a grant for the Global Environmental Research Coordination System from the Ministry of the Environment, Japan (MLIT1753). The research with GISS-E2.1 has been supported by the Aarhus University Interdisciplinary Centre for Climate Change (iClimate) OH fund (no. 2020-0162731), the FREYA project funded by the Nordic Council of Ministers (grant agreement nos. MST-227-00036 and MFVM-2019-13476), and the EVAM-SLCF funded by the Danish Environmental Agency (grant agreement no. MST-112-00298). Jesper Christensen (for DEHM model) received funding from the Danish Environmental Protection Agency (DANCEA funds for Environmental Support to the Arctic Region project; grant no. 2019-7975). Maria Sand has been supported by the Research Council of Norway (grant 315195, ACCEPT). Publisher's note: the article processing charges for this publication were not paid by a Russian or Belarusian institution.

Review statement

This paper was edited by Frank Dentener and reviewed by two anonymous referees.