ALBERT

Description: We present a method to infer spatially and spatio-temporally correlated emissions of greenhouse gases from atmospheric measurements and a chemical transport model. The method allows fast computation of spatial emissions using a hierarchical Bayesian framework as an alternative to Markov chain Monte Carlo algorithms. The spatial emissions follow a Gaussian process with a Matérn correlation structure which can be represented by a Gaussian Markov random field through a stochastic partial differential equation approach. The inference is based on an integrated nested Laplacian approximation (INLA) for hierarchical models with Gaussian latent fields. Combining an autoregressive temporal correlation and the Matérn field provides a full spatio-temporal correlation structure. We first demonstrate the method on a synthetic data example and follow this using a well-studied test case of inferring UK methane emissions from tall tower measurements of atmospheric mole fraction. Results from these two test cases show that this method can accurately estimate regional greenhouse gas emissions, accounting for spatio-temporal uncertainties that have traditionally been neglected in atmospheric inverse modelling.

Description: We report a 40-year history of SF6 atmospheric mole fractions measured at the Advanced Global Atmospheric Gases Experiment (AGAGE) monitoring sites, combined with archived air samples, to determine emission estimates from 1978 to 2018. Previously we reported a global emission rate of 7.3±0.6 Gg yr−1 in 2008 and over the past decade emissions have continued to increase by about 24 % to 9.04±0.35 Gg yr−1 in 2018. We show that changing patterns in SF6 consumption from developed (Kyoto Protocol Annex-1) to developing countries (non-Annex-1) and the rapid global expansion of the electric power industry, mainly in Asia, have increased the demand for SF6-insulated switchgear, circuit breakers, and transformers. The large bank of SF6 sequestered in this electrical equipment provides a substantial source of emissions from maintenance, replacement, and continuous leakage. Other emissive sources of SF6 occur from the magnesium, aluminium, and electronics industries as well as more minor industrial applications. More recently, reported emissions, including those from electrical equipment and metal industries, primarily in the Annex-1 countries, have declined steadily through substitution of alternative blanketing gases and technological improvements in less emissive equipment and more efficient industrial practices. Nevertheless, there are still demands for SF6 in Annex-1 countries due to economic growth, as well as continuing emissions from older equipment and additional emissions from newly installed SF6-insulated electrical equipment, although at low emission rates. In addition, in the non-Annex-1 countries, SF6 emissions have increased due to an expansion in the growth of the electrical power, metal, and electronics industries to support their continuing development. There is an annual difference of 2.5–5 Gg yr−1 (1990–2018) between our modelled top-down emissions and the UNFCCC-reported bottom-up emissions (United Nations Framework Convention on Climate Change), which we attempt to reconcile through analysis of the potential contribution of emissions from the various industrial applications which use SF6. We also investigate regional emissions in East Asia (China, S. Korea) and western Europe and their respective contributions to the global atmospheric SF6 inventory. On an average annual basis, our estimated emissions from the whole of China are approximately 10 times greater than emissions from western Europe. In 2018, our modelled Chinese and western European emissions accounted for ∼36 % and 3.1 %, respectively, of our global SF6 emissions estimate.

Description: High frequency measurements of trifluoromethane (HFC-23, CHF3), a potent hydrofluorocarbon greenhouse gas, largely emitted to the atmosphere as a by-product of the production of the hydrochlorofluorocarbon HCFC-22 (CHClF2), at five core stations of the Advanced Global Atmospheric Gases Experiment (AGAGE) network, combined with measurements on firn air, old Northern Hemisphere air samples and Cape Grim Air Archive (CGAA) air samples, are used to explore the current and historic changes in the atmospheric abundance of HFC-23. These measurements are used in combination with the AGAGE 2-D atmospheric 12-box model and a Bayesian inversion methodology to determine model atmospheric mole fractions and the history of global HFC-23 emissions. The global modelled annual mole fraction of HFC-23 in the background atmosphere was 28.9 ± 0.6 pmol mol−1 at the end of 2016, representing a 28 % increase from 22.6 ± 0.4 pmol mol−1 in 2009. Over the same time frame, the modelled mole fraction of HCFC-22 increased by 19 % from 199 ± 2 to 237 ± 2 pmol mol−1. However, unlike HFC-23, the annual average HCFC-22 growth rate slowed from 2009 to 2016 at an annual average rate of −0.5 pmol mol−1 yr−2. This slowing atmospheric growth is consistent with HCFC-22 moving from dispersive (high fractional emissions) to feedstock (low fractional emissions) uses, with HFC-23 emissions remaining as a consequence of incomplete mitigation from all HCFC-22 production.Our results demonstrate that, following a minimum in HFC-23 global emissions in 2009 of 9.6 ± 0.6, emissions increased to a maximum in 2014 of 14.5 ± 0.6 Gg yr−1 and then declined to 12.7 ± 0.6 Gg yr−1 (157 Mt CO2 eq. yr−1) in 2016. The 2009 emissions minimum is consistent with estimates based on national reports and is likely a response to the implementation of the Clean Development Mechanism (CDM) to mitigate HFC-23 emissions by incineration in developing (non-Annex 1) countries under the Kyoto Protocol. Our derived cumulative emissions of HFC-23 during 2010–2016 were 89 ± 2 Gg (1.1 ± 0.2 Gt CO2 eq.), which led to an increase in radiative forcing of 1.0 ± 0.1 mW m−2 over the same period. Although the CDM had reduced global HFC-23 emissions, it cannot now offset the higher emissions from increasing HCFC-22 production in non-Annex 1 countries, as the CDM was closed to new entrants in 2009. We also find that the cumulative European HFC-23 emissions from 2010 to 2016 were  ∼  1.3 Gg, corresponding to just 1.5 % of cumulative global HFC-23 emissions over this same period. The majority of the increase in global HFC-23 emissions since 2010 is attributed to a delay in the adoption of mitigation technologies, predominantly in China and East Asia. However, a reduction in emissions is anticipated, when the Kigali 2016 amendment to the Montreal Protocol, requiring HCFC and HFC production facilities to introduce destruction of HFC-23, is fully implemented.

Description: We present a method to derive atmospheric-observation-based estimates of carbon dioxide (CO2) fluxes at the national scale, demonstrated using data from a network of surface tall-tower sites across the UK and Ireland over the period 2013–2014. The inversion is carried out using simulations from a Lagrangian chemical transport model and an innovative hierarchical Bayesian Markov chain Monte Carlo (MCMC) framework, which addresses some of the traditional problems faced by inverse modelling studies, such as subjectivity in the specification of model and prior uncertainties. Biospheric fluxes related to gross primary productivity and terrestrial ecosystem respiration are solved separately in the inversion and then combined a posteriori to determine net ecosystem exchange of CO2. Two different models, Data Assimilation Linked Ecosystem Carbon (DALEC) and Joint UK Land Environment Simulator (JULES), provide prior estimates for these fluxes. We carry out separate inversions to assess the impact of these different priors on the posterior flux estimates and evaluate the differences between the prior and posterior estimates in terms of missing model components. The Numerical Atmospheric dispersion Modelling Environment (NAME) is used to relate fluxes to the measurements taken across the regional network. Posterior CO2 estimates from the two inversions agree within estimated uncertainties, despite large differences in the prior fluxes from the different models. With our method, averaging results from 2013 and 2014, we find a total annual net biospheric flux for the UK of 8±79 Tg CO2 yr−1 (DALEC prior) and 64±85 Tg CO2 yr−1 (JULES prior), where negative values represent an uptake of CO2. These biospheric CO2 estimates show that annual UK biospheric sources and sinks are roughly in balance. These annual mean estimates consistently indicate a greater net release of CO2 than the prior estimates, which show much more pronounced uptake in summer months.

Description: Atmospheric inverse modelling has become an increasingly useful tool for evaluating emissions of greenhouse gases including methane, nitrous oxide, and synthetic gases such as hydrofluorocarbons (HFCs). Atmospheric inversions for emissions of CO2 from fossil fuel combustion (ffCO2) are currently being developed. The aim of this paper is to investigate potential errors and uncertainties related to the spatial and temporal prior representation of emissions and modelled atmospheric transport for the inversion of ffCO2 emissions in the US state of California. We perform simulation experiments based on a network of ground-based observations of CO2 concentration and radiocarbon in CO2 (a tracer of ffCO2), combining prior (bottom-up) emission models and transport models currently used in many atmospheric studies. The potential effect of errors in the spatial and temporal distribution of prior emission estimates is investigated in experiments by using perturbed versions of the emission estimates used to create the pseudo-data. The potential effect of transport error was investigated by using three different atmospheric transport models for the prior and pseudo-data simulations. We find that the magnitude of biases in posterior total state emissions arising from errors in the spatial and temporal distribution in prior emissions in these experiments are 1 %–15 % of posterior total state emissions and are generally smaller than the 2σ uncertainty in posterior emissions. Transport error in these experiments introduces biases of −10 % to +6 % into posterior total state emissions. Our results indicate that uncertainties in posterior total state ffCO2 estimates arising from the choice of prior emissions or atmospheric transport model are on the order of 15 % or less for the ground-based network in California we consider. We highlight the need for temporal variations to be included in prior emissions and for continuing efforts to evaluate and improve the representation of atmospheric transport for regional ffCO2 inversions.

Description: Based on observations of the chlorofluorocarbons CFC-13 (chlorotrifluoromethane), ΣCFC-114 (combined measurement of both isomers of dichlorotetrafluoroethane), and CFC-115 (chloropentafluoroethane) in atmospheric and firn samples, we reconstruct records of their tropospheric histories spanning nearly 8 decades. These compounds were measured in polar firn air samples, in ambient air archived in canisters, and in situ at the AGAGE (Advanced Global Atmospheric Gases Experiment) network and affiliated sites. Global emissions to the atmosphere are derived from these observations using an inversion based on a 12-box atmospheric transport model. For CFC-13, we provide the first comprehensive global analysis. This compound increased monotonically from its first appearance in the atmosphere in the late 1950s to a mean global abundance of 3.18 ppt (dry-air mole fraction in parts per trillion, pmol mol−1) in 2016. Its growth rate has decreased since the mid-1980s but has remained at a surprisingly high mean level of 0.02 ppt yr−1 since 2000, resulting in a continuing growth of CFC-13 in the atmosphere. ΣCFC-114 increased from its appearance in the 1950s to a maximum of 16.6 ppt in the early 2000s and has since slightly declined to 16.3 ppt in 2016. CFC-115 increased monotonically from its first appearance in the 1960s and reached a global mean mole fraction of 8.49 ppt in 2016. Growth rates of all three compounds over the past years are significantly larger than would be expected from zero emissions. Under the assumption of unchanging lifetimes and atmospheric transport patterns, we derive global emissions from our measurements, which have remained unexpectedly high in recent years: mean yearly emissions for the last decade (2007–2016) of CFC-13 are at 0.48 ± 0.15 kt yr−1 (〉 15 % of past peak emissions), of ΣCFC-114 at 1.90 ± 0.84 kt yr−1 (∼ 10 % of peak emissions), and of CFC-115 at 0.80 ± 0.50 kt yr−1 (〉 5 % of peak emissions). Mean yearly emissions of CFC-115 for 2015–2016 are 1.14 ± 0.50 kt yr−1 and have doubled compared to the 2007–2010 minimum. We find CFC-13 emissions from aluminum smelters but if extrapolated to global emissions, they cannot account for the lingering global emissions determined from the atmospheric observations. We find impurities of CFC-115 in the refrigerant HFC-125 (CHF2CF3) but if extrapolated to global emissions, they can neither account for the lingering global CFC-115 emissions determined from the atmospheric observations nor for their recent increases. We also conduct regional inversions for the years 2012–2016 for the northeastern Asian area using observations from the Korean AGAGE site at Gosan and find significant emissions for ΣCFC-114 and CFC-115, suggesting that a large fraction of their global emissions currently occur in northeastern Asia and more specifically on the Chinese mainland.

Description: High-frequency, in situ global observations of HCFC-22 (CHClF2), HCFC-141b (CH3CCl2F), HCFC-142b (CH3CClF2) and HCFC-124 (CHClFCF3) and their main HFC replacements, HFC-134a (CH2FCF3), HFC-125 (CHF2CF3), HFC-143a (CH3CF3) and HFC-32 (CH2F2), have been used to determine their changing global growth rates and emissions in response to the Montreal Protocol and its recent amendments. Global mean mole fractions of HCFC-22, -141b, and -142b have increased throughout the observation period, reaching 234, 24.3 and 22.4 pmol mol−1, respectively, in 2015. HCFC-124 reached a maximum global mean mole fraction of 1.48 pmol mol−1 in 2007 and has since declined by 23 % to 1.14 pmol mol−1 in 2015. The HFCs all show increasing global mean mole fractions. In 2015 the global mean mole fractions (pmol mol−1) were 83.3 (HFC-134a), 18.4 (HFC-125), 17.7 (HFC-143a) and 10.5 (HFC-32). The 2007 adjustment to the Montreal Protocol required the accelerated phase-out of emissive uses of HCFCs with global production and consumption capped in 2013 to mitigate their environmental impact as both ozone-depleting substances and important greenhouse gases. We find that this change has coincided with a stabilisation, or moderate reduction, in global emissions of the four HCFCs with aggregated global emissions in 2015 of 449 ± 75 Gg yr−1, in CO2 equivalent units (CO2 eq.) 0.76 ± 0.1 Gt yr−1, compared with 483 ± 70 Gg yr−1 (0.82 ± 0.1 Gt yr−1 CO2 eq.) in 2010 (uncertainties are 1σ throughout this paper). About 79 % of the total HCFC atmospheric burden in 2015 is HCFC-22, where global emissions appear to have been relatively similar since 2011, in spite of the 2013 cap on emissive uses. We attribute this to a probable increase in production and consumption of HCFC-22 in Montreal Protocol Article 5 (developing) countries and the continuing release of HCFC-22 from the large banks which dominate HCFC global emissions. Conversely, the four HFCs all show increasing mole fraction growth rates with aggregated global HFC emissions of 327 ± 70 Gg yr−1 (0.65 ± 0.12 Gt yr−1 CO2 eq.) in 2015 compared to 240 ± 50 Gg yr−1 (0.47 ± 0.08 Gt yr−1 CO2 eq.) in 2010. We also note that emissions of HFC-125 and HFC-32 appear to have increased more rapidly averaged over the 5-year period 2011–2015, compared to 2006–2010. As noted by Lunt et al. (2015) this may reflect a change to refrigerant blends, such as R-410A, which contain HFC-32 and -125 as a 50 : 50 blend.

Description: We describe the motivation, design, and execution of the Greenhouse gAs Uk and Global Emissions (GAUGE) project. The overarching scientific objective of GAUGE was to use atmospheric data to estimate the magnitude, distribution, and uncertainty of the UK greenhouse gas (GHG, defined here as CO2, CH4, and N2O) budget, 2013–2015. To address this objective, we established a multi-year and interlinked measurement and data analysis programme, building on an established tall-tower GHG measurement network. The calibrated measurement network comprises ground-based, airborne, ship-borne, balloon-borne, and space-borne GHG sensors. Our choice of measurement technologies and measurement locations reflects the heterogeneity of UK GHG sources, which range from small point sources such as landfills to large, diffuse sources such as agriculture. Atmospheric mole fraction data collected at the tall towers and on the ships provide information on sub-continental fluxes, representing the backbone to the GAUGE network. Additional spatial and temporal details of GHG fluxes over East Anglia were inferred from data collected by a regional network. Data collected during aircraft flights were used to study the transport of GHGs on local and regional scales. We purposely integrated new sensor and platform technologies into the GAUGE network, allowing us to lay the foundations of a strengthened UK capability to verify national GHG emissions beyond the project lifetime. For example, current satellites provide sparse and seasonally uneven sampling over the UK mainly because of its geographical size and cloud cover. This situation will improve with new and future satellite instruments, e.g. measurements of CH4 from the TROPOspheric Monitoring Instrument (TROPOMI) aboard Sentinel-5P. We use global, nested, and regional atmospheric transport models and inverse methods to infer geographically resolved CO2 and CH4 fluxes. This multi-model approach allows us to study model spread in a posteriori flux estimates. These models are used to determine the relative importance of different measurements to infer the UK GHG budget. Attributing observed GHG variations to specific sources is a major challenge. Within a UK-wide spatial context we used two approaches: (1) Δ14CO2 and other relevant isotopologues (e.g. δ13CCH4) from collected air samples to quantify the contribution from fossil fuel combustion and other sources, and (2) geographical separation of individual sources, e.g. agriculture, using a high-density measurement network. Neither of these represents a definitive approach, but they will provide invaluable information about GHG source attribution when they are adopted as part of a more comprehensive, long-term national GHG measurement programme. We also conducted a number of case studies, including an instrumented landfill experiment that provided a test bed for new technologies and flux estimation methods. We anticipate that results from the GAUGE project will help inform other countries on how to use atmospheric data to quantify their nationally determined contributions to the Paris Agreement.

Description: As the second most populous country and third fastest growing economy, India has emerged as a global economic power. As such, its emissions of greenhouse and ozone-depleting gases are of global significance. However, unlike neighbouring China, the Indian sub-continent is very poorly monitored by atmospheric measurement networks. India's halocarbon emissions, here defined as chlorofluorocarbons (CFCs), hydrochlorofluorocarbons (HCFCs), hydrofluorocarbons (HFCs) and chlorocarbons, are not well-known. Previous measurements from the region have been obtained at observatories many hundreds of kilometres from source regions, or at high altitudes, limiting their value for the estimation of regional emission rates. Given the projected rapid growth in demand for refrigerants and solvents in India, emission estimates of these halocarbons are urgently needed to provide a benchmark against which future changes can be evaluated. In this study, we report atmospheric-measurement-derived halocarbon emissions from India. With the exception of dichloromethane, these top-down estimates are the first for India's halocarbons. Air samples were collected at low altitude during an aircraft campaign in June and July 2016, and emissions were derived from measurements of these samples using an inverse modelling framework. These results were evaluated to assess India's progress in phasing out ozone-depleting substances under the Montreal Protocol. India's combined CFC emissions are estimated to be 54 (27–86) Tg CO2 eq. yr−1 (5th and 95th confidence intervals are shown in parentheses). HCFC-22 emissions of 7.8 (6.0–9.9) Gg yr−1 are of similar magnitude to emissions of HFC-134a (8.2 (6.1–10.7) Gg yr−1). We estimate India's HFC-23 emissions to be 1.2 (0.9–1.5) Gg yr−1, and our results are consistent with resumed venting of HFC-23 by HCFC-22 manufacturers following the discontinuation of funding for abatement under the Clean Development Mechanism. We report small emissions of HFC-32 and HFC-143a and provide evidence to suggest that HFC-32 emissions were primarily due to fugitive emissions during manufacturing processes. A lack of significant correlation among HFC species and the small emissions derived for HFC-32 and HFC-143a indicate that in 2016, India's use of refrigerant blends R-410A, R-404A and R-507A was limited, despite extensive consumption elsewhere in the world. We also estimate emissions of the regulated chlorocarbons carbon tetrachloride and methyl chloroform from northern and central India to be 2.3 (1.5–3.4) and 0.07 (0.04–0.10) Gg yr−1 respectively. While the Montreal Protocol has been successful in reducing emissions of many ozone-depleting substances, growth in the global emission rates of the unregulated very short-lived substances poses an ongoing threat to the recovery of the ozone layer. Emissions of dichloromethane are found to be 96.5 (77.8–115.6) Gg yr−1, and our estimate suggests a 5-fold increase in emissions since the last estimate derived from atmospheric data in 2008. We estimate perchloroethene emissions from India and chloroform emissions from northern–central India to be 2.9 (2.5–3.3) and 32.2 (28.3–37.1) Gg yr−1 respectively. Given the rapid growth of India's economy and the likely increase in demand for halocarbons such as HFCs, the implementation of long-term atmospheric monitoring in the region is urgently required. Our results provide a benchmark against which future changes to India's halocarbon emissions may be evaluated.

Description: High frequency, in situ global observations of HCFC-22 (CHClF2), HCFC-141b (CH3CCl2F), HCFC-142b (CH3CClF2) and HCFC-124 (CHClFCF3) and their main HFC replacements HFC-134a (CH2FCF3), HFC-125 (CHF2CF3), HFC-143a (CH3CF3), and HFC-32 (CH2F2) have been used to determine their changing global growth rates and emissions in response to the Montreal Protocol and its recent amendments. The 2007 adjustment to the Montreal Protocol required the accelerated phase-out of HCFCs with global production and consumption capped in 2013, to mitigate their environmental impact as both ozone depleting substances and important greenhouse gases. We find that this change has coincided with a reduction in global emissions of the four HCFCs with aggregated global emissions in 2015 of 444 ± 75 Gg/yr, in CO2 equivalent units (CO2 e) 0.75 ± 0.1 Gt/yr, compared with 483 ± 70 Gg/yr (0.82 ± 0.1 Gt/yr CO2 e) in 2010. (All quoted uncertainties in this paper are 1 sigma). About 80 % of the total HCFC atmospheric burden in 2015 is HCFC-22, where global HCFC emissions appear to have been relatively constant in spite of the 2013 cap on global production and consumption. We attribute this to a probable increase in production and consumption of HCFC-22 in Montreal Protocol Article 5 (developing) countries and the continuing release of HCFC-22 from the large banks which dominate HCFC global emissions. Conversely, the four HFCs all show increasing annual growth rates with aggregated global HFCs emissions in 2015 of 329 ± 70 Gg/yr (0.65 ± 0.12 Gt/yr CO2 e) compared to 2010 with 240 ± 50 Gg/yr (0.47 ± 0.08 Gt/yr CO2 e). As HCFCs are replaced by HFCs we investigate the impact of the shift to refrigerant blends which have lower global warming potentials (GWPs). We also note that emissions of HFC-125 and HFC-32 appear to have increased more rapidly during the 2011–2015 5-yr period compared to 2006–2010.