# Evaluating modelled tropospheric columns of CH4, CO, and O3 in the Arctic using ground-based Fourier transform infrared (FTIR) measurements

Victoria A Flood, Kimberly Strong, Cynthia H Whaley, Kaley A Walker, Thomas Blumenstock, James W Hannigan, Johan Mellqvist, Justus Notholt, Mathias Palm, Amelie N Röhling, et al.

► To cite this version:

Victoria A Flood, Kimberly Strong, Cynthia H Whaley, Kaley A Walker, Thomas Blumenstock, et al. Evaluating modelled tropospheric columns of CH4, CO, and O3 in the Arctic using ground-based Fourier transform infrared (FTIR) measurements. Atmospheric Chemistry and Physics, 2024, 24, pp.1079-1118. 10.5194/acp-24-1079-2024. insu-04171742v2

HAL Id: insu-04171742

https://insu.hal.science/insu-04171742v2

Submitted on 26 Jan 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_1.png" /></figure>

<figure><img src="image_2.png" /></figure>

Evaluating modelled tropospheric columns of CH$_{4}$, CO,
and O$_{3}$ in the Arctic using ground-based Fourier
transform infrared (FTIR) measurements

Victoria A. Flood¹, Kimberly Strong¹, Cynthia H. Whaley², Kaley A. Walker¹, Thomas Blumenstock³,
James W. Hannigan⁴, Johan Mellqvist⁵, Justus Notholt⁶, Mathias Palm⁶, Amelie N. Röhling³,
Stephen Arnold⁷, Stephen Beagley⁸, Rong-You Chien⁹, Jesper Christensen¹⁰, Makoto Deushi¹¹,
Srdjan Dobricic¹², Xinyi Dong⁹, Joshua S. Fu⁹,¹³, Michael Gauss¹⁴, Wanmin Gong⁸, Joakim Langner¹⁵,
Kathy S. Law¹⁶, Louis Marelle¹⁶, Tatsuo Onishi¹⁶, Naga Oshima¹¹, David A. Plummer²,
Luca Pozzoli¹⁷,¹², Jean-Christophe Raut¹⁶, Manu A. Thomas¹⁷, Svetlana Tsyro¹⁴, and Steven Turnock¹⁸,⁷

¹Department of Physics, University of Toronto, Toronto, ON, Canada

²Canadian Centre for Climate Modelling and Analysis, Environment and Climate Change Canada,
Victoria, BC, Canada

³Institute of Meteorology and Climate Research (IMK-ASF), Karlsruhe Institute of Technology,
Karlsruhe, Germany

⁴Atmospheric Chemistry, Observations & Modeling, National Center for Atmospheric Research,
Boulder, CO, USA

⁵Earth and Space Sciences, Chalmers University of Technology, Gothenburg, Sweden

⁶Institute of Environmental Physics, University of Bremen, Bremen, Germany

⁷Institute for Climate and Atmospheric Science, School of Earth and Environment,
University of Leeds, Leeds, UK

⁸Air Quality Research Division, Environment and Climate Change Canada, Toronto, ON, Canada

⁹Department of Civil and Environmental Engineering, University of Tennessee, Knoxville, TN, USA

¹⁰Department of Environmental Science/Interdisciplinary Centre for Climate Change, Aarhus University,
Frederiksborgvej 399, Roskilde, Denmark

¹¹Meteorological Research Institute, Japan Meteorological Agency, Tsukuba, Japan

¹²Joint Research Centre, European Commission, Ispra, Italy

¹³Computational Earth Science Group, Oak Ridge National Laboratory, Oak Ridge, TN, USA

¹⁴Division for Climate Modelling and Air Pollution, Norwegian Meteorological Institute, Oslo, Norway

¹⁵Swedish Meteorological and Hydrological Institute, Norrköping, Sweden

¹⁶UVSQ, CNRS, LATMOS, Sorbonne Université, Paris, France

¹⁷FINCONS SPA, Via Torri Bianche 10, Vimercate, Italy

¹⁸Met Office Hadley Centre, Exeter, UK

**Correspondence:** Victoria A. Flood (v.flood@utoronto.com)

Received: 30 May 2023 – Discussion started: 3 July 2023

Revised: 27 October 2023 – Accepted: 20 November 2023 – Published: 24 January 2024

**Abstract.** This study evaluates tropospheric columns of methane, carbon monoxide, and ozone in the Arctic simulated by 11 models. The Arctic is warming at nearly 4 times the global average rate, and with changing emissions in and near the region, it is important to understand Arctic atmospheric composition and how it is changing. Both measurements and modelling of air pollution in the Arctic are difficult, making model validation with local measurements valuable. Evaluations are performed using data from five high-latitude ground-based Fourier transform infrared (FTIR) spectrometers in the Network for the Detection of Atmospheric Composition Change (NDACC). The models were selected as part of the 2021 Arctic Monitoring and Assessment Programme
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

Evaluating modelled tropospheric columns of CH$_{4}$, CO, and O$_{3}$ in the Arctic using ground-based Fourier transform infrared (FTIR) measurements

Victoria A. Flood¹, Kimberly Strong¹, Cynthia H. Whaley², Kaley A. Walker¹, Thomas Blumenstock³,
James W. Hannigan⁴, Johan Mellqvist⁵, Justus Notholt⁶, Mathias Palm⁶, Amelie N. Röhling³,
Stephen Arnold⁷, Stephen Beagley⁸, Rong-You Chien⁹, Jesper Christensen¹⁰, Makoto Deushi¹¹,
Srdjan Dobricic¹², Xinyi Dong⁹, Joshua S. Fu⁹,¹³, Michael Gauss¹⁴, Wanmin Gong⁸, Joakim Langner¹⁵,
Kathy S. Law¹⁶, Louis Marelle¹⁶, Tatsuo Onishi¹⁶, Naga Oshima¹¹, David A. Plummer²,
Luca Pozzoli¹⁷,¹², Jean-Christophe Raut¹⁶, Manu A. Thomas¹⁷, Svetlana Tsyro¹⁴, and Steven Turnock¹⁸,⁷

¹Department of Physics, University of Toronto, Toronto, ON, Canada

²Canadian Centre for Climate Modelling and Analysis, Environment and Climate Change Canada,
Victoria, BC, Canada

³Institute of Meteorology and Climate Research (IMK-ASF), Karlsruhe Institute of Technology,
Karlsruhe, Germany

⁴Atmospheric Chemistry, Observations & Modeling, National Center for Atmospheric Research,
Boulder, CO, USA

⁵Earth and Space Sciences, Chalmers University of Technology, Gothenburg, Sweden

⁶Institute of Environmental Physics, University of Bremen, Bremen, Germany

⁷Institute for Climate and Atmospheric Science, School of Earth and Environment,
University of Leeds, Leeds, UK

⁸Air Quality Research Division, Environment and Climate Change Canada, Toronto, ON, Canada

⁹Department of Civil and Environmental Engineering, University of Tennessee, Knoxville, TN, USA

¹⁰Department of Environmental Science/Interdisciplinary Centre for Climate Change, Aarhus University,
Frederiksborgvej 399, Roskilde, Denmark

¹¹Meteorological Research Institute, Japan Meteorological Agency, Tsukuba, Japan

¹²Joint Research Centre, European Commission, Ispra, Italy

¹³Computational Earth Science Group, Oak Ridge National Laboratory, Oak Ridge, TN, USA

¹⁴Division for Climate Modelling and Air Pollution, Norwegian Meteorological Institute, Oslo, Norway

¹⁵Swedish Meteorological and Hydrological Institute, Norrköping, Sweden

¹⁶UVSQ, CNRS, LATMOS, Sorbonne Université, Paris, France

¹⁷FINCONS SPA, Via Torri Bianche 10, Vimercate, Italy

¹⁸Met Office Hadley Centre, Exeter, UK

**Correspondence:** Victoria A. Flood (v.flood@utoronto.com)

Received: 30 May 2023 – Discussion started: 3 July 2023

Revised: 27 October 2023 – Accepted: 20 November 2023 – Published: 24 January 2024

**Abstract.** This study evaluates tropospheric columns of methane, carbon monoxide, and ozone in the Arctic simulated by 11 models. The Arctic is warming at nearly 4 times the global average rate, and with changing emissions in and near the region, it is important to understand Arctic atmospheric composition and how it is changing. Both measurements and modelling of air pollution in the Arctic are difficult, making model validation with local measurements valuable. Evaluations are performed using data from five high-latitude ground-based Fourier transform infrared (FTIR) spectrometers in the Network for the Detection of Atmospheric Composition Change (NDACC). The models were selected as part of the 2021 Arctic Monitoring and Assessment Programme
(AMAP) report on short-lived climate forcers. This work augments the model-measurement comparisons presented in that report by including a new data source: column-integrated FTIR measurements, whose spatial and temporal footprint is more representative of the free troposphere than in situ and satellite measurements. Mixing ratios of trace gases are modelled at 3-hourly intervals by CESM, CMAM, DEHM, EMEP MSC-W, GEM-MACH, GEOS-Chem, MATCH, MATCH-SALSA, MRI-ESM2, UKESM1, and WRF-Chem for the years 2008, 2009, 2014, and 2015. The comparisons focus on the troposphere (0–7 km partial columns) at Eureka, Canada; Thule, Greenland; Ny Ålesund, Norway; Kiruna, Sweden; and Harestua, Norway. Overall, the models are biased low in the tropospheric column, on average by –9.7% for CH₄, –21% for CO, and –18% for O₃. Results for CH₄ are relatively consistent across the 4 years, whereas CO has a maximum negative bias in the spring and minimum in the summer and O₃ has a maximum difference centered around the summer. The average differences for the models are within the FTIR uncertainties for approximately 15% of the model-location comparisons.

# 1 Introduction

Short-lived climate forcers (SLCFs) are a group of greenhouse gases and air pollutants with lifetimes less than 2 decades (IPCC, 2021). These include methane (CH₄), ozone (O₃), black carbon, halocarbons, sulfate, nitrate, and organic aerosols. The Intergovernmental Panel on Climate Change (IPCC) reports that in addition to radiative forcing, SLCFs have been found to have negative impacts on air quality, ecosystems, and human health. Due to their relatively short lifetimes, SLCFs are generally reflective of emission rates, meaning that mitigation can result in near-term impacts. Understanding the influences of SLCFs on the future climate will aid in policies and mitigation strategies to stay on track with the Paris Accord and its subsequent amendments. Reductions in SLCFs can be particularly beneficial in the Arctic because models have demonstrated a strong climate response in this region to local and remote forcing by SLCFs (Stohl et al., 2015).

The Arctic Monitoring and Assessment Programme (AMAP) was created by the Arctic Council to provide science-based analysis of Arctic pollution and climate change. AMAP has provided reports on SLCF impacts on the Arctic dating back to 2008. The 2021 AMAP SLCF assessment report assesses the impacts of black carbon, CH₄, O₃, and sulfate aerosols on the air quality, climate, and human health in the Arctic region (AMAP, 2021). A key difference from previous AMAP reports is the emphasis on air quality and human health. In addition to these SLCFs, the analysis includes the SLCF precursor gases carbon monoxide (CO), nitric oxide (NO), and nitrogen dioxide (NO₂). The report compares the output from 18 models with various historical measurements, including satellite, aircraft, ship, and in situ datasets. These observations are used to assess what processes need to be revised in the models and how these shortcomings impact the further application of the models, such as for climate and health predictions. Other chapters explore emissions, measurement advances, trends, climate air quality impacts, health ecosystem impacts, and next steps. A prominent theme in this report is the severity of change

happening in the Arctic. This includes the amplification of the pace of change in physical drivers, such as temperature and snow cover, and the frequency of extreme events, such as wildfires and incidents of rapid sea-ice loss. These factors contribute to ecosystem disruption, directly affecting local Arctic communities, in addition to having global repercussions. SLCF reductions are motivated by the near-term (20–30 years) benefits and by the goal of slowing the warming of the Arctic climate, which results in more wildfires and permafrost melt and, in turn, an increase in SLCF emissions and precursor gases (AMAP, 2021). The projections in this report provide guidance, objectives, and cautions for potential reduction implementation scenarios (AMAP, 2021). This study builds upon the model-measurement comparisons presented in the 2021 AMAP SLCF assessment report using an additional Arctic dataset that was not included in the original report.

The atmospheric measurements used to evaluate the models presented in this paper were made by Fourier transform infrared (FTIR) spectrometers that are contributing members of the Network for Detection of Atmospheric Composition Change (NDACC). NDACC has over 70 stations around the globe, collecting high-quality atmospheric composition measurements with ground-based, remote-sensing instruments (De Mazière et al., 2018). The network’s objective is to create a long-term database for various studies such as atmospheric trends, assessing links between climate and air quality, and as a resource for other atmospheric investigations such as satellite validation and model development. Atmospheric vertical profiles and trace gas columns are retrieved from high-resolution FTIR spectrometers that record solar spectra featuring characteristic atmospheric absorption lines. Five of the 28 NDACC FTIR stations are located at latitudes north of 60°N; for the purpose of this study, these will all be referred to as Arctic sites. The five sites are Eureka, Canada; Ny Ålesund, Norway; Thule, Greenland; Kiruna, Sweden; and Harestua, Norway. These high-latitude NDACC FTIR instruments provide a valuable set of long-term, measurements of multiple species of interest in the Arctic. Compared to surface in situ or satellite observations, the column-integrated
(AMAP) report on short-lived climate forcers. This work augments the model-measurement comparisons presented in that report by including a new data source: column-integrated FTIR measurements, whose spatial and temporal footprint is more representative of the free troposphere than in situ and satellite measurements. Mixing ratios of trace gases are modelled at 3-hourly intervals by CESM, CMAM, DEHM, EMEP MSC-W, GEM-MACH, GEOS-Chem, MATCH, MATCH-SALSA, MRI-ESM2, UKESM1, and WRF-Chem for the years 2008, 2009, 2014, and 2015. The comparisons focus on the troposphere (0–7 km partial columns) at Eureka, Canada; Thule, Greenland; Ny Ålesund, Norway; Kiruna, Sweden; and Harestua, Norway. Overall, the models are biased low in the tropospheric column, on average by –9.7% for CH$_4$, –21% for CO, and –18% for O$_3$. Results for CH$_4$ are relatively consistent across the 4 years, whereas CO has a maximum negative bias in the spring and minimum in the summer and O$_3$ has a maximum difference centered around the summer. The average differences for the models are within the FTIR uncertainties for approximately 15% of the model-location comparisons.

# 1 Introduction

Short-lived climate forcers (SLCFs) are a group of greenhouse gases and air pollutants with lifetimes less than 2 decades (IPCC, 2021). These include methane (CH$_4$), ozone (O$_3$), black carbon, halocarbons, sulfate, nitrate, and organic aerosols. The Intergovernmental Panel on Climate Change (IPCC) reports that in addition to radiative forcing, SLCFs have been found to have negative impacts on air quality, ecosystems, and human health. Due to their relatively short lifetimes, SLCFs are generally reflective of emission rates, meaning that mitigation can result in near-term impacts. Understanding the influences of SLCFs on the future climate will aid in policies and mitigation strategies to stay on track with the Paris Accord and its subsequent amendments. Reductions in SLCFs can be particularly beneficial in the Arctic because models have demonstrated a strong climate response in this region to local and remote forcing by SLCFs (Stohl et al., 2015).

The Arctic Monitoring and Assessment Programme (AMAP) was created by the Arctic Council to provide science-based analysis of Arctic pollution and climate change. AMAP has provided reports on SLCF impacts on the Arctic dating back to 2008. The 2021 AMAP SLCF assessment report assesses the impacts of black carbon, CH$_4$, O$_3$, and sulfate aerosols on the air quality, climate, and human health in the Arctic region (AMAP, 2021). A key difference from previous AMAP reports is the emphasis on air quality and human health. In addition to these SLCFs, the analysis includes the SLCF precursor gases carbon monoxide (CO), nitric oxide (NO), and nitrogen dioxide (NO$_2$). The report compares the output from 18 models with various historical measurements, including satellite, aircraft, ship, and in situ datasets. These observations are used to assess what processes need to be revised in the models and how these shortcomings impact the further application of the models, such as for climate and health predictions. Other chapters explore emissions, measurement advances, trends, climate air quality impacts, health ecosystem impacts, and next steps. A prominent theme in this report is the severity of change

happening in the Arctic. This includes the amplification of the pace of change in physical drivers, such as temperature and snow cover, and the frequency of extreme events, such as wildfires and incidents of rapid sea-ice loss. These factors contribute to ecosystem disruption, directly affecting local Arctic communities, in addition to having global repercussions. SLCF reductions are motivated by the near-term (20–30 years) benefits and by the goal of slowing the warming of the Arctic climate, which results in more wildfires and permafrost melt and, in turn, an increase in SLCF emissions and precursor gases (AMAP, 2021). The projections in this report provide guidance, objectives, and cautions for potential reduction implementation scenarios (AMAP, 2021). This study builds upon the model-measurement comparisons presented in the 2021 AMAP SLCF assessment report using an additional Arctic dataset that was not included in the original report.

The atmospheric measurements used to evaluate the models presented in this paper were made by Fourier transform infrared (FTIR) spectrometers that are contributing members of the Network for Detection of Atmospheric Composition Change (NDACC). NDACC has over 70 stations around the globe, collecting high-quality atmospheric composition measurements with ground-based, remote-sensing instruments (De Mazière et al., 2018). The network’s objective is to create a long-term database for various studies such as atmospheric trends, assessing links between climate and air quality, and as a resource for other atmospheric investigations such as satellite validation and model development. Atmospheric vertical profiles and trace gas columns are retrieved from high-resolution FTIR spectrometers that record solar spectra featuring characteristic atmospheric absorption lines. Five of the 28 NDACC FTIR stations are located at latitudes north of 60°N; for the purpose of this study, these will all be referred to as Arctic sites. The five sites are Eureka, Canada; Ny Ålesund, Norway; Thule, Greenland; Kiruna, Sweden; and Harestua, Norway. These high-latitude NDACC FTIR instruments provide a valuable set of long-term, measurements of multiple species of interest in the Arctic. Compared to surface in situ or satellite observations, the column-integrated
FTIR measurements have a spatial and temporal footprint that is more representative of the free troposphere. Performing model-measurement comparisons with partial column data supports and thus complements the assessments presented in the 2021 AMAP report. Previous studies have used FTIR data to examine model biases in the Arctic (e.g., Wespes et al., 2012; Zhou et al., 2019; Mahieu et al., 2021).

Measurements in the Arctic are difficult due to the harsh environment, remote locations, and high operating costs, resulting in a scarcity of monitoring stations and a limited representation of atmospheric vertical information. Using measurements to evaluate model simulations of the Arctic is important because the latter are used to project future changes in the Arctic, a region that is sensitive to climate change, warming at a rate 3 to 4 times the global average (Bush and Lemmen, 2019; Ballinger et al., 2020; AMAP, 2021; IPCC, 2021; Rantanen et al., 2022). These factors have led to initiatives like the AMAP SLCF assessment and the POLAR-CAT (Polar Study using Aircraft, Remote Sensing, Surface Measurements and Models, of Climate, Chemistry, Aerosols and Transport) Model Intercomparison Project (POLMIP), which, in part, aim to assess model performance in the Arctic region. POLMIP examined 11 atmospheric models in relation to a variety of Arctic observations taken as part of the International Polar Year in 2008 (Emmons et al., 2015). AMAP and POLMIP, in addition to the subsequent complementary publications (i.e., Wespes et al., 2012; Emmons et al., 2015; Monks et al., 2015; Whaley et al., 2022, 2023), provide a valuable point of reference for the modelling of CH₄, CO, and O₃ in the Arctic, which is explored in this paper. This allows for the findings presented here to be appraised relative to results from the same models compared to other instruments, with differing temporal frequency and altitude ranges (i.e., Whaley et al., 2022, 2023), with different simulations and Arctic FTIR measurements (i.e., Wespes et al., 2015), and to generally assess the similarities and differences that arise within Arctic SLCF modelling.

This project examines simulations from 11 models that were run for the 2021 AMAP SLCF assessment report, to assess the agreement between modelled trace gas concentrations and ground-based retrievals from high-latitude FTIR spectrometers. Specifically, this paper presents comparisons of CH₄, CO, and O₃ partial columns (from 0 to 7 km) for the years 2008, 2009, 2014, and 2015. The models examined are chemical transport and climate models: CESM, CMAM, DEHM, EMEP MSC-W, GEM-MACH, GEOS-CHEM, MATCH, MATCH-SALSA, MRI-ESM2, UKESM1, and WRF-CHEM. The objective is to utilize the high-quality, long-term Arctic FTIR datasets to assess how well the models perform. The remainder of this paper is organized as follows: Sect. 2 provides a description of the datasets used, Sect. 3 describes the analysis methodology, Sect. 4 examines the results and compares them with similar studies, and Sect. 5 presents the summary and conclusions.

# 2 Datasets

## 2.1 FTIR spectroscopy and retrievals

The FTIR measurement sites included in this study are summarized in Table 1, and the data are publicly available in the NDACC data repository (https://www-air.larc.nasa.gov/missions/ndacc/data.html, last access: 14 February 2023). These instruments require sunlight and a clear sight to the sun to make measurements, and so the high-latitude datasets are limited to the sunlit portion of the year at each location. To ensure high data quality and consistency between sites, NDACC has several specialized instrument and theme groups; the instruments used here are part of the Infrared Working Group (IRWG). The 10 standard gases reported by sites participating in the IRWG are C₂H₆, CH₄, CO, ClONO₂, HCl, HCN, HF, HNO₃, N₂O, and O₃, while several other gases are retrieved as research data products, including C₂H₂, CH₃OH, H₂CO, HCOOH, and OCS. The FTIR measurements cycle through a series of optical filters covering different spectral regions between approximately 650 and 4500 cm⁻¹ for the retrieval of multiple atmospheric gases. Atmospheric trace gas profiles and columns are retrieved with the SFIT4 algorithm, using optimal estimation to iteratively adjust an a priori profile to match a modelled spectrum to the measured spectrum within a defined convergence criterion (Rodgers, 2000; IRWG, 2020). The a priori profile information for the modelled spectra is provided by 40-year-average profiles from the Whole Atmosphere Community Climate Model (WACCM) (Marsh et al., 2013), with spectroscopic absorption parameters from the HITRAN 2008 line list (Rothman et al., 2009) and daily pressure and temperature profiles from the U.S. National Centers for Environmental Prediction (NCEP) (Kalnay et al., 1996). All sites included in this paper use SFIT4, except Kiruna, which uses a comparable retrieval code called PROFFIT, which has been shown to agree well with SFIT2 (which preceded SFIT4) (Hase et al., 2004). Primary references and further details of the sites are presented in Table 1.

The NDACC FTIR data files include the volume mixing ratio (VMR) in parts per million (ppm) and total columns and partial columns in molecules per centimeter squared (molec.cm⁻²). Other variables include altitude, date and time of observation, pressure, the a priori vertical profile, the averaging kernel (AVK) matrix, and retrieval uncertainties, both systematic and random. The random uncertainties are determined from the temperature, solar zenith angle, and measurement noise from the signal-to-noise ratio. Systematic uncertainties are determined from temperature and line parameters such as line strength and width.

The averaging kernel matrix represents the relationship between the retrieved state and the true atmospheric state at each altitude layer, and the sensitivity of a retrieval is calculated by taking the sum of the rows of the averaging kernel. This indicates how much of the information is coming
FTIR measurements have a spatial and temporal footprint that is more representative of the free troposphere. Performing model-measurement comparisons with partial column data supports and thus complements the assessments presented in the 2021 AMAP report. Previous studies have used FTIR data to examine model biases in the Arctic (e.g., Wespes et al., 2012; Zhou et al., 2019; Mahieu et al., 2021).

Measurements in the Arctic are difficult due to the harsh environment, remote locations, and high operating costs, resulting in a scarcity of monitoring stations and a limited representation of atmospheric vertical information. Using measurements to evaluate model simulations of the Arctic is important because the latter are used to project future changes in the Arctic, a region that is sensitive to climate change, warming at a rate 3 to 4 times the global average (Bush and Lemmen, 2019; Ballinger et al., 2020; AMAP, 2021; IPCC, 2021; Rantanen et al., 2022). These factors have led to initiatives like the AMAP SLCF assessment and the POLAR-CAT (Polar Study using Aircraft, Remote Sensing, Surface Measurements and Models, of Climate, Chemistry, Aerosols and Transport) Model Intercomparison Project (POLMIP), which, in part, aim to assess model performance in the Arctic region. POLMIP examined 11 atmospheric models in relation to a variety of Arctic observations taken as part of the International Polar Year in 2008 (Emmons et al., 2015). AMAP and POLMIP, in addition to the subsequent complementary publications (i.e., Wespes et al., 2012; Emmons et al., 2015; Monks et al., 2015; Whaley et al., 2022, 2023), provide a valuable point of reference for the modelling of CH₄, CO, and O₃ in the Arctic, which is explored in this paper. This allows for the findings presented here to be appraised relative to results from the same models compared to other instruments, with differing temporal frequency and altitude ranges (i.e., Whaley et al., 2022, 2023), with different simulations and Arctic FTIR measurements (i.e., Wespes et al., 2015), and to generally assess the similarities and differences that arise within Arctic SLCF modelling.

This project examines simulations from 11 models that were run for the 2021 AMAP SLCF assessment report, to assess the agreement between modelled trace gas concentrations and ground-based retrievals from high-latitude FTIR spectrometers. Specifically, this paper presents comparisons of CH₄, CO, and O₃ partial columns (from 0 to 7 km) for the years 2008, 2009, 2014, and 2015. The models examined are chemical transport and climate models: CESM, CMAM, DEHM, EMEP MSC-W, GEM-MACH, GEOS-CHEM, MATCH, MATCH-SALSA, MRI-ESM2, UKESM1, and WRF-CHEM. The objective is to utilize the high-quality, long-term Arctic FTIR datasets to assess how well the models perform. The remainder of this paper is organized as follows: Sect. 2 provides a description of the datasets used, Sect. 3 describes the analysis methodology, Sect. 4 examines the results and compares them with similar studies, and Sect. 5 presents the summary and conclusions.

# 2 Datasets

## 2.1 FTIR spectroscopy and retrievals

The FTIR measurement sites included in this study are summarized in Table 1, and the data are publicly available in the NDACC data repository (https://www-air.larc.nasa.gov/missions/ndacc/data.html, last access: 14 February 2023). These instruments require sunlight and a clear sight to the sun to make measurements, and so the high-latitude datasets are limited to the sunlit portion of the year at each location. To ensure high data quality and consistency between sites, NDACC has several specialized instrument and theme groups; the instruments used here are part of the Infrared Working Group (IRWG). The 10 standard gases reported by sites participating in the IRWG are C₂H₆, CH₄, CO, ClONO₂, HCl, HCN, HF, HNO₃, N₂O, and O₃, while several other gases are retrieved as research data products, including C₂H₂, CH₃OH, H₂CO, HCOOH, and OCS. The FTIR measurements cycle through a series of optical filters covering different spectral regions between approximately 650 and 4500 cm⁻¹ for the retrieval of multiple atmospheric gases. Atmospheric trace gas profiles and columns are retrieved with the SFIT4 algorithm, using optimal estimation to iteratively adjust an a priori profile to match a modelled spectrum to the measured spectrum within a defined convergence criterion (Rodgers, 2000; IRWG, 2020). The a priori profile information for the modelled spectra is provided by 40-year-average profiles from the Whole Atmosphere Community Climate Model (WACCM) (Marsh et al., 2013), with spectroscopic absorption parameters from the HITRAN 2008 line list (Rothman et al., 2009) and daily pressure and temperature profiles from the U.S. National Centers for Environmental Prediction (NCEP) (Kalnay et al., 1996). All sites included in this paper use SFIT4, except Kiruna, which uses a comparable retrieval code called PROFFIT, which has been shown to agree well with SFIT2 (which preceded SFIT4) (Hase et al., 2004). Primary references and further details of the sites are presented in Table 1.

The NDACC FTIR data files include the volume mixing ratio (VMR) in parts per million (ppm) and total columns and partial columns in molecules per centimeter squared (molec.cm⁻²). Other variables include altitude, date and time of observation, pressure, the a priori vertical profile, the averaging kernel (AVK) matrix, and retrieval uncertainties, both systematic and random. The random uncertainties are determined from the temperature, solar zenith angle, and measurement noise from the signal-to-noise ratio. Systematic uncertainties are determined from temperature and line parameters such as line strength and width.

The averaging kernel matrix represents the relationship between the retrieved state and the true atmospheric state at each altitude layer, and the sensitivity of a retrieval is calculated by taking the sum of the rows of the averaging kernel. This indicates how much of the information is coming
<table><caption>**Table 1.** Summary of NDACC FTIR sites used in this study.</caption><thead><tr><td>Site</td><td>Location</td><td>Key references</td><td>Operations</td></tr></thead><tbody><tr><td>Eureka, Canada</td><td>80.05° N, 86.42° W; 610 m a.s.l.</td><td>Batchelor et al. (2009)</td><td>Late February to mid-October Since 2006</td></tr><tr><td>Ny Ålesund, Norway</td><td>78.92° N, 11.93° E; 15 m a.s.l.</td><td>Notholt et al. (1997a, b, 2000)</td><td>Mid-March to September Since 1992</td></tr><tr><td>Thule, Greenland</td><td>76.53° N, 68.74° W; 225 m a.s.l.</td><td>Hannigan et al. (2009)</td><td>March to October Since 1999</td></tr><tr><td>Kiruna, Sweden</td><td>67.84° N, 20.41° E; 419 m a.s.l.</td><td>Blumenstock et al. (1997, 2009)</td><td>Mid-January to November Since 1996</td></tr><tr><td>Harestua, Norway</td><td>60.2° N, 10.8° E; 596 m a.s.l.</td><td>Galle et al. (1999)</td><td>All year Since 1994</td></tr></tbody></table>

from the a priori profile and how much comes from the measurement itself (Rodgers, 2000; Vigouroux et al., 2009). The degrees of freedom for signal (DOFS) are calculated by taking the trace of the averaging kernel; this indicates the number of independent pieces of information coming from each retrieval or, inversely, the number of components not constrained by the a priori profile information. The random and systematic FTIR partial column uncertainties are calculated using the error covariance matrices, following the method outlined in Vigouroux et al. (2009). The square root of the associated error is taken, and this is scaled to a percent uncertainty using the corresponding partial column sum. The mean systematic and random percent errors are added in quadrature to get the overall mean percent uncertainty for the species. The number of measurements, mean DOFS, and mean percent uncertainty in the 0–7 km partial columns of CH$_4$, CO, and O$_3$ for 2008, 2009, 2014, and 2015, for each station, are listed in Table 2. The mean partial column (0–7 and 7–20 km) and total column averaging kernels for CH$_4$, CO, and O$_3$ for 2008, 2009, 2014, and 2015, are shown in Fig. 1. The lowest-level difference between Kiruna and the other locations results from the use of a stronger constraint for the lowest level with the PROFFIT retrieval; however, retrieval error and noise indicate that the agreement between the averaging kernels (AVKs) is reasonable (Hase et al., 2004). The DOFS and averaging kernels are indicators of the vertical information within a retrieval. Figure 1 shows that the mean partial column averaging kernels for 0–7 and 7–20 km are distinguishable, with maxima at different altitudes. The mean total column averaging kernels for all three species appear smooth around 1.0, which indicates that contributions from all altitudes have similar weights in the total column. By altitude, the sensitivity of each species is > 0.5 in the partial column examined (not shown), meaning that more than half of the retrieved profile information comes from the measurement (Vigouroux et al., 2009). The average DOFS vary by species and station, given the reduced column height of 0–7 km; some of the values are less than 1, meaning the retrieval

is somewhat constrained by the a priori profile. However, it should be noted that the comparisons presented in this paper account for the vertical sensitivity of the FTIR measurements by smoothing the model data with the averaging kernels. This process is described in Sect. 3.

## 2.2 Atmospheric models

The models used in this study provide three-dimensional VMR fields on 3-hourly intervals for 2008, 2009, 2014, and 2015. These 4 years were selected for the 2021 AMAP SLCF assessment; 2008 and 2009 were previously evaluated in the 2015 AMAP report, and 2014 and 2015 were added to include more recent results from years for which Arctic measurements were available at the time (AMAP, 2021). The gases CH$_4$, CO, and O$_3$ were chosen for this study as model output for these species was available at 3-hourly intervals, and the FTIR measurements have good sensitivity for them throughout the 0–7 km with the FTIR, as discussed in the previous section. Note that not every model has provided all three gases; there are 3 which have CH$_4$, 9 with CO, and 11 with O$_3$ (see Table 3). The model simulations are the same as those discussed in Whaley et al. (2022, 2023) and the 2021 AMAP SLCF report; however, the analyses there were performed with the monthly mean output, while the analysis here is with the 3-hourly output, all of which is available at http://crd-data-donnees-rdc.ec.gc.ca/CCCMAP/products/AMAP/ (last access: 14 February 2023). While more models participated in the AMAP SLCF assessment (18 total) and other species were simulated, these were not included in the current study because either the models did not have 3-hourly outputs or the FTIR retrievals had insufficient tropospheric sensitivity (e.g., NO$_2$).

This set of models is a mix of Earth system models, chemical transport models, global transport models, and chemistry climate models. The models all used the same set of anthropogenic emissions from ECLIPSE v6b (Evaluating the Climate and Air Quality Impacts of Short-Lived Pollutants) by
<table><caption>**Table 1.** Summary of NDACC FTIR sites used in this study.</caption><thead><tr><td>Site</td><td>Location</td><td>Key references</td><td>Operations</td></tr></thead><tbody><tr><td>Eureka, Canada</td><td>80.05° N, 86.42° W; 610 m a.s.l.</td><td>Batchelor et al. (2009)</td><td>Late February to mid-October Since 2006</td></tr><tr><td>Ny Ålesund, Norway</td><td>78.92° N, 11.93° E; 15 m a.s.l.</td><td>Notholt et al. (1997a, b, 2000)</td><td>Mid-March to September Since 1992</td></tr><tr><td>Thule, Greenland</td><td>76.53° N, 68.74° W; 225 m a.s.l.</td><td>Hannigan et al. (2009)</td><td>March to October Since 1999</td></tr><tr><td>Kiruna, Sweden</td><td>67.84° N, 20.41° E; 419 m a.s.l.</td><td>Blumenstock et al. (1997, 2009)</td><td>Mid-January to November Since 1996</td></tr><tr><td>Harestua, Norway</td><td>60.2° N, 10.8° E; 596 m a.s.l.</td><td>Galle et al. (1999)</td><td>All year Since 1994</td></tr></tbody></table>

from the a priori profile and how much comes from the measurement itself (Rodgers, 2000; Vigouroux et al., 2009). The degrees of freedom for signal (DOFS) are calculated by taking the trace of the averaging kernel; this indicates the number of independent pieces of information coming from each retrieval or, inversely, the number of components not constrained by the a priori profile information. The random and systematic FTIR partial column uncertainties are calculated using the error covariance matrices, following the method outlined in Vigouroux et al. (2009). The square root of the associated error is taken, and this is scaled to a percent uncertainty using the corresponding partial column sum. The mean systematic and random percent errors are added in quadrature to get the overall mean percent uncertainty for the species. The number of measurements, mean DOFS, and mean percent uncertainty in the 0–7 km partial columns of CH$_4$, CO, and O$_3$ for 2008, 2009, 2014, and 2015, for each station, are listed in Table 2. The mean partial column (0–7 and 7–20 km) and total column averaging kernels for CH$_4$, CO, and O$_3$ for 2008, 2009, 2014, and 2015, are shown in Fig. 1. The lowest-level difference between Kiruna and the other locations results from the use of a stronger constraint for the lowest level with the PROFFIT retrieval; however, retrieval error and noise indicate that the agreement between the averaging kernels (AVKs) is reasonable (Hase et al., 2004). The DOFS and averaging kernels are indicators of the vertical information within a retrieval. Figure 1 shows that the mean partial column averaging kernels for 0–7 and 7–20 km are distinguishable, with maxima at different altitudes. The mean total column averaging kernels for all three species appear smooth around 1.0, which indicates that contributions from all altitudes have similar weights in the total column. By altitude, the sensitivity of each species is > 0.5 in the partial column examined (not shown), meaning that more than half of the retrieved profile information comes from the measurement (Vigouroux et al., 2009). The average DOFS vary by species and station, given the reduced column height of 0–7 km; some of the values are less than 1, meaning the retrieval

is somewhat constrained by the a priori profile. However, it should be noted that the comparisons presented in this paper account for the vertical sensitivity of the FTIR measurements by smoothing the model data with the averaging kernels. This process is described in Sect. 3.

## 2.2 Atmospheric models

The models used in this study provide three-dimensional VMR fields on 3-hourly intervals for 2008, 2009, 2014, and 2015. These 4 years were selected for the 2021 AMAP SLCF assessment; 2008 and 2009 were previously evaluated in the 2015 AMAP report, and 2014 and 2015 were added to include more recent results from years for which Arctic measurements were available at the time (AMAP, 2021). The gases CH$_4$, CO, and O$_3$ were chosen for this study as model output for these species was available at 3-hourly intervals, and the FTIR measurements have good sensitivity for them throughout the 0–7 km with the FTIR, as discussed in the previous section. Note that not every model has provided all three gases; there are 3 which have CH$_4$, 9 with CO, and 11 with O$_3$ (see Table 3). The model simulations are the same as those discussed in Whaley et al. (2022, 2023) and the 2021 AMAP SLCF report; however, the analyses there were performed with the monthly mean output, while the analysis here is with the 3-hourly output, all of which is available at http://crd-data-donnees-rdc.ec.gc.ca/CCCMAP/products/AMAP/ (last access: 14 February 2023). While more models participated in the AMAP SLCF assessment (18 total) and other species were simulated, these were not included in the current study because either the models did not have 3-hourly outputs or the FTIR retrievals had insufficient tropospheric sensitivity (e.g., NO$_2$).

This set of models is a mix of Earth system models, chemical transport models, global transport models, and chemistry climate models. The models all used the same set of anthropogenic emissions from ECLIPSE v6b (Evaluating the Climate and Air Quality Impacts of Short-Lived Pollutants) by
<table><caption>**Table 2.** Summary of FTIR measurement statistics.</caption><thead><tr><th rowspan="2">Site</th><th colspan="3">Number of measurements<br>(2008, 2009, 2014, 2015)</th><th colspan="3">Mean DOFS<br>(0–7 km)</th><th colspan="3">Mean percent uncertainty<br>(0–7 km)</th></tr><tr><th>CH<sub>4</sub></th><th>CO</th><th>O<sub>3</sub></th><th>CH<sub>4</sub></th><th>CO</th><th>O<sub>3</sub></th><th>CH<sub>4</sub></th><th>CO</th><th>O<sub>3</sub></th></tr></thead><tbody><tr><td>Eureka</td><td>754</td><td>736</td><td>684</td><td>0.84</td><td>1.1</td><td>0.80</td><td>4.6</td><td>3.9</td><td>8.2</td></tr><tr><td>Ny Ålesund</td><td>205</td><td>128</td><td>121</td><td>0.81</td><td>1.3</td><td>0.79</td><td>11.5</td><td>7.7</td><td>4.9</td></tr><tr><td>Thule</td><td>406</td><td>459</td><td>474</td><td>0.78</td><td>1.6</td><td>1.2</td><td>5.7</td><td>5.4</td><td>3.9</td></tr><tr><td>Kiruna</td><td>397</td><td>299</td><td>322</td><td>0.96</td><td>1.6</td><td>0.86</td><td>3.6</td><td>6.4</td><td>7.2</td></tr><tr><td>Harestua</td><td>151 (no 2008)</td><td>No CO</td><td>169 (no 2008)</td><td>0.78</td><td>No CO</td><td>1.12</td><td>5.2</td><td>No CO</td><td>4.1</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Mean 0–7 km partial column averaging kernels (lines with circle markers), mean 7–20 km partial column averaging kernels (dashed lines), and mean total column averaging kernels (solid lines), all in units of (molec. cm⁻² (molec. cm⁻²)⁻¹), by altitude, for (a) CH₄, (b) CO, and (c) O₃. Means are for 2008, 2009, 2014, and 2015 for all five FTIR sites except Harestua (no 2008 data).</figcaption></figure>

the IIASA GAINS (International Institute for Applied Systems Analysis – Greenhouse gas – Air pollution Interactions and Synergies) model (Amann et al., 2011; Klimont et al., 2017; Höglund-Isaksson et al., 2020). However, the models differ in their use of biogenic and volcanic emissions, tropospheric gas-phase chemistry complexity, and vertical and horizontal grids. Four of the 11 models simulate the stratosphere fully, one (GEOS-Chem) uses a simplified linearized stratospheric chemistry, one (GEM-MACH) only simulates the troposphere, and the rest use prescribed climatologies at the stratospheric boundary (Whaley et al., 2022). In total, 9 of the 11 models examined use the Global Fire Emissions Database (GFED; van der Werf et al., 2017) or GFED-based (CMIP6) forest fire emissions, and 9 of the 11 exclusively use ECLIPSEv6b for agricultural waste burning. A summary of the models is presented in Table 3, including which gases are examined in this study, their resolution, and to what degree stratospheric chemistry is considered. It should be noted that the CH₄ concentrations in these models have been prescribed (Whaley et al., 2022). The prescribed concentrations are input at the bottom model layer, and all come from the same dataset (Prather et al., 2012; Olivié et al., 2021), but the resulting CH₄ partial columns differ based on the processes within each model. For a full description of the models, see

Appendix A of Whaley et al. (2022) and the references in Table 3.

## 3 Methods

As mentioned, the models provided 3-hourly VMRs on model-specific pressure levels and latitude–longitude grids. The process of aligning the model output to FTIR data is described by the flow chart in Fig. 2.

This procedure modifies the model output to correspond to an FTIR measurement, making the resulting partial columns equivalent for further comparison. The date and time and volume-mixing-ratio profiles from the model output are extracted from the grid point that is closest to the FTIR location. The FTIR measurements are matched with the 3-hourly model measurement closest in time (± < 1.5 h); this is done to minimize the time difference between the two points, such that no measurement is greater than 1.5 h from a modelled output. If more than one FTIR measurement coincides with a model output (i.e., multiple measurements are within 1.5 h of the same model time), the FTIR measurements are averaged. After the model outputs are matched to the FTIR measurements, they are interpolated onto the pressure grid of the FTIR profile. Then, the model VMR profile is smoothed us-
<table><caption>**Table 2.** Summary of FTIR measurement statistics.</caption><thead><tr><th rowspan="2">Site</th><th colspan="3">Number of measurements<br>(2008, 2009, 2014, 2015)</th><th colspan="3">Mean DOFS<br>(0–7 km)</th><th colspan="3">Mean percent uncertainty<br>(0–7 km)</th></tr><tr><th>CH<sub>4</sub></th><th>CO</th><th>O<sub>3</sub></th><th>CH<sub>4</sub></th><th>CO</th><th>O<sub>3</sub></th><th>CH<sub>4</sub></th><th>CO</th><th>O<sub>3</sub></th></tr></thead><tbody><tr><td>Eureka</td><td>754</td><td>736</td><td>684</td><td>0.84</td><td>1.1</td><td>0.80</td><td>4.6</td><td>3.9</td><td>8.2</td></tr><tr><td>Ny Ålesund</td><td>205</td><td>128</td><td>121</td><td>0.81</td><td>1.3</td><td>0.79</td><td>11.5</td><td>7.7</td><td>4.9</td></tr><tr><td>Thule</td><td>406</td><td>459</td><td>474</td><td>0.78</td><td>1.6</td><td>1.2</td><td>5.7</td><td>5.4</td><td>3.9</td></tr><tr><td>Kiruna</td><td>397</td><td>299</td><td>322</td><td>0.96</td><td>1.6</td><td>0.86</td><td>3.6</td><td>6.4</td><td>7.2</td></tr><tr><td>Harestua</td><td>151 (no 2008)</td><td>No CO</td><td>169 (no 2008)</td><td>0.78</td><td>No CO</td><td>1.12</td><td>5.2</td><td>No CO</td><td>4.1</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Mean 0–7 km partial column averaging kernels (lines with circle markers), mean 7–20 km partial column averaging kernels (dashed lines), and mean total column averaging kernels (solid lines), all in units of (molec. cm⁻² (molec. cm⁻²)⁻¹), by altitude, for (a) CH₄, (b) CO, and (c) O₃. Means are for 2008, 2009, 2014, and 2015 for all five FTIR sites except Harestua (no 2008 data).</figcaption></figure>

the IIASA GAINS (International Institute for Applied Systems Analysis – Greenhouse gas – Air pollution Interactions and Synergies) model (Amann et al., 2011; Klimont et al., 2017; Höglund-Isaksson et al., 2020). However, the models differ in their use of biogenic and volcanic emissions, tropospheric gas-phase chemistry complexity, and vertical and horizontal grids. Four of the 11 models simulate the stratosphere fully, one (GEOS-Chem) uses a simplified linearized stratospheric chemistry, one (GEM-MACH) only simulates the troposphere, and the rest use prescribed climatologies at the stratospheric boundary (Whaley et al., 2022). In total, 9 of the 11 models examined use the Global Fire Emissions Database (GFED; van der Werf et al., 2017) or GFED-based (CMIP6) forest fire emissions, and 9 of the 11 exclusively use ECLIPSEv6b for agricultural waste burning. A summary of the models is presented in Table 3, including which gases are examined in this study, their resolution, and to what degree stratospheric chemistry is considered. It should be noted that the CH₄ concentrations in these models have been prescribed (Whaley et al., 2022). The prescribed concentrations are input at the bottom model layer, and all come from the same dataset (Prather et al., 2012; Olivié et al., 2021), but the resulting CH₄ partial columns differ based on the processes within each model. For a full description of the models, see

Appendix A of Whaley et al. (2022) and the references in Table 3.

## 3 Methods

As mentioned, the models provided 3-hourly VMRs on model-specific pressure levels and latitude–longitude grids. The process of aligning the model output to FTIR data is described by the flow chart in Fig. 2.

This procedure modifies the model output to correspond to an FTIR measurement, making the resulting partial columns equivalent for further comparison. The date and time and volume-mixing-ratio profiles from the model output are extracted from the grid point that is closest to the FTIR location. The FTIR measurements are matched with the 3-hourly model measurement closest in time (± < 1.5 h); this is done to minimize the time difference between the two points, such that no measurement is greater than 1.5 h from a modelled output. If more than one FTIR measurement coincides with a model output (i.e., multiple measurements are within 1.5 h of the same model time), the FTIR measurements are averaged. After the model outputs are matched to the FTIR measurements, they are interpolated onto the pressure grid of the FTIR profile. Then, the model VMR profile is smoothed us-
<table><caption>**Table 3.** Summary of the models used in this study.</caption><thead><tr><th>Model</th><th>Three-hourly outputs</th><th>Primary reference</th><th>Horizontal resolution and scale</th><th>Stratospheric chemistry</th></tr></thead><tbody><tr><td>CESM<br>Community Earth System Model</td><td>CO, O<sub>3</sub></td><td>Liu et al. (2016),<br>Danabasoglu et al. (2020)</td><td>1.9° × 2.5°<br>global</td><td>comprehensive</td></tr><tr><td>CMAM<br>Canadian Middle Atmosphere Model</td><td>CH<sub>4</sub>, CO, O<sub>3</sub></td><td>Jonsson et al. (2004),<br>Scinocca et al. (2008)</td><td>3.75° × 3.75°<br>global</td><td>comprehensive</td></tr><tr><td>DEHM<br>Danish Eulerian Hemispheric Model</td><td>O<sub>3</sub></td><td>Christensen (1997),<br>Brandt et al. (2012),<br>Massling et al. (2015)</td><td>50 km<br>polar stereographic</td><td>none</td></tr><tr><td>EMEP MSC-W<br>European Monitoring and Evaluation<br>System-Meteorological Synthesizing<br>Center - West</td><td>CO, O<sub>3</sub></td><td>Simpson et al. (2012, 2019)</td><td>0.5° × 0.5°<br>global</td><td>prescribed</td></tr><tr><td>GEM-MACH<br>Global Environmental Multiscale Model<br>- Modelling Air Quality and Chemistry</td><td>CO, O<sub>3</sub><br>(only 2015)</td><td>Gong et al. (2015),<br>Makar et al. (2015a, b),<br>Moran et al. (2018)</td><td>15 km<br>Arctic regional</td><td>none</td></tr><tr><td>GEOS-CHEM<br>Goddard Earth Observing System - Chemistry</td><td>CH<sub>4</sub>, CO, O<sub>3</sub></td><td>Bey et al. (2001)</td><td>2° × 2.5°<br>global</td><td>simplified</td></tr><tr><td>MATCH<br>Multi-Scale Atmospheric Transport Chemistry</td><td>CO, O<sub>3</sub></td><td>Robertson et al. (1999)</td><td>0.75° rotated<br>lat-long regional</td><td>prescribed</td></tr><tr><td>MATCH-SALSA<br>Multi-Scale Atmospheric Transport Chemistry -<br>Sectional Aerosol Module for Large<br>Scale Applications</td><td>CO, O<sub>3</sub></td><td>Robertson et al. (1999),<br>Andersson et al. (2007),<br>Kokkola et al. (2008)</td><td>0.75°<br>rotated lat-long<br>regional</td><td>prescribed</td></tr><tr><td>MRI-ESM2<br>Meteorological Research Institute –<br>Earth System Model Version 2</td><td>CH<sub>4</sub>, CO, O<sub>3</sub></td><td>Kawai et al. (2019),<br>Yukimoto et al. (2019),<br>Oshima et al. (2020)</td><td>chemistry: 280 km<br>general: 120 km<br>global</td><td>comprehensive</td></tr><tr><td>UKESM1<br>U.K. Earth System Model Version 1</td><td>O<sub>3</sub></td><td>Kuhlbrodt et al. (2018),<br>Williams et al. (2018),<br>Sellar et al. (2019)</td><td>140 km<br>global</td><td>comprehensive</td></tr><tr><td>WRF-CHEM<br>Weather Research and Forecasting<br>Model with Chemistry</td><td>CO, O<sub>3</sub><br>(only<br>2014–2015)</td><td>Marelle et al. (2017, 2018)</td><td>100 km<br>regional Arctic</td><td>prescribed</td></tr></tbody></table>

ing the respective FTIR measurement's averaging kernel and a priori profile. The purpose of smoothing the model data with the FTIR averaging kernel is to adjust the model to the vertical sensitivity of the FTIR measurement (Rodgers and Connor, 2003). The calculation for the smoothing is shown in Eq. (1), where $x_a$ is the FTIR a priori VMR vertical profile, A is the VMR averaging kernel matrix from the corresponding FTIR measurement, and $x_{\text{model}}$ is the modelled VMR vertical profile:

$$x_{\text{smooth}} = x_a + A \times [x_{\text{model}} - x_a]. \quad (1)$$

The model VMR profile is then transformed to a layer profile in units of molecules per centimeter squared using the ratio between the VMR and layer partial column (in molecules per centimeter squared) in the retrieved FTIR profile as the con-

version factor. At this point, the model output has the same altitude grid and units as the FTIR retrieval, which allows for partial columns to be summed. Partial columns from 0 to 7 km were calculated given AMAP's focus on SLCFs in the troposphere, with the cap at 7 km chosen to limit any stratospheric influence. Note that “0 km” is used as proxy for the minimum altitude, but this varies, based on location, with the altitude of each instrument listed in Table 1. The partial column examined here (0–7 km) encompasses 11 vertical layers for all sites, except Ny Ålesund, which has an additional (12th) layer given the lower altitude of its location (see Table 1).

To compare the model and FTIR partial columns, a model-measurement percent difference ($\Delta_i$) is calculated, as defined by Eq. (2) for a single model-measurement pair (i),
<table><caption>**Table 3.** Summary of the models used in this study.</caption><thead><tr><th>Model</th><th>Three-hourly outputs</th><th>Primary reference</th><th>Horizontal resolution and scale</th><th>Stratospheric chemistry</th></tr></thead><tbody><tr><td>CESM<br>Community Earth System Model</td><td>CO, O<sub>3</sub></td><td>Liu et al. (2016),<br>Danabasoglu et al. (2020)</td><td>1.9° × 2.5°<br>global</td><td>comprehensive</td></tr><tr><td>CMAM<br>Canadian Middle Atmosphere Model</td><td>CH<sub>4</sub>, CO, O<sub>3</sub></td><td>Jonsson et al. (2004),<br>Scinocca et al. (2008)</td><td>3.75° × 3.75°<br>global</td><td>comprehensive</td></tr><tr><td>DEHM<br>Danish Eulerian Hemispheric Model</td><td>O<sub>3</sub></td><td>Christensen (1997),<br>Brandt et al. (2012),<br>Massling et al. (2015)</td><td>50 km<br>polar stereographic</td><td>none</td></tr><tr><td>EMEP MSC-W<br>European Monitoring and Evaluation<br>System-Meteorological Synthesizing<br>Center - West</td><td>CO, O<sub>3</sub></td><td>Simpson et al. (2012, 2019)</td><td>0.5° × 0.5°<br>global</td><td>prescribed</td></tr><tr><td>GEM-MACH<br>Global Environmental Multiscale Model<br>- Modelling Air Quality and Chemistry</td><td>CO, O<sub>3</sub><br>(only 2015)</td><td>Gong et al. (2015),<br>Makar et al. (2015a, b),<br>Moran et al. (2018)</td><td>15 km<br>Arctic regional</td><td>none</td></tr><tr><td>GEOS-CHEM<br>Goddard Earth Observing System - Chemistry</td><td>CH<sub>4</sub>, CO, O<sub>3</sub></td><td>Bey et al. (2001)</td><td>2° × 2.5°<br>global</td><td>simplified</td></tr><tr><td>MATCH<br>Multi-Scale Atmospheric Transport Chemistry</td><td>CO, O<sub>3</sub></td><td>Robertson et al. (1999)</td><td>0.75° rotated<br>lat-long regional</td><td>prescribed</td></tr><tr><td>MATCH-SALSA<br>Multi-Scale Atmospheric Transport Chemistry -<br>Sectional Aerosol Module for Large<br>Scale Applications</td><td>CO, O<sub>3</sub></td><td>Robertson et al. (1999),<br>Andersson et al. (2007),<br>Kokkola et al. (2008)</td><td>0.75°<br>rotated lat-long<br>regional</td><td>prescribed</td></tr><tr><td>MRI-ESM2<br>Meteorological Research Institute –<br>Earth System Model Version 2</td><td>CH<sub>4</sub>, CO, O<sub>3</sub></td><td>Kawai et al. (2019),<br>Yukimoto et al. (2019),<br>Oshima et al. (2020)</td><td>chemistry: 280 km<br>general: 120 km<br>global</td><td>comprehensive</td></tr><tr><td>UKESM1<br>U.K. Earth System Model Version 1</td><td>O<sub>3</sub></td><td>Kuhlbrodt et al. (2018),<br>Williams et al. (2018),<br>Sellar et al. (2019)</td><td>140 km<br>global</td><td>comprehensive</td></tr><tr><td>WRF-CHEM<br>Weather Research and Forecasting<br>Model with Chemistry</td><td>CO, O<sub>3</sub><br>(only<br>2014–2015)</td><td>Marelle et al. (2017, 2018)</td><td>100 km<br>regional Arctic</td><td>prescribed</td></tr></tbody></table>

ing the respective FTIR measurement's averaging kernel and a priori profile. The purpose of smoothing the model data with the FTIR averaging kernel is to adjust the model to the vertical sensitivity of the FTIR measurement (Rodgers and Connor, 2003). The calculation for the smoothing is shown in Eq. (1), where $x_a$ is the FTIR a priori VMR vertical profile, A is the VMR averaging kernel matrix from the corresponding FTIR measurement, and $x_{\text{model}}$ is the modelled VMR vertical profile:

$$x_{\text{smooth}} = x_a + A \times [x_{\text{model}} - x_a]. \quad (1)$$

The model VMR profile is then transformed to a layer profile in units of molecules per centimeter squared using the ratio between the VMR and layer partial column (in molecules per centimeter squared) in the retrieved FTIR profile as the con-

version factor. At this point, the model output has the same altitude grid and units as the FTIR retrieval, which allows for partial columns to be summed. Partial columns from 0 to 7 km were calculated given AMAP's focus on SLCFs in the troposphere, with the cap at 7 km chosen to limit any stratospheric influence. Note that “0 km” is used as proxy for the minimum altitude, but this varies, based on location, with the altitude of each instrument listed in Table 1. The partial column examined here (0–7 km) encompasses 11 vertical layers for all sites, except Ny Ålesund, which has an additional (12th) layer given the lower altitude of its location (see Table 1).

To compare the model and FTIR partial columns, a model-measurement percent difference ($\Delta_i$) is calculated, as defined by Eq. (2) for a single model-measurement pair (i),
<figure><img src="image_4.png" /><figcaption>Figure 2. Flow chart depicting the process of matching model output to FTIR data.</figcaption></figure>

where PC<sub>M,i</sub> and PC<sub>F,i</sub> are the 0–7 km partial columns for the model and FTIR, respectively:

$$ \Delta_i = \left( \frac{\mathrm{PC}_{M,i} - \mathrm{PC}_{F,i}}{\mathrm{PC}_{F,i}} \right) \times 100. \qquad (2) $$

A regression line is fit to the raw scatter plot data of the model output versus FTIR measurements using all the available data points, where each plot includes the equation of this line and the correlation coefficient, R<sup>2</sup>. The normalized root mean square error (NRMSE), given by Eq. (3), is presented for each model and location, where N is the total number of model-measurement pairs (Kärnä and Baptista, 2016). The root mean square error is normalized to the standard deviation of the FTIR data (σ<sub>F</sub>) used in the respective analysis:

$$ \mathrm{NRMSE} = \frac{1}{\sigma_F} \sqrt{\left[ \sum_{i=1}^{N} (\mathrm{PC}_{M,i} - \mathrm{PC}_{F,i})^2 \right]}. \qquad (3) $$

In addition to evaluating the models using every available FTIR data point in the analysis years, the monthly mean annual cycles are also presented. The monthly mean partial columns (PC<sub>F,monthly,j</sub>) are calculated by taking the mean of every measurement in a given month (j), where N<sub>j</sub> is the number of points included in the month for all years considered. The monthly model mean partial columns (PC<sub>M,monthly,j</sub>) are made in the same manner, using only the smoothed partial columns that have a corresponding matching FTIR measurement, as defined above. Equation (4) outlines the calculation of a monthly mean partial column for month j for (a) the FTIRs (PC<sub>F,monthly,j</sub>) and (b) the models

(PC<sub>M,monthly,j</sub>):

$$ \mathrm{PC}_{M,\mathrm{monthly},j} = \frac{1}{N_j} \sum_{i=1}^{N_j} \mathrm{PC}_{M,i}, \qquad (4a) $$

$$ \mathrm{PC}_{F,\mathrm{monthly},j} = \frac{1}{N_j} \sum_{i=1}^{N_j} \mathrm{PC}_{F,i}. \qquad (4b) $$

The model-measurement monthly mean percent difference (Δ<sub>monthly,j</sub>), shown by Eq. (5), follows the same process as the monthly mean partial column and is the mean value from Eq. (2) for each month (j) across the years, where the error bars on the monthly mean plots represent the standard deviation of this mean:

$$ \Delta_{\mathrm{monthly},j} = \frac{1}{N_j} \sum_{i=1}^{N_j} \Delta_i. \qquad (5) $$

The mean of these monthly mean differences is used to calculate the overall mean percent difference (Δ<sub>O</sub>) for each model, sometimes referred to as model bias, where N<sub>months</sub> is the number of measurement months in a calendar year at that location (see Table 1), and the uncertainty given is the standard deviation of this mean:

$$ \Delta_O = \frac{1}{N_{\mathrm{months}}} \sum_{j=1}^{N_{\mathrm{months}}} \Delta_{\mathrm{monthly},j}. \qquad (6) $$

Finally, the monthly multi-model mean (MMM) partial column for month j (PC<sub>MMM,monthly,j</sub>) is calculated by taking the mean PC<sub>M,monthly,j</sub> for all models, at a given location, calculated with Eq. (4b), and the MMM monthly mean difference (Δ<sub>MMM,monthly,j</sub>) is the mean of Δ<sub>monthly,j</sub> for all models, at a given location, calculated with Eq. (5). The overall percent difference in the MMM measurement (Δ<sub>O,MMM</sub>) is given by Eq. (7):

$$ \Delta_{O,MMM} = \frac{1}{N_{\mathrm{months}}} \sum_{j=1}^{N_{\mathrm{months}}} \Delta_{MMM,\mathrm{monthly},j}. \qquad (7) $$

These steps are taken to establish the modelled seasonal cycles and quantify the differences between the models and measurements, by month and season. Further, assessing the MMM by month allows for a general overview of when and where models diverge from measurements and can help suggest shortcomings in the models. There are not enough measurements per day to evaluate a diurnal cycle, although it is expected to be small in the Arctic, and there are not enough years available in the 3-hourly dataset used here to examine long-term trends.

When discussing FTIR uncertainty, this refers to the mean uncertainty per gas and station, as listed in Table 2. When discussing the mean difference between the model and measurements, this refers to the overall mean difference (Δ<sub>O</sub>) as described by Eq. (6). In Sects. 4 and 5, these two parameters are used to assess model performance: if Δ<sub>O</sub> is within measurement (FTIR) uncertainty, the model can be considered in general agreement with the FTIR; if Δ<sub>O</sub>± the standard deviation of the mean is within the measurement uncertainty,
<figure><img src="image_5.png" /><figcaption>Figure 2. Flow chart depicting the process of matching model output to FTIR data.</figcaption></figure>

where PC<sub>M,i</sub> and PC<sub>F,i</sub> are the 0–7 km partial columns for the model and FTIR, respectively:

$$ \Delta_i = \left( \frac{\mathrm{PC}_{M,i} - \mathrm{PC}_{F,i}}{\mathrm{PC}_{F,i}} \right) \times 100. \qquad (2) $$

A regression line is fit to the raw scatter plot data of the model output versus FTIR measurements using all the available data points, where each plot includes the equation of this line and the correlation coefficient, R<sup>2</sup>. The normalized root mean square error (NRMSE), given by Eq. (3), is presented for each model and location, where N is the total number of model-measurement pairs (Kärnä and Baptista, 2016). The root mean square error is normalized to the standard deviation of the FTIR data (σ<sub>F</sub>) used in the respective analysis:

$$ \mathrm{NRMSE} = \frac{1}{\sigma_F} \sqrt{\left[ \sum_{i=1}^{N} (\mathrm{PC}_{M,i} - \mathrm{PC}_{F,i})^2 \right]}. \qquad (3) $$

In addition to evaluating the models using every available FTIR data point in the analysis years, the monthly mean annual cycles are also presented. The monthly mean partial columns (PC<sub>F,monthly,j</sub>) are calculated by taking the mean of every measurement in a given month (j), where N<sub>j</sub> is the number of points included in the month for all years considered. The monthly model mean partial columns (PC<sub>M,monthly,j</sub>) are made in the same manner, using only the smoothed partial columns that have a corresponding matching FTIR measurement, as defined above. Equation (4) outlines the calculation of a monthly mean partial column for month j for (a) the FTIRs (PC<sub>F,monthly,j</sub>) and (b) the models

(PC<sub>M,monthly,j</sub>):

$$ \mathrm{PC}_{M,\mathrm{monthly},j} = \frac{1}{N_j} \sum_{i=1}^{N_j} \mathrm{PC}_{M,i}, \qquad (4a) $$

$$ \mathrm{PC}_{F,\mathrm{monthly},j} = \frac{1}{N_j} \sum_{i=1}^{N_j} \mathrm{PC}_{F,i}. \qquad (4b) $$

The model-measurement monthly mean percent difference (Δ<sub>monthly,j</sub>), shown by Eq. (5), follows the same process as the monthly mean partial column and is the mean value from Eq. (2) for each month (j) across the years, where the error bars on the monthly mean plots represent the standard deviation of this mean:

$$ \Delta_{\mathrm{monthly},j} = \frac{1}{N_j} \sum_{i=1}^{N_j} \Delta_i. \qquad (5) $$

The mean of these monthly mean differences is used to calculate the overall mean percent difference (Δ<sub>O</sub>) for each model, sometimes referred to as model bias, where N<sub>months</sub> is the number of measurement months in a calendar year at that location (see Table 1), and the uncertainty given is the standard deviation of this mean:

$$ \Delta_O = \frac{1}{N_{\mathrm{months}}} \sum_{j=1}^{N_{\mathrm{months}}} \Delta_{\mathrm{monthly},j}. \qquad (6) $$

Finally, the monthly multi-model mean (MMM) partial column for month j (PC<sub>MMM,monthly,j</sub>) is calculated by taking the mean PC<sub>M,monthly,j</sub> for all models, at a given location, calculated with Eq. (4b), and the MMM monthly mean difference (Δ<sub>MMM,monthly,j</sub>) is the mean of Δ<sub>monthly,j</sub> for all models, at a given location, calculated with Eq. (5). The overall percent difference in the MMM measurement (Δ<sub>O,MMM</sub>) is given by Eq. (7):

$$ \Delta_{O,MMM} = \frac{1}{N_{\mathrm{months}}} \sum_{j=1}^{N_{\mathrm{months}}} \Delta_{MMM,\mathrm{monthly},j}. \qquad (7) $$

These steps are taken to establish the modelled seasonal cycles and quantify the differences between the models and measurements, by month and season. Further, assessing the MMM by month allows for a general overview of when and where models diverge from measurements and can help suggest shortcomings in the models. There are not enough measurements per day to evaluate a diurnal cycle, although it is expected to be small in the Arctic, and there are not enough years available in the 3-hourly dataset used here to examine long-term trends.

When discussing FTIR uncertainty, this refers to the mean uncertainty per gas and station, as listed in Table 2. When discussing the mean difference between the model and measurements, this refers to the overall mean difference (Δ<sub>O</sub>) as described by Eq. (6). In Sects. 4 and 5, these two parameters are used to assess model performance: if Δ<sub>O</sub> is within measurement (FTIR) uncertainty, the model can be considered in general agreement with the FTIR; if Δ<sub>O</sub>± the standard deviation of the mean is within the measurement uncertainty,
<table><caption>**Table 4.** The multi-model mean percent difference (Δ<sub>O,MMM</sub>) for each species at each location, including the overall average percent difference for each species and the standard deviation of the mean.</caption><thead><tr><th>Gas</th><th>Location</th><th>MMM percent difference</th></tr></thead><tbody><tr><td rowspan="5">CH<sub>4</sub></td><td>Eureka</td><td>-9.9 ± 0.7</td></tr><tr><td>Ny Ålesund</td><td>-10.2 ± 0.7</td></tr><tr><td>Thule</td><td>-7.5 ± 2.0</td></tr><tr><td>Kiruna</td><td>-11.6 ± 0.5</td></tr><tr><td>Harestua</td><td>-9.2 ± 1.4</td></tr><tr><td></td><td>Average</td><td>-9.7</td></tr><tr><td rowspan="5">CO</td><td>Eureka</td><td>-17.6 ± 5.6</td></tr><tr><td>Ny Ålesund</td><td>-16.7 ± 7.9</td></tr><tr><td>Thule</td><td>-24.4 ± 6.5</td></tr><tr><td>Kiruna</td><td>-23.7 ± 5.2</td></tr><tr><td>Average</td><td>-20.6</td></tr><tr><td rowspan="5">O<sub>3</sub></td><td>Eureka</td><td>-20.1 ± 10.2</td></tr><tr><td>Ny Ålesund</td><td>-28.5 ± 8.3</td></tr><tr><td>Thule</td><td>-17.6 ± 9.8</td></tr><tr><td>Kiruna</td><td>-14.6 ± 8.7</td></tr><tr><td>Harestua</td><td>-9.6 ± 9.5</td></tr><tr><td></td><td>Average</td><td>-18.1</td></tr></tbody></table>

then the model is sometimes in agreement with the measurements; and if the uncertainty and Δ<sub>O</sub> do not overlap, then the model and measurements do not agree.

# 4 Results and discussion

This section presents the analyses described above, for CH<sub>4</sub>, CO, and O<sub>3</sub>, and discusses the findings in the context of the 2021 AMAP SLCF assessment report and other related literature. Given the volume of data (three species, five locations, and 11 models), only selected plots are shown in the main text, with the remaining figures provided in Appendices A–C. These include plots for each location, showing the time series of the 0–7 km partial column for each measurement–model pair and the associated model–measurement percent difference, the equivalent plot reduced to monthly mean data (an individualized version of Figs. 3, 5, and 9), and the 0–7 km column of FTIR vs. smoothed model for the remaining locations (analogous to Figs. 4, 8, and 10). Figure 15 provides a summary of the overall differences for each model and location by species, as described by Eq. (6). Table 4 summarizes the overall MMM difference for each species at each location and the overall average for each species. All the comparisons shown are for a 0–7 km partial column, where the model output is smoothed as described by Eq. (1).

## 4.1 CH<sub>4</sub>

CH<sub>4</sub> is a powerful greenhouse gas (GHG), and its emissions are expected to increase in the Arctic due to melting per-

mafrost (IPCC, 2021). CH<sub>4</sub> is also involved in the formation of tropospheric O<sub>3</sub>, which is the third strongest anthropogenic GHG and an air pollutant at the surface. Therefore, it is important for both air quality and climate models to represent CH<sub>4</sub> accurately. The CH<sub>4</sub> plots for Ny Ålesund, Thule, Kiruna, and Harestua are provided in Appendix A, following the same order discussed here for Eureka.

Figure 3 shows the monthly mean 0–7 km partial column time series for the FTIR and models at each location (panels a–e), with the percent difference between the monthly mean model and monthly mean measurement for all locations shown in panel (f). This shows that, apart from a few outliers, the pattern of the seasonal cycle of CH<sub>4</sub> is consistent, although the amplitude is underestimated. The uniformity between the years (see Figs. A1–A5 for full data time series plots) and consistency of the model biases between sites is likely a consequence of CH<sub>4</sub> being prescribed in the models, in addition to the longer lifetime of CH<sub>4</sub>, relative to the other SLCFs. This is also seen in Fig. 4 (and Figs. A11–A14), where the model and FTIR columns are compared, with the line of best fit and R<sup>2</sup> are indicated in the legend.

A summary of the overall mean difference, R<sup>2</sup>, and the normalized root-mean-square error for each location is shown in Fig. 5. Across all three models, Arctic CH<sub>4</sub> is underpredicted compared to the FTIR measurements. The surface in situ CH<sub>4</sub> comparison in Whaley et al. (2022) showed that measured surface CH<sub>4</sub> VMRs are much more variable than the modelled VMRs. However, in the 0–7 km partial columns in this study, CH<sub>4</sub> is well-mixed and more homogeneous, resulting in better agreement between the models and the FTIR measurements. The low bias we find in this study for the Arctic sites is consistent with the global comparisons of these models to satellite measurements in Whaley et al. (2022), which found that some models did not distribute CH<sub>4</sub> with an accurate north–south gradient, resulting in low biases in the Arctic and high biases in lower latitudes. GEOS-Chem does simulate a north–south gradient, which is reflected in the smaller overall model–measurement percent difference, compared to other models, in all locations (note Fig. 6 in Whaley et al., 2022). However, the R<sup>2</sup> of GEOS-Chem vs. FTIR is smaller than that for the other models at some locations (Eureka and Kiruna), which can be attributed to the increase in variability the gradient introduces – including some instances of overestimation. The mean differences for each model across sites are relatively consistent, while the results vary more when comparing R<sup>2</sup> and NRMSE. Particularly, when comparing them between the same model, the R<sup>2</sup> for Ny Ålesund is the lowest and the NRMSE is the highest. The data from Ny Ålesund show less of a seasonal cycle than the other locations, and the FTIR uncertainty for CH<sub>4</sub> at Ny Ålesund is more than twice that of the other sites (see Fig. 15). The larger uncertainty may lead to reduced sensitivity to small changes and increased variability masking seasonal changes, which can contribute to the discrepancy between the models and observations. The mean difference
<table><caption>**Table 4.** The multi-model mean percent difference (Δ<sub>O,MMM</sub>) for each species at each location, including the overall average percent difference for each species and the standard deviation of the mean.</caption><thead><tr><th>Gas</th><th>Location</th><th>MMM percent difference</th></tr></thead><tbody><tr><td rowspan="5">CH<sub>4</sub></td><td>Eureka</td><td>-9.9 ± 0.7</td></tr><tr><td>Ny Ålesund</td><td>-10.2 ± 0.7</td></tr><tr><td>Thule</td><td>-7.5 ± 2.0</td></tr><tr><td>Kiruna</td><td>-11.6 ± 0.5</td></tr><tr><td>Harestua</td><td>-9.2 ± 1.4</td></tr><tr><td colspan="2">Average</td><td>-9.7</td></tr><tr><td rowspan="5">CO</td><td>Eureka</td><td>-17.6 ± 5.6</td></tr><tr><td>Ny Ålesund</td><td>-16.7 ± 7.9</td></tr><tr><td>Thule</td><td>-24.4 ± 6.5</td></tr><tr><td>Kiruna</td><td>-23.7 ± 5.2</td></tr><tr><td>Average</td><td>-20.6</td></tr><tr><td rowspan="5">O<sub>3</sub></td><td>Eureka</td><td>-20.1 ± 10.2</td></tr><tr><td>Ny Ålesund</td><td>-28.5 ± 8.3</td></tr><tr><td>Thule</td><td>-17.6 ± 9.8</td></tr><tr><td>Kiruna</td><td>-14.6 ± 8.7</td></tr><tr><td>Harestua</td><td>-9.6 ± 9.5</td></tr><tr><td colspan="2">Average</td><td>-18.1</td></tr></tbody></table>

then the model is sometimes in agreement with the measurements; and if the uncertainty and Δ<sub>O</sub> do not overlap, then the model and measurements do not agree.

# 4 Results and discussion

This section presents the analyses described above, for CH<sub>4</sub>, CO, and O<sub>3</sub>, and discusses the findings in the context of the 2021 AMAP SLCF assessment report and other related literature. Given the volume of data (three species, five locations, and 11 models), only selected plots are shown in the main text, with the remaining figures provided in Appendices A–C. These include plots for each location, showing the time series of the 0–7 km partial column for each measurement–model pair and the associated model–measurement percent difference, the equivalent plot reduced to monthly mean data (an individualized version of Figs. 3, 5, and 9), and the 0–7 km column of FTIR vs. smoothed model for the remaining locations (analogous to Figs. 4, 8, and 10). Figure 15 provides a summary of the overall differences for each model and location by species, as described by Eq. (6). Table 4 summarizes the overall MMM difference for each species at each location and the overall average for each species. All the comparisons shown are for a 0–7 km partial column, where the model output is smoothed as described by Eq. (1).

## 4.1 CH<sub>4</sub>

CH<sub>4</sub> is a powerful greenhouse gas (GHG), and its emissions are expected to increase in the Arctic due to melting per-

mafrost (IPCC, 2021). CH<sub>4</sub> is also involved in the formation of tropospheric O<sub>3</sub>, which is the third strongest anthropogenic GHG and an air pollutant at the surface. Therefore, it is important for both air quality and climate models to represent CH<sub>4</sub> accurately. The CH<sub>4</sub> plots for Ny Ålesund, Thule, Kiruna, and Harestua are provided in Appendix A, following the same order discussed here for Eureka.

Figure 3 shows the monthly mean 0–7 km partial column time series for the FTIR and models at each location (panels a–e), with the percent difference between the monthly mean model and monthly mean measurement for all locations shown in panel (f). This shows that, apart from a few outliers, the pattern of the seasonal cycle of CH<sub>4</sub> is consistent, although the amplitude is underestimated. The uniformity between the years (see Figs. A1–A5 for full data time series plots) and consistency of the model biases between sites is likely a consequence of CH<sub>4</sub> being prescribed in the models, in addition to the longer lifetime of CH<sub>4</sub>, relative to the other SLCFs. This is also seen in Fig. 4 (and Figs. A11–A14), where the model and FTIR columns are compared, with the line of best fit and R<sup>2</sup> are indicated in the legend.

A summary of the overall mean difference, R<sup>2</sup>, and the normalized root-mean-square error for each location is shown in Fig. 5. Across all three models, Arctic CH<sub>4</sub> is underpredicted compared to the FTIR measurements. The surface in situ CH<sub>4</sub> comparison in Whaley et al. (2022) showed that measured surface CH<sub>4</sub> VMRs are much more variable than the modelled VMRs. However, in the 0–7 km partial columns in this study, CH<sub>4</sub> is well-mixed and more homogeneous, resulting in better agreement between the models and the FTIR measurements. The low bias we find in this study for the Arctic sites is consistent with the global comparisons of these models to satellite measurements in Whaley et al. (2022), which found that some models did not distribute CH<sub>4</sub> with an accurate north–south gradient, resulting in low biases in the Arctic and high biases in lower latitudes. GEOS-Chem does simulate a north–south gradient, which is reflected in the smaller overall model–measurement percent difference, compared to other models, in all locations (note Fig. 6 in Whaley et al., 2022). However, the R<sup>2</sup> of GEOS-Chem vs. FTIR is smaller than that for the other models at some locations (Eureka and Kiruna), which can be attributed to the increase in variability the gradient introduces – including some instances of overestimation. The mean differences for each model across sites are relatively consistent, while the results vary more when comparing R<sup>2</sup> and NRMSE. Particularly, when comparing them between the same model, the R<sup>2</sup> for Ny Ålesund is the lowest and the NRMSE is the highest. The data from Ny Ålesund show less of a seasonal cycle than the other locations, and the FTIR uncertainty for CH<sub>4</sub> at Ny Ålesund is more than twice that of the other sites (see Fig. 15). The larger uncertainty may lead to reduced sensitivity to small changes and increased variability masking seasonal changes, which can contribute to the discrepancy between the models and observations. The mean difference
<figure><img src="image_5.png" /></figure>

**Figure 3.** (a–e) Monthly mean FTIR (black) and smoothed model (colour) 0–7 km partial columns of CH₄ (PCF,monthly, j and PCM,monthly, j, respectively), for each location, shown with the same y axis. Error bars represent the standard deviation of the monthly mean. (f) Mean model–measurement percent difference by month (Δ<sub>monthly, j</sub>) for each model (by colour) and location (by marker). Error bars represent the standard deviation of the monthly mean percent difference.

<figure><img src="image_6.png" /></figure>

**Figure 4.** Smoothed model vs. FTIR 0–7 km partial columns of CH₄ for Eureka, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and R² are noted in the legend. The 1 : 1 line is shown in light grey.

for GEOS-Chem is within the uncertainty in the FTIR mea-
surements for Ny Ålesund and Thule, as is the mean differ-
ence for MRI-ESM2 at Ny Ålesund; none of the other mod-
els are within the FTIR uncertainty at the given location (see
Fig. 15).

Figure 6 shows the multi-model mean (MMM) for each
location and the percent difference compared to the monthly
mean FTIR. The error bars and shading represent the stan-
dard deviation of the mean. The AMAP SLCF assessment
report compares the models with surface CH₄ measurements
and finds that the MMM bias for Arctic CH₄ is +1.3 %
(AMAP, 2021). When comparing them with 0–7 km FTIR
partial columns, the MMM bias ranges from –5 % to –15 %
(Fig. 6f) and unlike the results in the AMAP report, the com-
parisons are not improved by choosing a multi-model mean
because all three models have a negative bias. The FTIR re-
trievals show good sensitivity to tropospheric CH₄ (sensitiv-
ity > 0.5); however, as these column measurements average

out CH₄ biases over the tropospheric column, they are not
expected to exactly match the surface measurement compar-
isons. Furthermore, due to the sharp decrease in CH₄ above
the tropopause (Whaley et al., 2022), a poor representation of
the tropopause height may contribute to the low bias in the
modelled 0–7 km partial columns, as shown from O₃ data in
Whaley et al. (2023). The AMAP report also includes a com-
parison with upper-troposphere/lower-stratosphere (UTLS)
CH₄ VMRs as measured by the ACE-FTS (Atmospheric
Chemistry Experiment – Fourier Transform Spectrometer)
satellite instrument and finds that the models are biased low
by ~100 ppbv in the vicinity of the tropopause (300 hPa;
around ~8–9 km), indicating that the modelled tropopause
may be too low (Whaley et al., 2022). The results found here
are consistent with Whaley et al. (2022), in that the model
simulations of both the lower troposphere (0–7 km partial
columns) and the UTLS are biased low, and models with
north–south CH₄ gradients (here, only GEOS-Chem) have
<figure><img src="image_6.png" /></figure>

**Figure 3.** (a–e) Monthly mean FTIR (black) and smoothed model (colour) 0–7 km partial columns of CH₄ (PCF,monthly, j and PCM,monthly, j, respectively), for each location, shown with the same y axis. Error bars represent the standard deviation of the monthly mean. (f) Mean model–measurement percent difference by month (Δ<sub>monthly, j</sub>) for each model (by colour) and location (by marker). Error bars represent the standard deviation of the monthly mean percent difference.

<figure><img src="image_7.png" /></figure>

**Figure 4.** Smoothed model vs. FTIR 0–7 km partial columns of CH₄ for Eureka, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and R² are noted in the legend. The 1 : 1 line is shown in light grey.

for GEOS-Chem is within the uncertainty in the FTIR mea-
surements for Ny Ålesund and Thule, as is the mean differ-
ence for MRI-ESM2 at Ny Ålesund; none of the other mod-
els are within the FTIR uncertainty at the given location (see
Fig. 15).

Figure 6 shows the multi-model mean (MMM) for each
location and the percent difference compared to the monthly
mean FTIR. The error bars and shading represent the stan-
dard deviation of the mean. The AMAP SLCF assessment
report compares the models with surface CH₄ measurements
and finds that the MMM bias for Arctic CH₄ is +1.3 %
(AMAP, 2021). When comparing them with 0–7 km FTIR
partial columns, the MMM bias ranges from –5 % to –15 %
(Fig. 6f) and unlike the results in the AMAP report, the com-
parisons are not improved by choosing a multi-model mean
because all three models have a negative bias. The FTIR re-
trievals show good sensitivity to tropospheric CH₄ (sensitiv-
ity > 0.5); however, as these column measurements average

out CH₄ biases over the tropospheric column, they are not
expected to exactly match the surface measurement compar-
isons. Furthermore, due to the sharp decrease in CH₄ above
the tropopause (Whaley et al., 2022), a poor representation of
the tropopause height may contribute to the low bias in the
modelled 0–7 km partial columns, as shown from O₃ data in
Whaley et al. (2023). The AMAP report also includes a com-
parison with upper-troposphere/lower-stratosphere (UTLS)
CH₄ VMRs as measured by the ACE-FTS (Atmospheric
Chemistry Experiment – Fourier Transform Spectrometer)
satellite instrument and finds that the models are biased low
by ~100 ppbv in the vicinity of the tropopause (300 hPa;
around ~8–9 km), indicating that the modelled tropopause
may be too low (Whaley et al., 2022). The results found here
are consistent with Whaley et al. (2022), in that the model
simulations of both the lower troposphere (0–7 km partial
columns) and the UTLS are biased low, and models with
north–south CH₄ gradients (here, only GEOS-Chem) have
<figure><img src="image_7.png" /><figcaption>**Figure 5.** By model and location: (a) overall model-measurement mean percent difference for CH₄ 0–7 km partial columns (ΔO), with error bars that represent the standard deviation of the mean, as shown in the legend of Figs. A6–A10; (b) R² as shown in Figs. 4 and A11–A14; (c) normalized root-mean-square error.</figcaption></figure>

smaller biases than those that do not. Generally, the models can represent the temporal variability in the tropospheric column well, although they are biased low in magnitude, outside of the range of the FTIR uncertainty.

## 4.2 CO

Like CH₄, CO is involved in tropospheric O₃ formation in the presence of NOₓ. Thus, in order to properly simulate tropospheric O₃, it is important for models to accurately simulate CO. In the Arctic, CO is used as a tracer for identifying and quantifying influences from biomass burning and lower-latitude anthropogenic emissions (e.g., Fisher et al., 2010; Monks et al., 2015; Viatte et al., 2015; Lutsch et al., 2020).

Nine of the 11 models examined in this study provided 3-hourly outputs for CO; WRF-Chem only has outputs for 2014 and 2015, and GEM-MACH only has data for 2015 (Table 3). Seven of the nine CO models examined use GFED-based fire emissions. The remaining models are EMEP MSC-W, which uses FINN fire emissions, and GEM-MACH, which uses CFFEPS fire emissions (Whaley et al., 2022). Evidence of biomass burning events can be observed in the summer months when examining the CO seasonal cycle with all available measurement points, where there are sporadic increases in the measured CO (Figs. B1–B4). The CO time series data (i.e., Figs. B1–B4 and 7/B5–B8) indicates that the GFED-based models may overestimate CO from biomass burning as their bias shifts positively in the summertime relative to the rest of the time series. This feature is absent for GEM-MACH, which does not have a consistent trend between sites during the summer (although results are only available for 1 year), and for EMEP MSC-W, which shifts more negatively in the summertime. It is well known that the fire emissions inventories vary greatly from each other (AMAP, 2021), causing these differences in model results.

Figure 7 and Figs. B5–B8 show the monthly mean partial columns and percent differences between the models and the FTIR measurements. This allows for an overview of the mean percent difference and how the model biases change over the year. For example, MATCH exhibits a positive shift in bias from the end of summer to the fall in all locations. WRF-Chem is biased low in the spring and summer but agrees better with the observations from August onwards, in contrast to EMEP-MSC-W, which tends to diverge from the measurements in the mid- to late summer. GEM-MACH is the only model that has a positive mean difference in all locations. The year-round difference is likely due to the fact that this model used anthropogenic emissions produced locally for most of its regional domain, instead of the ECLIPSEv6B anthropogenic emissions that all of the other models used, and lateral regional boundary conditions provided from MOZART4 (Model for Ozone and Related Chemical Tracers, version 4) global simulations (Emmons et al., 2010; Gong et al., 2018; AMAP, 2021). Further, Fig. 8 and Figs. B9–B11 show the correlations between the modelled and FTIR partial columns, with the line of best fit and R² indicated in the legend. For many models, the 1 : 1 correlation and Figs. 8 and B9–B11 show that models have better agreement with the FTIR for low CO values and the disparity increases as CO increases; i.e., the line of best fit and 1 : 1 line diverge. The points with the maximum CO VMRs correspond to the FTIR springtime peak in the CO cycle (since wintertime CO measurements are not possible during polar night).

Figure 9 summarizes the overall model-measurement mean percent difference R² and normalized root-mean-square error for all locations. GEM-MACH has a mean percent difference that is within the FTIR uncertainty for Thule and Kiruna; EMEP MSC-W and MATCH are simulated within the mean FTIR uncertainty for Ny Ålesund (see Fig. 15). MATCH-SALSA and MRI-ESM2 exhibit high R² and low percent difference across all locations, relative to the other models’ values, although their columns do not fall within the FTIR uncertainties. GEM-MACH and MATCH have an NRMSE comparable to MATCH-SALSA and MRI-ESM2, despite generally lower R². WRF-Chem shows better agreement with the FTIR measurements from Eureka, where the NRMSE is comparable to CESM, CMAM, and GEOS-Chem. This is likely a result of the increased density of measurement points in August and September, when WRF-Chem exhibits a minimum bias compared to the FTIR data and because the comparison only includes data points from 2014 and 2015. The large negative biases earlier in the year lead to low R² and high NRMSE at all sites. This appears to be linked to negative biases in modelled surface CO over mid-latitude source regions and in the free troposphere compared to MOPITT data, as reported by Whaley et al. (2022). Overall, four model-location pairs have a mean difference within the average FTIR 0–7 km partial column uncertainty (see Ta-
<figure><img src="image_8.png" /><figcaption>**Figure 5.** By model and location: (a) overall model-measurement mean percent difference for CH₄ 0–7 km partial columns (ΔO), with error bars that represent the standard deviation of the mean, as shown in the legend of Figs. A6–A10; (b) R² as shown in Figs. 4 and A11–A14; (c) normalized root-mean-square error.</figcaption></figure>

smaller biases than those that do not. Generally, the models can represent the temporal variability in the tropospheric column well, although they are biased low in magnitude, outside of the range of the FTIR uncertainty.

## 4.2 CO

Like CH₄, CO is involved in tropospheric O₃ formation in the presence of NOₓ. Thus, in order to properly simulate tropospheric O₃, it is important for models to accurately simulate CO. In the Arctic, CO is used as a tracer for identifying and quantifying influences from biomass burning and lower-latitude anthropogenic emissions (e.g., Fisher et al., 2010; Monks et al., 2015; Viatte et al., 2015; Lutsch et al., 2020).

Nine of the 11 models examined in this study provided 3-hourly outputs for CO; WRF-Chem only has outputs for 2014 and 2015, and GEM-MACH only has data for 2015 (Table 3). Seven of the nine CO models examined use GFED-based fire emissions. The remaining models are EMEP MSC-W, which uses FINN fire emissions, and GEM-MACH, which uses CFFEPS fire emissions (Whaley et al., 2022). Evidence of biomass burning events can be observed in the summer months when examining the CO seasonal cycle with all available measurement points, where there are sporadic increases in the measured CO (Figs. B1–B4). The CO time series data (i.e., Figs. B1–B4 and 7/B5–B8) indicates that the GFED-based models may overestimate CO from biomass burning as their bias shifts positively in the summertime relative to the rest of the time series. This feature is absent for GEM-MACH, which does not have a consistent trend between sites during the summer (although results are only available for 1 year), and for EMEP MSC-W, which shifts more negatively in the summertime. It is well known that the fire emissions inventories vary greatly from each other (AMAP, 2021), causing these differences in model results.

Figure 7 and Figs. B5–B8 show the monthly mean partial columns and percent differences between the models and the FTIR measurements. This allows for an overview of the mean percent difference and how the model biases change over the year. For example, MATCH exhibits a positive shift in bias from the end of summer to the fall in all locations. WRF-Chem is biased low in the spring and summer but agrees better with the observations from August onwards, in contrast to EMEP-MSC-W, which tends to diverge from the measurements in the mid- to late summer. GEM-MACH is the only model that has a positive mean difference in all locations. The year-round difference is likely due to the fact that this model used anthropogenic emissions produced locally for most of its regional domain, instead of the ECLIPSEv6B anthropogenic emissions that all of the other models used, and lateral regional boundary conditions provided from MOZART4 (Model for Ozone and Related Chemical Tracers, version 4) global simulations (Emmons et al., 2010; Gong et al., 2018; AMAP, 2021). Further, Fig. 8 and Figs. B9–B11 show the correlations between the modelled and FTIR partial columns, with the line of best fit and R² indicated in the legend. For many models, the 1 : 1 correlation and Figs. 8 and B9–B11 show that models have better agreement with the FTIR for low CO values and the disparity increases as CO increases; i.e., the line of best fit and 1 : 1 line diverge. The points with the maximum CO VMRs correspond to the FTIR springtime peak in the CO cycle (since wintertime CO measurements are not possible during polar night).

Figure 9 summarizes the overall model-measurement mean percent difference R² and normalized root-mean-square error for all locations. GEM-MACH has a mean percent difference that is within the FTIR uncertainty for Thule and Kiruna; EMEP MSC-W and MATCH are simulated within the mean FTIR uncertainty for Ny Ålesund (see Fig. 15). MATCH-SALSA and MRI-ESM2 exhibit high R² and low percent difference across all locations, relative to the other models’ values, although their columns do not fall within the FTIR uncertainties. GEM-MACH and MATCH have an NRMSE comparable to MATCH-SALSA and MRI-ESM2, despite generally lower R². WRF-Chem shows better agreement with the FTIR measurements from Eureka, where the NRMSE is comparable to CESM, CMAM, and GEOS-Chem. This is likely a result of the increased density of measurement points in August and September, when WRF-Chem exhibits a minimum bias compared to the FTIR data and because the comparison only includes data points from 2014 and 2015. The large negative biases earlier in the year lead to low R² and high NRMSE at all sites. This appears to be linked to negative biases in modelled surface CO over mid-latitude source regions and in the free troposphere compared to MOPITT data, as reported by Whaley et al. (2022). Overall, four model-location pairs have a mean difference within the average FTIR 0–7 km partial column uncertainty (see Ta-
<figure><img src="image_8.png" /></figure>

**Figure 6.** (a–e) Monthly mean FTIR (black) and multi-model mean (coloured) 0–7 km partial columns of CH₄ (PCF,monthly,j and PC,MMM,monthly,j, respectively), with error bars and shaded areas representing the standard deviation of the mean. (f) Monthly mean percent difference in the MMM (ΔO,MMM) for all locations.

<figure><img src="image_9.png" /></figure>

**Figure 7.** (a–d) Monthly mean FTIR (black) and smoothed model (colour) 0–7 km partial columns of CO (PCF,monthly,j and PC,M,monthly,j, respectively), for each location, shown with the same y axis. Error bars represent the standard deviation of the monthly mean. (e) Model-measurement mean percent difference by month (Δmonthly,j) for each model (by colour) and location (by marker). Error bars represent standard deviation of the monthly mean percent difference.

ble 2) and when including the standard deviation of the mean difference, an additional 8 pairs out of 36 meet this criterion.

Figure 10 shows the monthly MMM for CO at each location, with the percent difference in the last panel (e). This highlights the general tendency of the models to underpredict tropospheric CO more in the spring than in the summer, which has been observed by other Arctic model-measurement comparison studies. The AMAP SLCF assessment report found that compared to CO from various sur-

face networks, the models had a greater bias than for the other SLFCs examined, underestimating CO in the spring and overestimating CO in the summer (AMAP, 2021). The same pattern was observed when comparing them with MO-PITT (Measurements of Pollution In The Troposphere) satellite CO in the free troposphere, at the 600 hPa level (Whaley et al., 2022). The change from a negative winter-spring bias to a positive summer bias was observed in model comparisons to surface CO measurements at two additional Arctic
<figure><img src="image_9.png" /><figcaption>**Figure 6.** (a–e) Monthly mean FTIR (black) and multi-model mean (coloured) 0–7 km partial columns of CH₄ (PCF,monthly,j and PC,MMM,monthly,j, respectively), with error bars and shaded areas representing the standard deviation of the mean. (f) Monthly mean percent difference in the MMM (ΔO,MMM) for all locations.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 7.** (a–d) Monthly mean FTIR (black) and smoothed model (colour) 0–7 km partial columns of CO (PCF,monthly,j and PC,M,monthly,j, respectively), for each location, shown with the same y axis. Error bars represent the standard deviation of the monthly mean. (e) Model-measurement mean percent difference by month (Δmonthly,j) for each model (by colour) and location (by marker). Error bars represent standard deviation of the monthly mean percent difference.</figcaption></figure>

ble 2) and when including the standard deviation of the mean difference, an additional 8 pairs out of 36 meet this criterion.

Figure 10 shows the monthly MMM for CO at each location, with the percent difference in the last panel (e). This highlights the general tendency of the models to underpredict tropospheric CO more in the spring than in the summer, which has been observed by other Arctic model-measurement comparison studies. The AMAP SLCF assessment report found that compared to CO from various sur-

face networks, the models had a greater bias than for the other SLFCs examined, underestimating CO in the spring and overestimating CO in the summer (AMAP, 2021). The same pattern was observed when comparing them with MO-PITT (Measurements of Pollution In The Troposphere) satellite CO in the free troposphere, at the 600 hPa level (Whaley et al., 2022). The change from a negative winter-spring bias to a positive summer bias was observed in model comparisons to surface CO measurements at two additional Arctic
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Smoothed model vs. FTIR 0–7 km partial column of CO for Eureka, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and $R^2$ are noted in the legend. The 1 : 1 line is shown in light grey.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9.** By model and location: (a) overall model-measurement mean percent difference for CO 0–7 km partial columns ($\Delta_O$), with error bars that represent the standard deviation of the mean, as shown in the legend of Figs. B5–B8. (b) $R^2$ as shown in Figs. 8 and B9–B11. (c) Normalized root-mean-square error.</figcaption></figure>

sites – Zeppelin, Norway, and Utqiaġvik (formerly Barrow), USA – with a –20% to 30% bias in the first 6 months of the year (Whaley et al., 2023), which is compatible with results shown in Fig. 10e.

In POLMIP, models were run for 2008 with a standardized emissions inventory; there is some overlap of models examined here, although a different emissions input was used

(see Emmons et al., 2015, for full project description). Similarly to the results presented here, the POLMIP study found that relative to surface, airborne, and satellite Arctic tropospheric measurements, CO was underpredicted by the models (MMM gross error 9%–12%), with a more negative bias in the winter and spring compared to the summer, although the models still broadly captured the seasonal cycle (Monks
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Smoothed model vs. FTIR 0–7 km partial column of CO for Eureka, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and $R^2$ are noted in the legend. The 1 : 1 line is shown in light grey.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 9.** By model and location: (a) overall model-measurement mean percent difference for CO 0–7 km partial columns ($\Delta_O$), with error bars that represent the standard deviation of the mean, as shown in the legend of Figs. B5–B8. (b) $R^2$ as shown in Figs. 8 and B9–B11. (c) Normalized root-mean-square error.</figcaption></figure>

sites – Zeppelin, Norway, and Utqiaġvik (formerly Barrow), USA – with a –20% to 30% bias in the first 6 months of the year (Whaley et al., 2023), which is compatible with results shown in Fig. 10e.

In POLMIP, models were run for 2008 with a standardized emissions inventory; there is some overlap of models examined here, although a different emissions input was used

(see Emmons et al., 2015, for full project description). Similarly to the results presented here, the POLMIP study found that relative to surface, airborne, and satellite Arctic tropospheric measurements, CO was underpredicted by the models (MMM gross error 9%–12%), with a more negative bias in the winter and spring compared to the summer, although the models still broadly captured the seasonal cycle (Monks
<figure><img src="image_12.png" /><figcaption>**Figure 10.** (a–d) Monthly mean FTIR (black) and multi-model mean (coloured) 0–7 km partial columns of CO (PC<sub>F,monthly,j</sub> and PC<sub>MMM,monthly,j</sub>, respectively), with error bars and shaded areas representing the standard deviation of the mean. (e) Monthly mean percent difference in the MMM (Δ<sub>O,MMM</sub>) for all locations.</figcaption></figure>

et al., 2015). Using an idealized tracer, POLMIP examined anthropogenic and biomass burning influences in Arctic regions, demonstrating a seasonal dependence of transport efficiency. It was shown that for anthropogenic emissions, Europe influences the surface CO, while Asia and North America have more influence higher in the troposphere (Monks et al., 2015). Furthermore, the tracer investigation in that study showed that OH differences account for more variability between the models than the transport mechanisms within the individual models. However, it can be noted that although models may reduce negative biases through better OH chemistry, this alone will not resolve the differences between the model and measurements (Monks et al., 2015).

The current study, the POLMIP study, and the AMAP report exhibit similarities in the model-measurement comparisons of CO. Most notably, all three studies show negative biases early in the year, which shift positively in the summer; the model-FTIR comparisons become less negative, while the AMAP-surface measurement comparisons change to a positive bias. Lutsch et al. (2020) also reported a low bias in GEOS-Chem lower-tropospheric CO columns compared with measurements from 10 FTIR stations, including four sites in this study, although they found a greater underestimation for Eureka and Thule in July and August due to transported boreal wildfire emissions not being fully captured by the model, particularly for years after 2015 not included in the present study. Previously published studies point to underestimated anthropogenic emissions as a source of the discrepancies (Monks et al., 2015; Whaley et al., 2022, 2023). The results of the model-FTIR comparisons presented here support this reasoning, as the only model with a positive bias (GEM-MACH) has additional local Arctic emissions (Gong et al., 2018). The models may be improved with more refined

OH chemistry, although it is unlikely to completely resolve the inconsistencies (Monks et al., 2015); improvements to long-range transport and biomass burning inventories could also reduce the differences between model results and measurements.

### 4.3 O<sub>3</sub>

Tropospheric O<sub>3</sub> is both a significant anthropogenic GHG and an air pollutant that has impacts on human health and ecosystems. In the troposphere, O<sub>3</sub> is a secondary pollutant, produced by photochemical oxidation of volatile organic compounds in the presence of NO<sub>x</sub>. In addition to atmospheric photochemistry, its production is highly sensitive to meteorological conditions. Diurnal impacts on O<sub>3</sub> production are minimal in the Arctic, relative to lower latitudes, due to the gradual and prolonged change in solar altitude and angle throughout the year. While O<sub>3</sub> processes are complex, O<sub>3</sub> is often quite well reproduced by models, possibly due to compensating biases in its precursors (Whaley et al., 2022). Although progress has been made, sparse observations, Arctic amplification, and a changing global climate hinder the understanding and modelling of O<sub>3</sub> in Arctic regions (Whaley et al., 2023). For a summary of the current understanding of Arctic tropospheric O<sub>3</sub>, see Whaley et al. (2023).

All 11 of the models examined in this study provide 3-hourly O<sub>3</sub> concentrations. The full data time series plots (Figs. C1–C5) demonstrate the variation between the models and throughout the year, which is likely a by-product of the complexity of modelling tropospheric O<sub>3</sub>. Figure 11 and Figs. C6–C10 show the monthly mean partial columns (panels a–e) and percent differences (panel f) to highlight the parts of the year which are over- or underpredicted. For
<figure><img src="image_13.png" /><figcaption>**Figure 10.** (a–d) Monthly mean FTIR (black) and multi-model mean (coloured) 0–7 km partial columns of CO (PC<sub>F,monthly,j</sub> and PC<sub>MMM,monthly,j</sub>, respectively), with error bars and shaded areas representing the standard deviation of the mean. (e) Monthly mean percent difference in the MMM (Δ<sub>O,MMM</sub>) for all locations.</figcaption></figure>

et al., 2015). Using an idealized tracer, POLMIP examined anthropogenic and biomass burning influences in Arctic regions, demonstrating a seasonal dependence of transport efficiency. It was shown that for anthropogenic emissions, Europe influences the surface CO, while Asia and North America have more influence higher in the troposphere (Monks et al., 2015). Furthermore, the tracer investigation in that study showed that OH differences account for more variability between the models than the transport mechanisms within the individual models. However, it can be noted that although models may reduce negative biases through better OH chemistry, this alone will not resolve the differences between the model and measurements (Monks et al., 2015).

The current study, the POLMIP study, and the AMAP report exhibit similarities in the model-measurement comparisons of CO. Most notably, all three studies show negative biases early in the year, which shift positively in the summer; the model-FTIR comparisons become less negative, while the AMAP-surface measurement comparisons change to a positive bias. Lutsch et al. (2020) also reported a low bias in GEOS-Chem lower-tropospheric CO columns compared with measurements from 10 FTIR stations, including four sites in this study, although they found a greater underestimation for Eureka and Thule in July and August due to transported boreal wildfire emissions not being fully captured by the model, particularly for years after 2015 not included in the present study. Previously published studies point to underestimated anthropogenic emissions as a source of the discrepancies (Monks et al., 2015; Whaley et al., 2022, 2023). The results of the model-FTIR comparisons presented here support this reasoning, as the only model with a positive bias (GEM-MACH) has additional local Arctic emissions (Gong et al., 2018). The models may be improved with more refined

OH chemistry, although it is unlikely to completely resolve the inconsistencies (Monks et al., 2015); improvements to long-range transport and biomass burning inventories could also reduce the differences between model results and measurements.

### 4.3 O<sub>3</sub>

Tropospheric O<sub>3</sub> is both a significant anthropogenic GHG and an air pollutant that has impacts on human health and ecosystems. In the troposphere, O<sub>3</sub> is a secondary pollutant, produced by photochemical oxidation of volatile organic compounds in the presence of NO<sub>x</sub>. In addition to atmospheric photochemistry, its production is highly sensitive to meteorological conditions. Diurnal impacts on O<sub>3</sub> production are minimal in the Arctic, relative to lower latitudes, due to the gradual and prolonged change in solar altitude and angle throughout the year. While O<sub>3</sub> processes are complex, O<sub>3</sub> is often quite well reproduced by models, possibly due to compensating biases in its precursors (Whaley et al., 2022). Although progress has been made, sparse observations, Arctic amplification, and a changing global climate hinder the understanding and modelling of O<sub>3</sub> in Arctic regions (Whaley et al., 2023). For a summary of the current understanding of Arctic tropospheric O<sub>3</sub>, see Whaley et al. (2023).

All 11 of the models examined in this study provide 3-hourly O<sub>3</sub> concentrations. The full data time series plots (Figs. C1–C5) demonstrate the variation between the models and throughout the year, which is likely a by-product of the complexity of modelling tropospheric O<sub>3</sub>. Figure 11 and Figs. C6–C10 show the monthly mean partial columns (panels a–e) and percent differences (panel f) to highlight the parts of the year which are over- or underpredicted. For
example, “springtime” (referred to here as when the sun rises, in approximately late February at the highest-latitude sites, until May) $O_3$ is of interest in the Arctic due to the springtime maximum in its seasonal cycle and the potential for both stratospheric ozone intrusions into the upper (mid-troposphere and surface $O_3$ depletion events (ODEs) due to bromine explosions and halogen chemistry. However, the FTIR 0–7 km partial-column $O_3$ seasonal cycle, shown here, is dominated by the free troposphere and stratospheric processes and does not have a springtime minimum from surface ODEs, as one might expect from surface measurements (Solberg et al., 1996; Berg et al., 2003; Skov et al., 2006; Eneroth et al., 2007; Whaley et al., 2023). The Arctic surface ODE features are primarily limited to the near-surface/lower boundary layer (< 2 km), whereas the 0–7 km partial column is dominated by the free troposphere (Zhao et al., 2016). It can be noted that all of the models in this study lack the necessary halogen chemistry needed to simulate ODEs in the high Arctic (Whaley et al., 2023). Figure 11 shows that across all locations, MATCH-SALSA overpredicts $O_3$ by 35%–75% in winter, which gradually declines until May, after which the bias becomes negative. GEM-MACH, GEOS-Chem, UKESM1, and WRF-Chem underestimate springtime $O_3$ most substantially across all sites. The discrepancies may arise from inaccuracies in model water vapor leading to an increase in $O_3$ destruction and/or a lack of $O_3$ transported from mid-latitudes, which is a substantial source of tropospheric $O_3$ in the Arctic (Hirdman et al., 2010; Whaley et al., 2023). In the case of the regional GEM-MACH model, low biases in $O_3$ or precursor species at the lateral boundary conditions may also be contributing. CESM, CMAM, DEHM, and MRI-ESM2 demonstrate reasonable agreement with measured springtime $O_3$ across locations, in addition to a smaller overall mean percent difference, relative to other models. EMEP MSC-W and WRF-Chem simulate springtime $O_3$ comparable to the aforementioned models, although negative biases later in the year lead to a larger overall mean percent difference. This may indicate that these models have too much photochemical $O_3$ loss in the summer months.

Figure 12 and Figs. C11–C14 show the model versus FTIR $O_3$ 0–7 km partial columns, with the line of best fit and $R^2$ shown in the legend, along with the 1:1 line. The general underprediction towards the largest values could be related to the underestimation in precursor species (such as CO or $NO_x$), a lack of long-range transport, an underestimation of ozone production in air masses during long-range transport to the Arctic, or a combination thereof. Using a MOZART-4 tagged tracer simulation of $O_3$, Wespes et al. (2012) examined source attributions of the tropospheric $O_3$ columns measured by the FTIR instruments at Thule and Eureka. Their analysis shows that the retrievals have a minimal contribution from the a priori profile (~1%), resulting in high vertical sensitivity throughout the troposphere. The tropospheric column source contributions were estimated, where over half was attributed to anthropogenic sources, followed by strato-

spheric influence and lastly lightning and biomass burning emissions (Wespes et al., 2012). The seasonal cycle of Arctic $O_3$ has been shown to vary based on geographical conditions, such as whether the site is coastal, inland, or at a high elevation (Whaley et al., 2023). Moreover, $O_3$ partial columns can be variable because they depend on the vertical distribution of $O_3$, which is determined by a combination of emissions, chemistry, dynamics, and radiation, all of which vary with altitude (Rap et al., 2015). Notably, Arctic $O_3$ columns have strong gradients in the influences on the vertical profile from mid-latitude regions (Europe, North America and Asia), which also vary with season (Monks et al., 2015). The combination of these factors leads to an increasingly complex series of model processes, which can also result in compounding errors. Without sensitivity simulations, like those carried out in Monks et al. (2015) and Rap et al. (2015), it is difficult to definitively say which of these processes are responsible for the underestimations found in this study.

Figure 13 shows the summary of $O_3$ mean percent differences, $R^2$, and normalized root-mean-square error. The model–FTIR comparisons reveal that the spatial resolution, and inclusion of stratospheric chemistry in the models do not necessarily improve results (refer to Table 3 for horizontal resolution and stratospheric chemistry). For example, WRF-Chem, EMEP MSC-W, and GEM-MACH show a low $R^2$ and higher NRMSE (varying between sites and models), although contributing to this for WRF-Chem and GEM-MACH could be the limited number of analysis years (two and one, respectively). These air-quality focused models have detailed chemistry and were run at a higher spatial resolutions, whereas for example CMAM, a climate-focused model, has a coarser resolution with simplified tropospheric chemistry and demonstrates larger $R^2$ and smaller mean percent differences (Fig. 13). However, when considering the stratosphere, CMAM, which includes comprehensive stratospheric chemistry, has comparable metrics in Fig. 13 to DEHM, which uses prescribed climatologies for the stratosphere. Similarly, Whaley et al. (2022) stated that the degree of stratospheric chemistry in the models did not reveal a consistent benefit or handicap when comparing the models with surface measurements. Here, the $O_3$ partial column comparisons show significant variation, although again models largely underpredict FTIR measurements. The $R^2$, mean percent difference, and NRMSE are relatively consistent, where models with a larger percent difference also have weaker correlations and higher NRMSEs. An exception to this is CESM, which has one of the smallest overall differences across the models and locations. However, in the model vs. FTIR plot(s) (and Figs. 12, C11–C14), CESM has considerable scatter above and below the line of best fit, resulting in a decreased mean difference, while also reducing $R^2$, unlike MRI-ESM2, which has a similar mean percent difference and NRMSE but a stronger linear correlation.

To supplement the aircraft and satellite campaigns undertaken for the POLARCAT study, daily mean $O_3$ measure-
example, “springtime” (referred to here as when the sun rises, in approximately late February at the highest-latitude sites, until May) O₃ is of interest in the Arctic due to the springtime maximum in its seasonal cycle and the potential for both stratospheric ozone intrusions into the upper (mid- ) troposphere and surface O₃ depletion events (ODEs) due to bromine explosions and halogen chemistry. However, the FTIR 0–7 km partial-column O₃ seasonal cycle, shown here, is dominated by the free troposphere and stratospheric processes and does not have a springtime minimum from surface ODEs, as one might expect from surface measurements (Solberg et al., 1996; Berg et al., 2003; Skov et al., 2006; Eneroth et al., 2007; Whaley et al., 2023). The Arctic surface ODE features are primarily limited to the near-surface/lower boundary layer (< 2 km), whereas the 0–7 km partial column is dominated by the free troposphere (Zhao et al., 2016). It can be noted that all of the models in this study lack the necessary halogen chemistry needed to simulate ODEs in the high Arctic (Whaley et al., 2023). Figure 11 shows that across all locations, MATCH-SALSA overpredicts O₃ by 35%–75% in winter, which gradually declines until May, after which the bias becomes negative. GEM-MACH, GEOS-Chem, UKESM1, and WRF-Chem underestimate springtime O₃ most substantially across all sites. The discrepancies may arise from inaccuracies in model water vapor leading to an increase in O₃ destruction and/or a lack of O₃ transported from mid-latitudes, which is a substantial source of tropospheric O₃ in the Arctic (Hirdman et al., 2010; Whaley et al., 2023). In the case of the regional GEM-MACH model, low biases in O₃ or precursor species at the lateral boundary conditions may also be contributing. CESM, CMAM, DEHM, and MRI-ESM2 demonstrate reasonable agreement with measured springtime O₃ across locations, in addition to a smaller overall mean percent difference, relative to other models. EMEP MSC-W and WRF-Chem simulate springtime O₃ comparable to the aforementioned models, although negative biases later in the year lead to a larger overall mean percent difference. This may indicate that these models have too much photochemical O₃ loss in the summer months.

Figure 12 and Figs. C11–C14 show the model versus FTIR O₃ 0–7 km partial columns, with the line of best fit and R² shown in the legend, along with the 1 : 1 line. The general underprediction towards the largest values could be related to the underestimation in precursor species (such as CO or NOₓ), a lack of long-range transport, an underestimation of ozone production in air masses during long-range transport to the Arctic, or a combination thereof. Using a MOZART-4 tagged tracer simulation of O₃, Wespes et al. (2012) examined source attributions of the tropospheric O₃ columns measured by the FTIR instruments at Thule and Eureka. Their analysis shows that the retrievals have a minimal contribution from the a priori profile (~ 1 %), resulting in high vertical sensitivity throughout the troposphere. The tropospheric column source contributions were estimated, where over half was attributed to anthropogenic sources, followed by strato-

spheric influence and lastly lightning and biomass burning emissions (Wespes et al., 2012). The seasonal cycle of Arctic O₃ has been shown to vary based on geographical conditions, such as whether the site is coastal, inland, or at a high elevation (Whaley et al., 2023). Moreover, O₃ partial columns can be variable because they depend on the vertical distribution of O₃, which is determined by a combination of emissions, chemistry, dynamics, and radiation, all of which vary with altitude (Rap et al., 2015). Notably, Arctic O₃ columns have strong gradients in the influences on the vertical profile from mid-latitude regions (Europe, North America and Asia), which also vary with season (Monks et al., 2015). The combination of these factors leads to an increasingly complex series of model processes, which can also result in compounding errors. Without sensitivity simulations, like those carried out in Monks et al. (2015) and Rap et al. (2015), it is difficult to definitively say which of these processes are responsible for the underestimations found in this study.

Figure 13 shows the summary of O₃ mean percent differences, R², and normalized root-mean-square error. The model–FTIR comparisons reveal that the spatial resolution, and inclusion of stratospheric chemistry in the models do not necessarily improve results (refer to Table 3 for horizontal resolution and stratospheric chemistry). For example, WRF-Chem, EMEP MSC-W, and GEM-MACH show a low R² and higher NRMSE (varying between sites and models), although contributing to this for WRF-Chem and GEM-MACH could be the limited number of analysis years (two and one, respectively). These air-quality focused models have detailed chemistry and were run at a higher spatial resolutions, whereas for example CMAM, a climate-focused model, has a coarser resolution with simplified tropospheric chemistry and demonstrates larger R² and smaller mean percent differences (Fig. 13). However, when considering the stratosphere, CMAM, which includes comprehensive stratospheric chemistry, has comparable metrics in Fig. 13 to DEHM, which uses prescribed climatologies for the stratosphere. Similarly, Whaley et al. (2022) stated that the degree of stratospheric chemistry in the models did not reveal a consistent benefit or handicap when comparing the models with surface measurements. Here, the O₃ partial column comparisons show significant variation, although again models largely underpredict FTIR measurements. The R², mean percent difference, and NRMSE are relatively consistent, where models with a larger percent difference also have weaker correlations and higher NRMSEs. An exception to this is CESM, which has one of the smallest overall differences across the models and locations. However, in the model vs. FTIR plot(s) (and Figs. 12, C11–C14), CESM has considerable scatter above and below the line of best fit, resulting in a decreased mean difference, while also reducing R², unlike MRI-ESM2, which has a similar mean percent difference and NRMSE but a stronger linear correlation.

To supplement the aircraft and satellite campaigns undertaken for the POLARCAT study, daily mean O₃ measure-
<figure><img src="image_13.png" /><figcaption>**Figure 11.** (a–e) Monthly mean FTIR (black) and smoothed model (colour) 0–7 km partial columns of O₃ (PC<sub>F,monthly,j</sub> and PC<sub>M,monthly,j</sub>, respectively), for each location, shown with the same y axis. Error bars represent the standard deviation of the monthly mean. (f) Model-measurement mean percent difference by month (Δ<sub>monthly,j</sub>) for each model (by colour) and location (by marker). Error bars represent standard deviation of the monthly mean percent difference.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 12.** Smoothed model vs. FTIR 0–7 km partial columns of O₃ for Eureka, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and R² are noted in the legend. The 1 : 1 line is shown in light grey.</figcaption></figure>

ments from the FTIR instruments at Eureka and Thule were
compared to MOZART-4 simulations in Wespes et al. (2012).
When examining a partial column from the ground to
300 hPa (approximately 9 km), the smoothed model showed
a bias of –15 % relative to the FTIR. This is consistent with
their analysis of aircraft observations, which revealed that
the model underestimated O₃ by 5 %–15 %. Results here are
similar to those presented in Wespes et al. (2012), where

across all the locations and models, 24 of the 55 model-
measurement mean percent differences were within ±15%
(see Fig. 15). The FTIR uncertainty for O₃ partial columns
ranges from 3.9% to 8.2%; the overall mean percent dif-
ference for MATCH-SALSA falls within these uncertainty
bounds for all locations, and CESM, DEHM, MATCH, and
MRI-ESM2 are within FTIR uncertainty for all locations but
Ny Ålesund.
<figure><img src="image_14.png" /><figcaption>**Figure 11.** (a–e) Monthly mean FTIR (black) and smoothed model (colour) 0–7 km partial columns of O₃ (PC<sub>F,monthly,j</sub> and PC<sub>M,monthly,j</sub>, respectively), for each location, shown with the same y axis. Error bars represent the standard deviation of the monthly mean. (f) Model-measurement mean percent difference by month (Δ<sub>monthly,j</sub>) for each model (by colour) and location (by marker). Error bars represent standard deviation of the monthly mean percent difference.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 12.** Smoothed model vs. FTIR 0–7 km partial columns of O₃ for Eureka, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and R² are noted in the legend. The 1 : 1 line is shown in light grey.</figcaption></figure>

ments from the FTIR instruments at Eureka and Thule were
compared to MOZART-4 simulations in Wespes et al. (2012).
When examining a partial column from the ground to
300 hPa (approximately 9 km), the smoothed model showed
a bias of –15 % relative to the FTIR. This is consistent with
their analysis of aircraft observations, which revealed that
the model underestimated O₃ by 5 %–15 %. Results here are
similar to those presented in Wespes et al. (2012), where

across all the locations and models, 24 of the 55 model-
measurement mean percent differences were within ±15%
(see Fig. 15). The FTIR uncertainty for O3 partial columns
ranges from 3.9% to 8.2%; the overall mean percent dif-
ference for MATCH-SALSA falls within these uncertainty
bounds for all locations, and CESM, DEHM, MATCH, and
MRI-ESM2 are within FTIR uncertainty for all locations but
Ny Ålesund.
<figure><img src="image_15.png" /><figcaption>**Figure 13.** By model and location: (a) overall model-measurement mean percent difference for O₃ 0–7 km partial columns (ΔO), with error bars that represent the standard deviation of the mean, as shown in the legend of Figs. C6–C10. (b) R² as shown in Figs. 12 and C11–C14. (c) Normalized root-mean-square error.</figcaption></figure>

The AMAP SLCF assessment report finds that the multi-model mean of Arctic O₃ has a bias of +11 ± 3% relative to surface measurements (AMAP, 2021). When partitioning results by region, all the models had positive biases when compared to the surface measurements in Alaska and negative biases in northern Europe, resulting in a relatively small mean bias across the Arctic as a whole (Whaley et al., 2022). Inaccuracies in long-range transport of O₃ and its precursors may have contributed to the increased discrepancy seen in the model-FTIR comparisons of the current study, particularly in partial columns with larger values. For example, the underestimation of CO may contribute to the negative bias in O₃ (see Figs. 9–10). Most models in AMAP (2021) show negative biases for Greenland and northern European locations, which would correspond closer geographically with the FTIR sites examined here. When comparing the AMAP models to TES (Tropospheric Emission Spectrometer) and ACE-FTS satellite O₃ measurements, the biases are negative at lower altitudes and become positive at higher altitudes (Whaley et al., 2022). AMAP model vs. ozonesonde comparisons showed similar elevated positive biases around 6–8 km of up to ±50%, again indicating that the models may produce too much O₃ from mid-latitude anthropogenic emissions or that there may be too much downward transport of O₃ from the stratosphere (Whaley et al., 2023). The best performance in that study came from the multi-model mean, which simulated O₃ within ±8% throughout the troposphere.

Figure 14 shows the monthly MMM for O₃ at all locations, along with the monthly mean FTIR and the associated percent difference. This shows that the models, as a

whole, have an increased negative bias in the middle of the year relative to the winter, while still exhibiting a negative bias overall. The longitudinal range of sites examined here may limit biases to be negative, not capturing the positive-negative gradient from west-east in O₃ found in the AMAP report (AMAP, 2021; Whaley et al., 2022). Nonetheless, the model-FTIR O₃ comparisons reflect the proclivity of the models to underpredict Arctic O₃ in the lower troposphere, as also found in the aforementioned studies. The results of this study agree with results from previous studies and suggest that improvements are still needed for accurate modelling of O₃ and CO in the Arctic (Whaley et al., 2023). Models still require improvements in their treatment of stratosphere-troposphere exchange and Arctic boundary layer processes to better simulate Arctic O₃. Also required are better understanding and implementation of processes influencing O₃ removal through dry deposition and O₃ photochemical production from anthropogenic, biomass burning, and natural sources in the lower and middle troposphere.

## 5 Conclusions

This study compares atmospheric models with data from five Arctic NDACC ground-based FTIR spectrometers. The models simulate SLCFs and precursor gases with 3-hourly outputs for the years 2008, 2009, 2014, and 2015. Here, a total of 3 models are evaluated for CH₄, 9 for CO, and 11 for O₃. The model simulations are compared with FTIR tropospheric partial column measurements to assess performance throughout the year and across locations.
<figure><img src="image_16.png" /><figcaption>**Figure 13.** By model and location: (a) overall model-measurement mean percent difference for O₃ 0–7 km partial columns (ΔO), with error bars that represent the standard deviation of the mean, as shown in the legend of Figs. C6–C10. (b) R² as shown in Figs. 12 and C11–C14. (c) Normalized root-mean-square error.</figcaption></figure>

The AMAP SLCF assessment report finds that the multi-model mean of Arctic O₃ has a bias of +11 ± 3% relative to surface measurements (AMAP, 2021). When partitioning results by region, all the models had positive biases when compared to the surface measurements in Alaska and negative biases in northern Europe, resulting in a relatively small mean bias across the Arctic as a whole (Whaley et al., 2022). Inaccuracies in long-range transport of O₃ and its precursors may have contributed to the increased discrepancy seen in the model-FTIR comparisons of the current study, particularly in partial columns with larger values. For example, the underestimation of CO may contribute to the negative bias in O₃ (see Figs. 9–10). Most models in AMAP (2021) show negative biases for Greenland and northern European locations, which would correspond closer geographically with the FTIR sites examined here. When comparing the AMAP models to TES (Tropospheric Emission Spectrometer) and ACE-FTS satellite O₃ measurements, the biases are negative at lower altitudes and become positive at higher altitudes (Whaley et al., 2022). AMAP model vs. ozonesonde comparisons showed similar elevated positive biases around 6–8 km of up to ±50%, again indicating that the models may produce too much O₃ from mid-latitude anthropogenic emissions or that there may be too much downward transport of O₃ from the stratosphere (Whaley et al., 2023). The best performance in that study came from the multi-model mean, which simulated O₃ within ±8% throughout the troposphere.

Figure 14 shows the monthly MMM for O₃ at all locations, along with the monthly mean FTIR and the associated percent difference. This shows that the models, as a

whole, have an increased negative bias in the middle of the year relative to the winter, while still exhibiting a negative bias overall. The longitudinal range of sites examined here may limit biases to be negative, not capturing the positive-negative gradient from west-east in O₃ found in the AMAP report (AMAP, 2021; Whaley et al., 2022). Nonetheless, the model-FTIR O₃ comparisons reflect the proclivity of the models to underpredict Arctic O₃ in the lower troposphere, as also found in the aforementioned studies. The results of this study agree with results from previous studies and suggest that improvements are still needed for accurate modelling of O₃ and CO in the Arctic (Whaley et al., 2023). Models still require improvements in their treatment of stratosphere-troposphere exchange and Arctic boundary layer processes to better simulate Arctic O₃. Also required are better understanding and implementation of processes influencing O₃ removal through dry deposition and O₃ photochemical production from anthropogenic, biomass burning, and natural sources in the lower and middle troposphere.

## 5 Conclusions

This study compares atmospheric models with data from five Arctic NDACC ground-based FTIR spectrometers. The models simulate SLCFs and precursor gases with 3-hourly outputs for the years 2008, 2009, 2014, and 2015. Here, a total of 3 models are evaluated for CH₄, 9 for CO, and 11 for O₃. The model simulations are compared with FTIR tropospheric partial column measurements to assess performance throughout the year and across locations.
<figure><img src="image_16.png" /><figcaption>**Figure 14.** (a–e) Monthly mean FTIR (black) and multi-model mean (coloured) 0–7 km partial columns of O₃ (PC<sub>F,monthly,j</sub> and PC<sub>MMM,monthly,j</sub>, respectively), with error bars and shaded areas representing the standard deviation of the mean. (f) Monthly mean percent difference in the MMM (Δ<sub>O,MMM</sub>) for all locations.</figcaption></figure>

Generally, across the five locations, the model simulations of 0–7 km partial columns of CH₄, CO, and O₃ are underestimated. There were no significant patterns in the biases identified between the sites, species, or models examined. Modelled CH₄ partial columns are relatively consistent across the year, broadly capturing seasonal cycles, with the exception of a few outliers. CO simulations are inconsistent in reproducing the seasonal cycle, underpredicting springtime partial columns compared to the rest of the year, and skewing differences to be more positive when there are enhancements due to biomass burning events. Similarly, the models underestimated O₃ maxima more than O₃ minima in the troposphere. The multi-model means are reflective of these trends, for which (ignoring outliers) the CH₄ mean percent difference is relatively consistent across the year, CO has a maximum difference in the spring and a minimum in the summer, and O₃ has maximum difference centered around the summer. The AMAP SLCF assessment report found the best results using a multi-model mean for all species when comparing it with surface measurements (AMAP, 2021; Whaley et al., 2022). However, here, the multi-model means of the tropospheric column for all species are biased low. The average MMM mean difference is approximately –10% for CH₄, –21% for CO, and –18% for O₃ (see Table 4), where the uncertainty in the FTIR 0–7 km partial column is on the order of 6% on average. When examining the models and location pairs individually, the mean difference (inclusive of standard deviation) is within the respective FTIR uncertainty, for 6 of 15 model–FTIR comparisons for CH₄, 12 of 34 comparisons for CO, and 25 of 55 comparisons for O₃ (see Fig. 15).

These evaluations show that models are lacking some degree of transport and/or emissions to accurately reproduce tropospheric columns and seasonal variability in the Arc-

<figure><img src="image_17.png" /><figcaption>**Figure 15.** Summary of model-measurement mean percent difference (Δ<sub>O</sub>) for each model and location by species. MMM is the multi-model mean (Δ<sub>O,MMM</sub>). The colour scale indicates the mean percent difference relative to the FTIR measurements, from blue (–50%) to red (+50%). A square marker indicates that the mean percent difference is within the FTIR uncertainty. A triangle marker indicates that the mean difference is within the FTIR uncertainty combined with the standard deviation of the monthly mean percent difference.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 14.** (a–e) Monthly mean FTIR (black) and multi-model mean (coloured) 0–7 km partial columns of O₃ (PC<sub>F,monthly,j</sub> and PC<sub>MMM,monthly,j</sub>, respectively), with error bars and shaded areas representing the standard deviation of the mean. (f) Monthly mean percent difference in the MMM (Δ<sub>O,MMM</sub>) for all locations.</figcaption></figure>

Generally, across the five locations, the model simulations of 0–7 km partial columns of CH₄, CO, and O₃ are underestimated. There were no significant patterns in the biases identified between the sites, species, or models examined. Modelled CH₄ partial columns are relatively consistent across the year, broadly capturing seasonal cycles, with the exception of a few outliers. CO simulations are inconsistent in reproducing the seasonal cycle, underpredicting springtime partial columns compared to the rest of the year, and skewing differences to be more positive when there are enhancements due to biomass burning events. Similarly, the models underestimated O₃ maxima more than O₃ minima in the troposphere. The multi-model means are reflective of these trends, for which (ignoring outliers) the CH₄ mean percent difference is relatively consistent across the year, CO has a maximum difference in the spring and a minimum in the summer, and O₃ has maximum difference centered around the summer. The AMAP SLCF assessment report found the best results using a multi-model mean for all species when comparing it with surface measurements (AMAP, 2021; Whaley et al., 2022). However, here, the multi-model means of the tropospheric column for all species are biased low. The average MMM mean difference is approximately –10% for CH₄, –21% for CO, and –18% for O₃ (see Table 4), where the uncertainty in the FTIR 0–7 km partial column is on the order of 6% on average. When examining the models and location pairs individually, the mean difference (inclusive of standard deviation) is within the respective FTIR uncertainty, for 6 of 15 model–FTIR comparisons for CH₄, 12 of 34 comparisons for CO, and 25 of 55 comparisons for O₃ (see Fig. 15).

These evaluations show that models are lacking some degree of transport and/or emissions to accurately reproduce tropospheric columns and seasonal variability in the Arc-

<figure><img src="image_18.png" /><figcaption>**Figure 15.** Summary of model-measurement mean percent difference (Δ<sub>O</sub>) for each model and location by species. MMM is the multi-model mean (Δ<sub>O,MMM</sub>). The colour scale indicates the mean percent difference relative to the FTIR measurements, from blue (–50%) to red (+50%). A square marker indicates that the mean percent difference is within the FTIR uncertainty. A triangle marker indicates that the mean difference is within the FTIR uncertainty combined with the standard deviation of the monthly mean percent difference.</figcaption></figure>
tic. Model evaluation can provide a valuable checkpoint to help improve the representation of the Arctic in atmospheric models. NDACC FTIR spectrometers were selected for this project because of the wide range of species measured, high spectral resolution, multiple high-latitude sites, and publicly available data; in addition, the column-integrated FTIR measurements used in this study have a spatial and temporal footprint that is more representative of the free troposphere than in situ and satellite measurements. Future work would benefit from the inclusion of sensitivity studies, furthering the model-measurement comparisons with mid-latitude NDACC FTIR sites and extending comparisons to a longer time frame, with some models and locations having data from as early as 1990.

## Appendix A: Additional figures for CH₄

<figure><img src="image_18.png" /><figcaption>**Figure A1.** (a) FTIR (black) and smoothed model (colour) 0–7 km partial columns of CH₄ by day of year, from Eureka. Model data are the nearest in time to each FTIR measurement. (b) Model-measurement percent difference (Δᵢ) from Eq. (2) by day of year. Each year is indicated by a different marker.</figcaption></figure>
tic. Model evaluation can provide a valuable checkpoint to help improve the representation of the Arctic in atmospheric models. NDACC FTIR spectrometers were selected for this project because of the wide range of species measured, high spectral resolution, multiple high-latitude sites, and publicly available data; in addition, the column-integrated FTIR measurements used in this study have a spatial and temporal footprint that is more representative of the free troposphere than in situ and satellite measurements. Future work would benefit from the inclusion of sensitivity studies, furthering the model-measurement comparisons with mid-latitude NDACC FTIR sites and extending comparisons to a longer time frame, with some models and locations having data from as early as 1990.

## Appendix A: Additional figures for CH$_4$

<figure><img src="image_19.png" /><figcaption>**Figure A1.** (a) FTIR (black) and smoothed model (colour) 0–7 km partial columns of CH$_4$ by day of year, from Eureka. Model data are the nearest in time to each FTIR measurement. (b) Model-measurement percent difference ($\Delta_i$) from Eq. (2) by day of year. Each year is indicated by a different marker.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>Figure A2. Same as Fig. A1 but for Ny Ålesund.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>Figure A3. Same as Fig. A1 but for Thule.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>Figure A2. Same as Fig. A1 but for Ny Ålesund.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>Figure A3. Same as Fig. A1 but for Thule.</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>Figure A4. Same as Fig. A1 but for Kiruna.</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>Figure A5. Same as Fig. A1 but for Harestua.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>Figure A4. Same as Fig. A1 but for Kiruna.</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>Figure A5. Same as Fig. A1 but for Harestua.</figcaption></figure>
<figure><img src="image_23.png" /></figure>

**Figure A6.** (a) Monthly mean FTIR (black) and smoothed model (colour) 0–7 km partial columns of CH₄ (PCF,monthly,j and PCM,monthly,j, respectively), from Eureka using model data that are the nearest in time to each FTIR measurement shown in Fig. A1. Error bars represent the standard deviation of the monthly mean. (b) Model-measurement mean percent difference by month (Δ<sub>monthly,j</sub>). Error bars represent standard deviation of the monthly mean percent difference. The legend in panel (b) shows the overall mean percent difference (Δ<sub>O</sub>) with the standard deviation of the overall mean percent difference.

<figure><img src="image_24.png" /></figure>

**Figure A7.** Same as Fig. A6 but for Ny Ålesund.

<figure><img src="image_25.png" /></figure>

**Figure A8.** Same as Fig. A6 but for Thule.
<figure><img src="image_24.png" /></figure>

**Figure A6.** (a) Monthly mean FTIR (black) and smoothed model (colour) 0–7 km partial columns of CH₄ (PCF,monthly,j and PCM,monthly,j, respectively), from Eureka using model data that are the nearest in time to each FTIR measurement shown in Fig. A1. Error bars represent the standard deviation of the monthly mean. (b) Model-measurement mean percent difference by month (Δ<sub>monthly,j</sub>). Error bars represent standard deviation of the monthly mean percent difference. The legend in panel (b) shows the overall mean percent difference (Δ<sub>O</sub>) with the standard deviation of the overall mean percent difference.

<figure><img src="image_25.png" /></figure>

**Figure A7.** Same as Fig. A6 but for Ny Ålesund.

<figure><img src="image_26.png" /></figure>

**Figure A8.** Same as Fig. A6 but for Thule.
<figure><img src="image_26.png" /><figcaption>**Figure A9.** Same as Fig. A6 but for Kiruna.</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>**Figure A10.** Same as Fig. A6 but for Harestua.</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>**Figure A11.** Smoothed model vs. FTIR 0–7 km partial column of CH₄ for Ny Ålesund, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and R² are noted in the legend. The 1 : 1 line is shown in light grey.</figcaption></figure>
<figure><img src="image_27.png" /><figcaption>**Figure A9.** Same as Fig. A6 but for Kiruna.</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>**Figure A10.** Same as Fig. A6 but for Harestua.</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>**Figure A11.** Smoothed model vs. FTIR 0–7 km partial column of CH₄ for Ny Ålesund, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and R² are noted in the legend. The 1 : 1 line is shown in light grey.</figcaption></figure>
<figure><img src="image_29.png" /><figcaption>**Figure A12.** Same as Fig. A9 but for Thule.</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>**Figure A13.** Same as Fig. A9 but for Kiruna.</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>**Figure A14.** Same as Fig. A9 but for Harestua.</figcaption></figure>
<figure><img src="image_30.png" /><figcaption>**Figure A12.** Same as Fig. A9 but for Thule.</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>**Figure A13.** Same as Fig. A9 but for Kiruna.</figcaption></figure>

<figure><img src="image_32.png" /><figcaption>**Figure A14.** Same as Fig. A9 but for Harestua.</figcaption></figure>
## Appendix B: Additional figures for CO

<figure><img src="image_32.png" /><figcaption>**Figure B1.** (a) FTIR (black) and smoothed model 0–7 km partial columns of CO by day of year, from Eureka. Model data are the nearest in time to each FTIR measurement. (b) Model–measurement percent difference ($\Delta_i$) from Eq. (2) by day of year. Each year is indicated by a different marker.</figcaption></figure>

<figure><img src="image_33.png" /><figcaption>**Figure B2.** Same as Fig. B1 but for Ny Ålesund.</figcaption></figure>
## Appendix B: Additional figures for CO

<figure><img src="image_33.png" /><figcaption>**Figure B1.** (a) FTIR (black) and smoothed model 0–7 km partial columns of CO by day of year, from Eureka. Model data are the nearest in time to each FTIR measurement. (b) Model–measurement percent difference ($\Delta_i$) from Eq. (2) by day of year. Each year is indicated by a different marker.</figcaption></figure>

<figure><img src="image_34.png" /><figcaption>**Figure B2.** Same as Fig. B1 but for Ny Ålesund.</figcaption></figure>
<figure><img src="image_34.png" /><figcaption>Figure B3. Same as Fig. B1 but for Thule.</figcaption></figure>

<figure><img src="image_35.png" /><figcaption>Figure B4. Same as Fig. B1 but for Kiruna.</figcaption></figure>
<figure><img src="image_35.png" /><figcaption>Figure B3. Same as Fig. B1 but for Thule.</figcaption></figure>

<figure><img src="image_36.png" /><figcaption>Figure B4. Same as Fig. B1 but for Kiruna.</figcaption></figure>
<figure><img src="image_36.png" /></figure>

**Figure B5.** (a) Monthly mean FTIR (black) and smoothed model 0–7 km partial columns of CO (PC<sub>F,monthly,j</sub> and PC<sub>M,monthly,j</sub>, respectively), from Eureka using model data that are the nearest in time to each FTIR measurement shown in Fig. B1. Error bars represent the standard deviation of the monthly mean. (b) Model-measurement mean percent difference by month (Δ<sub>monthly,j</sub>). Error bars represent standard deviation of the monthly mean percent difference. The legend in panel (b) shows the overall mean percent difference (Δ<sub>O</sub>) with the standard deviation of the overall mean percent difference.

<figure><img src="image_37.png" /></figure>

**Figure B6.** Same as Fig. B5 but for Ny Ålesund.

<figure><img src="image_38.png" /></figure>

**Figure B7.** Same as Fig. B5 but for Thule.
<figure><img src="image_37.png" /></figure>

**Figure B5.** (a) Monthly mean FTIR (black) and smoothed model 0–7 km partial columns of CO (PC<sub>F,monthly,j</sub> and PC<sub>M,monthly,j</sub>, respectively), from Eureka using model data that are the nearest in time to each FTIR measurement shown in Fig. B1. Error bars represent the standard deviation of the monthly mean. (b) Model-measurement mean percent difference by month (Δ<sub>monthly,j</sub>). Error bars represent standard deviation of the monthly mean percent difference. The legend in panel (b) shows the overall mean percent difference (Δ<sub>O</sub>) with the standard deviation of the overall mean percent difference.

<figure><img src="image_38.png" /></figure>

**Figure B6.** Same as Fig. B5 but for Ny Ålesund.

<figure><img src="image_39.png" /></figure>

**Figure B7.** Same as Fig. B5 but for Thule.
<figure><img src="image_39.png" /><figcaption>**Figure B8.** Same as Fig. B5 but for Kiruna.</figcaption></figure>

<figure><img src="image_40.png" /><figcaption>**Figure B9.** Smoothed model vs. FTIR 0–7 km partial columns of CO for Ny Ålesund, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and $R^2$ are noted in the legend. The 1:1 line is shown in light grey.</figcaption></figure>
<figure><img src="image_40.png" /><figcaption>**Figure B8.** Same as Fig. B5 but for Kiruna.</figcaption></figure>

<figure><img src="image_41.png" /><figcaption>**Figure B9.** Smoothed model vs. FTIR 0–7 km partial columns of CO for Ny Ålesund, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and $R^2$ are noted in the legend. The 1:1 line is shown in light grey.</figcaption></figure>
<figure><img src="image_41.png" /><figcaption>**Figure B10.** Same as Fig. B9 but for Thule.</figcaption></figure>

<figure><img src="image_42.png" /><figcaption>**Figure B11.** Same as Fig. B9 but for Kiruna.</figcaption></figure>
<figure><img src="image_42.png" /><figcaption>**Figure B10.** Same as Fig. B9 but for Thule.</figcaption></figure>

<figure><img src="image_43.png" /><figcaption>**Figure B11.** Same as Fig. B9 but for Kiruna.</figcaption></figure>
## Appendix C: Additional figures for O₃

<figure><img src="image_43.png" /><figcaption>**Figure C1.** (a) FTIR (black) and smoothed model partial columns of O₃ by day of year, from Eureka. Model data are the nearest in time to each FTIR measurement. (b) Model-measurement percent difference ($\Delta_i$) from Eq. (2) by day of year. Each year is indicated by a different marker.</figcaption></figure>

<figure><img src="image_44.png" /><figcaption>**Figure C2.** Same as Fig. C1 but for Ny Ålesund.</figcaption></figure>
## Appendix C: Additional figures for O₃

<figure><img src="image_44.png" /><figcaption>**Figure C1.** (a) FTIR (black) and smoothed model partial columns of O₃ by day of year, from Eureka. Model data are the nearest in time to each FTIR measurement. (b) Model-measurement percent difference ($\Delta_i$) from Eq. (2) by day of year. Each year is indicated by a different marker.</figcaption></figure>

<figure><img src="image_45.png" /><figcaption>**Figure C2.** Same as Fig. C1 but for Ny Ålesund.</figcaption></figure>
<figure><img src="image_45.png" /><figcaption>Figure C3. Same as Fig. C1 but for Thule.</figcaption></figure>

<figure><img src="image_46.png" /><figcaption>Figure C4. Same as Fig. C1 but for Kiruna.</figcaption></figure>
<figure><img src="image_46.png" /><figcaption>Figure C3. Same as Fig. C1 but for Thule.</figcaption></figure>

<figure><img src="image_47.png" /><figcaption>Figure C4. Same as Fig. C1 but for Kiruna.</figcaption></figure>
<figure><img src="image_47.png" /><figcaption>**Figure C5.** Same as Fig. C1 but for Harestua.</figcaption></figure>

<figure><img src="image_48.png" /><figcaption>**Figure C6.** (a) Monthly mean FTIR (black) and smoothed model 0–7 km partial columns of O₃ (PC<sub>F,monthly,j</sub> and PC<sub>M,monthly,j</sub>, respectively) from Eureka using model data that are the nearest in time to each FTIR measurement shown in Fig. C1. Error bars represent the standard deviation of the monthly mean. (b) Model-measurement mean percent difference by month (Δ<sub>monthly,j</sub>). Error bars represent standard deviation of the monthly mean percent difference. The legend in panel (b) shows the overall mean percent difference (Δ<sub>O</sub>) with the standard deviation of the overall mean percent difference.</figcaption></figure>
<figure><img src="image_48.png" /><figcaption>**Figure C5.** Same as Fig. C1 but for Harestua.</figcaption></figure>

<figure><img src="image_49.png" /><figcaption>**Figure C6.** (a) Monthly mean FTIR (black) and smoothed model 0–7 km partial columns of O₃ (PC<sub>F,monthly,j</sub> and PC<sub>M,monthly,j</sub>, respectively) from Eureka using model data that are the nearest in time to each FTIR measurement shown in Fig. C1. Error bars represent the standard deviation of the monthly mean. (b) Model–measurement mean percent difference by month (Δ<sub>monthly,j</sub>). Error bars represent standard deviation of the monthly mean percent difference. The legend in panel (b) shows the overall mean percent difference (Δ<sub>O</sub>) with the standard deviation of the overall mean percent difference.</figcaption></figure>
<figure><img src="image_49.png" /><figcaption>**Figure C7.** Same as Fig. C6 but for Ny Ålesund.</figcaption></figure>

<figure><img src="image_50.png" /><figcaption>**Figure C8.** Same as Fig. C6 but for Thule.</figcaption></figure>
<figure><img src="image_50.png" /><figcaption>**Figure C7.** Same as Fig. C6 but for Ny Ålesund.</figcaption></figure>

<figure><img src="image_51.png" /><figcaption>**Figure C8.** Same as Fig. C6 but for Thule.</figcaption></figure>
<figure><img src="image_51.png" /><figcaption>**Figure C9.** Same as Fig. C6 but for Kiruna.</figcaption></figure>

<figure><img src="image_52.png" /><figcaption>**Figure C10.** Same as Fig. C6 but for Harestua.</figcaption></figure>
<figure><img src="image_52.png" /><figcaption>**Figure C9.** Same as Fig. C6 but for Kiruna.</figcaption></figure>

<figure><img src="image_53.png" /><figcaption>**Figure C10.** Same as Fig. C6 but for Harestua.</figcaption></figure>
<figure><img src="image_53.png" /><figcaption>**Figure C11.** Model vs. FTIR 0–7 km partial columns of O₃ for Ny Ålesund, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and R² are noted in the legend. The 1 : 1 line is shown in light grey.</figcaption></figure>

<figure><img src="image_54.png" /><figcaption>**Figure C12.** Same as Fig. C11 but for Thule.</figcaption></figure>
<figure><img src="image_54.png" /><figcaption>**Figure C11.** Model vs. FTIR 0–7 km partial columns of O₃ for Ny Ålesund, showing all available model–FTIR corresponding data. The black line is the line of best fit, where the equation and R² are noted in the legend. The 1 : 1 line is shown in light grey.</figcaption></figure>

<figure><img src="image_55.png" /><figcaption>**Figure C12.** Same as Fig. C11 but for Thule.</figcaption></figure>
<figure><img src="image_55.png" /><figcaption>Figure C13. Same as Fig. C11 but for Kiruna.</figcaption></figure>

<figure><img src="image_56.png" /><figcaption>Figure C14. Same as Fig. C11 but for Harestua.</figcaption></figure>

**Data availability.** The FTIR data are publicly available in the NDACC data repository (https://www-air.larc.nasa.gov/missions/ndacc/data.html, NDACC, 2023), and the model data are available at https://doi.org/10.18164/e0a0ac5c-d851-45b9-b6d9-4abc29d7d419 (CCCma, 2023).

**Author contributions.** KAW, CHW, and KS conceived this project. VAF performed the formal analysis, including the comparisons between datasets, formation of plots and tables, and writing the paper. KS, CHW, and KAW provided scientific guidance and support throughout the work in addition to comments and edits to the paper. TB, JWH, JM, JN, MP, and ANR provided advice on the
<figure><img src="image_56.png" /><figcaption>Figure C13. Same as Fig. C11 but for Kiruna.</figcaption></figure>

<figure><img src="image_57.png" /><figcaption>Figure C14. Same as Fig. C11 but for Harestua.</figcaption></figure>

**Data availability.** The FTIR data are publicly available in the NDACC data repository (https://www-air.larc.nasa.gov/missions/ndacc/data.html, NDACC, 2023), and the model data are available at https://doi.org/10.18164/e0a0ac5c-d851-45b9-b6d9-4abc29d7d419 (CCCma, 2023).

**Author contributions.** KAW, CHW, and KS conceived this project. VAF performed the formal analysis, including the comparisons between datasets, formation of plots and tables, and writing the paper. KS, CHW, and KAW provided scientific guidance and support throughout the work in addition to comments and edits to the paper. TB, JWH, JM, JN, MP, and ANR provided advice on the
FTIR data and facilitated the operations and data management for their respective instruments. SA, SB, RYC, JC, MD, SD, XD, JSF, MG, WG, JL, KSL, LM, TO, NO, DAP, LP, JCR, MAT, SvT, and StT provided the model outputs for the AMAP report and guidance on the analysis of models. All co-authors provided feedback on the paper.

**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.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

**Acknowledgements.** The authors acknowledge the use of model datasets from the 2021 AMAP SLCF assessment report. The authors would like to thank all the people involved with the data collection and maintenance of the high-latitude NDACC FTIR instruments used in this study.

The Eureka FTIR measurements were made at the Polar Environment Atmospheric Research Laboratory (PEARL) by the Canadian Network for the Detection of Atmospheric Composition Change (CANDAC), which has been supported by the Atlantic Innovation Fund/Nova Scotia Research Innovation Trust, the Canada Foundation for Innovation, the Canadian Foundation for Climate and Atmospheric Sciences, the Canadian Space Agency, Environment and Climate Change Canada (ECCC), Government of Canada International Polar Year funding, the Natural Sciences and Engineering Research Council, the Northern Scientific Training Program, the Ontario Innovation Trust, the Polar Continental Shelf Program, and the Ontario Research Fund. We thank former CANDAC/PEARL PI James Drummond, PEARL Site Manager Pierre Fogal, CANDAC Data Manager Yan Tsehtik, the CANDAC operators, and the staff at ECCC's Eureka Weather Station for their contributions to data acquisition and for logistical and on-site support.

We gratefully acknowledge funding from the Transregional Collaborative Research Centre TR 172 - Arctic Amplification: Climate Relevant Atmospheric and Surface Processes (AC)³, project E02: Ny-Ålesund Column Thermodynamic Structure, Clouds, Aerosols, Trace Gases and Radiative Effects. We also thank the Senate of Bremer for financial support. The AWI Bremerhaven provided logistical support for measurements in Ny Ålesund.

Karlsruhe Institute of Technology would like to thank Uwe Rafalski from the Swedish Institute of Space Physics (IRF) for their continuing support of the NDACC FTIR site Kiruna.

Michael Gauss and Svetlana Tsyro received financial support from the Arctic Monitoring and Assessment Programme (AMAP).

Kathy S. Law, Jean-Christophe Raut, Louis Marelle, and Tatsuo Onishi (LATMOS) acknowledge support from 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) and for access to IDRIS HPC resources (GENCI allocation A009017141) as well as the IPSL mesoscale computing center (CICLAD: Calcul Intensif pour le CLimat, l'Atmosphère et la Dynamique) for model simulations. Kathy S. Law also acknowledges support from the French Space Agency (CNES) MERLIN (contract no. 7752).

Makoto Deushi and Naga Oshima were supported by the Environment Research and Technology Development Fund (grant nos. JPMEERF20202003 and JPMEERF20232001) of the Environmental Restoration and Conservation Agency Provided by the Ministry of Environment of Japan, the Arctic Challenge for Sustainability II (ArCS II), program grant no. JPMXD1420318865, and a grant for the Global Environmental Research Coordination System from the Ministry of the Environment, Japan (grant no. MLIT2253).

**Financial support.** Funding for this work was provided by the Canadian Space Agency under the Earth System Science Data Analyses Program (grant no. 21SUASABBC).

**Review statement.** This paper was edited by Manvendra Krishna Dubey and reviewed by two anonymous referees.

## References

Amann, M., Bertok, I., Borken-Kleefled, J., Cofala, J., Heyes, C., Höglund-Isaksson, L., Klimont, Z., Nguyen, B., Posch, M., Rafaj, P., Sandler, R., Schöpp, W., Wagner, F., and Winiwarter, W.: Cost-effective control of air quality and greenhouse gases in Europe: Modelling and policy applications, *Environ. Modell. Softw.*, 26, 1489–1501, https://doi.org/10.1016/j.envsoft.2011.07.012, 2011.

AMAP: AMAP Assessment 2021: Impacts of Short-lived Climate Forcers on Arctic Climate, Air Quality, and Human Health, Arctic Monitoring and Assessment Programme (AMAP), Tromsø, Norway, viii + 324 pp., https://www.amap.no/documents/doc/amap-assessment-2021-impacts-of-short-lived-climate-forcers-on-arctic-climate-air-quality-and-human-health/3614(last access: 14 February 2023), 2021.

Andersson, C., Langner, J., and Bergström, R.: Interannual variation and trends in air pollution over Europe due to climate variability during 1958–2001 simulated with a regional CTM coupled to the ERA40 reanalysis, *Tellus B*, 59, 77–98, https://doi.org/10.1111/j.1600-0889.2006.00231.x, 2007.

Ballinger, T. J., Overland, J. E., Wang, M., Bhatt, U. S., Hanna, E., Hanssen-Bauer, I., Kim, S.-J., Thoman, R. L., and Walsh, J. E.: Arctic Report Card 2020: Surface Air Temperature, Tech. rep., National Oceanic and Atmospheric Administration (NOAA), Office of Oceanic and Atmospheric Research, Pacific Marine Environmental Laboratory (U.S.), https://doi.org/10.25923/gcw8-2z06, 2020.

Batchelor, R. L., Strong, K., Lindenmaier, R., Mittermeier, R. L., Fast, H., Drummond, J. R., and Fogal, P. F.: A new Bruker IFS 125HR FTIR spectrometer for the Polar Environment Atmospheric Research Laboratory at Eureka, Nunavut, Canada: measurements and comparison with the existing Bomem
FTIR data and facilitated the operations and data management for their respective instruments. SA, SB, RYC, JC, MD, SD, XD, JSF, MG, WG, JL, KSL, LM, TO, NO, DAP, LP, JCR, MAT, SvT, and StT provided the model outputs for the AMAP report and guidance on the analysis of models. All co-authors provided feedback on the paper.

**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.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

**Acknowledgements.** The authors acknowledge the use of model datasets from the 2021 AMAP SLCF assessment report. The authors would like to thank all the people involved with the data collection and maintenance of the high-latitude NDACC FTIR instruments used in this study.

The Eureka FTIR measurements were made at the Polar Environment Atmospheric Research Laboratory (PEARL) by the Canadian Network for the Detection of Atmospheric Composition Change (CANDAC), which has been supported by the Atlantic Innovation Fund/Nova Scotia Research Innovation Trust, the Canada Foundation for Innovation, the Canadian Foundation for Climate and Atmospheric Sciences, the Canadian Space Agency, Environment and Climate Change Canada (ECCC), Government of Canada International Polar Year funding, the Natural Sciences and Engineering Research Council, the Northern Scientific Training Program, the Ontario Innovation Trust, the Polar Continental Shelf Program, and the Ontario Research Fund. We thank former CANDAC/PEARL PI James Drummond, PEARL Site Manager Pierre Fogal, CANDAC Data Manager Yan Tsehtik, the CANDAC operators, and the staff at ECCC's Eureka Weather Station for their contributions to data acquisition and for logistical and on-site support.

We gratefully acknowledge funding from the Transregional Collaborative Research Centre TR 172 - Arctic Amplification: Climate Relevant Atmospheric and Surface Processes (AC)³, project E02: Ny-Ålesund Column Thermodynamic Structure, Clouds, Aerosols, Trace Gases and Radiative Effects. We also thank the Senate of Bremeren for financial support. The AWI Bremerhaven provided logistical support for measurements in Ny Ålesund.

Karlsruhe Institute of Technology would like to thank Uwe Rafalski from the Swedish Institute of Space Physics (IRF) for their continuing support of the NDACC FTIR site Kiruna.

Michael Gauss and Svetlana Tsyro received financial support from the Arctic Monitoring and Assessment Programme (AMAP).

Kathy S. Law, Jean-Christophe Raut, Louis Marelle, and Tatsuo Onishi (LATMOS) acknowledge support from 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) and for access to IDRIS HPC resources (GENCI allocation A009017141) as well as the IPSL mesoscale computing center (CICLAD: Calcul Intensif pour le CLimat, l'Atmosphère et la Dynamique) for model simulations. Kathy S. Law also acknowledges support from the French Space Agency (CNES) MERLIN (contract no. 7752).

Makoto Deushi and Naga Oshima were supported by the Environment Research and Technology Development Fund (grant nos. JPMEERF20202003 and JPMEERF20232001) of the Environmental Restoration and Conservation Agency Provided by the Ministry of Environment of Japan, the Arctic Challenge for Sustainability II (ArCS II), program grant no. JPMXD1420318865, and a grant for the Global Environmental Research Coordination System from the Ministry of the Environment, Japan (grant no. MLIT2253).

**Financial support.** Funding for this work was provided by the Canadian Space Agency under the Earth System Science Data Analyses Program (grant no. 21SUASABBC).

**Review statement.** This paper was edited by Manvendra Krishna Dubey and reviewed by two anonymous referees.

## References

Amann, M., Bertok, I., Borken-Kleefled, J., Cofala, J., Heyes, C., Höglund-Isaksson, L., Klimont, Z., Nguyen, B., Posch, M., Rafaj, P., Sandler, R., Schöpp, W., Wagner, F., and Winiwarter, W.: Cost-effective control of air quality and greenhouse gases in Europe: Modelling and policy applications, *Environ. Modell. Softw.*, 26, 1489–1501, https://doi.org/10.1016/j.envsoft.2011.07.012, 2011.

AMAP: AMAP Assessment 2021: Impacts of Short-lived Climate Forcers on Arctic Climate, Air Quality, and Human Health, Arctic Monitoring and Assessment Programme (AMAP), Tromsø, Norway, viii + 324 pp., https://www.amap.no/documents/doc/amap-assessment-2021-impacts-of-short-lived-climate-forcers-on-arctic-climate-air-quality-and-human-health/3614(last access: 14 February 2023), 2021.

Andersson, C., Langner, J., and Bergström, R.: Interannual variation and trends in air pollution over Europe due to climate variability during 1958–2001 simulated with a regional CTM coupled to the ERA40 reanalysis, *Tellus B*, 59, 77–98, https://doi.org/10.1111/j.1600-0889.2006.00231.x, 2007.

Ballinger, T. J., Overland, J. E., Wang, M., Bhatt, U. S., Hanna, E., Hanssen-Bauer, I., Kim, S.-J., Thoman, R. L., and Walsh, J. E.: Arctic Report Card 2020: Surface Air Temperature, Tech. rep., National Oceanic and Atmospheric Administration (NOAA), Office of Oceanic and Atmospheric Research, Pacific Marine Environmental Laboratory (U.S.), https://doi.org/10.25923/gcw8-2z06, 2020.

Batchelor, R. L., Strong, K., Lindenmaier, R., Mittermeier, R. L., Fast, H., Drummond, J. R., and Fogal, P. F.: A new Bruker IFS 125HR FTIR spectrometer for the Polar Environment Atmospheric Research Laboratory at Eureka, Nunavut, Canada: measurements and comparison with the existing Bomem
DA8 spectrometer, J. Atmos. Ocean. Tech., 26, 1328–1340,
https://doi.org/10.1175/2009JTECHA1215.1, 2009.

Berg, T., Sekkesaeter, S., Steinnes, E., Valdal, A. K., and Wibetoe, G.: Springtime depletion of mercury in the European Arctic as observed at Svalbard, Sci. Total Environ., 304, 43–51, https://doi.org/10.1016/S0048-9697(02)00555-7, 2003.

Bey, I., Jacob, D. J., Yantosca, R. M., Logan, J. A., Field, B. D., Fiore, A. M., Li, Q., Liu, H. Y., Mickley, L. J., and Schultz, M. G.: Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation, J. Geophys. Res., 106, 23073–23095, https://doi.org/10.1029/2001JD000807, 2001.

Blumenstock, T., Fisher, H., Friedle, A., Hase, F., and Thomas, P.: Column Amounts of ClONO₂, HCl, HNO₃, and HF from Ground-Based FTIR Measurements Made Near Kiruna, Sweden, in Late Winter 1994, J. Atmos. Chem., 26, 311–321, https://doi.org/10.1023/A:1005728713762, 1997.

Blumenstock, T., Hase, F., Kramer, I., Mikuteit, S., Fischer, H., Goutail, F., and Raffalski, U.: Winter to winter variability of chlorine activation and ozone loss as observed by ground-based FTIR measurements at Kiruna since winter 1993/94, Int. J. Remote Sens., 30, 4055–4064, https://doi.org/10.1080/01431160902821916, 2009.

Brandt, J., Silver, J., Frohn, L. M., Geels, C., Gross, A., Hansen, A. B., Hansen, K. M., Hedegaard, G. B., Skjøth, C. A., Villadsen, H., Zare, A., and Christensen, J. H.: An integrated model study for Europe and North America using the Danish Eulerian Hemispheric Model with focus on intercontinental transport of air pollution, Atmos. Environ., 53, 156–176, https://doi.org/10.1016/j.atmosenv.2012.01.011, 2012.

Bush, E. and Lemmen, D. S. (Eds.): Canada's Changing Climate Report, Government of Canada, Ottawa, ON, 444 p., https://changingclimate.ca/CCCR2019/ (last access: 2 February 2023), 2019.

CCCma (Canadian Centre for Climate Modelling and analysis): AMAP SLCF models output in NetCDF format, CCCma [data set], https://doi.org/10.18164/e0a0ac5c-d851-45b9-b6d9-4abc29d7d419, 2023.

Christensen, J. H.: The Danish Eulerian hemispheric model – A three-dimensional air pollution model used for the Arctic, Atmos. Environ., 31, 4169–4191, https://doi.org/10.1016/S1352-2310(97)00264-1, 1997.

Danabasoglu, G., Lamarque, J., Bacmeister, J., Bailey, D. A., DuVivier, A. K., and Edwards, J.: The Community Earth System Model Version 2 (CESM2), J. Adv. Model. Earth Sy., 12, e2019MS001916, https://doi.org/10.1029/2019MS001916, 2020.

De Mazière, M., Thompson, A. M., Kurylo, M. J., Wild, J. D., Bernhard, G., Blumenstock, T., Braathen, G. O., Hannigan, J. W., Lambert, J.-C., Leblanc, T., McGee, T. J., Nedoluha, G., Petropavlovskikh, I., Seckmeyer, G., Simon, P. C., Steinbrecht, W., and Strahan, S. E.: The Network for the Detection of Atmospheric Composition Change (NDACC): history, status and perspectives, Atmos. Chem. Phys., 18, 4935–4964, https://doi.org/10.5194/acp-18-4935-2018, 2018.

Emmons, L. K., Walters, S., Hess, P. G., Lamarque, J.-F., Pfister, G. G., Fillmore, D., Granier, C., Guenther, A., Kinnison, D., Laepple, T., Orlando, J., Tie, X., Tyndall, G., Wiedinmyer, C., Baughcum, S. L., and Kloster, S.: Description and evaluation of the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4), Geosci. Model Dev., 3, 43–67, https://doi.org/10.5194/gmd-3-43-2010, 2010.

Emmons, L. K., Arnold, S. R., Monks, S. A., Huijnen, V., Tilmes, S., Law, K. S., Thomas, J. L., Raut, J.-C., Bouarar, I., Turquety, S., Long, Y., Duncan, B., Steenrod, S., Strode, S., Fleming, J., Mao, J., Langner, J., Thompson, A. M., Tarasick, D., Apel, E. C., Blake, D. R., Cohen, R. C., Dibb, J., Diskin, G. S., Fried, A., Hall, S. R., Huey, L. G., Weinheimer, A. J., Wisthaler, A., Mikoviny, T., Nowak, J., Peischl, J., Roberts, J. M., Ryerson, T., Warneke, C., and Helmig, D.: The POLARCAT Model Intercomparison Project (POLMIP): overview and evaluation with observations, Atmos. Chem. Phys., 15, 6721–6744, https://doi.org/10.5194/acp-15-6721-2015, 2015.

Eneroth, K., Holmén, K., Berg, T., Schmidbauer, N., and Solberg, S.: Springtime depletion of tropospheric ozone, gaseous elemental mercury and non-methane hydrocarbons in the European Arctic, and its relation to atmospheric transport, Atmos. Environ., 41, 8511–8526, https://doi.org/10.1016/j.atmosenv.2007.07.008, 2007.

Fisher, J. A., Jacob, D. J., Purdy, M. T., Kopacz, M., Le Sager, P., Carouge, C., Holmes, C. D., Yantosca, R. M., Batchelor, R. L., Strong, K., Diskin, G. S., Fuelberg, H. E., Holloway, J. S., Hyer, E. J., McMillan, W. W., Warner, J., Streets, D. G., Zhang, Q., Wang, Y., and Wu, S.: Source attribution and interannual variability of Arctic pollution in spring constrained by aircraft (ARCTAS, ARCPAC) and satellite (AIRS) observations of carbon monoxide, Atmos. Chem. Phys., 10, 977–996, https://doi.org/10.5194/acp-10-977-2010, 2010.

Galle, B., Mellqvist, J., Arlander, D. W., Fløisand, I., Chipperfield, M. P., and Lee, A. M.: Ground Based FTIR Measurements of stratospheric species from Harestua, Norway during SESAME and Comparison with Models, J. Atmos. Chem., 32, 147–164, https://doi.org/10.1023/A:1006137924562, 1999.

Gong, W., Makar, P. A., Zhang, J., Milbrandt, J., Gravel, S., Hayden, K. L., Macdonald, A. M., and Leaitch, W. R.: Modelling aerosol cloud meteorology interaction: A case study with a fully coupled air quality model GEM-MACH, Atmos. Environ., 115, 695–715, https://doi.org/10.1016/j.atmosenv.2015.05.062, 2015.

Gong, W., Beagley, S. R., Cousineau, S., Sassi, M., Munoz-Alpizar, R., Ménard, S., Racine, J., Zhang, J., Chen, J., Morrison, H., Sharma, S., Huang, L., Bellavance, P., Ly, J., Izdebski, P., Lyons, L., and Holt, R.: Assessing the impact of shipping emissions on air pollution in the Canadian Arctic and northern regions: current and future modelled scenarios, Atmos. Chem. Phys., 18, 16653–16687, https://doi.org/10.5194/acp-18-16653-2018, 2018.

Hannigan, J. W., Coffey, M. T., and Goldman, A.: Semiautonomous FTS Observation System for Remote Sensing of Stratospheric and Tropospheric Gases, J. Atmos. Ocean. Tech., 26, 1814–1828, https://doi.org/10.1175/2009JTECHA1230.1, 2009.

Hase, F., Hannigan, J. W., Coffey, M. T., Goldman, A., Höpfner, M., Jones, N. B., Rinsland, C. P., and Wood, S. W.: Intercomparison of retrieval codes used for the analysis of high-resolution, ground-based FTIR measurements, J. Quant. Spectrosc. Ra., 87, 25–52, https://doi.org/10.1016/j.jqsrt.2003.12.008, 2004.

Hirdman, D., Sodemann, H., Eckhardt, S., Burkhart, J. F., Jeffer- son, A., Mefford, T., Quinn, P. K., Sharma, S., Ström, J., and Stohl, A.: Source identification of short-lived air pollutants in the Arctic using statistical analysis of measurement data and par-
DA8 spectrometer, J. Atmos. Ocean. Tech., 26, 1328–1340,
https://doi.org/10.1175/2009JTECHA1215.1, 2009.

Berg, T., Sekkesaeter, S., Steinnes, E., Valdal, A. K., and Wibetoe, G.: Springtime depletion of mercury in the European Arctic as observed at Svalbard, Sci. Total Environ., 304, 43–51, https://doi.org/10.1016/S0048-9697(02)00555-7, 2003.

Bey, I., Jacob, D. J., Yantosca, R. M., Logan, J. A., Field, B. D., Fiore, A. M., Li, Q., Liu, H. Y., Mickley, L. J., and Schultz, M. G.: Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation, J. Geophys. Res., 106, 23073–23095, https://doi.org/10.1029/2001JD000807, 2001.

Blumenstock, T., Fisher, H., Friedle, A., Hase, F., and Thomas, P.: Column Amounts of ClONO₂, HCl, HNO₃, and HF from Ground-Based FTIR Measurements Made Near Kiruna, Sweden, in Late Winter 1994, J. Atmos. Chem., 26, 311–321, https://doi.org/10.1023/A:1005728713762, 1997.

Blumenstock, T., Hase, F., Kramer, I., Mikuteit, S., Fischer, H., Goutail, F., and Raffalski, U.: Winter to winter variability of chlorine activation and ozone loss as observed by ground-based FTIR measurements at Kiruna since winter 1993/94, Int. J. Remote Sens., 30, 4055–4064, https://doi.org/10.1080/01431160902821916, 2009.

Brandt, J., Silver, J., Frohn, L. M., Geels, C., Gross, A., Hansen, A. B., Hansen, K. M., Hedegaard, G. B., Skjøth, C. A., Villadsen, H., Zare, A., and Christensen, J. H.: An integrated model study for Europe and North America using the Danish Eulerian Hemispheric Model with focus on intercontinental air transport of air pollution, Atmos. Environ., 53, 156–176, https://doi.org/10.1016/j.atmosenv.2012.01.011, 2012.

Bush, E. and Lemmen, D. S. (Eds.): Canada's Changing Climate Report, Government of Canada, Ottawa, ON, 444 p., https://changingclimate.ca/CCCR2019/ (last access: 2 February 2023), 2019.

CCCma (Canadian Centre for Climate Modelling and analysis): AMAP SLCF models output in NetCDF format, CCCma [data set], https://doi.org/10.18164/e0a0ac5c-d851-45b9-b6d9-4abc29d7d419, 2023.

Christensen, J. H.: The Danish Eulerian hemispheric model – A three-dimensional air pollution model used for the Arctic, Atmos. Environ., 31, 4169–4191, https://doi.org/10.1016/S1352-2310(97)00264-1, 1997.

Danabasoglu, G., Lamarque, J., Bacmeister, J., Bailey, D. A., DuVivier, A. K., and Edwards, J.: The Community Earth System Model Version 2 (CESM2), J. Adv. Model. Earth Sy., 12, e2019MS001916, https://doi.org/10.1029/2019MS001916, 2020.

De Mazière, M., Thompson, A. M., Kurylo, M. J., Wild, J. D., Bernhard, G., Blumenstock, T., Braathen, G. O., Hannigan, J. W., Lambert, J.-C., Leblanc, T., McGee, T. J., Nedoluha, G., Petropavlovskikh, I., Seckmeyer, G., Simon, P. C., Steinbrecht, W., and Strahan, S. E.: The Network for the Detection of Atmospheric Composition Change (NDACC): history, status and perspectives, Atmos. Chem. Phys., 18, 4935–4964, https://doi.org/10.5194/acp-18-4935-2018, 2018.

Emmons, L. K., Walters, S., Hess, P. G., Lamarque, J.-F., Pfister, G. G., Fillmore, D., Granier, C., Guenther, A., Kinnison, D., Laepple, T., Orlando, J., Tie, X., Tyndall, G., Wiedinmyer, C., Baughcum, S. L., and Kloster, S.: Description and evaluation of the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4), Geosci. Model Dev., 3, 43–67, https://doi.org/10.5194/gmd-3-43-2010, 2010.

Emmons, L. K., Arnold, S. R., Monks, S. A., Huijnen, V., Tilmes, S., Law, K. S., Thomas, J. L., Raut, J.-C., Bouarar, I., Turquety, S., Long, Y., Duncan, B., Steenrod, S., Strode, S., Fleming, J., Mao, J., Langner, J., Thompson, A. M., Tarasick, D., Apel, E. C., Blake, D. R., Cohen, R. C., Dibb, J., Diskin, G. S., Fried, A., Hall, S. R., Huey, L. G., Weinheimer, A. J., Wisthaler, A., Mikoviny, T., Nowak, J., Peischl, J., Roberts, J. M., Ryerson, T., Warneke, C., and Helmig, D.: The POLARCAT Model Intercomparison Project (POLMIP): overview and evaluation with observations, Atmos. Chem. Phys., 15, 6721–6744, https://doi.org/10.5194/acp-15-6721-2015, 2015.

Eneroth, K., Holmén, K., Berg, T., Schmidbauer, N., and Solberg, S.: Springtime depletion of tropospheric ozone, gaseous elemental mercury and non-methane hydrocarbons in the European Arctic, and its relation to atmospheric transport, Atmos. Environ., 41, 8511–8526, https://doi.org/10.1016/j.atmosenv.2007.07.008, 2007.

Fisher, J. A., Jacob, D. J., Purdy, M. T., Kopacz, M., Le Sager, P., Carouge, C., Holmes, C. D., Yantosca, R. M., Batchelor, R. L., Strong, K., Diskin, G. S., Fuelberg, H. E., Holloway, J. S., Hyer, E. J., McMillan, W. W., Warner, J., Streets, D. G., Zhang, Q., Wang, Y., and Wu, S.: Source attribution and interannual variability of Arctic pollution in spring constrained by aircraft (ARCTAS, ARCPAC) and satellite (AIRS) observations of carbon monoxide, Atmos. Chem. Phys., 10, 977–996, https://doi.org/10.5194/acp-10-977-2010, 2010.

Galle, B., Mellqvist, J., Arlander, D. W., Flöisand, I., Chipperfield, M. P., and Lee, A. M.: Ground Based FTIR Measurements of stratospheric species from Harestua, Norway during SESAME and Comparison with Models, J. Atmos. Chem., 32, 147–164, https://doi.org/10.1023/A:1006137924562, 1999.

Gong, W., Makar, P. A., Zhang, J., Milbrandt, J., Gravel, S., Hayden, K. L., Macdonald, A. M., and Leaitch, W. R.: Modelling aerosol cloud meteorology interaction: A case study with a fully coupled air quality model GEM-MACH, Atmos. Environ., 115, 695–715, https://doi.org/10.1016/j.atmosenv.2015.05.062, 2015.

Gong, W., Beagley, S. R., Cousineau, S., Sassi, M., Munoz-Alpizar, R., Ménard, S., Racine, J., Zhang, J., Chen, J., Morrison, H., Sharma, S., Huang, L., Bellavance, P., Ly, J., Izdebski, P., Lyons, L., and Holt, R.: Assessing the impact of shipping emissions on air pollution in the Canadian Arctic and northern regions: current and future modelled scenarios, Atmos. Chem. Phys., 18, 16653–16687, https://doi.org/10.5194/acp-18-16653-2018, 2018.

Hannigan, J. W., Coffey, M. T., and Goldman, A.: Semiautonomous FTS Observation System for Remote Sensing of Stratospheric and Tropospheric Gases, J. Atmos. Ocean. Tech., 26, 1814–1828, https://doi.org/10.1175/2009JTECHA1230.1, 2009.

Hase, F., Hannigan, J. W., Coffey, M. T., Goldman, A., Höpfner, M., Jones, N. B., Rinsland, C. P., and Wood, S. W.: Intercomparison of retrieval codes used for the analysis of high-resolution, ground-based FTIR measurements, J. Quant. Spectrosc. Ra., 87, 25–52, https://doi.org/10.1016/j.jqsrt.2003.12.008, 2004.

Hirdman, D., Sodemann, H., Eckhardt, S., Burkhart, J. F., Jeffer- son, A., Mefford, T., Quinn, P. K., Sharma, S., Ström, J., and Stohl, A.: Source identification of short-lived air pollutants in the Arctic using statistical analysis of measurement data and par-
ticle dispersion model output, Atmos. Chem. Phys., 10, 669–693,
https://doi.org/10.5194/acp-10-669-2010, 2010.

Höglund-Isaksson, L., Gómez-Sanabria, A., Klimont, Z., Rafaj, P., and Schöpp, W.: Technical potentials and costs for reducing global anthropogenic methane emissions in the 2050 timeframe – results from the GAINS model, Environmental Research Communications, 2, 025004, https://doi.org/10.1088/2515-7620/ab7457, 2020.

IPCC: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Conners, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 2391 pp., https://doi.org/10.1017/9781009157896, 2021.

IRWG: Infrared Working Group Retrieval Code, SFIT – SFIT Spectral Data Analysis Model, https://wiki.ucar.edu/display/sfit4/ (last access: 2 February 2023), 2020.

Jonsson, A. I., de Grandpré, J., Fomichev, V. I., McConnell, J. C., and Beagley, S. R.: Doubled CO₂-induced cooling in the middle atmosphere: Photochemical analysis of the ozone radiative feedback, J. Geophys. Res., 109, D24103, https://doi.org/10.1029/2004JD005093, 2004.

Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., Saha, S., White, G., Woollen, J. and Zhu, Y., The NCEP/NCAR 40-year reanalysis project, B. Am. Meteorol. Soc., 77, 437–471, https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2, 1996.

Kärnä, T. and Baptista, A. M.: Evaluation of a long-term hindcast simulation for the Columbia River estuary, Ocean Model., 99, 1–14, https://doi.org/10.1016/j.ocemod.2015.12.007, 2016.

Kawai, H., Yukimoto, S., Koshiro, T., Oshima, N., Tanaka, T., Yoshimura, H., and Nagasawa, R.: Significant improvement of cloud representation in the global climate model MRI-ESM2, Geosci. Model Dev., 12, 2875–2897, https://doi.org/10.5194/gmd-12-2875-2019, 2019.

Klimont, Z., Kupiainen, K., Heyes, C., Purohit, P., Cofala, J., Rafaj, P., Borken-Kleefeld, J., and Schöpp, W.: Global anthropogenic emissions of particulate matter including black carbon, Atmos. Chem. Phys., 17, 8681–8723, https://doi.org/10.5194/acp-17-8681-2017, 2017.

Kokkola, H., Korhonen, H., Lehtinen, K. E. J., Makkonen, R., Asmi, A., Järvenoja, S., Anttila, T., Partanen, A.-I., Kulmala, M., Järvinen, H., Laaksonen, A., and Kerminen, V.-M.: SALSA – a Sectional Aerosol module for Large Scale Applications, Atmos. Chem. Phys., 8, 2469–2483, https://doi.org/10.5194/acp-8-2469-2008, 2008.

Kuhlbrodt, T., Jones, C. G., Sellar, A., Storkey, D., Blockley, E., Stringer, M., Hill, R., Graham, T., Ridley, J., Blaker, A., Calvert, D., Copsey, D., Ellis, R., Hewitt, H., Hyder, P., Ineson, S., Mulcahy, J., Siahaan, A., and Walton, J.: The Low-Resolution Version of HadGEM3 GC3.1: Development and Evaluation for Global Climate, J. Adv. Model. Earth Syst., 10, 2865–2888, https://doi.org/10.1029/2018MS001370, 2018.

Liu, X., Ma, P.-L., Wang, H., Tilmes, S., Singh, B., Easter, R. C., Ghan, S. J., and Rasch, P. J.: Description and evaluation of a

new four-mode version of the Modal Aerosol Module (MAM4) within version 5.3 of the Community Atmosphere Model, Geosci. Model Dev., 9, 505–522, https://doi.org/10.5194/gmd-9-505-2016, 2016.

Lutsch, E., Strong, K., Jones, D. B. A., Blumenstock, T., Conway, S., Fisher, J. A., Hannigan, J. W., Hase, F., Kasai, Y., Mahieu, E., Makarova, M., Morino, I., Nagahama, T., Notholt, J., Ortega, I., Palm, M., Poberovskii, A. V., Sussmann, R., and Warneke, T.: Detection and attribution of wildfire pollution in the Arctic and northern midlatitudes using a network of Fourier-transform infrared spectrometers and GEOS-Chem, Atmos. Chem. Phys., 20, 12813–12851, https://doi.org/10.5194/acp-20-12813-2020, 2020.

Mahieu, E., Fischer, E. V., Franco, B., Palm, M., Wizenberg, T., Smale, D., Clarisse, L., Clerbaux, C., Coheur, P-F., Hannigan, J. W., Lutsch, E., Notholt, J., Pardo Cantos, I., Prignon, M., Servais, C., and Strong, K.: First retrievals of peroxy-acetyl nitrate (PAN) from ground-based FTIR solar spectra recorded at remote sites, comparison with model and satellite data, Elementa: Science of the Anthropocene, 9, 00027, https://doi.org/10.1525/elementa.2021.00027, 2021.

Makar, P. A., Gong, W., Hogrefe, C., Zhang, Y., Curci, G., Zabkar, R., Milbrandt, J., Im, U., Balzarini, A., Baró, R., Bianconi, R., Cheung, P., Forkel, R., Gravel, S., Hirtl, M., Honzak, L., Hou, A., Jiménez-Guerrero, P., Langer, M., Moran, M., Pabla, B., Pérez, J., Pirovano, G., José, R. S., Tuccella, P., Werhahn, J., Zhang, J., and Galmarini, S.: Feedbacks between air pollution and weather, part 2: Effects on chemistry, Atmos. Environ., 115, 499–526, https://doi.org/10.1016/j.atmosenv.2014.10.021, 2015a.

Makar, P. A., Gong, W., Milbrandt, J., Hogrefe, C., Zhang, Y., Curci, G., Žabkar, R., Im, U., Balzarini, A., Baró, R., Bianconi, R., Cheung, P., Forkel, R., Gravel, S., Hirtl, M., Honzak, L., Hou, A., Jiménez-Guerrero, P., Langer, M., Moran, M., Pabla, B., Pérez, J., Pirovano, G., José, R. S., Tuccella, P., Werhahn, J., Zhang, J., and Galmarini, S.: Feedbacks between air pollution and weather, part 1: Effects on weather, Atmos. Environ., 115, 442–469, https://doi.org/10.1016/j.atmosenv.2014.12.003, 2015b.

Marelle, L., Raut, J.-C., Law, K. S., Berg, L. K., Fast, J. D., Easter, R. C., Shrivastava, M., and Thomas, J. L.: Improvements to the WRF-Chem 3.5.1 model for quasi-hemispheric simulations of aerosols and ozone in the Arctic, Geosci. Model Dev., 10, 3661–3677, https://doi.org/10.5194/gmd-10-3661-2017, 2017.

Marelle, L., Raut, J.-C., Law, K. S., and Duclaux, O.: Current and Future Arctic Aerosols and Ozone From Remote Emissions and Emerging Local Sources – Modeled Source Contributions and Radiative Effects, J. Geophys. Res.-Atmos., 123, 12942–12963, https://doi.org/10.1029/2018JD028863, 2018.

Marsh, D. R., Mills, M. J., Kinnison, D. E., Lamarque, J.-F., Calvo, N., and Polvani, L. M.: Climate Change from 1850 to 2005 Simulated in CESM1 (WACCM), J. Climate, 26, 7372–7391, https://doi.org/10.1175/JCLI-D-12-00558.1, 2013.

Massling, A., Nielsen, I. E., Kristensen, D., Christensen, J. H., Sørensen, L. L., Jensen, B., Nguyen, Q. T., Nøjgaard, J. K., Glaσsius, M., and Skov, H.: Atmospheric black carbon and sulfate concentrations in Northeast Greenland, Atmos. Chem. Phys., 15, 9681–9692, https://doi.org/10.5194/acp-15-9681-2015, 2015.

Monks, S. A., Arnold, S. R., Emmons, L. K., Law, K. S., Turquety, S., Duncan, B. N., Flemming, J., Huijnen, V., Tilmes, S., Langner, J., Mao, J., Long, Y., Thomas, J. L., Steenrod, S. D.,
ticle dispersion model output, Atmos. Chem. Phys., 10, 669–693,
https://doi.org/10.5194/acp-10-669-2010, 2010.

Höglund-Isaksson, L., Gómez-Sanabria, A., Klimont, Z., Rafaj, P., and Schöpp, W.: Technical potentials and costs for reducing global anthropogenic methane emissions in the 2050 timeframe – results from the GAINS model, Environmental Research Communications, 2, 025004, https://doi.org/10.1088/2515-7620/ab7457, 2020.

IPCC: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Conners, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 2391 pp., https://doi.org/10.1017/9781009157896, 2021.

IRWG: Infrared Working Group Retrieval Code, SFIT – SFIT Spectral Data Analysis Model, https://wiki.ucar.edu/display/sfit4/ (last access: 2 February 2023), 2020.

Jonsson, A. I., de Grandpré, J., Fomichev, V. I., McConnell, J. C., and Beagley, S. R.: Doubled CO₂-induced cooling in the middle atmosphere: Photochemical analysis of the ozone radiative feedback, J. Geophys. Res., 109, D24103, https://doi.org/10.1029/2004JD005093, 2004.

Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., Saha, S., White, G., Woollen, J. and Zhu, Y., The NCEP/NCAR 40-year reanalysis project, B. Am. Meteorol. Soc., 77, 437–471, https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2, 1996.

Kärnä, T. and Baptista, A. M.: Evaluation of a long-term hindcast simulation for the Columbia River estuary, Ocean Model., 99, 1–14, https://doi.org/10.1016/j.ocemod.2015.12.007, 2016.

Kawai, H., Yukimoto, S., Koshiro, T., Oshima, N., Tanaka, T., Yoshimura, H., and Nagasawa, R.: Significant improvement of cloud representation in the global climate model MRI-ESM2, Geosci. Model Dev., 12, 2875–2897, https://doi.org/10.5194/gmd-12-2875-2019, 2019.

Klimont, Z., Kupiainen, K., Heyes, C., Purohit, P., Cofala, J., Rafaj, P., Borken-Kleefeld, J., and Schöpp, W.: Global anthropogenic emissions of particulate matter including black carbon, Atmos. Chem. Phys., 17, 8681–8723, https://doi.org/10.5194/acp-17-8681-2017, 2017.

Kokkola, H., Korhonen, H., Lehtinen, K. E. J., Makkonen, R., Asmi, A., Järvenoja, S., Anttila, T., Partanen, A.-I., Kulmala, M., Järvinen, H., Laaksonen, A., and Kerminen, V.-M.: SALSA – a Sectional Aerosol module for Large Scale Applications, Atmos. Chem. Phys., 8, 2469–2483, https://doi.org/10.5194/acp-8-2469-2008, 2008.

Kuhlbrodt, T., Jones, C. G., Sellar, A., Storkey, D., Blockley, E., Stringer, M., Hill, R., Graham, T., Ridley, J., Blaker, A., Calvert, D., Copsey, D., Ellis, R., Hewitt, H., Hyder, P., Ineson, S., Mulcahy, J., Siahaan, A., and Walton, J.: The Low-Resolution Version of HadGEM3 GC3.1: Development and Evaluation for Global Climate, J. Adv. Model. Earth Syst., 10, 2865–2888, https://doi.org/10.1029/2018MS001370, 2018.

Liu, X., Ma, P.-L., Wang, H., Tilmes, S., Singh, B., Easter, R. C., Ghan, S. J., and Rasch, P. J.: Description and evaluation of a

new four-mode version of the Modal Aerosol Module (MAM4) within version 5.3 of the Community Atmosphere Model, Geosci. Model Dev., 9, 505–522, https://doi.org/10.5194/gmd-9-505-2016, 2016.

Lutsch, E., Strong, K., Jones, D. B. A., Blumenstock, T., Conway, S., Fisher, J. A., Hannigan, J. W., Hase, F., Kasai, Y., Mahieu, E., Makarova, M., Morino, I., Nagahama, T., Notholt, J., Ortega, I., Palm, M., Poberovskii, A. V., Sussmann, R., and Warneke, T.: Detection and attribution of wildfire pollution in the Arctic and northern midlatitudes using a network of Fourier-transform infrared spectrometers and GEOS-Chem, Atmos. Chem. Phys., 20, 12813–12851, https://doi.org/10.5194/acp-20-12813-2020, 2020.

Mahieu, E., Fischer, E. V., Franco, B., Palm, M., Wizenberg, T., Smale, D., Clarisse, L., Clerbaux, C., Coheur, P-F., Hannigan, J. W., Lutsch, E., Notholt, J., Pardo Cantos, I., Prignon, M., Servais, C., and Strong, K.: First retrievals of peroxy-acetyl nitrate (PAN) from ground-based FTIR solar spectra recorded at remote sites, comparison with model and satellite data, Elementa: Science of the Anthropocene, 9, 00027, https://doi.org/10.1525/elementa.2021.00027, 2021.

Makar, P. A., Gong, W., Hogrefe, C., Zhang, Y., Curci, G., Zabkar, R., Milbrandt, J., Im, U., Balzarini, A., Baró, R., Bianconi, R., Cheung, P., Forkel, R., Gravel, S., Hirtl, M., Honzak, L., Hou, A., Jiménez-Guerrero, P., Langer, M., Moran, M., Pabla, B., Pérez, J., Pirovano, G., José, R. S., Tuccella, P., Werhahn, J., Zhang, J., and Galmarini, S.: Feedbacks between air pollution and weather, part 2: Effects on chemistry, Atmos. Environ., 115, 499–526, https://doi.org/10.1016/j.atmosenv.2014.10.021, 2015a.

Makar, P. A., Gong, W., Milbrandt, J., Hogrefe, C., Zhang, Y., Curci, G., Žabkar, R., Im, U., Balzarini, A., Baró, R., Bianconi, R., Cheung, P., Forkel, R., Gravel, S., Hirtl, M., Honzak, L., Hou, A., Jiménez-Guerrero, P., Langer, M., Moran, M., Pabla, B., Pérez, J., Pirovano, G., José, R. S., Tuccella, P., Werhahn, J., Zhang, J., and Galmarini, S.: Feedbacks between air pollution and weather, part 1: Effects on weather, Atmos. Environ., 115, 442–469, https://doi.org/10.1016/j.atmosenv.2014.12.003, 2015b.

Marelle, L., Raut, J.-C., Law, K. S., Berg, L. K., Fast, J. D., Easter, R. C., Shrivastava, M., and Thomas, J. L.: Improvements to the WRF-Chem 3.5.1 model for quasi-hemispheric simulations of aerosols and ozone in the Arctic, Geosci. Model Dev., 10, 3661–3677, https://doi.org/10.5194/gmd-10-3661-2017, 2017.

Marelle, L., Raut, J.-C., Law, K. S., and Duclaux, O.: Current and Future Arctic Aerosols and Ozone From Remote Emissions and Emerging Local Sources – Modeled Source Contributions and Radiative Effects, J. Geophys. Res.-Atmos., 123, 12942–12963, https://doi.org/10.1029/2018JD028863, 2018.

Marsh, D. R., Mills, M. J., Kinnison, D. E., Lamarque, J.-F., Calvo, N., and Polvani, L. M.: Climate Change from 1850 to 2005 Simulated in CESM1 (WACCM), J. Climate, 26, 7372–7391, https://doi.org/10.1175/JCLI-D-12-00558.1, 2013.

Massling, A., Nielsen, I. E., Kristensen, D., Christensen, J. H., Sørensen, L. L., Jensen, B., Nguyen, Q. T., Nøjgaard, J. K., Glaσsius, M., and Skov, H.: Atmospheric black carbon and sulfate concentrations in Northeast Greenland, Atmos. Chem. Phys., 15, 9681–9692, https://doi.org/10.5194/acp-15-9681-2015, 2015.

Monks, S. A., Arnold, S. R., Emmons, L. K., Law, K. S., Turquety, S., Duncan, B. N., Flemming, J., Huijnen, V., Tilmes, S., Langner, J., Mao, J., Long, Y., Thomas, J. L., Steenrod, S. D.,