Description: We present a comprehensive study showing new results from a shallow gas seep area in approximate to 40 m water depth located in the North Sea, Netherlands sector B13 that we call Dutch Dogger Bank seep area. It has been postulated that methane presumably originating from a gas reservoir in approximate to 600 m depth below the seafloor is naturally leaking to the seafloor. Our ship-based subbottom echosounder data indicate that the migrating gas is trapped in numerous gas pockets in the shallow sediments. The gas pockets are located at the boundary between the top of the Late Pliocene section and overlying fine-grained sediments, which were deposited during the early Holocene marine transgression after the last glaciation. We mapped gas emissions during three R/V Heincke cruises in 2014, 2015, and 2016 and repeatedly observed up to 850 flares in the study area. Most of them (approximate to 80%) were concentrated at five flare clusters. Our repeated analysis revealed spatial similarities of seep clusters, but also heterogeneities in emission intensities. A first calculation of the methane released from these clusters into the water column revealed a flow rate of 277 L/min (SD=140), with two clusters emitting 132 and 142 L/min representing the most significant seepage sites. Above these two flare clusters, elevated methane concentrations were recorded in atmospheric measurements. Our results illustrate the effective transport of methane via gas bubbles through a approximate to 40 m water column, and furthermore provide an estimate of the emission rate needed to allow for a contribution to the atmospheric methane concentration.

Description: The Vestnesa Ridge comprises a 〉100 km long sediment drift located between the western continental slope of Svalbard and the Arctic mid-ocean ridges. It hosts a deep water (〉1000 m) gas hydrate and associated seafloor seepage system. Near-seafloor headspace gas compositions and its methane carbon isotopic signature along the ridge indicate a predominance of thermogenic gas sources feeding the system. Prediction of the base of the gas hydrate stability zone for theoretical pressure and temperature conditions and measured gas compositions results in an unusual underestimation of the observed bottom-simulating reflector (BSR) depth. The BSR is up to 60 m deeper than predicted for pure methane and measured gas compositions with 〉99% methane. Models for measured gas compositions with 〉4% higher-order hydrocarbons result in a better BSR approximation. However, the BSR remains 〉20 m deeper than predicted in a region without active seepage. A BSR deeper than predicted is primarily explained by unaccounted spatial variations in the geothermal gradient and by larger amounts of thermogenic gas at the base of the gas hydrate stability zone. Hydrates containing higher-order hydrocarbons form at greater depths and higher temperatures and contribute with larger amounts of carbons than pure methane hydrates. In thermogenic provinces, this may imply a significant upward revision (up to 50% in the case of Vestnesa Ridge) of the amount of carbon in gas hydrates.

Description: Estimating the spatial position of organisms is essential to quantify interactions between the organism and the characteristics of its surroundings, for example, predator–prey interactions, habitat selection, and social associations. Because marine mammals spend most of their time under water and may appear at the surface only briefly, determining their exact geographic location can be challenging. Here, we developed a photogrammetric method to accurately estimate the spatial position of marine mammals or birds at the sea surface. Digital recordings containing landscape features with known geographic coordinates can be used to estimate the distance and bearing of each sighting relative to the observation point. The method can correct for frame rotation, estimates pixel size based on the reference points, and can be applied to scenarios with and without a visible horizon. A set of R functions was written to process the images and obtain accurate geographic coordinates for each sighting. The method is applied to estimate the spatiotemporal fine-scale distribution of harbour porpoises in a tidal inlet. Video recordings of harbour porpoises were made from land, using a standard digital single-lens reflex (DSLR) camera, positioned at a height of 9.59 m above mean sea level. Porpoises were detected up to a distance of ~3136 m (mean 596 m), with a mean location error of 12 m. The method presented here allows for multiple detections of different individuals within a single video frame and for tracking movements of individuals based on repeated sightings. In comparison with traditional methods, this method only requires a digital camera to provide accurate location estimates. It especially has great potential in regions with ample data on local (a)biotic conditions, to help resolve functional mechanisms underlying habitat selection and other behaviors in marine mammals in coastal areas.

Description: Large internal wave breaking is observed exceeding a vertical array of 61 high-resolution temperature sensors at 1 m intervals between 7 and 67 m above the bottom. The array was moored for 5 days at 969 m of Opouawe Bank, New Zealand, a known methane seep area. As breaking internal waves dominate sediment resuspension above sloping topography in other ocean areas, they are expected to also influence methane transport. Despite being visible in single beam echosounder data, indications for turbulence due to rising gas bubbles are not found in the present 1 Hz sampled temperature records. Likely, the mooring was too far away from the very localized bubble release spot. Instead, the temperature sensors show detailed internal wave-turbulence transitions. Every tidal cycle, a solibore (a frontal turbulent bore with a train of trailing solitary waves) changes shape and intensity. These solibores are highly turbulent and they restratify the bottom boundary layer, thereby maintaining efficient mixing. Details of different turbulent bore developments are discussed. Averaged over a few tidal cycles and over the sensors range, mean vertical eddy diffusivity amounts to 3 ± 1 × 10−3 m2 s−1 and mean turbulent kinetic energy dissipation to 1.6 ± 0.7 × 10−7 W kg−1, with variations over 4 orders of magnitude. Such turbulence will affect the distribution of dissolved methane and other geochemical species in the lower 100–150 m above the bottom and their release from the bottom. The above mean values are remarkably similar to those found at various other sites in the NE Atlantic Ocean.

Description: We use new gas-hydrate geochemistry analyses, echosounder data, and three-dimensional P-Cable seismic data to study a gas-hydrate and free-gas system in 1200 m water depth at the Vestnesa Ridge offshore NW Svalbard. Geochemical measurements of gas from hydrates collected at the ridge revealed a thermogenic source. The presence of thermogenic gas and temperatures of similar to 3.3 degrees C result in a shallow top of the hydrate stability zone (THSZ) at similar to 340 m below sea level (mbsl). Therefore, hydrate-skinned gas bubbles, which inhibit gas-dissolution processes, are thermodynamically stable to this shallow water depth. This was confirmed by hydroacoustic observations of flares in 2010 and 2012 reaching water depths between 210 and 480 mbsl. At the seafloor, bubbles are released from acoustically transparent zones in the seismic data, which we interpret as regions where free gas is migrating through the hydrate stability zone (HSZ). These intrusions result in vertical variations in the base of the HSZ (BHSZ) of up to similar to 150 m, possibly making the shallow hydrate reservoir more susceptible to warming. Such Arctic gas-hydrate and free-gas systems are important because of their potential role in climate change and in fueling marine life, but remain largely understudied due to limited data coverage in seasonally ice-covered Arctic environments.

Description: We find that summer methane (CH4) release from seabed sediments west of Svalbard substantially increases CH4 concentrations in the ocean but has limited influence on the atmospheric CH4 levels. Our conclusion stems from complementary measurements at the seafloor, in the ocean, and in the atmosphere from land-based, ship and aircraft platforms during a summer campaign in 2014. We detected high concentrations of dissolved CH4 in the ocean above the seafloor with a sharp decrease above the pycnocline. Model approaches taking potential CH4 emissions from both dissolved and bubble-released CH4 from a larger region into account reveal a maximum flux compatible with the observed atmospheric CH4 mixing ratios of 2.4–3.8 nmol m−2 s−1. This is too low to have an impact on the atmospheric summer CH4 budget in the year 2014. Long-term ocean observatories may shed light on the complex variations of Arctic CH4 cycles throughout the year.

Description: Takahe seep, located on the Opouawe Bank, Hikurangi Margin, is characterized by a well-defined subsurface seismic chimney structure ca. 80,500 m2 in area. Sub-seafloor geophysical data based on acoustic anomaly layers indicated the presence of gas hydrate and free gas layers within the chimney structure. Reaction-transport modeling was applied to porewater data from 11 gravity cores to constrain methane turnover rates and benthic methane fluxes in the upper 10 m. Model results show that methane dynamics were highly variable due to transport and dissolution of ascending gas. The dissolution of gas (up to 3761 mmol m−2 yr−1) dwarfed the rate of methanogenesis within the simulated sediment column (2.6 mmol m−2 yr−1). Dissolved methane is mainly consumed by anaerobic oxidation of methane (AOM) at the base of the sulfate reduction zone and trapped by methane hydrate formation below it, with maximum rates in the central part of the chimney (946 and 2420 mmol m−2 yr−1, respectively). A seep-wide methane budget was constrained by combining the biogeochemical model results with geophysical data and led to estimates of AOM rates, gas hydrate formation and benthic dissolved methane fluxes of 3.68 × 104 mol yr−1, 73.85 × 104 mol yr−1and 1.19 × 104 mol yr−1, respectively. A much larger flux of methane probably escapes in gaseous form through focused bubble vents. The approach of linking geochemical model results with spatial geophysical data put forward here can be applied elsewhere to improve benthic methane turnover rates from limited single spot measurements to larger spatial scales.

Description: Quantifying marine methane fluxes of free gas (bubbles) from the seafloor into the water column is of importance for climate related studies, for example, in the Arctic, reliable methodologies are also of interest for studying man-made gas and oil leakage systems at hydrocarbon production sites. Hydroacoustic surveys with singlebeam and nowadays also multibeam systems have been proven to be a successful approach to detect bubble release from the seabed. A number of publications used singlebeam echosounder data to indirectly quantify free gas fluxes via empirical correlations between gas fluxes observed at the seafloor and the hydroacoustic response. Others utilize the hydroacoustic information in an inverse modeling approach to derive bubble fluxes. Here, we present an advanced methodology using data from splitbeam echosounder systems for analyzing gas release water depth (〉 100m). We introduce a new MATLAB-based software for processing and interactively editing data and we present how bubble-size distribution, bubble rising speed and the model used for calculating the backscatter response of single bubbles influence the final gas flow rate calculations. As a result, we highlight the need for further investigations on how large, wobbly bubbles, bubble clouds, and multi-scattering influence target strength. The results emphasize that detailed studies of bubble-size distributions and rising speeds need to be performed in parallel to hydroacoustic surveys to achieve realistic mediated methane flow rate and flux quantifications.

Description: Large reservoirs of methane present in Arctic marine sediments are susceptible to rapid warming, promoting increasing methane emissions. Gas bubbles in the water column can be detected, and flow rates can be quantified using hydroacoustic survey methods, making it possible to monitor spatiotemporal variability. We present methane (CH4) bubble flow rates derived from hydroacoustic data sets acquired during 11 research expeditions to the western Svalbard continental margin (2008-2014). Three seepage areas emit in total 725-1,125 t CH4/year, and bubble fluxes are up to 2 kg.m(-2).year (-1). Bubble fluxes vary between different surveys, but no clear trend can be identified. Flux variability analyses suggest that two areas are geologically interconnected, displaying alternating flow changes. Spatial migration of bubble seepage was observed to follow seasonal changes in the theoretical landward limit of the hydrate stability zone, suggesting that formation/dissociation of shallow hydrates, modulated by bottom water temperatures, influences seafloor bubble release. Plain Language Summary It has been speculated that the release of methane (a potent greenhouse gas) from the seafloor in some Arctic Ocean regions is triggered by warming seawater. Emissions of gas bubbles from the seafloor can be detected by ship-mounted sonars. In 2008, a methane seepage area west of Svalbard was hydroacoustically detected for the first time. This seepage was hypothesized to be caused by dissociation of hydrates (ice-like crystals consisting of methane and water) due to ocean warming. We present an analysis of sonar data from 11 surveys conducted between 2008 and 2014. This study is the first comparison of methane seepage-related hydroacoustic data over such a long period. The hydroacoustic mapping and quantification method allowed us to assess the locations and intensity of gas bubble release, and how these parameters change over time, providing necessary data for numerical flux and climate models. No trend of increasing gas flow was identified. However, we observed seasonal variations potentially controlled by seasonal formation and dissociation of shallow hydrates. The hydrate formation/dissociation process is likely controlled by changes of bottom water temperatures. Alternating gas emissions between two neighboring areas indicate the existence of fluid pathway networks within the sediments.

Description: Water column imaging multibeam echo sounder systems (MBESs) are a promising technology for quantitative estimates of the gas bubble volume flow within large gas seepage areas. Considerable progress has been made in recent years toward applicable calibration methods for MBESs as well as developing inversion models to convert acoustically measured backscattering cross sections to gas bubble volume flow. However, MBESs are still not commonly used for quantitative gas flow assessments. A reason for this is the absence of published processing methods that demonstrate how MBES data can be processed to quantitatively represent bubble streams. Here, we present a novel method (echo grid integration) that allows for assessing the aggregated backscattering cross section of targets within horizontal water layers. This derived value enables quantifying bubble stream gas flow rates using existing acoustic inversion methods. The presented method is based on averaging geo-referenced volume backscattering coefficients onto a high-resolution 3D voxel-grid. The results are multiplied with the voxel volume to represent measurements of the total backscattering cross-section within each voxel cell. Individual gridded values cannot be trusted because the beam pattern effects cause the values of individual targets to “smear” over multiple grid-cells. The true aggregated backscattering cross-section is thus estimated as the integral over the grid-cells affected by this smearing. Numerical simulation of MBES data acquisition over known targets assesses the method's validity and quantify it's uncertainty for different, realistic scenarios. The found low measurement bias (〈 1%), and dispersion (〈 5%) are promising for application in gas flow quantification methods.