The Korea-United States Air Quality (KORUS-AQ) field study was conducted during May–June 2016 to understand the factors controlling air quality in South Korea. Extensive aircraft and ground network observations from the campaign offer an opportunity to address issues in current air quality models and reduce model-observation disagreements. This study examines these issues using model evaluation against the KORUS-AQ observations and intercomparisons between models. Six regional and two global chemistry transport models using identical anthropogenic emissions participated in the model intercomparison study and were used to conduct air quality simulations focusing on ozone (O3), aerosols, and their precursors for the campaign. Using the KORUSv5 emissions inventory, which has been updated from KORUSv1, the models successfully reproduced observed nitrogen oxides (NOx) and volatile organic compounds mixing ratios in surface air, especially in the Seoul Metropolitan Area, but showed systematic low biases for carbon monoxide (CO), implying possible missing CO sources in the inventory in East Asia. Although the DC-8 aircraft-observed O3 precursor mixing ratios were well captured by the models, simulated O3 levels were lower than the observations in the free troposphere in part due to too low stratospheric O3 influxes, especially in regional models. During the campaign, the synoptic meteorology played an important role in determining the observed variability of PM2.5 (PM diameter ≤ 2.5 μm) concentrations in South Korea. The models successfully simulated the observed PM2.5 variability with significant inorganic sulfate-nitrate-ammonium aerosols contribution, but failed to reproduce that of organic aerosols, causing a large inter-model variability. From the model evaluation, we find that an ensemble of model results, incorporating individual models with differing strengths and weaknesses, performs better than most individual models at representing observed atmospheric compositions for the campaign. Ongoing model development and evaluation, in close collaboration with emissions inventory development, are needed to improve air quality forecasting.

## 1. Introduction

An international air quality field study, Korea-United States Air Quality (KORUS-AQ), which was jointly hosted by the Korean National Institute of Environmental Research (NIER) and U.S. National Aeronautics and Space Administration (NASA), occurred in South Korea during May–June 2016, to understand the factors controlling air quality across urban, rural, and coastal interfaces (Crawford et al., n.d.). During the campaign, extensive surface and aircraft observations of gas and aerosol species were conducted using various instruments (Crawford et al., n.d.). In addition, a number of 3-D chemistry transport models (CTMs) were used on a daily basis to produce up to 5-day air quality forecasts for planning aircraft observations in the peninsula and nearby oceans during the campaign.

Although the forecasts were valuable for identifying pollution plumes and other features targeted for observational sampling, the forecasts sometimes differed between models and failed to capture the observed magnitudes of aerosols, ozone (O3) and their precursors. A number of explanations are possible for the differences between models as well as between models and observations, including differences in the emissions used (or from reality), meteorology, representation of chemistry and aerosol formation, and other processes affecting atmospheric composition. In this study, we will address these issues using model evaluation against extensive surface and aircraft observations from the campaign as well as intercomparisons between models.

Model intercomparison studies of regional air quality models have been extensively performed in East Asia since Carmichael et al. (2001) initiated the Model Intercomparison Study of long-range transport and sulfur deposition in East Asia (MICS-Asia). The first intercomparison study mainly focused on regional model capabilities in simulating source–receptor (S–R) relationships for sulfur deposition in East Asia. MICS-Asia II was further expanded by including nitrogen compounds, O3, and aerosols as key species for regional acid deposition with the use of observations from the Acid Deposition Monitoring Network in East Asia (EANET) to evaluate participating models (Carmichael and Ueda, 2008). Recently, the MICS-Asia III updated the intercomparison of regional simulations for O3, aerosols, and their precursors in East Asia (Chen et al., 2019; Li et al., 2019) with mechanistic analyses to explain the diversity of model performances in reproducing observed species concentrations as well as deposition fluxes in East Asia (Akimoto et al., 2019; Itahashi et al., 2020; Tan et al., 2020; Tao et al., 2020).

The MICS-Asia initiative has provided a number of important scientific findings to enhance our understanding on regional and local air pollution problems in East Asia including S–R relationships of key species such as O3, aerosols, and their precursors. Scientific findings and quantitative analyses of models’ capabilities in reproducing observations from the MICS-Asia initiative have been critical to improving our ability to predict air quality by contributing to advances in models (Carmichael and Ueda, 2008). However, these previous studies focused on the evaluation of models using observations mostly from a surface network (EANET). Lack of observations in the free troposphere could be a missing component in model evaluation and intercomparison. Our intercomparison study focuses on model evaluation using extensive aircraft observations of major air pollutants and their precursors to address the formation of inorganic aerosols by gas/particle phase partitioning, chemical formation of organic aerosols, and several other issues raised from the campaign.

During the KORUS-AQ campaign, observations of major air pollutants in South Korea, including tropospheric O3, fine particulate matter (PM), and their precursors, showed several high episodes exceeding their air quality standards (NIER and NASA, 2017). For example, maximum 8-h average O3 mixing ratios in surface air frequently exceeded the Korean 8-h air quality standard of 60 ppbv during the campaign (Peterson et al., 2019). High O3 levels were observed in airborne and O3 sonde observations, showing campaign average values of 75.8 ppbv at 700 hPa (Miyazaki et al., 2019). These high O3 levels were accompanied with high precursor mixing ratios, especially aromatic volatile organic compounds (VOCs), which accounted for a large portion of total OH reactivity in the Seoul Metropolitan Area (SMA; Simpson et al., 2020).

PM2.5 (PM diameter ≤ 2.5 μm) concentrations in surface air are typically high in South Korea in winter (Jeong and Park, 2017). Therefore, high PM2.5 episodes and their driving factors were not the prime interest of the campaign. However, daily values at surface sites reached 90 μg/m3 during May 26–31, when frontal passages were associated with foreign inflow of pollutants to the Korean peninsula (Peterson et al., 2019). The exceedance of the 24-h national standard of 35 μg/m3 occurred not only during the transport period but also during several other episodes under stagnant conditions during the campaign.

Observed aerosol chemical composition averaged at six ground sites (Bangnyung, Bulkwang, Olympic park, Gwangju, Ulsan, and Jeju) showed that inorganic salts and carbonaceous aerosols contributed to PM2.5 by 53% and 47%, respectively. Similar chemical compositions of PM were also shown in airborne observations over the SMA within the boundary layer in which 55% and 43% of PM1 (PM diameter ≤ 1 μm) were composed of inorganic and carbonaceous aerosols, respectively (Nault et al., 2018). Organic aerosols were the main component of carbonaceous aerosols with a minor contribution of black carbon aerosols (<8%; Jordan et al., 2020).

Besides the abovementioned findings, there are many additional post-mission analyses that help improve our understanding of the various factors contributing to air quality in South Korea. With the basic understanding of the campaign provided by observations, we here use an updated KORUS emissions inventory from Woo et al. (n.d.) as an input to several air quality models to simulate gas and aerosol species and further evaluate the emissions, chemistry, and physical processes that affect model performance. The model evaluation and multi-model intercomparison gives us an opportunity to better understand the differences among models, which can be used as a stepping stone to improve our scientific understanding on the issues associated with air quality simulations in the Korean peninsula and broadly in East Asia.

## 2. Method

Six regional and two global CTMs participated in this MICS (Table 1) and were used to conduct air quality simulations focusing on O3, aerosols, and their precursor species for April–June, 2016. Out of the eight simulations, five correspond to unique modeling systems, with four simulations corresponding to different configurations of the WRF-Chem model. All simulations used identical anthropogenic emissions from the KORUSv5 inventory for East Asia, which was developed by Woo et al. (n.d.). For the model evaluation against surface and aircraft observations, we mapped model grid boxes to the corresponding 60-s averaged DC-8 flight tracks (KORUS-AQ data, 2019).

Table 1.

Summary of input options of participating models in the model intercomparison. DOI: https://doi.org/10.1525/elementa.2021.00139.t1

ModelInstitutionMeteorologyBoundary ConditionsFire EmissionsBiogenic EmissionsReferences
GEOS-Chem v12.7.2 Seoul Nat’l Univ. (SNU) GEOS-FP GEOS-Chem global run GFED4 MEGAN v2.1
CAM-Chem (CESM2) NCARa MERRA2 N/A FINN v1.5 MEGAN v2.1 Emmons et al. (2020)
CAMx v6.2 Ajou Univ. FNL 1° CMAQ Northeast Asia run N/A MEGAN v2.04
WRF-Chem v4.0 NCAR FNL 0.25° CAM-Chem global run FINN v1.5 MEGAN v2.04
WRF-Chem v4.0 Pusan Nat’l Univ. (PNU) FNL 0.25° Fixed climatology from NALROM2b FINN v1.5 MEGAN v2.04
WRF-Chem v3.6.1 Univ. of Iowa FNL 1° MOZARTc global run QFED v2.4 MEGAN v2.04 Saide et al. (2020)
WRF-Chem v3.6.1 Univ. of California, Los Angeles (UCLA) FNL 1° MOZART global run QFED v2.4 MEGAN v2.04 Saide et al. (2020)
CMAQ v5.0.1 Inha Univ. FNL 0.25° MACCd reanalysis GFED4 MEGAN v2.04
ModelInstitutionMeteorologyBoundary ConditionsFire EmissionsBiogenic EmissionsReferences
GEOS-Chem v12.7.2 Seoul Nat’l Univ. (SNU) GEOS-FP GEOS-Chem global run GFED4 MEGAN v2.1
CAM-Chem (CESM2) NCARa MERRA2 N/A FINN v1.5 MEGAN v2.1 Emmons et al. (2020)
CAMx v6.2 Ajou Univ. FNL 1° CMAQ Northeast Asia run N/A MEGAN v2.04
WRF-Chem v4.0 NCAR FNL 0.25° CAM-Chem global run FINN v1.5 MEGAN v2.04
WRF-Chem v4.0 Pusan Nat’l Univ. (PNU) FNL 0.25° Fixed climatology from NALROM2b FINN v1.5 MEGAN v2.04
WRF-Chem v3.6.1 Univ. of Iowa FNL 1° MOZARTc global run QFED v2.4 MEGAN v2.04 Saide et al. (2020)
WRF-Chem v3.6.1 Univ. of California, Los Angeles (UCLA) FNL 1° MOZART global run QFED v2.4 MEGAN v2.04 Saide et al. (2020)
CMAQ v5.0.1 Inha Univ. FNL 0.25° MACCd reanalysis GFED4 MEGAN v2.04

MEGAN = Model of Emissions of Gases and Aerosols from Nature.

aNational Center for Atmospheric Research.

bNOAA-Aeronomy Laboratory Regional Oxidation Model.

cModel for Ozone and Related chemical Tracers.

dMonitoring Atmospheric Composition and Climate.

Due to the different simulation domains from each model (Figure 1), available data coverage differed. GEOS-Chem used a nested domain with 0.25° × 0.3125° (latitude × longitude) spatial resolution for East Asia with boundary conditions from a 2° × 2.5° global simulation. CAM-Chem was used to conduct a global simulation with 0.47° × 0.63° spatial resolution. Also note that the Iowa and UCLA WRF-Chem used the smallest domain with the finest horizontal resolution (4 km), which is a nested one within the outer domain with a 20 km horizontal resolution, covering East Asia. CMAQ also used a nested simulation (9 km) with an outer domain of 27 km resolution covering East Asia. Hourly mean simulated gas and aerosol concentrations were used for the multi-model intercomparisons and comparisons with observations at surface sites.

Figure 1.

Simulation domain of each participating model for the model intercomparison. The domain for GEOS-Chem indicates the nested simulation domain. CAM-Chem was simulated in a global domain and therefore is not shown. Iowa and UCLA WRF-Chem used an outer domain with 20 km resolution covering 16.8–58.2°N and 68.2–148.7°E. CMAQ used an outer domain with 27 km resolution covering 19.9–49.2°N and 104.6–160.3°E. DOI: https://doi.org/10.1525/elementa.2021.00139.f1

Figure 1.

Simulation domain of each participating model for the model intercomparison. The domain for GEOS-Chem indicates the nested simulation domain. CAM-Chem was simulated in a global domain and therefore is not shown. Iowa and UCLA WRF-Chem used an outer domain with 20 km resolution covering 16.8–58.2°N and 68.2–148.7°E. CMAQ used an outer domain with 27 km resolution covering 19.9–49.2°N and 104.6–160.3°E. DOI: https://doi.org/10.1525/elementa.2021.00139.f1

### 2.1. Model description

Tables 1 and 2 summarize the participating models with brief descriptions of their driving meteorology, boundary conditions, wildfire emissions inventories, biogenic emissions schemes, and grid resolutions of the simulations. We allowed all models to choose their own options for meteorological fields and natural emissions. Therefore, despite the use of the same anthropogenic emissions, a difference in total emissions exists between models, which is discussed in Section 2.2.

Table 2.

Summary of model resolutions of participating models in the model intercomparison. DOI: https://doi.org/10.1525/elementa.2021.00139.t2

ModelHorizontal Resolution (Latitude × Longitude)Vertical Levels (<1.5 km)Model Top Height (km)Surface Level Height (m)
GEOS-Chem (SNU) 0.25° × 0.3125° 47 (11) 80 123
CAM-Chem (NCAR) 0.47° × 0.63° 56 (13) 45 120
CAMx (Ajou) 27 × 27 km 22 (12) 19 32
WRF-Chem (NCAR) 15 × 15 km 37 (13) 30 35
WRF-Chem (PNU) 27 × 27 km 31 (10) 19 38
WRF-Chem (Iowa) 4 × 4 km 53 (14) 20 47.5
WRF-Chem (UCLA) 4 × 4 km 53 (14) 20 47.5
CMAQ (Inha) 9 × 9 km 23 (10) 12 35
ModelHorizontal Resolution (Latitude × Longitude)Vertical Levels (<1.5 km)Model Top Height (km)Surface Level Height (m)
GEOS-Chem (SNU) 0.25° × 0.3125° 47 (11) 80 123
CAM-Chem (NCAR) 0.47° × 0.63° 56 (13) 45 120
CAMx (Ajou) 27 × 27 km 22 (12) 19 32
WRF-Chem (NCAR) 15 × 15 km 37 (13) 30 35
WRF-Chem (PNU) 27 × 27 km 31 (10) 19 38
WRF-Chem (Iowa) 4 × 4 km 53 (14) 20 47.5
WRF-Chem (UCLA) 4 × 4 km 53 (14) 20 47.5
CMAQ (Inha) 9 × 9 km 23 (10) 12 35

The third column indicates the total number of vertical levels and the number of levels below 1.5 km inside parentheses.

Table 3 summarizes gas-phase chemistry mechanisms, aerosol thermodynamics/microphysics modules, and secondary organic aerosol (SOA) schemes used in the models. As shown in the table, substantial differences in terms of the number of chemical species and complexity of chemical schemes exist between models. For example, the GEOS-Chem Tropchem mechanism has the largest number of kinetic reactions (713) and chemical species (220), while other schemes typically have 100–300 reactions, and some hydrocarbon mechanisms, such as the one used in Iowa WRF-Chem, were even more simplified for computational efficiency. Although complex gas phase chemistry is embedded within GEOS-Chem, the Tropchem mechanism, along with SAPRC99 used in CAMx, both lack detailed aromatic chemistry, compared to other hydrocarbon mechanisms such as the MOZART and RACM-based schemes, used in CAM-Chem and WRF-Chem (Carter, 2000; Ahmadov et al., 2012; Knote et al., 2014; Emmons et al., 2020).

Table 3.

Summary of chemistry and aerosol options used in the model intercomparison. DOI: https://doi.org/10.1525/elementa.2021.00139.t3

ModelChemistry MechanismAerosol Thermodynamics + MicrophysicsSOA Precursors + Scheme
GEOS-Chem (SNU) GEOS-Chem Tropchem (Mao et al., 2012; Fisher et al., 2016; Marais et al., 2016ISORROPIA II (Fountoukis and Nenes, 2007) + Bulk scheme Anthropogenic and biomass burning precursors scaled by CO emission + biogenic precursors scaled by isoprene, monoterpene emissions + simplified SOA with no VBS (Hodzic and Jimenez, 2011
CAM-Chem (NCAR) MOZART-T1 (Emmons et al., 2020N/A + MAM4 (Liu et al., 2016Aromatics, isoprene, terpenes, S/IVOCs + MAM4 with 5-bin VBS (Tilmes et al., 2019
CAMx (Ajou) SAPRC99 (Carter, 2000ISORROPIA I (Nenes et al., 1998) + CF scheme Aromatics, isoprene, terpenes + the original Secondary Organic Aerosol Partitioning (SOAP) (Strader et al., 1999
WRF-Chem (NCAR) MOZART-4 + updates from Knote et al. (2014)  MESA + MOSAIC (Zaveri et al., 2008Anthropogenic and biomass burning precursors scaled by CO emission + simplified SOA with no VBS (Hodzic and Jimenez, 2011) + 2-product approach for biogenic SOA (Shrivastava et al., 2011
WRF-Chem (PNU) RACM-ESRL (Ahmadov et al., 2015MARS (Grell et al., 2005) + MADE (Ackermann et al., 1998Alkanes, alkenes, aromatics, cresol, isoprene, terpenes + 4-bin VBS from Ahmadov et al. (2012)
WRF-Chem (Iowa) Reduced hydrocarbon mechanism from Pfister et al. (2014)  MESA + MOSAIC Anthropogenic and biomass burning precursors scaled by CO emission + simplified SOA with no VBS (Hodzic and Jimenez, 2011) + 2-product approach for biogenic SOA (Shrivastava et al., 2011
WRF-Chem (UCLA) RACM-ESRL MARS + MADE Alkanes, alkenes, aromatics, cresol, isoprene, terpenes + 4-bin VBS from Ahmadov et al. (2012)
CMAQ (Inha) CB05-TU (Whitten et al., 2010ISORROPIA I + AERO5 (Foley et al., 2010Alkanes, aromatics, isoprene, terpenes + 2-product approach using yields from Carlton et al. (2010)
ModelChemistry MechanismAerosol Thermodynamics + MicrophysicsSOA Precursors + Scheme
GEOS-Chem (SNU) GEOS-Chem Tropchem (Mao et al., 2012; Fisher et al., 2016; Marais et al., 2016ISORROPIA II (Fountoukis and Nenes, 2007) + Bulk scheme Anthropogenic and biomass burning precursors scaled by CO emission + biogenic precursors scaled by isoprene, monoterpene emissions + simplified SOA with no VBS (Hodzic and Jimenez, 2011
CAM-Chem (NCAR) MOZART-T1 (Emmons et al., 2020N/A + MAM4 (Liu et al., 2016Aromatics, isoprene, terpenes, S/IVOCs + MAM4 with 5-bin VBS (Tilmes et al., 2019
CAMx (Ajou) SAPRC99 (Carter, 2000ISORROPIA I (Nenes et al., 1998) + CF scheme Aromatics, isoprene, terpenes + the original Secondary Organic Aerosol Partitioning (SOAP) (Strader et al., 1999
WRF-Chem (NCAR) MOZART-4 + updates from Knote et al. (2014)  MESA + MOSAIC (Zaveri et al., 2008Anthropogenic and biomass burning precursors scaled by CO emission + simplified SOA with no VBS (Hodzic and Jimenez, 2011) + 2-product approach for biogenic SOA (Shrivastava et al., 2011
WRF-Chem (PNU) RACM-ESRL (Ahmadov et al., 2015MARS (Grell et al., 2005) + MADE (Ackermann et al., 1998Alkanes, alkenes, aromatics, cresol, isoprene, terpenes + 4-bin VBS from Ahmadov et al. (2012)
WRF-Chem (Iowa) Reduced hydrocarbon mechanism from Pfister et al. (2014)  MESA + MOSAIC Anthropogenic and biomass burning precursors scaled by CO emission + simplified SOA with no VBS (Hodzic and Jimenez, 2011) + 2-product approach for biogenic SOA (Shrivastava et al., 2011
WRF-Chem (UCLA) RACM-ESRL MARS + MADE Alkanes, alkenes, aromatics, cresol, isoprene, terpenes + 4-bin VBS from Ahmadov et al. (2012)
CMAQ (Inha) CB05-TU (Whitten et al., 2010ISORROPIA I + AERO5 (Foley et al., 2010Alkanes, aromatics, isoprene, terpenes + 2-product approach using yields from Carlton et al. (2010)

MOSAIC = Model for Simulating Aerosol Interactions and Chemistry; MADE = Modal Aerosol Dynamics Model for Europe; CF = Coarse/Fine; MAM4 = Modal Aerosol Model with 4 modes; VBS = volatility basis set.

As for the inorganic aerosol thermodynamics calculation, three different modules were used including ISORROPIA, Multicomponent Equilibrium Solver for Aerosols (MESA), and Model for an Aerosol Reacting System (MARS). ISOROPPIA computes gas-particle partitioning of nitric acid and ammonia (Fountoukis and Nenes, 2007) and MESA solves the solid-liquid equilibria within each sectional aerosol size bin (Zaveri et al., 2008). MARS calculates the thermodynamic equilibrium of sulfate, nitrate, and ammonium (Grell et al., 2005).

For aerosol microphysics simulations, most models except for GEOS-Chem used either a sectional or modal approach. WRF-Chem used either the Model for Simulating Aerosol Interactions and Chemistry (Zaveri et al., 2008) or Modal Aerosol Dynamics Model for Europe (Ackermann et al., 1998). CAMx employed the Coarse/Fine scheme with two static modes. CAM-Chem and CMAQ used the Modal Aerosol Model with 4 modes (Liu et al., 2016) and AERO5 (Foley et al., 2010), respectively, and GEOS-Chem used a bulk aerosol scheme.

Three types of SOA schemes were used in the models: a simplified scheme developed by Hodzic and Jimenez (2011), the 2-product approach (Odum et al., 1996), and the volatility basis set (VBS) approach (Donahue et al., 2006; Stanier et al., 2008). In the atmosphere, SOA can be formed by the oxidation of parent hydrocarbons, creating oxygenated semi-volatile compounds that either partition to the particle phase or undergo continuous oxidation. However, due to the complexity of various parent hydrocarbons, the simplified SOA scheme assumes that a lumped SOA precursor is emitted from the same sources as CO, and it is converted irreversibly to SOA with a fixed lifetime (approximately 1 day; Cubison et al., 2011; Hodzic and Jimenez, 2011; Hayes et al., 2015; Kim et al., 2015; Shrivastava et al., 2017). The 2-product approach, which was developed by Odum et al. (1996), requires speciated parent hydrocarbons that undergo oxidation and produces two SOA surrogates that partition (gas to particle) using yields from chamber experiments. Finally, the VBS approach (Donahue et al., 2006; Stanier et al., 2008) is the most complex scheme, which divides the semi-volatile oxidation products into several bins according to their volatilities to represent the continuous oxidation process of SOA.

In this intercomparison, we do not intend to judge individual schemes used in models or whether one scheme performs better than others. Instead, we try to understand the issues present with current schemes and configurations used in air quality simulations and gain some insights to improve them based on the findings.

### 2.2. Emissions

All models used an identical anthropogenic emissions inventory (KORUSv5), which was developed by Konkuk University for the campaign. Details on the inventory can be found elsewhere (Woo et al., n.d.). The inventory includes area, point, mobile, and ship emissions of species such as carbon monoxide (CO), nitrogen oxides (NOx = NO + NO2), ammonia (NH3), sulfur dioxide (SO2), and VOCs. Annual total emissions of individual species from the KORUS version 1 (v1) and KORUS version 5 (v5) inventories are summarized for East China and South Korea in Table 4. There are significant changes in the emissions of a few species from v1 to v5 inventories. For example, the anthropogenic CO and NOx emissions in South Korea were increased by factors of 2.5 and 1.4, respectively, but they were decreased by about 25% in East China. A more dramatic decrease up to 70% is shown for SO2 in East China. On the other hand, anthropogenic VOCs emissions were generally increased in both East China and South Korea from v1 to v5, especially for aromatic species, which were shown to be higher by factors of 1.5–2.3 in the v5 inventory compared to those of v1. One thing to note is that participating models used their own emission processors for chemical speciation and imposing diurnal variation of species emissions, which could result in differences of air quality simulations between the models despite the use of the same emissions inventory (Choi et al., 2019; Goldberg et al., 2019).

Table 4.

Total anthropogenic emissions (Tg/yr) in the KORUS version 1 and version 5 inventories in East China (27.7–40°N, 115–123°E) and South Korea (34.5–38°N, 126–130°E). DOI: https://doi.org/10.1525/elementa.2021.00139.t4

SpeciesKORUSv1 (Tg/yr)KORUSv5 (Tg/yr)
East ChinaSouth KoreaEast ChinaSouth Korea
CO 58.0 0.77 43.9 1.9
NOx 10.8 0.93 7.83 1.31
NH3 4.1 0.29 2.54 0.29
SO2 9.6 0.33 3.19 0.336
C2H6 0.25 0.02 0.568 0.029
C3H8 0.45 0.04 0.427 0.033
Alkanes 3.03 0.273 3.71 0.283
Alkenes 1.36 0.113 1.60 0.128
Aromatics 2.29 0.203 3.56 0.467
HCHO 0.10 0.008 0.04 0.009
pSO4 0.25 0.006 0.22 0.002
POA 1.02 0.02 0.43 0.05
BC 0.34 0.01 0.22 0.018
SpeciesKORUSv1 (Tg/yr)KORUSv5 (Tg/yr)
East ChinaSouth KoreaEast ChinaSouth Korea
CO 58.0 0.77 43.9 1.9
NOx 10.8 0.93 7.83 1.31
NH3 4.1 0.29 2.54 0.29
SO2 9.6 0.33 3.19 0.336
C2H6 0.25 0.02 0.568 0.029
C3H8 0.45 0.04 0.427 0.033
Alkanes 3.03 0.273 3.71 0.283
Alkenes 1.36 0.113 1.60 0.128
Aromatics 2.29 0.203 3.56 0.467
HCHO 0.10 0.008 0.04 0.009
pSO4 0.25 0.006 0.22 0.002
POA 1.02 0.02 0.43 0.05
BC 0.34 0.01 0.22 0.018

Individual emissions for NO and NO2 are provided in the inventories and the sum is used to represent NOx.

After the campaign, various modeling analyses using 3-D CTMs and 0-D box models, combined with remote sensing and in situ measurements, were applied to constrain the anthropogenic emissions in South Korea. Miyazaki et al. (2019) and Goldberg et al. (2019) used a top-down method accompanying satellite retrievals and CTMs, and Oak et al. (2019) conducted model evaluation against airborne data to examine bottom-up emission estimates for South Korea. These studies suggested 40–50% and 83% increases in NOx and CO emissions in South Korea, respectively. Updates in KORUSv5 emissions were in part based on these studies. All the rationales and supporting information for the updates of the KORUS inventory are available in Woo et al. (n.d.).

Participating models used either Model of Emissions of Gases and Aerosols from Nature (MEGAN) versions 2.04 or 2.1 (Guenther et al., 2012) for biogenic emissions of isoprene, terpenes, and other VOCs. However, each model used its own vegetation map and meteorology to determine the emissions, so there were differences in the resulting isoprene emissions among the models.

Different biomass burning emissions inventories were also used in the participating models. Wildfire emissions were generally much lower than anthropogenic emissions over South Korea during the campaign (Tang et al., 2019), but there existed a short period (May 18) when the influence of Siberian fires reached the Korean peninsula (Lamb et al., 2018; Peterson et al., 2019). Among the three different biomass burning inventories (GFED4, FINNv1.5, QFEDv2.4) used in the models, GFED4 had the largest CO emissions in South Korea (0.187, 0.039, and 0.068 Tg/yr, respectively). CO emissions from each inventory in eastern China and eastern Russia ranged from 0.8 to 1.6 Tg/yr and 33 to 74 Tg/yr, with GFED4 also being the largest and FINNv1.5 being the smallest in both regions. Emissions of primary organic aerosols (POA) and black carbon (BC) were also smallest in FINNv1.5, by factors of 2–10 for POA and 3–12 for BC, compared to those of GFED4 and QFEDv2.4 in the eastern Eurasia domain. POA and BC emissions from QFEDv2.4 were especially high compared to other inventories in eastern Russia. However, this difference had a minimal influence on the inter-model variability of PM simulations in South Korea among the participating models because of limited transboundary transport from Siberia into the peninsula.

### 2.3. Observations

We used observations from surface sites and the DC-8 aircraft during the campaign to evaluate the air quality simulations. Figure 2 shows site locations of the AirKorea network managed by NIER, affiliated with the Korean Ministry of Environment, which conducts regular measurements of six regulatory air pollutants including CO, O3, NO2, SO2, PM2.5, and PM10 (PM diameter ≤ 10 μm). This monitoring network provides hourly volume mixing ratios of gaseous pollutants and PM mass concentrations across the country, particularly focusing on urban air quality to which a large population is exposed. Therefore, most sites are situated in or around the major cities. However, because the AirKorea observations do not provide information on the chemical composition of PM, we used submicron particulate species (PM1) measurements conducted by the Korea Institute of Science and Technology (KIST), located in Seoul. Chemical compositions of non-refractive PM1 including sulfate (SO42–), nitrate (NO3), ammonium (NH4+), and organic aerosol (OA) were observed and analyzed using a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS; DeCarlo et al., 2006) and the positive matrix factorization (PMF) method by Kim et al. (2018). BC mass concentrations were also measured using a multi-angle absorption photometer (Petzold and Schönlinner, 2004) by KIST.

Figure 2.

Locations of 320 AirKorea (clustered in different colors for different regions) ground monitoring stations used in this study. Additional surface observations were conducted at KIST (37.601°N, 127.045°E), which is shown in red. Daesan (36.95°N, 126.41°E), where large power plant point sources are located, is also shown in black. DOI: https://doi.org/10.1525/elementa.2021.00139.f2

Figure 2.

Locations of 320 AirKorea (clustered in different colors for different regions) ground monitoring stations used in this study. Additional surface observations were conducted at KIST (37.601°N, 127.045°E), which is shown in red. Daesan (36.95°N, 126.41°E), where large power plant point sources are located, is also shown in black. DOI: https://doi.org/10.1525/elementa.2021.00139.f2

During the KORUS-AQ campaign, 20 research flights (RFs) of the NASA DC-8 aircraft were performed in the daytime (08–16 local standard time [LST]) and conducted extensive measurements of ambient concentrations of trace gases and aerosols over South Korea with high spatial and temporal resolutions (NIER and NASA, 2017; Crawford et al., n.d.). In order to focus on general characteristics during the campaign, we excluded observations from two flights (RF7 and RF18) that were designed for frontal cloud profiling and point source surveying. We also excluded observations near Daesan (Figure 2) from all flights to avoid influences from large power plant point sources (36.4–37.15°N, 126–126.88°E) in the comparisons below. Here we used 60-s averaged airborne observations of gaseous species including CO, O3, NOx, HNO3, H2O2, PAN, ANs (= peroxy + alkyl nitrates), VOCs, and HCHO for model evaluation. Table 5 summarizes instruments used for the measurement of each species.

Table 5.

Summary of observed gas species and instruments onboard DC-8. DOI: https://doi.org/10.1525/elementa.2021.00139.t5

SpeciesInstrumentInstitutionReferences
CO Differential Absorption CO Measurement (DACOM) NASA Langley Research Center Sachse et al. (1991)
O3, NO, NO2 NCAR NOxyO3 Chemiluminescence National Center for Atmospheric Research Ridley and Grahek (1990); Weinheimer et al. (1994)
HNO3 California Institute of Technology Chemical Ionization Mass Spectrometer (CIT-CIMS) California Institute of Technology Crounse et al. (2006); Paulot et al. (2009); St. Clair et al. (2010)
PAN Georgia Tech Chemical Ionization Mass Spectrometer (GTCIMS) Georgia Institute of Technology Slusher et al. (2004)
RO2NO2 (peroxy nitrates), RONO2 (alkyl nitrates) Thermal Dissociation-Laser Induced Fluorescence (TD-LIF) University of California, Berkeley Day et al. (2002); Wooldridge et al. (2010)
Toluene Whole Air Sampler (WAS) University of California, Irvine Colman et al. (2001)
Isoprene Proton Transfer Time of Flight Mass Spectrometer (PTR-ToF-MS) University of Oslo Müller et al. (2014)
HCHO Compact Atmospheric Multispecies Spectrometer (CAMS) University of Colorado Richter et al. (2015)
SpeciesInstrumentInstitutionReferences
CO Differential Absorption CO Measurement (DACOM) NASA Langley Research Center Sachse et al. (1991)
O3, NO, NO2 NCAR NOxyO3 Chemiluminescence National Center for Atmospheric Research Ridley and Grahek (1990); Weinheimer et al. (1994)
HNO3 California Institute of Technology Chemical Ionization Mass Spectrometer (CIT-CIMS) California Institute of Technology Crounse et al. (2006); Paulot et al. (2009); St. Clair et al. (2010)
PAN Georgia Tech Chemical Ionization Mass Spectrometer (GTCIMS) Georgia Institute of Technology Slusher et al. (2004)
RO2NO2 (peroxy nitrates), RONO2 (alkyl nitrates) Thermal Dissociation-Laser Induced Fluorescence (TD-LIF) University of California, Berkeley Day et al. (2002); Wooldridge et al. (2010)
Toluene Whole Air Sampler (WAS) University of California, Irvine Colman et al. (2001)
Isoprene Proton Transfer Time of Flight Mass Spectrometer (PTR-ToF-MS) University of Oslo Müller et al. (2014)
HCHO Compact Atmospheric Multispecies Spectrometer (CAMS) University of Colorado Richter et al. (2015)

For airborne observations of non-refractive PM1, we used 60-s averaged HR-ToF-AMS (DeCarlo et al., 2006) observations of particulate phase SO42–, NO3, NH4+, and OA conducted by University of Colorado at Boulder. Detailed descriptions of AMS measurement methods and PMF analysis used to determine the aerosol composition can be found in Nault et al. (2018). BC mass concentrations within the size range of 100–500 nm were measured using the Humidified Dual-Single Particle Soot Photometer (HD-SP2) by the NOAA Chemical Sciences Laboratory (Schwarz et al., 2013).

## 3. Model evaluation

We here conduct model evaluations by using various observations from the surface network and DC-8 aircraft during the KORUS-AQ campaign. This evaluation allows us to better understand issues with models used for the air quality simulations. Relative contributions to air quality in South Korea by local and transboundary pollution influences were predominantly determined by synoptic meteorology (Peterson et al., 2019). The campaign period can be grouped into four distinct synoptic regimes: dynamic weather (May 1–16), stagnant period (May 17–22), transport period (May 24–31), and blocking meteorology (June 1–10; Peterson et al., 2019), and model evaluations are conducted based on these regimes.

Model evaluations for surface air are based on hourly mean data including all hours during the campaign period, whereas model evaluations with the aircraft observations are based on daytime hours (08–16 LST) due to limited temporal coverage of the DC-8. For comparisons of species in surface air shown in Figure 3, all the simulated results with different spatial resolutions were first regridded to 0.5° resolution, then averaged to obtain a single 2-D gridded data set to represent the ensemble model.

Figure 3.

Model ensemble and observed mean CO, O3, NO2, SO2, PM2.5, and PM10 surface concentrations during the KORUS-AQ campaign (May 1–June 10). Simulated ensemble fields of the lowest model layer are shown in the background with observations overlaid in circles, which correspond to 320 AirKorea surface monitoring stations. The normalized mean bias (NMB) and Pearson correlation coefficient (R) are calculated using observations gridded to 0.5° horizontal resolution and the ensemble model. DOI: https://doi.org/10.1525/elementa.2021.00139.f3

Figure 3.

Model ensemble and observed mean CO, O3, NO2, SO2, PM2.5, and PM10 surface concentrations during the KORUS-AQ campaign (May 1–June 10). Simulated ensemble fields of the lowest model layer are shown in the background with observations overlaid in circles, which correspond to 320 AirKorea surface monitoring stations. The normalized mean bias (NMB) and Pearson correlation coefficient (R) are calculated using observations gridded to 0.5° horizontal resolution and the ensemble model. DOI: https://doi.org/10.1525/elementa.2021.00139.f3

For comparisons with the surface and aircraft observations in Figures 4 and 711, we used individual model results with different spatial resolutions by sampling simulated values from the closest grid box to the observed location. Model results were also sampled along the DC-8 flight track at the observed hour so that simulated values are coherent in time and space with observations. Here the ensemble model is the arithmetic mean of all models with no regridding or interpolation. For model evaluation, we use three statistical metrics: normalized mean bias (NMB), Pearson correlation coefficient (R), and root mean square error (RMSE), which are defined below.

(1) $NMB=∑i=1N(Simulationi−Observationi)∑i=1NObservationi,$

where N is sample size,

(2) $R=∑i=1N(Simulationi−S¯)(Observationi−O¯)∑i=1N(Simulationi−S¯)2∑i=1N(Observationi−O¯)2,$

where the overbar () is sample mean, and

(3) $RMSE=1N∑i=1N(Simulationi−Observationi)2.$

Figure 4.

Simulated and observed daytime (8-h) O3 (top) and daily PM2.5 (bottom) concentrations averaged at AirKorea surface sites located in Seoul, Busan, Incheon, Gwangju, Yeosu, and Gangwon. Observations are in black solid lines and individual models are in different colors with the ensemble model in red solid lines. Model ranges are shown in gray shades. Statistics (R, NMB, and RMSE) of each model and the ensemble model performances are denoted in the upper-left corners. DOI: https://doi.org/10.1525/elementa.2021.00139.f4

Figure 4.

Simulated and observed daytime (8-h) O3 (top) and daily PM2.5 (bottom) concentrations averaged at AirKorea surface sites located in Seoul, Busan, Incheon, Gwangju, Yeosu, and Gangwon. Observations are in black solid lines and individual models are in different colors with the ensemble model in red solid lines. Model ranges are shown in gray shades. Statistics (R, NMB, and RMSE) of each model and the ensemble model performances are denoted in the upper-left corners. DOI: https://doi.org/10.1525/elementa.2021.00139.f4

### 3.1. Evaluation for surface air

Figure 3 shows comparisons of simulated versus observed CO, O3, NO2, SO2, PM2.5, and PM10 concentrations in surface air averaged for the campaign period. Observed concentrations are from the AirKorea network and simulated values are from the ensemble of participating models. Observations show highest values in the SMA and industrial regions except for O3, which is strongly titrated by high NOx levels from traffic emissions (Figure 3c).

The spatial distributions of observed O3 mixing ratios are generally well captured by the model ensemble (R = .55), which clearly shows negative spatial correlations between simulated NO2 and O3. The models tend to underestimate the observations in part owing to the coarse spatial resolutions to simulate concentrated plumes from urban areas where most sites are located. The simulated low biases appear to be smaller for secondary pollutants relative to the primary; this is especially evident for CO, which was known to be severely underestimated in CTMs in South Korea (Huang et al., 2018; Tang et al., 2019; Gaubert et al., 2020; Lee et al., 2020). This low bias in simulated CO (NMB = –47%) indicates that there are additional missing sources of CO in the emissions inventory even though the domestic anthropogenic CO emissions have increased by a factor of 2.5 in the KORUSv5 inventory relative to that of the KORUSv1 inventory.

The model ensemble reproduced observed SO2 mixing ratios with fair spatial correlation (R = .74), capturing high values from major power plants in the Daesan, Yeosu, and Ulsan areas, but in general showed an underestimation (NMB = –27%). We can see large discrepancies in PM10 concentrations between the model ensemble and the observations, indicating a large uncertainty in simulating primary dust aerosols in models (Huneeus et al., 2011; Jeong and Park, 2018).

Figure 4 shows the time series of daytime O3 and 24-h mean PM2.5 concentrations averaged at AirKorea surface sites located in Seoul, Busan, Incheon, Gwangju, Yeosu, and Gangwon. The four different synoptic regimes are indicated in colored shadings. Simulated results from the participating models using KORUSv1 and KORUSv5 emissions are shown, although not all participating model results with the KORUSv1 emissions were available. From this comparison, a few notable points are found as follows:

1. In general, the models have capability in reproducing observed temporal variations of both O3 (R ≥ 0.53) and PM2.5 (R ≥ 0.69) in surface air in South Korea.

2. The model ensemble always outperforms individual models for O3 and PM2.5 during the KORUS-AQ.

3. The model performance in terms of R, NMB, and RMSE has improved using the KORUSv5 emissions relative to that of the KORUSv1 emissions, indicating that the campaign was very helpful to better understand the pollutant emissions in the peninsula.

4. The ensemble model with the KORUSv1 emissions significantly underestimated both observed O3 and PM2.5 concentrations by 18% and 47%, respectively, during May 19–24, due to underestimations of O3 precursors (i.e., NOx, aromatics) and OA. The models are much closer to observations with the KORUSv5 emissions, showing smaller biases (–7% and −22%, respectively) between the ensemble model and observations, although inter-model variability is larger with the KORUSv5 emissions relative to that with the KORUSv1 emissions.

### 3.2. Evaluation for aircraft observations

The DC-8 collected observations on 20 flights during the campaign. Each flight was designed based on the air quality forecasts to examine characteristics of air pollutants either emitted domestically or transported from outside South Korea. We divided individual flight observations and model results into four different synoptic regimes as discussed earlier to separate local and foreign influences. For model evaluation of vertical profiles (Figures 79) using DC-8 observations, we used data collected within the SMA region (37–37.6°N, 126.6–127.7°E), and statistics for individual models are summarized in Table 6. A total of 55 missed approaches were repeatedly conducted over this region, which provided more representative information on the atmospheric conditions during the campaign.

Table 6.

Summary of campaign average NMB (%) between models and observation for Figures 79. NMB below 1.5 km are presented inside parentheses. DOI: https://doi.org/10.1525/elementa.2021.00139.t6

SpeciesEnsembleGEOS-ChemCAM-ChemCAMxWRF-Chem (NCAR)WRF-Chem (PNU)WRF-Chem (Iowa)WRF-Chem (UCLA)CMAQ
CO –27 (–32) –25 (–29) –44 (–47) –43 (–49) –36 (–37) –42 (–47) –7 (–16) 0 (–8) –17 (–22)
O3 –16 (–18) –26 (–28) –9 (–15) –16 (–20) –10 (–15) –29 (–32) –11 (–14) –16 (–17) –9 (–7)
NOx 12 (11) –5 (–9) 11 (15) 18 (3) 56 (60) 19 (14) 3 (6) –7 (–3) 0 (4)
Toluene –9 (–9) –36 (–37) –3 (–0) –6 (–13) –30 (–28) –7 (–9) –15 (–12) 17 (21) 14 (15)
Isoprene –10 (–6) –28 (–33) –20 (–8) –37 (–37) –22 (–12) 14 (18) –14 (–2) –16 (–10) 46 (42)
HCHO –30 (–32) –33 (–38) –34 (–33) –41 (–53) –33 (–34) –24 (–28) –46 (–44) –26 (–19) –3 (–3)
BC 36 (36) 31 (28) –17 (–18) 12 (2) 67 (67) –21 (–24) 54 (60) 30 (34) 94 (94)
OA –46 (–46) –36 (–39) –38 (–39) –85 (–87) –33 (–27) –44 (–49) –54 (–54) –30 (–28) –36 (–36)
SO42– –4 (–10) –49 (–53) –35 (–44) –2 (–16) –13 (–13) –4 (–15) 8 (4) –5 (–9) 69 (71)
NO3 –33 (–37) 9 (–18) N/A –57 (–66) –12 (–5) –45 (–61) –2 (5) –29 (–21) –22 (–19)
NH4+ –26 (–29) –19 (–34) N/A –26 (–40) –14 (–9) –26 (–38) –16 (–17) –19 (–15) 19 (23)
PM1 –31 (–34) –24 (–34) –58 (–62) –54 (–62) –21 (–17) –33 (–42) –27 (–22) –23 (–20) –3 (–3)
NOx/NOy 0 (1) 5 (3) –24 (–14) 12 (6) –16 (–5) 2 (9) –10 (–9) –22 (–13) –29 (–17)
HNO3/NOy 74 (77) 12 (63) 187 (213) 94 (99) 128 (116) 117 (111) 41 (29) 69 (81) 29 (17)
pNO3/NOy 77 (54) 36 (–12) N/A –70 (–71) –38 (–33) –37 (–61) –4 (5) –43 (–35) –28 (–28)
PAN/NOy –34 (–29) –16 (–31) –32 (–39) –30 (–13) –35 (–42) –61 (–56) –19 (–16) 26 (28) –13 (22)
ANs/NOy –46 (–44) –88 (–88) –85 (–82) –66 (–64) –88 (–87) –74 (–72) –3 (2) –57 (–51) 136 (124)
SpeciesEnsembleGEOS-ChemCAM-ChemCAMxWRF-Chem (NCAR)WRF-Chem (PNU)WRF-Chem (Iowa)WRF-Chem (UCLA)CMAQ
CO –27 (–32) –25 (–29) –44 (–47) –43 (–49) –36 (–37) –42 (–47) –7 (–16) 0 (–8) –17 (–22)
O3 –16 (–18) –26 (–28) –9 (–15) –16 (–20) –10 (–15) –29 (–32) –11 (–14) –16 (–17) –9 (–7)
NOx 12 (11) –5 (–9) 11 (15) 18 (3) 56 (60) 19 (14) 3 (6) –7 (–3) 0 (4)
Toluene –9 (–9) –36 (–37) –3 (–0) –6 (–13) –30 (–28) –7 (–9) –15 (–12) 17 (21) 14 (15)
Isoprene –10 (–6) –28 (–33) –20 (–8) –37 (–37) –22 (–12) 14 (18) –14 (–2) –16 (–10) 46 (42)
HCHO –30 (–32) –33 (–38) –34 (–33) –41 (–53) –33 (–34) –24 (–28) –46 (–44) –26 (–19) –3 (–3)
BC 36 (36) 31 (28) –17 (–18) 12 (2) 67 (67) –21 (–24) 54 (60) 30 (34) 94 (94)
OA –46 (–46) –36 (–39) –38 (–39) –85 (–87) –33 (–27) –44 (–49) –54 (–54) –30 (–28) –36 (–36)
SO42– –4 (–10) –49 (–53) –35 (–44) –2 (–16) –13 (–13) –4 (–15) 8 (4) –5 (–9) 69 (71)
NO3 –33 (–37) 9 (–18) N/A –57 (–66) –12 (–5) –45 (–61) –2 (5) –29 (–21) –22 (–19)
NH4+ –26 (–29) –19 (–34) N/A –26 (–40) –14 (–9) –26 (–38) –16 (–17) –19 (–15) 19 (23)
PM1 –31 (–34) –24 (–34) –58 (–62) –54 (–62) –21 (–17) –33 (–42) –27 (–22) –23 (–20) –3 (–3)
NOx/NOy 0 (1) 5 (3) –24 (–14) 12 (6) –16 (–5) 2 (9) –10 (–9) –22 (–13) –29 (–17)
HNO3/NOy 74 (77) 12 (63) 187 (213) 94 (99) 128 (116) 117 (111) 41 (29) 69 (81) 29 (17)
pNO3/NOy 77 (54) 36 (–12) N/A –70 (–71) –38 (–33) –37 (–61) –4 (5) –43 (–35) –28 (–28)
PAN/NOy –34 (–29) –16 (–31) –32 (–39) –30 (–13) –35 (–42) –61 (–56) –19 (–16) 26 (28) –13 (22)
ANs/NOy –46 (–44) –88 (–88) –85 (–82) –66 (–64) –88 (–87) –74 (–72) –3 (2) –57 (–51) 136 (124)

Figure 5 shows simulated O3 mixing ratios from the model ensemble using KORUSv5 emissions for different synoptic regimes compared with DC-8 observations averaged below 1.5 km. DC-8 observations and the ensemble mean of the model results that were sampled along the DC-8 flight track were both regridded and interpolated to 0.25° horizontal resolution to avoid overlapping. Both observed and simulated O3 mixing ratios showed highest values for the transport period and relatively low values during the dynamic weather when frequent clouds and precipitation suppressed photochemical production. We can also see relatively high observed O3 levels during the stagnant and blocking periods throughout the peninsula. However, the ensemble model tends to underestimate the DC-8 observations by 24% and 7% during the two periods, although no significant biases were noted in the model evaluation against the surface observations above. We further discuss this issue in Section 3.2.1.

Figure 5.

Model ensemble (background) and observed (overlaid circles) O3 mixing ratios along the DC-8 flight track averaged below 1.5 km during different synoptic regimes. DC-8 observations and the ensemble model were gridded to 0.25° horizontal resolution. DOI: https://doi.org/10.1525/elementa.2021.00139.f5

Figure 5.

Model ensemble (background) and observed (overlaid circles) O3 mixing ratios along the DC-8 flight track averaged below 1.5 km during different synoptic regimes. DC-8 observations and the ensemble model were gridded to 0.25° horizontal resolution. DOI: https://doi.org/10.1525/elementa.2021.00139.f5

Figure 6 compares simulated PM1 concentrations from the model ensemble and DC-8 observations below 1.5 km averaged for each synoptic regime. As was discussed in Figure 5, we can also see clear differences in PM1 concentrations in South Korea depending on the synoptic meteorology such that the transport period showed the highest PM1 concentrations, whereas the lowest PM1 concentrations were shown in the dynamic weather period, suggesting that the meteorology plays an important role in determining air quality in South Korea. The ensemble model generally captures the spatial distribution of observed PM1 during the dynamic (R = .47) and transport (R = .53) periods. We see a strong west-east gradient of PM1 concentrations during the stagnant period, which is consistently shown in both DC-8 observations and the model ensemble. Predominant easterlies transported air pollutants to the Yellow Sea during this period. During the blocking period, we can see more concentrated PM1 concentrations mostly around the SMA with high emissions from mobile and industrial point sources. The ensemble model tends to pick up the local emissions to some extent, but generally underestimates observations, which calls for a more in-depth analysis on the chemical components of PM, which is discussed in Section 3.2.4.

Figure 6.

Model ensemble (background) and observed (overlaid circles) PM1 concentrations along the DC-8 flight track averaged below 1.5 km during different synoptic regimes. DC-8 observations and the ensemble model were gridded to 0.25° horizontal resolution. DOI: https://doi.org/10.1525/elementa.2021.00139.f6

Figure 6.

Model ensemble (background) and observed (overlaid circles) PM1 concentrations along the DC-8 flight track averaged below 1.5 km during different synoptic regimes. DC-8 observations and the ensemble model were gridded to 0.25° horizontal resolution. DOI: https://doi.org/10.1525/elementa.2021.00139.f6

#### 3.2.1. Gas profiles

Figure 7 compares simulated versus observed profiles of gas species mixing ratios in the SMA, during the KORUS-AQ campaign. A number of gas species were measured aboard the DC-8 during the campaign. Here we selected a few important species including CO, O3, NOx, toluene, isoprene, and HCHO, which are essential for understanding O3 chemistry in the peninsula (Schroeder et al., 2020). As discussed above, individual flight observations and coherently sampled model results were averaged for four periods based on the synoptic patterns during the campaign.

Figure 7.

Comparison of simulated and observed mean vertical profiles of gas species (CO, O3, NOx, toluene, isoprene, HCHO) in the SMA (37–37.6°N, 126.6–127.7°E) during different synoptic conditions. Boxplots and whiskers indicate the interquartile and 10–90 percentile ranges of observations. Vertical lines and closed black circles each indicate the median and mean values. Different colors represent the mean vertical profile of each model with solid lines in red indicating the model ensemble. Surface O3 mixing ratios from the AirKorea network (closed black circles) and the model ensembles (open red circles) are also shown for comparison. Observed and simulated profiles were binned with 0.5 km intervals and the number of data points used for calculating average values for each level is denoted in the right-most panel. DOI: https://doi.org/10.1525/elementa.2021.00139.f7

Figure 7.

Comparison of simulated and observed mean vertical profiles of gas species (CO, O3, NOx, toluene, isoprene, HCHO) in the SMA (37–37.6°N, 126.6–127.7°E) during different synoptic conditions. Boxplots and whiskers indicate the interquartile and 10–90 percentile ranges of observations. Vertical lines and closed black circles each indicate the median and mean values. Different colors represent the mean vertical profile of each model with solid lines in red indicating the model ensemble. Surface O3 mixing ratios from the AirKorea network (closed black circles) and the model ensembles (open red circles) are also shown for comparison. Observed and simulated profiles were binned with 0.5 km intervals and the number of data points used for calculating average values for each level is denoted in the right-most panel. DOI: https://doi.org/10.1525/elementa.2021.00139.f7

CO is a primary pollutant and a good indicator for pollutant emissions from fossil and biofuel burning. In the surface air evaluation, the ensemble model underestimated observed CO mixing ratios in surface air (Figure 3) in South Korea. Most models show an underestimation in comparison to the DC-8 observations, especially in the boundary layer (<1.5 km) regardless of synoptic patterns, indicating that the current emissions inventory may have a missing source of CO in the peninsula. This low bias in the boundary layer is largest in the transport period, when transboundary transport of Chinese emissions heavily influenced the peninsula, and the bias implies that the KORUSv5 emissions may also have some missing sources for CO in China. Gaubert et al. (2020) used MOPITT satellite retrievals to optimize emissions in the CAM-Chem model, starting with the KORUSv5 inventory, and suggested 33% and 80% increases in CO emissions for central and northern China, respectively.

Simulated ensemble CO mixing ratios in the free troposphere do not show any significant biases from the observations, but individual models show different signs in biases from the observed profile. We find that two WRF-Chem simulations by Iowa and UCLA, which used identical chemical boundary conditions generated from the MOZART-4 global model, show higher CO mixing ratios in the free troposphere relative to the observations. Considering the relatively long lifetime of CO, the large variability of simulated CO mixing ratios in the free troposphere could be related to lateral boundary conditions of CO, especially for regional models (Tang et al., 2007). CO profiles from the MOZART-4 boundary conditions were found to be higher than the boundary conditions used in other models by at least 40 ppbv in the free troposphere (Figure S1).

We also find that the planetary boundary layer (PBL) heights in the models could result in notable differences in simulated CO mixing ratios. Figure S2 compares the simulated versus the lidar-derived diurnally varying PBL heights at Seoul National University averaged for the campaign. Some models (CAM-Chem, CAMx, PNU WRF-Chem) overestimate daytime PBL heights relative to the lidar-derived value by 39–60% and simulate lower surface CO mixing ratios by approximately 20% compared to the models with low daytime PBL heights (NCAR WRF-Chem, Iowa WRF-Chem, UCLA WRF-Chem, CMAQ). In addition, biogenic and biomass burning sources and the chemical formation of CO along with different loss rates by the OH oxidation could contribute to the variability of simulated CO mixing ratios among the models.

Tropospheric O3, especially in surface air, is an important air pollutant. In East Asia, seasonal O3 levels are highest in spring (Li et al., 2007). During the campaign, we observed that synoptic meteorology had a distinct influence on O3 variations in the boundary layer with relatively low O3 levels during the dynamic weather and high O3 levels during the stagnant and transport periods (Peterson et al., 2019). Despite this variation, the observed mean O3 mixing ratios from the DC-8 are always higher than 60 ppbv, which is the 8-h average O3 air quality criteria in South Korea. This is somewhat different from what the surface observations showed in Figure 4, that the observed daytime O3 mixing ratios in surface air were often recorded higher than 60 ppbv but not for the whole time. We will discuss this issue of disparity in Section 4.

Figure 7 compares observed versus simulated O3 profiles up to 4 km averaged for different synoptic regimes. We find that the model ensemble is systematically lower than the DC-8 observations throughout the whole troposphere. The magnitudes of low biases in the model ensemble slightly vary with different periods but exist for the whole campaign period. We find that simulated O3 profiles are also dependent on the lateral boundary conditions. The PNU and UCLA WRF-Chem simulations used the same chemical mechanism (RACM-ESRL) with similar model configurations, but they tend to be on opposite extremes of free tropospheric O3, driven by different lateral O3 boundary conditions. O3 boundary condition profiles used in the UCLA WRF-Chem were approximately 10 ppbv higher throughout the troposphere than in the PNU WRF-Chem (Figure S1). Also, as the NCAR WRF-Chem used CAM-Chem results for their boundary conditions, similar profiles are shown in the free troposphere among the two models.

Figure 7 also shows comparisons of important O3 precursors including NOx and VOCs in South Korea. Simpson et al. (2020) showed using VOCs observations during the KORUS-AQ campaign and their OH reactivity that isoprene and toluene are two of the most important VOCs in the SMA. The comparisons of the model ensemble with observed profiles of O3 precursors, including NOx, toluene, and isoprene, during the KORUS-AQ campaign do not show any significant systematic biases (–10 to 12%) in the SMA (Figure 7, Table 6). However, a close investigation shows discrepancies between the models and the observations for the O3 precursors. We will address this issue in detail in Section 4.

#### 3.2.2. Aerosol profiles

Figure 8 shows simulated and observed profiles of PM chemical components, including carbonaceous aerosols (BC, OA), and inorganic sulfate (SO42–), nitrate (NO3), and ammonium (NH4+) aerosols, for four different periods during the campaign. Chemical components of PM1 were measured onboard the DC-8 using the AMS and HD-SP2 and used here for the model evaluation. We defined simulated and observed PM1 as the sum of the above five chemical components. Natural and anthropogenic dust and sea salt aerosols were not considered in this work. During the campaign, there were a few days with high soil dust concentrations transported from the Gobi Desert and nearby arid regions in China. However, not all models simulated natural soil dust aerosols, and therefore, this work does not include soil dust aerosols in the evaluation.

Figure 8.

Comparison of simulated and observed mean vertical profiles of PM1 and its components (BC, OA, SO42–, NO3, NH4+) in the SMA during different synoptic conditions. Boxplots and whiskers indicate the interquartile and 10–90 percentile ranges of observations. Vertical lines and closed black circles each indicate the median and mean values. Different colors represent the mean vertical profile of each model with solid lines in red indicating the model ensemble. Surface PM2.5 concentrations from the AirKorea network (closed black circles) and the model ensembles (open red circles) are also shown for comparison. Data sampling and averaging methods are identical to Figure 7. DOI: https://doi.org/10.1525/elementa.2021.00139.f8

Figure 8.

Comparison of simulated and observed mean vertical profiles of PM1 and its components (BC, OA, SO42–, NO3, NH4+) in the SMA during different synoptic conditions. Boxplots and whiskers indicate the interquartile and 10–90 percentile ranges of observations. Vertical lines and closed black circles each indicate the median and mean values. Different colors represent the mean vertical profile of each model with solid lines in red indicating the model ensemble. Surface PM2.5 concentrations from the AirKorea network (closed black circles) and the model ensembles (open red circles) are also shown for comparison. Data sampling and averaging methods are identical to Figure 7. DOI: https://doi.org/10.1525/elementa.2021.00139.f8

BC is one of the important PM components in East Asia and is mainly emitted from incomplete combustion (Bond et al., 2013). BC emissions from East Asia amount up to 30% of global anthropogenic BC emissions (Bond et al., 2004), with 91% originating from China (Zhang et al., 2009). Major sources of BC are fossil fuel and biofuel use in East Asia. Although the majority of Chinese BC emissions is from residential and industrial sectors, emissions from transportation are dominant in South Korea. In particular, diesel vehicle registrations in South Korea were 47.1% of total vehicles in 2014, and on-road diesel emissions accounted for roughly 60% of total BC emissions from transportation in 2015 (Anenberg et al., 2019). From KORUSv1 to KORUSv5, there was an 80% increase in domestic BC emissions and a 35% decrease in Chinese BC emissions.

We find that the model ensemble tends to overestimate observed BC concentrations, especially in the boundary layer for all four periods (18–44%). The discrepancy slightly varies with different synoptic patterns and is largest during the blocking period. The ensemble model using KORUSv1 emissions showed a 34% underestimation of BC in the boundary layer, which implies a possible overestimation of BC emissions in South Korea in the KORUSv5 inventory.

Simulated PBL heights also play an important role in simulated BC concentrations. In particular, the inter-model variability of BC concentrations in the boundary layer reflect that of PBL heights as shown in Figure S2. For example, the regional models simulating relatively high BC concentrations (NCAR WRF-Chem, Iowa WRF-Chem, UCLA WRF-Chem, CMAQ) show relatively lower PBL heights during daytimes compared to other models. Similarly, the models with high PBL heights, PNU WRF-Chem, CAM-Chem, and CAMx, simulate low BC concentrations. Although the variation of PBL heights among the models does not entirely explain the inter-model variability of simulated BC concentrations, the effect of PBL heights on BC is more noticeable than that of CO.

All models simulated both POA and SOA explicitly. The sum of POA and SOA concentrations was defined as OA mass concentrations, which was compared with the observed AMS OA collected by the DC-8. All models used the same anthropogenic POA emissions but different biomass burning emissions, which accounted for less than 20% of total South Korean POA emissions during the campaign. We find in the comparison that simulated total OA concentrations in most models fall within the observed 10–90 percentile range but with systematic low biases, which have been addressed in previous studies (Heald et al., 2005; Heald et al., 2011). Figure 8 reveals that the large OA variability among the models is mainly due to different treatments of SOA formation from its precursor species, which are summarized in Table 3.

Inorganic sulfate-nitrate-ammonium (SO42–-NO3-NH4+) aerosols were the largest contributor to PM2.5 concentrations during the campaign in South Korea (Jordan et al., 2020). It appears that the model ensemble generally well reproduces the profiles of observed inorganic aerosols for the stagnant and transport periods but shows overestimation during the dynamic weather period and underestimation during the blocking period. The discrepancy for the dynamic period could be associated with the treatment of wet scavenging processes in models (Itahashi et al., 2020). During the blocking period, observed SO42– is relatively well captured by the ensemble model (NMB = −8%), but NO3 and NH4+ aerosols are underestimated by 38% and 26%, respectively. One thing to note is that the discrepancy between the model ensemble and the observations is a lot smaller than the inter-model variability.

During the dynamic weather period, inorganic components of PM1 are overestimated and OA is underestimated by the ensemble model, resulting in a slight overestimation (1%) of observed PM1 concentration in the boundary layer. The ensemble model captures the enhancement of observed DC-8 PM1 concentrations during the transport period. From this, we can infer that the models were able to simulate the meteorological conditions that caused eastward transport of pollutants to the Korean peninsula.

Overall, the model performances for the stagnant and blocking periods show underestimations of PM1 (–33% and –41%, respectively) below 1.5 km. During the stagnant period, underestimation of PM1 is mostly due to low biases in simulated OA, whereas the underestimation of OA and inorganic aerosols both play a larger part during blocking. We also sampled AirKorea surface PM2.5 observations from the closest points along the DC-8 flight tracks below 1 km so that they are coherent with the DC-8 observations in both time and space. Figure 8d shows that collocated AirKorea ground observations of PM2.5 during the blocking period show a large difference from the DC-8 observations at 0.5 km, where the aircraft observation is much higher than that in surface air. This is likely caused by the mismatch between the aircraft flight tracks versus the PM2.5 monitoring stations (Figure 3e) located in major industrial source regions in the SMA, which is further discussed in Section 3.3.

#### 3.2.3. Evaluation for reactive nitrogen partitioning

Total reactive nitrogen (NOy), defined as the sum of NOx and its oxidation products (NOy = NOx + HNO3 + particulate nitrate + PAN + ANs + etc.), is also an important species in tropospheric chemistry. NOy partitioning between NOx and its reservoirs has an important implication for O3 and aerosol formation (Zellweger et al., 2003). Models were able to reproduce observed NOy mixing ratios but failed to accurately simulate its partitioning among reactive nitrogen compounds. In this section, we compare the observed and simulated vertical profiles of NOy species partitioning for different synoptic conditions during the campaign.

Figure 9 shows the simulated and observed profiles of each nitrogen compound ratio to NOy, including NOx, HNO3, particulate nitrate (pNO3), PAN, and ANs. First, we find that the model ensemble overestimates the observed NOy mixing ratios in the boundary layer by 7% during the whole campaign (Figure S3) but well captures the gradual increase of observed NOy from May 1 to June 10 (not shown).

Figure 9.

Comparison of simulated and observed mean vertical profiles of NOy (= NOx + HNO3 + particulate nitrate + PAN + ANs) partitioning in the SMA during different synoptic conditions. Boxplots and whiskers indicate the interquartile and 10–90 percentile ranges of observations. Vertical lines and closed black circles each indicate the median and mean values. Different colors represent the mean vertical profile of each model with solid lines in red indicating the model ensemble. Data sampling and averaging methods are identical to Figure 7. DOI: https://doi.org/10.1525/elementa.2021.00139.f9

Figure 9.

Comparison of simulated and observed mean vertical profiles of NOy (= NOx + HNO3 + particulate nitrate + PAN + ANs) partitioning in the SMA during different synoptic conditions. Boxplots and whiskers indicate the interquartile and 10–90 percentile ranges of observations. Vertical lines and closed black circles each indicate the median and mean values. Different colors represent the mean vertical profile of each model with solid lines in red indicating the model ensemble. Data sampling and averaging methods are identical to Figure 7. DOI: https://doi.org/10.1525/elementa.2021.00139.f9

The observed NOx/NOy ratio is in the range of 0.1–0.8, which is largest at the surface and generally decreases with altitude. The models relatively well capture the observed ratio in the low troposphere but generally underestimate it in the free troposphere (>3 km), implying too rapid conversion of NOx to its oxidation products, particularly HNO3.

We find that the discrepancy between the models and the observations is quite large for the HNO3/NOy ratio. Most models tend to overestimate the observed HNO3 fractions (0.1–0.3) by factors of 2–3. This high bias is in general large in the free troposphere, which was also found in the model evaluation against the observations from the TRACE-P and ACE-Asia campaigns in East Asia (Carmichael et al., 2003; Tang et al., 2004).

Figure S3 shows that simulated HNO3 mixing ratios show systematic high biases compared to HNO3 measured using the California Institute of Technology Chemical Ionization Mass Spectrometer (CIT-CIMS) instrument, used in Figure 9. The CIT-CIMS measures gas phase HNO3 through selective ion chemical ionization with an uncertainty of 40% (Crounse et al., 2006). In situ HNO3 was also measured by the University of New Hampshire Soluble Acidic Gases and Aerosols (SAGA) instrument onboard the DC-8 (Dibb et al., 2003). The SAGA HNO3 measurements are likely to have a considerable enhancement due to the contamination by submicron aerosol nitrate (McNaughton et al., 2009). When simulated HNO3 is compared to the SAGA measurements, the NMB of the ensemble model is reduced to –24%, relative to the comparison with the CIT-CIMS HNO3 (159%).

A large fraction of NOy is composed of organic nitrogen compounds (PAN and ANs) in the observations during the KORUS-AQ. Especially the observed PAN fraction amounts up to 0.5 above in the free troposphere, indicating large availability of oxidized VOCs. The model ensemble tends to underestimate the observed PAN fractions with large inter-model variability. We also find that the models in general underestimate the ANs fractions. The ensemble model underestimates ANs by 39%, which only make up 3% of total NOy, while the ANs contribute up to 6% of observed total NOy. In particular, models with relatively coarse resolutions showed large underestimations of ANs (Figure S3).

The underestimation of the organic nitrogen fraction may imply the importance of reactive VOCs in NOy partitioning and NOx-HOx recycling. For example, more detailed aromatic chemistry mechanisms increased the conversion of NOx to organic reservoirs such as PAN and ANs, by the oxidation of aromatic VOCs (Oak et al., 2019). As different VOC chemistry mechanisms are used in each model, large inter-model variabilities are noticeable in simulated organic nitrogen fractions.

### 3.3. Evaluation for aerosol compositions

Variations in the aerosol chemical composition during the campaign were distinct among different synoptic patterns, which were examined by previous studies to explain the elevations in PM loadings. Kim et al. (2018) compared the temporal variations of inorganic/organic aerosols and their precursor gases at the KIST ground site in Seoul. They found that the organic portion rapidly increased in May 20–24 under stagnant conditions with high temperature and strong radiation, and high levels of VOCs made a favorable environment for secondary formation of OA. They also concluded that during the high loading episode (May 26–31) when relative compositions of inorganics were dominant, direct eastward transport from SO2 source regions was combined with local secondary formation of nitrate.

Jordan et al. (2020) also provided a detailed analysis of aerosol observations and revealed that the majority of fine PM (PM2.5 or PM1) in South Korea was composed of secondary aerosols including inorganic sulfate-nitrate-ammonium and OA. The information of PM chemical composition is important to develop an efficient air quality policy to reduce PM concentrations in South Korea. Here we evaluate the models focusing on this aspect by comparing the simulations against the AMS observations aboard the DC-8 and in surface air during the campaign.

Figure 10 compares the simulated and observed PM1 chemical compositions along the DC-8 flight tracks (Figure 6) below 1.5 km over the SMA. We here defined the hydrocarbon-like component of OA as POA, and the oxidized component as SOA, based on PMF analysis by Nault et al. (2018). The observed PM1 mass concentrations vary with the meteorological regimes as we discussed earlier that the highest PM1 concentration occurred in the transport period and the lowest concentration was observed in the dynamic weather. Values are also high in the blocking pattern period due to the accumulation of local emissions from mobile and industrial sources in the Gyeonggi province (Figure S4).

Figure 10.

Comparison of simulated and observed mean PM1 chemical compositions below 1.5 km for different synoptic regimes in the SMA during the campaign. Each chemical component is indicated in different colors and the observations are from AMS aboard the DC-8 aircraft. Model results were sampled along the DC-8 flight tracks so that they are spatially and temporally coherent. DOI: https://doi.org/10.1525/elementa.2021.00139.f10

Figure 10.

Comparison of simulated and observed mean PM1 chemical compositions below 1.5 km for different synoptic regimes in the SMA during the campaign. Each chemical component is indicated in different colors and the observations are from AMS aboard the DC-8 aircraft. Model results were sampled along the DC-8 flight tracks so that they are spatially and temporally coherent. DOI: https://doi.org/10.1525/elementa.2021.00139.f10

This observed variability is generally captured by the model ensemble, but it fails to reproduce high levels of PM concentrations in the blocking period, during which DC-8 observations were considerably affected by domestic emissions from major industrial complexes located in the southern part of the SMA (Figure S4). The observed SOA concentrations and toluene, one of the main parent hydrocarbons, show relatively high values in this area, where PM1 concentrations are also high. Although these local sources are included in the present inventory, top-down estimations using HCHO show that toluene emissions in the KORUSv5 inventory are underestimated from these facilities (Fried et al., n.d.; Kwon et al., n.d.).

As to the chemical composition evaluation, we find a few issues as follows:

1. The models in general overestimate POA and underestimate SOA throughout the whole campaign (Figure 10).

2. Enhancements of the organic fraction during the stagnant and blocking periods are not captured by most models, especially for SOA.

3. Models that use the simplified SOA scheme (M1, M4, M6) that scales CO emissions to estimate SOA precursors generally underestimate local SOA production, which is in part due to CO underestimation.

4. Models that use the VBS scheme (M2, M5, M7) generally simulate higher SOA concentrations than those of other schemes, resulting in closer agreement with observations, but still underpredict SOA concentrations during the stagnant and blocking periods.

5. The models successfully reproduce increases of the inorganic fraction (especially sulfate) during the transport period.

6. When the models appear to be lower than the observations, it is mainly because of the OA underestimation.

We conducted a similar PM chemical composition evaluation using observations in surface air. Figure 11 compares the observed and simulated PM1 chemical composition at the KIST ground site (Figure 2) in Seoul during the campaign. For OA speciation, we defined POA as the sum of hydrocarbon-like OA and cooking OA, and SOA as the sum of oxidized OA, according to Kim et al. (2018). The observations show a similar temporal variability with aircraft observations, being higher in the transport period relative to other periods.

Figure 11.

Comparison of simulated and observed mean PM1 chemical compositions in surface air for different synoptic regimes during the campaign. Each chemical component is indicated in different colors and the observations are from AMS at the KIST ground site (37.601°N, 127.045°E). Model results from the lowest model layer were sampled coherently with the observations. DOI: https://doi.org/10.1525/elementa.2021.00139.f11

Figure 11.

Comparison of simulated and observed mean PM1 chemical compositions in surface air for different synoptic regimes during the campaign. Each chemical component is indicated in different colors and the observations are from AMS at the KIST ground site (37.601°N, 127.045°E). Model results from the lowest model layer were sampled coherently with the observations. DOI: https://doi.org/10.1525/elementa.2021.00139.f11

Underestimation of SOA during the stagnant period, when local secondary production was likely to dominate OA concentrations, is shown in both the airborne (Figure 10b) and surface (Figure 11b) data. Nault et al. (2018) and Kim et al. (2018) showed that local emissions and chemistry during the campaign resulted in intense formation of SOA in Seoul. The influence of Siberian wildfires on May 18 also brought a mixture of aged smoke plumes with additional SOA precursors to Korea (Peterson et al., 2019). This implies that the current SOA schemes in the models tend to underestimate SOA formation from precursor emissions.

One notable difference from the DC-8 observation is that the observed PM1 concentration is much lower in the blocking period than that of the DC-8 observation. Figure S5 shows that similar temporal variations of surface PM at the KIST site are also observed at the Olympic park ground site, which is located in southeast Seoul. As shown in Figure S4, the KIST and Olympic park sites are located in the northern part of the SMA, where the influence of industrial emissions is relatively small. The models tend to overestimate the ground observations during the blocking period, mainly driven by the overestimation of inorganic aerosols. We find in particular that the nitrate overestimation is the main reason for the model overestimation of the surface observations. The lack of nitrate and ammonium aerosols in CAM-Chem results in an overall underestimation of aerosol mass in comparison to both ground and aircraft observations, indicating the importance of including nitrate aerosol formation in air quality simulations.

## 4. Discussion

As was discussed in the surface air evaluation (Figure 4b), the models showed capability in reproducing observed daytime O3 levels during the campaign with a high correlation coefficient (R = .83) and no significant bias (NMB = 4%) between the model ensemble using the KORUSv5 emissions and the AirKorea observations. However, the model ensemble showed a systematic low bias in simulating DC-8 O3 profiles throughout the whole troposphere regardless of the synoptic patterns. This underestimation was also shown at the lowest altitude, which seems inconsistent with the result of the surface air evaluation. For comparison, we plotted averaged AirKorea surface O3 mixing ratios, which were coherently sampled with DC-8 observations below 1 km to examine the gradient of O3 mixing ratios from the surface to the lowest DC-8 observations. As shown in Figure 7, a sharp O3 gradient from the surface to the lowest DC-8 altitude is found in the observations and the model ensemble, although the simulated gradient is less than that of the observations.

Figure S6 compares O3 sonde measurements at Olympic park and Taehwa research forest, a downwind site located 30 km apart from Olympic park. O3 profiles from sonde measurements show sharp declines (approximately 30–40%) below 2 km to the surface, which were also reported in a case study during summer 2016 at Hangzhou, China, showing a 40% decrease within the boundary layer (Su et al., 2017). Observation-constrained box model simulations for the KORUS-AQ period show that O3 formation, especially through the reaction between HO2 and NO, decreases and the O3 destruction rate increases at lower altitudes near the surface driven by the concentrated NO and VOCs that react with O3, resulting in an overall decrease of net O3 production (Figure S6c). As the box model only considered photochemical loss of O3, combined effects of wet or dry deposition of O3 and its precursors could further lower surface O3 levels.

Among the models, we found that CAM-Chem reproduced observed O3 profiles relatively well with no significant low biases above the boundary layer, compared to other models (Table 6). In order to understand a possible cause for the low bias in the model ensemble, we used the CAM-Chem stratospheric O3 tracer to estimate stratospheric O3 contributions in the troposphere during the campaign (Figure S7). Stratospheric O3 influx varies with different synoptic regimes and shows generally high contributions during stagnant and blocking periods associated with persistent high pressure systems. We also find significant stratospheric O3 contributions even to the surface air (15–30 ppbv) from CAM-Chem simulations, which might be overestimated because this stratospheric tracer does not account for rapid titration by high NOx mixing ratios in surface air and dry deposition processes. Because the model top heights are approximately 20 km for the regional models, they did not simulate stratospheric chemistry, so the stratospheric influx of O3 was determined by the boundary and initial conditions used. Therefore, this additional O3 source may not be accurately represented in some models depending on the synoptic regimes, which show systematic low biases in simulated O3 mixing ratios compared with the observations in spring, when frequent intrusions of stratosphere air into the troposphere occur (Wang et al., 2020). Nonetheless, Kim and Lee (2010) suggested that springtime O3 maxima in Korea were more likely to be caused by photochemistry than stratospheric intrusions. We also found that the simulated stratospheric O3 using a linearized first-order reaction scheme (Linoz) in GEOS-Chem showed similar values to that of CAM-Chem, but showed lower levels in the troposphere, indicating that additional O3 sources were absent in the model.

Schroeder et al. (2020) used an observation-constrained photochemical box model to quantify O3 production sensitivities to various precursor VOCs in Korea. Souri et al. (2020) used a box model to evaluate the HCHO/NOx column ratio obtained from a remote sensing instrument during KORUS-AQ, which can be used as a proxy for classifying O3 production regimes. These studies emphasized the role of reactive VOCs, especially aromatics, in O3 production in the major metropolitan areas in South Korea, which were clearly found to be VOC-limited. Therefore, differences in the aromatic chemistry schemes among models are expected to result in considerable differences in O3 simulation.

Although no systematic biases in the model ensemble were found for reactive VOCs, such as toluene and isoprene, as shown in Figure 7, the simulated NOx concentration below 1.5 km was higher by 11% relative to the observations during the campaign. NOx was overestimated by 19% during the stagnant and blocking periods, when local chemistry had a larger influence on O3 production. Oak et al. (2019) investigated the O3 production sensitivity to aromatic chemistry in Korea during KORUS-AQ using a 3-D CTM and found that aromatic VOCs oxidation drives more NO to NO2 conversion, resulting in significant changes in O3 photochemistry in the boundary layer during 13–16 LST. They found that including detailed aromatic chemistry in the model resulted in a decrease of simulated NOx mixing ratios by 20% while increasing O3 by 13% below 1.5 km. Therefore, we expect that the biases below 1.5 km in the simulated NOx (11%) and O3 (–18%) in the model ensemble can be reduced by including detailed aromatic chemical mechanisms (Knote et al., 2014; Porter et al., 2017). Among the models, CAM-Chem and NCAR WRF-Chem used the most complex aromatic chemistry scheme (56 kinetic reactions) from MOZART-T1 described in Knote et al. (2014) and showed smaller model biases (–15%) compared to the ensemble model (–18%) for O3 below 1.5 km during the whole campaign.

## 5. Conclusions

An international air quality field study, KORUS-AQ, which was jointly hosted by the Korean NIER and the U.S. NASA, occurred during May–June 2016, to understand the factors controlling air quality in South Korea (Crawford et al., n.d.). Extensive aircraft and ground network observations in the peninsula and nearby oceans were available from the campaign for which a number of 3-D CTMs were also used on a daily basis to produce up to 5-day air quality forecasts for planning aircraft observations.

Although the forecasts were valuable for identifying pollution plumes and other features targeted for observational sampling, the forecasts sometimes differed between models and were not able to capture the observed magnitudes of aerosols, O3, and their precursors. This study addressed the issues associated with air quality simulations using model evaluation against extensive surface and aircraft observations from the campaign as well as intercomparisons between models.

Six regional and two global CTMs participated in the MICS and were used to conduct air quality simulations focusing on O3, aerosols, and their precursor species for the campaign using the KORUSv5 anthropogenic emissions inventory. The participating models chose their own options for emissions from biomass burning and biogenic sources.

Relative to the KORUSv1 inventory developed for the campaign, considerable changes were made in the latest KORUSv5 inventory. For example, the anthropogenic CO and NOx emissions in South Korea were increased by factors of 2.5 and 1.4, respectively, but they were decreased by about 25% in East China. A more dramatic decrease up to 70% was shown for SO2 in East China. On the other hand, anthropogenic VOCs emissions were generally increased in both East China and South Korea from v1 to v5, especially for aromatic species, which were shown to be higher by factors of 1.5–2.3 in v5 inventory compared to those of v1.

Our model evaluation focusing on surface air simulations revealed that the models, using the KORUSv5 emissions inventory, successfully reproduced the observed spatial and temporal variabilities of both O3 and PM2.5 concentrations in surface air in South Korea. However, we found a significant low bias for simulated CO mixing ratios in the peninsula and nearby oceans, implying possible missing CO sources in the inventory in East Asia.

Peterson et al. (2019) showed that synoptic meteorology played a critical role in determining characteristics of air quality in South Korea and broadly East Asia during the KORUS-AQ campaign. Following Peterson et al. (2019), we grouped DC-8 observations for four distinct synoptic regimes: dynamic weather (May1–16), stagnant period (May 17–22), transport period (May 24–31), and blocking meteorology (June 1–10), during which we conducted the model evaluation by comparing the simulated versus observed profiles of species concentrations primarily in the SMA.

Although the DC-8 aircraft observations of O3 precursors were well captured by the models, simulated O3 mixing ratios were lower than the observations in the troposphere regardless of synoptic regimes. The low O3 bias in the models was in part owing to too low stratospheric O3 influxes, especially in the middle and upper troposphere and in part owing to insufficient chemical production of O3 because of a simple representation of the chemistry of aromatic VOCs in most models.

Unlike O3, the synoptic meteorology played an important role in determining the observed variability of PM concentrations in South Korea. Highest PM1 concentrations observed onboard DC-8 occurred in South Korea during the transport period driven by the transboundary transport from China, whereas the lowest values were shown in the dynamic weather period. Despite the underestimation of PM1 concentrations in the models, we found that the models generally reproduced the observed variability of the PM1 mass concentrations depending on the synoptic regimes. Each chemical component comprising PM1 in the models, however, showed some discrepancies from the observations.

Observed chemical compositions of both PM1 and PM2.5 concentrations from the DC-8 and the surface network, respectively, showed a significant contribution by inorganic SO42–-NO3-NH4+ aerosols followed by carbonaceous aerosols including organic and BC aerosols. The large inorganic aerosol contribution was in general captured by the models. However, the ensemble of models tends to overestimate observed BC concentrations, especially in the boundary layer for all four synoptic periods (18–44%) with the KORUSv5 inventory, implying a possible overestimation of BC emissions in South Korea. We find in this study that the systematic low biases of simulated OA concentrations addressed in previous literature were also shown in the participating models with less degree of discrepancy from the observations, but with large OA variability among the models mainly due to different treatments of SOA formation from its precursor species.

Model deficiencies in simulating organic nitrogen compounds (e.g., ANs) in the boundary layer were also revealed in the model evaluation against DC-8 observations. The underestimation of the organic nitrogen fraction out of total reactive nitrogen (NOy) might imply the importance of reactive VOCs in NOy partitioning and O3 formation through NOx-HOx recycling.

From the model evaluation, we found that an ensemble of model results, incorporating individual models with differing strengths and weaknesses, performs better than most individual models at representing observed atmospheric composition for the campaign. Ongoing model development and evaluation, in close collaboration with emissions inventory development, are needed to improve air quality forecasting.

## Data accessibility statement

Observational data from KORUS-AQ used in this study can be downloaded in the International Consortium for Atmospheric Research on Transport and Transformation (ICARTT) format through the data archive website (https://www-air.larc.nasa.gov/cgi-bin/ArcView/korusaq). Hourly surface observations from the AirKorea network can be downloaded online through http://www.airkorea.or.kr/web. All data from the participating models for this study are also available through the KORUS-AQ data archive website in the ICARTT format for aircraft data and the netCDF format for hourly surface data.

## Supplemental files

Figure S1. Comparison of average vertical profiles of CO and O3 mixing ratios from the lateral boundary conditions used in each model.

Figure S2. Comparison of mean diurnal profiles of observed (lidar) and simulated PBL height at Seoul National University (SNU; 37.46°N, 126.95°E) during the whole campaign period.

Figure S3. Comparison of simulated and observed mean vertical profiles of NOy and its components (HNO3, PAN, ANs) in the SMA (37–37.6°N, 126.6–127.7°E) during different synoptic conditions.

Figure S4. Model ensemble (background) and observed (overlaid circles) a) PM1, b) secondary organic aerosol (SOA) concentrations, c) toluene, and d) HCHO mixing ratios along the DC-8 flight track averaged below 1.5 km in the SMA during the blocking period (June 1–10).

Figure S5. Comparison of simulated and observed mean PM2.5 chemical compositions in surface air at the Olympic park ground site (37.519°N, 127.122°E) for different synoptic regimes during the campaign.

Figure S6. Comparison of O3 sonde observations at a) Olympic park, b) Taehwa research forest (37.280°N, 127.227°E), and c) 0-D photochemical boxmodel simulations along the DC-8 flight track in the SMA below 2 km.

Figure S7. Simulated CAM-Chem (magenta) O3 mixing ratios during KORUS-AQ with stratospheric O3 (blue) contributions.

## Acknowledgments

We thank all members of the KORUS-AQ instrument team for their contributions during the field campaign and the agencies operating the measurements onboard DC-8 and ground sites. This work was supported by the Korean Government Ministry of Science and ICT (MSIT) and the National Research Foundation of Korea (NRF) No. 2018004494. This material is based upon work supported by the National Center for Atmospheric Research (NCAR), which is a major facility sponsored by the National Science Foundation (NSF) under Cooperative Agreement No. 1852977. The CESM project is supported primarily by the NSF. Computing and data storage resources, including the Cheyenne supercomputer (doi:10.5065/D6RX99HX), were provided by the Computational and Information Systems Laboratory (CISL) at NCAR.

## Funding

This study was supported and grant funded by the Korean Government Ministry of Science and ICT (MSIT), the National Research Foundation of Korea (NRF) No. 2018004494 (RJP, YJO), and the U.S. National Aeronautics and Space Administration (NASA) No. NNX16AD96G (LKE, GGP, BG, PES) and NNX15AU17G (GRC).

## Competing interests

The authors have declared that no competing interests exist.

## Author contributions

Contributed to conception and design: RJP, YJO, LKE.

Contributed to acquisition of data: All coauthors.

Contributed to analysis and interpretation of data: RJP, YJO, LKE, GGP, GRC, PES, BG.

Drafted and/or revised the article: RJP, YJO, LKE, GGP, PES, JHC, COS.

Approved the submitted version for publication: All co-authors.

## References

Ackermann
,
I
,
Hass
,
H
,
Memmesheimer
,
M
,
Ebel
,
A
,
Binkowski
,
F
,
Shankar
,
U
.
1998
.
Modal aerosol dynamics model for Europe: Development and first applications
.
Atmospheric Environment
17
:
2981
2999
.
,
R
,
McKeen
,
S
,
Robinson
,
A
,
Bahreini
,
R
,
Middlebrook
,
A
,
de Gouw
,
JA
,
Meagher
,
J
,
Hsie
,
EY
,
Edgerton
,
E
,
Shaw
,
S
,
Trainer
,
M
.
2012
.
A volatility basis set model for summertime secondary organic aerosols over the eastern United States in 2006
.
Journal of Geophysical Research (Atmospheres)
117
:
6301
. DOI: http://dx.doi.org/10.1029/2011JD016831.
,
R
,
McKeen
,
S
,
Trainer
,
M
,
Banta
,
R
,
Brewer
,
A
,
Brown
,
S
,
Edwards
,
PM
,
de Gouw
,
JA
,
Frost
,
GJ
,
Gilman
,
J
,
Helmig
,
D
,
Johnson
,
B
,
Karion
,
A
,
Koss
,
A
,
Langford
,
A
,
Lerner
,
B
,
Olson
,
J
,
Oltmans
,
S
,
Peischl
,
J
,
Pétron
,
G
,
Pichugina
,
Y
,
Roberts
,
JM
,
Ryerson
,
T
,
Schnell
,
R
,
Senff
,
C
,
Sweeney
,
C
,
Thompson
,
C
,
Veres
,
PR
,
Warneke
,
C
,
Wild
,
R
,
Williams
,
EJ
,
Yuan
,
B
,
Zamora
,
R
.
2015
.
Understanding high wintertime ozone pollution events in an oil- and natural gas-producing region of the western US
.
Atmospheric Chemistry and Physics
15
(
1
):
411
429
. DOI: http://dx.doi.org/10.5194/acp-15-411-2015.
Akimoto
,
H
,
Nagashima
,
T
,
Li
,
J
,
Fu
,
JS
,
Ji
,
D
,
Tan
,
J
,
Wang
,
Z
.
2019
.
Comparison of surface ozone simulation among selected regional models in MICS-Asia III – Effects of chemistry and vertical transport for the causes of difference
.
Atmospheric Chemistry and Physics
19
(
1
):
603
615
. DOI: http://dx.doi.org/10.5194/acp-19-603-2019.
Anenberg
,
S
,
Miller
,
J
,
Henze
,
D
,
Minjares
,
R
.
2019
.
A global snapshot of the air pollution-related health impacts of transportation sector emissions in 2010 and 2015
.
ICCT report
.
Washington, DC
.
Available at
https://www.theicct.org/publications/health-impacts-transport-emissions-2010-2015.
Accessed 10 September, 2020
.
Bond
,
TC
,
Doherty
,
SJ
,
Fahey
,
DW
,
Forster
,
PM
,
Berntsen
,
T
,
DeAngelo
,
BJ
,
Flanner
,
MG
,
Ghan
,
S
,
Kärcher
,
B
,
Koch
,
D
,
Kinne
,
S
,
Kondo
,
Y
,
Quinn
,
PK
,
Sarofim
,
MC
,
Schultz
,
MG
,
Schulz
,
M
,
Venkataraman
,
C
,
Zhang
,
H
,
Zhang
,
S
,
Bellouin
,
N
,
Guttikunda
,
SK
,
Hopke
,
PK
,
Jacobson
,
MZ
,
Kaiser
,
JW
,
Klimont
,
Z
,
Lohmann
,
U
,
Schwarz
,
JP
,
Shindell
,
D
,
Storelvmo
,
T
,
Warren
,
SG
,
Zender
,
CS
.
2013
.
Bounding the role of black carbon in the climate system: A scientific assessment
.
Journal of Geophysical Research: Atmospheres
118
(
11
):
5380
5552
. DOI: http://dx.doi.org/10.1002/jgrd.50171.
Bond
,
TC
,
Streets
,
DG
,
Yarber
,
KF
,
Nelson
,
SM
,
Woo
,
J-H
,
Klimont
,
Z
.
2004
.
A technology-based global inventory of black and organic carbon emissions from combustion
.
Journal of Geophysical Research: Atmospheres
109
(
D14
). DOI: http://dx.doi.org/10.1029/2003jd003697.
Carlton
,
AG
,
Bhave
,
PV
,
Napelenok
,
SL
,
Edney
,
EO
,
Sarwar
,
G
,
Pinder
,
RW
,
Pouliot
,
GA
,
Houyoux
,
M
.
2010
.
Model representation of secondary organic aerosol in CMAQv4.7
.
Environmental Science & Technology
44
(
22
):
8553
8560
. DOI: http://dx.doi.org/10.1021/es100636q.
Carmichael
,
GR
,
Hayami
,
H
,
Calori
,
G
,
Uno
,
I
,
Cho
,
SY
,
Engardt
,
M
,
Kim
,
S-B
,
Ichikawa
,
Y
,
Ikeda
,
Y
,
Ueda
,
H
,
Amann
,
M
.
2001
.
Model intercomparison study of long range transport and sulfur deposition in East Asia (MICS-ASIA)
.
Water, Air, and Soil Pollution
130
(
1
):
51
62
. DOI: http://dx.doi.org/10.1023/A:1012291200633.
Carmichael
,
GR
,
Tang
,
Y
,
Kurata
,
G
,
Uno
,
I
,
Streets
,
D
,
Woo
,
J-H
,
Huang
,
H
,
Yienger
,
J
,
Lefer
,
B
,
Shetter
,
R
,
Blake
,
D
,
Atlas
,
E
,
Fried
,
A
,
Apel
,
E
,
Eisele
,
F
,
Cantrell
,
C
,
Avery
,
M
,
Barrick
,
J
,
Sachse
,
G
,
Brune
,
W
,
Sandholm
,
S
,
Kondo
,
Y
,
Singh
,
H
,
Talbot
,
R
,
Bandy
,
A
,
Thorton
,
D
,
Clarke
,
A
,
Heikes
,
B
.
2003
.
Regional-scale chemical transport modeling in support of the analysis of observations obtained during the TRACE-P experiment
.
Journal of Geophysical Research: Atmospheres
108
(
D21
). DOI: http://dx.doi.org/10.1029/2002jd003117.
Carmichael
,
GR
,
Ueda
,
H
.
2008
.
MICS-Asia II: The model intercomparison study for Asia phase II
.
Atmospheric Environment
42
(
15
):
3465
3467
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2007.10.003.
Carter
,
W
.
2000
.
Documentation of the SAPRC-99 Chemical Mechanism for VOC Reactivity Assessment
.
Final Report to California Air Resources Board
.
Available at
https://intra.cert.ucr.edu/∼carter/absts.htm#saprc99.
Chen
,
L
,
Gao
,
Y
,
Zhang
,
MG
,
Fu
,
JS
,
Zhu
,
J
,
Liao
,
H
,
Li
,
JL
,
Huang
,
K
,
Ge
,
BZ
,
Wang
,
XM
,
Lam
,
YF
,
Lin
,
CY
,
Itahashi
,
S
,
Nagashima
,
T
,
Kajino
,
M
,
Yamaji
,
K
,
Wang
,
ZF
,
Kurokawa
,
J
.
2019
.
MICS-Asia III: Multi-model comparison and evaluation of aerosol over East Asia
.
Atmospheric Chemistry and Physics
19
(
18
):
11911
11937
. DOI: http://dx.doi.org/10.5194/acp-19-11911-2019.
Choi
,
J
,
Park
,
RJ
,
Lee
,
H-M
,
Lee
,
S
,
Jo
,
DS
,
Jeong
,
JI
,
Henze
,
DK
,
Woo
,
J-H
,
Ban
,
S-J
,
Lee
,
M-D
,
Lim
,
C-S
,
Park
,
M-K
,
Shin
,
HJ
,
Cho
,
S
,
Peterson
,
D
,
Song
,
C-K
.
2019
.
Impacts of local vs. trans-boundary emissions from different sectors on PM2.5 exposure in South Korea during the KORUS-AQ campaign
.
Atmospheric Environment
203
:
196
205
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2019.02.008.
Colman
,
JJ
,
Swanson
,
AL
,
Meinardi
,
S
,
Sive
,
BC
,
Blake
,
DR
,
Rowland
,
FS
.
2001
.
Description of the analysis of a wide range of volatile organic compounds in whole air samples collected during PEM-Tropics A and B
.
Analytical Chemistry
73
(
15
):
3723
3731
. DOI: http://dx.doi.org/10.1021/ac010027g.
Crawford
,
JH
,
Ahn
,
J-Y
,
,
J
,
Chang
,
L
,
Emmons
,
L
,
Kim
,
J
,
Lee
,
G
,
Park
,
J-H
,
Park
,
R
,
Woo
,
JH
,
Lefer
,
B
.
n.d.
The Korea-United States Air Quality (KORUS-AQ) Field Study
,
in press
.
Crounse
,
JD
,
McKinney
,
KA
,
Kwan
,
AJ
,
Wennberg
,
PO
.
2006
.
Measurement of gas-phase hydroperoxides by chemical ionization mass spectrometry
.
Analytical Chemistry
78
(
19
):
6726
6732
. DOI: http://dx.doi.org/10.1021/ac0604235.
Cubison
,
MJ
,
Ortega
,
AM
,
Hayes
,
PL
,
Farmer
,
DK
,
Day
,
D
,
Lechner
,
MJ
,
Brune
,
WH
,
Apel
,
E
,
Diskin
,
GS
,
Fisher
,
JA
,
Fuelberg
,
HE
,
Hecobian
,
A
,
Knapp
,
DJ
,
Mikoviny
,
T
,
Riemer
,
D
,
Sachse
,
GW
,
Sessions
,
W
,
Weber
,
RJ
,
Weinheimer
,
AJ
,
Wisthaler
,
A
,
Jimenez
,
JL
.
2011
.
Effects of aging on organic aerosol from open biomass burning smoke in aircraft and laboratory studies
.
Atmospheric Chemistry and Physics
11
(
23
):
12049
12064
. DOI: http://dx.doi.org/10.5194/acp-11-12049-2011.
Day
,
DA
,
Wooldridge
,
PJ
,
Dillon
,
MB
,
Thornton
,
JA
,
Cohen
,
RC
.
2002
.
A thermal dissociation laser-induced fluorescence instrument for in situ detection of NO2, peroxy nitrates, alkyl nitrates, and HNO3
.
Journal of Geophysical Research: Atmospheres
107
(
D6
):
ACH 4-1–ACH 4-14
. DOI: http://dx.doi.org/10.1029/2001JD000779.
DeCarlo
,
PF
,
Kimmel
,
JR
,
Trimborn
,
A
,
Northway
,
MJ
,
Jayne
,
JT
,
Aiken
,
AC
,
Gonin
,
M
,
Fuhrer
,
K
,
Horvath
,
T
,
Docherty
,
KS
,
Worsnop
,
DR
,
Jimenez
,
JL
.
2006
.
Field-deployable, high-resolution, time-of-flight aerosol mass spectrometer
.
Analytical Chemistry
78
(
24
):
8281
8289
. DOI: http://dx.doi.org/10.1021/ac061249n.
Dibb
,
JE
,
Talbot
,
RW
,
Scheuer
,
EM
,
Seid
,
G
,
Avery
,
MA
,
Singh
,
HB
.
2003
.
Aerosol chemical composition in Asian continental outflow during the TRACE-P campaign: Comparison with PEM-West B
.
Journal of Geophysical Research: Atmospheres
108
(
D21
). DOI: http://dx.doi.org/10.1029/2002jd003111.
Donahue
,
NM
,
Robinson
,
AL
,
Stanier
,
CO
,
Pandis
,
SN
.
2006
.
Coupled partitioning, dilution, and chemical aging of semivolatile organics
.
Environmental Science & Technology
40
(
8
):
2635
2643
. DOI: http://dx.doi.org/10.1021/es052297c.
Emmons
,
LK
,
Schwantes
,
RH
,
Orlando
,
JJ
,
Tyndall
,
G
,
Kinnison
,
D
,
Lamarque
,
J-F
,
Marsh
,
D
,
Mills
,
MJ
,
Tilmes
,
S
,
Bardeen
,
C
,
Buchholz
,
RR
,
Conley
,
A
,
Gettelman
,
A
,
Garcia
,
R
,
Simpson
,
I
,
Blake
,
DR
,
Meinardi
,
S
,
Pétron
,
G
.
2020
.
The chemistry mechanism in the community earth system model version 2 (CESM2)
.
Journal of Advances in Modeling Earth Systems
12
(
4
):
e2019MS001882
. DOI: http://dx.doi.org/10.1029/2019MS001882.
Fisher
,
JA
,
Jacob
,
DJ
,
Travis
,
KR
,
Kim
,
PS
,
Marais
,
EA
,
Miller
,
CC
,
Yu
,
K
,
Zhu
,
L
,
Yantosca
,
RM
,
Sulprizio
,
MP
,
Mao
,
J
,
Wennberg
,
PO
,
Crounse
,
JD
,
Teng
,
AP
,
Nguyen
,
TB
, St
Clair
,
JM
,
Cohen
,
RC
,
Romer
,
P
,
Nault
,
BA
,
Wooldridge
,
PJ
,
Jimenez
,
JL
,
Campuzano-Jost
,
P
,
Day
,
DA
,
Hu
,
W
,
Shepson
,
PB
,
Xiong
,
F
,
Blake
,
DR
,
Goldstein
,
AH
,
Misztal
,
PK
,
Hanisco
,
TF
,
Wolfe
,
GM
,
Ryerson
,
TB
,
Wisthaler
,
A
,
Mikoviny
,
T
.
2016
.
Organic nitrate chemistry and its implications for nitrogen budgets in an isoprene- and monoterpene-rich atmosphere: Constraints from aircraft (SEAC(4)RS) and ground-based (SOAS) observations in the Southeast US
.
Atmospheric Chemistry and Physics
16
(
9
):
5969
5991
. DOI: http://dx.doi.org/10.5194/acp-16-5969-2016.
Foley
,
KM
,
Roselle
,
SJ
,
Appel
,
KW
,
Bhave
,
PV
,
Pleim
,
JE
,
Otte
,
TL
,
Mathur
,
R
,
Sarwar
,
G
,
Young
,
JO
,
Gilliam
,
RC
,
Nolte
,
CG
,
Kelly
,
JT
,
Gilliland
,
AB
,
Bash
,
JO
.
2010
.
Incremental testing of the community multiscale air quality (CMAQ) modeling system version 4.7
.
Geoscientific Model Development
3
(
1
):
205
226
. DOI: http://dx.doi.org/10.5194/gmd-3-205-2010.
Fountoukis
,
C
,
Nenes
,
A
.
2007
.
ISORROPIA II: A computationally efficient thermodynamic equilibrium model for K+–Ca2+–Mg2+–NH4+–Na+–SO42−–NO3–Cl–H2O aerosols
.
Atmospheric Chemistry and Physics
7
(
17
):
4639
4659
. DOI: http://dx.doi.org/10.5194/acp-7-4639-2007.
Fried
,
A
,
Walega
,
J
,
Weibring
,
P
,
Richter
,
D
,
Simpson
,
IJ
,
Blake
,
DR
,
Blake
,
NJ
,
Meinardi
,
S
,
Barletta
,
B
,
Hughes
,
SC
,
Crawford
,
JH
,
Diskin
,
G
,
Barrick
,
J
,
Hair
,
J
,
Fenn
,
M
,
Wisthaler
,
A
,
Mikoviny
,
T
,
Woo
,
J-H
,
Kim
,
J
,
Min
,
K-E
,
Jeong
,
S
,
Wennberg
,
PO
,
Kim
,
MJ
,
Crounse
,
JD
,
Teng
,
AP
,
Benett
,
R
,
Yang-Martin
,
M
,
Shook
,
M
,
Huey
,
G
,
Tanner
,
D
,
Knote
,
C
,
Kim
,
JH
,
Park
,
RJ
,
Brune
,
W
.
n.d.
Airborne formaldehyde and VOC measurements over the daesan petrochemical complex on Korea’s northwest coast during the KORUS-AQ study: Estimation of emission fluxes and effects on air quality
,
in press
.
Gaubert
,
B
,
Emmons
,
LK
,
Raeder
,
K
,
Tilmes
,
S
,
Miyazaki
,
K
,
Arellano
Jr,
AF
,
Elguindi
,
N
,
Granier
,
C
,
Tang
,
W
,
Barré
,
J
,
Worden
,
HM
,
Buchholz
,
RR
,
Edwards
,
DP
,
Franke
,
P
,
Anderson
,
JL
,
Saunois
,
M
,
Schroeder
,
J
,
Woo
,
JH
,
Simpson
,
IJ
,
Blake
,
DR
,
Meinardi
,
S
,
Wennberg
,
PO
,
Crounse
,
J
,
Teng
,
A
,
Kim
,
M
,
Dickerson
,
RR
,
He
,
H
,
Ren
,
X
.
2020
.
Correcting model biases of CO in East Asia: Impact on oxidant distributions during KORUS-AQ
.
Atmospheric Chemistry and Physics Discuss
1
49
. DOI: http://dx.doi.org/10.5194/acp-2020-599.
Goldberg
,
DL
,
Saide
,
PE
,
Lamsal
,
LN
, de
Foy
,
B
,
Lu
Z
,
Woo
,
J-H
,
Kim
,
Y
,
Kim
,
J
,
Gao
,
M
,
Carmichael
,
G
,
Streets
,
DG
.
2019
.
A top-down assessment using OMI NO2 suggests an underestimate in the NOxemissions inventory in Seoul, South Korea, during KORUS-AQ
.
Atmospheric Chemistry and Physics
19
(
3
):
1801
1818
. DOI: http://dx.doi.org/10.5194/acp-19-1801-2019.
Grell
,
GA
,
Peckham
,
SE
,
Schmitz
,
R
,
McKeen
,
SA
,
Frost
,
G
,
Skamarock
,
WC
,
Eder
,
B
.
2005
.
Fully coupled “online” chemistry within the WRF model
.
Atmospheric Environment
39
(
37
):
6957
6975
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2005.04.027.
Guenther
,
AB
,
Jiang
,
X
,
Heald
,
CL
,
Sakulyanontvittaya
,
T
,
Duhl
,
T
,
Emmons
,
LK
,
Wang
,
X
.
2012
.
The model of emissions of gases and aerosols from nature version 2.1 (MEGAN2.1): An extended and updated framework for modeling biogenic emissions
.
Geoscientific Model Development
5
(
6
):
1471
1492
. DOI: http://dx.doi.org/10.5194/gmd-5-1471-2012.
Hayes
,
PL
,
Carlton
,
AG
,
Baker
,
KR
,
,
R
,
Washenfelder
,
RA
,
Alvarez
,
S
,
Rappenglück
,
B
,
Gilman
,
JB
,
Kuster
,
WC
,
de Gouw
,
JA
,
Zotter
,
P
,
Prévôt
,
ASH
,
Szidat
,
S
,
Kleindienst
,
TE
,
Offenberg
,
JH
,
Ma
,
PK
,
Jimenez
,
JL
.
2015
.
Modeling the formation and aging of secondary organic aerosols in Los Angeles during CalNex 2010
.
Atmospheric Chemistry and Physics
15
(
10
):
5773
5801
. DOI: http://dx.doi.org/10.5194/acp-15-5773-2015.
Heald
,
CL
,
Jacob
,
DJ
,
Park
,
RJ
,
Russell
,
LM
,
Huebert
,
BJ
,
Seinfeld
,
JH
,
Liao
,
H
,
Weber
,
RJ
.
2005
.
A large organic aerosol source in the free troposphere missing from current models
.
Geophysical Research Letters
32
(
18
). DOI: http://dx.doi.org/10.1029/2005gl023831.
Heald
,
CL
,
Coe
,
H
,
Jimenez
,
J
,
Weber
,
R
,
Bahreini
,
R
,
Middlebrook
,
A
,
Russell
,
L
,
Jolleys
,
M
,
Fu
,
T-M
,
Allan
,
J
,
Bower
,
K
,
Capes
,
G
,
Crosier
,
J
,
Morgan
,
W
,
Robinson
,
N
,
Williams
,
P
,
Cubison
,
M
,
DeCarlo
,
P
,
Dunlea
,
E
.
2011
.
Exploring the vertical profile of atmospheric organic aerosol: Comparing 17 aircraft field campaigns with a global model
.
Atmospheric Chemistry and Physics Discuss
11
:
25371
25425
. DOI: http://dx.doi.org/10.5194/acpd-11-25371-2011.
Hodzic
,
A
,
Jimenez
,
J
.
2011
.
Modeling anthropogenically controlled secondary organic aerosols in a megacity: A simplified framework for global and climate models
.
Geoscientific Model Development
4
:
901
917
. DOI: http://dx.doi.org/10.5194/gmd-4-901-2011.
Huang
,
M
,
Crawford
,
JH
,
Diskin
,
GS
,
Santanello
,
JA
,
Kumar
,
SV
,
Pusede
,
SE
,
Parrington
,
M
,
Carmichael
,
GR
.
2018
.
Modeling regional pollution transport events during KORUS-AQ: Progress and challenges in improving representation of land-atmosphere feedbacks
.
Journal of geophysical research Atmospheres: JGR
123
(
18
):
10732
10756
. DOI: http://dx.doi.org/10.1029/2018jd028554.
Huneeus
N
,
Schulz
M
,
Balkanski
Y
,
Griesfeller
J
,
Prospero
J
,
Kinne
,
S
,
Bauer
,
S
,
Boucher
,
O
,
Chin
,
M
,
Dentener
,
F
,
Diehl
,
T
,
Easter
,
R
,
Fillmore
,
D
,
Ghan
,
S
,
Ginoux
,
P
,
Grini
,
A
,
Horowitz
,
L
,
Koch
,
D
,
Krol
,
MC
,
Landing
,
W
,
Liu
,
X
,
Mahowald
,
N
,
Miller
,
R
,
Morcrette
,
JJ
,
Myhre
,
G
,
Penner
,
J
,
Perlwitz
,
J
,
Stier
,
P
,
Takemura
,
T
,
Zender
,
CS
.
2011
.
Global dust model intercomparison in AeroCom phase I
.
Atmospheric Chemistry and Physics
11
(
15
):
7781
7816
. DOI: http://dx.doi.org/10.5194/acp-11-7781-2011.
Itahashi
,
S
,
Ge
,
B
,
Sato
,
K
,
Fu
,
JS
,
Wang
,
X
,
Yamaji
,
K
,
Nagashima
,
T
,
Li
,
J
,
Kajino
,
M
,
Liao
,
H
,
Zhang
,
M
,
Wang
,
Z
,
Li
,
M
,
Kurokawa
,
J
,
Carmichael
,
GR
,
Wang
,
Z
.
2020
.
MICS-Asia III: Overview of model intercomparison and evaluation of acid deposition over Asia
.
Atmospheric Chemistry and Physics
20
(
5
):
2667
2693
. DOI: http://dx.doi.org/10.5194/acp-20-2667-2020.
Jeong
,
JI
,
Park
,
RJ
.
2017
.
Winter monsoon variability and its impact on aerosol concentrations in East Asia
.
Environmental Pollution
221
:
285
292
. DOI: http://dx.doi.org/10.1016/j.envpol.2016.11.075.
Jeong
,
JI
,
Park
,
RJ
.
2018
.
Efficacy of dust aerosol forecasts for East Asia using the adjoint of GEOS-Chem with ground-based observations
.
Environmental Pollution
234
:
885
893
. DOI: http://dx.doi.org/10.1016/j.envpol.2017.12.025.
Jordan
,
C
,
Crawford
,
J
,
Beyersdorf
,
A
,
Eck
,
T
,
Halliday
,
H
,
Nault
,
B
,
Chang
,
L-S
,
Park
,
J
,
Park
,
RJ
,
Lee
,
G
,
Kim
,
H
,
Ahn
,
J-y
,
Cho
,
S
,
Shin
,
H
,
Lee
,
J
,
Jung
,
J
,
Kim
,
D-S
,
Lee
,
M
,
Lee
,
T
,
Schwarz
,
J
.
2020
.
Investigation of factors controlling PM2.5 variability across the South Korean Peninsula during KORUS-AQ
.
Elementa: Science of the Anthropocene
8
:
28
. DOI: http://dx.doi.org/10.1525/elementa.424.
Kim
,
H
,
Zhang
,
Q
,
Heo
,
J
.
2018
.
Influence of intense secondary aerosol formation and long-range transport on aerosol chemistry and properties in the Seoul Metropolitan Area during spring time: Results from KORUS-AQ
.
Atmospheric Chemistry and Physics
18
(
10
):
7149
7168
. DOI: http://dx.doi.org/10.5194/acp-18-7149-2018.
Kim
,
JH
,
Lee
,
H
.
2010
.
What causes the springtime tropospheric ozone maximum over Northeast Asia?
27
(
3
):
543
551
. DOI: http://dx.doi.org/10.1007/s00376-009-9098-z.
Kim
,
PS
,
Jacob
,
DJ
,
Fisher
,
JA
,
Travis
,
K
,
Yu
,
K
,
Zhu
,
L
,
Yantosca
,
RM
,
Sulprizio
,
MP
,
Jimenez
,
JL
,
Campuzano-Jost
,
P
,
Froyd
,
KD
,
Liao
,
J
,
Hair
,
JW
,
Fenn
,
MA
,
Butler
,
CF
,
Wagner
,
NL
,
Gordon
,
TD
,
Welti
,
A
,
Wennberg
,
PO
,
Crounse
,
JD
,
St. Clair
,
JM
,
Teng
,
AP
,
Millet
,
DB
,
Schwarz
,
JP
,
Markovic
,
MZ
,
Perring
,
AE
.
2015
.
Sources, seasonality, and trends of southeast US aerosol: An integrated analysis of surface, aircraft, and satellite observations with the GEOS-Chem chemical transport model
.
Atmospheric Chemistry and Physics
15
(
18
):
10411
10433
. DOI: http://dx.doi.org/10.5194/acp-15-10411-2015.
Knote
,
C
,
Hodzic
,
A
,
Jimenez
,
JL
,
Volkamer
,
R
,
Orlando
,
JJ
,
Baidar
,
S
,
Brioude
,
J
,
Fast
,
J
,
Gentner
,
DR
,
Goldstein
,
AH
,
Hayes
,
PL
,
Knighton
,
WB
,
Oetjen
,
H
,
Setyan
,
A
,
Stark
,
H
,
Thalman
,
R
,
Tyndall
,
G
,
Washenfelder
,
R
,
Waxman
,
E
,
Zhang
,
Q
.
2014
.
Simulation of semi-explicit mechanisms of SOA formation from glyoxal in aerosol in a 3-D model
.
Atmospheric Chemistry and Physics
14
(
12
):
6213
6239
. DOI: http://dx.doi.org/10.5194/acp-14-6213-2014.
KORUS-AQ data
.
May
20
,
2019
. DOI: http://dx.doi.org/10.5067/Suborbital/KORUSAQ/DATA01.
Available at
https://www-air.larc.nasa.gov/missions/korus-aq/,
merged data version 6. Accessed 10 September 2020
.
Kwon
,
H
,
Park
,
RJ
,
Oak
,
YJ
,
Nowlan
,
CR
,
Janz
,
SJ
,
Kowalewski
,
MG
,
Fried
,
A
,
Walega
,
J
,
Bates
,
KH
,
Choi
,
J
,
Blake
,
DR
,
Wisthaler
,
A
,
Woo
J-H
.
n.d.
Top-down estimates of anthropogenic VOC emissions in South Korea using formaldehyde vertical column densities from aircraft platforms during the KORUS-AQ campaign
.
Elementa: Science of the Anthropocene
,
in press
.
Lamb
,
KD
,
Perring
,
AE
,
Samset
,
B
,
Peterson
,
D
,
Davis
,
S
,
Anderson
,
BE
,
Beyersdorf
,
A
,
Blake
,
DR
,
Campuzano-Jost
,
P
,
Corr
,
CA
,
Diskin
,
GS
,
Kondo
,
Y
,
Moteki
,
N
,
Nault
,
BA
,
Oh
,
J
,
Park
,
M
,
Pusede
,
SE
,
Simpson
,
IJ
,
Thornhill
,
KL
,
Wisthaler
,
A
,
Schwarz
,
JP
.
2018
.
Estimating source region influences on black carbon abundance, microphysics, and radiative effect observed over South Korea
.
Journal of Geophysical Research: Atmospheres
123
(
23
):
13,527–13,548
. DOI: http://dx.doi.org/10.1029/2018JD029257.
Lee
,
K
,
Yu
,
J
,
Lee
,
S
,
Park
,
M
,
Hong
,
H
,
Park
,
S
,
Choi
,
M
,
Kim
,
J
,
Kim
,
Y
,
Woo
,
J-H
,
Kim
,
S-W
,
Song
,
C
.
2020
.
Development of Korean air quality prediction system version 1 (KAQPS v1) with focuses on practical issues
.
Geoscientific Model Development
13
:
1055
1073
. DOI: http://dx.doi.org/10.5194/gmd-13-1055-2020.
Li
,
J
,
Nagashima
,
T
,
Kong
,
L
,
Ge
,
BZ
,
Yamaji
,
K
,
Fu
,
JS
,
Wang
,
XM
,
Fan
,
Q
,
Itahashi
,
S
,
Lee
,
HJ
,
Kim
,
CH
,
Lin
,
CY
,
Zhang
,
MG
,
Tao
,
ZN
,
Kajino
,
M
,
Liao
,
H
,
Li
,
M
,
Woo
,
JH
,
Kurokawa
,
J
,
Wang
,
Z
,
Wu
,
QZ
,
Akimoto
,
H
,
Carmichael
,
GR
,
Wang
,
ZF
.
2019
.
Model evaluation and intercomparison of surface-level ozone and relevant species in East Asia in the context of MICS-Asia Phase III - Part 1: Overview
.
Atmospheric Chemistry and Physics
19
(
20
):
12993
13015
. DOI: http://dx.doi.org/10.5194/acp-19-12993-2019.
Li
,
J
,
Wang
,
Z
,
Akimoto
,
H
,
Gao
,
C
,
Pochanart
,
P
,
Wang
,
X
.
2007
.
Modeling study of ozone seasonal cycle in lower troposphere over East Asia
.
Journal of Geophysical Research: Atmospheres
112
(
D22
). DOI: http://dx.doi.org/10.1029/2006JD008209.
Liu
,
X
,
Ma
,
PL
,
Wang
,
H
,
Tilmes
,
S
,
Singh
,
B
,
Easter
,
RC
,
Ghan
,
SJ
,
Rasch
,
PJ
.
2016
.
Description and evaluation of a new four-mode version of the Modal Aerosol Module (MAM4) within version 5.3 of the community atmosphere model
.
Geoscientific Model Development
9
(
2
):
505
522
. DOI: http://dx.doi.org/10.5194/gmd-9-505-2016.
Mao
,
J
,
Fan
,
S-M
,
Jacob
,
D
,
Travis
,
K
.
2012
.
Radical loss in the atmosphere from Cu-Fe redox coupling in aerosols
.
Atmospheric Chemistry & Physics Discussions
12
:
27053
27076
. DOI: http://dx.doi.org/10.5194/acpd-12-27053-2012.
Marais
,
EA
,
Jacob
,
DJ
,
Jimenez
,
JL
,
Campuzano-Jost
,
P
,
Day
,
DA
,
Hu
,
W
,
Krechmer
,
J
,
Zhu
,
L
,
Kim
,
PS
,
Miller
,
CC
,
Fisher
,
JA
,
Travis
,
K
,
Yu
,
K
,
Hanisco
,
TF
,
Wolfe
,
GM
,
Arkinson
,
HL
,
Pye
,
HOT
,
Froyd
,
KD
,
Liao
,
J
,
McNeill
,
VF
.
2016
.
Aqueous-phase mechanism for secondary organic aerosol formation from isoprene: Application to the southeast United States and co-benefit of SO2 emission controls
.
Atmospheric Chemistry and Physics
16
(
3
):
1603
1618
. DOI: http://dx.doi.org/10.5194/acp-16-1603-2016.
McNaughton
,
CS
,
Clarke
,
,
Kapustin
,
V
,
Shinozuka
,
Y
,
Howell
,
SG
,
Anderson
,
BE
,
,
E
,
Dibb
,
J
,
Scheuer
,
E
,
Cohen
,
RC
,
Wooldridge
,
P
,
Perring
,
A
,
Huey
,
LG
,
Kim
,
S
,
Jimenez
,
JL
,
Dunlea
,
EJ
,
DeCarlo
,
PF
,
Wennberg
,
PO
,
Crounse
,
JD
,
Weinheimer
,
AJ
,
Flocke
,
F
.
2009
.
Observations of heterogeneous reactions between Asian pollution and mineral dust over the Eastern North Pacific during INTEX-B
.
Atmospheric Chemistry and Physics
9
(
21
):
8283
8308
. DOI: http://dx.doi.org/10.5194/acp-9-8283-2009.
Miyazaki
,
K
,
Sekiya
,
T
,
Fu
,
D
,
Bowman
,
KW
,
Kulawik
,
SS
,
Sudo
,
K
,
Walker
,
T
,
Kanaya
,
Y
,
Takigawa
,
M
,
Ogochi
,
K
,
Eskes
,
H
,
Boersma
,
KF
,
Thompson
,
AM
,
Gaubert
,
B
,
Barre
,
J
,
Emmons
,
LK
.
2019
.
Balance of emission and dynamical controls on ozone during the Korea-United States air quality campaign from multiconstituent satellite data assimilation
.
Journal of Geophysical Research: Atmospheres
124
(
1
):
387
413
. DOI: http://dx.doi.org/10.1029/2018JD028912.
Müller
,
M
,
Mikoviny
,
T
,
Feil
,
S
,
Haidacher
,
S
,
Hanel
,
G
,
Hartungen
,
E
,
Jordan
,
A
,
Märk
,
L
,
Mutschlechner
,
P
,
Schottkowsky
,
R
,
Sulzer
,
P
,
Crawford
,
JH
,
Wisthaler
,
A
.
2014
.
A compact PTR-ToF-MS instrument for airborne measurements of volatile organic compounds at high spatiotemporal resolution
.
Atmospheric Measurement Techniques
7
(
11
):
3763
3772
. DOI: http://dx.doi.org/10.5194/amt-7-3763-2014.
National Institute of Environmental Research and National Aeronautics and Space Administration
.
2017
.
KORUS-AQ Rapid Science Synthesis Report
.
Available at
Accessed 10 September, 2020
.
Nault
,
BA
,
Campuzano-Jost
,
P
,
Day
,
DA
,
Schroder
,
JC
,
Anderson
,
B
,
Beyersdorf
,
AJ
,
Blake
,
DR
,
Brune
,
WH
,
Choi
,
Y
,
Corr
,
CA
,
de Gouw
,
JA
,
Dibb
,
J
,
DiGangi
,
JP
,
Diskin
,
GS
,
Fried
,
A
,
Huey
,
LG
,
Kim
,
MJ
,
Knote
,
CJ
,
Lamb
,
KD
,
Lee
,
T
,
Park
,
T
,
Pusede
,
SE
,
Scheuer
,
E
,
Thornhill
,
KL
,
Woo
,
JH
,
Jimenez
,
JL
.
2018
.
Secondary organic aerosol production from local emissions dominates the organic aerosol budget over Seoul, South Korea, during KORUS-AQ
.
Atmospheric Chemistry and Physics
18
(
24
):
17769
17800
. DOI: http://dx.doi.org/10.5194/acp-18-17769-2018.
Nenes
,
A
,
Pandis
,
SN
,
Pilinis
,
C
.
1998
.
ISORROPIA: A new thermodynamic equilibrium model for multiphase multicomponent inorganic aerosols
.
Aquatic Geochemistry
4
(
1
):
123
152
. DOI: http://dx.doi.org/10.1023/A:1009604003981.
Oak
,
YJ
,
Park
,
RJ
,
Schroeder
,
J
,
Crawford
,
J
,
Blake
,
D
,
Weinheimer
,
A
,
Woo
,
J-H
,
Kim
,
S-W
,
Yeo
,
H
,
Fried
,
A
,
Wisthaler
,
A
,
Brune
,
W
.
2019
.
Evaluation of simulated O3 production efficiency during the KORUS-AQ campaign: Implications for anthropogenic NOx emissions in Korea
.
Elementa: Science of the Anthropocene
7
:
56
. DOI: http://dx.doi.org/10.1525/elementa.394.
Odum
,
JR
,
Hoffmann
,
T
,
Bowman
,
F
,
Collins
,
D
,
Flagan
,
RC
,
Seinfeld
,
JH
.
1996
.
Gas/particle partitioning and secondary organic aerosol yields
.
Environmental Science & Technology
30
(
8
):
2580
2585
. DOI: http://dx.doi.org/10.1021/es950943+.
Paulot
,
F
,
Crounse
,
JD
,
Kjaergaard
,
HG
,
Kürten
,
A
,
St. Clair
,
JM
,
Seinfeld
,
JH
,
Wennberg
,
PO
.
2009
.
Unexpected epoxide formation in the gas-phase photooxidation of isoprene
.
Science
325
(
5941
):
730
. DOI: http://dx.doi.org/10.1126/science.1172910.
Peterson
,
D
,
Hyer
,
E
,
Han
,
S-O
,
Crawford
,
J
,
Park
,
R
,
Holz
,
R
,
Kuehn
,
R
,
Eloranta
,
E
,
Knote
,
C
,
Jordan
,
C
,
Lefer
,
B
.
2019
.
Meteorology influencing springtime air quality, pollution transport, and visibility in Korea
.
Elementa: Science of the Anthropocene
7
:
57
. DOI: http://dx.doi.org/10.1525/elementa.395.
Petzold
,
A
,
Schönlinner
,
M
.
2004
.
Multi-angle absorption photometry—A new method for the measurement of aerosol light absorption and atmospheric black carbon
.
Journal of Aerosol Science
35
(
4
):
421
441
. DOI: http://dx.doi.org/10.1016/j.jaerosci.2003.09.005.
Pfister
,
GG
,
Walters
,
S
,
Lamarque
,
J-F
,
Fast
,
J
,
Barth
,
MC
,
Wong
,
J
,
Done
,
J
,
Holland
,
G
,
Bruyère
,
CL
.
2014
.
Projections of future summertime ozone over the U.S
.
Journal of Geophysical Research: Atmospheres
119
(
9
):
5559
5582
. DOI: http://dx.doi.org/10.1002/2013jd020932.
Porter
,
WC
,
Safieddine
,
SA
,
Heald
,
CL
.
2017
.
Impact of aromatics and monoterpenes on simulated tropospheric ozone and total OH reactivity
.
Atmospheric Environment
169
:
250
257
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2017.08.048.
Richter
,
D
,
Weibring
,
P
,
Walega
,
JG
,
Fried
,
A
,
Spuler
,
SM
,
Taubman
,
MS
.
2015
.
Compact highly sensitive multi-species airborne mid-IR spectrometer
.
Applied Physics B
119
(
1
):
119
131
. DOI: http://dx.doi.org/10.1007/s00340-015-6038-8.
Ridley
,
BA
,
Grahek
,
FE
.
1990
.
A small, low flow, high sensitivity reaction vessel for NO chemiluminescence detectors
.
Journal of Atmospheric and Oceanic Technology
7
(
2
):
307
311
. DOI: http://dx.doi.org/10.1175/1520-0426(1990)007<0307:ASLFHS>2.0.CO;2.
Sachse
,
GW
,
Collins
,
JEJ
,
Hill
,
GF
,
,
LO
,
Burney
,
LG
,
Ritter
,
JA
.
1991
.
Airborne tunable diode laser sensor for high-precision concentration and flux measurements of carbon monoxide and methane
.
Proceedings SPIE 1433, Measurement of Atmospheric Gases
. DOI: http://dx.doi.org/10.1117/12.46162.
Saide
,
PE
,
Gao
,
M
,
Lu
,
Z
,
Goldberg
,
DL
,
Streets
,
DG
,
Woo
,
JH
,
Beyersdorf
,
A
,
Corr
,
CA
,
Thornhill
,
KL
,
Anderson
,
B
,
Hair
,
JW
,
Nehrir
,
AR
,
Diskin
,
GS
,
Jimenez
,
JL
,
Nault
,
BA
,
Campuzano-Jost
,
P
,
Dibb
,
J
,
Heim
,
E
,
Lamb
,
KD
,
Schwarz
,
JP
,
Perring
,
AE
,
Kim
,
J
,
Choi
,
M
,
Holben
,
B
,
Pfister
,
G
,
Hodzic
,
A
,
Carmichael
,
GR
,
Emmons
,
L
,
Crawford
,
JH
.
2020
.
Understanding and improving model representation of aerosol optical properties for a Chinese haze event measured during KORUS-AQ
.
Atmospheric Chemistry and Physics
20
(
11
):
6455
6478
. DOI: http://dx.doi.org/10.5194/acp-20-6455-2020.
Schroeder
,
J
,
Crawford
,
J
,
Ahn
,
J-Y
,
Chang
,
L
,
Fried
,
A
,
Walega
,
J
,
Weinheimer
,
A
,
Montzka
,
D
,
Hall
,
S
,
Ullmann
,
K
,
Wisthaler
,
A
,
Mikoviny
,
T
,
Chen
,
G
,
Blake
,
D
,
Blake
,
N
,
Hughes
,
S
,
Meinardi
,
S
,
Diskin
,
G
,
Digangi
,
J
,
Wennberg
,
P
.
2020
.
Observation-based modeling of ozone chemistry in the Seoul metropolitan area during the Korea-United States Air Quality Study (KORUS-AQ)
.
Elementa: Science of the Anthropocene
8
:
3
. DOI: http://dx.doi.org/10.1525/elementa.400.
Schwarz
,
JP
,
Samset
,
BH
,
Perring
,
AE
,
Spackman
,
JR
,
Gao
,
RS
,
Stier
,
P
,
Schulz
,
M
,
Moore
,
FL
,
Ray
,
EA
,
Fahey
,
DW
.
2013
.
Global-scale seasonally resolved black carbon vertical profiles over the Pacific
.
Geophysical Research Letters
40
(
20
):
5542
5547
. DOI: http://dx.doi.org/10.1002/2013GL057775.
Shrivastava
,
M
,
Cappa
,
CD
,
Fan
,
J
,
Goldstein
,
AH
,
Guenther
,
AB
,
Jimenez
,
JL
,
Kuang
,
C
,
,
A
,
Martin
,
ST
,
Ng
,
NL
,
Petaja
,
T
,
Pierce
,
JR
,
Rasch
,
PJ
,
Roldin
,
P
,
Seinfeld
,
JH
,
Shilling
,
J
,
Smith
,
JN
,
Thornton
,
JA
,
Volkamer
,
R
,
Wang
,
J
,
Worsnop
,
DR
,
Zaveri
,
RA
,
Zelenyuk
,
A
,
Zhang
,
Q
.
2017
.
Recent advances in understanding secondary organic aerosol: Implications for global climate forcing
.
Reviews of Geophysics
55
(
2
):
509
559
. DOI: http://dx.doi.org/10.1002/2016RG000540.
Shrivastava
,
M
,
Fast
,
J
,
Easter
,
R
,
Gustafson
Jr,
WI
,
Zaveri
,
RA
,
Jimenez
,
JL
,
Saide
,
P
,
Hodzic
,
A
.
2011
.
Modeling organic aerosols in a megacity: Comparison of simple and complex representations of the volatility basis set approach
.
Atmospheric Chemistry and Physics
11
(
13
):
6639
6662
. DOI: http://dx.doi.org/10.5194/acp-11-6639-2011.
Simpson
,
IJ
,
Blake
,
DR
,
Blake
,
NJ
,
Meinardi
,
S
,
Barletta
,
B
,
Hughes
,
SC
,
Fleming
,
LT
,
Crawford
,
JH
,
Diskin
,
GS
,
Emmons
,
LK
,
Fried
,
A
,
Guo
,
H
,
Peterson
,
DA
,
Wisthaler
,
A
,
Woo
,
J-H
,
Barré
,
J
,
Gaubert
,
B
,
Kim
,
J
,
Kim
,
MJ
,
Kim
,
Y
,
Knote
,
C
,
Mikoviny
,
T
,
Pusede
,
SE
,
Schroeder
,
JR
,
Wang
,
Y
,
Wennberg
,
PO
,
Zeng
,
L
.
2020
.
Characterization, sources and reactivity of volatile organic compounds (VOCs) in Seoul and surrounding regions during KORUS-AQ
.
Elementa: Science of the Anthropocene
8
:
37
. DOI: http://dx.doi.org/10.1525/elementa.434.
Slusher
,
DL
,
Huey
,
LG
,
Tanner
,
DJ
,
Flocke
,
FM
,
Roberts
,
JM
.
2004
.
A thermal dissociation–chemical ionization mass spectrometry (TD-CIMS) technique for the simultaneous measurement of peroxyacyl nitrates and dinitrogen pentoxide
.
Journal of Geophysical Research: Atmospheres
109
(
D19
). DOI: http://dx.doi.org/10.1029/2004JD004670.
Souri
,
AH
,
Nowlan
,
CR
,
Wolfe
,
GM
,
Lamsal
,
LN
,
Chan Miller
,
CE
,
,
GG
,
Janz
,
SJ
,
Fried
,
A
,
Blake
,
DR
,
Weinheimer
,
AJ
,
Diskin
,
GS
,
Liu
,
X
,
Chance
,
K
.
2020
.
Revisiting the effectiveness of HCHO/NO2 ratios for inferring ozone sensitivity to its precursors using high resolution airborne remote sensing observations in a high ozone episode during the KORUS-AQ campaign
.
Atmospheric Environment
224
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2020.117341.
St. Clair
,
JM
,
McCabe
,
DC
,
Crounse
,
JD
,
Steiner
,
U
,
Wennberg
,
PO
.
2010
.
Chemical ionization tandem mass spectrometer for the in situ measurement of methyl hydrogen peroxide
.
Review of Scientific Instruments
81
(
9
):
094102
. DOI: http://dx.doi.org/10.1063/1.3480552.
Stanier
,
CO
,
Donahue
,
N
,
Pandis
,
SN
.
2008
.
Parameterization of secondary organic aerosol mass fractions from smog chamber data
.
Atmospheric Environment
42
(
10
):
2276
2299
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2007.12.042.
,
R
,
Lurmann
,
F
,
Pandis
,
SN
.
1999
.
Evaluation of secondary organic aerosol formation in winter
.
Atmospheric Environment
33
(
29
):
4849
4863
. DOI: http://dx.doi.org/10.1016/S1352-2310(99)00310-6.
Su
,
W
,
Liu
,
C
,
Hu
,
Q
,
Fan
,
G
,
Xie
,
Z
,
Huang
,
X
,
Zhang
,
T
,
Chen
,
Z
,
Dong
,
Y
,
Ji
,
X
,
Liu
,
H
,
Wang
,
Z
,
Liu
,
J
.
2017
.
Characterization of ozone in the lower troposphere during the 2016 G20 conference in Hangzhou
.
Scientific Reports
7
(
1
):
17368
. DOI: http://dx.doi.org/10.1038/s41598-017-17646-x.
Tan
,
J
,
Fu
,
JS
,
Carmichael
,
GR
,
Itahashi
,
S
,
Tao
,
Z
,
Huang
,
K
,
Dong
,
X
,
Yamaji
,
K
,
Nagashima
,
T
,
Wang
,
X
,
Liu
,
Y
,
Lee
,
HJ
,
Lin
,
CY
,
Ge
,
B
,
Kajino
,
M
,
Zhu
,
J
,
Zhang
,
M
,
Liao
,
H
,
Wang
,
Z
.
2020
.
Why do models perform differently on particulate matter over East Asia? A multi-model intercomparison study for MICS-Asia III
.
Atmospheric Chemistry and Physics
20
(
12
):
7393
7410
. DOI: http://dx.doi.org/10.5194/acp-20-7393-2020.
Tang
,
Y
,
Carmichael
,
GR
,
Seinfeld
,
JH
,
Dabdub
,
D
,
Weber
,
RJ
,
Huebert
,
B
,
Clarke
,
,
Guazzotti
,
SA
,
Sodeman
,
DA
,
Prather
,
KA
,
Uno
,
I
,
Woo
,
J-H
,
Yienger
,
JJ
,
Streets
,
DG
,
Quinn
,
PK
,
Johnson
,
JE
,
Song
,
C-H
,
Grassian
,
VH
,
Sandu
,
A
,
Talbot
,
RW
,
Dibb
,
JE
.
2004
.
Three-dimensional simulations of inorganic aerosol distributions in east Asia during spring 2001
.
Journal of Geophysical Research: Atmospheres
109
(
D19
). DOI: http://dx.doi.org/10.1029/2003jd004201.
Tang
,
Y
,
Carmichael
,
GR
,
Thongboonchoo
,
N
,
Chai
,
T
,
Horowitz
,
LW
,
Pierce
,
RB
,
,
JA
,
Pfister
,
G
,
Vukovich
,
JM
,
Avery
,
MA
,
Sachse
,
GW
,
Ryerson
,
TB
,
Holloway
,
JS
,
Atlas
,
EL
,
Flocke
,
FM
,
Weber
,
RJ
,
Huey
,
LG
,
Dibb
,
JE
,
Streets
,
DG
,
Brune
,
WH
.
2007
.
Influence of lateral and top boundary conditions on regional air quality prediction: A multiscale study coupling regional and global chemical transport models
.
Journal of Geophysical Research: Atmospheres
112
(
D10
). DOI: http://dx.doi.org/10.1029/2006jd007515.
Tang
,
W
,
Emmons
,
LK
,
Arellano
Jr,
AF
,
Gaubert
,
B
,
Knote
,
C
,
Tilmes
,
S
,
Buchholz
,
RR
,
Pfister
,
GG
,
Diskin
,
GS
,
Blake
,
DR
,
Blake
,
NJ
,
Meinardi
,
S
,
DiGangi
,
JP
,
Choi
,
Y
,
Woo
,
JH
,
He
,
C
,
Schroeder
,
JR
,
Suh
,
I
,
Lee
,
HJ
,
Jo
,
HY
,
Kanaya
,
Y
,
Jung
,
J
,
Lee
,
Y
,
Kim
,
D
.
2019
.
Source contributions to carbon monoxide concentrations during KORUS-AQ based on CAM-chem model applications
.
Journal of Geophysical Research: Atmospheres
124
(
5
):
2796
2822
. DOI: http://dx.doi.org/10.1029/2018jd029151.
Tao
,
Z
,
Chin
,
M
,
Gao
,
M
,
Kucsera
,
T
,
Kim
,
D
,
Bian
,
H
,
Kurokawa
,
J
,
Wang
,
Y
,
Liu
,
Z
,
Carmichael
,
GR
,
Wang
,
Z
,
Akimoto
,
H
.
2020
.
Evaluation of NU-WRF model performance on air quality simulation under various model resolutions – An investigation within the framework of MICS-Asia Phase III
.
Atmospheric Chemistry and Physics
20
(
4
):
2319
2339
. DOI: http://dx.doi.org/10.5194/acp-20-2319-2020.
Tilmes
,
S
,
Hodzic
,
A
,
Emmons
,
LK
,
Mills
,
MJ
,
Gettelman
,
A
,
Kinnison
,
DE
,
Park
,
M
,
Lamarque
,
J-F
,
Vitt
,
F
,
Shrivastava
,
M
,
Campuzano-Jost
,
P
,
Jimenez
,
JL
,
Liu
,
X
.
2019
.
Climate forcing and trends of organic aerosols in the community earth system model (CESM2)
.
Journal of Advances in Modeling Earth Systems
11
(
12
):
4323
4351
. DOI: http://dx.doi.org/10.1029/2019ms001827.
Wang
,
Y
,
Wang
,
H
,
Wang
,
W
.
2020
.
A stratospheric intrusion-influenced ozone pollution episode associated with an intense horizontal-trough event
.
Atmosphere
11
:
164
. DOI: http://dx.doi.org/10.3390/atmos11020164.
Weinheimer
,
AJ
,
Walega
,
JG
,
Ridley
,
BA
,
Gary
,
BL
,
Blake
,
DR
,
Blake
,
NJ
,
Rowland
,
FS
,
Sachse
,
GW
,
Anderson
,
BE
,
Collins
,
JE
.
1994
.
Meridional distributions of NOx, NOy, and other species in the lower stratosphere and upper troposphere during AASE II
.
Geophysical Research Letters
21
(
23
):
2583
2586
. DOI: http://dx.doi.org/10.1029/94GL01897.
Whitten
,
GZ
,
Heo
,
G
,
Kimura
,
Y
,
McDonald-Buller
,
E
,
Allen
,
DT
,
Carter
,
WPL
,
Yarwood
,
G
.
2010
.
A new condensed toluene mechanism for Carbon Bond: CB05-TU
.
Atmospheric Environment
44
(
40
):
5346
5355
. DOI: http://dx.doi.org/10.1016/j.atmosenv.2009.12.029.
Woo
,
J-H
,
Kim
,
Y
,
Kim
,
J
.
n.d. in press
.
Wooldridge
,
PJ
,
Perring
,
AE
,
Bertram
,
TH
,
Flocke
,
FM
,
Roberts
,
JM
,
Singh
,
HB
,
Huey
,
LG
,
Thornton
,
JA
,
Wolfe
,
GM
,
Murphy
,
JG
,
Fry
,
JL
,
Rollins
,
AW
,
LaFranchi
,
BW
,
Cohen
,
RC
.
2010
.
Total Peroxy Nitrates (ΣPNs) in the atmosphere: The Thermal Dissociation-Laser Induced Fluorescence (TD-LIF) technique and comparisons to speciated PAN measurements
.
Atmospheric Measurement Techniques
3
(
3
):
593
607
. DOI: http://dx.doi.org/10.5194/amt-3-593-2010.
Zaveri
,
RA
,
Easter
,
RC
,
Fast
,
JD
,
Peters
,
LK
.
2008
.
Model for simulating aerosol interactions and chemistry (MOSAIC)
.
Journal of Geophysical Research: Atmospheres
113
(
D13
). DOI: http://dx.doi.org/10.1029/2007JD008782.
Zellweger
,
C
,
Forrer
,
J
,
Hofer
,
P
,
Nyeki
,
S
,
Schwarzenbach
,
B
,
Weingartner
,
E
,
Ammann
,
M
,
Baltensperger
,
U
.
2003
.
Partitioning of reactive nitrogen (NOy) and dependence on meteorological conditions in the lower free troposphere
.
Atmospheric Chemistry and Physics
3
(
3
):
779
796
. DOI: http://dx.doi.org/10.5194/acp-3-779-2003.
Zhang
,
Q
,
Streets
,
DG
,
Carmichael
,
GR
,
He
,
KB
,
Huo
,
H
,
Kannari
,
A
,
Klimont
,
Z
,
Park
,
IS
,
Reddy
,
S
,
Fu
,
JS
,
Chen
,
D
,
Duan
,
L
,
Lei
,
Y
,
Wang
,
LT
,
Yao
,
ZL
.
2009
.
Asian emissions in 2006 for the NASA INTEX-B mission
.
Atmospheric Chemistry and Physics
9
(
14
):
5131
5153
. DOI: http://dx.doi.org/10.5194/acp-9-5131-2009.

How to cite this article: Park, RJ, Oak, YJ, Emmons, LK, Kim, CH, Pfister, GG, Carmichael, GR, Saide, PE, Cho, SY, Kim, S, Woo, JH, Crawford, JH, Gaubert, B, Lee, HJ, Park, SY, Jo, YJ, Gao, M, Tang, B, Stanier, CO, Shin, SS, Park, HY, Bae, C, Kim, E. 2021. Multi-model intercomparisons of air quality simulations for the KORUS-AQ campaign. Elementa Science of the Anthropocene 9(1). DOI: https://doi.org/10.1525/elementa.2021.00139.

Domain Editor-in-Chief: Detlev Helmig, Boulder A.I.R. LLC, Boulder, CO, USA

Associate Editor: Alex Guenther, Atmospheric Chemistry Division, National Center for Atmospheric Research, Boulder, CO, USA

Knowledge Domain: Atmospheric Science

Part of an Elementa Special Feature: KORUS-AQ

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. See http://creativecommons.org/licenses/by/4.0/.