Ozone Pollution in the North China Plain during the 2016 Air Chemistry Research in Asia (ARIAs) Campaign: Observations and a Modeling Study

: To study air pollution in the North China Plain (NCP), the Air Chemistry Research in Asia (ARIAs) campaign conducted airborne measurements of air pollutants in spring 2016. High pollutant concentrations, with O 3 > 100 ppbv, CO > 500 ppbv, and NO 2 > 10 ppbv, were observed. CMAQ simulations with the 2010 EDGAR emissions capture the spatial and temporal variations in ozone and its major precursors such as NO 2 and VOCs, with significant underestimation. Differences between CMAQ simulations and satellite observations reflect changes in anthropogenic emissions, decreased NO x emissions in megacities such as Beijing, but slight increases in other cities and rural areas. CMAQ also underestimates HCHO and CO, suggesting adjustments of the 2010 EDGAR emissions are necessary. HCHO/NO 2 column ratios derived from OMI measurements and CMAQ simulations show that VOC-sensitive chemistry dominates the ozone photochemical production in eastern China, suggesting the importance of tightening regulations on anthropogenic VOC emissions. After adjusting emissions based on satellite observations, better model performance was achieved. Because of the VOC-sensitive environment in ozone chemistry over the NCP, the underestimation of anthropogenic emissions could be important for CMAQ simulations, while future study and regulations should focus on VOC emissions with continuous controls on NO x emissions in China.

Elevated ozone concentrations have adverse impacts on both human health [29][30][31] and the ecosystem [32][33][34].Tropospheric ozone absorbs thermal radiation and acts as the third most important anthropogenic contribution to radiative forcing of climate [35][36][37].In the lower troposphere, the photolysis of ozone is an important source of atmospheric hydroxyl (OH) radicals that control the lifetimes of atmospheric species such as CO and volatile organic compounds (VOCs) [38][39][40].Tropospheric ozone has a relatively long lifetime of several days to weeks [41,42], leading to significant long-range transport of ozone and its precursors [43][44][45].Thus, investigation of ozone pollution in China is essential to support national and international policy decisions for air quality and the climate.
Tropospheric ozone is produced through complex photochemical reactions of precursors including nitrogen oxides (NO x = NO + NO 2 ) and VOCs in the presence of sunlight [46][47][48][49].In China, the sectors of power generation, industry, and transportation are the dominant sources of NO x emissions [50][51][52].Before 2010, NO x emissions in China increased substantially [53,54].Analysis of satellite data revealed that recently, NO x emissions have started decreasing in highly developed regions such as the Pearl River Delta (PRD), but still increased in other regions [55][56][57].Anthropogenic VOC emissions had a similar increasing trend in the past decades [58][59][60][61] and are projected to increase in the future [62].Therefore, the recent increase in tropospheric ozone in China could likely be explained by the enhanced anthropogenic emissions of ozone precursors.
Due to the complex O 3 -NO x -VOCs chemistry, we need to understand the photochemical regime for local ozone production, i.e., whether it is NO x -sensitive or VOCsensitive [63,64].Duncan et al. [65] used the ratio of tropospheric columns of formaldehyde (HCHO) and nitrogen dioxide (NO 2 ) observed by the National Aeronautics and Space Administration (NASA) Ozone Monitoring Instrument (OMI) onboard the Aura satellite to characterize ozone sensitivity.Studies show that a NO x -sensitive regime dominates in the United States, except in megacities such as Los Angles and New York City, where the local ozone production is in VOC-sensitive or transition regimes [65][66][67].However, studies reported that VOC-sensitive and transition regimes for ozone photochemical production exist ubiquitously in China due to the large amount of NO x emissions, especially over the NCP [67][68][69][70].As such, although the current regulations in China focus only on reduction of NO x emissions [71][72][73], air quality can also benefit from VOC controls [74,75].
Aircraft measurements are essential to study the precursor emissions, photochemical production, and transport of ozone pollution at regional scale, while airborne campaigns are sparse in China [76][77][78][79][80].To better understand the characteristics of ozone pollution, the Air Chemistry Research in Asia (ARIAs) aircraft campaign was conducted in Hebei Province in the NCP during May-June 2016, and was affiliated with the Aerosol Atmosphere Boundary-Layer Cloud (A 2 BC) experiment [81][82][83][84].Concentrations of major air pollutants in the lower atmosphere were measured during 11 research flights over the NCP, which were conducted in association with NASA's Korea U.S. Air Quality Study (KORUS-AQ) campaign in downwind South Korea.Measurements collected by the ARIAs research flights and the A 2 BC surface observations can provide a comprehensive dataset to thoroughly study the tropospheric ozone pollution and emissions of its precursors in China.
In this study, we evaluated anthropogenic emissions and ozone pollution in the NCP using a combination of aircraft measurements, satellite observations, and modeling results.The U.S. Environmental Protection Agency (EPA) Community Multiscale Air Quality (CMAQ) model was used to simulate the atmospheric chemistry for the ARIAs campaign.We evaluated the emission data through comparison with aircraft measurements and satellite products, and adjusted emissions to improve the CMAQ performance.Lastly, we investigated the sensitivity of ozone production derived from CMAQ simulations and OMI observations and discussed the future ozone pollution in China.

Aircraft Campaign in the NCP
With more than 250 million tons of iron and steel produced in 2016 (data from http://data.stats.gov.cn,accessed in 10 September 2018), Hebei Province in the NCP is the most industrialized area in China.Due to its high emissions and proximity to megacities Beijing and Tianjin, the Beijing-Tianjin-Hebei area had severe air pollution in the past decade [85,86].In May and June 2016, the ARIAs aircraft campaign was conducted over Hebei Province to investigate emissions, chemical evolution, and transport of air pollutants.The airborne campaign was coordinated with the A 2 BC field campaign in Xingtai (XT, 37.18 • N, 114.36 • E, 182 m above sea level, ASL) and the NASA KORUS-AQ campaign to expand the study to East Asia [81,87].A Harbin Y12 research airplane (similar to the de Havilland Twin Otter) was employed to measure concentrations of air pollutants including O 3 , carbon monoxide (CO), CO 2 , and NO 2 , and aerosol optical properties [81,88,89].The research airplane was based in Luancheng airport (hereafter referred to as LC, 37.91 • N, 114.59 • E, 58 m ASL), south of Shijiazhuang, the capital city of Hebei province with 10 million population.Eleven research flights were conducted during the ARIAs campaign (Figure S1 in the Supplementary Materials).Vertical profiles of air pollutants from near-surface (~100 m above ground level, AGL) to the free troposphere (>3000 m) were conducted over LC, XT (the supersite of the A 2 BC campaign), Julu (JL, 37.22 The airborne measurements of ozone were conducted using a commercially available analyzer (Model 49C, Thermo Environmental Instruments, TEI, Franklin, MA, USA) [90].NO 2 was measured using a modified commercially available cavity ring-down spectroscopy (CRDS) detector [91,92].Nitrogen oxide (NO) and reactive nitrogen compound (NO y ) concentrations were analyzed using a commercially available NO analyzer (Model 42C, Thermo Environmental Instruments) with a hot molybdenum convertor working at 375 • C [93,94].Ambient gas input was switched with and without the converter frequently to measure NO and NO y simultaneously.However, due to high power demand of the instrument and converter, NO and NO y were only measured during some research flights.Concentrations of CO and CO 2 were monitored with a 4-channel Picarro CRDS instrument (Model G2401-m, Picarro Inc., Santa Clara, CA, USA), calibrated with CO/CO 2 standards certified at the National Institute of Standards and Technology [95].All the instruments were routinely serviced, calibrated and used for airborne measurements in the United States and China [79,90,[95][96][97][98][99].Detailed information about the instrumentation including the sampling frequency, precisions and accuracies is listed in Table S1 of the Supplementary Materials.Measurements of ambient air pollutants were logged at 1 Hz frequency, but the average times for different instruments were different, as shown in Table S1.All measurements were synchronized based on the Picarro measurements of CO 2 and CO with time, geolocation, and altitude from the Global Positioning System (GPS).The delay and lag time of each instrument was considered during the post-processing of observation data and averaged to a 1-minute record for further analysis and model evaluation.
In the ARIAs research flights, 28 whole air samples (WAS) were collected in vertical spirals at different altitudes from ~400 m to ~3500 m.The WAS were analyzed using gas chromatography (GC) with Flame Ionization Detection (FID) and Mass Spectroscopy (MS) by the College of Environment Sciences and Engineering at Peking University.Seventyfour species of alkanes, alkenes/alkynes, aromatics, and halocarbons were identified and quantified for a study on ozone photochemical chemistry (see details in [88,100]).Detection limits for the compounds ranged from 2 to 60 pptv.Hourly surface O 3 , CO, and NO 2 concentrations for May and June 2016 were obtained from the public website of the China Ministry of Ecology and Environment (CMEE, available at https://quotsoft.net/air/ (accessed on 15 February 2020).The CMEE network had ~1450 monitoring sites in 2016 (Figure S2 in the Supplementary Materials).The units of CMEE measurements for O 3 /NO 2 and CO are µg/m 3 and mg/m 3 , respectively.We converted them to parts per billion by volume (ppbv) and parts per million by volume (ppmv) at standard temperature and pressure (STP) for model evaluation.

Satellite Products
To evaluate the emissions and atmospheric chemistry in the NCP and greater East Asia, we used satellite observations of CO, NO 2 , and HCHO for May and June 2016.The Measurements of Pollution In the Troposphere (MOPITT) instrument onboard the NASA Terra satellite retrieved CO column contents with ~10:30 a.m.local overpass time [101].We used the latest version 7 MOPITT Level 3 daily gridded average products (1 • × 1 • spatial resolution, available at https://eosweb.larc.nasa.gov/project/mopitt/mop03j_v007,accessed on 15 February 2020) for the ARIAs campaign period [102].MOPITT thermal-infrared and near-infrared (TIR + NIR) products show improved sensitivity to near-surface CO in China [103].We used MOPITT near-surface CO (~900 hPa) products and related averaging kernels (AKs) to evaluate the CMAQ results [104].
OMI, onboard the NASA Aura satellite, is a UV/Vis solar backscatter spectrometer in a polar sun-synchronous orbit with a ~1:35 p.m. local overpass time.With high spatial resolution (13 km × 24 km for the center at nadir) and nearly daily coverage, OMI provided monitoring of trace gases and aerosol properties from 2005 [105].The Version 3 OMI Level 2 NO 2 products (https://disc.gsfc.nasa.gov/datasets/OMNO2_V003/summary,accessed on 15 February 2020) [106] were used to evaluate the emissions and atmospheric chemistry in East Asia.Under clear sky, tropospheric NO 2 columns from OMI have precision of ~0.5 × 10 16 molecules cm −2 and an accuracy of ±30% [107].OMI HCHO Version 3 Smithsonian Astronomical Observatory (SAO) (https://disc.gsfc.nasa.gov/datasets/OMHCHO_V003/summary, accessed on 15 February 2020) Level 2 products were used in this study [108,109].The precision of column HCHO is ~1.0 × 10 16 molecules cm −2 , and SAO products have an accuracy of ±25-30% without clouds [110,111].Data in OMI pixels affected by the row anomaly and contaminated by clouds were filtered out using quality flags for both NO 2 and HCHO columns.

Model Set-Up
We used CMAQ version 5.2 [112] to simulate atmospheric chemistry for the ARIAs campaign.The Weather Research and Forecasting (WRF) model Version 3.8.1 [113] was driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim products (ds627.0,https://rda.ucar.edu/datasets/ds627.0, accessed on 15 February 2020) [114] to generate meteorological fields.Two domains with spatial resolution of 36 km and 12 km (Figure 1) were used to cover East Asia, with 35 layers from the surface to 50 hPa and ~20 layers in the lowest 2 km.Major physical options in WRF include the Rapid Radiative Transfer Model (RRTM) radiation scheme [115], the Pleim-Xiu surface layer and land surface model [116,117], the Asymmetric Convective Model (ACM2) boundary layer scheme [118], the Kain-Fritsch cumulus scheme [119], and the WRF Single-Moment 6 (WSM-6) microphysics [120].The National Centers for Environmental Prediction (NCEP) ADP Global Surface and Upper Air Observational Weather Data (ds461.0and ds351.0,https://rda.ucar.edu,accessed on 15 February 2020) were used to perform observational and analysis nudging on all domains following the method developed for NASA aircraft campaigns [97,121].WRF outputs were processed by the EPA Meteorology-Chemistry Interface Processor Version 4.3 (MCIP v4.3, released in November 2015) for emission processing and CMAQ simulations.We used the latest version 7 MOPITT Level 3 daily gridded average products (1° × 1° spatial resolution, available at https://eosweb.larc.nasa.gov/project/mopitt/mop03j_v007,accessed on 15 February 2020) for the ARIAs campaign period [102].MOPITT thermal-infrared and near-infrared (TIR + NIR) products show improved sensitivity to near-surface CO in China [103].We used MOPITT near-surface CO (~900 hPa) products and related averaging kernels (AKs) to evaluate the CMAQ results [104].OMI, onboard the NASA Aura satellite, is a UV/Vis solar backscatter spectrometer in a polar sun-synchronous orbit with a ~1:35 p.m. local overpass time.With high spatial resolution (13 km × 24 km for the center at nadir) and nearly daily coverage, OMI provided monitoring of trace gases and aerosol properties from 2005 [105].The Version 3 OMI Level 2 NO2 products (https://disc.gsfc.nasa.gov/datasets/OMNO2_V003/summary,accessed on 15 February 2020) [106] were used to evaluate the emissions and atmospheric chemistry in East Asia.Under clear sky, tropospheric NO2 columns from OMI have precision of ~0.5 × 10 16 molecules cm −2 and an accuracy of ±30% [107].OMI HCHO Version 3 Smithsonian Astronomical Observatory (SAO) (https://disc.gsfc.nasa.gov/datasets/OMHCHO_V003/summary,accessed on 15 February 2020) Level 2 products were used in this study [108,109].The precision of column HCHO is ~1.0 × 10 16 molecules cm −2 , and SAO products have an accuracy of ±25-30% without clouds [110,111].Data in OMI pixels affected by the row anomaly and contaminated by clouds were filtered out using quality flags for both NO2 and HCHO columns.

Model Set-Up
We used CMAQ version 5.2 [112] to simulate atmospheric chemistry for the ARIAs campaign.The Weather Research and Forecasting (WRF) model Version 3.8.1 [113] was driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim products (ds627.0,https://rda.ucar.edu/datasets/ds627.0, accessed on 15 February 2020) [114] to generate meteorological fields.Two domains with spatial resolution of 36 km and 12 km (Figure 1) were used to cover East Asia, with 35 layers from the surface to 50 hPa and ~20 layers in the lowest 2 km.Major physical options in WRF include the Rapid Radiative Transfer Model (RRTM) radiation scheme [115], the Pleim-Xiu surface layer and land surface model [116,117], the Asymmetric Convective Model (ACM2) boundary layer scheme [118], the Kain-Fritsch cumulus scheme [119], and the WRF Single-Moment 6 (WSM-6) microphysics [120].The National Centers for Environmental Prediction (NCEP) ADP Global Surface and Upper Air Observational Weather Data (ds461.0and ds351.0,https://rda.ucar.edu,accessed on 15 February 2020) were used to perform observational and analysis nudging on all domains following the method developed for NASA aircraft campaigns [97,121].WRF outputs were processed by the EPA Meteorology-Chemistry Interface Processor Version 4.3 (MCIP v4.3, released in November 2015) for emission processing and CMAQ simulations.Anthropogenic emissions data were from the Emissions Database for Global Atmospheric Research Version 4.2 (EDGAR v4.2, 0.1 • × 0.1 • resolution) for year 2010, which are widely used for chemical transport modeling [122].We used the EPA Sparse Matrix Operator Kernel Emissions (SMOKE) modeling system Version 4.5 [123] to project EDGAR emissions to the modeling domain.Emissions of air pollutants were speciated into Carbon Bond 05 chemical mechanism [124] and updated AERO6 aerosol module [125].The EDGAR v4.2 inventory has emissions for energy, industry, residential, and transport sectors.Without stack height information for power plants in the energy sector, we followed the approach developed in He et al. [99] to locate these anthropogenic emissions at ~200 m above the surface as an approximation for averaged stack height and plume rise.We used the United States Geological Survey (USGC) 24 category land use dataset combined with the Biogenic Emission Inventory System (BEIS) emission factors table to generate the input files for the CMAQ inline biogenic emissions modeling.Biogenic emissions were estimated using the BEIS module inline in CMAQ [112].
CMAQ v5.2 uses the updated Carbon Bond 6 (CB6r3) chemical mechanism [126] including improved chemistry mechanism for organic nitrates and peroxyacyl nitrates (PAN) chemistry and leads to better performance in simulating Secondary Organic Aerosols (SOA) and tropospheric ozone in the United States [127].CMAQ was run with a coarse domain and a nested domain (Figure 1).Chemical initial and boundary conditions for the coarse domain were obtained from the default concentration profiles built in CMAQ [112].Results from the CMAQ coarse domain were used to generate boundary conditions for the nested domain.The WRF-CMAQ system was run from mid-April to June, with the first 2 weeks as spin-up.Hourly concentrations of air pollutants were saved for further analysis and model evaluation.

Air Pollution in the NCP and CMAQ Performance
Figure 2 summarizes all aircraft measurements of O 3 , NO 2 , CO, and CO 2 over the NCP from 11 research flights.Generally, we observed high concentrations of air pollutants, with maximum values as >100 ppbv of O 3 , >20 ppbv of NO 2 , >500 ppbv of CO, and >450 ppmv of CO 2 , in the aircraft campaign area (defined as 36.5-38.5 • N, 114.0-115.5 • E hereafter).We conducted vertical spirals over XT (the A 2 BC supersite), LC (the airport in south of Shijiazhuang), and two rural areas (JL and QZ) during the ARIAs research flights.Figure 3a summarizes vertical distributions and the mean profiles of air pollutants over XT, with mean O 3 concentrations of 80 ppbv in the lower atmosphere.We observed isolated plumes with >10 ppbv of NO 2 , >1000 ppbv of CO, and >440 ppmv of CO 2 over XT, usually with a secondary maximum between 800 and 1200 m (a sample plume is presented in Figure S3 of the Supplementary Materials).These plumes aloft can play an important role in long-range transport of air pollutants to downwind regions.Profiles over LC (Figure 3b) show higher O 3 concentrations (>100 ppbv) and relatively moderate NO 2 (~3 ppbv) and CO (~250 ppbv).The rural areas, JL and QZ, have relatively clean environments with <80 ppbv of O 3 , <2 ppbv of NO 2 , and <300 ppbv of CO (Figure 3c,d).Even the concentrations of air pollutants over the rural region in the NCP are comparable to or higher than values in urban areas in North America and Europe.These results suggest both south-north and east-west gradients of air pollution in the campaign region, i.e., higher concentrations of air pollutants in the west XT-LC corridor near the mountain as compared with east side of JL-QZ, and higher concentrations in the north LC as compared with the south XT.The ARIAs research flights had good coverage of regions with both high and moderate concentrations of air pollutants and can fairly represent the regional nature of air pollution over the NCP.To evaluate the model performance in eastern China, we compared the mean concentrations of daily 8-hour maximum (MDA8) ozone in May and June 2016 from CMEE observations and CMAQ simulations driven by the EDGAR inventory (named the CMAQ_baseline run hereafter).The CMAQ_baseline run slightly underestimates concentrations of surface ozone in eastern China but cannot reproduce high surface ozone in some cities in the NCP (Figure 4a), and the overall underestimate is moderate with a linear regression slope of 0.86 (Figure 4b).For CO and NO 2 , two important ozone precursors, significant underestimations exist ubiquitously in eastern China (Figure S4 in the Supplementary Materials).These results suggest that the underestimation of ozone precursors in CMAQ could lead to poor model performance in simulating tropospheric ozone.
Similar analyses were conducted to investigate air pollutant concentrations in the lower troposphere over the NCP observed by the aircraft.A research flight on 11 June 2016 (Figure S5 in the Supplementary Materials) shows that CMAQ well captures the vertical gradient of air pollutants, while substantially underestimating concentrations of all trace gases except NO y .Since CMAQ generated hourly outputs, to alleviate the uncertainty of comparing 1-minute aircraft data and hourly model simulations, we used 10-minute averaged aircraft measurements that were matched to the closest hourly model output following the approach described in Goldberg et al. [128].Figure 5 shows similar underestimation (50% to 75% for all air pollutants) as compared with surface measurements (Figure 4 and Figure S4 in the Supplementary Materials).CMAQ overestimates NO y but substantially underestimates NO and NO 2 , which suggests that a significant amount of reactive nitrogen compounds could exist in the form of organic nitrates or nitrate aerosols in the model.Figure 6 compares total VOC concentrations from WAS samples and CMAQ simulations, indicating that VOC levels are also significantly underestimated by 80%.The model evaluation with surface and aircraft measurements suggests that ozone pollution in the NCP has been significantly underestimated in the CMAQ_baseline run, which could be due to the uncertainty introduced by using the 2010 EDGAR emissions to simulate the 2016 ARIAs campaign period.Thus, we need to evaluate the emissions inventory data to improve the CMAQ performance and investigate the sensitivity of ozone production.trations of surface ozone in eastern China but cannot reproduce high surface ozone in some cities in the NCP (Figure 4a), and the overall underestimate is moderate with a linear regression slope of 0.86 (Figure 4b).For CO and NO2, two important ozone precursors, significant underestimations exist ubiquitously in eastern China (Figure S4 in the Supplementary Materials).These results suggest that the underestimation of ozone precursors in CMAQ could lead to poor model performance in simulating tropospheric ozone.Similar analyses were conducted to investigate air pollutant concentrations in the lower troposphere over the NCP observed by the aircraft.A research flight on 11 June 2016 (Figure S5 in the Supplementary Materials) shows that CMAQ well captures the vertical gradient of air pollutants, while substantially underestimating concentrations of all trace gases except NOy.Since CMAQ generated hourly outputs, to alleviate the uncertainty of comparing 1-minute aircraft data and hourly model simulations, we used 10-minute averaged aircraft measurements that were matched to the closest hourly model output following the approach described in Goldberg et al. [128].Figure 5 shows similar underesti-

Evaluation of Emissions Inventory in the NCP
The EDGAR v4.2 emission inventory in East Asia was created based on the 2010 MIX emission inventory [129], so substantial changes were anticipated when used for the ARIAs campaign in 2016.Anthropogenic emission inventories are usually based on the "bottom-up" approach, which relies on the statistics of fossil fuel usage and emission factors (EFs) for each sector, defined as the ratio of the amount of air pollutants released by a unit of CO2 emissions, e.g., CO/CO2 and NOx/CO2.The 2010 EDGAR inventory has emissions for four sectors: Energy, Industry, Transportation, and Residential.We calculated the CO/CO2, NOx/CO2, and NOx/CO ratios through averaging the EFs from these four sectors (Figure 7).

Evaluation of Emissions Inventory in the NCP
The EDGAR v4.2 emission inventory in East Asia was created based on the 2010 MIX emission inventory [129], so substantial changes were anticipated when used for the ARIAs campaign in 2016.Anthropogenic emission inventories are usually based on the "bottom-up" approach, which relies on the statistics of fossil fuel usage and emission factors (EFs) for each sector, defined as the ratio of the amount of air pollutants released by a unit of CO 2 emissions, e.g., CO/CO 2 and NO x /CO 2 .The 2010 EDGAR inventory has emissions for four sectors: Energy, Industry, Transportation, and Residential.We calculated the CO/CO 2 , NO x /CO 2 , and NO x /CO ratios through averaging the EFs from these four sectors (Figure 7).

Evaluation of Emissions Inventory in the NCP
The EDGAR v4.2 emission inventory in East Asia was created based on the 2010 MIX emission inventory [129], so substantial changes were anticipated when used for the ARIAs campaign in 2016.Anthropogenic emission inventories are usually based on the "bottom-up" approach, which relies on the statistics of fossil fuel usage and emission factors (EFs) for each sector, defined as the ratio of the amount of air pollutants released by a unit of CO2 emissions, e.g., CO/CO2 and NOx/CO2.The 2010 EDGAR inventory has emissions for four sectors: Energy, Industry, Transportation, and Residential.We calculated the CO/CO2, NOx/CO2, and NOx/CO ratios through averaging the EFs from these four sectors (Figure 7).To evaluate the emission inventory data in the NCP, we used a plume recognition method to calculate the emission enhancements (EEs) from Y12 observations.We first selected 60 1-second aircraft measurements with a 60-second moving window.Then, we conducted linear regression of observed air pollutant (CO, NOx, etc.) concentrations vs. CO2 concentrations in each 60-second window and calculated the slope (i.e., ΔCO/ΔCO2 and ΔNOx/ΔCO2) and correlation (R).The slope is defined as EEs in each window, standing for a 'plume' tested in the 60-second window.Lastly, we only selected EEs that were within the PBL (below 1.5 km AGL in this study) and statistically significant (R 2 > 0.6), so the values of these selected EEs could act as a proxy of EFs in the air mass observed.The detailed information about this plume recognition method can be found in Halliday et al. [87].
EEs observed during the research flights have a broad range of values.ΔCO/ΔCO2 ranges from below 1%, a typical value of modern automobile emissions, to higher than To evaluate the emission inventory data in the NCP, we used a plume recognition method to calculate the emission enhancements (EEs) from Y12 observations.We first selected 60 1-second aircraft measurements with a 60-second moving window.Then, we conducted linear regression of observed air pollutant (CO, NO x , etc.) concentrations vs. CO 2 concentrations in each 60-second window and calculated the slope (i.e., ∆CO/∆CO 2 and ∆NO x /∆CO 2 ) and correlation (R).The slope is defined as EEs in each window, standing for a 'plume' tested in the 60-second window.Lastly, we only selected EEs that were within the PBL (below 1.5 km AGL in this study) and statistically significant (R 2 > 0.6), so the values of these selected EEs could act as a proxy of EFs in the air mass observed.The detailed information about this plume recognition method can be found in Halliday et al. [87].
EEs observed during the research flights have a broad range of values.∆CO/∆CO 2 ranges from below 1%, a typical value of modern automobile emissions, to higher than 10%, a value indicating fossil fuel combustion with high emissions such as biomass burning (Figure 7a).The mean of observed EE for CO (3.7%) was close to that calculated from the EDGAR inventory (4.0%) in the aircraft campaign area.Observed ∆NO x /∆CO 2 ratios also had isolated high values (>0.1%), with a mean value of 0.05%, which is substantially higher than the EF (~0.03%) derived from the EDGAR inventory (Figure 7b).Since estimation of anthropogenic CO 2 flux in an urban/suburban area is challenging [130,131], the underestimation of CO and NO x in the NCP could be caused by either underestimated EFs or uncertainty in anthropogenic CO 2 emission data used in the 'bottom-up' approach.
To further investigate the characteristics of air pollutant emissions in the NCP, we conducted a similar analysis of ∆NO x /∆CO, which are usually co-emitted in combustion processes.Since around half of the CO and NO x are from mobile sources in the EDGAR emission inventory, this ratio can approximately represent the emission characteristic of mobile sources in the NCP.The mean observed ∆NO x /∆CO ratio was ~1.3%, significantly lower than the 5.8% based on the EDGAR emission inventory (Figure 7c).These results suggest that the EDGAR emission inventory substantially overestimates the ratios of NO x /CO, while the automobile emissions over the NCP in 2016 have been greatly improved due to recent regulations, i.e., EDGAR overestimates the contribution from combustion with high emission factors.It is worth noting that we only evaluated the emission ratios (EEs or EFs) in the EDGAR inventory during the ARIAs campaign, while the underestimation of CO and NO x emissions could be caused by inaccurate CO 2 emissions, which is beyond the scope of this study.

Evaluation of CO, NO x , and VOC Emissions Using Satellite Data
Satellite observations are widely used to evaluate the anthropogenic emissions in East Asia, sometimes supplemented by model simulations, e.g., CO emissions using the MOPITT CO products [132,133], anthropogenic NO x emissions using OMI NO 2 products [134][135][136], and VOC emissions using OMI HCHO products [137].In this study, we used measurements from multiple satellite instruments to evaluate the CMAQ performance of NO 2 , HCHO, and CO.Since NO 2 and HCHO can be treated as proxies of NO x and VOC emissions, we can further improve the 2010 EDGAR emissions over the NCP based on satellite data.
We followed the approach developed in Canty et al. [138] to compare the tropospheric column contents of NO 2 from OMI products and CMAQ simulations.Level 2 OMI NO 2 swatch information including row anomaly and quality flags was used to sample NO 2 vertical profiles from CMAQ outputs, and then the CMAQ NO 2 column was calculated using the OMI averaging kernel (AK).Lastly, we averaged OMI and CMAQ NO 2 column contents to create daily 0.25 • × 0.25 • Level 3 products (see details in [138]).A similar approach was used to integrate HCHO column contents from CMAQ simulations based on OMI HCHO retrievals (see details in [66]) and construct daily 0.25 • × 0.25 • Level 3 HCHO products.For tropospheric CO, we selected the CO concentrations at ~900 hPa in CMAQ and averaged them to 1.0 • × 1.0 • daily products using the MOPITT CO averaging kernel [102].All gridded daily data of satellite and CMAQ were averaged in May and June 2016 for comparison.
We plotted OMI NO 2 columns in eastern China (Figure 8a), which showed strong signals over the NCP.NO 2 columns from OMI and CMAQ generally agreed well over eastern China, but large discrepancies existed in some regions (Figure 9a).For the aircraft campaign area, CMAQ underestimates NO 2 columns (slope = 0.95 and mean ratio = 0.81, i.e., only predicts 81% of OMI NO 2 column) with uncertainties relatively smaller within the campaign period (Figure 9b).However, in urban regions such as Beijing, the Yangtze River Delta (YRD), and the PRD, CMAQ substantially overestimates column NO 2 by up to 30% (Figure 8a).Because the CMAQ_baseline run used the 2010 anthropogenic emission data, these differences should represent the changes in NO x emissions due to recent air pollution regulations.The comparison of NO 2 columns suggests that NO x pollution of megacities in China has been substantially improved since 2010, while NO x pollution in smaller cities and rural area has worsened, consistent with results from independent studies using OMI products [57,139].OMI HCHO retrievals also show high values over the NCP in spring, when plants' photosynthetic activity is relatively low, reflecting the dominance of anthropogenic VOC emissions in northern China [61,140].CMAQ has good agreement with OMI HCHO within the aircraft campaign area (<20% underestimation), but substantially underestimates HCHO columns in southern China, where biogenic VOC emissions dominate (Figure 8b).The MOPPIT products show high near-surface CO concentrations over eastern China (Figure 8c), while the CMAQ_baseline run substantially underestimates CO over northern China and only predicts 42% of the CO over the aircraft campaign area.
smaller cities and rural area has worsened, consistent with results from independent studies using OMI products [57,139].OMI HCHO retrievals also show high values over the NCP in spring, when plants' photosynthetic activity is relatively low, reflecting the dominance of anthropogenic VOC emissions in northern China [61,140].CMAQ has good agreement with OMI HCHO within the aircraft campaign area (<20% but substantially underestimates HCHO columns in southern China, where biogenic VOC emissions dominate (Figure 8b).The MOPPIT products show high near-surface CO concentrations over eastern China (Figure 8c), while the CMAQ_baseline run substantially underestimates CO over northern China and only predicts 42% of the CO over the aircraft campaign area.Using NO2 and HCHO as proxies of NOx and VOC emissions, the comparison of satellite observations and the CMAQ_baseline simulations suggests that both NOx and VOC emissions in the aircraft campaign area need to be adjusted for a better simulation of tropospheric ozone.Also, the underestimation of CO, as an important precursor, can lead to underprediction of tropospheric ozone.We calculated the model/satellite ratios of Using NO 2 and HCHO as proxies of NO x and VOC emissions, the comparison of satellite observations and the CMAQ_baseline simulations suggests that both NO x and VOC emissions in the aircraft campaign area need to be adjusted for a better simulation of tropospheric ozone.Also, the underestimation of CO, as an important precursor, can lead to underprediction of tropospheric ozone.We calculated the model/satellite ratios of NO x , Air 2024, 2 194 HCHO, and CO in East Asia (Figure S6 in the Supplementary Materials) and used these ratios to adjust their anthropogenic emissions in CMAQ.The results will be discussed in Section 3.5.

Tropospheric Ozone Production Sensitivity from OMI and CMAQ
Photochemical production of tropospheric ozone is highly non-linear and dependent on concentrations of NO x and VOCs [64,141,142].A maximum rate of ozone production can be achieved with an optimal VOCs/NO x ratio.With other VOCs/NO x ratios, ozone production can be either in the VOC-sensitive regime (the rate of ozone production is controlled by VOC concentrations) or in the NO x -sensitive regime (the rate of production is controlled by NO x concentrations).Different pollution control strategies can be implemented to reduce the tropospheric ozone levels in these two regimes.For instance, in a VOC-sensitive environment, reducing NO x emissions will lead to limited effects until the ozone production has been changed to a NO x -sensitive environment with the continuous removal of NO x from the atmosphere.Duncan et al. [65] developed an approach using the OMI HCHO/NO 2 column ratio to estimate the ozone production sensitivity as follows: (1) HCHO/NO 2 < 1: VOC-sensitive regime; (2) HCHO/NO 2 1~2: transition regime; (3) HCHO/NO 2 > 2: NO x -sensitive regime.Studies show that urban areas in the U.S. such as Los Angeles, New York City and Houston are in VOC-sensitive or transition regimes, which lead to difficulty in local regulation of air quality [65,66,121].Recent studies suggest new threshold values of HCHO/NO 2 ratios between VOC-sensitive, transition, and NO x -sensitive regimes in the U.S. [67,143].
Using the Duncan et al. [65] approach, results from OMI products suggest large areas of eastern China are either in a VOC-sensitive regime (mostly megacities such as Beijing) or in a transition regime [67,68,144].We followed the approach described in Ring et al. [66] to calculate the column HCHO/NO 2 ratios from OMI observations and CMAQ simulations for East Asia.OMI column HCHO/NO 2 ratios suggest that the ozone photochemical production is in a VOC-sensitive or transition regime over the NCP and other major urban areas such as YRD and PRD (Figure 10a) if the Duncan et al. [65] approach is applicable for these areas.CMAQ successfully captured the spatial distribution of the regional nature of ozone production sensitivity in eastern China but predicted that the rate of ozone production is controlled more by VOCs with the CMAQ HCHO/NO 2 ratio lower than 1.0 in Beijing, YRD, and PRD (Figure 10b).The VOC-sensitive environment from both OMI observations and CMAQ simulations suggests the rate of ozone photochemical production in the NCP is controlled not only by NO x emissions, but also by VOC emissions, for which regulations are currently lacking in China.With continuous reduction of anthropogenic NO x emissions in China, VOC controls might be efficient in these VOC-sensitive regions.

Improvements in Tropospheric Ozone Simulation Using Satellite Products
Results of the previous two sections show that the CMAQ_baseline run substantially underestimates the concentrations of ozone and its major precursors in the NCP.Independent studies using KORUS-AQ observations and satellite products suggested that major ozone precursor emissions such as CO and NO x could have large discrepancies as compared with emission inventory in East Asia [145,146].To identify the individual and combined effects of emission discrepancies for major ozone precursors in the NCP, we designed a series of sensitivity experiments with emissions adjusted to satellite observations.Unlike the top-down approach using global chemical transport models such as GEOS-Chem [135,147], here we simply applied the ratios of air pollutant column contents from satellite observations and CMAQ simulations on each 0.25 degree grid (Figure S6 in the Supplementary Materials) as CO CMAQ /CO MOPITT , NO 2CMAQ /NO 2OMI , and HCHO CMAQ /HCHO OMI ratios for anthropogenic CO, NO x , and VOC emissions, respectively.To estimate the contribution from biogenic VOC emissions, we conducted one more run with the in-line BEIS module turned off.

Improvements in Tropospheric Ozone Simulation Using Satellite Products
Results of the previous two sections show that the CMAQ_baseline run substantially underestimates the concentrations of ozone and its major precursors in the NCP.Independent studies using KORUS-AQ observations and satellite products suggested that major ozone precursor emissions such as CO and NOx could have large discrepancies as compared with emission inventory in East Asia [145,146].To identify individual and combined effects of emission discrepancies for major ozone precursors in the NCP, we designed a series of sensitivity experiments with emissions adjusted to satellite observations.Unlike the top-down approach using global chemical transport models such as GEOS-Chem [135,147], here we simply applied the ratios of air pollutant column contents from satellite observations and CMAQ simulations on each 0.25 degree grid (Figure S6 in the Table 1 shows the emission adjustments for the five sensitivity experiments.CMAQ was run for the nested 12 km domain (D02) with the same meteorology, initial conditions, and boundary conditions derived from the coarse domain results of the CMAQ_baseline run.We conducted a similar analysis of Figure 8 to investigate the discrepancy between satellite products and CMAQ simulations over eastern China (Figure S7 in the Supplementary Materials).The CMAQ_all run successfully reproduced the column NO 2 and near-surface CO in eastern China and column HCHO in the NCP.In southern China, where biogenic VOCs dominate, adjusting anthropogenic VOC emissions showed limited improvements on column HCHO simulations.Figure 11 compares the mean MDA8 ozone concentrations from CMEE observations and the CMAQ_all run, and similar analyses for the other four CMAQ sensitivity experiments are presented in Figure S8 of the Supplementary Materials.CMAQ sensitivity experiments driven by modified emissions showed better results in simulating the surface ozone in eastern China, while the CMAQ_all run has successfully produced the high surface ozone in the NCP.In Table 2, we summarize the statistics from linear regression analyses of CMEE observations and all CMAQ sensitivity experiments in eastern China.Generally, the adjustment of anthropogenic emissions has improved CMAQ simulations of O 3 , CO, and NO 2 .Improvements in surface ozone simulations are sensitive to adjustments of both NO x and VOC emissions; improvements in modeled CO are only sensitive to the adjustment of CO emissions; improvements in surface NO 2 results are sensitive to VOC emission adjustments.These results are consistent with our previous conclusion that the VOC-sensitive environment dominates in eastern China.The ARIAs flights covered a large area (~10 4 km 2 ) in Hebei Province representing the regional nature of air pollution over the NCP.A case comparison of CMAQ_All run and Y12 measurements on 11 June 2016 (Figure 12) showed better results in both concentrations and vertical gradient of air pollutants (compared with Figure S5 in the Supplementary Materials), indicating the effectiveness of improving the emission inventories based on satellite observations.Table 3 summarizes the model performance of CMAQ as compared with aircraft measurements, and scatter plots for each CMAQ sensitivity experiment are shown in Figure S9 of the Supplementary Materials.The adjustments of the EDGAR emissions with satellite observations substantially improved simulations of ozone pollution, with the root mean square error (RSME) decreasing from 25.1 ppbv (the CMAQ_baseline run) to 21.2 ppbv (the CMAQ_All run), and the mean ratio of CMAQ simulations to aircraft observations increasing from 0.75 to 0.82.The model performance of CO has also been improved, with the RMSE decreasing from 247.0 ppbv to 203.6 ppbv, and the mean ratio increasing from 0.40 to 0.66.For nitrogen compounds including NO 2 , NO, and NO y , the adjustments of EDGAR emissions had small impacts on improving the CMAQ performance.The reason could be that the ozone photochemistry is mainly VOC-sensitive over the NCP, so the adjustments of NO x emissions have limited impacts close to sources.
CMAQ_baseline run) to 21.2 ppbv (the CMAQ_All run), and the mean ratio of CMAQ simulations to aircraft observations increasing from 0.75 to 0.82.The model performance of CO has also been improved, with the RMSE decreasing from 247.0 ppbv to 203.6 ppbv, and the mean ratio increasing from 0.40 to 0.66.For nitrogen compounds including NO2, NO, and NOy, the adjustments of EDGAR emissions had small impacts on improving the CMAQ performance.The reason could be that the ozone photochemistry is mainly VOCsensitive over the NCP, so the adjustments of NOx emissions have limited impacts close to sources.

Conclusions and Discussion
The ARIAs campaign conducted aircraft measurements over the NCP and observed high concentrations of air pollutants, including O 3 , CO, and nitrogen compounds.CMAQ simulations driven by the 2010 EDGAR emissions substantially underestimated levels of ozone and its precursors in the campaign region.Analysis of emission enhancements of CO and NO x with respect to concurrent CO 2 measurements suggested that the usage of the 2010 EDGAR emissions for the 2016 ARIAs campaign could introduce substantial uncertainty due to the recent changes in anthropogenic emissions in China.Comparison of atmospheric columns of NO 2 from CMAQ simulations and satellite observations showed that NO x emissions decreased in megacities such as Beijing and Shanghai but increased in rural areas from 2010 to 2016.Similar analysis of HCHO and CO showed that the EDGAR VOC and CO emissions could be also underestimated in the NCP.The HCHO/NO 2 column ratio from OMI observations indicated tropospheric ozone production was mainly in the VOC-sensitive regime in the NCP, which was also confirmed by CMAQ simulations.To improve the model performance, we adjusted the EDGAR emissions over East Asia based on satellite observations.Better performance in simulating ozone and its precursors was achieved, as compared with surface and airborne observations.The ARIAs aircraft measurements revealed substantial amounts of tropospheric ozone observed in the free troposphere (i.e., >80 ppbv with 2-3 km AGL in Figure 2), which agrees with the global studies [148][149][150][151][152]. The high ozone concentrations aloft can be transported for a long distance and mixed downwind through the PBL development to impact surface ozone concentrations, as we observed in the U.S. [97].Our modeling approach used the default built-in profiles as boundary conditions, so the transport from regions outside China was not taken into consideration, which might have overestimated the local photochemical production of ozone in the NCP.more rigorous modeling study using dynamic boundary conditions downscaled from a global chemical transport model and advanced numerical tools such as the source apportionment technique is needed for future studies.
In the NCP, stratospheric intrusions could also have an important effect on surface ozone, with a typical enhancement of 20 ppbv [153].Our aircraft measurements did not observe high ozone events with low CO concentrations (Figure 3), which are used as an indicator of stratospheric intrusions to the surface (SITS).However, due to limited flights within the 2-month campaign, these SITS could be missed by the research flights.The CMAQ model used in this study does not have the stratospheric chemistry component, so the potential impacts from the ozone-rich stratospheric airmass on the surface ozone were also not simulated.
Recent reports [154] suggest that emissions of nitrous acid (HONO) from fertilized soils can contribute to ozone formation.The version of CMAQ used in this study does not consider this source, so inclusion of this process could further increase ozone production rates and sensitivity to VOCs in rural areas.
Both satellite observations and CMAQ simulations indicate that the VOC-sensitive chemistry dominates the ozone photochemical production in eastern China, so the rate of local ozone production is mainly controlled by the VOC emissions.In the past few years, despite implementation of control measures mainly on SO 2 and NO x , ozone concentrations have increased in China.Our study indicated that high NO x concentrations were pervasive in the PBL over rural areas of the NCP, where anthropogenic VOCs were also abundant.Reducing NOx emissions is essential to control ozone on the regional scale, but our model simulations indicated that reducing VOC emissions can lower the rate of photochemical smog production.
Currently, studies and regulations on anthropogenic VOC emissions are lacking, so with the expectation of further decreasing NO x emissions in China, more severe ozone pollution could be anticipated.Even though hydrocarbon controls can have a beneficial impact on the local rate of ozone production in the VOC-sensitive regime, the ozone levels will not decrease until NO x emissions are substantially lower, i.e., regulations on VOCs are needed in addition to continued controls on NOx emissions in China.These results can also partially explain why ozone pollution intensified in the past few years while PM 2.5 pollution abated.New datasets such as the updated 'bottom-up' emissions inventory for East Asia and high-resolution satellite observations such as TROPOMI products are needed to improve the modeling of ozone pollution in China, which can provide scientific evidence for future national and international regulations on air quality.

Figure 1 .
Figure 1.CMAQ domains for the ARIAs campaign (the proximate aircraft campaign area is demonstrated by a red dot).The coarse domain (d01, 36 km resolution) covers East Asia, and the nested domain (d02, 12 km resolution) focuses on eastern China.

Figure 4 .
Figure 4. Comparison of mean MDA8 ozone concentrations in May and June 2016 over eastern China.(a) Contour plot; the background stands for mean MDA8 ozone concentrations from the CMAQ_baseline run; the dots represent observed values from the CMEE network.(b) Scatter plot; blue line is the linear regression fitting.

Figure 4 .Figure 5 .
Figure 4. Comparison of mean MDA8 ozone concentrations in May and June 2016 over eastern China.(a) Contour plot; the background stands for mean MDA8 ozone concentrations from the CMAQ_baseline run; the dots represent observed values from the CMEE network.(b) Scatter plot; blue line is the linear regression fitting.

Figure 6 .
Figure 6.Comparison of total VOC concentrations from WAS samples and CMAQ simulations.Values are in units of parts per billion Carbon (ppbC).Black line shows the 1:1 ratio; red line stands for the linear regression fitting line.M_Diff: mean difference; R: correlation; NMB: normalized mean bias; NME: normalized mean error; RMSE: root mean square error; M_Ratio: mean ratio.

Figure 6 .
Figure 6.Comparison of total VOC concentrations from WAS samples and CMAQ simulations.Values are in units of parts per billion Carbon (ppbC).Black line shows the 1:1 ratio; red line stands for the linear regression fitting line.M_Diff: mean difference; R: correlation; NMB: normalized mean bias; NME: normalized mean error; RMSE: root mean square error; M_Ratio: mean ratio.

Air 2024, 2 , 13 Figure 6 .
Figure 6.Comparison of total VOC concentrations from WAS samples and CMAQ simulations.Values are in units of parts per billion Carbon (ppbC).Black line shows the 1:1 ratio; red line stands for the linear regression fitting line.M_Diff: mean difference; R: correlation; NMB: normalized mean bias; NME: normalized mean error; RMSE: root mean square error; M_Ratio: mean ratio.

Figure 7 .
Figure 7.Comparison of emission enhancements (EEs) from the ARIAs campaign and emission factors (EFs) from the EDGAR emission inventory.(a) ∆CO/∆CO 2 , (b) ∆NO x /∆CO 2 , (c) ∆NO x /∆CO.Blue histogram shows the distribution of EEs observed by the Y12 aircraft; red line shows the ratio calculated through averaging EFs from four sectors of the EDGAR anthropogenic emissions inventory.

Figure 9 .
Figure 9.Comparison of OMI NO2 and CMAQ NO2 columns averaged in May and June 2016.(a) Scatter plot of NO2 columns over eastern China.(b) Scatter plot of NO2 columns over the campaign area (error bars were calculated as the standard deviation of daily OMI products and daily CMAQ simulations during the 2-month campaign).

Figure 9 .
Figure 9.Comparison of OMI NO 2 and CMAQ NO 2 columns averaged in May and June 2016.(a) Scatter plot of NO 2 columns over eastern China.(b) Scatter plot of NO 2 columns over the campaign area (error bars were calculated as the standard deviation of daily OMI products and daily CMAQ simulations during the 2-month campaign).

Figure 10 .
Figure 10.Column HCHO/NO2 ratios over East Asia in 2016.(a) Ratio derived from collocated OMI HCHO and NO2 observation.(b) Ratio calculated from CMAQ HCHO and NO2 simulations with OMI quality information and averaging kernel (AK).

Figure 10 .
Figure 10.Column HCHO/NO 2 ratios over East Asia in 2016.(a) Ratio derived from collocated OMI HCHO and NO 2 observation.(b) Ratio calculated from CMAQ HCHO and NO 2 simulations with OMI quality information and averaging kernel (AK).

Figure 11 .
Figure 11.Similar to Figure 4, but compares MDA8 ozone concentrations from CMEE observations and the CMAQ_all run.(a) Contour plot; the background stands for mean MDA8 ozone

Figure 11 .
Figure 11.Similar to Figure 4, but compares MDA8 ozone concentrations from CMEE observations and the CMAQ_all run.(a) Contour plot; the background stands for mean MDA8 ozone concentrations from the CMAQ_baseline run; the dots represent observed values from the CMEE network.(b) Scatter plot; blue line is the linear regression fitting.

:
ARIAs flights over the NCP.Eleven research flights were conducted from May to mid-June 2016; Figure S2: CMEE Air quality monitoring sites in 2016; Figure S3: A plume observed over Xingtai during the flight on 11 June 2016; Figure S4: Comparison of mean CO and NO 2 concentrations in May and June 2016 over eastern China; Figure S5: A case study comparing aircraft observations and the CMAQ_baseline run results on 11 June 2016; Figure S6: Ratios of column contents of the CMAQ_baseline simulations and satellite observations; Figure S7: Similar to Figure 8, but shows differences between satellite observations and CMAQ simulations in May and June 2016; Figure S8: Similar to Figure 11, but compares CMEE observations and other CMAQ runs; Figure S9: Similar to Figure 5, but shows results from other CMAQ sensitivity experiments in the campaign region; Table • N, 115.02 • E, 30 m ASL), and Quzhou (QZ, 36.76 • N, 114.96 • E, 40 m ASL).

Table 1 .
List of CMAQ simulations with adjusted emissions based on satellite observations.Anthropogenic NO x , VOC, and CO emissions were adjusted using OMI NO 2 , OMI HCHO, and MOPITT CO, observations.

Table 2 .
Statistics of CMAQ performance in six sensitivity experiments compared with CMEE observations in eastern China.Slope and Correlation R (Corr.R) are from the linear regression analysis; NMB = normalized mean bias; RMSE = root mean square error.

Table 3 .
Statistics of CMAQ performance in six sensitivity experiments compared with ARIAs aircraft measurements over the NCP.