Journal topic
Atmos. Chem. Phys., 20, 3987–3998, 2020
https://doi.org/10.5194/acp-20-3987-2020
Atmos. Chem. Phys., 20, 3987–3998, 2020
https://doi.org/10.5194/acp-20-3987-2020

Research article 02 Apr 2020

Research article | 02 Apr 2020

# Using ship-borne observations of methane isotopic ratio in the Arctic Ocean to understand methane sources in the Arctic

Using ship-borne observations of methane isotopic ratio in the Arctic Ocean to understand methane sources in the Arctic
Antoine Berchet1, Isabelle Pison1, Patrick M. Crill2, Brett Thornton2, Philippe Bousquet1, Thibaud Thonat1, Thomas Hocking1, Joël Thanwerdas1, Jean-Daniel Paris1, and Marielle Saunois1 Antoine Berchet et al.
• 1Laboratoire des Sciences du Climat et de l'Environnement, CEA-CNRS-UVSQ, IPSL, Gif-sur-Yvette, France
• 2Department of Geological Sciences, Stockholm University, 10691 Stockholm, Sweden

Correspondence: Antoine Berchet (antoine.berchet@lsce.ipsl.fr)

Abstract

Characterizing methane sources in the Arctic remains challenging due to the remoteness, heterogeneity and variety of such emissions. In situ campaigns provide valuable datasets to reduce these uncertainties. Here we analyse data from the summer 2014 SWERUS-C3 campaign in the eastern Arctic Ocean, off the shore of Siberia and Alaska. Total concentrations of methane, as well as relative concentrations of 12CH4 and 13CH4, were measured continuously during this campaign for 35 d in July and August. Using a chemistry-transport model, we link observed concentrations and isotopic ratios to regional emissions and hemispheric transport structures. A simple inversion system helped constrain source signatures from wetlands in Siberia and Alaska, and oceanic sources, as well as the isotopic composition of lower-stratosphere air masses. The variation in the signature of lower-stratosphere air masses, due to strongly fractionating chemical reactions in the stratosphere, was suggested to explain a large share of the observed variability in isotopic ratios. These results point towards necessary efforts to better simulate large-scale transport and chemistry patterns to make relevant use of isotopic data in remote areas. It is also found that constant and homogeneous source signatures for each type of emission in a given region (mostly wetlands and oil and gas industry in our case at high latitudes) are not compatible with the strong synoptic isotopic signal observed in the Arctic. A regional gradient in source signatures is highlighted between Siberian and Alaskan wetlands, the latter having lighter signatures (more depleted in 13C). Finally, our results suggest that marine emissions of methane from Arctic continental-shelf sources are dominated by thermogenic-origin methane, with a secondary biogenic source as well.

1 Introduction

Methane (CH4) is both a potent greenhouse gas and a precursor of ozone with very diverse sources and sinks in the atmosphere . The wide variety of CH4 sources and their spatial and temporal heterogeneity make the uncertainties in CH4 budgets very large on both regional and global scales . This impairs our understanding of the variations in atmospheric concentrations, particularly of which sources of methane and/or regions are causing these variations, which have been rapid in recent decades .

In the Arctic, major CH4 sources are natural wetlands, inland waters (lakes, streams, deltas, estuaries), leaks from oil and gas extraction and transport, wildfires, seabeds, and geological seepage. The magnitude of all these sources suffers with very high uncertainties . The large areas of wetlands above 50 N, and the high sensitivity of their CH4 emissions to the changing climate, make this zone a key region for the global CH4 budget. The present uncertainties in CH4 sources and sinks in the Arctic are very large due to the complexity of the involved processes and the difficult access to these remote regions (e.g. Thornton et al.2016b; Bohn et al.2015). Moreover, in addition to increased CH4 emissions from wetlands and thawing permafrost, increasing ocean temperatures could lead to the destabilization of methane hydrates on the Arctic continental shelf, potentially emitting large quantities of CH4. For instance, significant point emissions have been detected along the East Siberian Arctic Shelf , taking the shape of CH4 flaring from the seafloor and extending up to the surface. However, upscaling point measurements of “hotspots” proves difficult, and there is no proof that such methane hydrate emissions currently reach the atmosphere in large quantities . Other potential Arctic seafloor sources of CH4 include emissions from degrading subsea permafrost , leakage from natural gas reservoirs and degrading terrestrial organic carbon transported onto the continental shelf . CH4 emissions from the Arctic would then have a positive feedback on climate change. Better knowledge of Arctic CH4 emissions would reduce uncertainties in its global budget and help to better quantify the sensitivity of Arctic regional sources and sinks to climate change.

For more than 10 years, atmospheric measurements of methane concentrations have been performed in the Arctic at surface stations , during mobile field campaigns such as the YAK-AEROSIB (Airborne Extensive Regional Observations in SIBeria) aircraft campaigns and the TROICA train campaign (Tarasova et al., 2006, 2009), or during oceanographic campaigns . In the present work, we analyse data from the SWERUS-C3 campaign aboard a ship in the Arctic Ocean during summer 2014 . Such short-term mobile campaigns are necessary for complementing the limited number of long-term fixed, mostly coastal stations currently available. In particular, oceanic campaigns are expected to provide information not only on oceanic sources but also on land sources located upwind. However, CH4 from various sources is being mixed during the atmospheric transport of the air masses, which makes it difficult to separate them without resorting to numerical modelling .

Atmospheric inversions merge together observations, numerical modelling and emission datasets to attribute the observed variability in CH4 concentrations to emitting regions and thus optimize the CH4 budget. Such methods were successfully applied in the Arctic using in situ fixed stations as well as satellites when available . But despite technical progress in numerical modelling and inversion methods, it is hardly feasible to separate co-located emissions from different emitting sectors upwind observation sites based on observations of CH4 concentrations alone. Observations of methane isotopic ratios could help in separating emission sectors, as the main emission processes are isotopically fractionating, causing significantly different isotopic source signatures. For example, high-latitude wetlands were attributed signatures in a range of −80 ‰ to −55 ‰ . The δ13CCH4 signature of atmospheric CH4 above the Arctic Ocean has been previously reported in the range of −50 ‰ to −47 ‰ . Isotopes have already been used to characterize the origin of air masses in the Arctic , though these studies concluded that refinements in qualifying source emission isotopic signatures are required.

In the following, we explore the potential of using observations of isotopic ratios in the Arctic Ocean together with total CH4 concentrations to separate pan-Arctic emission sources. We further analyse emission isotopic signatures in the Arctic from integrated atmospheric observations. We base our analysis on the unique observation set collected during the ship-based campaign SWERUS-C3 during summer 2014 in the Arctic Ocean. By comparing measurements to simulations of total CH4 and the isotopic ratio, we analyse the extent to which the observable signal in the Arctic Ocean is exploitable in a numerical inversion system. In Sect. 2, we explain our inversion approach alongside giving details on the SWERUS-C3 observation campaign and on the model CHIMERE used in our study. In Sect. 3, we compare observations to simulations to assess the main contributions to the signal variability and then implement a simplified inversion system to quantify isotopic emission signatures from various emission sectors around the Arctic.

2 Methods

Figure 1Path of the icebreaker Oden during the SWERUS-C3 campaign and domain of simulations. (a) The ship positions are represented by grey and brown lines, with brown parts corresponding to locations where isotopic observations where carried out. The area delimited by coloured lines is the domain of CHIMERE simulations used for this study (see Sect. 2.2). The shaded areas and associated numbers correspond to the regions and their IDs used to separate contributions from remote emissions to the observed signal, as detailed in Sect. 2.3. ESAS is the East Siberian Arctic Shelf. CHIMERE boundary conditions are split along the four sides of the domain as indicated by the coloured lines. (b) Zoomed-in view of the area covered by the campaign. The icebreaker's locations are coloured based on their corresponding dates. Ship positions with a black edge are locations where isotopic observations where carried out. More details on the campaign in . The shaded area corresponds to the ESAS emission region used in our simulation set-up.

## 2.1 Campaign and instrument description

Observations were carried out during the SWERUS-C3 campaign aboard the Swedish icebreaker Oden between 14 July and 26 September 2014. The cruise path went through the central and outer Laptev and East Siberian seas and finally the Chukchi Sea to Point Barrow, Alaska, in a first leg (see Fig. 1). A second leg of the cruise headed north from Point Barrow back through the Chukchi Sea and into the Arctic Ocean. As shown in Fig. S1 in the Supplement, sea ice cover was present during a large portion of the campaign. Regions known to have active seafloor gas seeps (see Thornton et al.2016a) occurred in both ice-free (in the Laptev Sea) and ice-covered (in the East Siberian Sea) regions.

Concentrations of total CH4 were measured during the whole campaign using an off-axis cavity ring-down laser spectrometer from Los Gatos Research (LGR), Inc. (Model 0010, FGGA-24EP, Mountain View, California, USA). Air inlets were located at 9, 15, 20 and 35 m above the sea surface; air was pulled through all inlets continuously and analysed from one inlet at a time for 2 min before switching to the next inlet. Data were filtered using wind speed and direction to avoid contamination from the ship exhaust. As no local sources influenced our measurements, concentrations are similar at all levels. We concatenate measurements from all inlets indifferently for our study. The spectrometer was calibrated every 2 h using two synthetic air target gases; the target gases themselves were calibrated before, during and after the cruise to two NOAA Earth System Research Laboratory-certified standards for CH4. The reported precision was 0.5 ppb. Further details on the campaign conditions and instrument configuration are available in .

Isotopic ratios were measured only during the first leg of the campaign, from 14 July to 26 August (see Fig. 1), using an Aerodyne Research, Inc. (Billerica, MA, USA), direct absorption interband cascade laser spectrometer. This spectrometer measured the concentrations of the CH4 isotopologues 12CH4, 13CH4 and CH3D, the last of which is not discussed in the current paper. The more common isotope ratio mass spectrometry methods directly provide (as their name implies) an isotope ratio. In contrast, because the Aerodyne spectrometer measures the individual isotopologues, they must be individually calibrated before converting to 13CCH4 values; this method is described in .

## 2.2 Model description

The Eulerian model CHIMERE was run to simulate total concentrations of CH4 as well as partial 12CH4 and 13CH4 concentrations to compute CH4 isotopic ratios afterwards using the following formula:

$\begin{array}{}\text{(1)}& \mathit{\delta }{}^{\mathrm{13}}\mathrm{C}=\frac{{\left(\frac{\left[{}^{\mathrm{13}}\mathrm{C}\right]}{\left[{}^{\mathrm{12}}\mathrm{C}\right]}\right)}_{\text{sim}}}{{\left(\frac{\left[{}^{\mathrm{13}}\mathrm{C}\right]}{\left[{}^{\mathrm{12}}\mathrm{C}\right]}\right)}_{\text{ref}}}-\mathrm{1},\end{array}$

with ${\left(\frac{\left[{}^{\mathrm{13}}\mathrm{C}\right]}{\left[{}^{\mathrm{12}}\mathrm{C}\right]}\right)}_{\text{ref}}=\mathrm{0.0112372}$ being the reference ratio from .

The domain of simulations spans over most of the Northern Hemisphere, with a horizontal resolution of ∼100km in order to include most contributions from distant sources (see Fig. 1). Similarly, the model uses 34 vertical levels from the surface up to 150 hPa to represent stratosphere-to-troposphere intrusions. A spin-up period of 6 months prior to the campaign was used to properly assess the impact of air masses transported for long periods before reaching the Arctic Ocean. The chemical sink of CH4 by OH radicals is explicitly computed in CHIMERE using pre-computed fixed OH fields from the chemical model LMDz-INCA (Interaction with Chemistry and Aerosols; Hauglustaine et al.2004; Folberth et al.2006).

CHIMERE runs use the following input data streams: (i) meteorological fields downloaded from the European Centre for Medium Range Weather Forecasts (https://www.ecmwf.int) at 0.5 resolution every 3 h, (ii) anthropogenic emissions aggregated at the CHIMERE resolution from the EDGARv4.3.2 database at 0.1 horizontal resolution , (iii) wetland emissions interpolated from the model ORCHIDEE at 0.5 horizontal resolution , (iv) boundary CH4 concentration fields extracted from the general circulation model LMDz (these global simulations include both the chemical sinks of OH and chlorine as well as their impact on the isotopic ratios, and Cl and OH fields are prescribed offline from the chemical model LMDz-INCA) and (v) and isotopic signatures of the different sources chosen from .

The chemical sink by chlorine is not included in our set-up to keep simulations as light as possible. This sink can be separated into two main contributions: the upper stratosphere and the Arctic Ocean boundary layer. The upper stratosphere is not included in our model of simulation, but chlorine sink (and isotope fractionation) is explicitly accounted for in global LMDz simulations used as boundary conditions in our set-up. Regarding the Arctic Ocean boundary layer, the set-up by was adapted to our case, including the boundary layer Cl sink using pre-computed fields from the model LMDz-INCA. It resulted in differences of concentrations lower than 1 ppb over the Arctic Ocean and less than 0.02 ‰ for the isotopic ratio of air masses, which is negligible compared to the signal we are inquiring into.

Other fluxes not included in our set-up play a significant role in the regional pan-Arctic budget, such as inland water bodies, wildfires and the sink in soil, but have limited impact on our observations. These fluxes were tested in our case and were quantified to cause differences in simulated concentrations lower than 2 ppb and less than 0.01 ‰ in simulated isotopic ratios at the locations sampled during the SWERUS-C3 campaign.

## 2.3 Atmospheric inversion of isotopic signature

Usually observations of δ13CCH4 are used to help with constraining methane fluxes and differentiating between different sources with known signatures. However, the intrinsic spatial and temporal variability in source isotopic signatures limits the robustness of this approach (e.g. Fisher et al.2017, as illustrated in Sect. 3.1). Here, we conversely assume that total CH4 is properly simulated by our model (as confirmed by the good performance of the model to reproduce total CH4 concentrations, highlighted in Sect. 3.1) and that the relative contributions of various sources from various regions are correct. Thus we use δ13CCH4 observations to help reduce uncertainties in source isotopic signatures: we test the ability of the ship-based measurements to help constrain the isotopic signature of remote sources, such as wetland sources and oceanic emissions from the Laptev, East Siberian and Chukchi seas, dominant in the region explored during the campaign.

To do so, δ13CCH4 observations are implemented into a classical analytical Bayesian framework . The designed inversion system optimizes source signatures from different source types and different regions. At every time step when an isotopic observation is available, the system fits observations of isotopic ratios by altering the isotopic ratio in air masses coming from relevant source types and regions. Thus, the control vector contains one isotopic ratio value to optimize for each time step, each sector and each region, as detailed in Eq. (3) below.

The isotopic ratios of wetlands, solid fossil fuels, oil and gas, other anthropogenic sources from various land regions, and a potential variety of marine sources (gas field leaks, decomposing hydrates, degrading permafrost, etc.) from the East Siberian Arctic Shelf (ESAS) as well as from air masses coming from the sides and roof of our domain of simulations are optimized in the system. Apart from the ESAS, emissions are spatially differentiated into 23 geographical regions (see Fig. 1). Contributions from different regions and sectors are differentiated by computing so-called response functions by region, emission type and boundary side. That is to say, we carry out individual CHIMERE chemistry-transport simulations for every region, every type of emission and every side of the domain, with all the other emissions and boundary conditions being switched off, resulting in an ensemble of 98 response functions (23 regions × 4 sectors + ESAS + 4 sides + top).

The simulated isotopic final composition y(t) at every given time step t when an observation is available is retrieved by scaling relative contributions according to assumed source signatures (or original average composition for boundary conditions) as follows:

$\begin{array}{}\text{(2)}& \mathbit{y}\left(t\right)=\sum _{r\in \text{regions}}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}\sum _{s\in \text{sectors}}{\mathit{\alpha }}_{r,s}\left(t\right)×{\mathit{\delta }}_{r,s}\left(t\right),\end{array}$

with r and s varying over all available regions and sectors respectively; ${\mathit{\alpha }}_{r,s}^{t}$ ($\mathrm{0}<{\mathit{\alpha }}_{r,s}^{t}<\mathrm{1}$) is the relative contribution of the sector s from region r at time t, and ${\mathit{\delta }}_{r,s}^{t}$ is the signature (in ‰) of the sector s from region r at time t.

This linear relationship allows us to define the control vector x and the observation operator, linking the control vector to observations of isotopic ratios, to easily compute and scale the simulated isotopic composition:

Given the prior control vector xb containing assumed source signatures before inversion, the observation vector yo and the observation operator H, optimized signatures are obtained by solving the Bayesian problem equation:

$\begin{array}{}\text{(4)}& {\mathbit{x}}^{\mathrm{a}}={\mathbit{x}}^{\mathrm{b}}+\mathbf{K}\left({\mathbit{y}}^{\mathrm{o}}-\mathbf{H}{\mathbit{x}}^{\mathrm{b}}\right),\end{array}$

with $\mathbf{K}={\mathbf{P}}^{\mathrm{b}}{\mathbf{H}}^{\mathrm{T}}\left(\mathbf{R}+{\mathbf{HP}}^{\mathrm{b}}{\mathbf{H}}^{\mathrm{T}}{\right)}^{-\mathrm{1}}$ being the Kalman matrix.

The matrix R represents uncertainties in the observations and in the capability of the model to reproduce them. In our case, we set them uniformly to 1.5 ‰ (1 ‰ from observation errors and 0.5 ‰ from simulation errors). The matrix Pb represents uncertainties and covariances in the prior knowledge we have on source signatures. We build the matrix Pb following the values in Table 1, deduced from and . Ranges and prior signatures for boundary conditions are deduced from global simulations with the model LMDz. Observation time steps are not optimized separately. Instead, we use temporal correlations in the Pb between different time steps. We represent temporal correlations between two time steps ti and tj as

$\begin{array}{}\text{(5)}& r=\mathrm{exp}\left(-\frac{|{t}_{i}-{t}_{j}|}{\mathit{\tau }}\right),\end{array}$

with τ being the temporal correlation scale of Table 1.

As shown in Table 1, the values of source signatures are not well known, and a very large range of signatures are available in the literature. To account for this large variety of realistic signatures, we carry out a Monte Carlo ensemble of 8000 inversions with varying prior signatures and uncertainties instead of running one single inversion. Prior signatures are sampled following a normal distribution with average and standard deviation from Table 1; the standard deviation is chosen as half of the min–max range. Uncertainties are sampled following a uniform distribution spanning over $\left[{\mathit{\sigma }}_{\text{ref}}/\mathrm{2},{\mathit{\sigma }}_{\text{ref}}\right]$, with σref equaling half of the min–max range of Table 1.

In the end, we obtain hourly posterior signatures for each simulated sector and region for each of the 8000 inversions. Even though posterior signatures are available for each region and each sector at each observation time step, we do not inquire into the temporal variability in sources, as constraints provided by the SWERUS observations are very heterogeneous in time and space. Instead, we compute overall posterior distributions for each simulated sector and region based on an ensemble of 8 000 000 (8000 inversions × 1000 hourly observations). To minimize the impact of control vector components that are ill-constrained by the inversion, all data points are not evenly counted in posterior distributions. Posterior distributions of signatures are computed, accounting for all the Monte Carlo samples and weighted by the corresponding values of the sensitivity matrix KH , which gives an indicator of how much observations constrain one component of the control vector. The posterior optimal signature for each region and sector is computed as the maximum of the probability distribution.

Table 1Isotopic signatures for the inputs in the CHIMERE model. The min–max range is deduced from existing literature . The prior signature is computed as the centre of the min–max range.

3 Results and discussion

## 3.1 Forward modelling of total methane and isotopic ratio

Figure 2 shows observations of total CH4 and of isotopic ratios as measured during the campaign and compared to simulations. The model CHIMERE reproduces most of the variability in the total CH4 signal well. The average bias over the period is lower than 5 ppb, with a correlation of 0.66 between observations and simulations on an hourly basis. Most peaks spanning more than 1 d are properly represented in the model, proving the capability of the model to reproduce the synoptic variability in the observations. Smaller peaks are missed by the model, in particular on 5, 12 and 15 August, indicating that some local sources are not included in the model or are dispersed too quickly in the numerical realm. These could be local intense seeps that met along the ship's track or onshore wetlands not well represented with the model ORCHIDEE at 0.5 horizontal resolution. We do not investigate further missing emissions, as most peaks are well explained by the model, which we assume is sufficient to carry out an inversion of isotopic signatures as described in Sect. 3.2.

When computing the intersect with the y axis of the linear fit between δ13CCH4 and total CH4 (see Keeling plots in the Supplement), the observed isotope ratios point to an average generic Arctic source of −63.0 ‰, consistent with dominant biogenic sources in Arctic regions. The model reproduces this average signature well at −59.5 ‰. Observations highlight a strong synoptic variability in isotopic ratios in the Arctic, with a standard deviation of 0.50 ‰ and a range of 2 ‰. Most of this is missed by the model (see Fig. 2a, prior simulation). Simulated ratios with fixed (temporally and spatially) isotopic signatures for the emission sectors detailed in Sect. 2.3 barely exhibit any variations. The prior standard deviation is 0.22 ‰ (or 0.12 ‰ when removing the wetland event on 21 August), with a range of 1.5 ‰ (or 0.5 ‰). Considering the good fit of simulations to observations of total CH4, the missing variability indicates that the classical assumption of uniform signatures for given sectors and regions is not valid in the Arctic, consistent with and . Contributions to modelled concentrations from different regions of a given emission sector can change much more than the variability in total CH4, as indicated in Fig. 2. For instance, on 22 July, contributions from wetlands turn from a dominating Siberian influence to a North American one, causing a change of ∼30ppb in the signal. Differences in the average wetland source signatures between these two regions of ∼20 ‰ (as suggested by Ganesan et al.2018) would thus translate into ∼0.3 ‰ in measured isotopic ratio, partly explaining the corresponding observed event (see Fig. 2b).

Figure 2(a) Observed and simulated total CH4 concentrations. (b) Simulated contributions to total CH4 concentrations. Individual regions simulated by the model (see Fig. 1) are aggregated into two main continental components: North America (NA) and Eurasia. Light green areas depict Eurasian (mostly Siberian) wetlands, while dark green ones are North American wetlands. Shaded blue areas represent contributions from the sides of the CHIMERE simulation domain (see Fig. 1). Orange shades represent minor anthropogenic contribution. So-called “top” line gives the simulated concentrations originating from the lower stratosphere (i.e. from the top of CHIMERE simulation domain). Please note the gap in y axis scale at 1700 ppb, highlighted by the dashed line. (c) Observed and simulated isotopic ratios before (prior) and after (posterior) inversion.

Figure 3(a) Map of regions constrained by the observations in the inversion. For land regions, only wetlands are constrained. The green region corresponds to oceanic sources from the East Siberian Arctic Shelf. The ship path is indicated in black, with red points highlighting locations with available isotope observations. (b) Distribution of posterior hourly signatures as deduced by the inversion for regions constrained by the observations for the ensemble of 8000 Monte Carlo inversions (see details in Sect. 2.3). Prior signature distributions (dashed lines) are those of Table 1. The optimal posterior signatures, defined as the maximum of the posterior distribution, are highlighted by plain horizontal lines.

Still, more critical for the composition of air masses are the changes in very large-scale hemispheric contributions. As indicated by the blue shading in Fig. 2b, depending on the dominant large-scale transport patterns, contributions from the stratosphere and from the model lateral sides (located in the Tropics) can vary by more than 400 ppb within a few days. This corresponds to dominantly updraught or downdraught transport patterns, as illustrated by Fig. S2 in the Supplement. These very strong variations in total CH4 enhance the impact of uncertainties in the vertical and horizontal distribution of isotopic ratios at the hemispheric scale. First, tropical air masses are influenced by tropical wetlands and anthropogenic emissions, causing a spatial and temporal variability in tropical isotopic ratio of up to 1 ‰, which is not accounted for in our CHIMERE set-up with fixed isotopic ratios at the simulation domain sides (see Sect. 2.2). Second, the vertical profiles of isotopic ratios in the Arctic (see simulated example from the global transport model LMDz in Fig. S3 in the Supplement) are very steep. Such gradients are poorly represented in most global models due to issues in the representation of the vertical transport or to the insufficiently quantified fractionating OH and chlorine sinks in the stratosphere and upper troposphere. These two sources of uncertainties in chemistry-transport models coupled with the strong real-world variations in stratospheric and tropospheric contributions could explain why the regional model CHIMERE does not reproduce the strong synoptic variability in δ13CCH4 observed during the SWERUS-C3 campaign. In particular, for the above-mentioned event on 22 July, contributions from the domain sides vary by more than 300 ppb. Such a variability in CH4 contributions, associated with differences of a few percent per mille between the isotopic ratios of lower-stratosphere air masses and mid-latitude and low-latitude air masses, could explain the observed event.

Thus, the first-order variability in atmospheric isotopic ratios is due to a balance between non-regional transport-related hemispheric features and regional contributions of wetland, ocean and anthropogenic emissions.

## 3.2 Optimization of Arctic source signatures

Assuming that the mix of CH4 sources is correct, we now attempt to separate hemispheric and regional contributions by optimizing source signatures for a set of geographical regions and different emission sectors in the Arctic as detailed in Sect. 2.3. Posterior isotopic ratios in Fig. 2c follow most of the variability in observations, indicating that the inverse method does fit the observations in a satisfying way. The rest of the signal is within the observation uncertainties of 0.1 ‰. This proves that even though the model is not perfect in representing the transport, it is reasonable to use simulated contributions to optimize isotopic signatures.

Figure 3 shows the posterior signature distributions as deduced from the 8000 Monte Carlo inversions for the four regions that are the most constrained by the observations, i.e. those weighted by the sensitivity matrix as detailed in Sect. 2.3. Accounting for the sensitivity matrix, it appears that only the roof boundary conditions (i.e. air masses from the lower stratosphere), ESAS emissions (i.e. emissions from the Laptev, East Siberian and Chukchi seas) and wetland regions on the shores of the Arctic Ocean are reasonably constrained by the SWERUS-C3 ship-based campaign. Even though anthropogenic emissions were optimized in our system, only the wetland emission sector is significantly constrained for land regions (Fig. 3). The lower-stratosphere signatures are in the short range of 48.5 to −46.5 ‰. Wetlands are suggested to have a heavier signature in Canada (optimal signature of −69.9 ‰) than in eastern Siberia (optimal signature of −65.9 ‰, with a node of similar importance at −55 ‰), consistent with and the compilation by . Wetlands in Alaska exhibit a narrow posterior distribution at −51.3 ‰, with a secondary mode at −75 ‰. Alaska is thus well constrained by the inversion. However, the final value may suggest that the inversion has difficulties in differentiating collocated emissions and mixes the signal due to thermogenic sources with co-located wetland emissions, as is the case in Alaska with extensive extraction of raw oil and gas.

Posterior ESAS signatures are significantly shifted by more than 5 ‰ to −49.5 ‰ from the prior signature towards lighter values. This compares with previous studies and points towards a mix of different processes taking place in the Arctic Shelf, such as inputs from the seabed . The posterior signature could thus be explained by mixed biogenic and thermogenic sources, confirming that ESAS emissions, possibly including a hydrate contribution, are not as depleted as wetland sources .

Overall, the approach developed here reveals that the spatial and temporal variations in isotopic source signatures must be accounted for in order to properly represent δ13CCH4 observations. Such an approach does not allow us to reach definitive conclusions when considering the spread of the inferred regional isotopic signatures. However, it is crucial to account for isotopic ratios to avoid misallocating methane flux variations in methane inversions. We also show that atmospheric δ13CCH4 signals can be significant (larger than observation errors), indicating a good potential for the use of isotopic observations based on an oceanic campaign to improve our knowledge of the Arctic methane cycle. Finally, the weight of the boundary conditions in the signal points towards necessary progress in global simulations (including fractionating chemical reactions in the stratosphere) of CH4 atmospheric isotopic ratios.

4 Conclusions

Observations of total atmospheric methane and isotopic ratio were carried out in summer 2014 in the Arctic Ocean during the SWERUS-C3 campaign aboard the Swedish icebreaker Oden. A unique continuous dataset of 45 d of atmospheric isotopic ratios over the Arctic Ocean is available from this campaign. Consistent with other campaigns in the region collecting flasks, the synoptic variability in atmospheric isotopic ratios in the Arctic is very strong, at ∼2 ‰, largely above observation error. Using forward simulations, we confirmed that the assumption of uniform isotopic signatures to represent emission sectors is invalid in the Arctic areas dominated by natural sources. We also exhibited the strong dependency of atmospheric isotopic ratios on large-scale changes in air mass origin (lateral boundaries of our simulation domain, corresponding to mid-latitude and low-latitude air masses; top boundaries corresponding to lower-stratosphere air masses). Based on a simplified inversion framework, the SWERUS-C3 data were used to infer isotopic source signatures of the Arctic regions and emission sectors. Due to the limited number of available observations and the important distance between sources and observations, our system was not able to provide any significant constraints on anthropogenic emissions and could optimize signatures from the ESAS and wetlands near the Arctic Ocean shores only. Wetland and oceanic ESAS source signatures were found to span a very wide range with a multimodal distribution for wetlands. The inversion also indicated that CH4 emissions from the ESAS are composed of a mixture of dominant thermogenic methane, complemented by some biogenic methane.

Overall, only a strong spatial and temporal variability in emission signatures and in stratospheric isotopic ratios can explain the variability in observations. Therefore, our study points towards necessary improvements in simulating the first-order transport and chemistry of methane and its isotopes to reproduce large-scale hemispheric features, especially stratosphere to troposphere exchanges. This makes it necessary to improve (i) the quality of continuous isotopic measurements to capture the synoptic signal with even higher confidence, (ii) numerical chemistry-transport models so that the uncertainties in the first-order processes are at least 1 order of magnitude smaller than the regional signal, which is not the case in our study, and (iii) the mapping of isotopic emission signatures used as priors in inversions as initiated by .

Data availability
Data availability.

All observational data used in this work are publicly available on the Bolin Centre Database (https://bolin.su.se/data/?s=SWERUS-C3, last access: 30 March 2020) alongside complementary data during the SWERUS-C3 campaign. They are divided into three separate sub-datasets: (i) atmospheric measurements for total CH4 during the first leg of the SWERUS-C3 campaign (https://bolin.su.se/data/thornton-2016; ), (ii) atmospheric measurements for total CH4 during the second leg of the campaign (https://bolin.su.se/data/swerus-2014-ghg; ), and (iii) carbon isotopic ratio in atmospheric methane (https://bolin.su.se/data/swerus-2014-d13c; ).

Supplement
Supplement.

Author contributions
Author contributions.

AB, TT and IP designed the simulation experiments. AB and TT developed the code and performed the CHIMERE numerical simulations. TH and JT ran global simulations. PMC and BT designed, carried out and provided observation data from the SWERUS-C3 campaign. PB, MS and JDP contributed to the scientific analysis of this work. AB prepared the paper, with contributions from all co-authors.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

We thank the crew of I/B Oden, who made the SWERUS-C3 expedition possible. SWERUS-C3 funding was provided by the Knut and Alice Wallenberg Foundation, Vetenskapsrådet (Swedish Research Council), Stockholm University, Swedish Polar Research Secretariat, and the Bolin Centre for Climate Research. δ13CCH4 observations from SWERUS-C3 campaign will be posted to the Bolin Centre database at https://bolin.su.se/data/. This work has been supported by the Swedish Research Council VR through a Franco-Swedish project called IZOMET-FS: “Distinguishing Arctic CH4 sources to the atmosphere using inverse analysis of high-frequency CH4, δ13CCH4 and CH3D measurements” (grant no. VR 2014-6584). The study extensively relies on the meteorological data provided by the ECMWF. Calculations were performed using the computing resources of LSCE, maintained by François Marabelle and the LSCE IT team.

Financial support
Financial support.

This research has been supported by the Franco-Swedish project IZOMET-FS (grant no. VR 2014-6584).

Review statement
Review statement.

This paper was edited by Eliza Harris and reviewed by two anonymous referees.

References

Arora, V. K., Berntsen, T., Biastock, A., Bousquet, P., Bruhwiler, L., Bush, E., Chan, E., Christensen, T. R., Dlugokencky, E., Fisher, R. E., France, J., Gauss, M., Höglund-Isaksson, L., Houweling, S., Huissteden, K., and Janssens-Maenhout, G.: AMAP, 2015. AMAP Assessment 2015: Methane as an Arctic climate forcer, Tech. Rep., Arctic Monitoring and Assessment Programme (AMAP), Oslo, Norway, available at: https://www.amap.no/documents/doc/amap-assessment-2015-methane-as-an-arctic-climate-forcer/1285 (last access: 9 March 2020), 2015. a

Arshinov, M. Y., Belan, B. D., Davydov, D. K., Inouye, G., Krasnov, O. A., Maksyutov, S., Machida, T., Fofonov, A. V., and Shimoyama, K.: Spatial and temporal variability of CO2 and CH4 concentrations in the surface atmospheric layer over West Siberia, Atmos. Ocean Opt., 22, 84–93, https://doi.org/10.1134/S1024856009010126, 2009. a

Berchet, A., Pison, I., Chevallier, F., Paris, J.-D., Bousquet, P., Bonne, J.-L., Arshinov, M. Y., Belan, B. D., Cressot, C., Davydov, D. K., Dlugokencky, E. J., Fofonov, A. V., Galanin, A., Lavriĉ, J., Machida, T., Parker, R., Sasakawa, M., Spahni, R., Stocker, B. D., and Winderlich, J.: Natural and anthropogenic methane fluxes in Eurasia: a mesoscale quantification by generalized atmospheric inversion, Biogeosciences, 12, 5393–5414, https://doi.org/10.5194/bg-12-5393-2015, 2015. a, b

Berchet, A., Bousquet, P., Pison, I., Locatelli, R., Chevallier, F., Paris, J.-D., Dlugokencky, E. J., Laurila, T., Hatakka, J., Viisanen, Y., Worthy, D. E. J., Nisbet, E., Fisher, R., France, J., Lowry, D., Ivakhov, V., and Hermansen, O.: Atmospheric constraints on the methane emissions from the East Siberian Shelf, Atmos. Chem. Phys., 16, 4147–4157, https://doi.org/10.5194/acp-16-4147-2016, 2016. a, b, c, d

Bohn, T. J., Melton, J. R., Ito, A., Kleinen, T., Spahni, R., Stocker, B. D., Zhang, B., Zhu, X., Schroeder, R., Glagolev, M. V., Maksyutov, S., Brovkin, V., Chen, G., Denisov, S. N., Eliseev, A. V., Gallego-Sala, A., McDonald, K. C., Rawlins, M. A., Riley, W. J., Subin, Z. M., Tian, H., Zhuang, Q., and Kaplan, J. O.: WETCHIMP-WSL: intercomparison of wetland methane emissions models over West Siberia, Biogeosciences, 12, 3321–3349, https://doi.org/10.5194/bg-12-3321-2015, 2015. a

Cardinali, C., Pezzulli, S., and Andersson, E.: Influence-matrix diagnostic of a data assimilation system, Q. J. R. Meteorol. Soc., 130, 2767–2786, https://doi.org/10.1256/qj.03.205, 2004. a

Charkin, A. N., Dudarev, O. V., Semiletov, I. P., Kruhmalev, A. V., Vonk, J. E., Sánchez-García, L., Karlsson, E., and Gustafsson, Ö.: Seasonal and interannual variability of sedimentation and organic matter distribution in the Buor-Khaya Gulf: the primary recipient of input from Lena River and coastal erosion in the southeast Laptev Sea, Biogeosciences, 8, 2581–2594, https://doi.org/10.5194/bg-8-2581-2011, 2011. a

Craig, H.: Isotopic standards for carbon and oxygen and correction factors for mass-spectrometric analysis of carbon dioxide, Geochim. Cosmochim. Acta, 12, 133–149, https://doi.org/10.1016/0016-7037(57)90024-8, 1957. a

Cramer, B., Poelchau, H. S., Gerling, P., Lopatin, N. V., and Littke, R.: Methane released from groundwater: the source of natural gas accumulations in northern West Siberia, Mar. Petrol. Geol., 16, 225–244, https://doi.org/10.1016/S0264-8172(98)00085-3, 1999. a

Crippa, M., Janssens-Maenhout, G., Dentener, F., Guizzardi, D., Sindelarova, K., Muntean, M., Van Dingenen, R., and Granier, C.: Forty years of improvements in European air quality: regional policy-industry interactions with global impacts, Atmos. Chem. Phys., 16, 3825–3841, https://doi.org/10.5194/acp-16-3825-2016, 2016. a

Dlugokencky, E. J., Bruhwiler, L., White, J. W. C., Emmons, L. K., Novelli, P. C., Montzka, S. A., Masarie, K. A., Lang, P. M., Crotwell, A. M., Miller, J. B., and Gatti, L. V.: Observational constraints on recent increases in the atmospheric CH4 burden, Geophys. Res. Lett., 36, L18803, https://doi.org/10.1029/2009GL039780, 2009. a

Dlugokencky, E. J., Crotwell, A. M., Lang, P. M., and Masarie, K. A.: Atmospheric methane dry air mole Fractions from quasi-continuous measurements at Barrow, Alaska and Mauna Loa, Hawaii, 1986–2013, available at: ftp://ftp.cmdl.noaa.gov/data/trace_gases/ch4/in-situ/surface/ (last access: 9 March 2020), 2014. a

Dmitrenko, I. A., Kirillov, S. A., Tremblay, L. B., Kassens, H., Anisimov, O. A., Lavrov, S. A., Razumov, S. O., and Grigoriev, M. N.: Recent changes in shelf hydrography in the Siberian Arctic: Potential for subsea permafrost instability, J. Geophys. Res.-Oceans, 116, C10027, https://doi.org/10.1029/2011JC007218, 2011. a

Fisher, R. E., Sriskantharajah, S., Lowry, D., Lanoisellé, M., Fowler, C. M. R., James, R. H., Hermansen, O., Lund Myhre, C., Stohl, A., Greinert, J., Nisbet-Jones, P. B. R., Mienert, J., and Nisbet, E. G.: Arctic methane sources: Isotopic evidence for atmospheric inputs, Geophys. Res. Lett., 38, L21803, https://doi.org/10.1029/2011GL049319, 2011. a, b

Fisher, R. E., France, J. L., Lowry, D., Lanoisellé, M., Brownlow, R., Pyle, J. A., Cain, M., Warwick, N., Skiba, U. M., Drewer, J., Dinsmore, K. J., Leeson, S. R., Bauguitte, S. J.-B., Wellpott, A., O'Shea, S. J., Allen, G., Gallagher, M. W., Pitt, J., Percival, C. J., Bower, K., George, C., Hayman, G. D., Aalto, T., Lohila, A., Aurela, M., Laurila, T., Crill, P. M., McCalley, C. K., and Nisbet, E. G.: Measurement of the 13C isotopic signature of methane emissions from northern European wetlands, Global Biogeochem. Cy., 31, 605–623, https://doi.org/10.1002/2016GB005504, 2017. a, b

Folberth, G. A., Hauglustaine, D. A., Lathière, J., and Brocheton, F.: Interactive chemistry in the Laboratoire de Météorologie Dynamique general circulation model: model description and impact analysis of biogenic hydrocarbons on tropospheric chemistry, Atmos. Chem. Phys., 6, 2273–2319, https://doi.org/10.5194/acp-6-2273-2006, 2006. a

Ganesan, A. L., Stell, A. C., Gedney, N., Comyn-Platt, E., Hayman, G., Rigby, M., Poulter, B., and Hornibrook, E. R. C.: Spatially Resolved Isotopic Source Signatures of Wetland Methane Emissions, Geophys. Res. Lett., 45, 3737–3745, https://doi.org/10.1002/2018GL077536, 2018. a, b, c, d, e

Hauglustaine, D. A., Hourdin, F., Jourdain, L., Filiberti, M.-A., Walters, S., Lamarque, J.-F., and Holland, E. A.: Interactive chemistry in the Laboratoire de Météorologie Dynamique general circulation model: Description and background tropospheric chemistry evaluation, J. Geophys. Res.-Atmos., 109, D04314, https://doi.org/10.1029/2003JD003957, 2004. a

Ishizawa, M., Chan, D., Worthy, D., Chan, E., Vogel, F., and Maksyutov, S.: Analysis of atmospheric CH4 in Canadian Arctic and estimation of the regional CH4 fluxes, Atmos. Chem. Phys., 19, 4637–4658, https://doi.org/10.5194/acp-19-4637-2019, 2019. a, b

James, R. H., Bousquet, P., Bussmann, I., Haeckel, M., Kipfer, R., Leifer, I., Niemann, H., Ostrovsky, I., Piskozub, J., Rehder, G., Treude, T., Vielstädte, L., and Greinert, J.: Effects of climate change on methane emissions from seafloor sediments in the Arctic Ocean: A review, Limnol. Oceanogr., 61, S283–S299, https://doi.org/10.1002/lno.10307, 2016. a

Kirschke, S., Bousquet, P., Ciais, P., Saunois, M., Canadell, J. G., Dlugokencky, E. J., Bergamaschi, P., Bergmann, D., Blake, D. R., Bruhwiler, L., Cameron-Smith, P., Castaldi, S., Chevallier, F., Feng, L., Fraser, A., Heimann, M., Hodson, E. L., Houweling, S., Josse, B., Fraser, P. J., Krummel, P. B., Lamarque, J.-F., Langenfelds, R. L., Le Quéré, C., Naik, V., O'Doherty, S., Palmer, P. I., Pison, I., Plummer, D., Poulter, B., Prinn, R. G., Rigby, M., Ringeval, B., Santini, M., Schmidt, M., Shindell, D. T., Simpson, I. J., Spahni, R., Steele, L. P., Strode, S. A., Sudo, K., Szopa, S., van der Werf, G. R., Voulgarakis, A., van Weele, M., Weiss, R. F., Williams, J. E., and Zeng, G.: Three decades of global methane sources and sinks, Nat. Geosci., 6, 813–823, https://doi.org/10.1038/ngeo1955, 2013. a

Lorenson, T.: Gas composition and isotopic geochemistry of cuttings, core, and gas hydrate from the JAPEX/JNOC/GSC Mallik 2L-38 gas hydrate research well, Bulletin of the Geological Survey of Canada, p. 21, available at: http://pubs.er.usgs.gov/publication/70021662 (last access: 9 March 2020), 1999. a

McCalley, C. K., Woodcroft, B. J., Hodgkins, S. B., Wehr, R. A., Kim, E.-H., Mondav, R., Crill, P. M., Chanton, J. P., Rich, V. I., Tyson, G. W., and Saleska, S. R.: Methane dynamics regulated by microbial community response to permafrost thaw, Nature, 514, 478–481, https://doi.org/10.1038/nature13798, 2014. a

McGuire, A. D., Anderson, L. G., Christensen, T. R., Dallimore, S., Guo, L., Hayes, D. J., Heimann, M., Lorenson, T. D., Macdonald, R. W., and Roulet, N.: Sensitivity of the carbon cycle in the Arctic to climate change, Ecol. Monogr., 79, 523–555, https://doi.org/10.1890/08-2025.1, 2009. a

Menut, L., Bessagnet, B., Khvorostyanov, D., Beekmann, M., Blond, N., Colette, A., Coll, I., Curci, G., Foret, G., Hodzic, A., Mailler, S., Meleux, F., Monge, J.-L., Pison, I., Siour, G., Turquety, S., Valari, M., Vautard, R., and Vivanco, M. G.: CHIMERE 2013: a model for regional atmospheric composition modelling, Geosci. Model Dev., 6, 981–1028, https://doi.org/10.5194/gmd-6-981-2013, 2013. a

Nisbet, E. G., Dlugokencky, E. J., Manning, M. R., Lowry, D., Fisher, R. E., France, J. L., Michel, S. E., Miller, J. B., White, J. W. C., Vaughn, B., Bousquet, P., Pyle, J. A., Warwick, N. J., Cain, M., Brownlow, R., Zazzeri, G., Lanoisellé, M., Manning, A. C., Gloor, E., Worthy, D. E. J., Brunke, E.-G., Labuschagne, C., Wolff, E. W., and Ganesan, A. L.: Rising atmospheric methane: 2007–2014 growth and isotopic shift, Global Biogeochem. Cy., 30, 1356–1370, https://doi.org/10.1002/2016GB005406, 2016. a

Nisbet, E. G., Manning, M. R., Dlugokencky, E. J., Fisher, R. E., Lowry, D., Michel, S. E., Myhre, C. L., Platt, S. M., Allen, G., Bousquet, P., Brownlow, R., Cain, M., France, J. L., Hermansen, O., Hossaini, R., Jones, A. E., Levin, I., Manning, A. C., Myhre, G., Pyle, J. A., Vaughn, B. H., Warwick, N. J., and White, J. W. C.: Very Strong Atmospheric Methane Growth in the 4 Years 2014–2017: Implications for the Paris Agreement, Global Biogeochem. Cy., 33, 318–342, https://doi.org/10.1029/2018GB006009, 2019. a

Pankratova N., Skorokhod A., Belikov I., Elansky N., Rakitin V., Shtabkin Y., and Berezina, E.: Evidence of atmospheric response to methane emissions from the east siberian arctic shelf. geography, environment, sustainability, 11, 85–92, https://doi.org/10.24057/2071-9388-2018-11-1-85-92, 2018. a

Pankratova, N., Belikov, I., Skorokhod, A., Belousov, V., Artamonov, A., Repina, I., and Shishov, E.: Measurements and data processing of atmospheric CO2, CH4, H2O and δ13C-CH4 mixing ratio during the ship campaign in the East Arctic and the Far East seas in autumn 2016, IOP Conf. Ser.: Earth Environ. Sci., 231, 012041, https://doi.org/10.1088/1755-1315/231/1/012041, 2019. a, b

Paris, J. D., Ciais, P., Nédélec, P., Stohl, A., Belan, B. D., Arshinov, M. Y., Carouge, C., Golitsyn, G. S., and Granberg, I. G.: New insights on the chemical composition of the Siberian air shed from the YAK-AEROSIB aircraft campaigns, B. Am. Meteorol. Soc., 91, 625–641, 2010. a

Pisso, I., Myhre, C. L., Platt, S. M., Eckhardt, S., Hermansen, O., Schmidbauer, N., Mienert, J., Vadakkepuliyambatta, S., Bauguitte, S., Pitt, J., Allen, G., Bower, K. N., O'Shea, S., Gallagher, M. W., Percival, C. J., Pyle, J., Cain, M., and Stohl, A.: Constraints on oceanic methane emissions west of Svalbard from atmospheric in situ measurements and Lagrangian transport modeling, J. Geophys. Res.-Atmos., 121, 025590, https://doi.org/10.1002/2016JD025590, 2016. a, b

Ringeval, B., de Noblet-Ducoudré, N., Ciais, P., Bousquet, P., Prigent, C., Papa, F., and Rossow, W. B.: An attempt to quantify the impact of changes in wetland extent on methane emissions on the seasonal and interannual time scales, Global Biogeochem. Cy., 24, GB2003, https://doi.org/10.1029/2008GB003354, 2010. a

Ruppel, C. D. and Kessler, J. D.: The interaction of climate change and methane hydrates, Rev. Geophys., 55, 126–168, https://doi.org/10.1002/2016RG000534, 2017. a

Sapart, C. J., Shakhova, N., Semiletov, I., Jansen, J., Szidat, S., Kosmach, D., Dudarev, O., van der Veen, C., Egger, M., Sergienko, V., Salyuk, A., Tumskoy, V., Tison, J.-L., and Röckmann, T.: The origin of methane in the East Siberian Arctic Shelf unraveled with triple isotope analysis, Biogeosciences, 14, 2283–2292, https://doi.org/10.5194/bg-14-2283-2017, 2017. a, b

Sasakawa, M., Shimoyama, K., Machida, T., Tsuda, N., Suto, H., Arshinov, M., Davydov, D., Fofonov, A., Krasnov, O., Saeki, T., Koyama, Y., and Maksyutov, S.: Continuous measurements of methane from a tower network over Siberia, Tellus B, 62, 403–416, https://doi.org/10.1111/j.1600-0889.2010.00494.x, 2010. a

Saunois, M., Bousquet, P., Poulter, B., Peregon, A., Ciais, P., Canadell, J. G., Dlugokencky, E. J., Etiope, G., Bastviken, D., Houweling, S., Janssens-Maenhout, G., Tubiello, F. N., Castaldi, S., Jackson, R. B., Alexe, M., Arora, V. K., Beerling, D. J., Bergamaschi, P., Blake, D. R., Brailsford, G., Brovkin, V., Bruhwiler, L., Crevoisier, C., Crill, P., Covey, K., Curry, C., Frankenberg, C., Gedney, N., Höglund-Isaksson, L., Ishizawa, M., Ito, A., Joos, F., Kim, H.-S., Kleinen, T., Krummel, P., Lamarque, J.-F., Langenfelds, R., Locatelli, R., Machida, T., Maksyutov, S., McDonald, K. C., Marshall, J., Melton, J. R., Morino, I., Naik, V., O'Doherty, S., Parmentier, F.-J. W., Patra, P. K., Peng, C., Peng, S., Peters, G. P., Pison, I., Prigent, C., Prinn, R., Ramonet, M., Riley, W. J., Saito, M., Santini, M., Schroeder, R., Simpson, I. J., Spahni, R., Steele, P., Takizawa, A., Thornton, B. F., Tian, H., Tohjima, Y., Viovy, N., Voulgarakis, A., van Weele, M., van der Werf, G. R., Weiss, R., Wiedinmyer, C., Wilton, D. J., Wiltshire, A., Worthy, D., Wunch, D., Xu, X., Yoshida, Y., Zhang, B., Zhang, Z., and Zhu, Q.: The global methane budget 2000–2012, Earth Syst. Sci. Data, 8, 697–751, https://doi.org/10.5194/essd-8-697-2016, 2016. a, b

Saunois, M., Bousquet, P., Poulter, B., Peregon, A., Ciais, P., Canadell, J. G., Dlugokencky, E. J., Etiope, G., Bastviken, D., Houweling, S., Janssens-Maenhout, G., Tubiello, F. N., Castaldi, S., Jackson, R. B., Alexe, M., Arora, V. K., Beerling, D. J., Bergamaschi, P., Blake, D. R., Brailsford, G., Bruhwiler, L., Crevoisier, C., Crill, P., Covey, K., Frankenberg, C., Gedney, N., Höglund-Isaksson, L., Ishizawa, M., Ito, A., Joos, F., Kim, H.-S., Kleinen, T., Krummel, P., Lamarque, J.-F., Langenfelds, R., Locatelli, R., Machida, T., Maksyutov, S., Melton, J. R., Morino, I., Naik, V., O'Doherty, S., Parmentier, F.-J. W., Patra, P. K., Peng, C., Peng, S., Peters, G. P., Pison, I., Prinn, R., Ramonet, M., Riley, W. J., Saito, M., Santini, M., Schroeder, R., Simpson, I. J., Spahni, R., Takizawa, A., Thornton, B. F., Tian, H., Tohjima, Y., Viovy, N., Voulgarakis, A., Weiss, R., Wilton, D. J., Wiltshire, A., Worthy, D., Wunch, D., Xu, X., Yoshida, Y., Zhang, B., Zhang, Z., and Zhu, Q.: Variability and quasi-decadal changes in the methane budget over the period 2000–2012, Atmos. Chem. Phys., 17, 11135–11161, https://doi.org/10.5194/acp-17-11135-2017, 2017. a

Shakhova, N., Semiletov, I., Salyuk, A., Yusupov, V., Kosmach, D., and Gustafsson, Ö.: Extensive Methane Venting to the Atmosphere from Sediments of the East Siberian Arctic Shelf, Science, 327, 1246–1250, https://doi.org/10.1126/science.1182221, 2010. a

Shakhova, N., Semiletov, I., Leifer, I., Sergienko, V., Salyuk, A., Kosmach, D., Chernykh, D., Stubbs, C., Nicolsky, D., Tumskoy, V., and Gustafsson, Ö.: Ebullition and storm-induced methane release from the East Siberian Arctic Shelf, Nat. Geosci., 7, 64–70, https://doi.org/10.1038/ngeo2007, 2014. a

Sherwood, O. A., Schwietzke, S., Arling, V. A., and Etiope, G.: Global Inventory of Gas Geochemistry Data from Fossil Fuel, Microbial and Burning Sources, version 2017, Earth Syst. Sci. Data, 9, 639–656, https://doi.org/10.5194/essd-9-639-2017, 2017. a, b, c

Skorokhod, A. I., Pankratova, N. V., Belikov, I. B., Thompson, R. L., Novigatsky, A. N., and Golitsyn, G. S.: Observations of atmospheric methane and its stable isotope ratio (δ13C) over the Russian Arctic seas from ship cruises in the summer and autumn of 2015, Dokl. Earth Sc., 470, 1081–1085, https://doi.org/10.1134/S1028334X16100160, 2016. a

Tan, Z., Zhuang, Q., Henze, D. K., Frankenberg, C., Dlugokencky, E., Sweeney, C., Turner, A. J., Sasakawa, M., and Machida, T.: Inverse modeling of pan-Arctic methane emissions at high spatial resolution: what can we learn from assimilating satellite retrievals and using different process-based wetland and lake biogeochemical models?, Atmos. Chem. Phys., 16, 12649–12666, https://doi.org/10.5194/acp-16-12649-2016, 2016. a

Tarantola, A.: Inverse problem theory and methods for model parameter estimation, Vol. 89, siam, 2005. a

Tarasova, O., Brenninkmeijer, C., Assonov, S., Elansky, N., Röckmann, T., and Brass, M.: Atmospheric CH4 along the Trans-Siberian railroad (TROICA) and river Ob: Source identification using stable isotope analysis, Atmos. Environ., 40, 5617–5628, https://doi.org/10.1016/j.atmosenv.2006.04.065, 2006. a, b

Tarasova, O. A., Houweling, S., Elansky, N., and Brenninkmeijer, C. A. M.: Application of stable isotope analysis for improved understanding of the methane budget: comparison of TROICA measurements with TM3 model simulations, J. Atmos. Chem., 63, 49–71, https://doi.org/10.1007/s10874-010-9157-y, 2009. a

Thompson, R. L., Sasakawa, M., Machida, T., Aalto, T., Worthy, D., Lavric, J. V., Lund Myhre, C., and Stohl, A.: Methane fluxes in the high northern latitudes for 2005–2013 estimated using a Bayesian atmospheric inversion, Atmos. Chem. Phys., 17, 3553–3572, https://doi.org/10.5194/acp-17-3553-2017, 2017. a

Thonat, T., Saunois, M., Bousquet, P., Pison, I., Tan, Z., Zhuang, Q., Crill, P. M., Thornton, B. F., Bastviken, D., Dlugokencky, E. J., Zimov, N., Laurila, T., Hatakka, J., Hermansen, O., and Worthy, D. E. J.: Detectability of Arctic methane sources at six sites performing continuous atmospheric measurements, Atmos. Chem. Phys., 17, 8371–8394, https://doi.org/10.5194/acp-17-8371-2017, 2017. a

Thornton, B. F., Geibel, M. C., Crill, P. M., Humborg, C., and Mörth, C.-M.: Methane fluxes from the sea to the atmosphere across the Siberian shelf seas, Geophys. Res. Lett., 43, 5869–5877, https://doi.org/10.1002/2016GL068977, 2016a. a, b, c, d, e

Thornton, B. F., Wik, M., and Crill, P. M.: Double-counting challenges the accuracy of high-latitude methane inventories, Geophys. Res. Lett., 43, 12569–12577, https://doi.org/10.1002/2016GL071772, 2016b. a, b, c

Thornton, B. F., Prytherch, J., Andersson, K., Brooks, I. M., Salisbury, D., Tjernström, M., and Crill, P. M.: Shipborne eddy covariance observations of methane fluxes constrain Arctic sea emissions, Sci. Adv., 6, eaay7934, https://doi.org/10.1126/sciadv.aay7934, 2020. a, b

Thornton, B., Geibel, M., Crill, P., and Humborg, C.: Atmospheric methane and carbon dioxide and surface water methane from the SWERUS-C3 Arctic Ocean expedition in 2014, leg 1, Dataset version 2.0, Bolin Centre Database, https://doi.org/10.17043/thornton-2016, 2020a. a

Thornton, B., Andersson, K., and Crill, P.: Atmospheric methane and carbon dioxide from the SWERUS-C3 Arctic Ocean expedition in 2014, leg 2, Dataset version 1.0, Bolin Centre Database, https://doi.org/10.17043/swerus-2014-ghg, 2020b. a

Thornton, B., Andersson, K., and Crill, P.: Carbon isotope ratio (δ13C) in atmospheric methane from the SWERUS-C3 Arctic Ocean expedition in 2014, Dataset version 1.0, Bolin Centre Database, https://doi.org/10.17043/swerus-2014-d13c, 2020c.  a

Turner, A. J., Frankenberg, C., and Kort, E. A.: Interpreting contemporary trends in atmospheric methane, P. Natl. Acad. Sci. USA, 116, 2805–2813, https://doi.org/10.1073/pnas.1814297116, 2019. a

Warwick, N. J., Cain, M. L., Fisher, R., France, J. L., Lowry, D., Michel, S. E., Nisbet, E. G., Vaughn, B. H., White, J. W. C., and Pyle, J. A.: Using δ13C-CH4 and δD-CH4 to constrain Arctic methane emissions, Atmos. Chem. Phys., 16, 14891–14908, https://doi.org/10.5194/acp-16-14891-2016, 2016. a

Yu, J., Xie, Z., Sun, L., Kang, H., He, P., and Xing, G.: δ13C-CH4 reveals CH4 variations over oceans from mid-latitudes to the Arctic, Sci. Rep., 5, 13760, https://doi.org/10.1038/srep13760, 2015. a, b