NSF NCAR Atmospheric Chemistry Observations & Modeling

Abstract. The Copernicus Atmosphere Monitoring Service (CAMS) reanalysis is the latest global reanalysis dataset of atmospheric composition produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), consisting of three-dimensional time-consistent atmospheric composition fields, including aerosols and chemical species. The dataset currently covers the period 2003–2016 and will be extended in the future by adding 1 year each year. A reanalysis for greenhouse gases is being produced separately. The CAMS reanalysis builds on the experience gained during the production of the earlier Monitoring Atmospheric Composition and Climate (MACC) reanalysis and CAMS interim reanalysis. Satellite retrievals of total column CO; tropospheric column NO2; aerosol optical depth (AOD); and total column, partial column and profile ozone retrievals were assimilated for the CAMS reanalysis with ECMWF's Integrated Forecasting System. The new reanalysis has an increased horizontal resolution of about 80 km and provides more chemical species at a better temporal resolution (3-hourly analysis fields, 3-hourly forecast fields and hourly surface forecast fields) than the previously produced CAMS interim reanalysis. The CAMS reanalysis has smaller biases compared with most of the independent ozone, carbon monoxide, nitrogen dioxide and aerosol optical depth observations used for validation in this paper than the previous two reanalyses and is much improved and more consistent in time, especially compared to the MACC reanalysis. The CAMS reanalysis is a dataset that can be used to compute climatologies, study trends, evaluate models, benchmark other reanalyses or serve as boundary conditions for regional models for past periods.

Abstract. The evolution of tropospheric ozone from 1850 to 2100 has been studied using data from Phase 6 of the Coupled Model Intercomparison Project (CMIP6). We evaluate long-term changes using coupled atmosphere–ocean chemistry–climate models, focusing on the CMIP Historical and ScenarioMIP ssp370 experiments, for which detailed tropospheric-ozone diagnostics were archived. The model ensemble has been evaluated against a suite of surface, sonde and satellite observations of the past several decades and found to reproduce well the salient spatial, seasonal and decadal variability and trends. The multi-model mean tropospheric-ozone burden increases from 247 ± 36 Tg in 1850 to a mean value of 356 ± 31 Tg for the period 2005–2014, an increase of 44 %. Modelled present-day values agree well with previous determinations (ACCENT: 336 ± 27 Tg; Atmospheric Chemistry and Climate Model Intercomparison Project, ACCMIP: 337 ± 23 Tg; Tropospheric Ozone Assessment Report, TOAR: 340 ± 34 Tg). In the ssp370 experiments, the ozone burden increases to 416 ± 35 Tg by 2100. The ozone budget has been examined over the same period using lumped ozone production (PO3) and loss (LO3) diagnostics. Both ozone production and chemical loss terms increase steadily over the period 1850 to 2100, with net chemical production (PO3-LO3) reaching a maximum around the year 2000. The residual term, which contains contributions from stratosphere–troposphere transport reaches a minimum around the same time before recovering in the 21st century, while dry deposition increases steadily over the period 1850–2100. Differences between the model residual terms are explained in terms of variation in tropopause height and stratospheric ozone burden.

Abstract. This study presents airborne in situ and satellite remote sensing climatologies of cirrus clouds and humidity. The climatologies serve as a guide to the properties of cirrus clouds, with the new in situ database providing detailed insights into boreal midlatitudes and the tropics, while the satellite-borne data set offers a global overview. To this end, an extensive, quality-checked data archive, the Cirrus Guide II in situ database, is created from airborne in situ measurements during 150 flights in 24 campaigns. The archive contains meteorological parameters, ice water content (IWC), ice crystal number concentration (Nice), ice crystal mean mass radius (Rice), relative humidity with respect to ice (RHice), and water vapor mixing ratio (H2O) for each of the flights. Depending on the parameter, the database has been extended by about a factor of 5–10 compared to earlier studies. As one result of our investigation, we show that the medians of Nice, Rice, and RHice have distinct patterns in the IWC–T parameter space. Lookup tables of these variables as functions of IWC and T can be used to improve global model cirrus representation and remote sensing retrieval methods. Another outcome of our investigation is that across all latitudes, the thicker liquid-origin cirrus predominate at lower altitudes, while at higher altitudes the thinner in situ-origin cirrus prevail. Further, examination of the radiative characteristics of in situ-origin and liquid-origin cirrus shows that the in situ-origin cirrus only slightly warm the atmosphere, while liquid-origin cirrus have a strong cooling effect. An important step in completing the Cirrus Guide II is the provision of the global cirrus Nice climatology, derived by means of the retrieval algorithm DARDAR-Nice from 10 years of cirrus remote sensing observations from satellite. The in situ measurement database has been used to evaluate and improve the satellite observations. We found that the global median Nice from satellite observations is almost 2 times higher than the in situ median and increases slightly with decreasing temperature. Nice medians of the most frequently occurring cirrus sorted by geographical regions are highest in the tropics, followed by austral and boreal midlatitudes, Antarctica, and the Arctic. Since the satellite climatologies enclose the entire spatial and temporal Nice occurrence, we could deduce that half of the cirrus are located in the lowest, warmest (224–242 K) cirrus layer and contain a significant amount of liquid-origin cirrus. A specific highlight of the study is the in situ observations of cirrus and humidity in the Asian monsoon anticyclone and the comparison to the surrounding tropics. In the convectively very active Asian monsoon, peak values of Nice and IWC of 30 cm−3 and 1000 ppmv are detected around the cold point tropopause (CPT). Above the CPT, ice particles that are convectively injected can locally add a significant amount of water available for exchange with the stratosphere. We found IWCs of up to 8 ppmv in the Asian monsoon in comparison to only 2 ppmv in the surrounding tropics. Also, the highest RHice values (120 %–150 %) inside of clouds and in clear sky are observed around and above the CPT. We attribute this to the high H2O mixing ratios (typically 3–5 ppmv) observed in the Asian monsoon compared to 1.5 to 3 ppmv found in the tropics. Above the CPT, supersaturations of 10 %–20 % are observed in regions of weak convective activity and up to about 50 % in the Asian monsoon. This implies that the water available for transport into the stratosphere might be higher than the expected saturation value.

We use a large laboratory, modeling, and field dataset to investigate the isoprene + O3 reaction, with the goal of better understanding the fates of the C1 and C4 Criegee intermediates in the atmosphere. Although ozonolysis can produce several distinct Criegee intermediates, the C1 stabilized Criegee (CH2OO, 61 ± 9%) is the only one observed to react bimolecularly. We suggest that the C4 Criegees have a low stabilization fraction and propose pathways for their decomposition. Both prompt and non-prompt reactions are important in the production of OH (28% ± 5%) and formaldehyde (81% ± 16%). The yields of unimolecular products (OH, formaldehyde, methacrolein (42 ± 6%) and methyl vinyl ketone (18 ± 6%)) are fairly insensitive to water, i.e., changes in yields in response to water vapor (≤4% absolute) are within the error of the analysis. We propose a comprehensive reaction mechanism that can be incorporated into atmospheric models, which reproduces laboratory data over a wide range of relative humidities. The mechanism proposes that CH2OO + H2O (k(H2O)∼ 1 × 10(-15) cm(3) molec(-1) s(-1)) yields 73% hydroxymethyl hydroperoxide (HMHP), 6% formaldehyde + H2O2, and 21% formic acid + H2O; and CH2OO + (H2O)2 (k(H2O)2∼ 1 × 10(-12) cm(3) molec(-1) s(-1)) yields 40% HMHP, 6% formaldehyde + H2O2, and 54% formic acid + H2O. Competitive rate determinations (kSO2/k(H2O)n=1,2∼ 2.2 (±0.3) × 10(4)) and field observations suggest that water vapor is a sink for greater than 98% of CH2OO in a Southeastern US forest, even during pollution episodes ([SO2] ∼ 10 ppb). The importance of the CH2OO + (H2O)n reaction is demonstrated by high HMHP mixing ratios observed over the forest canopy. We find that CH2OO does not substantially affect the lifetime of SO2 or HCOOH in the Southeast US, e.g., CH2OO + SO2 reaction is a minor contribution (<6%) to sulfate formation. Extrapolating, these results imply that sulfate production by stabilized Criegees is likely unimportant in regions dominated by the reactivity of ozone with isoprene. In contrast, hydroperoxide, organic acid, and formaldehyde formation from isoprene ozonolysis in those areas may be significant.

Abstract. Atmospheric carbon monoxide (CO) concentrations have been decreasing since 2000, as observed by both satellite- and ground-based instruments, but global bottom-up emission inventories estimate increasing anthropogenic CO emissions concurrently. In this study, we use a multi-species atmospheric Bayesian inversion approach to attribute satellite-observed atmospheric CO variations to its sources and sinks in order to achieve a full closure of the global CO budget during 2000–2017. Our observation constraints include satellite retrievals of the total column mole fraction of CO, formaldehyde (HCHO), and methane (CH4) that are all major components of the atmospheric CO cycle. Three inversions (i.e., 2000–2017, 2005–2017, and 2010–2017) are performed to use the observation data to the maximum extent possible as they become available and assess the consistency of inversion results to the assimilation of more trace gas species. We identify a declining trend in the global CO budget since 2000 (three inversions are broadly consistent during overlapping periods), driven by reduced anthropogenic emissions in the US and Europe (both likely from the transport sector), and in China (likely from industry and residential sectors), as well as by reduced biomass burning emissions globally, especially in equatorial Africa (associated with reduced burned areas). We show that the trends and drivers of the inversion-based CO budget are not affected by the inter-annual variation assumed for prior CO fluxes. All three inversions contradict the global bottom-up inventories in the world's top two emitters: for the sign of anthropogenic emission trends in China (e.g., here -0.8±0.5 % yr−1 since 2000, while the prior gives 1.3±0.4 % yr−1) and for the rate of anthropogenic emission increase in South Asia (e.g., here 1.0±0.6 % yr−1 since 2000, smaller than 3.5±0.4 % yr−1 in the prior inventory). The posterior model CO concentrations and trends agree well with independent ground-based observations and correct the prior model bias. The comparison of the three inversions with different observation constraints further suggests that the most complete constrained inversion that assimilates CO, HCHO, and CH4 has a good representation of the global CO budget, and therefore matches best with independent observations, while the inversion only assimilating CO tends to underestimate both the decrease in anthropogenic CO emissions and the increase in the CO chemical production. The global CO budget data from all three inversions in this study can be accessed from https://doi.org/10.6084/m9.figshare.c.4454453.v1 (Zheng et al., 2019).

The concentration of nitrogen oxides (NOx) plays a central role in controlling air quality. On a global scale, the primary sink of NOx is oxidation to form HNO3. Gas-phase HNO3 photolyses slowly with a lifetime in the troposphere of 10 days or more. However, several recent studies examining HONO chemistry have proposed that particle-phase HNO3 undergoes photolysis 10–300 times more rapidly than gas-phase HNO3. We present here constraints on the rate of particle-phase HNO3 photolysis based on observations of NOx and HNO3 collected over the Yellow Sea during the KORUS-AQ study in summer 2016. The fastest proposed photolysis rates are inconsistent with the observed NOx to HNO3 ratios. Negligible to moderate enhancements of the HNO3 photolysis rate in particles, 1–30 times faster than in the gas phase, are most consistent with the observations. Small or moderate enhancement of particle-phase HNO3 photolysis would not significantly affect the HNO3 budget but could help explain observations of HONO and NOx in highly aged air.

Abstract Nocturnal dinitrogen pentoxide (N 2 O 5 ) heterogeneous chemistry impacts regional air quality and the distribution and lifetime of tropospheric oxidants. Formed from the oxidation of nitrogen oxides, N 2 O 5 is heterogeneously lost to aerosol with a highly variable reaction probability, γ (N 2 O 5 ), dependent on aerosol composition and ambient conditions. Reaction products include soluble nitrate (HNO 3 or NO 3 − ) and nitryl chloride (ClNO 2 ). We report the first‐ever derivations of γ (N 2 O 5 ) from ambient wintertime aircraft measurements in the critically important nocturnal residual boundary layer. Box modeling of the 2015 Wintertime INvestigation of Transport, Emissions, and Reactivity (WINTER) campaign over the eastern United States derived 2,876 individual γ (N 2 O 5 ) values with a median value of 0.0143 and range of 2 × 10 −5 to 0.1751. WINTER γ (N 2 O 5 ) values exhibited the strongest correlation with aerosol water content, but weak correlations with other variables, such as aerosol nitrate and organics, suggesting a complex, nonlinear dependence on multiple factors, or an additional dependence on a nonobserved factor. This factor may be related to aerosol phase, morphology (i.e., core shell), or mixing state, none of which are commonly measured during aircraft field studies. Despite general agreement with previous laboratory observations, comparison of WINTER data with 14 literature parameterizations (used to predict γ (N 2 O 5 ) in chemical transport models) confirms that none of the current methods reproduce the full range of γ (N 2 O 5 ) values. Nine reproduce the WINTER median within a factor of 2. Presented here is the first field‐based, empirical parameterization of γ (N 2 O 5 ), fit to WINTER data, based on the functional form of previous parameterizations.

Anthropogenic mercury (Hg(0)) emissions oxidize to gaseous Hg(II) compounds, before deposition to Earth surface ecosystems. Atmospheric reduction of Hg(II) competes with deposition, thereby modifying the magnitude and pattern of Hg deposition. Global Hg models have postulated that Hg(II) reduction in the atmosphere occurs through aqueous-phase photoreduction that may take place in clouds. Here we report that experimental rainfall Hg(II) photoreduction rates are much slower than modelled rates. We compute absorption cross sections of Hg(II) compounds and show that fast gas-phase Hg(II) photolysis can dominate atmospheric mercury reduction and lead to a substantial increase in the modelled, global atmospheric Hg lifetime by a factor two. Models with Hg(II) photolysis show enhanced Hg(0) deposition to land, which may prolong recovery of aquatic ecosystems long after Hg emissions are lowered, due to the longer residence time of Hg in soils compared with the ocean. Fast Hg(II) photolysis substantially changes atmospheric Hg dynamics and requires further assessment at regional and local scales.

Abstract. We have compared a suite of recent global CO2 atmospheric inversion results to independent airborne observations and to each other, to assess their dependence on differences in northern extratropical (NET) vertical transport and to identify some of the drivers of model spread. We evaluate posterior CO2 concentration profiles against observations from the High-Performance Instrumented Airborne Platform for Environmental Research (HIAPER) Pole-to-Pole Observations (HIPPO) aircraft campaigns over the mid-Pacific in 2009–2011. Although the models differ in inverse approaches, assimilated observations, prior fluxes, and transport models, their broad latitudinal separation of land fluxes has converged significantly since the Atmospheric Carbon Cycle Inversion Intercomparison (TransCom 3) and the REgional Carbon Cycle Assessment and Processes (RECCAP) projects, with model spread reduced by 80 % since TransCom 3 and 70 % since RECCAP. Most modeled CO2 fields agree reasonably well with the HIPPO observations, specifically for the annual mean vertical gradients in the Northern Hemisphere. Northern Hemisphere vertical mixing no longer appears to be a dominant driver of northern versus tropical (T) annual flux differences. Our newer suite of models still gives northern extratropical land uptake that is modest relative to previous estimates (Gurney et al., 2002; Peylin et al., 2013) and near-neutral tropical land uptake for 2009–2011. Given estimates of emissions from deforestation, this implies a continued uptake in intact tropical forests that is strong relative to historical estimates (Gurney et al., 2002; Peylin et al., 2013). The results from these models for other time periods (2004–2014, 2001–2004, 1992–1996) and re-evaluation of the TransCom 3 Level 2 and RECCAP results confirm that tropical land carbon fluxes including deforestation have been near neutral for several decades. However, models still have large disagreements on ocean–land partitioning. The fossil fuel (FF) and the atmospheric growth rate terms have been thought to be the best-known terms in the global carbon budget, but we show that they currently limit our ability to assess regional-scale terrestrial fluxes and ocean–land partitioning from the model ensemble.

Atmospheric iodine causes tropospheric ozone depletion and aerosol formation, both of which have significant climate impacts, and is an essential dietary element for humans. However, the evolution of atmospheric iodine levels at decadal and centennial scales is unknown. Here, we report iodine concentrations in the RECAP ice-core (coastal East Greenland) to investigate how atmospheric iodine levels in the North Atlantic have evolved over the past 260 years (1750-2011), this being the longest record of atmospheric iodine in the Northern Hemisphere. The levels of iodine tripled from 1950 to 2010. Our results suggest that this increase is driven by anthropogenic ozone pollution and enhanced sub-ice phytoplankton production associated with the recent thinning of Arctic sea ice. Increasing atmospheric iodine has accelerated ozone loss and has considerably enhanced iodine transport and deposition to the Northern Hemisphere continents. Future climate and anthropogenic forcing may continue to amplify oceanic iodine emissions with potentially significant health and environmental impacts at global scale.

Laboratory studies of atmospheric chemistry characterize the nature of atmospherically relevant processes down to the molecular level, providing fundamental information used to assess how human activities drive environmental phenomena such as climate change, urban air pollution, ecosystem health, indoor air quality, and stratospheric ozone depletion. Laboratory studies have a central role in addressing the incomplete fundamental knowledge of atmospheric chemistry. This article highlights the evolving science needs for this community and emphasizes how our knowledge is far from complete, hindering our ability to predict the future state of our atmosphere and to respond to emerging global environmental change issues. Laboratory studies provide rich opportunities to expand our understanding of the atmosphere via collaborative research with the modeling and field measurement communities, and with neighboring disciplines.

) in early stage (<3 h) wildfire plume evolution. These results highlight the role of HONO as a major component of reactive nitrogen emissions from wildfires and the main driver of initial photochemical oxidation.

, secondary organic aerosols SOA) in different parts of the world in response to modified emissions of primary pollutants during the COVID-19 pandemic. We quantify the respective effects of the reductions in NOx and in volatile organic carbon (VOC) emissions, which, in most cases, affect oxidants in opposite ways. Using model simulations, we show that the level of NOx has been reduced by typically 40% in China during February 2020 and by similar amounts in many areas of Europe and North America in mid-March to mid-April 2020, in good agreement with space and surface observations. We show that, relative to a situation in which the emission reductions are ignored and despite the calculated increase in hydroxyl and peroxy radicals, the ozone concentration increased only in a few NOx-saturated regions (northern China, northern Europe, and the US) during the winter months of the pandemic when the titration of this molecule by NOx was reduced. In other regions, where ozone is NOx-controlled, the concentration of ozone decreased. SOA concentrations decrease in response to the concurrent reduction in the NOx and VOC emissions. The model also shows that atmospheric meteorological anomalies produced substantial variations in the concentrations of chemical species during the pandemic. In Europe, for example, a large fraction of the ozone increase in February 2020 was associated with meteorological anomalies, while in the North China Plain, enhanced ozone concentrations resulted primarily from reduced emissions of primary pollutants.

Abstract We calculate ozone radiative forcing (RF) and stratospheric temperature adjustments for the period 1850–2014 using the newly available Coupled Model Intercomparison Project phase 6 (CMIP6) ozone data set. The CMIP6 total ozone RF (1850s to 2000s) is 0.28 ± 0.17 W m −2 (which is 80% higher than our CMIP5 estimation), and 0.30 ± 0.17 W m −2 out to the present day (2014). The total ozone RF grows rapidly until the 1970s, slows toward the 2000s, and shows a renewed growth thereafter. Since the 1990s the shortwave RF exceeds the longwave RF. Global stratospheric ozone RF is positive between 1930 and 1970 and then turns negative but remains positive in the Northern Hemisphere throughout. Derived stratospheric temperature changes show a localized cooling in the subtropical lower stratosphere due to tropospheric ozone increases and cooling in the upper stratosphere due to ozone depletion by more than 1 K already prior to the satellite era (1980) and by more than 2 K out to the present day (2014).

Abstract Tropospheric O 3 has been decreasing across much of the eastern U.S. but has remained steady or even increased in some western regions. Recent increases in VOC and NO x emissions associated with the production of oil and natural gas (O&amp;NG) may contribute to this trend in some areas. The Northern Front Range of Colorado has regularly exceeded O 3 air quality standards during summertime in recent years. This region has VOC emissions from a rapidly developing O&amp;NG basin and low concentrations of biogenic VOC in close proximity to urban‐Denver NO x emissions. Here VOC OH reactivity (OHR), O 3 production efficiency (OPE), and an observationally constrained box model are used to quantify the influence of O&amp;NG emissions on regional summertime O 3 production. Analyses are based on measurements acquired over two summers at a central location within the Northern Front Range that lies between major regional O&amp;NG and urban emission sectors. Observational analyses suggest that mixing obscures any OPE differences in air primarily influenced by O&amp;NG or urban emission sector. The box model confirms relatively modest OPE differences that are within the uncertainties of the field observations. Box model results also indicate that maximum O 3 at the measurement location is sensitive to changes in NO x mixing ratio but also responsive to O&amp;NG VOC reductions. Combined, these analyses show that O&amp;NG alkanes contribute over 80% to the observed carbon mixing ratio, roughly 50% to the regional VOC OHR, and approximately 20% to regional photochemical O 3 production.

Abstract. The chemistry–climate model ECHAM-HAMMOZ contains a detailed representation of tropospheric and stratospheric reactive chemistry and state-of-the-art parameterizations of aerosols using either a modal scheme (M7) or a bin scheme (SALSA). This article describes and evaluates the model version ECHAM6.3-HAM2.3-MOZ1.0 with a focus on the tropospheric gas-phase chemistry. A 10-year model simulation was performed to test the stability of the model and provide data for its evaluation. The comparison to observations concentrates on the year 2008 and includes total column observations of ozone and CO from IASI and OMI, Aura MLS observations of temperature, HNO3, ClO, and O3 for the evaluation of polar stratospheric processes, an ozonesonde climatology, surface ozone observations from the TOAR database, and surface CO data from the Global Atmosphere Watch network. Global budgets of ozone, OH, NOx, aerosols, clouds, and radiation are analyzed and compared to the literature. ECHAM-HAMMOZ performs well in many aspects. However, in the base simulation, lightning NOx emissions are very low, and the impact of the heterogeneous reaction of HNO3 on dust and sea salt aerosol is too strong. Sensitivity simulations with increased lightning NOx or modified heterogeneous chemistry deteriorate the comparison with observations and yield excessively large ozone budget terms and too much OH. We hypothesize that this is an impact of potential issues with tropical convection in the ECHAM model.

The 2015 and 2020 ozone holes set record sizes in October-December. We show that these years, as well as other recent large ozone holes, still adhere to a fundamental recovery metric: the later onset of early spring ozone depletion as chlorine and bromine diminishes. This behavior is also captured in the Whole Atmosphere Chemistry Climate Model. We quantify observed recovery trends of the onset of the ozone hole and in the size of the September ozone hole, with good model agreement. A substantial reduction in ozone hole depth during September over the past decade is also seen. Our results indicate that, due to dynamical phenomena, it is likely that large ozone holes will continue to occur intermittently in October-December, but ozone recovery will still be detectable through the later onset, smaller, and less deep September ozone holes: metrics that are governed more by chemical processes.

Understanding the efficiency and variability of photochemical ozone (O3) production from western wildfire plumes is important to accurately estimate their influence on North American air quality. A set of photochemical measurements were made from the NOAA Twin Otter research aircraft as a part of the Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) experiment. We use a zero-dimensional (0-D) box model to investigate the chemistry driving O3 production in modeled plumes. Modeled afternoon plumes reached a maximum O3 mixing ratio of 140 ± 50 ppbv (average ± standard deviation) within 20 ± 10 min of emission compared to 76 ± 12 ppbv in 60 ± 30 min in evening plumes. Afternoon and evening maximum O3 isopleths indicate that plumes were near their peak in NOx efficiency. A radical budget describes the NOx volatile - organic compound (VOC) sensitivities of these plumes. Afternoon plumes displayed a rapid transition from VOC-sensitive to NOx-sensitive chemistry, driven by HOx (=OH + HO2) production from photolysis of nitrous acid (HONO) (48 ± 20% of primary HOx) and formaldehyde (HCHO) (26 ± 9%) emitted directly from the fire. Evening plumes exhibit a slower transition from peak NOx efficiency to VOC-sensitive O3 production caused by a reduction in photolysis rates and fire emissions. HOx production in evening plumes is controlled by HONO photolysis (53 ± 7%), HCHO photolysis (18 ± 9%), and alkene ozonolysis (17 ± 9%).

Abstract We examine in detail a 1 year global reanalysis of carbon monoxide (CO) that is based on joint assimilation of conventional meteorological observations and Measurement of Pollution in The Troposphere (MOPITT) multispectral CO retrievals in the Community Earth System Model (CESM). Our focus is to assess the impact to the chemical system when CO distribution is constrained in a coupled full chemistry‐climate model like CESM. To do this, we first evaluate the joint reanalysis (MOPITT Reanalysis) against four sets of independent observations and compare its performance against a reanalysis with no MOPITT assimilation (Control Run). We then investigate the CO burden and chemical response with the aid of tagged sectoral CO tracers. We estimate the total tropospheric CO burden in 2002 (from ensemble mean and spread) to be 371 ± 12% Tg for MOPITT Reanalysis and 291 ± 9% Tg for Control Run. Our multispecies analysis of this difference suggests that (a) direct emissions of CO and hydrocarbons are too low in the inventory used in this study and (b) chemical oxidation, transport, and deposition processes are not accurately and consistently represented in the model. Increases in CO led to net reduction of OH and subsequent longer lifetime of CH 4 (Control Run: 8.7 years versus MOPITT Reanalysis: 9.3 years). Yet at the same time, this increase led to 5–10% enhancement of Northern Hemisphere O 3 and overall photochemical activity via HO x recycling. Such nonlinear effects further complicate the attribution to uncertainties in direct emissions alone. This has implications to chemistry‐climate modeling and inversion studies of longer‐lived species.

Abstract. Wildfire smoke is one of the most significant concerns of human and environmental health, associated with its substantial impacts on air quality, weather, and climate. However, biomass burning emissions and smoke remain among the largest sources of uncertainties in air quality forecasts. In this study, we evaluate the smoke emissions and plume forecasts from 12 state-of-the-art air quality forecasting systems during the Williams Flats fire in Washington State, US, August 2019, which was intensively observed during the Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) field campaign. Model forecasts with lead times within 1 d are intercompared under the same framework based on observations from multiple platforms to reveal their performance regarding fire emissions, aerosol optical depth (AOD), surface PM2.5, plume injection, and surface PM2.5 to AOD ratio. The comparison of smoke organic carbon (OC) emissions suggests a large range of daily totals among the models, with a factor of 20 to 50. Limited representations of the diurnal patterns and day-to-day variations of emissions highlight the need to incorporate new methodologies to predict the temporal evolution and reduce uncertainty of smoke emission estimates. The evaluation of smoke AOD (sAOD) forecasts suggests overall underpredictions in both the magnitude and smoke plume area for nearly all models, although the high-resolution models have a better representation of the fine-scale structures of smoke plumes. The models driven by fire radiative power (FRP)-based fire emissions or assimilating satellite AOD data generally outperform the others. Additionally, limitations of the persistence assumption used when predicting smoke emissions are revealed by substantial underpredictions of sAOD on 8 August 2019, mainly over the transported smoke plumes, owing to the underestimated emissions on 7 August. In contrast, the surface smoke PM2.5 (sPM2.5) forecasts show both positive and negative overall biases for these models, with most members presenting more considerable diurnal variations of sPM2.5. Overpredictions of sPM2.5 are found for the models driven by FRP-based emissions during nighttime, suggesting the necessity to improve vertical emission allocation within and above the planetary boundary layer (PBL). Smoke injection heights are further evaluated using the NASA Langley Research Center's Differential Absorption High Spectral Resolution Lidar (DIAL-HSRL) data collected during the flight observations. As the fire became stronger over 3–8 August, the plume height became deeper, with a day-to-day range of about 2–9 km a.g.l. However, narrower ranges are found for all models, with a tendency of overpredicting the plume heights for the shallower injection transects and underpredicting for the days showing deeper injections. The misrepresented plume injection heights lead to inaccurate vertical plume allocations along the transects corresponding to transported smoke that is 1 d old. Discrepancies in model performance for surface PM2.5 and AOD are further suggested by the evaluation of their ratio, which cannot be compensated for by solely adjusting the smoke emissions but are more attributable to model representations of plume injections, besides other possible factors including the evolution of PBL depths and aerosol optical property assumptions. By consolidating multiple forecast systems, these results provide strategic insight on pathways to improve smoke forecasts.