2.1 Measurement sites

The details of the observational sites used, such as their coordinates, their altitude, the air inlet height above ground, and the height of each site in the different transport model versions (Sect. 2.4), are summarized in Table 1. Their location can be seen in Fig. 1. Since the main goal of this study is to quantify the differences between low- (7 km) and high-resolution inversions (1 km) in Switzerland, the observational sites chosen should be sensitive to Swiss emissions. Most of the Swiss F-gas emissions can be expected to originate from the region called the Swiss Plateau. It is located north of the Alps, covering about $\mathrm{1}/\mathrm{3}$ of the area of Switzerland and including about $\mathrm{2}/\mathrm{3}$ of the population of Switzerland. The biggest cities of Switzerland are located in this region and most of the industrial activity takes place here as well.

Table 1Details of the observational sites used in the study, including the location, altitude, and height of the model topography in the different FLEXPART model versions.

Download Print Version | Download XLSX

Figure 1Measurement locations (black crosses) of sites used in the inversions as well as COSMO and IFS model domains (polygons) used as input to FLEXPART (a). Swiss measurement locations (black crosses) and major cities (red circles) on top of topographic relief including major rivers and lakes (Swiss coordinate system, LV03) (b).

The Beromünster (BRM) tall tower site (Table 1) is located in the middle of the Swiss Plateau on a hill with an elevation of about 800 m a.s.l. GHG measurements at the tower were established in 2012 (Berhanu et al., 2016; Oney et al., 2015) and since 2016 the site has been part of the Swiss air quality observing network (NABEL). The area surrounding Beromünster is mainly rural and used for agricultural activities. BRM is sensitive to emissions from most of the Swiss Plateau, as can be seen in Fig. 2. The closest city to Beromünster is Lucerne (urban area population of 220 000), located 20 km south of the site, whereas the Zurich urban area (approximately 1.3 million) is 40 km to the east. The measurements at the site were taken on a tall tower with a height of 217 m a.g.l. at a sample inlet height of 212 m a.g.l.

Figure 2Simulated total surface sensitivity (footprints) for Beromünster (a, b) and for Sottens (c, d) for the duration of the measurement campaigns (1 September 2019–31 August 2020 for BRM and 5 March 2021–24 October 2021 for SOT) as obtained from FLEXPART-COSMO-7 (a, c) and FLEXPART-COSMO-1 (b, d). The surface sensitivity is given as particle residence time per air density. The locations of the sites are indicated with black crosses, and major cities are indicated with black circles.

The Sottens (SOT, Table 1) tall tower is located in the western part of the Swiss Plateau in the canton of Vaud at an altitude of about 800 m a.s.l. The area surrounding SOT is also rural. Measurements in SOT are sensitive to emissions from the western and center parts of the Swiss Plateau as well as from the deep Alpine Rhone valley (canton of Valais, Fig. 2). The closest larger city is Lausanne (urban area population of 430 000), located 20 km southwest of the site, while Geneva is about 80 km west-southwest of the site (urban area population of 500 000). The measurements on the tall tower were taken at an inlet height of 120 m a.g.l.

The Jungfraujoch (JFJ, Table 1) site is a high-altitude observatory located in the Bernese Alps on the border between the cantons of Valais and Bern. The observatory is located at a steep mountain saddle connecting two major mountains (Jungfrau at 4158 m a.s.l. and Mönch at 4099 m a.s.l.). JFJ is part of the AGAGE network and has been measuring halocarbons since 2000. Although JFJ is representative of lower-free-tropospheric conditions in the winter, it frequently receives fresh boundary layer pollution during the summer months from both the Swiss Plateau and the south of the Alps, but also from more distant sources throughout central Europe (Henne et al., 2010; Herrmann et al., 2015).

Two additional sites were used in the inversions conducted within this study: Mace Head (MHD) and Tacolneston (TAC), as shown in Table 1. MHD is located in County Galway on the west coast of Ireland. Its exposure to the North Atlantic Ocean makes it an ideal location for background observations due to the dominating westerly flow. TAC is located 150 km to the northeast of London on the east coast of England in southern Norfolk. Its location is optimal to constrain emissions from the UK and partly from the Benelux region, which is one of the regions with the highest emission density in Europe (Manning et al., 2021). These sites were also used in a previous study to constrain Swiss halocarbon emissions (Rust et al., 2022). Adding these sites outside Switzerland allows the inversion to constrain larger-scale European emissions. Leaving these unconstrained by observations may have led to biased estimates for the Swiss domain.

Two more AGAGE sites were employed in sensitivity inversions to explore any further impact of additional observations on Swiss emissions: Monte Cimone (CMN) and Taunus Observatory (TOB), as shown in Table 1. CMN is a high-altitude observatory on the highest peak of the northern Apennines in Italy. Its remote location, high altitude, and large distance from big cities and hence major emission sources make it representative of the free troposphere and background values in southern Europe and the northern Mediterranean basin, but it can also occasionally receive pollution events from the Po Valley (Bonasoni et al., 2000). TOB is located on the second-highest peak in the Taunus mountain region in central Germany. Its close proximity to major emission sources (Frankfurt and Mainz) and its location in central Europe make the site well suited for European air pollution studies.

2.2 Observational data

In this study, we used observational data on 1,1,1,2-tetrafluoroethane (HFC-134a), 1,1,1,2,2-pentafluoroethane (HFC-125), difluoromethane (HFC-32), and sulfur hexafluoride (SF₆), which together with 1,1,1-trifluoroethane (HFC-143a, not reliably measured from BRM and SOT) account for more than 80 % of total Swiss halocarbon emissions in terms of CO₂ equivalents (Reimann et al., 2021). The observational data come from two extended measurement campaigns at BRM and SOT, conducted within the project IHALOME (Rust et al., 2022), and from the AGAGE monitoring network (Prinn et al., 2018). The measurement campaigns at BRM and SOT were performed to explore the impact of a denser measurement network on the Swiss national emission estimates by top-down methods. Semicontinuous air samples were taken with a frequency of approximately two ambient air measurements within 3 h using a Medusa pre-concentration unit, coupled to gas chromatography (GC) and mass spectrometry (MS) (Miller et al., 2008). In total, about 60 fully calibrated halocarbons were measured with atmospheric abundances in the dry-air mole fraction range of parts per quadrillion (femtomol mol⁻¹) to parts per trillion (picomol mol⁻¹; Rust et al., 2022). Measurements of HFC-143a suffered from instrumental problems and were therefore not used in the present analysis. Measurements at MHD and TAC were also performed with a Medusa GC instrument, while those at TOB and CMN were performed with different pre-concentration and GC-MS instruments (Schuck et al., 2018; Maione et al., 2013). The measurements used in the present study are based on fully intercalibrated reference standards. The measurements are based on the Scripps Institution of Oceanography (SIO) primary calibration scales SIO-05 for HFC-134a and SF₆, SIO-14 for HFC-125, and SIO-07 for HFC-32.

The observational data employed for the inversions cover the period from August 2019 to October 2021. Data from TAC, MHD, and JFJ were used for the whole period, while BRM and SOT data were available only during the field campaigns. The campaign in BRM lasted from August 2019 to September 2020 and in SOT from March 2021 to October 2021. There was no temporal overlap because the same instrument had to be used at both locations. CMN and TOB observations were employed in sensitivity inversions for the BRM campaign period only. Measurements from TOB come from flask samples, which are collected weekly for offline analysis.

2.3 Baseline

To run our inversions, 24-hourly (3-hourly for sensitivity inversion) mean values were produced from the available observations of the abovementioned sites. To correctly infer regional emissions from a limited model domain, accurate knowledge of the so-called background (or baseline) mole fraction of a compound is needed. An observed mole fraction of a compound can be decomposed into a baseline fraction, y_o,b, and the contribution due to recent emissions, as targeted by the regional simulation, y_o,p:

$$\begin{array}{}\text{(1)}& y_{\mathrm{o}}=y_{\mathrm{o},\mathrm{b}}+y_{\mathrm{o},\mathrm{p}}.\end{array}$$

An underestimation of the baseline will magnify an emission event, whereas an overestimation will reduce the intensity of an emission event. We estimate our baseline mole fractions by using the robust extraction of the baseline signal (REBS) method (Ruckstuhl et al., 2012). The REBS method is an iterative filter, which assumes that the mean of the baseline can be approximated by a smooth curve and its uncertainty distribution can be given by a Gaussian distribution with a constant (in time) standard deviation. The smooth curve is estimated by piecewise local weighted linear regressions. The first weight function acts to decrease the impact of observations in proportion to their distance from the point of the time series that is to be estimated, t₀. The second reduces the influence of the data according to their distance from the expected value of the baseline. In each iteration, an updated estimate of the baseline value and the baseline standard deviation is calculated, and the method usually converges after 5–10 iterations. For our baseline estimates, a tuning factor of b=3.5, a temporal window width of 60 d, and a maximum of 10 iterations were used. In our inversions, we used the baseline estimated from JFJ for all Swiss sites, the one estimated for MHD for the sites on the British Isles and TOB, and the one estimated for CMN for CMN itself. This selection is motivated by the fact that the REBS method works best for sites that are mostly sampling background, whereas for typical continental boundary layer sites (like BRM, SOT, TAC) only a few “pure” background observations exist throughout the year, and hence REBS-estimated baselines tend to overestimate true baselines. All baselines were updated as part of the emission inversion step for each site individually.

Other statistical methods for baseline estimation have been applied to greenhouse gas observations (e.g., Thoning et al., 1989; El Yazidi et al., 2018) some of which use additional transport model information (trajectories or footprints) to select background sectors (O'Doherty et al., 2001). Differences between estimated background conditions are often small or limited to certain events or situations. There is no consensus regarding which of these methods is most robust under all circumstances.

2.4 Transport models

The inversion system utilized for this study is comprised of an atmospheric transport model, which relates the spatial emissions, x, of the compound of interest to the mole fractions measured at the receptor site, y_o, via a linear mapping y_o=ℋx. Here, the LPDM FLEXPART (Pisso et al., 2019) was driven by the meteorological fields from two Eulerian NWP models: the limited-area NWP model COSMO and the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). Simulations with IFS were used to extend FLEXPART-COSMO simulations beyond the COSMO model domain (Katharopoulos et al., 2022).

2.4.2 FLEXPART LPDM

FLEXPART was initially designed to estimate the mesoscale and synoptic dispersion of radio-nuclei from point sources, such as releases during a nuclear accident like Chernobyl. Nowadays, FLEXPART (Stohl et al., 2005; Pisso et al., 2019) and other LPDMs like STILT, NAME, and HYSPLIT (Lin et al., 2003; Jones et al., 2007; Stein et al., 2015b) are utilized for a large variety of tracer transport problems, simulating the transport, diffusion, conversion, and deposition of various compounds ranging from inert GHGs to aerosol particles.

One of FLEXPART's major applications is in inverse modeling studies for the estimation of regional- and continental-scale emissions of atmospheric compounds (Fang and Michalak, 2015; Henne et al., 2016; Brunner et al., 2012; Stohl et al., 2010). This is due to FLEXPART's ability for both forward- and backward-in-time simulations. For backward simulations, particle trajectories are integrated backward in time using a negative time step. The final product is an estimate of the sensitivity of a concentration measured at the receptor y_i to an emission source x_i, called the source–receptor relationship (Seibert and Frank, 2004). Source–receptor relationships derived from FLEXPART are linear since all atmospheric processes considered during the transport of the tracers are linear (advection, diffusion, convective mixing). The compounds we are interested in possess very long atmospheric lifetimes (5 years or more), so for the regional-scale transport (less than 10 d) we can assume these to be inert. Thus, linear relationships, m_i,l, in units of s m³ kg⁻¹ mol mol⁻¹ with i referring to different grid cells and l referring to different receptors, can be easily derived from FLEXPART. If the spatial distribution of emissions E_i is multiplied by source sensitivities, m_i,l, the product yields the mixing ratio increment, y_l, of the tracer at the receptor site, l, resulting from emissions in the considered domain and time window,

$$\begin{array}{}\text{(2)}& y_l=\sum _i{m}_{i,l}{E}_i,\end{array}$$

to which the baseline concentration y_b,l needs to be added to obtain the absolute mixing ratio. In our case, y_b,l was estimated from observations using the REBS method (Sect. 2.3).

Here, we utilize two versions of FLEXPART in backward mode: FLEXPART-COSMO and FLEXPART-IFS. FLEXPART-COSMO is a version of FLEXPART adapted to the COSMO model (Henne et al., 2016). The meteorological fields driving FLEXPART are directly used in the hybrid-z coordinate system of COSMO with no additional interpolation. FLEXPART-IFS interpolates the meteorological fields from the hybrid pressure coordinate system of ECMWF-IFS to a terrain-following z-based system (Stohl et al., 2005). The meteorological fields employed in FLEXPART simulations are some of the driving NWP's prognostic variables (winds, temperature, pressure, etc.) and accumulated fluxes (precipitation, surface heat, momentum, moisture).

FLEXPART-COSMO is employed at two different spatial resolutions, 7 and 1 km (Sect. 2.4). When we refer to spatial resolution, we always mean the resolution of the driving NWP, here COSMO, and not the LPDM. In the Lagrangian framework, there is no discretization of space and the frame of reference is centered on each particle following its trajectory in the space–time continuum.

Receptor-oriented FLEXPART simulations were carried out by releasing 50 000 particles at each different receptor continuously over 3 h periods. Particles were then traced back for 8 d for FLEXPART-COSMO-7 and for 4 and 8 d for FLEXPART-COSMO-1 coupled (see below) to FLEXPART-IFS, respectively. Source sensitivities were stored on two different output domains for FLEXPART-COSMO simulations: a larger domain (main, 0.16^∘ × 0.12^∘ horizontal resolution) covering a similar area as the COSMO-7 simulations and a smaller and finer domain (nest, 0.02^∘ × 0.015^∘ horizontal resolution) focusing on Switzerland. For simulations with FLEXPART-IFS, the output grid was at a lower horizontal resolution grid for the whole of Europe (0.1^∘ × 0.1^∘).

If FLEXPART is driven only by COSMO-1 fields, source sensitivities can only be produced for the limited COSMO-1 domain, and any European contributions from larger distances (as from the COSMO-7 domain) would be neglected. To account for this limitation, we offline-nest FLEXPART-COSMO-1 to FLEXPART-IFS in order to continue the integration of the particles in Europe once they leave the COSMO-1 domain (Katharopoulos et al., 2022).

Particle transport in FLEXPART is modeled by a simple zero acceleration scheme,

$$\begin{array}{}\text{(3)}& {\displaystyle }\mathit{X}(t+\mathrm{d}t)=\mathit{X}\left(t\right)+\mathit{u}(\mathit{X},t)\mathrm{d}t,\text{(4)}& {\displaystyle }\mathit{u}(\mathit{X},t)={\mathit{u}}_{\mathrm{g}}(\mathit{X},t)+{\mathit{u}}^{\mathbf{\prime }}(\mathit{X},t)+{\mathit{u}}_{\mathrm{m}}(\mathit{X},t),\end{array}$$

where X​​​​​​​ is the particle's position, and u is the wind vector at the particle's location comprised of three components. The term u_g is the average wind vector at the particle location (in our case taken from the COSMO model), u^′ is the fluctuation from the mean wind representing the turbulence in the atmosphere (modeled as a stochastic Markov chain process), and u_m represents additional mesoscale wind variations (Stohl et al., 2005).

For simulations with FLEXPART-COSMO-7, we use the original turbulence parameterization of FLEXPART, the Hanna turbulence scheme (Stohl et al., 2005), while for FLEXPART-COSMO-1 we utilize the novel scheme introduced by Katharopoulos et al. (2022). We recently showed that since the Hanna scheme is developed to parameterize the whole turbulence spectrum and COSMO-1 wind fields explicitly resolve part of the turbulence spectrum, eddies of the size of the model grid and bigger can be represented by the wind fields – that leads to duplication of parts of the turbulence spectrum in the model, and, as a result, to increased diffusion.

2.5 Inversion framework

As we have already mentioned, FLEXPART was utilized in the backward mode to produce source sensitivities, M, which translate spatial emissions, x, to mole fractions, y, at the receptor site:

$$\begin{array}{}\text{(5)}& \mathit{y}=\mathbf{M}\mathit{x}.\end{array}$$

The state vector x, $\mathit{x}=(x_{\mathrm{1}},\mathrm{\dots },x_k{)}^T$ contained K elements, which correspond to the sum of the total number of grid cells, N_E, in our inversion grid and the total number of baseline nodes, N_B, to be optimized by the inversion. The baseline nodes are baseline factors at discrete time intervals, since we do not optimize the baseline at every time step to reduce the size of the matrix, M, and also avoid ending up with an underdetermined system of equations. Here, the time interval between our baseline nodes was set to τ_b=30 d. The estimation of the a priori baseline is described in Sect. 2.3. The rectangular matrix M (size L×K) is a column block matrix with two blocks, M_E and M_B, representing the sensitivity of the observations to emissions for each grid cell and the baseline mole fractions, respectively. The mole fractions at the receptor sites are the product of the sensitivity matrix with the state vector, and its length is equal to the number of observations at all receptor sites, $\mathit{y}=(y_{\mathrm{1}},\mathrm{\dots },y_L{)}^T$.

If we used the complete output grid of our transport model as the inversion grid, then the size of our sensitivity matrix would be too large to be computationally manageable and the solution probably would be underdetermined depending on the spatial correlation lengths. Fine grids with negligible source sensitivities and very low a priori emissions are also more prone to be assigned negative emissions in typical dipole patterns since we assume Gaussian-distributed errors. To reduce the size of the inversion problem, an irregularly sized inversion grid is introduced that assigns finer (lower) grid cells in areas with larger (smaller) average source sensitivities (Henne et al., 2016). The reduced-resolution grid serves two purposes. On the one hand, it reduces the number of state vector elements, removing many elements with very little sensitivity. On the other hand, it helps to smooth out transport model errors, which tend to grow with distance from the point of observation. The number of grid cells in our inversions varies from 1000 to 2500 depending on the number of observations available for different inversions.

Bayesian inverse modeling is employed to statistically optimize the estimates of the variables of interest, x, by constraining them with the observational data, y⁰ (top-down constraint), and with the prior estimate of the variables of interest, x_b (bottom-up constraint). Gaussian-distributed errors are always assumed between the observations and the simulated mole fractions and between the a priori and the a posteriori emissions,

$$\begin{array}{}\text{(6)}& {\displaystyle }P\left(\mathit{x}\right)={\displaystyle \frac{\mathrm{1}}{\sqrt{\mathrm{2}\mathit{\pi }\left|\mathbf{B}\right|}}}{e}^{-\frac{\mathrm{1}}{\mathrm{2}}{\left(\mathit{x}-{\mathit{x}}_{\mathrm{b}}\right)}^T{\mathbf{B}}^{-\mathrm{1}}\left(\mathit{x}-{\mathit{x}}_{\mathrm{b}}\right)},\text{(7)}& {\displaystyle }P\left({\mathit{y}}^{\mathrm{0}}\right|\mathit{x})={\displaystyle \frac{\mathrm{1}}{\sqrt{\mathrm{2}\mathit{\pi }\left|\mathbf{R}\right|}}}{e}^{-\frac{\mathrm{1}}{\mathrm{2}}{\left({\mathit{y}}^{\mathrm{0}}-\mathbf{M}\mathit{x}\right)}^T{\mathbf{R}}^{-\mathrm{1}}\left({\mathit{y}}^{\mathrm{0}}-\mathbf{M}\mathit{x}\right)},\end{array}$$

where B is the a priori error covariance matrix, and R is the observational error covariance matrix. The construction of these matrices is discussed in Sect. 2.5.1. By applying Bayes' theorem, $P\left(\mathit{x}\right|\mathit{y})\approx P(\mathit{y}\left|\mathit{x}\right)P\left(\mathit{x}\right)$, we obtain the a posteriori Gaussian probability distribution function for the error of the emissions. The cost function in Eq. (8) is the negative logarithm of P(x|y). We minimize the cost function to find the value of the state of the emissions, x, that minimizes the observational and a priori error:

$$\begin{array}{}\text{(8)}& \begin{array}{rl}J\left(\mathit{x}\right)& ={\left({\mathit{y}}^{\mathrm{0}}-\mathbf{M}\mathit{x}\right)}^T{\mathbf{R}}^{-\mathrm{1}}\left({\mathit{y}}^{\mathrm{0}}-\mathbf{M}\mathit{x}\right)\\ & +{\left(\mathit{x}-{\mathit{x}}_{\mathrm{b}}\right)}^T{\mathbf{B}}^{-\mathrm{1}}\left(\mathit{x}-{\mathit{x}}_{\mathrm{b}}\right).\end{array}\end{array}$$

The minimization problem can be solved analytically, since the sensitivity matrix, Mx, is a linear mapping, $\mathit{x}\in {\mathbb{R}}^K\to \mathbf{M}\mathit{x}\in {\mathbb{R}}^L$. Major advantages of the analytical approach are (1) the complete characterization of the a posteriori error (Brasseur and Jacob, 2017) as part of the solution and (2) that it can be fast and well suited for a plethora of sensitivity inversions. The minimization of the cost function yields the solution

$$\begin{array}{}\text{(9)}& \widehat{\mathit{x}}={\mathit{x}}_{\mathrm{b}}+\mathbf{G}\left(\mathit{y}-\mathbf{M}{\mathit{x}}_{\mathrm{b}}\right),\end{array}$$

where G is the gain matrix,

$$\begin{array}{}\text{(10)}& \mathbf{G}={\mathbf{BM}}^T{\left({\mathbf{MBM}}^T+\mathbf{R}\right)}^{-\mathrm{1}},\end{array}$$

giving the sensitivity of the optimal state to the observations. In the analytical inversion, the a posteriori error covariance matrix can be directly calculated as

$$\begin{array}{}\text{(11)}& \widehat{\mathbf{S}}=\mathbf{B}-\mathbf{BMG},\end{array}$$

describing the uncertainty of the posterior estimate.

2.5.1 Covariance matrices

Our design of the error covariance matrices, B and R, follows a maximum likelihood approach for which initial estimates of the matrices are needed (Henne et al., 2016). Both covariance matrices are symmetric block matrices. The observational error covariance matrix, $\mathbf{R}=\left[{\mathit{ϵ}}_{\mathrm{o}}{\mathit{ϵ}}_{\mathrm{o}}^T\right]$, contains contributions from the instrument error, ϵ_I, the representation error, ϵ_R, and the model error, ϵ_M. These errors are assumed to be uncorrelated, so the covariance matrix can be calculated as the sum of squares of the individual covariance matrices for each source of error, ${\mathit{ϵ}}_{\mathrm{o}}^{\mathrm{2}}={\mathit{ϵ}}_{\mathrm{I}}^{\mathrm{2}}+{\mathit{ϵ}}_{\mathrm{R}}^{\mathrm{2}}+{\mathit{ϵ}}_{\mathrm{M}}^{\mathrm{2}}$.

The block matrix R is a row block matrix, containing a number of blocks equal to the number of different receptors. Diagonal elements of R are estimated as follows:

$$\begin{array}{}\text{(12)}& R_{i,i}={\mathit{ϵ}}_{\mathrm{I}}^{\mathrm{2}}+\mathit{\alpha }+\mathit{\beta }{\mathit{\chi }}_{\mathrm{p},i}^{\mathrm{2}}.\end{array}$$

Representation and model errors are considered to be a single error, increasing linearly with a-priori-simulated mixing ratios, χ_p,i. The factors α and β are determined by the log-likelihood approach. For each block matrix representing an individual receptor, temporal correlation in the error is added to the covariance matrix by setting the non-diagonal entries to

$$\begin{array}{}\text{(13)}& R_{i,j}=e^{-\frac{T_{i,j}}{{\mathit{\tau }}_{\mathrm{0}}}}\sqrt{R_{i,i}}\sqrt{R_{j,j}}i\ne j.\end{array}$$

The factor T_i,j is the time difference between measurements, and τ₀ is the temporal correlation length, here set to a very small value of 0.01 d, meaning that we assume almost zero auto-correlation (independence) between daily average observation–model errors. Although this may underestimate the true error correlation in some situations, in our experience it allows capturing pronounced pollution events more realistically in the a posteriori simulations. Additionally, error correlation between different sites is neglected.

The matrix B consists of two block matrices. The first corresponds to the emissions, B^E, and the second to the baseline, B^B. The diagonal elements of matrix B^E are proportional (factor f_E) to the a priori emissions, while the off-diagonal elements are spatially correlated. The correlation fades as an exponential function of their distance, d_i,j, scaled by a correlation length scale, L:

$$\begin{array}{}\text{(14)}& {\displaystyle }{B}_{i,i}^{\mathrm{E}}={\left(f_{\mathrm{E}}{x}_{\mathrm{b},i}\right)}^{\mathrm{2}},\text{(15)}& {\displaystyle }{B}_{i,j}^{\mathrm{E}}=e^{-\frac{d_{i,j}}{L}}\sqrt{B_{i,i}^{\mathrm{E}}}\sqrt{B_{j,j}^{\mathrm{E}}}i\ne j.\end{array}$$

The diagonal values of block matrix B^B are proportional (factor f_b) to the baseline error, while the non-diagonal elements are set to be correlated in time. The correlation fades as an exponential function of the time difference between baseline nodes, T_i,j, scaled by a temporal correlation length, τ_b:

$$\begin{array}{}\text{(16)}& {\displaystyle }{B}_{i,i}^{\mathrm{B}}={\left(f_{\mathrm{b}}{\mathit{\sigma }}_{\mathrm{b}}\right)}^{\mathrm{2}},\text{(17)}& {\displaystyle }{B}_{i,i}^{\mathrm{B}}=e^{-\frac{T_{i,j}}{{\mathit{\tau }}_{\mathrm{b}}}}\sqrt{B_{i,i}^{\mathrm{B}}}\sqrt{B_{j,j}^{\mathrm{B}}}i\ne j.\end{array}$$

2.5.2 Maximum likelihood

Accurate knowledge of the a priori and observational error covariance matrices, B and R, is, in general, unavailable and often “expert judgments” are used to estimate or set the parameters describing the matrices. Similarly, our initial values of the covariance matrices are a mix of expert judgment and methods used in the literature (Henne et al., 2016). To overcome the partial subjectivity of the construction of the covariance matrices, we employ a maximum likelihood optimization step (Michalak et al., 2005). The parameters that we optimize in the maximum likelihood optimization are the correlation length, L, the factor f_E, which gives the variance of the emissions at each grid cell relative to the prior emissions, and the temporal correlation length of the baseline, τ_B. Additionally, for each different receptor, the factors f_b, α, and β are optimized. The maximum likelihood estimate of the covariance parameters is obtained by minimizing Eq. (18) with respect to the covariance parameters (Michalak et al., 2005),

$$\begin{array}{}\text{(18)}& \begin{array}{rl}{L}_{\mathit{\theta }}& ={\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}\ln \left|{\mathbf{MBM}}^T+\mathbf{R}\right|\\ & +{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}{\left(\mathit{y}-\mathcal{H}{\mathit{x}}_{\mathrm{b}}\right)}^T{\left({\mathbf{MBM}}^T+\mathbf{R}\right)}^{-\mathrm{1}}\left(\mathit{y}-\mathbf{M}{\mathit{x}}_{\mathrm{b}}\right).\end{array}\end{array}$$

2.6 Sensitivity tests

The main focus of this study is to assess the impact of high-resolution FLEXPART-COSMO-1 simulations on the emission estimates of halocarbons. The transport model resolution is one of the factors which can influence the total inverse emission estimates, their spatial distribution, and their uncertainty. The kind of analytical inversion used here to optimize the emissions was shown to likely underestimate the uncertainty of the a posteriori state vector (e.g., Berchet et al., 2015). To capture the whole range of uncertainty of our a posteriori, we conducted additional sensitivity tests (Table 2), in which we vary different parameters and aspects of our inversion (transport model, inversion grid, spatial distribution of a priori emissions, optimization of different covariance parameters during the maximum likelihood estimation step, temporal resolution of the assimilated observations, sensitivity of the inversion to the inclusion of observations from additional sites, seasonality of emissions). In the following, our BASE inversion (Table 2) corresponds to inversions for which BRM, SOT, and JFJ are employed as the observational sites in Switzerland. TAC and MHD are used in this setup as additional non-Swiss observational sites. Furthermore, the maximum likelihood step is calculated for all covariance parameters except L and f_b, which are fixed to specific values for each compound. The observations are aggregated over 24 h intervals and the irregular grid size is increased to approximately 2000 grid cells, compared to the inversions with fewer cells in Rust et al. (2022). This setup is used for comparisons across inversions with different transport model resolutions for the same tracer (BASE1 and BASE7, Table 2). All the different sensitivity tests described in the following sections are summarized in Table 2.

Table 2Different groups of inversions conducted in this study.

^∗ Two different approaches for setting covariance parameters were used for the seasonal inversions. See text for details.

Download Print Version | Download XLSX

2.6.2 Spatial distribution of a priori emissions

We conducted inversions using three different spatial distributions of the a priori emissions, x_b, in order to test the sensitivity of the a posteriori estimated vector and its uncertainty to different a priori choices (Table 2). Please note that independent of the spatial distribution of the emissions, the probability distribution of each element of the state vector x_b always follows a Gaussian distribution in all our analytical inversions. The a priori emissions for individual countries were taken from the annual national inventory reports (NIR) to the United Nations Framework Convention on Climate Change (UNFCCC) for the reference year 2018 as reported in April 2020.

For some widely used substances, such as HFC-134a​​​​​​​, we expect that the usage and hence the emissions mostly follow proxies like population and traffic (HFC-134a in mobile air-conditioning). For other compounds, such as SF₆, used as insulator gas in high-voltage installations, the choice of the a priori is not as obvious since their emissions may be more dominated by individual emission hotspots. Nevertheless, a population-based a priori is generally still meaningful (Hu et al., 2023). For these substances, using different a priori fields allows for illustrating how strongly the inverse solution is guided by the a priori and reveals if the higher-resolution transport model inhibits a larger potential to localize emissions independent of the a priori.

Figure 3Different a priori spatial emission distributions – presented on the irregular grid – utilized in the inversions. The population-based a priori can be seen in panel (a), the elevation-dependent a priori in panel (b), and the uniformly distributed emissions per country in panel (c).

The different a priori fields used in this study consist of a population-based a priori, a uniform-per-country a priori, and an elevation-dependent a priori. In the population-based a priori, an emission factor represents the average emissions for each person in the country, and the emissions are given by the emission factor multiplied by the number of residents in each grid cell. In the uniform-per-country case, the emissions are distributed uniformly in the whole country, while in the elevation-dependent a priori, the emissions are distributed uniformly per country below an elevation threshold of 1000 m, whereas above that threshold the emissions were set to 5 % of the low-elevation value. Above the elevation threshold, population densities are usually low in the Alps and very few industrial installations are present, suggesting very limited emissions of the current substances of interest. The spatial distribution of the different a priori emissions can be seen in Fig. 3.

3.1 Preliminary screening tests

Swiss halocarbon emissions using the low-resolution model were estimated after the measurement campaign in Beromünster in 2020 and the results are summarized in Rust et al. (2022). The inversions conducted for their study included observations from two sites in Switzerland (JFJ and BRM) and two sites on the British Isles (MHD and TAC) to constrain the European emissions. Here, we first examine the impact of the additional observational sites (CMN and TOB) on the Swiss national emission estimates. Sensitivity inversions were conducted for both the high- and low-resolution models and for HFC-134a including observations from CMN and TOB (PREL_SITEXT7 and PREL_SITEXT1, Table 2). Only results from the high-resolution model are discussed in the following. In Fig. 4, the spatial distributions of the a posteriori emissions of HFC-134a are displayed for the inversion excluding (PREL_SITRED1) (Fig. 4a) and including CMN and TOB (PREL_SITEXT1) (Fig. 4b), as is the resulting difference (Fig. 4c). For HFC-134a no large spatial differences between the a posteriori emissions of the two inversions can be seen (Fig. 4c). The differences in the total Swiss emissions for the two inversions are also not significant: 308 ± 48 Mg yr⁻¹ for the BASE inversion and 312 ± 50 Mg yr⁻¹ when including CMN and TOB.

Figure 4Spatial distribution of Swiss HFC-134a a posteriori emissions for the inversion PREL_SITRED1 (a) and for PREL_SITEXT1 (b). Panel (c) shows the a posteriori emission differences between panels (a) and (b), while in panel (d) the a posteriori emission differences between panel (a) and BASE1 inversion are shown. In all cases, results from the high-resolution transport model are given.

The same cannot be claimed for observational sites in Switzerland though, since the PREL_SITRED1 inversion changes significantly in terms of both spatial distribution and total emissions when SOT (BASE1) is included (Sect. 3.2, Fig. 4d). This is an expected result since SOT adds sensitivity to regions where BRM is not very sensitive (Fig. 2). Additional observational sites are also essential, since they allow for sampling emissions from the same region at different sites and hence under different atmospheric conditions (advection direction, turbulence regime), thereby improving the representation of dispersion. This happens because turbulent dispersion behaves differently in the near and the far field. In the near field, dispersion approaches isotropy mainly at large scales, meaning that the diffusion in the near field is independent of the size of the eddies. During that phase, the size of the plume is much smaller compared to the larger turbulent eddies, and turbulence acts more like a mean transport mechanism (Csanady, 1973).

Concerning the covariance parameters which were excluded from the maximum likelihood step (Sect. 2.6.3), f_b was initially set to 1, meaning that the baseline is assigned an uncertainty equal to the uncertainty calculated in the REBS method (Sect. 2.3). The latter sometimes leads to unrealistically large adjustments in the baseline and usually underestimation of the emissions, since most of these adjustments tend to increase the baseline considerably. Different sensitivity tests with different values of the factor f_b were conducted to find a representative value for each different receptor and for each different compound (PREL_COV1 and PREL_COV7). Then, for all the inversions for this compound, the value of the factor f_b was fixed and not further optimized in the maximum likelihood step.

Estimating the baseline concentration purely from observations and optimizing it by site may not be the best solution to the baseline problem. Alternatively, baseline observations and transport model information can be used to reconstruct a spatially and temporally resolved baseline concentration at the domain boundaries, from which, again with the transport model information, a baseline concentration for each site and time can be sampled (e.g., Manning et al., 2021; Hu et al., 2023). Common baseline concentrations at the domain boundary, instead of individual baseline concentrations at the sites, may then be included as part of the state vector.

Another factor that is poorly constrained by the maximum likelihood approach is the correlation length, L. Values pointing to overfitting were also obtained from the maximum likelihood, mainly for the low-resolution model. A series of sensitivity runs were deployed for the estimation of a meaningful correlation length, which afterward was used as a fixed value in the inversion for both model resolutions without being further optimized by the maximum likelihood (PREL_COV1 and PREL_COV7). The total emission estimates were not sensitive to small to medium deviations of the correlation length from the chosen value.

Furthermore, we investigated whether we obtain additional information from the high-resolution inversions if we use 3-hourly observation aggregates to drive the inversion instead of 24-hourly aggregates, as used for the low-resolution inversions and in previous studies (Rust et al., 2022) (PREL_AGR1 and PREL_AGR7). There was no significant difference when the 3-hourly aggregates were used, so we maintained the 24-hourly aggregates for the inversions in this study since they result in considerably reduced computational costs.

Finally, we fixed the parameters, which influence the resolution of the inversion grid, to values that lead to similar inversion grids for both FLEXPART model resolutions (PREL_NCEL1 and PREL_NCEL7). However, the total country or total inversion emissions did not show sensitivity to the resolution of the inversion grid (< 1 % differences across the two models) within the range of tested resolutions.

3.2 Emissions of HFC-134a

Currently, HFC-134a is the most often used HFC in Switzerland with reported emissions of 455 Mg yr⁻¹ for 2019 and 415 Mg yr⁻¹ for 2020 (FOEN, 2022). HFC-134a is employed as a refrigerant in mobile air-conditioning (i.e., road traffic) and in stationary refrigeration systems. It is also used as a foam blowing agent. Its 100-year GWP is 1430 and its atmospheric lifetime is approximately 14 years (Engel et al., 2018).

Figure 5Spatial distribution of Swiss HFC-134a a posteriori emissions for the BASE inversion with the 7 km model (a, c, e) and the 1 km model (b, d, f) starting from a population-based a priori (a, b). The same plots (c, d) starting from a spatially uniform a priori and from an elevation-dependent a priori (e, f).

Figure 6A posteriori minus a priori emission differences for HFC-134a for the BASE inversion with the 7 km model (a, c, e) and the 1 km model (b, d, f) starting from a population-based a priori (a, b), a spatially uniform a priori (c, d), and an elevation-based a priori (e, f).

BASE7 inversion leads to an a posteriori estimate of annual Swiss emissions of 260 ± 49 Mg yr⁻¹ (Fig. 9). The a posteriori distribution of the emissions for this inversion can be seen in Fig. 5a and the difference between the a posteriori values and the a priori in Fig. 6a. Compared to the UNFCCC reported bottom-up emissions, there is a significant reduction in HFC-134a emissions almost everywhere in Switzerland except for the region south of SOT and in the canton of Valais, where the emissions are increased compared to the a priori. In Rust et al. (2022) the Swiss emission estimate for HFC-134a – using observations from BRM, JFJ, MHD, TAC, a population-based a priori, and inversions with the low-resolution model (PREL_SITRED1) – was 274 ± 67 Mg yr⁻¹ (2 standard deviations) for 2019–2020.

The same a posteriori emissions distributions but obtained with the high-resolution transport model can be seen in Fig. 5b and the difference from the a priori emissions in Fig. 6b. The total emissions estimate for BASE1 inversion is 351 ± 44 Mg yr⁻¹ (Fig. 9), which is 35 % higher compared to the BASE7 estimate. The BASE1 estimate is closer to the value of the inventory, while BASE7 gives a 40 % lower estimate than the inventory. The relative uncertainty in the estimate is higher for BASE7 at 18.8 % compared to the BASE1 inversion at 12.5 %. This is an indication of improved use of the information content of the observations by the high-resolution model due to the improved representation of the atmospheric flow.

If we consider the difference between the a posteriori emissions of the high- and low-resolution model inversions (BASE1 and BASE7), the BASE1 inversion enhances the emissions in all big cities of Switzerland, in the regions with the most industrial activity (canton of Aargau, west-northwest of Zurich), and along the traffic network (Fig. 8). The arc with increased emissions from Zurich to Bern in Fig. 8 could point to emissions from the main highway of Switzerland (A1) forming the main west–east transport route and connecting two of the biggest cities of Switzerland. To evaluate the connection between HFC-134a emissions and traffic, we calculated the correlation between a posteriori emissions and CO₂ traffic emissions, as taken from the spatially resolved Swiss emission inventory (Heldstab et al., 2021), for the two different inversions. The a posteriori emissions from the BASE1 inversion show a higher correlation of r=0.6 with the traffic CO₂ emissions compared to the BASE7 inversion emissions of r=0.35. Since we start from a population-based a priori, the correlation between the a posteriori emissions and the population was additionally estimated. A posteriori emissions from the BASE1 inversion possess a correlation of r=0.96 with the population, while the emissions from the BASE7 inversion are slightly less correlated at r=0.8. Hence, the high-resolution model inversion stays closer to the a priori distribution compared to the low-resolution model.

Table 3Statistical measures used to assess the reliability of inversion for different compounds, different transport model resolutions, and different a priori emissions. The table displays the reduced χ² index, degrees of freedom (DOF), root mean squared error (RMSE), and correlations of simulated compound values against observations for Beromünster (BRM) and Sottens (SOT), as well as a posteriori emissions, E_a, for Switzerland.

Download Print Version | Download XLSX

Figure 7 shows the mole fraction time series of daily averaged HFC-134a for both low- and high-resolution BASE inversions and the observations at BRM (Fig. 7a) and SOT (Fig. 7b). Both inversions represent the observed variability closely; however, some distinct differences do exist. To assess the performance of our inversions, we used the following statistical measures: reduced χ² index, which is a measure of the normalized variance between the observations and the simulated values; the degrees of freedom, which is a measure of the relative uncertainty reduction between the a priori and the a posteriori; the correlation coefficient, r; and the root mean square error (RMSE) of the simulated versus the observed mole fractions (Table 3). From these statistics for HFC-134a, we conclude that both BASE1 and BASE7 inversions are reliable, but the BASE1 inversions show improved performance at the receptors in Switzerland. A general performance improvement of FLEXPART-COSMO-1 versus FLEXPART-COSMO-7 can already be seen in the a priori simulations (see Table S1​​​​​​​ in the Supplement), where it is solely due to improvements in the transport description and flow in complex terrain, as emission distribution and baseline values were the same for both sets of simulations.

Figure 7Time series of observed (blue lines) and simulated (red and green lines) HFC-134a mole fractions at BRM (a) and SOT (b). A posteriori simulations for the BASE inversion and with the low-resolution (C7, green lines) and high-resolution (C1, red lines) transport model are given.

Download

Figure 8A posteriori emission differences between the high- and low-resolution model inversions with population-based a priori for HFC-134a.

Figure 9A posteriori emissions for all the substances utilized in this study for the different model resolutions and the different a priori choices. All results correspond to the BASE inversions.

Download

For the inversions with uniformly distributed a priori emissions by country (BASE_ED; Figs. 5c, d–6c, d), the results are inferior to the results obtained with the population-based a priori. The BASE_ED inversions tend to retain and cannot completely remove the emissions from the Alpine region, while the distribution of emissions in the Swiss Plateau does not look very plausible, at least for the low-resolution model. The inversion with the high-resolution model seems to be able to locate the emission hotspots north of BRM and northwest of Zurich. These results highlight the importance of the a priori spatial distribution. If the inversion is initialized with a highly unrealistic a priori and there is an insufficient observational constraint, the inversion may not converge to a realistic state.

In Figs. 5 and 6, panels (e) and (f) correspond to inversions with an elevation-dependent a priori (BASE_ED). The total Swiss emission estimate is 318 ± 62 Mg yr⁻¹ for the BASE_ED1 and 217 ± 46 Mg yr⁻¹ for the BASE_ED7 (Fig. 9). For the inversions with the high-resolution model (Figs. 5f and 6f), we can see that the inversion converges again towards a population-based distribution, especially close to the observational sites. The hotspots of emissions in the cantons of Zurich and Aargau are reconstructed by the inversion, although not as sharply as for a population-based a priori, along with the hotspots in the Lausanne and Geneva regions. However, the inversion using the low-resolution transport model cannot recover the population-based distribution to the same degree (Fig. 5c). The emissions from Zurich in particular seem to be allocated too far to the west at a closer distance to BRM. Similarly, emissions from Geneva are not indicated as prominently. These observations are corroborated by the correlation between the a posteriori emissions and population, which was r=0.56 for the high-resolution model and only r=0.31 for the low-resolution model. Since we are highly confident that the emissions of this substance should be correlated with population density, the BASE_ED inversions show that the high-resolution model inversions are much more accurate in reconstructing the true distribution.

Moreover, the HFC-134a inversions with seasonally variable emissions (SEAS) reveal the existence of a seasonality pattern in the emissions in Switzerland. In both types of seasonal inversions with the high-resolution model (see Sect. 2.6.6) there is a clear annual variability of HFC-134a with the peak during the summer months of June–August (JJA) – 433 ± 94 Mg yr⁻¹ for SEAS1 and 364 ± 92 Mg yr⁻¹ for the SEAS2 inversion – and the minimum during the winter – 238 ± 98 Mg yr⁻¹ for the SEAS1 239 ± 100 Mg yr⁻¹ for the SEAS2 inversion. This corresponds to a seasonal amplitude of approximately 1.3, which is similar to seasonal amplitudes obtained by Hu et al. (2017) for HFC emissions in North America. Given the relatively large a posteriori uncertainties in seasonal emissions, summer and winter emissions are significantly different at the 95 % confidence level for the SEAS1 inversions, but not for the SEAS2 inversions. The spatial distribution of the emissions for the different seasons is similar, pointing to the conclusion that the emissions in all different seasons have the same sources, but the leakage of HFC-134a from refrigeration systems is higher when they are in use. According to the statistical measures used, the SEAS2 inversion is superior to the SEAS1, possessing a better correlation of the simulated values when compared to the observations and reduced χ² much closer to 1. The total annual emission estimate for the two inversions is 320 ± 50 Mg yr⁻¹ for SEAS2 and 342 ± 48 Mg yr⁻¹ for the SEAS1 inversion, close to the BASE estimate. Hence, there seems to be no big gain when we consider inversions with seasonality when the main target is the validation of annual total emissions. However, these simulations could help improve our understanding of the release mechanisms of these compounds. Additionally, the a posteriori model performance slightly increased for all seasonal inversions compared to the annual mean population-based inversions (see Table S2). For the RMSE, the performance increase was of the order of 10 % but was less pronounced for the correlation coefficient. SEAS2 inversions achieved a larger degree of freedom, but also revealed a χ² index considerably larger than 1, indicating some degree of overfitting. For the low-resolution model, a much smaller (insignificant) seasonal amplitude was obtained in the a posteriori emissions. Whether this is due to the reduced ability of the model to realistically reproduce the diurnal mole fraction variability compared to the high-resolution model (Katharopoulos et al., 2022) or to potential seasonal transport biases will need to be investigated in future studies.

Based on the analysis in this section, we can claim that the high-resolution inversions reconstruct the spatial distribution of the HFC-134a emissions in Switzerland better and with more detail than the low-resolution inversions. The total Swiss emission estimates between the two models differ significantly, with the high-resolution model predicting values closer to those in the inventory. For the remaining halocarbon inversions in this work, we present only the results from population-based a priori and elevation-dependent a priori since the elevation-dependent a priori can be seen as an improved version of the uniform-by-country a priori. Figures for the simulations with a uniform spatial a priori distribution can be found in the Supplement.

3.3 Emissions of HFC-125

HFC-125 is the second most abundant HFC in Switzerland, with reported emissions of 122 Mg yr⁻¹ for 2019 and 2020 (FOEN, 2022). HFC-125 is employed as a refrigerant mainly in stationary refrigeration systems, and as a result, its emissions are expected to be from static sources. It is also used as a fire suppression agent in fire extinguishers, but this use is forbidden in Switzerland. Although, on a mass basis, HFC-125 emissions are lower compared to HFC-134a, their impact as a GHG is higher since the 100-year GWP is 3500, almost 3 times that of HFC-134a. The setup used to estimate the Swiss emissions of HFC-125 is identical to the setup used for HFC-134a including the two Swiss sites (BRM, SOT).

Figure 10Spatial distribution of a posteriori emission for HFC-125 (a–d), SF₆ (e–h), and HFC-32 (i–l) for the high- and low-resolution inversions as well as the population-based (rows 1–2) and elevation-dependent a priori (rows 3–4).

Figure 11Spatial distribution of a posteriori minus a priori emission difference for HFC-125 (a–d), HFC-32 (e–h), and SF₆ (i–l) for the high- and low-resolution inversions as well as the population-based (rows 1–2) and elevation-dependent a priori (rows 3–4).

The BASE7 inversion yields Swiss a posteriori emissions of 78 ± 20 Mg yr⁻¹ (Figs. 9, 10, and 11)​​​​​​​. In Rust et al. (2022) the Swiss emission estimate for HFC-125 with PREL_SITRED7 was 107 ± 28 Mg yr⁻¹ (2 standard deviations) for 2019–2020 (Figs. S1–S2). The estimate from a second top-down method used in Rust et al. (2022) to estimate the emissions from BRM (tracer-ratio method) was 94 ± 19 Mg yr⁻¹. The addition of SOT yields approximately 20 % lower annual Swiss emissions estimates. The a posteriori minus a priori emission difference figures for the two cases (not shown) depict the fact that the PREL_SITRED7 inversion increases the emissions of HFC-125 compared to the a priori north and northwest of BRM, in the Valais region, and in the regions north and south of SOT, whereas the BASE7 inversion including both BRM and SOT increases the emissions compared to the a priori only in a small radius around BRM and in the canton of Valais. In all other regions of the Swiss Plateau, a significant decrease in emissions is observed.

For the high-resolution BASE inversion, BASE1, the a posteriori emissions can be seen in Fig. 10b and the difference from the a priori emissions in Fig. 11b. The total emissions estimate in this case is 101 ± 21 Mg yr⁻¹. This number is 22 % higher compared to the low-resolution model estimate. Similar to HFC-134a, the a posteriori uncertainty is higher (although only slightly) in BASE7 at 26 % compared to the BASE1 inversion at 21 %. The BASE1 estimate is closer to the value of the inventory, whereas the inversion estimate with the BASE7 corresponds to about $\mathrm{2}/\mathrm{3}$ of the inventory value. The BASE1 inversion increases the emissions in the region to the north of BRM, in Basel, and in eastern Switzerland close to the borders with Austria and Germany (Fig. 11b). In contrast to the BASE7 inversion, the BASE1 inversion increases the emissions for all the big cities of Switzerland and in the industrial region ranging from Zurich to Basel (Fig. S3).

Additionally, in Figs. 10c–d and 11c–d, the a posteriori emissions and the differences from the a priori for HFC-125 can be seen starting from an elevation-dependent a priori (BASE_ED). As with HFC-134a, the inversions converge again towards a population-based a priori that is especially close to the observational sites. The hotspots of emissions in the cantons of Zurich and Aargau are reconstructed by the inversion, along with the hotspots in the Lausanne and Geneva regions. The BASE_ED1 again tends to produce an a posteriori distribution closer to the population distribution compared to the BASE_ED7 inversion. The total Swiss emission estimate for the BASE_ED1 inversion with the elevation-dependent a priori is 90 ± 18 Mg yr⁻¹, while that for the BASE_ED7 inversion is 83 ± 20 Mg yr⁻¹. The results for the inversions starting from uniformly distributed emissions (BASE_UNI) can be seen in the Supplement (Figs. S14–S17).

Both the high- and low-resolution inversions are reliable (Table 3), since they present reasonable reduced χ² and they both lower the uncertainty of the a priori emissions (DOF). The high-resolution inversion for HFC-125 possesses a slightly higher correlation and slightly smaller RMSE of simulated versus observed values at the Swiss receptors (Table 3). Since HFC-125 is used in stationary air-conditioning systems, its usage should be concentrated in the big cities and in industrial areas and should partially follow a population-based distribution. Hence, the increase in emissions in the areas west of Zurich, north of BRM, and south of Basel looks reasonable.

3.4 Emissions of HFC-32

Difluoromethane or HFC-32 is the fourth most emitted HFC in Switzerland with reported emissions of 57 Mg for 2020 and an increasing emission trend. HFC-32 is employed as a refrigerant for the same purposes as HFC-134a. Hence, we expect the spatial distribution of its emissions to be similar to HFC-134a. Its lifetime is only 5 years and its GWP is correspondingly low (705). Thus, it is a relatively low-risk choice among HFC refrigerants.

The BASE7 and BASE1 inversions for the period 2019–2021 yield a posteriori annual emissions of 38 ± 8 and 50 ± 8 Mg yr⁻¹, respectively (Figs. 9, 10, and 11)​​​​​​​. In Rust et al. (2022) the Swiss emissions estimate with PREL_SITRED7 inversion with population-based a priori emissions for HFC-32 was 44 ± 12 Mg for 2019–2020 (Figs. S4–S5).

Comparing panel (e) with panel (f) in Fig. 11 again reveals significant differences between the two model versions. While both increase the emissions in the canton of Valais, the Lausanne region, and around SOT, there is a significant difference for the rest of the Swiss Plateau, where the BASE7 inversion decreases the emissions, whereas the BASE1 inversion mostly increases the emissions. In Fig. S6 the a posteriori emission differences between the high- and low-resolution inversions are shown. These results are very similar to the results for HFC-125. The latter, together with the resemblance of the a posteriori emissions between HFC-32, HFC-125, and HFC-134a for all model resolutions, verifies our prior assumption that the emissions for these substances have similar sources.

In panels (g) and (h) in Figs. 10 and 11 the a posteriori emissions and a posteriori minus a priori emission differences for HFC-32 are depicted starting from an elevation-dependent a priori (BASE_ED). Again the results are very similar to those of HFC-134a, and the a posteriori emissions again reconstruct a population-based distribution. The hotspots of emissions in the cantons of Zurich and Aargau are reconstructed by the inversion, along with the hotspots in the Lausanne and Geneva regions. The total Swiss emission estimate for the BASE_ED1 inversion with the elevation-dependent a priori is 47 ± 5 Mg yr⁻¹, while for the BASE_ED7 inversion it is 35 ± 4 Mg yr⁻¹. The results for the inversions starting from a uniform distribution by country (BASE_UNI) can be seen in the Supplement (Figs. S18–S21). The statistical measures assessing the reliability and performance of the results are summarized in Table 3, confirming the generally improved performance of the high-resolution model at the receptor sites for all a priori distributions.

3.5 Emissions of SF₆

Sulfur hexafluoride or SF₆ is also an F-gas that is mainly used (80 % of its emissions, Simmonds et al., 2020) in the electrical power industry as a gaseous dielectric medium. Another utilization of SF₆ is in semiconductor manufacturing and as inert gas for the casting of magnesium. The Swiss inventory value for SF₆ emissions was 6.7 and 6.0 Mg yr⁻¹ for 2019 and 2020, respectively. Although its emissions are lower compared to those of the HFCs presented before, its very large GWP (23 500) makes SF₆ the fourth most important contributor of the F-gases to anthropogenic warming. Since SF₆ is used as an electrical insulator by the electrical industry, its emissions may be concentrated on point sources, and the choice of a population-based a priori may not be as obvious as with the HFCs.

The BASE7 and BASE1 inversions for the period 2019–2021 yield a posteriori annual emissions of 7.6 ± 1.1 and 9.0 ± 1.1 Mg yr⁻¹, respectively (Figs. 9, 10, and 11). Sensitivity tests using both uniform and population-based a priori (Rust et al., 2022) showed that the population-based a priori is still the best option for this compound. The Swiss emissions estimate for SF₆ – using only observations from BRM (PREL_SITRED7) and a population-based a priori – was 9 ± 2 Mg yr⁻¹ for 2019–2020 (Figs. S7–S9). These results are consistently larger than the inventory values. A potential cause of this discrepancy could be the strong emission sources in southwestern Germany (as indicated in our a posteriori results; also compare Simmonds et al., 2020), which may have potentially been incorrectly attributed to Switzerland as well.

Comparing panels (i) with (j) in Fig. 11 the distribution of the emissions for both model resolutions reveal very similar patterns of increases and decreases in emissions compared to the a priori. The difference of 15 % in total emissions between the two inversions comes from the areas of Lausanne, Geneva, and north and northwest of BRM. The BASE1 inversion increases the emissions more in these regions than the BASE7 inversion (Fig. S9).

In panels (k) and (l) in Figs. 10 and 11 the a posteriori emissions and a posteriori minus a priori emission differences for SF₆ are summarized starting from an elevation-dependent a priori (BASE_ED). The BASE_ED1 inversion does a better job of reconstructing the hotspots of emissions north of BRM and around Lake Geneva. The total Swiss emission estimate for the BASE_ED1 inversion is 7 ± 1.2 Mg yr⁻¹, while for the BASE_ED7 inversion it is 5 ± 1 Mg yr⁻¹. The results for the BASE_UNI inversions can be seen in the Supplement (Figs. S22–S25), whereas the statistical measures assessing the reliability and performance of the results are summarized in Table 3. In contrast to the other compounds, the spatial distribution of the emissions for SF₆ – when a uniform a priori is used – seems more reasonable and highlights the same emission regions as the inversions with the population- and elevation-dependent a priori.