Journal cover Journal topic
Atmospheric Chemistry and Physics An interactive open-access journal of the European Geosciences Union
Journal topic
ACP | Articles | Volume 20, issue 8
Atmos. Chem. Phys., 20, 4787–4807, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
Atmos. Chem. Phys., 20, 4787–4807, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 24 Apr 2020

Research article | 24 Apr 2020

Trends and emissions of six perfluorocarbons in the Northern Hemisphere and Southern Hemisphere

Trends and emissions of six perfluorocarbons in the Northern Hemisphere and Southern Hemisphere
Elise S. Droste1, Karina E. Adcock1, Matthew J. Ashfold2, Charles Chou3, Zoë Fleming4,a, Paul J. Fraser5, Lauren J. Gooch1, Andrew J. Hind1, Ray L. Langenfelds5, Emma Leedham Elvidge1, Norfazrin Mohd Hanif1,6, Simon O'Doherty7, David E. Oram1,8, Chang-Feng Ou-Yang9, Marios Panagi4, Claire E. Reeves1, William T. Sturges1, and Johannes C. Laube1,10 Elise S. Droste et al.
  • 1Centre for Ocean and Atmospheric Sciences, School of Environmental Sciences, University of East Anglia, Norwich, NR4 7TJ, UK
  • 2School of Environmental and Geographical Sciences, University of Nottingham Malaysia, 43500 Semenyih, Malaysia
  • 3Research Center for Environmental Changes, Academia Sinica, Taipei 11529, Taiwan
  • 4National Centre for Atmospheric Science (NCAS), Department of Chemistry, University of Leicester, Leicester, UK
  • 5Commonwealth Scientific and Industrial Research Organisation, Oceans and Atmosphere, Climate Science Centre, Aspendale, Australia
  • 6Department of Earth Sciences and Environment, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia
  • 7Department of Chemistry, University of Bristol, Bristol, UK
  • 8National Centre for Atmospheric Science, School of Environmental Sciences, University of East Anglia, Norwich, NR4 7TJ, UK
  • 9Department of Atmospheric Sciences, National Central University, Taoyuan, Taiwan
  • 10Institute of Energy and Climate Research – Stratosphere (IEK-7), Forschungszentrum Jülich GmbHJ, Jülich, Germany
  • anow at: Center for Climate and Resilience Research (CR2), University of Chile, Santiago, Chile

Correspondence: Elise S. Droste (

Back to toptop

Perfluorocarbons (PFCs) are potent greenhouse gases with global warming potentials up to several thousand times greater than CO2 on a 100-year time horizon. The lack of any significant sinks for PFCs means that they have long atmospheric lifetimes of the order of thousands of years. Anthropogenic production is thought to be the only source for most PFCs. Here we report an update on the global atmospheric abundances of the following PFCs, most of which have for the first time been analytically separated according to their isomers: c-octafluorobutane (c-C4F8), n-decafluorobutane (n-C4F10), n-dodecafluoropentane (n-C5F12), n-tetradecafluorohexane (n-C6F14), and n-hexadecafluoroheptane (n-C7F16). Additionally, we report the first data set on the atmospheric mixing ratios of perfluoro-2-methylpentane (i-C6F14). The existence and significance of PFC isomers have not been reported before, due to the analytical challenges of separating them. The time series spans a period from 1978 to the present. Several data sets are used to investigate temporal and spatial trends of these PFCs: time series of air samples collected at Cape Grim, Australia, from 1978 to the start of 2018; a time series of air samples collected between July 2015 and April 2017 at Tacolneston, UK; and intensive campaign-based sampling collections from Taiwan. Although the remote “background” Southern Hemispheric Cape Grim time series indicates that recent growth rates of most of these PFCs are lower than in the 1990s, we continue to see significantly increasing mixing ratios that are between 6 % and 27 % higher by the end of 2017 compared to abundances measured in 2010. Air samples from Tacolneston show a positive offset in PFC mixing ratios compared to the Southern Hemisphere baseline. The highest mixing ratios and variability are seen in air samples from Taiwan, which is therefore likely situated much closer to PFC sources, confirming predominantly Northern Hemispheric emissions for most PFCs. Even though these PFCs occur in the atmosphere at levels of parts per trillion molar or less, their total cumulative global emissions translate into 833 million metric tonnes of CO2 equivalent by the end of 2017, 23 % of which has been emitted since 2010. Almost two-thirds of the CO2 equivalent emissions within the last decade are attributable to c-C4F8, which currently also has the highest emission rates that continue to grow. Sources of all PFCs covered in this work remain poorly constrained and reported emissions in global databases do not account for the abundances found in the atmosphere.

1 Introduction
Back to toptop

Perfluorocarbons (PFCs) are fully fluorinated hydrocarbon chemicals that occur ubiquitously in the atmosphere. Even though the magnitudes of their mixing ratios are of the order of parts per trillion molar (ppt) or less, monitoring their atmospheric concentrations is highly warranted by their strong infra-red absorption properties and inertness. As a result of these physical and chemical properties, PFCs can have global warming potentials (GWPs) of several thousands on a 100-year time horizon and extremely long atmospheric lifetimes (Bravo et al.2010; Hartmann et al.2013) (Table 1). Even though their low abundances mean that their contribution to radiative forcing is relatively minor at present, their continued accumulation in the atmosphere might result in a significant contribution to climate change by the end of the century (Zhang et al.2011). Previous studies have shown the emergence of atmospheric PFCs in the 1960s based on firn air samples and a subsequent accelerated increase in the 1990s (Worton et al.2007; Mühle et al.2010; Ivy et al.2012a; Laube et al.2012; Oram et al.2012; Trudinger et al.2016). With the exception of tetrafluoromethane (CF4) (Harnisch et al.1996; Deeds et al.2008; Mühle et al.2010), all known PFCs are thought to be almost exclusively anthropogenic in origin (EDGAR2014). For the lowest molecular weight PFCs, such as CF4 and hexafluoroethane (C2F6), emissions from aluminium production are one of the largest sources. Other sources of PFCs include the use in chemical vapour deposition (CVD) chamber cleaning, plasma etching, refrigerants, fire protection, fluoropolymer production, and – especially for the higher molecular weight PFCs, such as n-C6F14 and n-C7F16 – use as heat transfer fluids and solvents (Robin and Iikubo1992; Mazurin et al.1994; Beu2005; EPA2008; Hartmann et al.2013; EDGAR2014; Mühle et al.2019). Atmospheric trends vary, but mixing ratios are increasing for all known PFCs. Previously reported PFC trends in the atmosphere include CF4, C2F6, octafluoropropane (C3F8), cyclic-octafluorobutane (c-C4F8), decafluorobutane (C4F10), dodecafluoropentane (C5F12), tetradecafluorohexane (C6F14), hexadecafluoroheptane (C7F16), and octadecafluorooctane (C8F18) (Laube et al.2012; Ivy et al.2012a; Oram et al.2012; Mühle et al.2010, 2019). No analytically separated isomers of these PFC species have been reported in the atmosphere before. As the GWPs and atmospheric lifetimes can differ substantially between isomers (Bravo et al.2010), it is important to ascertain the full analytical separation and quantification of such isomers.

Table 1Atmospheric lifetimes and global warming potentials (GWPs) on a 100-year time horizon as reported in the 2013 Intergovernmental Panel on Climate Change (IPCC) report for all PFCs discussed here (Myhre et al.2013).

a The lifetime for i-C6F14 in this work is here assumed to be the same as for n-C6F14. b The GWP for i-C6F14 is estimated based on the radiative efficiency reported in Bravo et al. (2010).

Download Print Version | Download XLSX

The current work focuses on six PFC isomers: c-C4F8, n-C4F10, n-C5F12, n-C6F14, i-C6F14 (perfluoro-2-methylpentane), and n-C7F16. Even though all but one of these PFCs have been reported before, this is the first time that we report on their separated isomers, except for c-C4F8 and n-C5F12, for which we are updating and expanding our previously reported data set (Laube et al.2012; Oram et al.2012). The GWPs on a 100-year time horizon of all six compounds studied in this paper range between 7370 (i-C6F14) and 9540 (n-C4F10) (Table 1) (Bravo et al.2010; Hartmann et al.2013). To assess the recent changes in atmospheric abundances of these PFCs, this work aims to achieve the following objectives.

  1. To extend the 1978–2010 PFC records from Cape Grim, Australia (Oram et al.2012; Laube et al.2012), into early 2018. The air sampled at Cape Grim under baseline conditions has typically travelled over the open ocean many days prior and is known to be very well-mixed (Fraser et al.1999). The Cape Grim time series is therefore an important data set used to study Southern Hemisphere baseline trends of various trace gases (Oram et al.1995; Langenfelds et al.1996; Fraser et al.1999; Reeves et al.2005; Laube et al.2012; Adcock et al.2018). The global trend between 1978 and 2010 for c-C4F8 has been discussed in Oram et al. (2012) and most recently in Mühle et al. (2019), while those for n-C4F10, n-C5F12, n-C6F14, and n-C7F16 have previously been investigated in Laube et al. (2012) and Ivy et al. (2012a). Data acquired since then allow for a re-analysis of the conclusions, given the most recent trends.

  2. To transfer some measurements to a new calibration scale. Re-calibration of the PFCs was necessary, as a number of PFC isomers have recently been identified in our laboratory at the University of East Anglia (UEA), including two isomers of C4F10 (n-CF3(CF2)2CF3 and i-CF(CF3)3), two isomers of C6F14 (n-CF3(CF2)4CF3 and i-CF3(CF2)2CF(CF3)2), and three isomers of C7F16, of which one is the n-isomer (CF3(CF2)5CF3). The chemical structure of the other two isomers is still unknown. A suspicion that the chromatogram peaks previously reported for C6F14 and C7F16 consisted of more than one isomer had already been expressed by Laube et al. (2012). Since then, a new gas-chromatography method reported here (Sect. 2.2) has been developed that led to the separation of these isomers, consequently increasing the accuracy of the measurements.

  3. To study the atmospheric trend and abundance of i-C6F14 for the first time. The separation of isomers allows for more accurate estimations of the impact on the atmospheric radiative forcing, especially since i-C6F14 has a different radiative efficiency than n-C6F14 (Bravo et al.2010).

  4. To compare the abundance and variability of PFCs between the well-mixed, unpolluted air in the Southern Hemisphere, which are well-represented by the Cape Grim data, and the Northern Hemisphere, where most of the emissions occur. Regarding the latter, this study presents recent observations from air samples collected in Tacolneston, United Kingdom, and in Taiwan. Neither of these sites necessarily represent average abundances for the Northern Hemisphere, but the variability in the observations gives information on the geographical distance of these sites to potential sources. This leads to the final objective.

  5. To estimate the global emissions of these six PFCs and compare them to emissions reported in the EDGARv4.2 database (EDGAR2014). PFCs are a part of the 1997 Kyoto Protocol and are covered in the 2014 revised F-Gas regulation of the European Union (EU), which will be a part of the Nationally Determined Contribution of the EU to the Paris Agreement. However, substantial discrepancies have been found between emissions derived from atmospheric background measurements and reported emissions (specifically for the compounds focused on in this study; see Laube et al.2012; Oram et al.2012; Ivy et al.2012b; Mühle et al.2019). The gas-chromatographic separation of isomers is a useful development that may give some insight into, for example, the identity and potential commonality of sources and possibly a better understanding of the changes in sources with time.

The methods are discussed in Sect. 2, including sampling, analysis, emission modelling, uncertainties, and air mass trajectory modelling. The results are split up into three sections. The first section (Sect. 3.1) presents the Cape Grim observations and trends since 1978 for each PFC, with a focus on data since 2010. These are then compared to measurements from Tacolneston and Taiwan. The second section (Sect. 3.2) discusses the emissions based on the emission modelling and compares these to the emissions from the emission report-based EDGARv4.2 database. The third section of the results (Sect. 3.3) covers potential source regions and types based on air mass trajectories simulated by the NAME model for the Taiwan samples. Implications of the findings are outlined in Sect. 4. Section 5 summarises the main conclusions.

2 Methods
Back to toptop

2.1 Sampling and measurement sites

The surface air samples from Cape Grim, Tacolneston, and two sites in Taiwan have been collected over different periods of time and at different frequencies (Table 2).

  1. Air samples collected at Cape Grim Baseline Air Pollution Station, Tasmania, since 1978 have been analysed regularly at UEA since the early 1990s. Between 1978 and 1994 and again since 2012, samples were collected as sub-samples from parent air archive tanks. Between 1995 and 2011, the majority of samples were collected directly in UEA flasks (3 L), and only six samples were sub-sampled from parent air archives. Sampling procedures for this period are outlined in Langenfelds et al. (1996), Oram et al. (1995), and Fraser et al. (1999), and confidence in the stability of the PFCs in the archive is justified in Oram et al. (2012). A selection of these samples has been re-analysed using an updated method in order to revise the n-C4F10, n-C6F14, and n-C7F16 mixing ratios previously published in Laube et al. (2012) (see Sect. 2.2). Samples for UEA have been collected since 1994 in electro-polished stainless steel canisters (Rasmussen) until 2000 and in Silcosteel-treated stainless steel canisters (Restek Corp.) since 2001 (Laube et al.2012; Oram et al.2012).

  2. Samples from Tacolneston, United Kingdom, were collected on a near-bi-weekly basis between July 2015 and April 2017 (Stanley et al.2018; Adcock et al.2018). Samples were collected in Silcosteel-treated stainless steel canisters (Restek Corp.).

  3. In Taiwan, air samples were collected between 2013 and 2016, in March and April each year. The collection site alternated each year between a site in the north (Cape Fuguei) and in the south of Taiwan (Hengchun), starting in the south in 2013 (Laube et al.2016; Oram et al.2017; Adcock et al.2018). Sample collection was aimed at the winter monsoon flow from the north. Even though Hengchun is in the south and air at this time of year travels primarily from the north over Taiwan, the air at Hengchun has a similar trace gas composition to air in Cape Fuguei. This is because Taiwan's topography re-directs the northern winds along the eastern side of the mountains where the population is sparse, and thus it does not pick up additional anthropogenic PFC emissions before it arrives at Hengchun. Due to local PFC contamination issues during sampling, measurements taken in 2013 have been excluded.

Table 2Details on locations of sampling sites and dates, gas-chromatography columns used, and number of samples of which the measurements for each PFC compound are reported here. Differences among the number of samples per compound are due to a combination of the length of the period during which the compound was measured, the chromatography method, and the precision of the measurement affected by baseline distortions.

Download Print Version | Download XLSX

2.2 Analytical methods

All air samples were analysed using an Agilent 6890 Gas Chromatograph coupled to a Waters AutoSpec magnetic sector mass spectrometer (AutoSpec GC-MS). An in-house built cryo-trapping system pre-concentrates 200–300 mL of the air sample at −78C in a HayeSep D-filled loop (80/100 mesh) after passing through a magnesium perchlorate (Mg(ClO4)2) drier to remove any water vapour. Two columns are routinely used as part of AutoSpec analysis at UEA: a GS-GasPro column (length ∼30–50 m; ID: 0.32 mm) and a KCl-passivated CP-PLOT Al2O3 column (length: 50 m; ID: 0.32 mm). As many samples as possible are analysed on both columns to provide a large suite (50+ compounds) of trace gas measurements (e.g. Adcock et al.2018; Laube et al.2016; Leedham Elvidge et al.2018) (Table 2). For the current work, the CP-PLOT column is the most relevant, as it has the ability to separate the higher molecular weight PFCs and their isomers based on their boiling point and symmetry. For compounds analysed on the CP-PLOT column, an Ascarite (NaOH-coated silica) trap is included before the magnesium perchlorate trap to remove CO2, which can otherwise distort or reduce the signal of a number of compounds, including those discussed in this study. The sample is injected into the GC at 100 C. Instrumental drift was accounted for by bracketing samples with calibrated working standards (see Sect. 2.3).

The Taiwan samples collected in 2014 were analysed on the GS-GasPro column only, using a method optimised for a particular set of trace gases that does not include isomers of C6F14 and C7F16. Taiwan samples collected in 2015 and 2016, as well as all Tacolneston samples, are measured on the CP-PLOT column.

The dry air mole fractions are measured and the mixing ratios (in units of parts per trillion, ppt) reported in this study are used as an equivalent to picomole per mole. The absolute calibration method has been described in detail in Laube et al. (2010) and has been described specifically in Oram et al. (2012) for c-C4F8 and in Laube et al. (2012) for C4F10, C5F12, C6F14 and C7F16, including linearity of the detector response and identification. For a detailed description and assessment of the calibration procedure and system, we refer the reader to the Supplement of Laube et al. (2010). The only procedural difference for the calibrations here is the temperature of the system (100 C instead of 80 C). The ions (mz) used for the quantification of c-C4F8, n-C4F10, and n-C5F12 are 131.0, 119.0, and 119.0, respectively. i-C6F14, n-C6F14, and n-C7F16 have been analysed on both quantifying ions (mz) 169.0 and 219.0. Ion mz 169 is used for trend analysis, unless a baseline distortion occurred in the chromatogram, in which case the 219 ion is used. Other studies that used the same set-up and method include Laube et al. (2016) and Adcock et al. (2018).

Even though two isomers of C4F10, two isomers of C6F14, and three isomers of C7F16 have been identified with the current analytical method, the current study only focuses on a selection of these (n-C4F10, i-C6F14, n-C6F14, n-C7F16) for the following reasons. The main signal of the quantifying ion for i-C4F10 (mz 131.0) is not well-separated from the larger n-C4F10 peak on that same ion. Nevertheless, it is still possible to determine the n-isomer for C4F10 on a different quantifying ion, which we confirm to be unaffected by the i-isomer (see Sect. 2.3). Regarding the isomers of C7F16, two out of the three isomers have very small signals resulting in bad precisions. Additionally, these two smallest peaks are not well-separated from each other (but are well-separated from the main isomer peak), and the number of possible isomers with very similar mass spectra is too high to allow for a high-confidence identification. Therefore, their current quantification is too inaccurate and imprecise. However, it is possible to quantify the n-isomer of C7F16 with sufficient precision, because it has a larger signal. For the same reasons, both isomers of C6F14 are quantified here.

Measurement errors (indicated by the error bars on the observational data points) consist of a combination of the 1σ standard deviations of the working standard and sample replicates on the same analysis day. For samples that have been analysed against an older, previously used working standard, the uncertainties include an internal conversion accounted for as the 1σ standard deviation of the peak ratio on inter-comparison days (n=8).

2.3 Calibrations

The mixing ratios in the samples are determined based on a secondary calibration standard, which we refer to as the working standard. Our working standard consists of clean Northern Hemispheric air, which therefore contains all relevant PFCs. For many gases, the mixing ratios in the working standard have been calibrated by NOAA. As certified values do not exist for the relevant PFCs, we calibrate them at UEA using an independent calibration scale (see more details in the Supplement).

The previously reported calibration scales (Laube et al.2012) (now referred to as UEA2010) were revised to accommodate the separated isomers. The new calibration scale is referred to as UEA2018. The calibration procedure has been described in Laube et al. (2010) and has undergone little alteration. High-purity compounds for n-C4F10, n-C6F14, i-C6F14, and n-C7F16 were diluted to ppt levels in a two-step static dilution series using Oxygen Free Nitrogen (OFN) gas (British Oxygen Company). This calibration procedure was repeated three times for each compound and for a range of dilutions (Table S1 in the Supplement). All pure compounds had a 98 % purity or higher at time of acquisition and were subsequently further purified in our lab by subjecting each to three repeated freeze–heating cycles. We used the same compounds as in Laube et al. (2012) for n-C4F10 and n-C6F14, for which we now confirm their isomeric purity. For i-C6F14 and n-C7F16 we used newly purchased pure compounds, which were also confirmed to be isomerically pure. In addition, for n-C7F16 this represents a significant improvement to Laube et al. (2012), as the latter used a technical mixture with an 85 % n-isomer content.

The OFN was filtered to remove small amounts of trace gas contaminations by flowing the gas through 60 cm 14′′ stainless steel tubing filled with HayeSep D, which was submersed in an ethanol-dry-ice cold trap (−78C). The calibration system operates under low pressures (<300 mbar) and high temperatures (100 C). This results in insignificant virial coefficients, and thus the dry-air mole fractions could be calculated using the ideal gas law (Laube et al.2010).

n-C4F10 was re-calibrated to rule out any bias from the observed i-C4F10 isomer (a small side peak which was observed in both the working standard and atmospheric samples) and to confirm the leak tightness of the calibration system (Laube et al.2016). The combined influence of the i-isomer and the leak tightness of the calibration system is <2.8 %. This value has been determined by calculating the relative difference of the n-C4F10 mixing ratio in the working standard using the two (UEA2010 versus UEA2018) calibration scales (see more details in the Supplement). Due to the leak tightness of the system and a lack of observed isomers of c-C4F8 and n-C5F12 in the working standard and atmospheric samples, a revision of the calibration scales was deemed unnecessary for these two compounds. Note that the data for c-C4F8, n-C4F10, and n-C5F12 in Oram et al. (2012) and Laube et al. (2012) are included in the current work.

The confidence in the complete separation of all isomers for C6F14 and C7F16 is based on the ratio of the peak areas of the two most abundant ions into which these compounds are ionised during mass spectrometry: mass-to-charge (mz) ratios 219.0 and 169.0. A comparison of the 219 : 169 ion ratio is made between the Cape Grim air samples and the calibrations done (Fig. S1). Since the calibration samples are based on dilutions of high-purity isomer compounds (>98 %), a significant deviation from the ion ratios in the air samples compared to the ion ratios in the calibration samples might suggest that the peak measured in the chromatogram for the air samples actually consists of more than one isomer. In turn, if the ion ratios measured in the air samples are similar to the ion ratios measured in the calibration samples, then confidence can be attributed within the uncertainties that the signal measured in the air samples is for one particular isomer only. For i-C6F14, n-C6F14, and n-C7F16, the ion ratio in the Cape Grim samples did not differ significantly from the ion ratio in the calibrations within 1σ measurement uncertainty (Fig. S1). It should be noted that even if other isomers do exist and co-elute, they (a) likely do not have a trend that is substantially different from any of the C6F14 and C7F16 isomers discussed here and (b) currently have mixing ratios that are extremely small compared to those of the main isomers.

Table 3Mixing ratios determined in the working standard, average analytical precisions of the calibrations, and measures of accuracies of CFC-11 calibrations expressed as the average relative differences compared to NOAA scales ((UEA  NOAA)  UEA). The analytical precision is the relative 1σ standard deviation of the compound signal in multiple working standard runs analysed on the day of each calibration. The ranges of analytical precisions and differences to NOAA are shown here for each compound in brackets and include all three calibrations done in this work for that compound. For details on the calibration of c-C4F8 and n-C5F12, consult Oram et al. (2012) and Laube et al. (2012).

a New calibration (UEA2018). b Old calibration (UEA2010). c Oram et al. (2012). d Laube et al. (2012).

Download Print Version | Download XLSX

Calibrated mixing ratios for the working standard (i.e. unpolluted Northern Hemispheric air from 2017) are given in Table 3, along with the average analytical precision per compound. The latter is the relative 1σ standard deviation of the compound signal in multiple working standard runs analysed on the day of each calibration. To obtain a measure of accuracy of the calibrations (as in Laube et al.2010, 2012), trichlorofluoromethane (CFC-11) was diluted along with the pure PFCs as a reference compound, as our working tank has been calibrated by the globally recognised GMD (Global Monitoring Division) of the NOAA-ESRL (National Oceanic and Atmospheric Administration – Earth System Research Laboratory) for CFC-11. For all UEA calibrations, the average offset with values determined by NOAA shows that the UEA-calibrated concentrations for CFC-11 are consistently slightly lower by on average 4.2 %±0.3 %. To put the accuracy of these calibrations into perspective, the accuracies of calibrations for the PFCs in common with Ivy et al. (2012a) and Laube et al. (2012) are between 4.0 % and 7.8 % and between −5.5 % and 2.8 %, respectively.

The improvement of the new calibration scale by separating the isomers of C6F14 and C7F16 is substantial. Even though the n-isomer of both C6F14 and C7F16 is dominant for both PFCs, the old calibration scale would have overestimated the mixing ratios of these n-isomers by 20 % and 11 %, respectively. These levels of improvement in the UEA2018 calibration scale are determined based on the relative difference in concentration determined for the respective n-isomers in the working standard between using the UEA2010 and UEA2018 calibration scales. The new UEA calibration scale thus allows for a more accurate analysis of atmospheric mixing ratios of these PFCs.

2.4 Emission modelling

The annual global emissions are derived using a two-dimensional global atmospheric chemistry-transport model (run using the Facsimile, version 7, software; also see Laube et al.2012; Oram et al.2012; Reeves et al.2005; Fraser et al.1999). The model domain consists of 24 equal-area, zonally averaged latitudinal bands, which each have 12 atmospheric layers between the surface and the maximum altitude represented in the model. Each layer is attributed with a height of 2 km, resulting in a total altitude of 24 km.

Due to the high stability of the C–F bonds, PFCs are considered to be chemically inert, and thus the model is set to have no chemical or photolytic loss for these compounds. Photolytic loss would only be relevant at mesospheric altitudes (Morris et al.1995), which are not represented in the model domain. The only sink for PFCs in the model is the diffusive loss at the top boundary of the model domain, where the diffusive loss is controlled by a fixed ratio of the mixing ratio of the PFC in the top atmospheric layer to the layer above the model domain, thereby creating a gradient. If there is no gradient, i.e. the ratio is equal to 1, then there is no diffusive loss. The upper boundary ratio used for each modelled gas is set such that the diffusive loss replicates the lifetime of the gas in the upper atmosphere (i.e. above the model domain). For the PFCs, this is effectively their atmospheric lifetime. As PFCs have lifetimes that are much longer than the time period studied, the gradient out of the model is set very low (i.e. ratio of 0.997; Laube et al.2012).

The emission distribution is based on the global distribution of reported C6F14 emissions in the EDGARv4.2 data set in 2005. This set-up is kept the same as in Laube et al. (2012), because the EDGARv4.2 emission distribution for C6F14 between 2005 and 2010 has not changed significantly and also does not significantly affect the simulated mixing ratios. The reason the EDGARv4.2 data for C6F14 are used for all PFCs is that the discrepancy between atmospherically derived and reported emissions is smallest for C6F14 (Laube et al.2012). In the model, PFC emissions are set to occur 99 % in the northern mid-latitudes. This is a realistic assumption according to reported emissions recorded in the EDGARv4.2 database, which indicate that 98 % to 100 % of the PFC emissions occur in the Northern Hemisphere, depending on the PFC compound (EDGAR2014).

Model runs start in 1934 and end in 2018. The annual emissions are iteratively altered to obtain a best fit of the modelled mixing ratios to the observed mixing ratios at Cape Grim (located within the latitudinal band spanning 35.7 to 41.8 S). The earliest annual PFC emissions set in the model start in the 1950s, when certain PFCs were first detected above detection limits in firn air (Laube et al.2012). The best fit is determined by minimising the sum of least squares between the modelled and observed Cape Grim mixing ratios (Table S3).

2.5 Trend and emission uncertainties

Two versions of uncertainties have been computed for the model-derived PFC trends for Cape Grim and Tacolneston data (visualised as the shaded envelopes in the figures). The first version will be referred to as the “trend uncertainty”, which consists of the annually averaged analytical uncertainty, the averaged model-fit uncertainty of the Cape Grim trend, and a 5 % modelling uncertainty (Table S2). The analytical uncertainty is described in Sect. 2.2. The model-fit uncertainty is the average relative difference between the observed mixing ratio and the simulated mixing ratio. The modelling uncertainty accounts for uncertainties regarding the model transport scheme and is the error in simulating the concentration at Cape Grim of a long-lived gas that is emitted primarily in the Northern Hemisphere at a well-known rate (Reeves et al.2005). Lifetime uncertainties are not included, because they are too long for the time period studied here and would not have a significant effect on the results. Note that the trend uncertainty does not include the calibration uncertainty. Even though it is important to consider the calibration uncertainty in order to evaluate the complete uncertainty of the model simulation, it does not affect the trend. However, to give an overview of all uncertainties, the calibration is included in the “total uncertainty”, which is the second type of uncertainty referred to in this work (Table S2). An additional uncertainty for the trends of c-C4F8 and n-C5F12 is added in the total uncertainty, which is the error in the conversion of the mixing ratio between internal working standards, as these two PFCs have not been re-calibrated for the current study (see Sect. 2.3).

The “trend uncertainty” and the “total uncertainty” for each PFC were used to determine their respective uncertainty bands around the observed mixing ratio trend. The model was then re-run to find the emissions that fitted the simulated maximum and minimum mixing ratios best. The result is the minimum and maximum emissions according to both types (“trend” and “total”) of uncertainties.

2.6 NAME modelling

This work used the UK Met Office's Langrangian particle dispersion model, Numerical Atmospheric Modelling Environment (NAME) (Jones et al.2007), for tracking and understanding the origin of air masses arriving at Taiwan. The model was run in the backward mode for 12 d long simulations to generate the footprints of where the air sampled during the campaigns had previously been close to the Earth's surface (see Fig. 1). The analysis begins by releasing batches of 30 000 inert backward particles over a 3 h period encompassing the collection time of each sample. Over the course of the 12 d travel time, the locations of all particles within the lowest 100 m of the model atmosphere were recorded every 15 min on a grid with a resolution of 0.25 longitude and 0.25 latitude. The trajectories were calculated using three-dimensional meteorological fields produced by the UK Met Office's Numerical Weather Prediction tool, the Unified Model (UM). These fields have a horizontal grid resolution of 0.23 longitude by 0.16 latitude and 59 vertical levels below ∼30 km.

Figure 1Example of the particle footprints generated by the NAME model indicating where air sampled in Taiwan came from. The colour scale is logarithmic and represents the calculated time-integrated particle density (g s m−3) within the surface layer (0–100 m) during 12 d prior to the sampling days given a point release at Taiwan of g s−1. Darker colours indicate a greater influence of the region on the chemical composition of the air sampled in Taiwan compared to the lighter colours.

Figure 2Regions used in the NAME model to quantify the simulated contributions to particle densities in Taiwan.

In order to quantify the contribution of various regions to each footprint, the domain was divided into 15 regions using shape files produced by ArcGIS, a geographic information system (GIS) for working with maps and geographic information (Fig. 2). The segregation of China from the East Asia region category has enabled a more detailed analysis to be conducted to determine which specific Chinese sub-regions could contribute to the variations in the PFC mixing ratios in Taiwan. The contribution of each region is quantified by summing the particle concentration (g s m−3) in each grid cell within each shape file (see Fleming et al.2012, and O'Shea et al.2017).

For the purpose of identifying the possible sources of PFCs, the NAME footprints were used to calculate the mixing ratio of a tracer at the Taiwan measurement site, given emissions from a range of different emission sectors. Carbon monoxide (CO) was chosen as the tracer because it is emitted from a number of different emission sources and because widely used and tested bottom-up emission estimates of CO are available from the Representative Concentration Pathway 8.5 (RCP8.5) inventory (Riahi et al.2011; Van Vuuren et al.2011) (, last access: 26 October 2018). Moreover, the RCP8.5 CO emissions are divided into various emission sectors (industry, power plants, solvents, agricultural waste burning, waste, forest burning, grassland burning, residential, international shipping, surface transportation, and agriculture), allowing us to test whether there was a correlation between the observed PFC mixing ratios and whether the air modelled to have arrived at the measurement site had likely been subject to emissions from any particular sector.

The NAME footprints were combined individually with the distribution of CO emissions from each sector to calculate a modelled mixing ratio of the emitted species at the measurement site (see Supplement). Emissions were taken for the year 2010, for the timescale of each footprint (12 d) (more details can be found in Oram et al.2017). A correlation analysis was performed between the modelled CO from each sector and the measured PFC mixing ratios. Industry (combustion and processing) and solvent applications are expected to show some correlations with PFC mixing ratios, as they are most closely associated with PFC sources, and so we would expect their emissions to be co-located. Results of these analyses are found in Sect. 3.3.

3 Results and discussion
Back to toptop

3.1 Atmospheric abundances

3.1.1 c-C4F8 trends

Out of all PFC compounds reported in this work, c-C4F8 is the most abundant. c-C4F8 seems to still have substantial Northern Hemispheric emissions, for which the measured abundances provide three potential pieces of evidence: (1) a considerable, and even accelerating, rate of increase in concentration over time, (2) the existence of a relatively large interhemispheric gradient given its long atmospheric lifetime, and (3) large variations above background levels at Northern Hemispheric sites. Each of these pieces of evidence will be illustrated with the Cape Grim, Tacolneston, and Taiwan observations, respectively.

Figure 3(a) Mixing ratios at Cape Grim of c-C4F8 between 1978 and 2018 (blue circles). Data prior to 2010 are shown as empty markers (Oram et al.2012), while data from samples collected after 2010 are shown as filled markers to illustrate the part of the time series that is extended in the current work. The atmospheric trend simulated by the model is represented by the red line. Total uncertainties (light blue) and trend uncertainties (dark blue) are indicated with the shaded areas along the trend line. Results for the c-C4F8 trend published in Mühle et al. (2019) are shown in the dotted black line. (b) Mixing ratios after 2010 for Cape Grim (circles), for Tacolneston (diamonds), and Taiwan (squares). The red line indicates the modelled Southern Hemisphere baseline trend at Cape Grim, while the dashed red line indicates the modelled trend for mixing ratios at Tacolneston. Total and trend uncertainties are indicated with the shaded grey areas along the Tacolneston dashed-trend line.


The c-C4F8 mixing ratio in the Southern Hemisphere increased from 0.31 ppt in 1978 to 1.52 ppt by the end of 2017 (Fig. 3a). As reported previously by Oram et al. (2012), the mixing ratios seemed to stabilise in the 1990s, but picked up again in the 2000s. The extended time series since 2008 indicate a continuing increasing trend similar to the one pre-1990, totalling a 27 % enhancement in the last 8 years (2010–2018). In fact, atmospheric c-C4F8 abundances show accelerating growth in recent years, potential reasons for which are explored in Sect. 3.2.1. Mühle et al. (2019) reported a slightly higher mixing ratio of 1.61 ppt at the end of 2017 in the high-latitude Southern Hemisphere. Since Mühle et al. (2019) already compared their Southern Hemisphere time series to ours until 2008 (published in Oram et al.2012), we here focus our comparison on the extended part of the record, i.e. from 2008 until 2017. The average difference of annual mixing ratios for this part of the time series is 0.09±0.01 ppt (1σ). This difference, which is equivalent to an average of 7 %, is to some degree likely due to the two different independent absolute calibration scales. However, since the difference is clearly not constant over time (see Fig. 3), which has also been discussed in Mühle et al. (2019), no simple conversion factor can be derived for the two data sets.

The mixing ratios at Tacolneston, UK, seem to have a consistent offset compared to the observations at Cape Grim on average of 0.07 ppt (based on monthly averages; Fig. 3b), equivalent to an interhemispheric ratio of 1.05 (average enhancement of Tacolneston mixing ratio measurements against the background trend at Cape Grim). Note that mixing ratios at Tacolneston have not been shown to be representative of Northern Hemisphere background mixing ratios. However, the lack of variability in the data suggests that Tacolneston is not in close proximity to any major sources of c-C4F8. Such a proximity of sources should cause multiple occurrences of substantially higher mixing ratios, as can be seen in the Taiwan data (discussed below). When considering the simulated atmospheric trend from the emission model for Tacolneston, it can be noted that it overestimates the observed concentrations, despite the well-fitted Southern Hemispheric trend. The difficulty in validating trace gas trends for Tacolneston is 3-fold. (1) The time series for Tacolneston only covers the time period between 2015 and 2017, which does not enable a validation of the atmospheric trends over the period during which c-C4F8 has been emitted. (2) The model is two-dimensional and as such represents zonal averages. In reality, emissions will occur at specific longitudes, which leads to heterogeneity in atmospheric concentrations at the latitudes of the emissions, even for these very long-lived compounds. This is clearly illustrated by the Taiwan data (see below). (3) The mixing ratios simulated by the model are dependent on a given global distribution of emissions, which harbours a substantial amount of uncertainty. Varying the distribution of the 99 % of global emissions within the Northern Hemisphere does not significantly affect the simulated mixing ratios for Cape Grim, but naturally it does for Northern Hemispheric sites, such as Tacolneston. It is likely that the distribution for c-C4F8 in the Northern Hemisphere requires revision to include, for example, a distribution that changes over time rather than a fixed distribution (e.g. see a similar discussion for CFC-113a in Adcock et al.2018). It is not unrealistic to consider a Northern Hemispheric distribution that moves from higher to lower latitudes over time, as the industries that utilise PFCs and other halocarbons, such as the electronics industry, have moved from more northern latitudes to regions such as South-East Asia (Montzka et al.2009, for hydrochlorofluorocarbons). However, available information on regions that are sources of c-C4F8 (and other PFCs in general) provide insufficient evidence to base a well-founded emission distribution on (see Sect. 3.2). The global emission distribution for c-C4F8 is thus kept consistent with that for the other PFCs discussed in this work, for which the simulated Tacolneston trend compares reasonably well with the observations, i.e. within the measurement and modelling uncertainties.

The c-C4F8 mixing ratios observed in samples collected in Taiwan are contrasting to those at Cape Grim and Tacolneston because of their large day-to-day variability (Fig. 3b). The maximum measured mixing ratios for each year in Taiwan are on average 17 % higher than the average simulated mixing ratio in the respective year. This strongly suggests a close proximity to one or more PFC sources. Considering that much of the air sampled in Taiwan during the campaigns has been transported from the direction of China (Laube et al.2016; Adcock et al.2018) and that Taiwan itself has a major semi-conductor industry as well (Saito et al.2010), it is not surprising to find such high levels of c-C4F8 at these sites.

3.1.2 n-C4F10 and n-C5F12 trends

Between 1978 and 2018, the background mixing ratios for n-C4F10 and n-C5F12 at Cape Grim have increased from 0.02 to 0.18 ppt and from 0.02 to 0.15 ppt, respectively (Fig. 4a). In contrast to c-C4F8, n-C4F10 and n-C5F12 seem to have much smaller Northern Hemispheric sources at the present time. First of all, the Southern Hemispheric mixing ratio growth rate has declined since ∼2000. n-C4F10 and n-C5F12 mixing ratios at Cape Grim only increased by 9 % and 6 % between 2010 and 2018, respectively. Even though growth rates are currently slower than they were in the 1990s, the continuing increasing mixing ratios shown in the extended time series indicate that sources still exist. The similarity in the trends of n-C4F10 and n-C5F12 suggests that these sources might be similar or common. The ratio between n-C4F10 and n-C5F12 remained relatively stable between 1978 and 2018 at an average of about 1.3±0.07.

Figure 4(a) Mixing ratios at Cape Grim of n-C4F10 (blue circles) and n-C5F12 (green circles) between 1978 and 2018. Data prior to 2010 are shown as empty symbols (Laube et al.2012), while data from samples collected after 2010 are shown as filled symbols to illustrate the part of the time series that is extended in the current work. Note that the data from Laube et al. (2012) have been converted to a new and improved calibration scale. The atmospheric trends simulated by the model are represented by the red line for n-C4F10 and by the magenta line for n-C5F12. Total uncertainties (light blue for n-C4F10, light green for n-C5F12) and trend uncertainties (dark blue for n-C4F10, dark green for n-C5F12) are indicated with the shaded areas along the trend line. (b) n-C4F10 mixing ratios after 2010 for Cape Grim (blue circles), for Tacolneston (orange diamonds), and Taiwan (light blue squares); n-C5F12 mixing ratios after 2010 for Cape Grim (green circles), for Tacolneston (light green diamonds), and Taiwan (magenta squares). The red and magenta lines indicate the modelled Southern Hemisphere baseline trend at Cape Grim for n-C4F10 and n-C5F12, respectively, while the dashed red and magenta lines indicate the modelled trend for mixing ratios at Tacolneston for n-C4F10 and n-C5F12, respectively. Total and trend uncertainties are indicated with the shaded grey areas along the Tacolneston dashed-trend lines for both n-C4F10 and n-C5F12.


Second of all, as would be expected from a slowly increasing trace gas – and based on the samples collected at Tacolneston and at Cape Grim (Fig. 4b) – there does not seem to be an interhemispheric gradient for these PFCs that is discernible given the measurement uncertainties. This is also captured in the simulated mixing ratios for both sites by the model. If the interhemispheric ratios of the mixing ratios of n-C4F10 and n-C5F12 were to be similar to that of c-C4F8, then the differences between the Cape Grim and Tacolneston data for these two PFCs would have to be at least 0.01 ppt. This is not the case, but even so, these differences in mixing ratios are too small to discern given the measurement uncertainty.

With the current measurement precisions and low growth rates, the data allow no detailed conclusions on the appropriateness of the global emission distributions used for the model simulations of these two gases. As the Cape Grim observations are representative of the well-mixed Southern Hemispheric air, the similarity in abundances between the Tacolneston and Cape Grim observations might mean that the Tacolneston site is situated relatively far away from n-C4F10 and n-C5F12 sources. Additionally, it indicates that the emissions of n-C4F10 and n-C5F12 occur at a low rate, which means that these PFCs will be relatively well-mixed in the Northern Hemisphere. Note that the n-C4F10 data at Tacolneston do display some scatter, but it is well within the uncertainties of the Cape Grim record.

Third of all, smaller Northern Hemispheric sources of n-C4F10 and n-C5F12 compared to c-C4F8 are reflected by the smaller elevations of their mixing ratios measured in Taiwan relative to those of c-C4F8, despite their enhancement compared to Cape Grim measurements (Fig. 4b). Maximum values are around 0.25 ppt for n-C4F10 and 0.20 ppt for n-C5F12. On average, the mixing ratios measured in Taiwan are 5 % higher than the average mixing ratios simulated in the respective year for Cape Grim for both n-C4F10 and n-C5F12. Even though this provides some evidence that there are sources of n-C4F10 and n-C5F12 in East and South-East Asia, these sources do not appear to have a large impact on the interhemispheric ratio between the unpolluted Northern Hemisphere and Southern Hemisphere sites. Much lower mixing ratio elevations for both n-C4F10 and n-C5F12 were observed in Taiwan in 2015 compared to 2014 and 2016. This is likely related to different air mass origins of samples collected in this particular year (see also Oram et al.2017; Adcock et al.2018).

3.1.3 i-C6F14 and n-C6F14 trends

The trends of i-C6F14 and n-C6F14 are somewhat similar to each other. n-C6F14 mixing ratios at Cape Grim have increased from 0.01 ppt in 1978 to 0.22 ppt in 2018 (Fig. 5a). The fastest increase in atmospheric abundance occurred in the 1990s and has since slowed down. Mixing ratios appear to be approaching stabilisation, but have still increased by about 9 % since 2010. Observed long-term trends are similar for i-C6F14, although its atmospheric abundance is much lower (Fig. 5a). This points towards commonality or co-location in sources with n-C6F14. i-C6F14 mixing ratios were <0.01 ppt before 1987, but reached 0.07 ppt by 2018 (note that the earliest observation for i-C6F14 shown within the Cape Grim archive is in 1987, due to limited precision of measurements on samples collected prior to that year). In the current work, it is the PFC with the lowest atmospheric abundance throughout the Cape Grim record. Despite the similarities between these two isomers, a notable difference is observed in their recent rate of change: the relative growth rate of i-C6F14 since 2010 (∼19 %) is double that of n-C6F14. Given their similar physico-chemical properties it is likely that these two isomers are emitted to the atmosphere by the same anthropogenic processes. If their production has not changed much over time, the i-C6F14 : n-C6F14 ratio is expected to remain relatively stable over that same time period. However, the i-C6F14 : n-C6F14 ratio has increased from an average of 0.26±0.01 (1σ standard deviation) between 2003 and 2008 to an average of 0.3±0.01 between 2013 and 2018 (1σ) (Fig. S2). This might indicate a shift in production method or new sources that emit these isomers in different ratios.

Figure 5(a) Mixing ratios at Cape Grim of n-C6F14 (blue circles, left axis) between 1978 and 2018 and i-C6F14 (green circles, right axis) between 1987 and 2018. Data prior to 2010 are shown as empty symbols (Laube et al.2012), while data from samples collected after 2010 are shown as filled symbols to illustrate the part of the time series that is extended in the current work. Note that the data from Laube et al. (2012) have been converted to a new and improved calibration scale. The atmospheric trends simulated by the model are represented by the red line for n-C6F14 and by magenta line for i-C6F14. Total uncertainties (light blue for n-C6F14, light green for i-C6F14) and trend uncertainties (dark blue for n-C6F14, dark green for i-C6F14) are indicated with the shaded areas along the trend line. Note that there is only a very small difference between the trend and total uncertainties, and thus the latter is difficult to distinguish. (b) n-C6F14 mixing ratios after 2010 (left axis) for Cape Grim (blue circles, left axis), for Tacolneston (orange diamonds), and Taiwan (light blue squares); i-C6F14 mixing ratios (right axis) after 2010 for Cape Grim (green circles, right axis), for Tacolneston (light green diamonds), and Taiwan (magenta squares). The red and magenta lines indicate the modelled Southern Hemisphere baseline trend at Cape Grim for n-C6F14 and i-C6F14, respectively, while the dashed red and magenta lines indicate the modelled trend for mixing ratios at Tacolneston for n-C6F14 and i-C6F14, respectively. Total and trend uncertainties are indicated with the shaded grey areas along the Tacolneston dashed-trend lines for both n-C6F14 and i-C6F14.


Observed Tacolneston mixing ratios for n-C6F14 and i-C6F14 in 2015 and 2016 compared to those at Cape Grim in the same years show no clear interhemispheric gradient (Fig. 5b). The average monthly differences between the observations at these two sites are 0.008 and 0.005 ppt for n-C6F14 and i-C6F14, respectively, which approaches or even exceeds the differences in mixing ratio one would expect if the interhemispheric ratio was at least that of c-C4F8 (0.01 ppt for n-C6F14 and 0.003 for i-C6F14). However, such differences in mixing ratios are too small compared to the measurement uncertainties, and thus no interhemispheric gradient is discernible based on these data.

The mixing ratios at Tacolneston for the isomers of C6F14 have greater variability in observed mixing ratios compared to those of the other PFCs discussed in this work. Mixing ratios occasionally exceed the modelled uncertainty envelopes for Tacolneston, such as in 2015, but none of these excursions exceeded the 2σ measurement uncertainties. The simulated trend for the Tacolneston data corresponds well to the lower end of the variability in the observations.

No measurements of either n-C6F14 or i-C6F14 were made for Taiwan samples taken in 2014, as these samples were analysed on a GC column that did not allow for a separation of these isomers (see Sect. 2.2). Mixing ratios in Taiwan in 2015 and 2016 are again extremely variable on short timescales and reach higher abundances than at Cape Grim or Tacolneston, ranging between 0.21 and 0.47 ppt for n-C6F14 and between 0.06 and 0.13 ppt for i-C6F14 (Fig. 5b). On average, the mixing ratios measured for both n-C6F14 and i-C6F14 are 9 % higher in Taiwan than at Cape Grim.

Even though substantial variability of mixing ratios above the Cape Grim baseline for both C6F14 isomers is observed at Tacolneston and in Taiwan, the low absolute mixing ratios in both hemispheres combined with the measurement uncertainties do not allow for any further conclusions on the hemispheric distribution of their emissions.

3.1.4 n-C7F16 trends

Similarly to c-C4F8, n-C7F16 shows substantial sources in the Northern Hemisphere, as is suggested by its rate of change, interhemispheric gradient, and large variations above background levels in Taiwan (Fig. 6a, b). Mixing ratios of n-C7F16 at Cape Grim increased from 0.01 ppt in 1978 to 0.11 ppt in 2018. Its trend is different from those of the other PFC compounds reported here, as its mixing ratios have been continuously increasing at an approximately constant rate since 1985. Between 2010 and 2018, atmospheric background levels in the Southern Hemisphere have increased by 21 % and show no signs of slowing down. Two statistical outliers appear in measurements on samples collected in 1988 and 1992.

Figure 6(a) Mixing ratios at Cape Grim of n-C7F16 between 1978 and 2018 (blue circles). Data prior to 2010 are shown as empty markers, while data from samples collected after 2010 are shown as filled markers to illustrate the part of the time series that is extended in the current work. Note that the data from Laube et al. (2012) have been converted to a new and improved calibration scale. The atmospheric trend simulated by the model is represented by the red line. Total uncertainties (light blue) and trend uncertainties (dark blue) are indicated with the shaded areas along the trend line. (b) Mixing ratios after 2010 for Cape Grim (circles), for Tacolneston (diamonds), and Taiwan (squares). The red line indicates the modelled Southern Hemisphere baseline trend at Cape Grim, while the dashed red line indicates the modelled trend for mixing ratios at Tacolneston. Total and trend uncertainties are indicated with the shaded grey areas along the Tacolneston dashed-trend line.


Observations at Tacolneston and Cape Grim together display an interhemispheric gradient for n-C7F16 of 1.04, which is captured by the model. However, even though this interhemispheric gradient is comparable to that for c-C4F8, the measurement uncertainties are too large to be able to conclude based on these data that the interhemispheric gradient is clearly discernible. The Taiwan measurements range between 0.10 and 0.22 ppt, which is consistent with the large variability seen at Taiwan's measurement sites for all other PFCs analysed in this work. Mixing ratios in Taiwan are on average 15 % higher than at Cape Grim. Having discussed the atmospheric trends, the next step is to investigate how global emissions changed over time and what their current status is.

3.2 Global emissions

3.2.1 c-C4F8 emissions

The model-derived global c-C4F8 emissions have changed substantially over time (Fig. 7). An increase in emission rates occurred in the 1980s, peaking around 1.7 Gg yr−1 (gigagrams per year) in 1986. A possible contributing source to these emissions might have been the increasing use of liquid PFC coolants in temperature control units within semiconductor processing. The design of these temperature control units, which was initially made to be used with water and glycol, was inadequate to prevent any leakage of the replacement PFCs from the pumps and seals (EPA2008). The following rapid decrease in emissions between 1985 and 1995 despite the growing demand in the semiconductor industry remains unexplained (Oram et al.2012), especially because it precedes the Kyoto Protocol and any formal initiative taken by the semiconductor industry to reduce PFC emissions.

Figure 7Estimated global emission rates for c-C4F8 in this study (full black line), in Oram et al. (2012) (dotted magenta line), in Mühle et al. (2019) (dotted green line), and as reported in the EDGARv4.2 database (dashed black line). Note that the EDGARv4.2 emission rates are plotted on the secondary axis. Uncertainty envelopes include contributions from measurements, modelling, and calibrations.


Emissions have subsequently continued to increase from the 1990s minimum. The global emission trend is consistent up to 2006 with our previously reported data by Oram et al. (2012) and agrees well on average with the most recent results of Mühle et al. (2019), although our data show larger variability for some periods. However, the extended time series reveals that emission rates do not show any signs of stabilisation and have instead continued to increase, and even accelerate after ∼2012. Annual emissions are now approaching rates of 2.0 Gg yr−1, which are, within uncertainties, comparable to rates determined for the mid-1980s. This compares quite well with the estimated global emissions of 2.2 Gg yr−1 in 2017 by the recent work of Mühle et al. (2019). The average difference between Mühle et al. (2019)'s and our emission estimates for 2008–2017 is 0.17±0.09 Gg yr−1 (1σ). However, the two data sets agree within the uncertainty ranges. Emission rates have increased by approximately 50 % since 2010. The increasing emissions are interesting especially in the light of the substantial PFC emission-reduction efforts by the semiconductor industry since 1999 by means of process optimisation, alternative chemistries, PFC recovery, and abatement (Beu2005). The report by the International Sematech Manufacturing Initiative (ISMI) refers to effective and successful reduction efforts for all PFCs by the semiconductor industry, explicitly describing those for CF4, C2F6, C3F8, and c-C4F8 (Beu2005). Giving credibility to these reduction efforts by the World Semiconductor Council (WSC), it is perhaps unexpected that c-C4F8 emissions are increasing at their current rate. However, the list of WSC members is not exclusive to all major semiconductor manufacturers in the world, and while technology is becoming more efficient, demand within the electronics industry may be offsetting PFC reduction efforts. Additionally, c-C4F8 is one of the candidates as a replacement for lower molecular weight PFCs, such as CF4 and C2F6, which have longer atmospheric lifetimes and higher GWPs (Tsai et al.2002). Finally, the results of Mühle et al. (2019) point toward fluoropolymer production as another major source of atmospheric c-C4F8.

An enormous discrepancy is evident between c-C4F8 emissions based on reported production (bottom-up) in the EDGARv4.2 database and emissions derived based on the model simulations of the atmospheric mixing ratios (top-down). This finding continues a trend noted for our earlier measurements (Oram et al.2012). The c-C4F8 emission rates from the EDGARv4.2 database are lower than model-derived values by 2 orders of magnitude. This reveals the continued lack of reporting by nations on PFC production and emissions. However, a fraction of the discrepancies between atmospherically derived and report-derived emissions is possibly attributable to inadvertent PFC production and emission, despite efforts to reduce PFC release into the atmosphere through leakage (Beu2005). Some of these inadvertent PFC sources may remain unidentified, complicating attempts to locate source regions and source types. Even the recent evidence towards fluoropolymer production sources by Mühle et al. (2019), which is consistent with their much larger observational data set for c-C4F8, still remains somewhat speculative on a global scale.

3.2.2 n-C4F10 and n-C5F12 emissions

Emission rates for n-C4F10 and n-C5F12 are comparable to each other in terms of both trend and magnitude, which is consistent with previous work (Laube et al.2012; Ivy et al.2012b) (Fig. 8). Emissions of both compounds initially rise steadily to peak at around 0.30 Gg yr−1 in the mid-1990s, after which they decrease. n-C4F10 and n-C5F12 are used in refrigeration technology and fire extinguishing methods (Robin and Iikubo1992; Mazurin et al.1994; EDGAR2014). The upward emission trend prior to the mid-1990s may be attributed to the growth of the electronics industry, which required cooling technology, and the substitution of ozone-depleting chlorofluorocarbons (CFCs) in refrigeration and fire extinguishing applications. Subsequent replacement strategies of PFC applications after the signing of the Kyoto Protocol in 1997, such as substitution by alternative chemistries, may have triggered the decline in these emissions.

Figure 8Estimated global emission rates for n-C4F10 and n-C5F12 (full black and red lines, respectively) in this study, global emission rates for C4F10 as reported in the EDGARv4.2 database (dashed black line; no distinction between isomers), modelled emissions by Ivy et al. (2012b) for C4F10 and C5F12 (dash-dotted blue and brown lines, respectively; no distinction between isomers), and modelled emissions by Laube et al. (2012) for C4F10 and C5F12 (dotted black and red lines, respectively; no distinction between isomers). Note that the maximum EDGARv4.2 emission rate for n-C5F12 is less than 0.6×10-4 and is therefore not plotted here. Shadings illustrate the uncertainties.


Maximum emission rates reported in Laube et al. (2012) are 0.27 Gg yr−1 for n-C4F10 and 0.31 Gg yr−1 n-C5F12, which agrees within the uncertainties of the current work. A similar conclusion is valid for the maximum emissions rates in Ivy et al. (2012b) (Fig. 8). In order to fit the observed mixing ratios, emission rates are required to have stabilised in the last decade. In the current work, n-C4F10 stabilises at lower emission rates than in the late 1970s. These trends are also qualitatively comparable to the work by Ivy et al. (2012b). Laube et al. (2012) already suspected a stabilisation in emission rates, which is supported by the extended data set presented in the current study. This low emission rate is consistent with (1) the low NH–SH gradient seen in the Cape Grim and Tacolneston observations and (2) the Taiwan data only showing moderately elevated mixing ratios.

As with c-C4F8, it is apparent that the top-down and bottom-up approaches estimate very different annual global emissions for n-C4F10 (Fig. 8) and n-C5F12 (emissions <5.3×10-5 Gg yr−1). It is important to note that data reported in the EDGARv4.2 database do not distinguish between different isomers. Hence, the reported emissions in the EDGARv4.2 database are likely a combination of various isomers. Nevertheless, it is still worthwhile comparing these data on reported emissions to the modelled emissions for the n-isomers in the current study, because the abundance of other isomers besides the n-isomer is very low for C4F10 and not detectable for C5F12. The discrepancy between the modelled emissions for the isomers in the current study and the total reported emissions in the EDGARv4.2 database is around 2 orders of magnitude for C4F10 and 4 orders of magnitude for C5F12, with maximum EDGARv4.2 emission rates of 0.02 Gg yr−1 (C4F10) and of 5.26×10-5 Gg yr−1 (C5F12) (EDGARv4.2 data for n-C5F12 are not shown) (Fig. 8).

3.2.3 i-C6F14 and n-C6F14 emissions

One of the most striking features of the top-down derived emission trend for n-C6F14 is the sudden increase in emissions in the mid-1990s (Fig. 9). Emission rates are estimated to be constant at less than 0.2 Gg yr−1 up until ∼1994, but then increase by a factor of 6 to reach 1.21±0.10 Gg yr−1 by ∼1997. n-C6F14 is liquid at room temperature and thus has widely used applications as a heat transfer fluid. As described in Sect. 3.2.1, the switch to using PFCs as heat transfer fluids was rapid due to their effectiveness at regulating heat, which was necessary to cope with the increasing demand of the semiconductor manufacturing industry. This development could be an explanation for the rapid onset of n-C6F14 emissions. Additionally, the onset roughly follows the signing of the Montreal Protocol in 1987 and may have triggered the replacement of some CFCs, such as CFC-113, which also had applications in equipment cooling (EPA2008). However, the Cape Grim data in this time period are relatively scarce, and thus the exact timing of the increase in emissions is not well-constrained. By 2013, emission rates have decreased to below 0.2 Gg yr−1 again. This reduction in emissions may partly be attributable to the substitution of n-C6F14 as heat-transfer fluid by hydrofluoroethers, which have some superior properties compared to most liquid PFCs (Tuma and Tousignant2001).

Figure 9Estimated global emission rates for i-C6F14 and n-C6F14 in this study (full red and black lines, respectively) and global emission rates for C6F14 as reported in the EDGARv4.2 database (dashed black line; no distinction between isomers), by Ivy et al. (2012b) (dash-dotted blue line; no distinction between isomers) and by Laube et al. (2012) (dotted black line; no distinction between isomers). The sum of the global emissions for the isomers of C6F14, as reported in the current work, is illustrated by the dotted green line. Shadings illustrate uncertainties.


The effect of separating the C6F14 isomers on the estimated emissions is apparent when comparing them to emission estimates in Laube et al. (2012), where the isomers are not separated (Fig. 9); the emission trend for i-C6F14 is disaggregated from that of the actual n-C6F14 emission trend. Observations presented in Laube et al. (2012) led n-C6F14 emission rate estimates to increase earlier and more gradually and to decline more gradually as well, underestimating the maximum emissions in the late 1990s compared to those in EDGARv4.2 (although they are comparable within the uncertainties) (Fig. 9).

Emission rates for i-C6F14 are estimated to have started increasing in 1992. This is in contrast to 1994, which is what was estimated for the onset of significant n-C6F14 emissions. However, the measurement uncertainties and especially the sparsity of the early data set do not allow for any further conclusions on the exact timing of the onset of emission increases for both isomers. i-C6F14 emissions increase at a slower rate than those for the n-isomer and reach a maximum of 0.25±0.02 Gg yr−1 in 1996–1997, which is when n-C6F14 emissions reach their maximum values as well. i-C6F14 emissions subsequently decrease gradually and stabilise at about 0.09±0.007 Gg yr−1 by 2005. Interestingly, global emissions for i-C6F14 have been constant since ∼2004, but emissions for n-C6F14 seem to continue to decrease until about 2012 (Fig. 5). This is consistent with the increasing i-C6F14 : n-C6F14 ratio in the last two decades (Sect. 3.1.3).

The current study has shown that the i-isomer of C6F14 is clearly present and increasing in the atmosphere (Sect. 3.1.3). The emissions reported in Laube et al. (2012) and Ivy et al. (2012a) (which are based on dilutions of the n-isomer calibrated against a working tank containing compressed, unpolluted air that inevitably contains both isomers) cannot be directly compared to the emissions derived here for the isomers separately. Similarly, due to the fact that EDGARv4.2 does not distinguish PFC isomers, caution must be taken when comparing the emission estimates for the C6F14 isomers from the model to the emissions reported for the C6F14 compound in EDGARv4.2.

Thus, it would be more accurate to compare the sum of the observation-derived emissions of the separate isomers of C6F14 to the reported EDGARv4.2 data. This results in three major improvements. One, it reduces the gap between derived emissions and EDGARv4.2 emissions prior to 1995, as seen for n-C6F14 emission estimates. Two, peak emission rates of the sum of the C6F14 isomers have increased from 1.21±0.1 Gg yr−1 for n-C6F16 alone to 1.46±0.12 Gg yr−1 for the two isomers combined and thereby exceed EDGARv4.2 values. However, given that reported emission rates in EDGARv4.2 inherently also have an uncertainty that is not shown here and that it is more likely that nations under-report emissions, this increase in peak emissions is probably realistic. Three, the estimated emission rates of the sum of both C6F14 isomers agree exceedingly well (relative to all other PFCs reported here) with the EDGARv4.2-reported values from 1999 onwards.

3.2.4 n-C7F16 emissions

Emission rate estimates for n-C7F16 increase from 0.05±0.003 Gg yr−1 in 1980 to a maximum of 0.18±0.01 Gg yr−1 in 1985. Emissions remain stable at this emission level until 2017 (Fig. 10), despite emission-reduction efforts for higher molecular weight PFCs as heat transfer fluids by leak proofing pumps and applying alternative chemistry (EPA2008). This suggests that n-C7F16 either has properties and/or applications that are challenging to substitute and/or unknown atmospheric sources.

Figure 10Estimated global emission rates for n-C7F16 in this study (full black line) and global emission rates for C7F16 as reported in the EDGARv4.2 database (dashed black line; no distinction between isomers), as reported in Ivy et al. (2012b) (dash-dotted blue line; no distinction between isomers) and as reported in Laube et al. (2012) (dotted black line; no distinction between isomers).


Similarly to Laube et al. (2012), model-based emission estimates are higher compared to the EDGARv4.2 reported values before 1990 and after 1999 and are lower in between. Maximum emission rates for n-C7F16 derived by Laube et al. (2012) are 0.23±0.1 Gg yr−1. This slightly higher maximum as compared to the results in this work is likely due to calibration differences as (a) at least two minor isomers of C7F16 are present in the atmosphere and were not chromatographically separated by Laube et al. (2012) and (b) the n-C7F16 calibration in Laube et al. (2012) was not based on a pure isomer, but a technical mixture of isomers. The uncertainties are larger in the work by Laube et al. (2012), due to the additional uncertainty that had to be attributed based on the limited purity of C7F16 used for the calibration.

The estimated emissions have a distinctly different trend compared to those in Ivy et al. (2012b), which seem to be in better agreement with the peak seen in the 1990s in the EDGARv4.2 data. Their larger data set, which also included time series in both the Northern Hemisphere and Southern Hemisphere, is able to better constrain their model. This results in more detailed features in the simulated atmospheric trends. However, to assume that the measured signal for C7F16 in previous work only consists of the n-C7F16 isomer would overestimate n-C7F16 mixing ratios by an unknown amount (see Sect. 2.3), because the signal for C7F16 in the working standard tanks and in the air samples is not exclusive to the n-isomer, as isomers were not separated in the earlier studies.

Observations of the continuing increasing background mixing ratios in the Southern Hemisphere strongly suggest that global emissions of n-C7F16 persist at an approximately constant rate. These conclusions remain consistent with those made in Laube et al. (2012) and Ivy et al. (2012b) for C7F16.

3.3 Possible source regions and source types

The high variability of the PFC mixing ratios in Taiwan described in this work suggests that these measurement sites are close to major PFC sources. The fact that these sampled air masses are not yet well-mixed offers the possibility of investigating the following questions. (1) Which PFCs are likely co-emitted and have similar sources? (2) What are likely source regions for these PFCs within Asia? (3) And what are likely source types of the PFCs measured in Taiwan? To pursue the answers to these questions, the air measurements are combined with the NAME model data, as described in Sect. 2.4.

Table 4An overview of the best Spearman rank correlation values (squared) for all PFCs measured in Taiwan and reported here, including C2F6 and C3F8, with PFC species, regions, and sources. EC: East China, ECS: East China Sea, and SCC: South-central China. A complete overview of correlation coefficients can be found in the Supplement.

Download Print Version | Download XLSX

The similarity of sources among the PFCs can be studied by looking at the inter-species correlations (Fig. S3). These results show that all PFCs, including C2F6 and C3F8 (CF4 was not measured), are significantly correlated with each other, with squared Spearman coefficients ranging from 0.14 to 0.68 (Table S4). The best correlations are observed between C2F6 and C3F8 (in contrast to what was reported by Zhang et al.2017) and between c-C4F8 and n-C4F10 (Table 4). i-C6F14 and n-C6F14 are not very well-correlated, which supports the hypothesis of at least partly independent sources. Some correlation plots seem to show bifurcation within the scattered data (most notably in the correlations of most PFCs with C3F8), suggesting that multiple source signatures may exist (Fig. S3).

When PFC mixing ratios are correlated with the 15 regions used to quantify sources of particle densities used in the dispersion modelling, the R-squared values are generally much smaller, although sample numbers limit the statistics (Table S5). The best correlations are observed for East China for all compounds with the highest squared Spearman rank coefficients found for n-C7F16 (0.49) and n-C5F12 (0.47) (Table 4). Significant correlations are also found for all PFC species with the East China Sea, which can be explained by the fact that air masses from East China travel over the East China Sea before they reach Taiwan.

When using simulated CO mixing ratios as a tracer for source types, all PFCs correlate well with energy industry and domestic sources, which is linked to population density (Table S6). Other best coefficients are found either with the signature of power plants (C2F6, C3F8, c-C4F8, i-C6F14, n-C6F14) or solvents (n-C4F10, n-C5F12, n-C7F16) (Table 4).

The combination of the NAME modelling work and air measurements shows potential to identify and quantify regional PFC sources and source types. However, the interpretation of the results may be challenging, such as the possibility of the heterogeneity in emissions over time (e.g. discontinuous emissions may lead to two similar particle density signatures linked to two very different atmospheric concentrations measured in Taiwan). The collection of air samples at the sites in Taiwan is ongoing. New data in addition to the current data set can be used to inspect PFC source regions and types more rigorously, including the significance of the initial signs of any bifurcation in the correlations between PFCs and CO concentrations that share source types. Nevertheless, this regional study on Taiwan is indicative of the relation between measured atmospheric PFC concentrations and important economic sectors, which feeds back to the implications for rising global PFC trends.

4 Implications
Back to toptop

The changing trends and emissions of six PFCs have been thoroughly discussed in Sects. 3.1 and 3.2, but since these greenhouse gases only occur at very low ppt and sub-ppt levels in our atmosphere, the significance of their influence on our climate system needs further clarification. In order to gain more insight into this matter, the annual global emission rates of all six PFCs covered in this work have been converted to the equivalent of CO2 by using the updated GWPs from the latest IPCC report (Myhre et al.2013) (Table 1). As isomers have the same molecular mass, the GWP for i-C6F14 is estimated based on the GWP reported for n-C6F14 (7910; (Myhre et al.2013) and scaled with the ratio of the radiative efficiencies of i-C6F14 (0.41 W m−2 ppbv−1; Bravo et al.2010) and n-C6F14 (0.44 W m−2 ppbv−1; Myhre et al.2013). PFCs have extremely long atmospheric lifetimes, and thus their CO2 equivalents are calculated here as accumulating over time (Fig. 11) (Laube et al.2012). By 2017, the cumulative global emissions of c-C4F8, n-C4F10, n-C5F12, i-C6F14, n-C6F14, and n-C7F16 amounted to 833 million metric tonnes of CO2 equivalents. This represents an increase in the total global cumulative CO2 equivalent of 23 % between the beginning of 2010 and the end of 2017; 61 % of the total CO2 equivalent is attributable to c-C4F8. The difference between the total CO2 equivalent by the end of 2009 reported in Laube et al. (2012) (∼750 million metric tonnes CO2 equivalent) and in the current work (∼678 million metric tonnes CO2 equivalent) is partly the result of the separation of isomers, but is mostly due to the updated global warming potentials (Myhre et al.2013).

Figure 11Cumulative CO2 equivalent (million metric tonnes) of c-C4F8, n-C4F10, n-C5F12, i-C6F14, n-C6F14, and n-C7F16 based on emissions between 1978 and 2017.


The importance of PFCs extends beyond the scope of climate in the troposphere. The inert nature and continuously increasing concentrations in the troposphere of PFCs over time make them interesting as potential new age-of-air (AoA) tracers. AoA tracers are crucial for research on stratospheric circulation (e.g. Engel et al.2017) and troposphere–stratosphere fluxes (e.g. Bönisch et al.2009), which enables further understanding of threats to the ozone layer. SF6 and CO2 are among the most commonly used AoA tracers, but have limitations that significantly compromise conclusions on dynamic and chemical stratospheric processes (Stiller et al.2012; Ray et al.2017; Leedham Elvidge et al.2018). Some of the lower molecular weight PFCs, such as CF4, C2F6, and C3F8, have already been shown to be suitable age tracers (Leedham Elvidge et al.2018). Considering all PFCs presented in the current work, c-C4F8 has a particular potential to be an excellent age tracer for three reasons. (1) It has a well-constrained time series from 1978 to the present. (2) It still has a sufficiently fast tropospheric growth rate (3.2±0.2 % on average between 2010 and 2018). (3) Observations for background c-C4F8 levels in the Southern Hemisphere have an average measurement uncertainty of about 1 % (1.0±0.7 %). The other PFCs presented in this work are currently not suitable as AoA tracers as their annual growth rates since 2010 (∼0.8–2.9 %) do not exceed their average measurement precisions (∼3.1–7.3 %).

5 Conclusions
Back to toptop

This work has extended an existing time series of observations at Cape Grim, Australia, for atmospheric mixing ratios of c-C4F8, n-C4F10, n-C5F12, n-C6F14, and n-C7F16 from 2010 to 2018 and converted it onto an improved calibration scale. The calibration scale is improved due to the separation of isomers, which leads to a more accurate analysis of tropospheric trends. Moreover, the background trend in i-C6F14 mixing ratios in the Southern Hemisphere is reported here for the first time, and some of the other PFCs have been measured of the first time as discrete isomers.

The background mixing ratios for all six PFCs continue to increase and reached 1.51 ppt (c-C4F8), 0.18 ppt (n-C4F10), 0.15 ppt (n-C5F12), 0.07 ppt (i-C6F14), 0.22 ppt (n-C6F14), and 0.11 ppt (n-C7F16) by the end of 2017. An increasing trend is observed most clearly for atmospheric background concentrations of c-C4F8, i-C6F14, and n-C7F16, which have increased by 27 %, 19 %, and 21 % in the last 8 years, respectively. Atmospheric mixing ratios for n-C4F10, n-C5F12, and n-C6F14 have increased by 9 %, 6 %, and 9 %, respectively.

Southern Hemispheric background trends are compared to mixing ratios in the Northern Hemisphere. Air samples from the Tacolneston site in the UK show an interhemispheric gradient for c-C4F8 and n-C7F16, indicating substantial emissions within the Northern Hemisphere. This is consistent with their faster (as compared to the other four PFCs) increasing atmospheric levels observed in unpolluted, Southern Hemispheric air.

All PFC mixing ratios in Taiwan air samples are extremely variable and show frequent enhancements. These observations show that the Taiwan measurement sites are in close proximity to substantial PFC sources, which is not surprising given the extensive presence of the electronic manufacturing industry in East Asia, which is believed to be one of the main atmospheric sources of these gases (Beu2005; Saito et al.2010).

Finally, emission rate model estimates generally agree well with previous modelling work. Most noteworthy is the continuation of the increasing trend for c-C4F8 emissions. This increase reflects an acceleration of mixing ratio increases in recent years. Our work provides an independent verification of the recent trend of that particular gas and is largely in agreement with the findings of the more extensive c-C4F8-focused work of Mühle et al. (2019). For the longer-chain PFCs, differences with results by Laube et al. (2012) are mainly due to the improved calibration scales. Especially the comparison for C6F14 has improved as its two main isomers are now separated.

The current study demonstrates that for most of the six PFCs in this work, the emissions determined with the bottom-up approach of the EDGARv4.2 database are much lower than those determined by the top-down approach using observations of mixing ratios. This suggests that production and emissions of PFCs are going unreported, which demonstrates the need to better determine the locations of major PFC sources. Analyses of Taiwan PFC mixing ratios and simulated particle dispersion and CO concentrations by the NAME model indicate the potential of such regional studies in understanding the global trends. This analysis will gain statistical weight with more measurements added to the data set each year and is therefore considered ongoing. Voluntary PFC reduction plans exist among semiconductor industry associations (Beu2005). Even though these may be adhered to, PFC emissions are not reported by most nations, including those in East Asia (Saito et al.2010). Determining the contribution of various regions to rising PFC mixing ratios thus remains a challenge.

Monitoring and regulating PFC emissions are relevant to future climate, as illustrated by the cumulative CO2 equivalents of over 830 million tonnes for all PFCs reported here (an increase of 155 million tonnes between the beginning of 2010 and the end of 2017) and the potential to be applied in other areas of atmospheric research, such as circulation and chemistry changes in the stratosphere.

Data availability
Back to toptop
Data availability. 

The measurement and modelling data presented in this work have been made available on Zenodo and can be publicly accessed using (Droste2019) or by contacting

Back to toptop

The supplement related to this article is available online at:

Author contributions
Back to toptop
Author contributions. 

ESD led the manuscript writing process and carried out the emission modelling and some of the measurements. KEA, LJG, AJH, ELE, and JCL also contributed to the measurements, while CER, MJA, ZF, NMH, and MP contributed to the modelling parts of the study. CC, PJF, RLF, SO'D, DEO, CFOY, and WTS contributed through the coordination and execution of the various sampling activities. JCL developed the concept for this study and all the authors contributed to developing it further as well as to the manuscript.

Competing interests
Back to toptop
Competing interests. 

The authors declare that they have no conflict of interest.

Back to toptop

This work was supported by the European Research Council's (ERC) funding for the EXC3ITE project (EXploring stratospheric Chemistry, Composition, and Circulation using Innovative TEchniques). The collection and curation of the Cape Grim Air Archive are jointly funded by CSIRO, the Bureau of Meteorology (BoM), and Refrigerant Reclaim Australia; BoM/CGBAPS staff at Cape Grim were/are largely responsible for the collection of archive samples and UEA flask air samples; the original (mid-1990s) subsampling of the archive for UEA was funded by AFEAS and CSIRO, with ongoing subsampling by CSIRO. Operation of Tacolneston is funded by the Department of Business, Energy & Industrial Strategy (BEIS) through contract TN 1537/06/2018. We would specifically like to thank Stephen Humphrey and Andy MacDonald for all their work in collecting the samples at the Tacolneston site and transporting them to UEA. Taiwan-related work was supported by NERC IOF award and the Academia Sinica. This work used the NAME atmospheric dispersion model, for which the UK Met Office provided the NWP meteorological data sets. Karina E. Adcock's PhD is supported by the Natural Environment Research Council through the EnvEast Doctoral Training Partnership. Norfazrin Mohd Hanif's PhD is supported by the Ministry of Education Malaysia (MOE) and Universiti Kebangsaan Malaysia (UKM). Johannes C. Laube and Lauren J. Gooch received funding from the UK Natural Environment Research Council (Research Fellowship and studentship) and David E. Oram from the National Centre for Atmospheric Science. Lastly, we would like to thank our anonymous referees, whose detailed and constructive comments helped to substantially improve this paper.

Financial support
Back to toptop
Financial support. 

This research has been supported by the European Research Council (grant no. EXC3ITE-678904-ERC-2015-STG), the UK National Environment Research Council (NERC IOF Award NE/J016012/1, Research Fellowship NE/I021918/1, and PhD studentships NE/1210143 and NE/L002582//1), and the Taiwanese Academia Sinica (grant no. AS-KPQ-106-DDPP).

Review statement
Back to toptop
Review statement. 

This paper was edited by Drew Gentner and reviewed by three anonymous referees.

Back to toptop

Adcock, K. E., Reeves, C. E., Gooch, L. J., Leedham Elvidge, E. C., Ashfold, M. J., Brenninkmeijer, C. A. M., Chou, C., Fraser, P. J., Langenfelds, R. L., Mohd Hanif, N., O'Doherty, S., Oram, D. E., Ou-Yang, C.-F., Phang, S. M., Samah, A. A., Röckmann, T., Sturges, W. T., and Laube, J. C.: Continued increase of CFC-113a (CCl3CF3) mixing ratios in the global atmosphere: emissions, occurrence and potential sources, Atmos. Chem. Phys., 18, 4737–4751,, 2018. a, b, c, d, e, f, g, h

Beu, L.: Reductions in Perfluorocompound (PFC) Emissions: 2005 State-of-the-Technology Report, International SEMATECH Manufacturing Initiative, Austin, TX, USA, 2005. a, b, c, d, e, f

Bönisch, H., Engel, A., Curtius, J., Birner, Th., and Hoor, P.: Quantifying transport into the lowermost stratosphere using simultaneous in-situ measurements of SF6 and CO2, Atmos. Chem. Phys., 9, 5905–5919,, 2009. a

Bravo, I., Aranda, A., Hurley, M. D., Marston, G., Nutt, D. R., Shine, K. P., Smith, K., and Wallington, T. J.: Infrared absorption spectra, radiative efficiencies, and global warming potentials of perfluorocarbons: Comparison between experiment and theory, J. Geophys. Res.-Atmos., 115, D24317,, 2010. a, b, c, d, e, f

Deeds, D. A., Vollmer, M. K., Kulongoski, J. T., Miller, B. R., Mühle, J., Harth, C. M., Izbicki, J. A., Hilton, D. R., and Weiss, R. F.: Evidence for crustal degassing of CF4 and SF6 in Mojave Desert groundwaters, Geochim. Cosmochim. Ac., 72, 999–1013,, 2008. a

Droste, E. S.: Trends and Emissions of Six Perfluorocarbons in the Northern and Southern Hemispheres (Version v1), Data set, Zenodo,, 2019. a

EDGAR: EDGAR: Emission Database for Global Atmospheric Research (EDGAR), release version 4.2, European Commission, Joint Research Centre (JRC)/Netherlands Environmental Assessment Agency (PBL), available at: (last access: 31 August 2018), 2014. a, b, c, d, e

Engel, A., Bönisch, H., Ullrich, M., Sitals, R., Membrive, O., Danis, F., and Crevoisier, C.: Mean age of stratospheric air derived from AirCore observations, Atmos. Chem. Phys., 17, 6825–6838,, 2017. a

EPA: Uses and Emissions of Liquid PFC Heat Transfer Fluids from the Electronics Sector, available at: (last access: 24 October 2018), 2008. a, b, c, d

Fleming, Z. L., Monks, P. S., and Manning, A. J.: Untangling the influence of air-mass history in interpreting observed atmospheric composition, Atmos. Res., 104, 1–39,, 2012. a

Fraser, P., Oram, D., Reeves, C., Penkett, S., and McCulloch, A.: Southern Hemispheric halon trends (1978–1998) and global halon emissions, J. Geophys. Res.-Atmos., 104, 15985–15999,, 1999. a, b, c, d

Harnisch, J., Borchers, R., Fabian, P., and Maiss, M.: Tropospheric trends for CF4 and C2F6 since 1982 derived from SF6 dated stratospheric air, Geophys. Res. Lett., 23, 1099–1102,, 1996. a

Hartmann, D., Klein Tank, A., Rusticucci, M., Alexander, L., Brönnimann, S., Charabi, Y., Dentener, F., Dlugokencky, E., Easterling, D., Kaplan, A., Soden, B., Thorne, P., Wild, M., and Zhai, P.: Observations: atmosphere and surface, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, chap. 2, Cambridge University Press, 159–254,, 2013. a, b, c

Ivy, D. J., Arnold, T., Harth, C. M., Steele, L. P., Mühle, J., Rigby, M., Salameh, P. K., Leist, M., Krummel, P. B., Fraser, P. J., Weiss, R. F., and Prinn, R. G.: Atmospheric histories and growth trends of C4F10, C5F12, C6F14, C7F16 and C8F18, Atmos. Chem. Phys., 12, 4313–4325,, 2012a. a, b, c, d, e

Ivy, D. J., Rigby, M., Baasandorj, M., Burkholder, J. B., and Prinn, R. G.: Global emission estimates and radiative impact of C4F10, C5F12, C6F14, C7F16 and C8F18, Atmos. Chem. Phys., 12, 7635–7645,, 2012b. a, b, c, d, e, f, g, h, i

Jones, A., Thomson, D., Hort, M., and Devenish, B.: The UK Met Office's next-generation atmospheric dispersion model, NAME III, in: Air pollution modeling and its application XVII, Springer, 580–589,, 2007. a

Langenfelds, R., Fraser, P., Francey, R., Steele, L., Porter, L., and Allison, C.: The Cape Grim air archive: The first seventeen years, 1978–1995, in: Baseline Atmospheric Program Australia, edited by: Francey, R. J., Dick, A. L., and Derek, N., Bureau of Meteorology and CSIRO Division of Atmospheric Research, Melbourne, Australia, 53–70,, 1996. a, b

Laube, J. C., Martinerie, P., Witrant, E., Blunier, T., Schwander, J., Brenninkmeijer, C. A. M., Schuck, T. J., Bolder, M., Röckmann, T., van der Veen, C., Bönisch, H., Engel, A., Mills, G. P., Newland, M. J., Oram, D. E., Reeves, C. E., and Sturges, W. T.: Accelerating growth of HFC-227ea (1,1,1,2,3,3,3-heptafluoropropane) in the atmosphere, Atmos. Chem. Phys., 10, 5903–5910,, 2010. a, b, c, d, e

Laube, J. C., Hogan, C., Newland, M. J., Mani, F. S., Fraser, P. J., Brenninkmeijer, C. A. M., Martinerie, P., Oram, D. E., Röckmann, T., Schwander, J., Witrant, E., Mills, G. P., Reeves, C. E., and Sturges, W. T.: Distributions, long term trends and emissions of four perfluorocarbons in remote parts of the atmosphere and firn air, Atmos. Chem. Phys., 12, 4081–4090,, 2012. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w, x, y, z, aa, ab, ac, ad, ae, af, ag, ah, ai, aj, ak, al, am, an, ao, ap, aq, ar, as, at, au

Laube, J. C., Mohd Hanif, N., Martinerie, P., Gallacher, E., Fraser, P. J., Langenfelds, R., Brenninkmeijer, C. A. M., Schwander, J., Witrant, E., Wang, J.-L., Ou-Yang, C.-F., Gooch, L. J., Reeves, C. E., Sturges, W. T., and Oram, D. E.: Tropospheric observations of CFC-114 and CFC-114a with a focus on long-term trends and emissions, Atmos. Chem. Phys., 16, 15347–15358,, 2016. a, b, c, d, e

Leedham Elvidge, E., Bönisch, H., Brenninkmeijer, C. A. M., Engel, A., Fraser, P. J., Gallacher, E., Langenfelds, R., Mühle, J., Oram, D. E., Ray, E. A., Ridley, A. R., Röckmann, T., Sturges, W. T., Weiss, R. F., and Laube, J. C.: Evaluation of stratospheric age of air from CF4, C2F6, C3F8, CHF3, HFC-125, HFC-227ea and SF6; implications for the calculations of halocarbon lifetimes, fractional release factors and ozone depletion potentials, Atmos. Chem. Phys., 18, 3369–3385,, 2018. a, b, c

Mazurin, I. M., Stolyarevsky, A. Y., Doronin, A. S., and Shevtsov, A. G.: Ozone-Friendly Working Mixture For Refrigeration Equipment, european Patent 0676460A1, European Patent Office, Torino, Italy, 1994. a, b

Montzka, S., Hall, B., and Elkins, J.: Accelerated increases observed for hydrochlorofluorocarbons since 2004 in the global atmosphere, Geophys. Res. Lett., 36, L03804,, 2009. a

Morris, R. A., Miller, T. M., Viggiano, A., Paulson, J. F., Solomon, S., and Reid, G.: Effects of electron and ion reactions on atmospheric lifetimes of fully fluorinated compounds, J. Geophys. Res.-Atmos., 100, 1287–1294,, 1995. a

Mühle, J., Ganesan, A. L., Miller, B. R., Salameh, P. K., Harth, C. M., Greally, B. R., Rigby, M., Porter, L. W., Steele, L. P., Trudinger, C. M., Krummel, P. B., O'Doherty, S., Fraser, P. J., Simmonds, P. G., Prinn, R. G., and Weiss, R. F.: Perfluorocarbons in the global atmosphere: tetrafluoromethane, hexafluoroethane, and octafluoropropane, Atmos. Chem. Phys., 10, 5145–5164,, 2010. a, b, c

Mühle, J., Trudinger, C. M., Western, L. M., Rigby, M., Vollmer, M. K., Park, S., Manning, A. J., Say, D., Ganesan, A., Steele, L. P., Ivy, D. J., Arnold, T., Li, S., Stohl, A., Harth, C. M., Salameh, P. K., McCulloch, A., O'Doherty, S., Park, M.-K., Jo, C. O., Young, D., Stanley, K. M., Krummel, P. B., Mitrevski, B., Hermansen, O., Lunder, C., Evangeliou, N., Yao, B., Kim, J., Hmiel, B., Buizert, C., Petrenko, V. V., Arduini, J., Maione, M., Etheridge, D. M., Michalopoulou, E., Czerniak, M., Severinghaus, J. P., Reimann, S., Simmonds, P. G., Fraser, P. J., Prinn, R. G., and Weiss, R. F.: Perfluorocyclobutane (PFC-318, c-C4F8) in the global atmosphere, Atmos. Chem. Phys., 19, 10335–10359,, 2019. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o

Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang, J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., Nakajima, T., Robock, A., Stephens, G., Takemura, T., and Zhang, H.: Anthropogenic and natural radiative forcing, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, chap. 8, Cambridge University Press, Cambridge, UK, and New York, USA, 658–740, 2013. a, b, c, d, e

Oram, D., Reeves, C., Penkett, S., and Fraser, P.: Measurements of HCFC-142b and HCFC-141b in the Cape Grim air Archive: 1978–1993, Geophys. Res. Lett., 22, 2741–2744,, 1995. a, b

Oram, D. E., Mani, F. S., Laube, J. C., Newland, M. J., Reeves, C. E., Sturges, W. T., Penkett, S. A., Brenninkmeijer, C. A. M., Röckmann, T., and Fraser, P. J.: Long-term tropospheric trend of octafluorocyclobutane (c-C4F8 or PFC-318), Atmos. Chem. Phys., 12, 261–269,, 2012. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t

Oram, D. E., Ashfold, M. J., Laube, J. C., Gooch, L. J., Humphrey, S., Sturges, W. T., Leedham-Elvidge, E., Forster, G. L., Harris, N. R. P., Mead, M. I., Samah, A. A., Phang, S. M., Ou-Yang, C.-F., Lin, N.-H., Wang, J.-L., Baker, A. K., Brenninkmeijer, C. A. M., and Sherry, D.: A growing threat to the ozone layer from short-lived anthropogenic chlorocarbons, Atmos. Chem. Phys., 17, 11929–11941,, 2017. a, b, c

O'Shea, S. J., Choularton, T. W., Flynn, M., Bower, K. N., Gallagher, M., Crosier, J., Williams, P., Crawford, I., Fleming, Z. L., Listowski, C., Kirchgaessner, A., Ladkin, R. S., and Lachlan-Cope, T.: In situ measurements of cloud microphysics and aerosol over coastal Antarctica during the MAC campaign, Atmos. Chem. Phys., 17, 13049–13070,, 2017. a

Ray, E. A., Moore, F. L., Elkins, J. W., Rosenlof, K. H., Laube, J. C., Röckmann, T., Marsh, D. R., and Andrews, A. E.: Quantification of the SF6 lifetime based on mesospheric loss measured in the stratospheric polar vortex, J. Geophys. Res.-Atmos., 122, 4626–4638,, 2017. a

Reeves, C. E., Sturges, W. T., Sturrock, G. A., Preston, K., Oram, D. E., Schwander, J., Mulvaney, R., Barnola, J.-M., and Chappellaz, J.: Trends of halon gases in polar firn air: implications for their emission distributions, Atmos. Chem. Phys., 5, 2055–2064,, 2005. a, b, c

Riahi, K., Rao, S., Krey, V., Cho, C., Chirkov, V., Fischer, G., Kindermann, G., Nakicenovic, N., and Rafaj, P.: RCP 8.5 – A scenario of comparatively high greenhouse gas emissions, Climatic Change, 109, 33–57,, 2011. a

Robin, L. R. and Iikubo, Y.: Fire Extinguishing Methods Utilizing Perfluorocarbons, US Patent 5117917, US Patent Office, USA, 1992. a, b

Saito, T., Yokouchi, Y., Stohl, A., Taguchi, S., and Mukai, H.: Large emissions of perfluorocarbons in East Asia deduced from continuous atmospheric measurements, Environ. Sci. Technol., 44, 4089–4095,, 2010. a, b, c

Stanley, K. M., Grant, A., O'Doherty, S., Young, D., Manning, A. J., Stavert, A. R., Spain, T. G., Salameh, P. K., Harth, C. M., Simmonds, P. G., Sturges, W. T., Oram, D. E., and Derwent, R. G.: Greenhouse gas measurements from a UK network of tall towers: technical description and first results, Atmos. Meas. Tech., 11, 1437–1458,, 2018. a

Stiller, G. P., von Clarmann, T., Haenel, F., Funke, B., Glatthor, N., Grabowski, U., Kellmann, S., Kiefer, M., Linden, A., Lossow, S., and López-Puertas, M.: Observed temporal evolution of global mean age of stratospheric air for the 2002 to 2010 period, Atmos. Chem. Phys., 12, 3311–3331,, 2012. a

Trudinger, C. M., Fraser, P. J., Etheridge, D. M., Sturges, W. T., Vollmer, M. K., Rigby, M., Martinerie, P., Mühle, J., Worton, D. R., Krummel, P. B., Steele, L. P., Miller, B. R., Laube, J., Mani, F. S., Rayner, P. J., Harth, C. M., Witrant, E., Blunier, T., Schwander, J., O'Doherty, S., and Battle, M.: Atmospheric abundance and global emissions of perfluorocarbons CF4, C2F6 and C3F8 since 1800 inferred from ice core, firn, air archive and in situ measurements, Atmos. Chem. Phys., 16, 11733–11754,, 2016.  a

Tsai, W.-T., Chen, H.-P., and Hsien, W.-Y.: A review of uses, environmental hazards and recovery/recycle technologies of perfluorocarbons (PFCs) emissions from the semiconductor manufacturing processes, J. Loss Prevent. Proc., 15, 65–75,, 2002. a

Tuma, P. E. and Tousignant, L.: Reducing emissions of PFC heat transfer fluids, in: SEMI Technical Symposium, 3M Specialty Materials, 16 July 2001, San Francisco, CA, USA, 1–8, 2001. a

Van Vuuren, D. P., Edmonds, J., Kainuma, M., Riahi, K., Thomson, A., Hibbard, K., Hurtt, G. C., Kram, T., Krey, V., Lamarque, J.-F., Masui, T., Meinshausen, M., Nakicenovic, N., Smith, S. J., and Rose, S. K.: The representative concentration pathways: an overview, Climatic Change, 109, 5–31,, 2011. a

Worton, D. R., Sturges, W. T., Gohar, L. K., Shine, K. P., Martinerie, P., Oram, D. E., Humphrey, S. P., Begley, P., Gunn, L., Barnola, J.-M., Schwander, J., and Mulvaney, R.: Atmospheric trends and radiative forcings of CF4 and C2F6 inferred from firn air, Environ. Sci. Technol., 41, 2184–2189,, 2007. a

Zhang, G., Yao, B., Vollmer, M. K., Montzka, S. A., Mühle, J., Weiss, R. F., O'Doherty, S., Li, Y., Fang, S., and Reimann, S.: Ambient mixing ratios of atmospheric halogenated compounds at five background stations in China, Atmos. Environ., 160, 55–69,, 2017. a

Zhang, H., Wu, J., and Shen, Z.: Radiative forcing and global warming potential of perfluorocarbons and sulfur hexafluoride, Science China Earth Sciences, 54, 764–772,, 2011. a

Publications Copernicus
Short summary
We update the tropospheric trends and emissions of six perfluorocarbon (PFC) gases, including separate isomers. Trends for these strong greenhouse gases are still increasing, but at slower rates than previously. The lack of natural sinks results in the global accumulation of 833 million metric tonnes of CO2 equivalent for these six PFCs by 2017. Modelling results indicate potential source regions and types in East Asia, but we find that many emissions are unaccounted for in emission reports.
We update the tropospheric trends and emissions of six perfluorocarbon (PFC) gases, including...