Articles | Volume 18, issue 1
Research article
 | Highlight paper
18 Feb 2022
Research article | Highlight paper |  | 18 Feb 2022

Autonomous methane seep site monitoring offshore western Svalbard: hourly to seasonal variability and associated oceanographic parameters

Knut Ola Dølven, Bénédicte Ferré, Anna Silyakova, Pär Jansson, Peter Linke, and Manuel Moser

Improved quantification techniques of natural sources are needed to explain variations in atmospheric methane. In polar regions, high uncertainties in current estimates of methane release from the seabed remain. We present unique 10- and 3-month time series of bottom water measurements of physical and chemical parameters from two autonomous ocean observatories deployed at separate intense seabed methane seep sites (91 and 246 m depth) offshore western Svalbard from 2015 to 2016. Results show high short-term (100–1000 nmol L−1 within hours) and seasonal variation, as well as higher (2–7 times) methane concentrations compared to previous measurements. Rapid variability is explained by uneven distribution of seepage and changing ocean current directions. No overt influence of tidal hydrostatic pressure or water temperature variations on methane concentration was observed, but an observed negative correlation with temperature at the 246 m site fits with hypothesized seasonal blocking of lateral methane pathways in the sediments. Negative correlation between bottom water methane concentration (and variability) and wind forcing, concomitant with signs of weaker water column stratification, indicates increased potential for methane release to the atmosphere in fall and winter. We present new information about short- and long-term methane variability and provide a preliminary constraint on the uncertainties that arise in methane inventory estimates from this variability.

1 Introduction

Unexplained changes in atmospheric methane (CH4) mole fraction motivate research in understanding and quantifying non-anthropogenic sources (Saunois et al.2020). The atmospheric forcing of CH4 is particularly sensitive to changes in emission rates due to its high warming potential and short lifetime. Improved knowledge about atmospheric CH4 fluxes is therefore crucial to constrain future climate projections (Pachauri and Meyer2014; Myhre et al.2016b). These properties of atmospheric CH4 also make reducing anthropogenic CH4 emissions a potential solution for rapid climate change mitigation (Saunois et al.2016). A global effort to cut greenhouse gas emissions through international agreements is, however, dependent on precise estimates of sources and sinks to verify contributions from different nations.

Seabed seepage is considered a minor source of atmospheric CH4, but there is high uncertainty in current and predicted emission estimates (Saunois et al.2016). Current estimates suggest a total contribution of 7 (5–10) Tg yr−1 (Etiope et al.2019; Saunois et al.2020), which is  1 % of the total CH4 emissions to the atmosphere. Methane is released from the seabed as free gas (bubbles) and dissolved gas in sediment pore water. Bubbles rise quickly towards the sea surface, but most CH4 dissolves near the seafloor because of gas exchange across the bubble rims and bubble dissolution (McGinnis et al.2006; Jansson et al.2019a). Dissolved CH4 is dispersed and advected by ocean currents (Silyakova et al.2020) and is continuously transformed to carbon dioxide (CO2) by bacterial aerobic oxidation (Hanson and Hanson1996; Reeburgh2007). These processes significantly limit the lifetime of CH4 in the water column, and the amount of CH4 that can reach the atmosphere is highly dependent on the depth where the seepage occurs (McGinnis et al.2006; Graves et al.2015). Intense CH4 seepage at shallow depths in coastal areas and on continental shelves is therefore the main potential source of seabed CH4 to the atmosphere.

The shallow continental margins of the Arctic Ocean store large amounts of CH4 as free gas, gas dissolved in pore water fluid, and gas hydrates (James et al.2016; Ruppel and Kessler2017), i.e., clathrate structures composed of water trapped by hydrocarbon molecules formed and kept stable at low temperature and high pressure (Sloan1998). Increasing bottom water temperature has the potential to liberate methane from these reservoirs via various mechanisms, potentially resulting in a positive climate feedback loop (Westbrook et al.2009; Shakhova et al.2010; James et al.2016).

Studies on CH4 inventory, distribution and release in the Arctic Ocean are mainly based on research cruise data from late spring to early fall, when ice and weather conditions allow fieldwork in the region (Gentz et al.2014; Sahling et al.2014; Mau et al.2017), whereas winter data are sparse. Bottom water temperature (Westbrook et al.2009; Reagan et al.2011; Ferré et al.2012; Braga et al.2020), water mass origins (Steinle et al.2015), micro-seismicity (Franek et al.2017) and hydrostatic pressure (Linke et al.2009; Römer et al.2016) have all been proposed to be linked with sources and sinks of CH4 in the water column. These processes act on a wide range of timescales, from hours (e.g., hydrostatic pressure) to decades (bottom water temperature). Without a better understanding of the spatial and temporal variability of CH4 in Arctic seep sites, it is challenging to untangle these processes. Unconstrained local variability in CH4 seepage and concentration also imposes a high degree of uncertainty on CH4 inventory estimates (Saunois et al.2020). The combination of climate-sensitive CH4 storages, vast shallow ocean regions and limited data availability highlights the need for more understanding of seabed CH4 seepage on Arctic shelves.

To assess the aforementioned challenges, we have obtained, analyzed and compared two unique long-term underwater multi-parameter time series from seafloor observatories deployed at two distinct intense CH4 seep sites on the western Svalbard continental shelf (Fig. 1) where no CH4 measurements have previously been done in winter season. We combine high-frequency physical (ocean currents, temperature, salinity, pressure) and chemical (O2, CO2, CH4) data to perform hypothesis testing and provide new insights on CH4 distribution, content, and variability on short (minutes) and long (seasonal) timescales and potential implications.

1.1 Regional settings

Two observatories (O91 and O246) were deployed from June 2015 (CAGE 15-3 cruise) to May 2016 (CAGE 16-4 cruise) from R/V Helmer Hanssen at the inter-trough shelf region between Isfjorden and Kongsfjorden, west of Prins Karls Forland. The O91 observatory was deployed at 91 m water depth on the continental shelf (78.561 N, 10.142 E), and the O246 observatory was deployed at 246 m water depth further offshore close to the shelf break (78.655 N, 9.433 E, Fig. 1).

Figure 1Bathymetry of the study area with the locations of the observatories O91 and O246 offshore western Svalbard. Flares detected by single-beam echo sounder survey prior to recovering the observatories (May 2016, cruise CAGE 16-4) are indicated with red dots, and ship tracks are shown as brown lines. The inset map shows the working area (red square) offshore Svalbard. WSC and CC refer to the warm West Spitsbergen Current and cold Coastal Current, respectively.

Both sites were located in areas with thousands of previously mapped CH4 gas seeps (e.g., Sahling et al.2014; Veloso-Alarcón et al.2019; Silyakova et al.2020; this work; see Fig. 1), often referred to as “flares” due to the appearance of bubble streams in echo sounder data. Nonetheless, atmospheric sampling in this region suggests that any emissions to the atmosphere are small (Platt et al.2018). Gas accumulation at the O246 seep site has been suggested to be a result of gas migration in permeable layers within the seabed from deeper free gas or hydrate reservoirs (Rajan et al.2012; Sarkar et al.2012; Veloso-Alarcón et al.2019), while seepage at site O91 has been attributed to thawing sub-sea permafrost due to ice sheet retreat at the end of the last glaciation (Sahling et al.2014; Portnov et al.2016). Water sampling has indicated high temporal variability, with bottom water concentrations (average) changing from 200 nmol L−1 within 1 week in July 2014 at O91 (Myhre et al.2016a) and  80 nmol L−1 within 20 h (two single-point measurements) at O246 in August 2010 (Gentz et al.2014). A consistent pattern of decreasing concentrations from the seafloor to the sea surface at both sites (400 to < 8 nmol L−1 at O91, Myhre et al.2016a, and from > 500 to < 20 nmol L−1 at O246, Gentz et al.2014) has also been observed. Further offshore, continuous measurements from a towed fast-response underwater laser spectrometer also revealed very high spatial CH4 variability (Jansson et al.2019b).

The local water masses are characterized by exchange and convergence of warm, saline Atlantic water (e.g., temperature T>3C and salinity SA>34.9; Swift and Aagaard1981) in the West Spitsbergen Current and colder, fresher Arctic water (e.g., T<0C, 34.3<SA<34.8; Loeng1991) in the Coastal Current combined with seasonal cooling, ice formation and freshwater input from land (Nilsen et al.2016) (Fig. 1). Local mixing rates can be strongly affected by synoptic-scale weather systems, causing upwelling and disruption of the front between the two ocean currents (Saloranta and Svendsen2001; Cottier et al.2007). Freshwater input in summer stratifies the water column, while cooling, storm activity and sea ice formation can facilitate vertical mixing in winter (Saloranta and Svendsen2001; Nilsen et al.2016).

2 Methods

The “K-Lander” ocean observatories were designed to monitor CH4 release and associated physical and chemical parameters in challenging environments (see Appendix A). A launcher equipped with camera and telemetry allowed for safe deployment at a site selected by visual control. Observatory O91 recorded data from 2 July 2015 to 6 May 2016, while O246 recorded data from 1 July until 3 October 2015, when data recording ceased due to an electrical malfunction.

Both observatories were equipped with an acoustic Doppler current profiler (ADCP), a conductivity–temperature–depth (CTD) profiler with an oxygen optode, and Contros HydroC CO2 II and HydroC Plus CH4 sensors (Fig. A1a; details about the instrumentation are provided in Appendix B). The deployed HydroC CH4, being a younger iteration of the sensor, relies on a tunable diode laser absorption spectrometry (TDLAS) detector (rather than non-dispersive infrared spectrometry (NDIR)), while the CO2 sensors use NDIR detectors. Both sensors were equipped with polydimethylsiloxane (PDMS) membranes and Seabird SBE 5M pumps (see Appendix B).

Table 1Correlation coefficients between variables at O91. The terms RTC CH4 and raw CH4 refer to response-time-corrected and untreated CH4 data, respectively (Sect. 2 and Appendix B).

Download Print Version | Download XLSX

The high power consumption of the Contros HydroC CH4 and CO2 sensors required a power cycling mode to allow for long-term monitoring while simultaneously capturing rapid short-term variability. Partial pressure of CH4 and CO2 was therefore measured continuously for 24 h every 21 d and for 1 h every day (see Table B1). Methane concentration data were corrected for slow response time onto a 3 min interval grid and converted to absolute concentration following Dølven et al. (2021), which is the default “CH4 concentration” discussed and described in this text (see Appendix B). Faulty pumps in the CO2 sensors ambiguously increased the response time, which prevented response time correction, making CO2 data suitable only for long-term qualitative analysis.

Uncertainty ranges for the CH4 sensor data are reported as 95  % confidence intervals and typically vary between 5  % and 20  % (Fig. B1b). We did not perform any post-validation and/or intermittent validation. Although always an advantage for all sensors in long-term deployments, such validation is not a requirement for the TDLAS-based sensor (as opposed to NDIR), due to its high long-term stability. Standard post-processing (e.g., inspection of metadata such as internal pressure and temperature) and evaluation of fit residuals in the response time correction procedure (see Appendix B and Dølven et al., 2021) also indicated consistent sensor behavior throughout the deployments. It is also worth noting that the current paper is concerned with large changes and high concentrations, and we are confident that the quality of the response-time-corrected Contros HydroC CH4 data is sufficient to support the inferences described herein.

We calculated correlation coefficient (R) matrices to give a first-order overview of the linear relationships between the measured parameters. We mapped the flares in the area using single-beam echo sounder data collected during the observatory recovery cruise in 2016 (CAGE 16-4, Fig. 1) and estimated gas flow rates using the FlareHunter software (Veloso et al.2015). Additionally, we obtained 10 m wind reanalysis data from the ERA-Interim database.

We calculated seawater density (McDougall and Barker2011) and CH4 solubility (Kossel et al.2013) using the CTD data. A CTD cast (SBE plus 24 Hz) prior to the O91 recovery (6 May 2016) showed a salinity drift in the conductivity sensor of around 0.4 (here and elsewhere in the paper, salinity values are practical salinity). Post-calibration inspection of the conductivity signal and potential water mass mixing end-members indicates that this might have been caused by mud pollution occurring in late 2015 or early 2016.

3 Results

3.1 Time series at site O91

Dissolved CH4 concentration at site O91 ranged from 5 ± 3 nmol L−1 (6 December in 2015) to 1748 ± 142 nmol L−1 (20 August in 2015) (Fig. 2a and Appendix C), with 2.5 and 97.5 percentiles of 16 and 785 nmol L−1. The data follow a nearly lognormal distribution, with a mean and median of 227 and 165 nmol L−1, respectively, and interquartile range of 88–334 nmol L−1. Large variations (> 100 up to almost 1000 nmol L−1) in CH4 concentration occurred on short timescales (< 1 h) throughout the measurement period (see Fig. 2a, d and all 24 h periods in Appendix C), with an average range for all 24 h periods of 840 nmol L−1 and a median rate of change (ROC) of 3.2 nmol L−1 min−1. We also observe a long-term trend of decreasing running median (2-week window) concentrations towards winter, from 495 nmol L−1 in July–August 2015 to 53 nmol L−1 in January 2016 (Fig. 2). There was a relatively weak but significant negative correlation between the wind speed and CH4 concentration (RRTC=-0.33) but otherwise weak to non-existent linear relationships between CH4 concentration and the measured ocean parameters (Table 1).

CO2 averaged 403 µatm with an increase towards mid-November 2015 ( 410 µatm) and then a decrease until 6 May ( 391 µatm) in 2016 (Fig. 2a). CO2 dropped to  305 µatm on 24 August, concurrent with a rapid decrease in salinity (0.5), increase in temperature and oxygen, and high CH4 concentration. The increase in oxygen rules out methanogenesis. Instead, there might be at least two explanations for the reduction of CO2 and enrichment of CH4: (i) water column mixing, which brings oxygen-rich, warm and fresh surface water to deeper depths (and with it CO2-depleted water), or (ii) methane enrichment by zooplankton following the summer bloom.

Figure 2Time series from (a) O91 and (b) O246 showing response-time-corrected (see Appendix B) CH4, CO2, temperature, salinity, pressure, CH4 solubility, oxygen and wind speed (10 m) data. The O246 data are truncated due to an electrical malfunction in the system on 3 October. (c) A 2 d histogram and 1 SD variance ellipse of bottom current velocity (81 m depth at O91 and 236 m depth at O246) and (d) an example of 24 and 1 h (2 and 3 August) CH4 concentration measurement period from O91 (green box). All 24 h measurement periods are shown in Appendix C. Note the different scales between O91 and O246.


Bottom water temperature increased steadily from  3 in July 2015 to  5.5 C in October–November 2015, with occasional sharp shifts (T±1C) occurring within hours to days (Fig. 2b). Temperature then decreased from the beginning of December 2015 to  1.8 C at the end of the deployment in May 2016, showing more frequent and stronger episodes of rapid temperature shifts (T±2C also occurring over hours or days). Despite uncertainty in salinity data, it is worth noting that these rapid shifts in temperature and salinity were reproduced by the Svalbard 800 model in the same area (Silyakova et al.2020) by eddy activity.

Hydrostatic pressure was mostly governed by tides (94.5 % of variance) with a dominant semi-diurnal M2 tide (M2 refers to a tidal constituent with period 12.42 h; see, e.g., Gerkema2019). Amplitudes varied from  1.2 to 1.5 m during neap and spring cycles (Fig. 2c).

The calculated CH4 solubility decreased from 0.016 mol L−1 in July 2015 to 0.015 mol L−1 at the end of November 2015 and increased to almost 0.017 mol L−1 in May 2016 (Fig. 2c). This long-term trend was mainly caused by temperature variability (R=-0.99), while tidal pressure changes caused a semi-diurnal variation of ±0.005 mol L−1.

Dissolved O2 decreased from  385 µmol L−1 in July 2015 to  350 µmol L−1 at the beginning of December 2015 and increased to  400 µmol L−1 towards 6 May 2016 (Fig. 2d) and followed temperature inversely (R=-0.94), with similar long- and short-term variability.

The averaged bottom water current (81 m above the seafloor) was 4 cm s−1 in a northwestward direction (321 N) (Fig. 2c). The current usually had one counterclockwise rotation every 23.93 h period, corresponding to the diurnal K1 tidal constituent (tide with period 23.93 h; see Gerkema2019) with a secondary semi-diurnal (M2) modulation.

3.2 Time series at site O246

CH4 concentration at site O246 ranged from 10 ± 3 nmol L−1 on 21 September 2015 to 2727 ± 182 nmol L−1 on 18 August 2015, with 2.5 and 97.5 percentiles of 107 and 1374 nmol L−1. The data approximately follow lognormal distribution with average and median of 577 and 600 nmol L−1, respectively, and an interquartile range of 293–721 nmol L−1. The median ROC of CH4 (31 nmol L−1 min−1) was almost 20 times higher than site O91 (Fig. 2b and Appendix C). There was also clear diurnal periodicity in CH4 concentration at O246. The long-term trend (2-week running mean) shows decreasing concentrations until 3 October 2015 (end of the measuring period, Fig. 2b). Dissolved O2 decreased from  380 to  300 µmol L−1 and was negatively correlated with water temperature (R=-0.61; see Table 2 for the complete correlation matrix).

Table 2Correlation coefficients between variables at O91. The terms RTC CH4 and raw CH4 refer to response-time-corrected and untreated CH4 (see Sect. 2 and Appendix B).

Download Print Version | Download XLSX

Temperature and salinity increased from  2.5 to  4.0 C and  34.85 up to  35.0, respectively, from the deployment until October 2015 (Fig. 2b), with Atlantic water being dominant throughout the measuring period. Rapid shifts of around ± 1 C and 0.05 salinity occurred occasionally over a period of hours to days.

Variance in hydrostatic pressure was mainly explained by the tides (95.2 %), which were mainly governed by the semi-diurnal M2 tide, with weaker diurnal and fortnightly modulation (Fig. 2b). Changes in pressure varied from  1.2 to  1.5 m during periods of neap and spring tide.

Being governed mainly by temperature (R=-0.99), CH4 solubility dropped from 0.042 to 0.040 mol L−1 from the deployment in July until October 2015, with a semi-diurnal variation of  0.005 mol L−1 due to tidal changes in hydrostatic pressure.

The averaged current was  10 cm s−1 northward (7 N) (Fig. 2c). Variability in the along-slope current (10 N direction) was strongly related to the semi-diurnal M2 tidal component, while the cross-slope currents were governed by the diurnal K1 frequency. The bottom water current rotated counterclockwise with a period of 23.93 h (K1 tidal constituent) and semi-diurnal modulation in the along-slope component. Dissolved CH4 concentration was weakly anti-correlated with wind speed (R=-0.29), temperature (R=-0.31), and salinity (R=-0.24) and positively correlated with CH4 solubility (R=0.33) and oxygen (R=0.3).

4 Discussion

4.1 CH4 variability

Combining mapped flares and flow rates from the recovery cruise (May 2016) with bottom water current velocity (9 m above the seafloor) reveals that CH4 concentration was strongly affected by whether water was advected from areas where we mapped strong or weak seepage in May 2016 (Fig. 3). Strong seeps (flow rate > 200 mL−1 min−1) were mainly located between  30 and 80 m to the north or northeast of site O91, and only weak and more distant seepage was observed southwest of the observatory (Fig. 3a). Consequently, the averaged CH4 concentration from water coming from northeast was  440 nmol L−1, while water from the southwest averaged  100 nmol L−1. Similarly, a strong CH4 seep (flow rate  1200 mL min−1) was mapped  40 m north of site O246, making water advected from this direction highly elevated in CH4, with an average of  1400 nmol L−1 compared to the overall average of 577 nmol L−1 (Fig. 3b). The rapid changes in dissolved CH4 can to a high degree be explained by this relationship, due to the high variability in ocean current velocity. That this relationship holds for most of the measuring period also shows that even though observed average concentrations are lower in winter months, the seep configuration did not change significantly from July 2015 to May 2016, and dissolved CH4 was efficiently dispersed in relatively high concentrations in the whole seepage area.

Furthermore, daily CH4 concentrations at site O91 were higher on average than the 24 h measurements (313 vs. 200 nmol L−1). This can be explained by the comparable measurement periodicity (24 h) and tidal periodicity (23.93 h) in the ocean currents, resulting in predominantly eastward advection during daily measurements, thus systematically transferring water from a weak seepage area (Fig. 3). We did not observe this effect at site O246, most likely due to less tidal variance in the current direction (Fig. 2b). Nonetheless, this systematic tide-induced bias on the daily measurements at site O91 highlights the importance of taking the oceanographic conditions into account to avoid misinterpretation of variability.

Figure 3O91 (a) and O246 (b) location (yellow dot) and flow rates from flares mapped in their vicinity during CAGE 16-4 (color scale). Background color (green and blue) illustrates seafloor bathymetry. The compass diagram shows the relationship between ocean current direction (angle) and CH4 concentration (distance from center; black is response-time-corrected (RTC) data, and raw data are in blue).


Since currents are mostly northward and seepage is mostly located to the north of both observatories, averaged measured CH4 concentrations are likely lower than the average over the immediate surrounding area (Fig. 3). Despite this, the observatory data show higher average CH4 concentrations than previously reported. In the area surrounding site O91, Silyakova et al. (2020) reported average concentrations of 92, 70 and 61 nmol L−1 in June 2014, July 2015 and May 2016, respectively, based on discrete water sampling. Averaged CH4 concentrations measured at site O91 in July 2015 and May 2016 were 566 and 110 nmol L−1, respectively, i.e., around 8 and 2 times higher than values reported by Silyakova et al. (2020), respectively. The maximum CH4 concentration at O91 of 1748 ± 142 nmol L−1 on 20 August 2015 also significantly exceeds the previously maximum recorded concentration in the area of 480 nmol L−1 (July 2014, Silyakova et al., 2020). At site O246 the August 2016 average (564 nmol L−1) was 8 times higher than what Gentz et al. (2014) found in August 2010 (70 nmol L−1) using an altimeter-controlled CTD towed at 2 m above the seafloor. The maximum concentration in August 2016 also significantly exceeded previous observations, with 2661 ± 163 nmol L−1 compared to 524 nmol L−1 measured by Gentz et al. (2014).

These differences could be a result of temporal, local or regional differences in CH4 concentration. However, strong vertical gradients in dissolved CH4 are well documented at both seep sites (Gentz et al.2014), and our sensors measured closer to the seafloor (1.2 m above seafloor) compared to Gentz et al. (2014) (2 m above seafloor) and Silyakova et al. (2020) (5 to 15 m above seafloor). Additionally, the observatories were deployed close to seeps using a launcher, as opposed to “blind” water sampling from a shipborne rosette. Methane was also measured in situ, thereby avoiding potential CH4 outgassing after retrieval of water samples (Schlüter et al.1998).

Dissolved CH4 within shallow seep sites where gas can bypass the oceanic sinks often present heterogeneous distribution and rapid temporal variability (Gentz et al.2014; Myhre et al.2016a). Our results show that the temporal variability at the two seep sites is higher than previously reported and that changing ocean currents and configuration of nearby seeps are major contributors. This high short-term variability introduces a conceptual error in studies relying on discrete water sampling (e.g., to calculate inventories) because the time required to conduct the survey (measured in days) is much longer than large temporal variations in concentration (up to an order of 103 nmol L−1 within hours).

Figure 4Relative standard error of the mean for different numbers of samples N for O91 24 h data, data presented in Silyakova et al. (2020) (“June-14”, “July-15”, “May-16”), and O246 24 h data (in black). Relative standard deviation (corresponding to the standard error with N=1) is given in the legend (σrel). * Data from Silyakova et al. (2020) calculated assuming that the sample distribution resembles the underlying distribution (see Appendix D).


We can obtain a first-order constraint on errors caused by short-term variability in a hypothetical water sampling survey using the 24 h time series from the observatories. We assume the hypothetical survey seeks to find the average concentration in the bottom layer of the seep site. The expected error can then be found by calculating the standard error of the mean (SEM) for a given number of samples N using the 24 h time series as an underlying distribution representing the sub-daily variability of the seep site (Fig. 4; Appendix D contains a detailed outline of the methodology). Even though surveys often require more than 24 h to complete (2–3 d in Silyakova et al., 2020), a majority of processes causing short-term variability have periods below or at  24 h (for instance tides and many turbulent eddies; see, e.g., Sect. 3.2 and 3.1 and Talley et al.2011), likely making the daily distribution relevant also for surveys with longer duration. We compared SEM calculations based on the observatory 24 h time series with SEM calculations for the bottom water ( 5 m above the seafloor) discrete water sample data used for average and inventory estimates of the O91 seep site in Silyakova et al. (2020) (also included in Fig. 4).

The absolute SEM (in nmol L−1) is generally higher for time series with higher averaged concentrations, making the relative SEM cluster well, with gradually diminishing range for increasing N (an inherent property of the SEM, e.g., 12 %–45 % for N=10, 9 %–30 % for N=30, Fig. 4). The SEM of the data from Silyakova et al. (2020) is similar to the SEM of the 24 h time series, with a common range of 5 %–15 % expected error for surveys with N∼60 samples (N=64, 62, and 63 in Silyakova et al. 2020). It should be noted that the comparison with data from Silyakova et al. (2020) has caveats, e.g., that the observatory data do not contain errors due to spatial variability and an assumption of representative short-term temporal variability at the observatory sites (see also Appendix D).

Evidently, detailed surveys of individual seep sites, such as the study by Silyakova et al. (2020), can provide reasonable estimates of local inventories (< 15 % uncertainty) despite high short-term temporal variability. However, it is important to note that the area investigated in Silyakova et al. (2020) was densely mapped and homogeneous in the sense that it is an area where seepage is well documented (Silyakova et al.2020). Interpolation or averaging across larger regions where the amount of seepage is mostly unknown can result in considerable errors due to false interpolation assumptions and amplification of individual measurement errors, which can be large (expected errors up to  140 % for single measurements; see listed standard deviations in Fig. 4). These effects can potentially explain some of the discrepancies in estimates of oceanic CH4 inventories and fluxes.

Our findings stress the importance of sufficiently dense mapping and knowledge about the underlying seep condition when collecting water samples for inventory estimates. They also highlight the advantage of towed or autonomous instrumentation capable of providing continuous CH4 data, giving considerably better coverage and representation of the CH4 distribution in less time (e.g., Sommer et al.2015; Grilli et al.2018; Canning et al.2021). Assuming a distribution that better reflects the uneven spread of CH4 when applying interpolation and extrapolation techniques could also limit estimation errors. Future studies should investigate how initial errors due to short-term and small-scale variability propagate via different upscaling techniques and how these errors can be mitigated.

4.2 Hydrostatic pressure

Tidal changes in hydrostatic pressure can trigger CH4 release by build-up of CH4 in sediment pore water at rising tide and subsequent release when pore pressure decreases at falling tide, as observed at the Hikurangi Margin (Linke et al.2009) and Clayoquot Slope (Römer et al.2016). Our study sites differ from these sites in depth (they are >600 m) and in tidal amplitude (4 m at Clayoquot Slope compared to 1.5 offshore Prins Karls Forland). Linke et al. (2009) and Römer et al. (2016) also observed bubbles hydro-acoustically, while we measure dissolved CH4, which is strongly affected by the (also tidally dependent) current direction (Fig. 3).

To evaluate the effect that hydrostatic pressure changes have on the in situ concentration, we need to constrain the variance caused by changing current directions (since they operate in the same frequency domain). To do this, we first binned the CH4 concentration data into overlapping bins defined by the current direction at the time when the measurement was obtained and calculated standard scores (the number of standard deviations each value deviates from the sample mean; see, e.g., Kreyszig1979) for the data in each bin. We used larger current direction intervals for O246 due to the shorter data set, with a 12 window for O91 and a 30 window for O246. This resulted in a data set (i.e., the standard scores from all bins) effectively unrelated to the current direction. We then binned all the standard-scored CH4 data according to when the data were collected in relation to the M2 governed tidal cycle peak using overlapping 30 min bins (the M2 tide explains 79.2 % and 80.3 % of the pressure variance at O91 and O246, respectively). Average and median values were calculated for each bin, giving the averaged or median normalized dissolved CH4 value (standard score) for each current velocity defined data bin as a function of the M2 tidal cycle (Fig. 5). This partial decoupling of variability in hydrostatic pressure and current direction was possible since the bottom water current and hydrostatic pressure changes had different dominant tidal constituents; i.e., the current was mainly dominated by the diurnal K1 constituent ( 23.91 h period), while the M2 tide is semi-diurnal (12.42 h period).

A strong effect of the hydrostatic pressure on local seepage should elevate the standard scores at decreasing pressure (from 0 to 6.2 h, i.e., in the right half of Fig. 5), which we observe at both observatories. However, we observe stronger peaks at increasing hydrostatic pressure (3 h) at site O91 and at the M2 peak (0 h) at site O246, which contradicts this hypothesis. This does not mean that there is no effect of hydrostatic pressure changes but rather that the seepage in the area is widespread at both falling and rising tide conditions. The high variability caused by the strong effect of current direction also makes it particularly challenging to detect moderate changes in seepage intensity.

Figure 5Median and averaged standard scores of CH4 binned according to bottom water current direction and according to where the data were sampled in relation to the phase of the M2 pressure tide.


4.3 Bottom water temperature

Bottom water temperature can affect CH4 release by altering hydrate stability and CH4 solubility in pore water and the water column (Sloan1998; Jansson et al.2019a). Seasonal CH4 release variability resulting from temperature variations in the bottom water has been linked to migration of the gas hydrate stability zone (GHSZ) and hydrate dissociation further offshore at  390 m water depth (Berndt et al.2014; Ferré et al.2020). Our observatories were deployed in areas too shallow for gas hydrate to form. However, inversely varying seepage intensity between seepage at the GHSZ depth (390 m) and site O246 can suggest that these areas are fed by the same hydrocarbon source and that hydrates seasonally block the lateral pathways between these seep sites (Veloso-Alarcón et al.2019). This is in agreement with the observed long-term ( 3 months) negative correlation between bottom water temperature and dissolved CH4 at site O246 (R=-0.31). It should be noted that the same relationship is observed at O91; however, no geophysical data are available from this area due to the shallow depth.

Tidal pressure variations can affect CH4 release via pore water solubility (Sect. 4.2), but on longer timescales CH4 solubility is almost exclusively a function of water temperature. Higher CH4 solubility implies more CH4 dissolved in pore water and within bubble streams, potentially increasing the amount of CH4 dissolved in bottom water. A small but significant (R=0.33) positive correlation between CH4 solubility and concentration at site O246 and site O91 (considering the same time period, i.e., until 3 October in 2015) could indicate such an effect. This is also an alternative explanation for the negative correlation between temperature and CH4 concentration at site O246.

4.4 Pore water seepage

Short-term temperature increase further offshore (390 m depth) has been linked with release of warm, CH4-rich fluids from the sediments triggered by short-duration seismic events (Franek et al.2017). This means that increased CH4 concentration should be accompanied by increased water temperature and reduced salinity due to admixture of warmer, less saline pore water. We compared short-term anomalies (i.e., deviations from daily means) in these three variables in the 24 h data sets at both seep sites but found no corroborating evidence for this hypothesis. Instead, the covariance between current velocity and temperature and salinity anomalies indicates that short-term variability is mainly caused by cross-shelf exchange of Atlantic water in the West Spitsbergen Current and the colder, fresher Arctic water in the Coastal Current due to eddies (Hattermann et al.2016). It also indicates that CH4 release comes mainly from bubble dissolution and not from pore water seepage.

Figure 6(a) Bulk velocity shear (ΔH=8 m) and two-dimensional correlation with wind stress (contours). Relationships between 5 d accumulated wind stress and median (b) and median absolute deviation (c) of CH4 concentration for 24 h data periods. Persistent wind events with speeds of more than 10 m s−1 in periods over 6 h are indicated with blue stars along the x axis of panel (a). Blue highlights fall and winter water column conditions as described in the text.


4.5 Seasonal variation of CH4 distribution at site O91

Low release of CH4 to the atmosphere from the O91 seep area during summer despite high seabed influx has been explained by suppression of vertical mixing by strong stratification (Myhre et al.2016a) or absence of mechanical forcing such as wind stress (Silyakova et al.2020). However, in fall and winter, the water column offshore Prins Karls Forland is expected to have more horizontal and vertical mixing due to weaker stratification from cooling or sea ice formation (Tverberg et al.2014), baroclinic instability in the frontal structures of the West Spitsbergen Current (von Appen et al.2016; Hattermann et al.2016), and more frequent storms (Nilsen et al.2016).

We expect lower CH4 variability and lower CH4 concentration during periods of high mixing and dispersion due to weaker horizontal and vertical gradients and more efficient dispersion of CH4 away from sources. We use three sets of parameters to evaluate long-term changes in the amount of mixing in the water column (see Appendix E): (i) the 4-week averaged bulk velocity shear (Sb), (ii) the two-dimensional correlation between wind stress and current velocity (RWC), and (iii) the number of stormy days as defined by persistent winds > 11 ms−1 lasting longer than 6 h (Fig. 6). Calm weather and low Sb and RWC until mid-September 2015 indicate a stable water column with limited mixing in the bottom waters. From mid-September, Sb increased and stayed high until mid-November, together with a gradual increase in RWC that can be attributed to a gradual breakdown of stratification and increasing number of storm events (Fig. 6a). RWC remained high (RWC>0.5 at 60 m depth) until March 2016, indicating a significant effect of wind forcing in the water column. From March until observatory retrieval, RWC decreased to <0.2 below 50 m depth while Sb increased below 60 m depth, indicating available energy for mixing in the bottom waters.

We quantified CH4 variability during the 24 h measurements using the median absolute deviation (MAD) and used the median as a measure of the amount of dissolved CH4. The three 24 h periods collected during the calmer period prior to mid-September had high median concentration (> 300 nmol L−1) and the overall highest variability (MAD > 160 nmol L−1), as expected for low mixing conditions (Fig. 6b and c). From mid-September until the end of March (i.e., fall and winter seasons), the 24 h CH4 concentration time series had generally lower MAD and median concentration. In this period, CH4 variability and median also showed a good statistical relationship with the 5 d accumulated wind stress (R=-0.82 for MAD and R=-0.61 for median concentration), indicating that wind forcing has a deep impact on mixing and redistribution of CH4 in the water column (which also fits well with a high RWC). The two last 24 h CH4 time series (10 April and 1 May) had low median concentration, which could be explained by the absence of stratification (Silyakova et al.2020) and generation of mixing from the observed increase in Sb.

Accumulated wind stress, Sb and RWC are only limited indicators of water column dispersion and mixing. Nonetheless, the relationship between these parameters and the MAD and medians of the 24 h period CH4 time series gives a good indication of the seasonal cycle of distribution and vertical transport of CH4: strong stratification, less wind forcing and eddy activity in summer limit mixing and prevent CH4 from reaching the atmosphere. However, in fall and winter, reduced stratification makes the water column more prone to mixing, and distribution of CH4 seems to be strongly linked with wind forcing from September to April.

5 Conclusions

Time series of dissolved CH4 at both lander locations show considerably higher CH4 concentrations (up to 1748 ± 142 nmol L−1 at O91 and 2727 ± 182 nmol L−1 at O246) than previously found in ship-based water sampling surveys (maximum of 482 near O91 and of 564 near O246). The time series also uncover high CH4 variability (up to  1000 nmol L−1) within short timescales (<24 h), highlighting the potential uncertainty of flux and inventory estimates based on interpolation and extrapolation techniques relying on, e.g., an assumption of linearity. We calculated the standard error of a mean estimate based on a hypothetical discrete water sampling survey based on a range of samples by using the 24 h time series as the underlying distribution. The results aligned well with previous discrete water sampling surveys in the area, giving a standard error of the mean of 5 %–15 % for  60 samples.

Variability can be linked to directional ocean current variations occurring at tidal timescales, which shows the importance of taking the current direction and seep locations into account when interpreting intense seep site observations. The persistent relationship between current direction and location of seeps during recovery shows that there was seepage throughout the year and that the seep configuration was relatively constant.

We did not observe a direct effect of tidal pressure variations on CH4 release, but this could be hidden by the strong effect of variations in current direction. A negative (long-term) correlation between temperature and dissolved CH4 at O246 is in agreement with the hypothesized seasonal blocking of lateral CH4 pathways in the sediments (Veloso-Alarcón et al.2019) but could also be explained by increased CH4 solubility in the water column.

Short-term, small-scale variations in temperature and salinity were not linked with increased amounts of dissolved CH4, but rather with cross-frontal exchange of water masses due to eddies.

We observed a seasonal cycle in the characteristics of the 24 h time series that fits with seasonal changes in dispersion and mixing characteristics of the water column. Higher CH4 concentration and variability in early fall, when stratification was strong, was followed by lower median concentrations and variability in late fall and winter when the water column was more affected by mixing. In late fall and winter, wind forcing was statistically coupled to the concentration and variability of CH4, probably due to weaker water column stratification.

When estimating the atmospheric impact of a particular CH4 source based on sparse measurements, it is crucial to have some constraints on the temporal and spatial variability. These constraints can either be direct knowledge about variability itself or how inventory and fluxes are affected by related physical and/or chemical parameters. We observed considerable temporal and spatial variability at the two seep sites that need to be taken into account to obtain meaningful estimates of CH4 fluxes or inventories. That no strong direct link was found with other oceanographic parameters illustrates the nonlinearity of the system, making careful interpretation of measurements important. Future studies should aim to identify the errors that arise via different upscaling and interpolation techniques and how these errors can be mitigated. Based on our observations, we suggest that uncertainties in CH4 inventory and seep estimates can be mitigated by taking the local seep configuration, ocean currents and mixing rates into account and employing autonomous instrumentation capable of resolving the steep horizontal gradients in dissolved CH4. This, alongside direct measurements of seepage by, e.g., acoustic instrumentation, can help constrain future estimates of CH4 flux to the atmosphere from seabed seepage.

Appendix A: The K-Lander

Figure A1(a) The K-Lander is a 1.6 m high and 3.6 m wide trawl-proof stainless steel frame with multiple instrument mounts and batteries. The side panels are perforated to allow unobstructed water flow to the instruments inside the structure. See Appendix B for details on instrumentation. (b) One of the K-Landers during deployment with a launcher mounted on top and camera system mounted on a boom for visual control of landing area. (c) The two K-landers before deployment.


Appendix B: Measurement intervals, general post-processing and data

The CTD (oxygen sensor) and ADCP conducted measurements every 4 and 9 min, respectively, during the continuous monitoring of CH4 and CO2 measurements and every 21 and 29 min, respectively, during the rest of the deployment period (see Table B1 for acronyms, description and measurement accuracy). Salinity was measured on the practical salinity scale.

The upward-mounted ADCP measured ocean currents in 1 m bins with a bottom 7 m blank distance, where the topmost 20 % of the water column was disregarded due to side lobe interference. The high resolution, relatively short ensemble time (1 min) and potential presence of CH4 bubbles in the water resulted in noisy data. We dampened the noise by first removing any data points with error velocities exceeding one short-term (1 week) standard deviation, smoothed the data using a second-order Butterworth low-pass filter with a 3 h cutoff period and a spatial (i.e., vertical) moving average filter with a 5 m Hann window (increasing the blank distance to 10 m). The accuracy of the ADCP data is therefore not explicitly constrained and is based on comparing current velocity frequency spectra before and after filtering, combined with averaged error velocity of the raw data (Table B1).

Since sensors were recording at different frequencies, chronological alignment of the data was carried out by identifying nearest neighbor data points or by resampling. For correlation coefficients, histograms and Fourier analysis, the data sets were resampled to a uniform 15 min or 1 h measuring interval depending on the sample frequency of the raw data, using a poly-phase anti-aliasing filter. Due to the power-cycling mode of the CH4 and CO2 sensors and differing sampling frequencies, some statistics were based on more data points than others (outlined in Table B1). Daily measurements of CH4 were excluded from these statistics due to the high probability of systematic errors induced by periodic diurnal effects.

Harmonic analysis of hydrostatic pressure and ocean currents was done using t_tide (see Pawlowicz et al.2002) and the fast Fourier transform.

We calculated the rate of change (ROC) in CH4 concentration using the response-time-corrected CH4 data and the absolute value of the three-point (9 min) finite differences to limit the effect of noise on the calculation.

The absolute concentration of CH4 in the water (nmol L−1) was estimated from the partial pressure of CH4, pressure, temperature and salinity using Henry's law and Henry constants obtained from Harvey (1996) and the practical molar volume and gamma term from Duan and Mao (2006).

The CH4 sensors were calibrated to relevant water temperatures prior to deployment. The TDLAS detectors (Contros GmbH2018) provide measurements with good selectivity (fit for purpose) and high long-term stability (intermittent calibration not necessary) and are unaffected by dissolved oxygen content (unless via complete depletion). Biofouling was also minimal at retrieval (due to the cold water and local setting), and the PDMS membranes are almost unaffected by cold water. Generally, we did no observations indicating issues with any of the sensors except for what has already been mentioned regarding the conductivity probe and electrical malfunction of O246. Furthermore, we discarded all data recorded during instrument warm-up (i.e., when internal temperature was below correct operating temperature) before the individual measurement periods (the instruments were turned on  35 min prior to recording the data used in the analysis).

In the Contros HydroC CH4 and CO2 sensors, dissolved gases diffuse through a hydrophobic membrane into a gas chamber and equilibrate with the ambient environment. This results in the slow response time (e.g., τ6350 min under certain conditions for our membrane and pump setup for the CH4 sensor) and poor representation of the rapid changes in CH4 we expected in our study area (Gentz et al.2014; Myhre et al.2016a). We therefore performed a response time correction for the dissolved CH4 data following the methodology presented in Dølven et al. (2021), modulating the response time using the temperature data (effects of salinity on membrane permeability were not taken into account since these are negligible for the local ranges; see Robb1968). The CO2 sensors had a faulty pump, which ambiguously increased the response time of the sensors, making response time correction impossible.

The response time correction was performed for each period individually (1 and 24 h, i.e., 377 periods), using the stated measurement accuracy of the instrument (2 µatm or 3 % of measured value, whichever is higher) as input uncertainty. We first identified the ideal Δt according to the maximum curvature point in the L curves of the 24 h measurement periods. These varied slightly between each measurement period but averaged close to 180 s (176.4 s). To keep the same measuring interval for all the CH4 data, we therefore corrected all the data with a specified Δt of 180 s, which falls well within the bend of the L curve and should therefore safeguard a good balance between noise and model error (Fig. B1a). Inspection of model fit residuals showed a slight modulation following the variance in the signal, which is explained by our choice to use the same 3 min measurement grid across a relatively wide variance range, but the residuals were otherwise Gaussian. Although this is expected, it indicates that errors might be slightly overestimated for low-variance sections of the time series and vice versa for high-variance sections.

The uncertainty estimate varies depending on the amount of CH4 measured by the TDLAS unit in the measurement chamber of the instrument. The distribution of the uncertainty estimates is shown as percentages in Fig. B1b. Estimated uncertainty ranged from 3 to 205 nmol L−1 (95 % confidence, high for high concentrations in the measurement chamber and vice versa) or usually between 5 % and 20 %, although there were some outliers when the concentration was low and uncertainty estimate was high (Fig. B1b).

Table B1Instruments mounted on O91 and O246 (see Fig. A1), measured parameters, height in meters above seafloor (masf), and stated accuracy. ADCP stands for acoustic Doppler current profiler. N shows the number of data points used for later multi-variable analysis of O91 and O246.

* The Contros HydroC CH4 outputs partial pressure from the internal gas chamber. ** We report absolute concentration in seawater (nmol L−1) using Henry's law. *** We report accuracy only for response-time-corrected (RTC) concentrations (see Fig. B1) since the accuracy for untreated CH4 concentration data is ambiguous due to the slow response time. n/a: not applicable.

Download Print Version | Download XLSX

Figure B1(a) L curve for response time correction of CH4 data showing the location of the chosen Δt (180 s) for 6 May at O91(b) Estimated relative (%) uncertainty for response-time-corrected CH4 data (both observatories).


Appendix C: The 24 h measurements of CH4

Figure C1All 24 h periods of CH4 concentration at O91 using response-time-corrected data (black) with the uncertainty estimate (grey shading, 95 % confidence) and raw data (blue) from O91.


Figure C2All 24 h periods of CH4 concentration at O246 using response-time-corrected data (black) with the uncertainty estimate (grey shading, 95 % confidence) and raw data (blue) from O246.


Appendix D: Standard error of mean estimate due to temporal variability

To obtain the (theoretical) true dissolved CH4 average or inventory for an area requires a known concentration everywhere at a single point T0 in time. Considering a hypothetical ship-based discrete water sampling survey, any small-scale spatial variability not resolved by the sampling grid or localized (not seep-wide) short-term temporal variability occurring during the survey time can be considered a measurement errors for the purpose of the survey. Assuming that the water samples are sufficiently spaced out to be considered independent samples, the estimated average concentration from N samples in a particular depth layer at a seep site can be expressed as follows:

(D1) E m , ϵ t , ϵ s = n = 1 N m + ϵ t n + ϵ s n N ,

where m is the average of the seep site at T0, ϵt is error due to temporal short-term deviation from m at sampling time T0t and ϵs is spatial deviation in concentration from m. The expected standard error of E(m,ϵt,ϵs) from the short-term temporal and spatial variability is then given by

(D2) σ E m , ϵ t , ϵ s = σ N ,

where σ is the standard deviation of the distribution we sample from Ayyub and McCuen (2011). From Eqs. (D1) and (D2) we obtain

(D3) σ E m , ϵ t , ϵ s = σ E m , ϵ t + σ E m , ϵ s = σ t N + σ s N ,

where σt and σs are the distribution standard deviation related to temporal (ϵt) and spatial (ϵs) variability and σE(m,ϵt) and σE(m,ϵs) are the corresponding contributions to the standard error of the mean. Assuming the daily variance at the observatory is representative for the seep site, we can describe the expected error caused by sub-daily variability (all ϵt) in a scenario where a seep site is being sampled N times using the 24 h time series as the underlying distribution. In essence, we treat every measurement as having an associated probability distribution that is represented by the 24 h time series (which gives the sub-daily variability).

In the discrete water sample data presented in Silyakova et al. (2020), the underlying distribution is unknown, and we can only assume that the sample distribution resembles the underlying distribution, i.e., that

(D4) σ E m , ϵ t , ϵ s σ ^ E ( m , ϵ t , ϵ s ) = σ sampled N ,

where σ^E(m,ϵt,ϵs) is the standard error estimate of the mean based on the sample distribution and σsampled is the standard deviation of the measurements. All three data sets, i.e., “June-14” (N=64), “July-15” (N=62) and “May-16” (N=63), have a similarly skewed distribution compared to what is found in the observatory data (see Fig. D1), which supports this assumption. The survey in Silyakova et al. (2020) required 2–3 d to complete, while the observatory data only concern sub-daily variability (24 h time series). Nonetheless, we believe the comparison is valid, since the known major contributors to short-term (timescales shorter than weeks) variability act on sub-daily (or at least daily) scales, such as the dominant frequencies in the ocean currents and pressure changes.

There is a clear relationship of increasing σE(m,ϵt,ϵs) with increasing daily average, making relative σE(m,ϵt,ϵs) a meaningful quantity to use, as opposed to absolute σE(m,ϵt,ϵs). Additionally, for simplicity, we have not differentiated in the notation of the standard error of the mean (SEM) in the main text of the paper, referring to it as simply SEM in all situations.

It is also enlightening to consider the distribution of average estimates and how the skewed underlying distribution affects the distribution of average estimate errors for smaller N. We did this by simulating hypothetical surveys by random sampling from the 24 h data sets (Fig. D2), which shows the elevated probability of underestimating the average for estimates based on few samples (N30), i.e., that the median error is smaller than the average error. This is caused by an inheritance of the skewed underlying distribution in the CH4 concentration data (see Fig. D2a). This also allows for severe overestimates due to the long right-hand-side tail of the distribution. For larger values of N (N30), average estimates tend towards being normally distributed, thus avoiding these effects (see Fig. D2b).

Error estimates of more complicated properties, such as the total CH4 content in a volume of water based on interpolation techniques, require an assessment of the individual uncertainties of each measurement and how these errors propagate via, e.g., linear interpolation in the spatial domain. While not being explicitly applicable to inventory estimates, the σE still describes how random errors cancel out for larger values of N in evenly sampled grids, assuming this variability is representative for the seep site.

Figure D1Distribution of CH4 concentrations from the (a) June 2014, (b) July 2015 and (c) May 2016 data in Silyakova et al. (2020) and (d) from the 24 h data (all periods) at O91. Note the different scale for the y axis between (a)(c) and (d).


Figure D2Histograms of simulated average estimates based on N=10 (a) and N=30 (b) samples from the 24 h data set from 23 August at O91 showing the median and mean as vertical lines.


Appendix E: Bulk velocity shear and wind stress correlation

We calculated bulk wind stress using 10 m above sea level ERA-interim reanalysis wind data (Dee et al.2011; Large and Pond1981). Water column bulk velocity shear Sb (see, e.g., Lincoln et al.2016) was calculated as follows:

(E1) S b 2 = u u - u l h diff 2 + v u - v l h diff 2 ,

where uu, ul, vu and vl refer to the eastward and northward ADCP velocity components in the upper (subscript u) and lower (subscript l) layer and hdiff is the vertical distance between layers. The direct effect of wind stress is usually confined to surface water, although indirect effects such as Ekman transport or overturning and the formation of eddies can facilitate currents and mixing at deeper depths (Cushman-Roisin and Beckers2011). The two-dimensional correlation coefficient RWC between the wind and ocean currents was calculated using Kundu (1976) and the complex representations of the wind stress and de-tided current velocity vectors (τc and uc, respectively) as follows:

(E2) R WC = τ c * u c τ c * τ c 1 2 u c * u c 1 2 ,

where 〈…〉 gives the normalized inner product of the vectors and * annotates the complex conjugate. We allow time lags of up to 15 h to account for the gradual and indirect effects of wind stress on the ocean currents. Both properties were estimated throughout the valid current velocity profile, but this was only done only down to 80 m depth due to the 8 m vertical distance between the defined layers used in the bulk velocity shear calculation.

Code and data availability

All data presented in this paper can be obtained upon request to the authors and are also available in the platform Open research Data at the University of Tromsø – The Arctic University of Norway (, Dølven2022). All computer code used can be obtained upon request to the corresponding author.

Author contributions

The study was conceptualized by KOD, BF, AS, PL, and PJ. The data were curated by KOD, BF, and MM. Formal analysis was done by KOD, BF, and MM. Funding was acquired by BF. Investigation was done by KOD, BF, and AS. Methodology was implemented and developed by KOD and BF. The project was administrated by BF and AS. Resources was obtained by BF. Any software used was developed by KOD. The project was supervised by BF. Visualizations were made by KOD and MM. KOD wrote the manuscript (original draft), and KOD, BF, AS, PL, PJ, and MM contributed to reviewing and editing the manuscript.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We thank the crew of R/V Helmer Hanssen for their assistance during the deployment (CAGE 15-3) and recovery (CAGE 16-4) cruises. This study is a part of CAGE (Centre for Arctic Gas Hydrate, Environment and Climate), Norwegian Research Council (grant no. 223259). We thank Nicholas Warner for proofreading the article. We the thank two anonymous reviewers and the editor for their constructive feedback on the manuscript.

Financial support

This research has been supported by the Norwegian Research Council (grant no. 223259).

Review statement

This paper was edited by Mario Hoppema and reviewed by two anonymous referees.


Ayyub, B. M. and McCuen, R. H.: Probability, Statistics, and Reliability for Engineers and Scientists, Chapman & Hall/CRC, 3rd Edn., CRC Press, p. 409, ISBN 9781439809518, 2011. a

Berndt, C., Feseker, T., Treude, T., Krastel, S., Liebetrau, V., Niemann, H., Bertics, V. J., Dumke, I., Dünnbier, K., Ferré, B., Graves, C., Gross, F., Hissmann, K., Hühnerbach, V., Krause, S., Lieser, K., Schauer, J., and Steinle, L.: Temporal Constraints on Hydrate-Controlled Methane Seepage off Svalbard, Science, 343, 284–287,, 2014. a

Braga, R., Iglesias, R., Romio, C., Praeg, D., Miller, D., Viana, A., and Ketzer, J.: Modelling methane hydrate stability changes and gas release due to seasonal oscillations in bottom water temperatures on the Rio Grande cone, offshore southern Brazil, Mar. Petrol. Geol., 112, 104071,, 2020. a

Canning, A., Fietzek, P., Rehder, G., and Körtzinger, A.: Technical note: Seamless gas measurements across the land–ocean aquatic continuum – corrections and evaluation of sensor data for CO2, CH4 and O2 from field deployments in contrasting environments, Biogeosciences, 18, 1351–1373,, 2021. a

Contros GmbH: CONTROS HydroC™ CH4 Sensor for dissolved methane, available at: (last access: 5 January 2022), 2018. a

Cottier, F., Nilsen, F., Inall, M. E., Gerland, S., Tverberg, V., and Svendsen, H.: Wintertime warming of an Arctic shelf in response to large-scale atmospheric circulation, Geophys. Res. Lett., 34, L10607,, 2007. a

Cushman-Roisin, B. and Beckers, J.-M.: Introduction to Geophysical Fluid Dynamics, Elsevier Academic Press, 2nd Edn., ISBN 9780120887590, 2011. a

Dee, D., Uppala, S., Simmons, A., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hersbach, H., Hólm, E., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A., Monge-Sanz, B., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597,, 2011. a

Dølven, K. O.: Replication data for Autonomous methane seep site monitoring offshore Western Svalbard: Hourly to seasonal variability and associated oceanographic parameters, V1, DataverseNO [data set],, 2022. a

Dølven, K. O., Vierinen, J., Grilli, R., Triest, J., and Ferré, B.: Response time correction of slow response sensor data by deconvolution of the growth-law equation, Geosci. Instrum. Method. Data Syst. Discuss. [preprint],, in review, 2021. a, b

Duan, Z. and Mao, S.: A thermodynamic model for calculating methane solubility, density and gas phase composition of methane-bearing aqueous fluids from 273 to 523 K and from 1 to 2000 bar, Geochim. Cosmochim. Ac., 70, 3369–3386,, 2006. a

Etiope, G., Ciotoli, G., Schwietzke, S., and Schoell, M.: Gridded maps of geological methane emissions and their isotopic signature, Earth Syst. Sci. Data, 11, 1–22,, 2019. a

Ferré, B., Mienert, J., and Feseker, T.: Ocean temperature variability for the past 60 years on the Norwegian-Svalbard margin influences gas hydrate stability on human time scales, J. Geophys. Res.-Ocean., 117, C10017,, 2012. a

Ferré, B., Jansson, P., Moser, M., Portnov, A., Graves, C., Panieri, G., Gründger, F., Berndt, C., Lehmann, M., and Niemann, H.: Reduced methane seepage from Arctic sediments during cold bottom-water conditions, Nat. Geosci., 13, 144–148,, 2020. a

Franek, P., Plaza-Faverola, A., Mienert, J., Buenz, S., Ferré, B., and Hubbard, A.: Microseismicity Linked to Gas Migration and Leakage on the Western Svalbard Shelf, Geochem. Geophy. Geosy., 18, 4623–4645,, 2017. a, b

Gentz, T., Damm, E., von Deimling, J. S., Mau, S., McGinnis, D. F., and Schlüter, M.: A water column study of methane around gas flares located at the West Spitsbergen continental margin, Cont. Shelf Res., 72, 107–118,, 2014. a, b, c, d, e, f, g

Gerkema, T.: Tidal Constituents and the Harmonic Method, in: Introduction to Tides, Cambridge University Press, 1st Edn., 60–86, ISBN 9781108474269,, 2019. a, b

Graves, C. A., Lea, S., Gregor, R., Niemann, H., Connely, D. P., Lowry, D., Fisher, R. E., Stott, A. W., Sahling, H., and James, R. H.: Fluxes and fate of dissolved methane released at the seafloor at the landward limit of the gas hydrate stability zone offshore western Svalbard, J. Geophys. Res.-Ocean., 120, 6185–6201,, 2015. a

Grilli, R., Triest, J., Chappellaz, J., Calzas, M., Desbois, T., Jansson, P., Guillerm, C., Ferré, B., Lechevallier, L., Ledoux, V., and Romanini, D.: Sub-Ocean: Subsea Dissolved Methane Measurements Using an Embedded Laser Spectrometer Technology, Environ. Sci. Technol., 52, 10543–10551,, 2018. a

Hanson, R. S. and Hanson, T. E.: Methanotrophic bacteria, Microbiol. Rev., 60, 439–471,, 1996. a

Harvey, A. H.: Semiempirical correlation for Henry's constants over large temperature ranges, AIChE J., 42, 1491–1494,, 1996. a

Hattermann, T., Erik, I. P., Wilken Jon, A., Jon, A., and Arild, S.: Eddy-driven recirculation of Atlantic Water in Fram Strait, Geophys. Res. Lett., 43, 3406–3414,, 2016. 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,, 2016. a, b

Jansson, P., Ferré, B., Silyakova, A., Dølven, K. O., and Omstedt, A.: A new numerical model for understanding free and dissolved gas progression toward the atmosphere in aquatic methane seepage systems, Limnol. Oceanogr.-Method., 17, 223–239,, 2019a. a, b

Jansson, P., Triest, J., Grilli, R., Ferré, B., Silyakova, A., Mienert, J., and Chappellaz, J.: High-resolution underwater laser spectrometer sensing provides new insights into methane distribution at an Arctic seepage site, Ocean Sci., 15, 1055–1069,, 2019b. a

Kossel, E., Bigalke, N., Piñero, E., and Haeckel, M.: The SUGAR Toolbox, PANGAEA,, 2013. a

Kreyszig, E.: Advanced Engineering Mathematics, Wiley, 4 Edn., John Wiley and Sons Ltd., ISBN 9780471042716, 1979. a

Kundu, P. K.: Ekman Veering Observed near the Ocean Bottom, J. Phys. Oceanogr., 6, 238–242,<0238:EVONTO>2.0.CO;2, 1976. a

Large, W. G. and Pond, S.: Open Ocean Momentum Flux Measurements in Moderate to Strong Winds, J. Phys. Oceanogr., 11, 324–336,<0324:OOMFMI>2.0.CO;2, 1981. a

Lincoln, B. J., Rippeth, T. P., and Simpson, J. H.: Surface mixed layer deepening through wind shear alignment in a seasonally stratified shallow sea, J. Geophys. Res.-Ocean., 121, 6021–6034,, 2016. a

Linke, P., Sommer, S., Rovelli, L., and McGinnis, D. F.: Physical limitations of dissolved methane fluxes: The role of bottom-boundary layer processes, Mar. Geol., 272, 209–222,, 2009. a, b, c

Loeng, H.: Features of the physical oceanographic conditions of the Barents Sea, Polar Res., 10, 5–18,, 1991. a

Mau, S., Romer, M., Torres, M. E., Bussmann, I., Pape, T., Damm, E., Geprags, P., Wintersteller, P., Hsu, C.-W., Loher, M., and Bohrmann, G.: Widespread methane seepage along the continental margin off Svalbard – from Bjørnøya to Kongsfjorden, Sci. Rep., 7, 42997,, 2017. a

McDougall, T. J. and Barker, P. M.: Getting started with TEOS-10 and the Gibbs Seawater (GSW) Oceanographic Toolbox, SCOR/IAPSO WG127, 22 pp., ISBN 9780646556215, 2011. a

McGinnis, D. F., Greinert, J., Artemov, Y., Beaubien, S. E., and Wüest, A.: Fate of rising methane bubbles in stratified waters: How much methane reaches the atmosphere?, J. Geophys. Res.-Ocean., 111, C09007,, 2006. a, b

Myhre, C. L., Ferré, B., Platt, S. M., Silyakova, A., Hermansen, O., Allen, G., Pisso, I., Schmidbauer, N., Stohl, A., Pitt, J., Jansson, P., Greinert, J., Percival, C., Fjaeraa, A. M., O'Shea, S. J., Gallagher, M., Le Breton, M., Bower, K. N., Bauguitte, S. J. B., Dalsøren, S., Vadakkepuliyambatta, S., Fisher, R. E., Nisbet, E. G., Lowry, D., Myhre, G., Pyle, J. A., Cain, M., and Mienert, J.: Extensive release of methane from Arctic seabed west of Svalbard during summer 2014 does not influence the atmosphere, Geophys. Res. Lett., 43, 4624–4631,, 2016a. a, b, c, d, e

Myhre, C. L., Hermansen, O., Fiebig, M., Lunder, C., Fjæraa, A. M., Svendby, T., Platt, M., Hansen, G., Scmidbauer, N., and T., K.: Monitoring of greenhouse gases and aerosols at Svalbard and Birkenes in 2015 – Annual report, Norwegian Institute for Air Research (NILU), NILU report, 31/2016, 2016b. a

Nilsen, F., Skogseth, R., Vaardal-Lunde, J., and Inall, M.: A Simple Shelf Circulation Model: Intrusion of Atlantic Water on the West Spitsbergen Shelf, J. Phys. Oceanogr., 46, 1209–1230,, 2016. a, b, c

Pachauri, R. K. and Meyer, L. A. (Eds.): IPCC, 2014: Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, IPCC, Geneva, Switzerland, 151 pp., 2014. a

Pawlowicz, R., B., B., and Lentz, S.: Classical Tidal Harmonic Analysis Including Error Estimates in MATLAB using ttide, Comput. Geosci., 28, 929–937, 2002. a

Platt, S. M., Eckhardt, S., Ferré, B., Fisher, R. E., Hermansen, O., Jansson, P., Lowry, D., Nisbet, E. G., Pisso, I., Schmidbauer, N., Silyakova, A., Stohl, A., Svendby, T. M., Vadakkepuliyambatta, S., Mienert, J., and Lund Myhre, C.: Methane at Svalbard and over the European Arctic Ocean, Atmos. Chem. Phys., 18, 17207–17224,, 2018. a

Portnov, A., Vadakkepuliyambatta, S., Mienert, J., and Hubbard, A.: Ice-sheet-driven methane storage and release in the Arctic, Nat. Commun., 7, 10314,, 2016. a

Rajan, A., Mienert, J., and Bünz, S.: Acoustic evidence for a gas migration and release system in Arctic glaciated continental margins offshore NW-Svalbard, Mar. Petrol. Geol., 32, 36–49,, 2012. a

Reagan, M. T., Moridis, G. J., Elliott, S. M., and Maltrud, M.: Contribution of oceanic gas hydrate dissociation to the formation of Arctic Ocean methane plumes, J. Geophys. Res.-Ocean., 116, C09014,, 2011. a

Reeburgh, W. S.: Oceanic Methane Biogeochemistry, Chem. Rev., 107, 486–513,, 2007. a

Robb, W. L.: Thin silicone membranes – Their permeation properties and some applications, Ann. NY Acad. Sci., 146, 119–137,, 1968. a

Römer, M., Riedel, M., Scherwath, M., Heesemann, M., and Spence, G. D.: Tidally controlled gas bubble emissions: A comprehensive study using long-term monitoring data from the NEPTUNE cabled observatory offshore Vancouver Island, Geochem. Geophy. Geosy., 17, 3797–3814,, 2016. a, b, c

Ruppel, C. and Kessler, J.: The interaction of climate change and methane hydrates, Rev. Geophys., 55, 126–168,, 2017. a

Sahling, H., Römer, M., Pape, T., Bergès, B., dos Santos Fereirra, C., Boelmann, J., Geprägs, P., Tomczyk, M., Nowald, N., Dimmler, W., Schroedter, L., Glockzin, M., and Bohrmann, G.: Gas emissions at the continental margin west of Svalbard: mapping, sampling, and quantification, Biogeosciences, 11, 6029–6046,, 2014. a, b, c

Saloranta, T. M. and Svendsen, H.: Across the Arctic front west of Spitsbergen: high-resolution CTD sections from 1998–2000, Polar Res., 20, 177–184, 2001. a, b

Sarkar, S., Berndt, C., Minshull, T. A., Westbrook, G. K., Klaeschen, D., Masson, D. G., Chabert, A., and Thatcher, K. E.: Seismic evidence for shallow gas-escape features associated with a retreating gas hydrate zone offshore west Svalbard, J. Geophys. Res.-Sol. Ea., 117, B09102,, 2012. a

Saunois, M., Jackson, R. B., Bousquet, P., Poulter, B., and Canadell, J. G.: The growing role of methane in anthropogenic climate change, Environ. Res. Lett., 11, 120207,, 2016. a, b

Saunois, M., R. Stavert, A., Poulter, B., Bousquet, P., G. Canadell, J., B. Jackson, R., A. Raymond, P., J. Dlugokencky, E., Houweling, S., K. Patra, P., Ciais, P., K. Arora, V., Bastviken, D., Bergamaschi, P., R. Blake, D., Brailsford, G., Bruhwiler, L., M. Carlson, K., Carrol, M., Castaldi, S., Chandra, N., Crevoisier, C., M. Crill, P., Covey, K., L. Curry, C., Etiope, G., Frankenberg, C., Gedney, N., I. Hegglin, M., Höglund-Isaksson, L., Hugelius, G., Ishizawa, M., Ito, A., Janssens-Maenhout, G., M. Jensen, K., Joos, F., Kleinen, T., B. Krummel, P., L. Langenfelds, R., G. Laruelle, G., Liu, L., MacHida, T., Maksyutov, S., C. McDonald, K., McNorton, J., A. Miller, P., R. Melton, J., Morino, I., Müller, J., Murguia-Flores, F., Naik, V., Niwa, Y., Noce, S., O'Doherty, S., J. Parker, R., Peng, C., Peng, S., P. Peters, G., Prigent, C., Prinn, R., Ramonet, M., Regnier, P., J. Riley, W., A. Rosentreter, J., Segers, A., J. Simpson, I., Shi, H., J. Smith, S., Paul Steele, L., F. Thornton, B., Tian, H., Tohjima, Y., N. Tubiello, F., Tsuruta, A., Viovy, N., Voulgarakis, A., S. Weber, T., Van Weele, M., R. Van Der Werf, G., F. Weiss, R., Worthy, D., Wunch, D., Yin, Y., Yoshida, Y., Zhang, W., Zhang, Z., Zhao, Y., Zheng, B., Zhu, Q., Zhu, Q., and Zhuang, Q.: The global methane budget 2000–2017, Earth Syst. Sci. Data, 12, 1561–1623,, 2020. a, b, c

Schlüter, M., Linke, P., and Suess, E.: Geochemistry of a sealed deep-sea borehole on the Cascadia Margin, Mar. Geol., 148, 9–20,, 1998. a

Shakhova, N., Semiletov, I., Leifer, I., Salyuk, A., Rekant, P., and Kosmach, D.: Geochemical and geophysical evidence of methane release over the East Siberian Arctic Shelf, J. Geophys. Res.-Ocean., 115, C08007,, 2010. a

Silyakova, A., Jansson, P., Serov, P., Ferré, B., Pavlov, A. K., Hattermann, T., Graves, C. A., Platt, S. M., Myhre, C. L., Gründger, F., and Niemann, H.: Physical controls of dynamics of methane venting from a shallow seep area west of Svalbard, Cont. Shelf Res., 194, 104030,, 2020. a, b, c, d, e, f, g

Sloan, E. D.: Physical/chemical properties of gas hydrates and application to world margin stability and climatic change, Geol. Soc. Lond. Sp. Publ., 137, 31–50, 1998. a, b

Sommer, S., Schmidt, M., and Linke, P.: Continuous inline mapping of a dissolved methane plume at a blowout site in the Central North Sea UK using a membrane inlet mass spectrometer – Water column stratification impedes immediate methane release into the atmosphere, Mar. Petrol. Geol., 68, 766–775,, 2015. a

Steinle, L., Graves, C., Treude, T., Ferre, B., Biastoch, A., Bussmann, I., Berndt, C., Krastel, S., James, R., Behrens, E., Böning, C., Greinert, J., Sapart, C., Scheinert, M., Sommer, S., Lehmann, M., and Niemann, H.: Water column methanotrophy controlled by a rapid oceanographic switch, Nat. Geosci., 8, 378–382,, 2015. a

Swift, J. H. and Aagaard, K.: Seasonal transitions and water mass formation in the Iceland and Greenland seas, Deep-Sea Res. Pt. A., 28, 1107–1129,, 1981.  a

Talley, L. D., Pickard, G. L., Emery, W. J., and Swift, J. H.: Chapter 1 – Introduction to Descriptive Physical Oceanography, in: Descriptive Physical Oceanography, 6th Edn., edited by: Talley, L. D., Pickard, G. L., Emery, W. J., and Swift, J. H., Academic Press, Boston, 1–6,, 2011. a

Tverberg, V., Nøst, O. A., Lydersen, C., and Kovacs, K. M.: Winter sea ice melting in the Atlantic Water subduction area, Svalbard Norway, J. Geophys. Res.-Ocean., 119, 5945–5967,, 2014. a

Veloso, M., Greinert, J., Mienert, J., and Batist, M.: A new methodology for quantifying bubble flow rates in deep water using splitbeam echosounders: Examples from the Arctic offshore NW-Svalbard, Limnol. Oceanogr.-Method., 13, 267–287, 2015. a

Veloso-Alarcón, M. E., Jansson, P., Batist, M. D., Minshull, T. A., Westbrook, G. K., Pälike, H., Bünz, S., Wright, I., and Greinert, J.: Variability of Acoustically Evidenced Methane Bubble Emissions Offshore Western Svalbard, Geophys. Res. Lett., 46, 9072–9081,, 2019. a, b, c, d

von Appen, W.-J., Schauer, U., Hattermann, T., and Beszczynska-Möller, A.: Seasonal Cycle of Mesoscale Instability of the West Spitsbergen Current, J. Phys. Oceanogr., 46, 1231–1254,, 2016. a

Westbrook, G. K., Thatcher, K. E., Rohling, E. J., Piotrowski, A. M., Pälike, H., Osborne, A. H., Nisbet, E. G., Minshull, T. A., Lanoisellé, M., James, R. H., Hühnerbach, V., Green, D., Fisher, R. E., Crocker, A. J., Chabert, A., Bolton, C., Beszczynska-Möller, A., Berndt, C., and Aquilina, A.: Escape of methane gas from the seabed along the West Spitsbergen continental margin, Geophys. Res. Lett., 36, L15608,, 2009. a, b

Short summary
Natural sources of atmospheric methane need to be better described and quantified. We present time series from ocean observatories monitoring two seabed methane seep sites in the Arctic. Methane concentration varied considerably on short timescales and seasonal scales. Seeps persisted throughout the year, with increased potential for atmospheric release in winter due to water mixing. The results highlight and constrain uncertainties in current methane estimates from seabed methane seepage.