Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Increase in CFC-11 emissions from eastern China based on atmospheric observations


The recovery of the stratospheric ozone layer relies on the continued decline in the atmospheric concentrations of ozone-depleting gases such as chlorofluorocarbons1. The atmospheric concentration of trichlorofluoromethane (CFC-11), the second-most abundant chlorofluorocarbon, has declined substantially since the mid-1990s2. A recently reported slowdown in the decline of the atmospheric concentration of CFC-11 after 2012, however, suggests that global emissions have increased3,4. A concurrent increase in CFC-11 emissions from eastern Asia contributes to the global emission increase, but the location and magnitude of this regional source are unknown3. Here, using high-frequency atmospheric observations from Gosan, South Korea, and Hateruma, Japan, together with global monitoring data and atmospheric chemical transport model simulations, we investigate regional CFC-11 emissions from eastern Asia. We show that emissions from eastern mainland China are 7.0 ± 3.0 (±1 standard deviation) gigagrams per year higher in 2014–2017 than in 2008–2012, and that the increase in emissions arises primarily around the northeastern provinces of Shandong and Hebei. This increase accounts for a substantial fraction (at least 40 to 60 per cent) of the global rise in CFC-11 emissions. We find no evidence for a significant increase in CFC-11 emissions from any other eastern Asian countries or other regions of the world where there are available data for the detection of regional emissions. The attribution of any remaining fraction of the global CFC-11 emission rise to other regions is limited by the sparsity of long-term measurements of sufficient frequency near potentially emissive regions. Several considerations suggest that the increase in CFC-11 emissions from eastern mainland China is likely to be the result of new production and use, which is inconsistent with the Montreal Protocol agreement to phase out global chlorofluorocarbon production by 2010.

Access options

Rent or Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: Observations of atmospheric CFC-11 at Gosan and Hateruma.
Fig. 2: CFC-11 emissions derived from atmospheric observations.
Fig. 3: Spatial distribution of the derived CFC-11 fluxes in the NAME-HB inversion.

Data availability

Data from remote AGAGE stations and Gosan data are available from the AGAGE website ( Hateruma data are available at World Data Centre for Greenhouse Gases ( NOAA data are taken as the ‘combined set’ data record (a best-estimate record based on flask samples and in situ measurements), available from the NOAA Global Monitoring Division data server (, or more specifically

Code availability

The inversion models of NAME-HB, NAME-InTEM, FLEXPART-MIT and FLEXPART-Empa are used to determine the regional CFC-11 emissions. Enquiries about the model codes should be directed to M.R., A.J.M., S.R. and R.G.P. Licences to use NAME and InTEM are available for research purposes via a request to the UK Met Office or on request from A.J.M. and A.L.R. The code for the NAME-based hierarchical Bayesian inversion (NAME-HB) is available on request from M.R. and L.M.W. The code of the dispersion model FLEXPART is available from The code for the FLEXPART-based Bayesian inversion (FLEXPART-MIT) is available on request from X.F. The inversion code used by Empa is available from or on request from S.H.


  1. 1.

    Harris, N. R. P. & Wuebbles, D. J. in Scientific Assessment of Ozone Depletion: 2014 Ch. 5, 416 (World Meteorological Organization, 2014).

  2. 2.

    Carpenter, L. J. & Reimann, S. in Scientific Assessment of Ozone Depletion: 2014 Ch. 1 (World Meteorological Organization, 2014).

  3. 3.

    Montzka, S. A. et al. An unexpected and persistent increase in global emissions of ozone-depleting CFC-11. Nature 557, 413–417 (2018).

    CAS  ADS  Article  Google Scholar 

  4. 4.

    Prinn, R. G. et al. History of chemically and radiatively important atmospheric gases from the Advanced Global Atmospheric Gases Experiment (AGAGE). Earth Syst. Sci. Data 10, 985–1018 (2018).

    ADS  Article  Google Scholar 

  5. 5.

    Yokouchi, Y. et al. High frequency measurements of HFCs at a remote site in east Asia and their implications for Chinese emissions. Geophys. Res. Lett. 33, L21814 (2006).

  6. 6.

    Hu, L. et al. Considerable contribution of the Montreal Protocol to declining greenhouse gas emissions from the United States. Geophys. Res. Lett. 44, 8075–8083 (2017).

    CAS  ADS  Article  Google Scholar 

  7. 7.

    Fang, X. et al. Changes in emissions of ozone-depleting substances from China due to implementation of the Montreal Protocol. Environ. Sci. Technol. 52, 11359–11366 (2018).

    CAS  ADS  Article  Google Scholar 

  8. 8.

    Palmer, P. I. et al. Eastern Asian emissions of anthropogenic halocarbons deduced from aircraft concentration data. J. Geophys. Res. Atmos. 108, 4753 (2003).

    ADS  Google Scholar 

  9. 9.

    An, X. et al. Estimating emissions of HCFC-22 and CFC-11 in China by atmospheric observations and inverse modeling. Sci. China Chem. 55, 2233–2241 (2012).

    CAS  Article  Google Scholar 

  10. 10.

    Li, S. et al. Emissions of halogenated compounds in East Asia determined from measurements at Jeju Island, Korea. Environ. Sci. Technol. 45, 5668–5675 (2011).

    CAS  ADS  Article  Google Scholar 

  11. 11.

    Kim, J. et al. Regional atmospheric emissions determined from measurements at Jeju Island, Korea: halogenated compounds from China. Geophys. Res. Lett. 37, L12801 (2010).

    ADS  Google Scholar 

  12. 12.

    Vollmer, M. K. et al. Emissions of ozone-depleting halocarbons from China. Geophys. Res. Lett. 36, L15823 (2009).

    ADS  Article  Google Scholar 

  13. 13.

    Saito, T. et al. Extraordinary halocarbon emissions initiated by the 2011 Tohoku earthquake. Geophys. Res. Lett. 42, 2500–2507 (2015).

    CAS  ADS  Article  Google Scholar 

  14. 14.

    Manning, A. J., O’Doherty, S., Jones, A. R., Simmonds, P. G. & Derwent, R. G. Estimating UK methane and nitrous oxide emissions from 1990 to 2007 using an inversion modeling approach. J. Geophys. Res. 116, D02305 (2011).

    ADS  Article  Google Scholar 

  15. 15.

    Stohl, A., Forster, C., Frank, A., Seibert, P. & Wotawa, G. Technical note: the Lagrangian particle dispersion model FLEXPART version 6.2. Atmos. Chem. Phys. 5, 2461–2474 (2005).

    CAS  ADS  Article  Google Scholar 

  16. 16.

    Lunt, M. F., Rigby, M., Ganesan, A. L. & Manning, A. J. Estimation of trace gas fluxes with objectively determined basis functions using reversible-jump Markov chain Monte Carlo. Geosci. Model Dev. 9, 3213–3229 (2016).

    CAS  ADS  Article  Google Scholar 

  17. 17.

    Arnold, T. et al. Inverse modelling of CF4 and NF3 emissions in East Asia. Atmos. Chem. Phys. 18, 13305–13320 (2018).

    CAS  ADS  Article  Google Scholar 

  18. 18.

    Fang, X. et al. Rapid increase in ozone-depleting chloroform emissions from China. Nat. Geosci. 12, 89–93 (2019).

    CAS  ADS  Article  Google Scholar 

  19. 19.

    Henne, S. et al. Validation of the Swiss methane emission inventory by atmospheric observations and inverse modelling. Atmos. Chem. Phys. 16, 3683–3710 (2016).

    CAS  ADS  Article  Google Scholar 

  20. 20.

    Ko, M. K. W., Newman, P. A., Reimann, S. & Strahan, S. E. in SPARC Report No. 6: Lifetimes of Stratospheric Ozone-Depleting Substances, Their Replacements, and Related Species (eds Ko, M. et al.) (Stratospheric Processes And their Role in Climate, 2013).

  21. 21.

    Duan, H. et al. Chilling prospect: climate change effects of mismanaged refrigerants in China. Environ. Sci. Technol. 52, 6350–6356 (2018).

    CAS  ADS  Article  Google Scholar 

  22. 22.

    The Intergovernmental Panel on Climate Change. Safeguarding the Ozone Layer and the Global Climate System: Issues Related to Hydrofluorocarbons and Perfluorocarbons (SROC) (IPCC/TEAP, Cambridge Univ. Press, 2005).

  23. 23.

    Lunt, M. F. et al. Continued emissions of the ozone-depleting substance carbon tetrachloride from eastern Asia. Geophys. Res. Lett. 45, 11423–11430 (2018).

    CAS  ADS  Article  Google Scholar 

  24. 24.

    Miller, B. R. et al. Medusa: a sample preconcentration and GC/MS detector system for in situ measurements of atmospheric trace halocarbons, hydrocarbons, and sulfur compounds. Anal. Chem. 80, 1536–1545 (2008).

    CAS  Article  Google Scholar 

  25. 25.

    Enomoto, T., Yokouchi, Y., Izumi, K. & Inagaki, T. Development of an analytical method for atmospheric halocarbons and its application to airborne observation. J. Jpn. Soc. Atmospheric Environ. 40, 1–8 (2005).

    CAS  Google Scholar 

  26. 26.

    Cullen, M. J. P. The Unified Forecast/Climate Model. Meteorol. Mag. 122, 81–94 (1993).

    Google Scholar 

  27. 27.

    Chipperfield, M. P. et al. Multimodel estimates of atmospheric lifetimes of long-lived ozone-depleting substances: present and future. J. Geophys. Res. 119, 2555–2573 (2014).

    ADS  Article  Google Scholar 

  28. 28.

    Ganesan, A. L. et al. Characterization of uncertainties in atmospheric trace gas inversions using hierarchical Bayesian methods. Atmos. Chem. Phys. 14, 3855–3864 (2014).

    ADS  Article  Google Scholar 

  29. 29.

    Green, P. J. Reversible jump Markov chain Monte Carlo computation and Bayesian model determination. Biometrika 82, 711–732 (1995).

    MathSciNet  Article  Google Scholar 

  30. 30.

    Stohl, A. et al. An analytical inversion method for determining regional and global emissions of greenhouse gases: sensitivity studies and application to halocarbons. Atmos. Chem. Phys. 9, 1597–1620 (2009).

    CAS  ADS  Article  Google Scholar 

  31. 31.

    O’Doherty, S. et al. In situ chloroform measurements at Advanced Global Atmospheric Gases Experiment atmospheric research stations from 1994 to 1998. J. Geophys. Res. 106, 20429–20444 (2001).

    ADS  Article  Google Scholar 

  32. 32.

    Rigby, M. et al. Re-evaluation of the lifetimes of the major CFCs and CH3CCl3 using atmospheric trends. Atmos. Chem. Phys. 13, 2691–2702 (2013).

    ADS  Article  Google Scholar 

  33. 33.

    Rigby, M. et al. Recent and future trends in synthetic greenhouse gas radiative forcing. Geophys. Res. Lett. 41, 2623–2630 (2014).

    CAS  ADS  Article  Google Scholar 

Download references


We are indebted to the site operators who oversee the day-to-day running of the AGAGE, NIES and NOAA stations. We particularly thank the NASA Upper Atmosphere Research Program for its continuing support of AGAGE, including providing modelling, field station and instrumentation support though grant NNX16AC98G to MIT, and supporting the overall experimental programme through grants NNX16AC96G and NNX16AC97G to SIO. Observations at Cape Grim are supported largely by the Australian Bureau of Meteorology, CSIRO, the Australian Department of the Environment and Energy (DoEE) and Refrigerant Reclaim Australia (RRA). Mace Head, Ireland, is supported by the Department for Business, Energy & Industrial Strategy (BEIS, UK, formerly the Department of Energy and Climate Change (DECC)) contract 1028/06/2015 to the University of Bristol and the UK Meteorological Office. Ragged Point, Barbados is supported by the National Oceanic and Atmospheric Administration (NOAA, USA), contract RA-133-R15-CN-0008 to the University of Bristol. M.R., L.M.W., M.F.L. and R.L.T. were supported by Natural Environment Research Council grants NE/I021365/1, NE/I027282/1, NE/M014851/1, NE/L013088/1 and NE/N016548/1. A.L.G. is supported by NERC Independent Research Fellowship NE/L010992/1. S.P., T.L., S.L., M.-K.P. and K.-R.K. were supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (no. NRF-2016R1A2B2010663). Observations at Hateruma were partly supported by the Ministry of the Environment of Japan. S.A.M., B.D.H. and G.S.D. are indebted to C. Siso, D. Mondeel, J. D. Nance, F. Moore, J.W. Elkins, B. Vasel, C. Schultz, R. Schnell and J. H. Butler for discussions and assistance with the NOAA measurements considered here, which were made possible in part with support from the NOAA Climate Program Office’s AC4 program.

Reviewer information

Nature thanks Andreas Stohl, Guus Velders and the other anonymous reviewer(s) for their contribution to the peer review of this work.

Author information




S.P., T.S., G.S.D., P.J.F., B.D.H., C.M.H., J.K., K.-R.K., P.B.K., T.L., S.L., S.A.M., J.M., S.O., M.-K.P., P.K.S., P.S., R.F.W., Y.Y. and D.Y. contributed observational data. L.M.W., A.L.R., X.F. and S.H. carried out atmospheric model simulations and inverse analysis with support from R.L.T., M.F.L. and A.L.G., and under the supervision of M.R., A.J.M., S.R. and R.G.P. All authors wrote the manuscript.

Corresponding author

Correspondence to S. Park.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 CFC-11 mole fractions measured at different sites in the AGAGE network and at affiliated sites.

Box plots indicate the 25th and 75th of the individual measurement data (approximately 2-hourly), with the median shown as a horizontal line within each box. The whiskers show the 10th and 90th percentiles. The lower percentiles are typically representative of baseline (‘unpolluted’) mole fractions, and the difference between the lower and higher percentiles indicates the magnitude of above-baseline events due to the interception of air masses containing recently emitted CFC-11.

Extended Data Fig. 2 NAME-derived mean sensitivities of atmospheric mole fractions at Gosan and Hateruma to potential emission.

This mean sensitivity map was derived for sampling events at these two sites during 2008 to 2017. Black triangle and circle indicate the Gosan and Hateruma stations, respectively. Thin grey lines show the boundaries of the provinces within the region we denote eastern mainland China. This region contains the provinces Anhui, Beijing, Hebei, Jiangsu, Liaoning, Shandong, Shanghai, Tianjin and Zhejiang. The west of the grey line transecting Japan, containing the regions Chūgoku, Kansai, Kyūshū & Okinawa, and Shikoku, is the region we denote western Japan.

Extended Data Fig. 3 Simulated and observed CFC-11 mole fractions at Gosan and Hateruma.

Left, a comparison of the simulated CFC-11 mole fractions from the four different inversion analyses and those that were measured at Gosan and Hateruma. Right, residuals between the simulated and observed mole fractions (data minus model). Shading denotes 1 s.d. of the model-data mismatch uncertainties assumed in the inversions. Simulated mole fractions are derived from the a posteriori emissions. For the NAME-InTEM, FLEXPART-EMPA, NAME-HB and FLEXPART-MIT inversions, 2-hourly, 3-hourly, 24-hourly and 24-hourly averaging was applied to the model and data, respectively.

Extended Data Fig. 4 Derived emissions from different inverse models considering subsets of measurement data.

Shown are emission estimates derived using the Gosan CFC-11 measurements (green dashed line), the Hateruma CFC-11 measurements (red dotted line), or both records in the inversion analysis (pink solid line). ad, Estimates shown use the NAME-HB (a), NAME-InTEM (b), FLEXPART-MIT (c) and FLEXPART-Empa (d) inversion techniques.

Extended Data Fig. 5 Derived emissions for eastern mainland China for different a priori emission magnitudes and different spatial distributions.

The inversions are as described in Methods (‘a priori’, pink solid line), with a priori emissions twice as high (‘2× a priori’, green dashed line), or the same magnitude but distributed in space according to population density (‘population weighted’, red dotted line). ad, Estimates shown use the NAME-HB (a), NAME-InTEM (b), FLEXPART-MIT (c) and FLEXPART-Empa (d) inversion methods.

Extended Data Fig. 6 Top-down CFC-11 emissions estimates in eastern Asia.

ac, Emission estimates are shown for western Japan (a), South Korea (b) and North Korea (c) using the NAME-HB (yellow lines), NAME-InTEM (blue lines), FLEXPART-MIT (pink lines) and FLEXPART-Empa (grey lines) inverse frameworks described in the Methods. See Extended Data Fig. 2 for the definition of the western Japan region. All lines and symbols are the a posteriori mean, and shading denotes s.d. uncertainty.

Extended Data Fig. 7 Maps of mean CFC-11 emission fluxes from different models and time periods, and the differences between time periods.

ai, CFC-11 emissions for NAME-InTEM (ac), FLEXPART-MIT (df) and FLEXPART-Empa (gi) using the inversion framework described in the Methods, and similar to Fig. 3. a, d, g, Average spatial emissions are shown for the period 2008–2012. b, e, h, Average spatial emissions are shown for the period 2014–2017. c, f, i, The difference in emissions from the 2008–2012 period to the 2013–2017 period, using NAME-InTEM, FLEXPART-MIT and FLEXPART-Empa, respectively.

Extended Data Fig. 8 The implied bank release fraction of CFC-11 from eastern mainland China.

The implied release fraction of banks of CFC-11 from the eastern mainland China region defined in Extended Data Fig. 2, assuming that there is no non-reported production after 2008. This calculation assumes that the 2008 emissions estimates are 5% of the bank size (estimated to be towards the upper limit of expected global bank fractional release rate3). The coloured lines show the implied bank release fraction for each year after 2008, required to sustain the mean a posteriori emissions from the NAME-HB (yellow), NAME-InTEM (blue), FLEXPART-MIT (pink) and FLEXPART-Empa (grey) inversions. The black line shows the bank release fraction from a bottom-up study for the whole of China7.

Extended Data Fig. 9 CFC-12 mole fractions at the Gosan measurement site.

Measured mole fractions of CFC-12 at Gosan, South Korea, from 2008 to 2017. Box plots are as defined in Extended Data Fig. 1.

Supplementary information

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Rigby, M., Park, S., Saito, T. et al. Increase in CFC-11 emissions from eastern China based on atmospheric observations. Nature 569, 546–550 (2019).

Download citation

Further reading


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links