We use a series of chemical transport model and chemistry climate model simulations to investigate the observed negative trends in MOPITT CO over several regions of the world, and to examine the consistency of time-dependent emission inventories with observations. We find that simulations driven by the MACCity inventory, used for the Chemistry Climate Modeling Initiative (CCMI), reproduce the negative trends in the CO column observed by MOPITT for 2000–2010 over the eastern United States and Europe. However, the simulations have positive trends over eastern China, in contrast to the negative trends observed by MOPITT. The model bias in CO, after applying MOPITT averaging kernels, contributes to the model–observation discrepancy in the trend over eastern China. This demonstrates that biases in a model's average concentrations can influence the interpretation of the temporal trend compared to satellite observations. The total ozone column plays a role in determining the simulated tropospheric CO trends. A large positive anomaly in the simulated total ozone column in 2010 leads to a negative anomaly in OH and hence a positive anomaly in CO, contributing to the positive trend in simulated CO. These results demonstrate that accurately simulating variability in the ozone column is important for simulating and interpreting trends in CO.
Carbon monoxide (CO) is an air pollutant that contributes to ozone formation and affects the oxidizing capacity of the troposphere (Thompson, 1992; Crutzen, 1973). Its primary loss is through reaction with OH, which leads to a lifetime of 1–2 months (Bey et al., 2001) and makes CO an excellent tracer of long-range transport. Both fossil fuel combustion and biomass burning are major sources of CO. The biomass burning source shows large interannual variability (van der Werf et al., 2010), while fossil fuel emissions typically change more gradually. The time-dependent MACCity inventory (Granier et al., 2011) shows decreases in CO emissions from the United States and Europe from 2000 to 2010 due to increasing pollution controls but increases in emissions from China. MACCity emissions for years after 2000 are based on the Representative Concentration Pathway (RCP) 8.5 (Riahi et al., 2007). The REAS (Kurokawa et al., 2013) and EDGAR4.2 (EC-JRC/PBL, 2011) inventories also show increasing CO emissions from China. The bottom-up inventory of Zhang et al. (2009) shows an 18 % increase in CO emissions from China from 2001 to 2006, and Zhao et al. (2012) estimate a 6 % increase between 2005 and 2009. However, there is considerable uncertainty in bottom-up inventories, and comparison of model hindcast simulations driven by bottom-up inventories with observations provides an important test of the time-dependent emission estimates. Space-based observations of CO are now available for over a decade and show trends at both hemispheric and regional scales. Warner et al. (2013) found significant negative trends in both background CO and recently emitted CO at 500 hPa over southern hemispheric oceans and northern hemispheric land and ocean in Atmospheric Infrared Sounder (AIRS) data. Worden et al. (2013) calculated trends in the CO column from several thermal infrared (TIR) instruments including MOPITT and AIRS. They found statistically significant negative trends over Europe, the eastern United States, and China for 2002–2012. He et al. (2013) also report a negative trend in MOPITT near-surface CO over western Maryland.
Surface concentrations of CO show downward trends over the United States driven by emission reductions (EPA, 2011), consistent with the space-based trends. Decreases in the partial column of CO from FTIR stations in Europe also show decreases from 1996 to 2006, consistent with emissions decreases (Angelbratt et al., 2011). Yoon and Pozzer (2014) found that a model simulation of 2001 to 2010 reproduced negative trends in surface CO over the eastern United States and western Europe, but showed a positive trend in surface CO over southern Asia.
The cause of the negative trend over China seen in MOPITT and AIRS data is
uncertain. The trend is consistent with the results of Li and Liu (2011),
who found decreases in surface CO measurements in Beijing, and with
decreases in CO emissions in 2008 inferred from the correlation of CO with
CO
This study examines whether global hindcast simulations can reproduce the trends and variability in carbon monoxide seen in the MOPITT record. We examine the role of averaging kernels and the contribution of trends at different altitudes to the trends observed by MOPITT. We then examine the impact of OH variability on the simulated trends in CO.
The MOPITT instrument onboard the Terra Satellite provides the longest satellite-based record of atmospheric CO, with observations available from March 2000 to present. It provides nearly global coverage every 3 days (Edwards et al., 2004). We use the monthly Level 3 daytime column data from the Version 5 TIR product, which has negligible drift in the bias over time (Deeter et al., 2013). The Level 3 data are a gridded product and include the a priori and averaging kernel for each grid box. Supplemental Fig. S1 shows the MOPITT column averaging kernels averaged over four regions. The column averaging kernels depend on the observed scene, and vary year to year as well as seasonally. The dependence of the column averaging kernels on the CO mixing ratio profile (Deeter, 2009) explains the high values in the lower troposphere over eastern China in winter.
We calculate trends and deseasonalized anomalies for the eastern United States,
Europe, and eastern China regions described by Worden et al. (2013). Trends
that differ from zero by more than the 2
We use a suite of chemistry climate model (CCM) and chemical transport model (CTM) simulations to interpret the observed trends. The Global Modeling Initiative (GMI) CTM includes both tropospheric (Duncan et al., 2007) and stratospheric (Strahan et al., 2007) chemistry, including over 400 reactions and 124 chemical species. Meteorology for the GMI simulations comes from the Modern-Era Retrospective Analysis for Research and Applications (MERRA) (Rienecker et al., 2011). The GEOS-5 Chemistry Climate Model (GEOSCCM) (Oman et al., 2011) incorporates the GMI chemical mechanism into the GEOS-5 atmospheric general circulation model (AGCM). The GEOSCCM simulations are forced by observed sea surface temperatures (SSTs) from Reynolds et al. (2002).
The Community Earth System Model, CESM1 CAM4-chem, includes 191 chemical tracers and over 400 reactions for both troposphere and stratosphere (Tilmes et al., 2016). The model can be run fully coupled to a free-running ocean, with prescribed SSTs, or with nudged meteorology from GEOS-5 or MERRA analysis. CESM1 CAM4-chem is further coupled to the land model, providing biogenic emissions from the Model of Emissions and Aerosols from Nature (MEGAN), version 2.1 (Guenther et al., 2012).
Several simulations were conducted as part of the Chemistry-Climate Model
Initiative (CCMI) project (Eyring et al., 2013). These include the Ref-C1
simulation of the GEOSCCM and a Ref-C1 CESM1 CAM4-Chem simulation, hereafter
called G-Ref-C1 and C-Ref-C1, respectively, and the Ref-C1-SD simulation of
the GMI CTM. Both the Ref-C1 and the Ref-C1-SD simulations use
time-dependent anthropogenic and biomass burning emissions from the MACCity
inventory (Granier et al., 2011), but the Ref-C1-SD simulations use
specified meteorology while the Ref-C1 simulations run with prescribed SSTs.
The MACCity inventory linearly interpolates the decadal anthropogenic
emissions from the ACCMIP inventory (Lamarque et al., 2010) for 2000, and
the RCP8.5 emissions for 2005 and 2010, to each year in between. The MACCity
biomass burning emissions have year-to-year variability based on the GFED-v2
(van der Werf et al., 2006) inventory. From 2000 to 2010, CO emissions in
the MACCity inventory decreased from 31 to 11 Tg yr
Given the uncertainty in CO emissions, we conduct a GMI CTM simulation using
an alternative time-dependent emissions scenario, called AltEmis. This
simulation is described in detail in Strode et al. (2015b). Briefly,
anthropogenic emissions include time dependence based on EPA (
Trends in the CO emissions used in the Ref-C1 and Ref-C1-SD simulations (blue bars) and AltEmis simulation (purple bars) over 2000–2010 for the United States, Europe, China, and the world.
Description of simulations.
We regrid the model output to the MOPITT grid and convolve the simulated CO
with the MOPITT averaging kernels and a priori in order to compare the
simulated and observed CO columns. The averaging kernels are space- and time-dependent. We use the following equation from Deeter et al. (2013):
We deseasonalize the simulated CO columns and calculate their linear trend following the same procedure that we applied to the MOPITT CO. Months that do not have MOPITT data (June–July 2001 and August–September 2009) are excluded from the analysis of the model trends as well.
The Ref-C1 and Ref-C1-SD simulations requested by CCMI extend until 2010. However, the MACCity biomass burning emissions extend only until 2008. CAM4-Chem therefore repeated the biomass burning emissions for 2008 for years 2009–2010. In contrast, the GEOSCCM Ref-C1 and GMI Ref-C1-SD simulations used emissions from GFED3 (van der Werf et al., 2010) for years after 2008. Some simulations were available through 2011, while others ended in 2010. We therefore report results for 2000–2010, but note that extending the analysis through 2011 does not alter the conclusions.
Regional trends and correlations:
The hindcast simulations driven by MACCity emissions (G-Ref-C1, Ref-C1-SD, and C-Ref-C1) show negative trends in CO over the US and Europe that agree with the observed slope from MOPITT within the uncertainty (Fig. 2, Table 2). The MOPITT trends for both regions are statistically significant for both regions, as shown by Worden et al. (2013). These results are consistent with the findings of Yin et al. (2015), whose inversion of MOPITT data showed a posteriori trends in CO emissions over the US and western Europe that were consistent with but slightly larger than the a priori trends. The EmFix hindcast shows a positive, though non-significant, trend for both regions, indicating that the decrease in CO emissions is necessary for reproducing the downward trend in the CO column. The AltEmis simulation fails to produce the negative trends, despite including negative trends in regional emissions for both the US and Europe. The impact of these negative regional trends is insufficient to overcome the positive global emission trend in the AltEmis scenario (Fig. 1), leading to positive trends in CO.
The time series and trends (left column) and deseasonalized
monthly anomalies (right column) of the CO column from MOPITT (black), the
MOPITT a priori (gray), and simulated by G-Ref-C1 (red), Ref-C1-SD (blue),
EmFix (green), C-Ref-C1 (orange), and AltEmis (purple) for 2000–2010. The
regions shown are
Figure 2 also reveals a negative bias in the simulated CO column between the models and MOPITT. A low bias in simulated CO at northern latitudes is often present in global models (Naik et al., 2013) and may indicate a high bias in northern hemispheric OH (Strode et al., 2015a) or CO dry deposition (Stein et al., 2014), as well as an underestimate of CO emissions.
The deseasonalized anomalies in the MOPITT and simulated CO columns are shown in Fig. 2b and d; the correlation coefficients between the observed and simulated monthly anomalies are presented in Table 2b. The highest correlations are for the AltEmis and Ref-C1-SD simulations of the GMI CTM. This result is consistent with the use of year-specific meteorology, which we expect to better match the transport of particular years. The lowest correlations are for the EmFix simulation. This is expected since the EmFix simulation does not include inter-annual variability (IAV) in biomass burning. The IAV in biomass burning makes a large contribution to the IAV of CO (Voulgarakis et al., 2015).
The role of biomass burning in driving the CO variability is even more
evident at the hemispheric scale. Figure 2g and h show the anomalies in MOPITT
and the simulations for the Northern Hemisphere (0–60
Observations from MOPITT show a negative trend in the CO column over eastern China for 2002–2012 (Worden et al., 2013). The negative trend for the years 2000–2014 exceeds that for 2000–2010 (Table 2), showing that it is not driven solely by temporary emission reductions in 2008. Our simulations do not reproduce this trend, and instead show increases in the CO column (Fig. 2e), which is expected given that CO emissions from China increase in four of the five simulations. The anomalies (Fig. 2f) show that the discrepancy in the simulated versus observed trends is driven largely by the failure of the simulations to capture the 2008 dip in the CO column, leading to an overestimate that continues through 2010. This suggests emission reductions in China during this time period are not adequately captured by the emission inventories. However, the good agreement between the observed and simulated decreases in CO for the Northern Hemisphere as a whole (Fig. 2g, h) suggest that on a global scale, the emission time series is reasonable. Consequently, we examine several other factors that may contribute to the difference in sign between the MOPITT and simulated CO trends.
Regional trends in CO are expected to vary with altitude, with surface concentrations most heavily influenced by local emissions. MOPITT TIR retrievals have higher sensitivity to CO in the mid-troposphere than at the surface (Deeter et al., 2004), so the trend in the MOPITT CO column will be weighted towards the trends in free tropospheric CO rather than near-surface CO. We quantify this impact on our Ref-C1-SD CO column trends by comparing the trend in the pure-model CO column with that of the simulated column convolved with the MOPITT averaging kernels.
The simulated CO trend over eastern China for 2000–2010 is positive (but not
significant) both with and without the averaging kernels, but application of
the MOPITT kernels increases the positive trend from 1.
Yoon et al. (2013) show that since the averaging kernels vary over time, a
bias between the true atmosphere and the a priori assumed by MOPITT can lead
to an artificial trend in the retrieved CO. Similarly, the bias between the
average simulated CO concentrations and the MOPITT a priori, evident in
Fig. 2, can lead to an artifact in the simulated CO trend when the
simulation is convolved with the MOPITT averaging kernels. This is due to
the changing contribution of the a priori when the vertical sensitivity
(averaging kernel) is varying in time. MOPITT vertical sensitivity varies
with time due to instrument degradation as well as the change in CO
abundance. The bias in CO varies with altitude, so if the vertical
sensitivity described by the averaging kernel changes, this will change the
value of the convolved CO column even if there were no changes in the CO
profile. Furthermore, changes in the averaging kernel result in more or less
weight placed on the a priori versus the CO simulated by the model. Thus, a
difference between the a priori and the model means that placing more (or
less) weight on the a priori will change the resulting value of
Deseasonalized monthly anomalies in the total ozone column (left), mean tropospheric OH (center), and CO column (right) from the EmFix simulation as a function of latitude and month.
The second term on the right-hand side shows that the larger the bias between the modeled CO and the a priori, the larger the impact of the changing averaging kernel.
We quantify this effect by convolving the simulated CO for each year with
the MOPITT averaging kernels for the year 2008, thus removing the effect of
the time dependence of the averaging kernels. The resulting trend,
0.
Figure 2 and Table 2 also show a positive trend in the GMI EmFix simulation
for eastern China. This larger trend in the EmFix simulation than the
Ref-C1-SD simulation indicates that the net decrease in emissions
contributes to decreasing CO over eastern China, consistent with the
observed negative trend, but other factors in the model cause an increase in
CO over eastern China even when all emissions are constant. Subtracting the
EmFix trend from the Ref-C1-SD trend shows that the changing emissions
contribute a CO trend of
Since the EmFix simulation shows a positive trend in the Northern
Hemisphere, we next examine the variability in the CO sink, OH. We also
examine variability in the total ozone column, since overhead ozone is a
major driver of OH variability (Duncan and Logan, 2008). Figure 3 shows the
variability in CO and OH in the EmFix simulation. The positive and negative
anomalies in CO correspond with the negative and positive anomalies,
respectively, in OH. The anomalies in OH are in turn inversely related to
anomalies in the total ozone column. The correlation coefficient between OH
and column ozone is
Monthly ozone column
The large anomaly in the simulated total ozone column in 2010 is
overestimated compared to observations. Figure 4 shows the time dependence
of the total ozone column from 30 to 60
While the impact of OH interannual variability on the apparent trend in CO
is clear in the EmFix simulation, this source of variability is partially
masked by large interannual variability in CO emissions in the other
simulations. We examine the correlation between the detrended and
deseasonalized CO anomalies from 10
We conducted a series of multi-year simulations to analyze the causes of the negative trends in MOPITT CO reported by Worden et al. (2013). Both CTM and CCM simulations driven by the MACCity emissions reproduce the observed trends over the eastern United States and Europe, providing confidence in the regional emission trends.
Monthly simulated CO column anomalies from the Ref-C1-SD
simulation as a function of CO emissions for 10
None of the simulations reproduce the observed negative trend over eastern
China. This negative trend persists even with the MOPITT data extended out
to 2014. The MOPITT averaging kernels are weighted towards the free
troposphere, where the relative importance of hemispheric versus local
trends is greater. However, our simulations indicate that this effect is
insufficient to explain the negative trends over China. Indeed, the negative
trend in MOPITT CO over eastern China (
Variability in emissions is the primary driver of year-to-year variability in simulated CO, but OH variability also plays a role. The simulated OH is anti-correlated with both CO and the total ozone column, highlighting the importance of realistic overhead ozone columns for accurately simulating CO variability and trends. In addition, further work is needed to understand recent changes in CO emissions from China.
This work was supported by NASA's Modeling, Analysis, and Prediction Program and computing resources from the NASA High-End Computing Program. We thank Bruce Van Aartsen for contributing to the GMI simulations. The CESM project is supported by the National Science Foundation and the Office of Science (BER) of the US Department of Energy. The MOPITT project is supported by the NASA Earth Observing System (EOS) Program. The National Center for Atmospheric Research (NCAR) is sponsored by the National Science Foundation.Edited by: A. Gettelman