Search

To quantify changes in air pollution over Europe at the 2050 horizon, we designed a comprehensive modelling system that captures the external factors considered to be most relevant, and that relies on up-to-date and consistent sets of air pollution and climate policy scenarios. Global and regional climate as well as global chemistry simulations are based on the recent representative concentration pathways (RCP) produced for the Fifth Assessment Report (AR5) of the IPCC (Intergovernmental Panel on Climate Change) whereas regional air quality modelling is based on the updated emissions scenarios produced in the framework of the Global Energy Assessment. We explored two diverse scenarios: a reference scenario where climate policies are absent and a mitigation scenario which limits global temperature rise to within 2 °C by the end of this century. This first assessment of projected air quality and climate at the regional scale based on CMIP5 (5th Coupled Model Intercomparison Project) climate simulations is in line with the existing literature using CMIP3. The discrepancy between air quality simulations obtained with a climate model or with meteorological reanalyses is pointed out. Sensitivity simulations show that the main factor driving future air quality projections is air pollutant emissions, rather than climate change or intercontinental transport of pollution. Whereas the well documented "climate penalty" that weights upon ozone (increase of ozone pollution with global warming) over Europe is confirmed, other features appear less robust compared to the literature, such as the impact of climate on PM2.5. The quantitative disentangling of external factors shows that, while several published studies focused on the climate penalty bearing upon ozone, the contribution of the global ozone burden is somewhat overlooked in the literature.

We examine the response of Arctic gas and aerosol concentrations to perturbations in pollutant emissions from Europe, East and South Asia, and North America using results from a coordinated model intercomparison. These sensitivities to regional emissions (mixing ratio change per unit emission) vary widely across models and species. Intermodel differences are systematic, however, so that the relative importance of different regions is robust. North America contributes the most to Arctic ozone pollution. For aerosols and CO, European emissions dominate at the Arctic surface but East Asian emissions become progressively more important with altitude, and are dominant in the upper troposphere. Sensitivities show strong seasonality: surface sensitivities typically maximize during boreal winter for European and during spring for East Asian and North American emissions. Mid-tropospheric sensitivities, however, nearly always maximize during spring or summer for all regions. Deposition of black carbon (BC) onto Greenland is most sensitive to North American emissions. North America and Europe each contribute ~40% of total BC deposition to Greenland, with ~20% from East Asia. Elsewhere in the Arctic, both sensitivity and total BC deposition are dominated by European emissions. Model diversity for aerosols is especially large, resulting primarily from differences in aerosol physical and chemical processing (including removal). Comparison of modeled aerosol concentrations with observations indicates problems in the models, and perhaps, interpretation of the measurements. For gas phase pollutants such as CO and O3, which are relatively well-simulated, the processes contributing most to uncertainties depend on the source region and altitude examined. Uncertainties in the Arctic surface CO response to emissions perturbations are dominated by emissions for East Asian sources, while uncertainties in transport, emissions, and oxidation are comparable for European and North American sources. At higher levels, model-to-model variations in transport and oxidation are most important. Differences in photochemistry appear to play the largest role in the intermodel variations in Arctic ozone sensitivity, though transport also contributes substantially in the mid-troposphere.

Fifteen chemistry-transport models are used to quantify, for the first time, the export of oxidised nitrogen (NOy) to and from four regions (Europe, North America, South Asia, and East Asia), and to estimate the uncertainty in the results. Between 12 and 24% of the NOx emitted is exported from each region annually. The strongest impact of each source region on a foreign region is: Europe on East Asia, North America on Europe, South Asia on East Asia, and East Asia on North America. Europe exports the most NOy, and East Asia the least. East Asia receives the most NOy from the other regions. Between 8 and 15% of NOx emitted in each region is transported over distances larger than 1000 km, with 3–10% ultimately deposited over the foreign regions.

Global tropospheric ozone distributions, budgets, and radiative forcings from an ensemble of 26 state-of-the-art atmospheric chemistry models have been intercompared and synthesized as part of a wider study into both the air quality and climate roles of ozone. Results from three 2030 emissions scenarios, broadly representing “optimistic,” “likely,” and “pessimistic” options, are compared to a base year 2000 simulation. This base case realistically represents the current global distribution of tropospheric ozone. A further set of simulations considers the influence of climate change over the same time period by forcing the central emissions scenario with a surface warming of around 0.7K. The use of a large multimodel ensemble allows us to identify key areas of uncertainty and improves the robustness of the results. Ensemble mean changes in tropospheric ozone burden between 2000 and 2030 for the 3 scenarios range from a 5% decrease, through a 6% increase, to a 15% increase. The intermodel uncertainty (±1 standard deviation) associated with these values is about ±25%. Model outliers have no significant influence on the ensemble mean results. Combining ozone and methane changes, the three scenarios produce radiative forcings of −50, 180, and 300 mW m−2, compared to a CO2 forcing over the same time period of 800–1100 mW m−2. These values indicate the importance of air pollution emissions in short- to medium-term climate forcing and the potential for stringent/lax control measures to improve/worsen future climate forcing. The model sensitivity of ozone to imposed climate change varies between models but modulates zonal mean mixing ratios by ±5 ppbv via a variety of feedback mechanisms, in particular those involving water vapor and stratosphere-troposphere exchange. This level of climate change also reduces the methane lifetime by around 4%. The ensemble mean year 2000 tropospheric ozone budget indicates chemical production, chemical destruction, dry deposition and stratospheric input fluxes of 5100, 4650, 1000, and 550 Tg(O3) yr−1, respectively. These values are significantly different to the mean budget documented by the Intergovernmental Panel on Climate Change (IPCC) Third Assessment Report (TAR). The mean ozone burden (340 Tg(O3)) is 10% larger than the IPCC TAR estimate, while the mean ozone lifetime (22 days) is 10% shorter. Results from individual models show a correlation between ozone burden and lifetime, and each model's ozone burden and lifetime respond in similar ways across the emissions scenarios. The response to climate change is much less consistent. Models show more variability in the tropics compared to midlatitudes. Some of the most uncertain areas of the models include treatments of deep tropical convection, including lightning NOx production; isoprene emissions from vegetation and isoprene's degradation chemistry; stratosphere-troposphere exchange; biomass burning; and water vapor concentrations.

We analyze present-day and future carbon monoxide (CO) simulations in 26 state-of-the-art atmospheric chemistry models run to study future air quality and climate change. In comparison with near-global satellite observations from the MOPITT instrument and local surface measurements, the models show large underestimates of Northern Hemisphere (NH) extratropical CO, while typically performing reasonably well elsewhere. The results suggest that year-round emissions, probably from fossil fuel burning in east Asia and seasonal biomass burning emissions in south-central Africa, are greatly underestimated in current inventories such as IIASA and EDGAR3.2. Variability among models is large, likely resulting primarily from intermodel differences in representations and emissions of nonmethane volatile organic compounds (NMVOCs) and in hydrologic cycles, which affect OH and soluble hydrocarbon intermediates. Global mean projections of the 2030 CO response to emissions changes are quite robust. Global mean midtropospheric (500 hPa) CO increases by 12.6 ± 3.5 ppbv (16%) for the high-emissions (A2) scenario, by 1.7 ± 1.8 ppbv (2%) for the midrange (CLE) scenario, and decreases by 8.1 ± 2.3 ppbv (11%) for the low-emissions (MFR) scenario. Projected 2030 climate changes decrease global 500 hPa CO by 1.4 ± 1.4 ppbv. Local changes can be much larger. In response to climate change, substantial effects are seen in the tropics, but intermodel variability is quite large. The regional CO responses to emissions changes are robust across models, however. These range from decreases of 10–20 ppbv over much of the industrialized NH for the CLE scenario to CO increases worldwide and year-round under A2, with the largest changes over central Africa (20–30 ppbv), southern Brazil (20–35 ppbv) and south and east Asia (30–70 ppbv). The trajectory of future emissions thus has the potential to profoundly affect air quality over most of the world's populated areas.

Within ACCENT, a European Network of Excellence, eighteen atmospheric models from the U.S., Europe, and Japan calculated present (2000) and future (2030) concentrations of ozone at the Earth's surface with hourly temporal resolution. Comparison of model results with surface ozone measurements in 14 world regions indicates that levels and seasonality of surface ozone in North America and Europe are characterized well by global models, with annual average biases typically within 5–10 nmol/mol. However, comparison with rather sparse observations over some regions suggest that most models overestimate annual ozone by 15–20 nmol/mol in some locations. Two scenarios from the International Institute for Applied Systems Analysis (IIASA) and one from the Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios (IPCC SRES) have been implemented in the models. This study focuses on changes in near-surface ozone and their effects on human health and vegetation. Different indices and air quality standards are used to characterise air quality. We show that often the calculated changes in the different indices are closely inter-related. Indices using lower thresholds are more consistent between the models, and are recommended for global model analysis. Our analysis indicates that currently about two-thirds of the regions considered do not meet health air quality standards, whereas only 2–4 regions remain below the threshold. Calculated air quality exceedances show moderate deterioration by 2030 if current emissions legislation is followed and slight improvements if current emissions reduction technology is used optimally. For the "business as usual" scenario severe air quality problems are predicted. We show that model simulations of air quality indices are particularly sensitive to how well ozone is represented, and improved accuracy is needed for future projections. Additional measurements are needed to allow a more quantitative assessment of the risks to human health and vegetation from changing levels of surface ozone.

This paper investigates the acetone variability in the upper troposphere (UT) as sampled during the CARIBIC airborne experiment and simulated by the LMDz- INCA global chemistry climate model. The aim is to (1) describe spatial distribution and temporal variability of acetone; (2) propose benchmarks deduced from the observed data set; and (3) investigate the representativeness of the observational data set. According to the model results, South Asia (including part of the Indian Ocean, all of India, China, and the Indochinese peninsula) and Europe (including Mediterranean Sea) are net source regions of acetone, where nearly 25% of North Hemispheric (NH) primary emissions and nearly 40% of the NH chemical production of acetone take place. The impact of these net source regions on continental upper tropospheric acetone is studied by analysing CARIBIC observations of 2006 and 2007 when most flight routes stretched between Frankfurt (Germany) and Manila (Philippines), and by focussing over 3 sub-regions where acetone variability is strong: Europe-Mediterranean, Central South China and South China Sea. Important spatial variability was observed over different scales: (1) east-west positive gradient of annually averaged acetone vmr in UT over the Eurasian continent, namely a factor two increase from east to west; (2) ocean/continent contrast with 50% enhancement over the continents; (3) the acetone volume mixing ration (vmr) may vary in summer by more than 1000 pptv within only 5 latitude-longitude degrees; (4) the standard deviation for measurements acquired during a short flight sequence over a sub-region may reach 40 %. Temporal variability is also important: (1) the acetone volume mixing ratio (vmr) in the UT varies with the season, increasing from winter to summer by a factor 2 to 4; (2) a difference as large as 200 pptv may be observed between successive inbound and outbound flights over the same sub-region due to different flight specifications (trajectory in relation to the plume, time of day). A satisfactory agreement for the abundance of acetone is found between model results and observations, with e.g. only 30% overestimation of the annual average over Central-South China and the South China Sea (between 450 and 600 pptv), and an underestimation by less than 20% over Europe-Mediterranean (around 800 pptv). Consequently, annual budget terms could be computed with LMDz-INCA, yielding a global atmospheric burden of 7.2 Tg acetone, a 127 Tg yr−1 global source/sink strength, and a 21-day mean residence time. Moreover the study shows that LMDz-INCA can reproduce the impact of summer convection over China when boundary layer compounds are lifted to cruise altitude of 10– 11 km and higher. The consequent enhancement of acetone vmr during summer is reproduced by LMDz-INCA, to reach agreement on an observed maximum of 970±400 pptv (average during each flight sequence over the defined zone ± standard deviation). The summer enhancement of acetone is characterized by a high spatial and temporal heterogeneity, showing the necessity to increase the airborne measurement frequency over Central-South China and the South China Sea in August and September, when the annual maximum is expected (daily average model values reaching potentially 3000 pptv). In contrast, the annual cycle in the UT over Europe-Mediterranean is not reproduced by LMDz-INCA, in particular the observed summer enhancement of acetone to 1400±400 pptv after long-range transport of free tropospheric air masses over North Atlantic Ocean is not reproduced. In view of the agreement on the acetone annual cycle at surface level, this disagreement in UT over Europe indicates misrepresentation of simulated transport of primary acetone or biased spatial distribution of acetone chemical sinks and secondary sources. The sink and source budget in long-range transported free tropospheric air masses may be studied by analysing atmospheric chemical composition observed by CARIBIC in summer flights between North America and Europe.

As part of the Hemispheric Transport of Air Pollution (HTAP; http:// www.htap.org) project, we analyze results from 15 global and 1 hemispheric chemical transport models and compare these to Clean Air Status and Trends Network (CASTNet) observations in the United States (US) for 2001. Using the policy-relevant maximum daily 8-h average ozone (MDA8 O3) statistic, the multi-model ensemble represents the observations well (mean r2=0.57, ensemble bias = +4.1 ppbv for all US regions and all seasons) despite a wide range in the individual model results. Correlations are strongest in the northeastern US during spring and fall (r2=0.68); and weakest in the midwestern US in summer (r2=0.46). However, large positive mean biases exist during summer for all eastern US regions, ranging from 10–20 ppbv, and a smaller negative bias is present in the western US during spring (~3 ppbv). In nearly all other regions and seasons, the biases of the model ensemble simulations are ≤5 ppbv. Sensitivity simulations in which anthropogenic O3-precursor emissions (NOx + NMVOC + CO + aerosols) were decreased by 20% in four source regions: East Asia (EA), South Asia (SA), Europe (EU) and North America (NA) show that the greatest response of MDA8 O3 to the summed foreign emissions reductions occurs during spring in the West (0.9 ppbv reduction due to 20% emissions reductions from EA + SA + EU). East Asia is the largest contributor to MDA8 O3 at all ranges of the O3 distribution for most regions (typically ~0.45 ppbv) followed closely by Europe. The exception is in the northeastern US where emissions reductions in EU had a slightly greater influence than EA emissions, particularly in the middle of the MDA8 O3 distribution (response of ~0.35 ppbv between 35–55 ppbv). EA and EU influences are both far greater (about 4x) than that from SA in all regions and seasons. In all regions and seasons O3-precursor emissions reductions of 20% in the NA source region decrease MDA8 O3 the most – by a factor of 2 to nearly 10 relative to foreign emissions reductions. The O3 response to anthropogenic NA emissions is greatest in the eastern US during summer at the high end of the O3 distribution (5–6 ppbv for 20% reductions). While the impact of foreign emissions on surface O3 in the US is not negligible – and is of increasing concern given the recent growth in Asian emissions – domestic emissions reductions remain a far more effective means of decreasing MDA8 O3 values, particularly those above 75 ppb (the current US standard).

This study describes a simple parameterization to estimate regionally averaged changes in surface ozone due to past or future changes in anthropogenic precursor emissions based on results from 14 global chemistry transport models. The method successfully reproduces the results of full simulations with these models. For a given emission scenario it provides the ensemble mean surface ozone change, a regional source attribution for each change, and an estimate of the associated uncertainty as represented by the variation between models. Using the Representative Concentration Pathway (RCP) emission scenarios as an example, we show how regional surface ozone is likely to respond to emission changes by 2050 and how changes in precursor emissions and atmospheric methane contribute to this. Surface ozone changes are substantially smaller than expected with the SRES A1B, A2 and B2 scenarios, with annual global mean reductions of as much as 2 ppb by 2050 vs. increases of 4–6 ppb under SRES, and this reflects the assumptions of more stringent precursor emission controls under the RCP scenarios. We find an average difference of around 5 ppb between the outlying RCP 2.6 and RCP 8.5 scenarios, about 75% of which can be attributed to differences in methane abundance. The study reveals the increasing importance of limiting atmospheric methane growth as emissions of other precursors are controlled, but highlights differences in modelled ozone responses to methane changes of as much as a factor of two, indicating that this remains a major uncertainty in current models.

Understanding the surface O3 response over a “receptor” region to emission changes over a foreign “source” region is key to evaluating the potential gains from an international approach to abate ozone (O3) pollution. We apply an ensemble of 21 global and hemispheric chemical transport models to estimate the spatial average surface O3 response over east Asia (EA), Europe (EU), North America (NA), and south Asia (SA) to 20% decreases in anthropogenic emissions of the O3 precursors, NOx, NMVOC, and CO (individually and combined), from each of these regions. We find that the ensemble mean surface O3 concentrations in the base case (year 2001) simulation matches available observations throughout the year over EU but overestimates them by >10 ppb during summer and early fall over the eastern United States and Japan. The sum of the O3 responses to NOx, CO, and NMVOC decreases separately is approximately equal to that from a simultaneous reduction of all precursors. We define a continental-scale “import sensitivity” as the ratio of the O3 response to the 20% reductions in foreign versus “domestic” (i.e., over the source region itself) emissions. For example, the combined reduction of emissions from the three foreign regions produces an ensemble spatial mean decrease of 0.6 ppb over EU (0.4 ppb from NA), less than the 0.8 ppb from the reduction of EU emissions, leading to an import sensitivity ratio of 0.7. The ensemble mean surface O3 response to foreign emissions is largest in spring and late fall (0.7–0.9 ppb decrease in all regions from the combined precursor reductions in the three foreign regions), with import sensitivities ranging from 0.5 to 1.1 (responses to domestic emission reductions are 0.8–1.6 ppb). High O3 values are much more sensitive to domestic emissions than to foreign emissions, as indicated by lower import sensitivities of 0.2 to 0.3 during July in EA, EU, and NA when O3 levels are typically highest and by the weaker relative response of annual incidences of daily maximum 8-h average O3 above 60 ppb to emission reductions in a foreign region (

We have analysed time-slice simulations from 17 global models, participating in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), to explore changes in present-day (2000) hydroxyl radical (OH) concentration and methane (CH4) lifetime relative to preindustrial times (1850) and to 1980. A comparison of modeled and observation-derived methane and methyl chloroform lifetimes suggests that the present-day global multi-model mean OH concentration is overestimated by 5 to 10% but is within the range of uncertainties. The models consistently simulate higher OH concentrations in the Northern Hemisphere (NH) compared with the Southern Hemisphere (SH) for the present-day (2000; inter-hemispheric ratios of 1.13 to 1.42), in contrast to observation-based approaches which generally indicate higher OH in the SH although uncertainties are large. Evaluation of simulated carbon monoxide (CO) concentrations, the primary sink for OH, against ground-based and satellite observations suggests low biases in the NH that may contribute to the high north–south OH asymmetry in the models. The models vary widely in their regional distribution of present-day OH concentrations (up to 34%). Despite large regional changes, the multi-model global mean (mass-weighted) OH concentration changes little over the past 150 yr, due to concurrent increases in factors that enhance OH (humidity, tropospheric ozone, nitrogen oxide (NOx) emissions, and UV radiation due to decreases in stratospheric ozone), compensated by increases in OH sinks (methane abundance, carbon monoxide and non-methane volatile organic carbon (NMVOC) emissions). The large inter-model diversity in the sign and magnitude of preindustrial to present-day OH changes (ranging from a decrease of 12.7% to an increase of 14.6%) indicate that uncertainty remains in our understanding of the long-term trends in OH and methane lifetime. We show that this diversity is largely explained by the different ratio of the change in global mean tropospheric CO and NOx burdens (ΔCO/ΔNOx, approximately represents changes in OH sinks versus changes in OH sources) in the models, pointing to a need for better constraints on natural precursor emissions and on the chemical mechanisms in the current generation of chemistry-climate models. For the 1980 to 2000 period, we find that climate warming and a slight increase in mean OH (3.5 ± 2.2%) leads to a 4.3 ± 1.9% decrease in the methane lifetime. Analysing sensitivity simulations performed by 10 models, we find that preindustrial to present-day climate change decreased the methane lifetime by about four months, representing a negative feedback on the climate system. Further, we analysed attribution experiments performed by a subset of models relative to 2000 conditions with only one precursor at a time set to 1860 levels. We find that global mean OH increased by 46.4 ± 12.2% in response to preindustrial to present-day anthropogenic NOx emission increases, and decreased by 17.3 ± 2.3%, 7.6 ± 1.5%, and 3.1 ± 3.0% due to methane burden, and anthropogenic CO, and NMVOC emissions increases, respectively.

We present a systematic comparison of tropospheric NO2 from 17 global atmospheric chemistry models with three state-of-the-art retrievals from the Global Ozone Monitoring Experiment (GOME) for the year 2000. The models used constant anthropogenic emissions from IIASA/EDGAR3.2 and monthly emissions from biomass burning based on the 1997–2002 average carbon emissions from the Global Fire Emissions Database (GFED). Model output is analyzed at 10:30 local time, close to the overpass time of the ERS-2 satellite, and collocated with the measurements to account for sampling biases due to incomplete spatiotemporal coverage of the instrument. We assessed the importance of different contributions to the sampling bias: correlations on seasonal time scale give rise to a positive bias of 30–50% in the retrieved annual means over regions dominated by emissions from biomass burning. Over the industrial regions of the eastern United States, Europe and eastern China the retrieved annual means have a negative bias with significant contributions (between –25% and +10% of the NO2 column) resulting from correlations on time scales from a day to a month. We present global maps of modeled and retrieved annual mean NO2 column densities, together with the corresponding ensemble means and standard deviations for models and retrievals. The spatial correlation between the individual models and retrievals are high, typically in the range 0.81–0.93 after smoothing the data to a common resolution. On average the models underestimate the retrievals in industrial regions, especially over eastern China and over the Highveld region of South Africa, and overestimate the retrievals in regions dominated by biomass burning during the dry season. The discrepancy over South America south of the Amazon disappears when we use the GFED emissions specific to the year 2000. The seasonal cycle is analyzed in detail for eight different continental regions. Over regions dominated by biomass burning, the timing of the seasonal cycle is generally well reproduced by the models. However, over Central Africa south of the Equator the models peak one to two months earlier than the retrievals. We further evaluate a recent proposal to reduce the NOx emission factors for savanna fires by 40% and find that this leads to an improvement of the amplitude of the seasonal cycle over the biomass burning regions of Northern and Central Africa. In these regions the models tend to underestimate the retrievals during the wet season, suggesting that the soil emissions are higher than assumed in the models. In general, the discrepancies between models and retrievals cannot be explained by a priori profile assumptions made in the retrievals, neither by diurnal variations in anthropogenic emissions, which lead to a marginal reduction of the NO2 abundance at 10:30 local time (by 2.5–4.1% over Europe). Overall, there are significant differences among the various models and, in particular, among the three retrievals. The discrepancies among the retrievals (10–50% in the annual mean over polluted regions) indicate that the previously estimated retrieval uncertainties have a large systematic component. Our findings imply that top-down estimations of NOx emissions from satellite retrievals of tropospheric NO2 are strongly dependent on the choice of model and retrieval.

Acid-driven multiphase chemistry of isoprene epoxydiols (IEPOX), key isoprene oxidation products, with inorganic sulfate aerosol yields substantial amounts of secondary organic aerosol (SOA) through the formation of organosulfur compounds. The extent and implications of inorganic-to-organic sulfate conversion, however, are unknown. In this article, we demonstrate that extensive consumption of inorganic sulfate occurs, which increases with the IEPOX-to-inorganic sulfate concentration ratio (IEPOX/Sulfinorg), as determined by laboratory measurements. Characterization of the total sulfur aerosol observed at Look Rock, Tennessee, from 2007 to 2016 shows that organosulfur mass fractions will likely continue to increase with ongoing declines in anthropogenic Sulfinorg, consistent with our laboratory findings. We further demonstrate that organosulfur compounds greatly modify critical aerosol properties, such as acidity, morphology, viscosity, and phase state. These new mechanistic insights demonstrate that changes in SO2 emissions, especially in isoprene-dominated environments, will significantly alter biogenic SOA physicochemical properties. Consequently, IEPOX/Sulfinorg will play an important role in understanding the historical climate and determining future impacts of biogenic SOA on the global climate and air quality.

Methane is an important greenhouse gas, responsible for about 20% of the warming induced by long-lived greenhouse gases since pre-industrial times. By reacting with hydroxyl radicals, methane reduces the oxidizing capacity of the atmosphere and generates ozone in the troposphere. Although most sources and sinks of methane have been identified, their relative contributions to atmospheric methane levels are highly uncertain. As such, the factors responsible for the observed stabilization of atmospheric methane levels in the early 2000s, and the renewed rise after 2006, remain unclear. Here, we construct decadal budgets for methane sources and sinks between 1980 and 2010, using a combination of atmospheric measurements and results from chemical transport models, ecosystem models, climate chemistry models and inventories of anthropogenic emissions. The resultant budgets suggest that data-driven approaches and ecosystem models overestimate total natural emissions. We build three contrasting emission scenarios-which differ in fossil fuel and microbial emissions-to explain the decadal variability in atmospheric methane levels detected, here and in previous studies, since 1985. Although uncertainties in emission trends do not allow definitive conclusions to be drawn, we show that the observed stabilization of methane levels between 1999 and 2006 can potentially be explained by decreasing-to-stable fossil fuel emissions, combined with stable-to-increasing microbial emissions. We show that a rise in natural wetland emissions and fossil fuel emissions probably accounts for the renewed increase in global methane levels after 2006, although the relative contribution of these two sources remains uncertain. © 2013 Macmillan Publishers Limited.