ALBERT

Description: We describe the setup and first results of an inverse modelling system for atmospheric N2O, based on a four-dimensional variational (4DVAR) technique and the atmospheric transport zoom model TM5. We focus in this study on the European domain, utilizing a comprehensive set of quasi-continuous measurements over Europe, complemented by N2O measurements from the Earth System Research Laboratory of the National Oceanic and Atmospheric Administration (NOAA/ESRL) cooperative global air sampling network. Despite ongoing measurement comparisons among networks parallel measurements at a limited number of stations show that significant offsets exist among the different laboratories. Since the spatial gradients of N2O mixing ratios are of the same order of magnitude as these biases, the direct use of these biased datasets would lead to significant errors in the derived emissions. Therefore, in order to also use measurements with unknown offsets, a new bias correction scheme has been implemented within the TM5-4DVAR inverse modelling system, thus allowing the simultaneous assimilation of observations from different networks. The N2O bias corrections determined in the TM5-4DVAR system agree within ~0.1 ppb (dry-air mole fraction) with the bias derived from the measurements at monitoring stations where parallel NOAA discrete air samples are available. The N2O emissions derived for the northwest European and east European countries for 2006 show good agreement with the bottom-up emission inventories reported to the United Nations Framework Convention on Climate Change (UNFCCC). Moreover, the inverse model can significantly narrow the uncertainty range reported in N2O emission inventories for these countries, while the lack of measurements does not allow to reduce the uncertainties of emission estimates in southern Europe. Several sensitivity experiments were performed to test the robustness of the results. It is shown that also inversions without detailed a priori spatio-temporal emission distributions are capable to reproduce major regional emission patterns within the footprint of the existing atmospheric network, demonstrating the strong constraints of the atmospheric observations on the derived emissions.

Description: Seasonal cycles in the mixing ratios of tropospheric nitrous oxide (N2O) are derived by detrending long-term measurements made at sites across four global surface monitoring networks. The detrended monthly data display large interannual variability, which at some sites challenges the concept of a "mean" seasonal cycle. In the Northern Hemisphere, correlations between polar winter lower stratospheric temperature and detrended N2O data, around the month of the seasonal minimum, provide empirical evidence for a stratospheric influence, which varies in strength from year to year and can explain much of the interannual variability in the surface seasonal cycle. Even at sites where a strong, competing, regional N2O source exists, such as from coastal upwelling at Trinidad Head, California, the stratospheric influence must be understood to interpret the biogeochemical signal in monthly mean data. In the Southern Hemisphere, detrended surface N2O monthly means are correlated with polar spring lower stratospheric temperature in months preceding the N2O minimum, providing empirical evidence for a coherent stratospheric influence in that hemisphere as well, in contrast to some recent atmospheric chemical transport model (ACTM) results. Correlations between the phasing of the surface N2O seasonal cycle in both hemispheres and both polar lower stratospheric temperature and polar vortex break-up date provide additional support for a stratospheric influence. The correlations discussed above are generally more evident in high-frequency in situ data than in data from weekly flask samples. Furthermore, the interannual variability in the N2O seasonal cycle is not always correlated among in situ and flask networks that share common sites, nor do the mean seasonal amplitudes always agree. The importance of abiotic influences such as the stratospheric influx and tropospheric transport on N2O seasonal cycles suggests that, at sites remote from local sources, surface N2O mixing ratio data by themselves are unlikely to provide information about seasonality in surface sources, e.g., for atmospheric inversions, unless the ACTMs employed in the inversions accurately account for these influences. An additional abioitc influence is the seasonal ingassing and outgassing of cooling and warming surface waters, which creates a thermal signal in tropospheric N2O that is of particular importance in the extratropical Southern Hemisphere, where it competes with the biological ocean source signal.

Description: The recent increase of atmospheric methane is investigated by using two atmospheric inversions to quantify the distribution of sources and sinks for the 2006–2008 period, and a process-based model of methane emissions by natural wetland ecosystems. Methane emissions derived from the two inversions are consistent at a global scale: emissions are decreased in 2006 (−7 Tg) and increased in 2007 (+21 Tg) and 2008 (+18 Tg), as compared to the 1999–2006 period. The agreement on the latitudinal partition of the flux anomalies for the two inversions is fair in 2006, good in 2007, and not good in 2008. In 2007, a positive anomaly of tropical emissions is found to be the main contributor to the global emission anomalies (~60–80%) for both inversions, with a dominant share attributed to natural wetlands (~2/3), and a significant contribution from high latitudes (~25%). The wetland ecosystem model produces smaller and more balanced positive emission anomalies between the tropics and the high latitudes for 2006, 2007 and 2008, mainly due to precipitation changes during these years. At a global scale, the agreement between the ecosystem model and the inversions is good in 2008 but not satisfying in 2006 and 2007. Tropical South America and Boreal Eurasia appear to be major contributors to variations in methane emissions consistently in the inversions and the ecosystem model. Finally, changes in OH radicals during 2006–2008 are found to be less than 1% in inversions, with only a small impact on the inferred methane emissions.

Description: Air was sampled from the porous firn layer at the NEEM site in Northern Greenland. We use an ensemble of ten reference tracers of known atmospheric history to characterise the transport properties of the site. By analysing uncertainties in both data and the reference gas atmospheric histories, we can objectively assign weights to each of the gases used for the depth-diffusivity reconstruction. We define an objective root mean square criterion that is minimised in the model tuning procedure. Each tracer constrains the firn profile differently through its unique atmospheric history and free air diffusivity, making our multiple-tracer characterisation method a clear improvement over the commonly used single-tracer tuning. Six firn air transport models are tuned to the NEEM site; all models successfully reproduce the data within a 1σ Gaussian distribution. A comparison between two replicate boreholes drilled 64 m apart shows differences in measured mixing ratio profiles that exceed the experimental error. We find evidence that diffusivity does not vanish completely in the lock-in zone, as is commonly assumed. The ice age- gas age difference (Δage) at the firn-ice transition is calculated to be 182+3−9 yr. We further present the first intercomparison study of firn air models, where we introduce diagnostic scenarios designed to probe specific aspects of the model physics. Our results show that there are major differences in the way the models handle advective transport. Furthermore, diffusive fractionation of isotopes in the firn is poorly constrained by the models, which has consequences for attempts to reconstruct the isotopic composition of trace gases back in time using firn air and ice core records.

Description: Atmospheric radiocarbon (14C) represents an important observational constraint on emissions of fossil-fuel derived carbon into the atmosphere due to the absence of 14C in fossil fuel reservoirs. The high sensitivity and precision that accelerator mass spectrometry (AMS) affords in atmospheric 14C analysis has greatly increased the potential for using such measurements to evaluate bottom-up emissions inventories of fossil fuel CO2 (CO2ff), as well as those for other co-emitted species. Here we use observations of 14CO2 and a series of primary hydrocarbons and combustion tracers from discrete air samples collected between June 2009 and September 2010 at the National Oceanic and Atmospheric Administration Boulder Atmospheric Observatory (BAO; Lat: 40.050° N, Lon: 105.004° W) to derive emission ratios of each species with respect to CO2ff. The BAO tower is situated at the boundary of the Denver metropolitan area to the south and a large industrial and agricultural region to the north and east, making it an ideal location to study the contrasting mix of emissions from the activities in each region. The species considered in this analysis are carbon monoxide (CO), methane (CH4), acetylene (C2H2), benzene (C6H6), and C3–C5 alkanes. We estimate emissions for a subset of these species by using the Vulcan high resolution CO2ff emission data product as a reference. We find that CO is overestimated in the 2008 National Emissions Inventory (NEI08) by a factor of ~2. A close evaluation of the inventory suggests that the ratio of CO emitted per unit fuel burned from on-road gasoline vehicles is likely over-estimated by a factor of 2.5. Using a wind-directional analysis of the data, we find enhanced concentrations of CH4, relative to CO2ff, in air influenced by emissions to the north and east of the BAO tower when compared to air influenced by emissions in the Denver metro region to the south. Along with enhanced CH4, the strongest enhancements of the C3–C5 alkanes are also found in the north and east wind sector, suggesting that both the alkane and CH4 enhancements are sourced from oil and gas fields located to the northeast, though it was not possible to rule out the contribution of non oil and gas CH4 sources.

Description: We describe an assimilation system for atmospheric methane (CH4), CarbonTracker-CH4, and demonstrate the diagnostic value of global or zonally averaged CH4 abundances for evaluating the results. We show that CarbonTracker-CH4 is able to simulate the observed zonal average mole fractions and capture inter-annual variability in emissions quite well at high northern latitudes (53–90° N). In contrast, CarbonTracker-CH4 is less successful in the tropics where there are few observations and therefore misses significant variability and is more influenced by prior flux estimates. CarbonTracker-CH4 estimates of total fluxes at high northern latitudes are about 81 ± 7 Tg CH4 yr−1, about 12 Tg CH4 yr−1 (13%) lower than prior estimates, a result that is consistent with other atmospheric inversions. Emissions from European wetlands are decreased by 30%, a result consistent with previous work by Bergamaschi et al. (2005); however, unlike their results, emissions from wetlands in boreal Eurasia are increased relative to the prior estimate. Although CarbonTracker-CH4 does not estimate an increasing trend in emissions from high northern latitudes for 2000 through 2010, significant inter-annual variability in high northern latitude fluxes is recovered. Exceptionally warm growing season temperatures in the Arctic occurred in 2007, a year that was also anonymously wet. Estimated emissions from natural sources were greater than the decadal average by 4.4 ± 3.8 Tg CH4 yr−1 in 2007. CarbonTracker-CH4 estimates for temperate latitudes are only slightly increased over prior estimates, but about 10 Tg CH4 yr−1 is redistributed from Asia to North America. This difference exceeds the estimated uncertainty for North America (±3.5 Tg CH4 yr−1). We used time invariant prior flux estimates, so for the period from 2000 to 2006, when the growth rate of global atmospheric CH4 was very small, the assimilation does not produce increases in natural or anthropogenic emissions in contrast to bottom-up emission data sets. After 2006, when atmospheric CH4 began its recent increases, CarbonTracker-CH4 allocates some of the increases to anthropogenic emissions at temperate latitudes, and some to tropical wetland emissions. For temperate North America the prior flux increases by about 4 Tg CH4 yr−1 during winter when biogenic emissions are small. Examination of the residuals at some North American observation sites suggests that increased gas and oil exploration may play a role since sites near fossil fuel production are particularly hard for the inversion to fit and the prior flux estimates at these sites are apparently lower and lower over time than what the atmospheric measurements imply. The tropics are not currently well resolved by CarbonTracker-CH4 due to sparse observational coverage and a short assimilation window. However, there is a small uncertainty reduction and posterior emissions are about 18% higher than prior estimates. Most of this increase is allocated to tropical South America rather than being distributed among the global tropics. Our estimates for this source region are about 32 ± 4 Tg CH4 yr−1, in good agreement with the analysis of Melack et al. (2004) who obtained 29 Tg CH4 yr−1 for the most productive region, the Amazon Basin.

Description: A synthesis inversion based on the atmospheric zoom model TM5 is used to derive top-down estimates of CH4 emissions from individual European countries for the year 2001. We employ a model zoom over Europe with 1° × 1° resolution that is two-way nested into the global model domain (with resolution of 6° × 4°. This approach ensures consistent boundary conditions for the zoom domain and thus European top-down estimates consistent with global CH4 observations. The TM5 model, driven by ECMWF analyses, simulates synoptic scale events at most European and global sites fairly well, and the use of high-frequency observations allows exploiting the information content of individual synoptic events. A detailed source attribution is presented for a comprehensive set of 56 monitoring sites, assigning the atmospheric signal to the emissions of individual European countries and larger global regions. The available observational data put significant constraints on emissions from different regions. Within Europe, in particular several Western European countries are well constrained. The inversion results suggest up to 50-90% higher anthropogenic CH4 emissions in 2001 for Germany, France and UK compared to reported UNFCCC values (EEA, 2003). A recent revision of the German inventory, however, resulted in an increase of reported CH4 emissions by 68.5% (EEA, 2004), being now in very good agreement with our top-down estimate. The top-down estimate for Finland is distinctly smaller than the a priori estimate, suggesting much smaller CH4 emissions from Finnish wetlands than derived from the bottom-up inventory. The EU-15 totals are relatively close to UNFCCC values (within 4-30%) and appear very robust for different inversion scenarios.

Description: We use an ensemble Kalman filter (EnKF), together with the GEOS-Chem chemistry transport model, to estimate regional monthly methane (CH4) fluxes for the period June 2009–December 2010 using proxy dry-air column-averaged mole fractions of methane (XCH4) from GOSAT (Greenhouse gases Observing SATellite) and/or NOAA ESRL (Earth System Research Laboratory) and CSIRO GASLAB (Global Atmospheric Sampling Laboratory) CH4 surface mole fraction measurements. Global posterior estimates using GOSAT and/or surface measurements are between 510–516 Tg yr−1, which is less than, though within the uncertainty of, the prior global flux of 529 ± 25 Tg yr−1. We find larger differences between regional prior and posterior fluxes, with the largest changes in monthly emissions (75 Tg yr−1) occurring in Temperate Eurasia. In non-boreal regions the error reductions for inversions using the GOSAT data are at least three times larger (up to 45%) than if only surface data are assimilated, a reflection of the greater spatial coverage of GOSAT, with the two exceptions of latitudes 〉60° associated with a data filter and over Europe where the surface network adequately describes fluxes on our model spatial and temporal grid. We use CarbonTracker and GEOS-Chem XCO2 model output to investigate model error on quantifying proxy GOSAT XCH4 (involving model XCO2) and inferring methane flux estimates from surface mole fraction data and show similar resulting fluxes, with differences reflecting initial differences in the proxy value. Using a series of observing system simulation experiments (OSSEs) we characterize the posterior flux error introduced by non-uniform atmospheric sampling by GOSAT. We show that clear-sky measurements can theoretically reproduce fluxes within 10% of true values, with the exception of tropical regions where, due to a large seasonal cycle in the number of measurements because of clouds and aerosols, fluxes are within 15% of true fluxes. We evaluate our posterior methane fluxes by incorporating them into GEOS-Chem and sampling the model at the location and time of surface CH4 measurements from the AGAGE (Advanced Global Atmospheric Gases Experiment) network and column XCH4 measurements from TCCON (Total Carbon Column Observing Network). The posterior fluxes modestly improve the model agreement with AGAGE and TCCON data relative to prior fluxes, with the correlation coefficients (r2) increasing by a mean of 0.04 (range: −0.17 to 0.23) and the biases decreasing by a mean of 0.4 ppb (range: −8.9 to 8.4 ppb).

Description: We present atmospheric sulfur hexafluoride (SF6) mole fractions and emissions estimates from the 1970s to 2008. Measurements were made of archived air samples starting from 1973 in the Northern Hemisphere and from 1978 in the Southern Hemisphere, using the Advanced Global Atmospheric Gases Experiment (AGAGE) gas chromatographic-mass spectrometric (GC-MS) systems. These measurements were combined with modern high-frequency GC-MS and GC-electron capture detection (ECD) data from AGAGE monitoring sites, to produce a unique 35-year atmospheric record of this potent greenhouse gas. Atmospheric mole fractions were found to have increased by more than an order of magnitude between 1973 and 2008. The 2008 growth rate was the highest recorded, at 0.29 ± 0.02 pmolmol−1 yr−1. A three-dimensional chemical transport model and a minimum variance Bayesian inverse method was used to estimate annual emission rates using the measurements, with a priori estimates from the Emissions Database for Global Atmospheric Research (EDGAR, version 4). Consistent with the mole fraction growth rate maximum, global emissions during 2008 were also the highest in the 1973–2008 period, reaching 7.4 ± 0.6 Gg yr−1 (1-σ uncertainties) and surpassing the previous maximum in 1995. The 2008 values follow an increase in emissions of 48 ± 20% since 2001. A second global inversion which also incorporated National Oceanic and Atmospheric Administration (NOAA) flask measurements and in situ monitoring site data agreed well with the emissions derived using AGAGE measurements alone. By estimating continent-scale emissions using all available AGAGE and NOAA surface measurements covering the period 2004–2008, with no pollution filtering, we find that it is likely that much of the global emissions rise during this five-year period originated primarily from Asian developing countries that do not report detailed, annual emissions to the United Nations Framework Convention on Climate Change (UNFCCC). We also find it likely that SF6 emissions reported to the UNFCCC were underestimated between at least 2004 and 2005.