Hydroacoustic Research Articles

Noise sensors - Part 1: Usage for the redesign of neighbourhoods

An interdisciplinary research project at the Stuttgart University of Applied Sciences with participation of the departments of urban planning, surveying and geoinformatics, business psychology and acoustics, present alternative urban scenarios. Different aspects for the redesign of the district, especially measures to reduce the individual motorised transport, were investigated. Virtual reality and augmented reality experiments are used in public spaces, to show redesign possibilities. Surveys of people, who had the experience of the future scenarios via the AR application and VR application, were carried out. Additionally, acoustical surveys and measurements in this research project were performed using so called noise sensors. These sensors based on MEMS microphones were developed for a simple and cost-effective measurement of external noise in a citizen sciences project. The measurement results of the noise sensors are saved online and can be accessed on the homepage "Sensor.Community". One of the noise sensors was used in the project area over a longer period. Measurement values were generated in third octave bands at intervals of 5 s and are available for further evaluations, e.g. for the calculation of indoor noise levels. The paper presents the research project and reports about the application of the noise sensor.

Room acoustic survey of open and enclosed learning spaces in Finland

Use of open learning spaces for education has increased in schools while enclosed learning spaces are still used in most schools in overall. There is some debate that these two learning space types might differ in acoustic perceptions. If they do, room acoustic quality might be one explaining factor. Therefore, our purpose was to determine, whether these learning space types differ in room acoustic properties. We studied 10 schools with mostly enclosed and another 10 schools with open or flexible learning spaces in Finland. Altogether 73 learning spaces were investigated among the 20 schools. Measurements concerned reverberation time, Speech Transmission Index, and background noise level. Most of enclosed learning spaces (63%) met the new requirements of Finnish building regulations for enclosed learning spaces. Opposite to that, most of open learning spaces (91%) did not meet the regulations for open learning spaces. The results can be used in the root cause analysis of a teacher survey conducted in the same 20 schools. Teacher survey will be reported in another study. Parallel analysis of objective and subjective results is necessary to judge whether the current regulations are justified.

Underwater sound produced by operational floating wind turbines

Deployment of offshore wind turbines is growing rapidly. Conventional "fixed bottom" wind turbines are limited to water depths of ~60 m, so deployment of offshore turbines to greater water depths requires floating systems. Both fixed bottom and floating turbines produce operational noise that impacts the marine environment. While the operational noises from fixed bottom turbines are relatively well understood, floating turbines are comparatively novel technologies and there is uncertainty as to the noise they will produce and its impact on the marine environment. Results from acoustic surveys of two floating wind farms off the Scottish Coast are presented; Kincardine wind farm which uses semi-submersible floating systems; and Hywind which uses a floating spar structure. The acoustic signature of each consists of continuous tonal noise related to drivetrain vibration and transient broadband events related to the mooring systems. A machine learning algorithm was used to discriminate between drivetrain noise and mooring noise. Sound propagation modelling was used to determine the source levels for turbines at each wind farm allowing a comparison of the noise signature of different floating structures.

Acoustic study of urban bat diversity in Veliko Tarnovo, Bulgaria

Urbanisation presents both challenges and opportunities for bat populations. Understanding the composition and distribution of urban bat species is the first step toward their effective conservation. However, comprehensive studies on urban bat diversity remain scarce in Bulgaria. To address this gap, we conducted a year-long acoustic survey on five different habitats in the urbanised territory of Veliko Tarnovo. Our study identified the presence of eight bat species and other distinct acoustic groups within the urban environment. Bats were found to utilise various habitats across the urban landscape. Although the activity was highest during May and September, bats remained active all year round, including the winter. Our study provides insights into synanthropic bats’ seasonal and microhabitat-dependent use of urban areas. These findings will establish a baseline for future research focusing on habitat preferences, population dynamics, and interactions with urbanisation. Such research is essential for the effective management and conservation of bat populations in urban environments.

Double tombolo formation by regressive barrier widening and landside submergence: The case of Orbetello, Italy

The double tombolo of Orbetello, in Italy, has formed during the Holocene around an ancient central tombolo. Earlier models consider that its sand barriers formed as sand spits that stretched from the mainland to a coastal island before enlarging seawards. This evolution, however, remains speculative. In order to test these models, we conducted the first study of a double tombolo that combines coring of its sand barriers and comprehensive imaging of its internal structure using sub-bottom acoustic surveys offshore and in the back-barrier. Sediment ages were constrained by 14C, luminescence, and U/Th dating. Acoustic images below the lagoon show that the barriers are in fact broad regressive strandplains that initiated on the flanks of the preexisting central isthmus when sea level was −7 ± 1 m lower than today. The strandplains then rose upwards and outwards, tracking sea level rise over the past 7 kyr. The oldest and lowest parts of the strandplains were flooded into the shallow intervening lagoon. The central isthmus is composed of regressive sand barriers accreted around a MIS 5.5 core during subsequent stages MIS 5.3 and MIS 5.1. The emplacement of the isthmus interrupted longshore drift between the mainland and the coastal island, converting the flanks of the initial tombolo into terminal sinks in which sand accretion accelerated, spurring early and rapid regression during the Holocene. A review of the environmental parameters conducive to double tombolo formation suggests that double tombolos may represent a frequent, albeit short-lived stage during the enlargement of single tombolos.

Fourier-transform near-infrared spectroscopy first application to age determination in European fish species: the case of the Atlantic horse mackerel from the central Mediterranean Sea

Fourier-transform near-infrared spectroscopy (FT-NIRS) has recently been used to reduce the processing time for estimates of annual age in several fish species. The present study represents the first application of this technique to marine organisms from the European waters. Atlantic horse mackerel (Trachurus trachurus) from the central Mediterranean Sea was selected for its ecological role, its socioeconomic value, and because its age is regularly estimated by otolith reading under a stereomicroscope for stock assessment purposes. Absorption spectra of the whole otoliths were acquired by FT-NIRS across a multiyear dataset, obtained during acoustic surveys carried out in different regions of the central Mediterranean Sea. The acquired spectra were processed to optimize calibration models to predict age. The best linear models obtained by the optimizing procedure predicted age successfully with a coefficient of determination of 0.95–0.96, mean squared error of 0.5 years, and bias <0.03 years. Although regional calibration models were also obtained separately, the models with all regions combined appeared much more robust and encompassing the whole age range, suggesting its usefulness for the entire spatial distribution range considered in the preset study. The use of FT-NIRS allowed to successfully predict age for horse mackerel between the ages 1 and 14, with similar or better precision (and bias) of the traditional methods, also increasing the standardization by repeatability of age assignments. The preliminary results obtained by this study encourage further effort to fine-tune the calibration model for new species, in order to apply a precise, rapid, cost-effective method for stock assessment purposes of Mediterranean fish resources.

A Three–Year Comparison of Fluctuations in the Occurrence of the Giant Jellyfish (Nemopilema nomurai)

In this study, acoustic, sighting, trawl, and marine environmental surveys were used to determine the vertical distribution and density of giant jellyfish that have been observed in Korean waters over the past 3 years. From 2020 to 2022, annual surveys were conducted in May and July in the East China Sea and waters adjacent to South Korea. The acoustic data were processed by identifying and eliminating all signals considered as noise while excluding those suspected to be jellyfish signals. Subsequently, a single target detection method was employed. Giant jellyfish are distributed mostly in the middle and low layers. In May 2021, the average population density of giant jellyfish was recorded as 11.6 (ind./ha), which was the highest density. In July 2022, this value decreased to 1.7 (ind./ha), marking the lowest density. The sighting survey, which allows for the identification of jellyfish distributed in the surface layer, exhibited a difference of approximately 0.13 times compared to the acoustic survey conducted in the middle and low layers in 2020. In 2021 and 2022, this difference was approximately 0.11 times and 0.24 times, respectively. The average of this difference was 0.16 times or greater.

Passive acoustic surveys demonstrate high densities of sperm whales off the mid-Atlantic coast of the USA in winter and spring

Oceans are increasingly crowded by anthropogenic activities yet the impact on Outer Continental Shelf (OCS) marine life remains largely unquantified. The MAPS (Marine Mammal Acoustic and Spatial Ecology) study of 2019 included passive acoustic and visual vessel surveys over the Mid-Atlantic OCS of the USA to address data gaps in winter/spring for deep-diving cetaceans, including sperm whales. Echolocation clicks were used to derive slant ranges to sperm whales for design- and model-based density estimates. Although more survey effort was realised in the spring, high densities of whales were identified in both winter and spring (10.46 and 8.89 per 1000 km2 respectively). The spring model-based abundance estimate of 1587 whales (CI 946–2663) was considered the most representative figure, in part due to lower coefficients of variation. Modelling suggested that high densities of whales were associated with warm core rings, eddies and edges. As OCS waters provide an important foraging habitat for North Atlantic sperm whales, appropriate mitigation is required to ensure commercial pressures to develop offshore energy do not negatively affect this endangered species.

Vertical distribution of Lake Superior cisco (Coregonus artedi) spawning aggregations and implications for population monitoring

Cisco (Coregonus artedi) support an evolving commercial roe fishery in Wisconsin waters of Lake Superior. To monitor trends in spawning cisco abundance, fishery managers recently began estimating adult biomass and exploitation using fall hydroacoustic surveys, which were combined with gill net surveys to inform apportionments of acoustic data. The gill net survey design consisted of paired top-suspended and bottom-set gill nets, but only the sex ratios from top nets are currently used with the hydroacoustic surveys due to an assumption that cisco in Lake Superior are pelagic spawners. However, the vertical sex distribution of cisco during spawning aggregations has been described as dynamic, with males becoming more bottom-oriented throughout the spawning season. We used multilevel aggregated binomial regressions to: 1) determine if there is bias between top and bottom gill net catches of cisco for either sex and if it changes throughout the spawning season, 2) evaluate how the vertical distribution of males and females may create bias in sex ratios used to estimate exploitation, and 3) explore the effect that maturity (i.e., gonadal development) has on vertical distribution during spawning aggregations. We identified sex-specific bias in vertical catch location that has the potential to bias estimates of sex ratio, and the source of this bias may be attributable to maturity driven changes in behavior. These findings highlight a need for caution when relying on gill nets to apportion cisco sex ratios during spawning aggregations and provide support for a non-pelagic alternative hypothesis of spawning behavior.

Anomalous ocean currents and European anchovy dispersal in the Iberian ecosystem

Unlike other upwelling areas where sardine and anchovy species dominate the pelagic ecosystems, the Western Iberian ecosystem has been consistently dominated by European sardine Engraulis encrasicolus, while anchovy had a residual presence from the start of acoustic surveys, in 1989, to 2014. Since 2015, the abundance of anchovy in the Western Iberian margins has sharply increased and continues to show an increasing trend as of 2023. It is unclear if this increase is a result of dispersal from nearby recruitment areas, higher survival rates of early life stages due to favorable environmental conditions, or both. We used a set of different models to simulate the dispersion and survival of anchovy early life stages in the Iberian region for the years preceding the increase in anchovy abundance. An ocean model simulation with the model CROCO provided the fields used as background for Lagrangian simulations coupled to an individual-based model of anchovy eggs and larvae. We simulated the years 2013-2015, and the results show that in 2014 and 2015, anomalous upper-ocean circulation patterns with strong and persistent westward currents transported a large number of eggs and larvae from the Bay of Biscay (BoB) westward along the Northern Iberian margin. The maximum transport occurred in June and July 2015, when 8 and 4%, respectively, of the eggs spawned in the BoB potentially reached the Iberian west coast as larvae. This process might explain the increase in anchovy abundance in the Western Iberian ecosystem. The results of the study show that episodes of anomalous intense ocean currents, when coincident with high presence of eggs, can lead to the colonization of new areas, and connectivity between areas varies dramatically with time.

Island area and diet predict diversity and distribution of bats in a Pacific Northwest archipelago

Abstract The island biogeography theory predicts that species richness in islands and island-like systems is the ultimate result of island isolation and area. Species with high dispersal capabilities are predicted to be less affected by these factors because of their capacity to move more efficiently between islands or habitats, and here we test this idea in bats, the only mammals capable of flight. We conducted mist net and acoustic surveys across 21 islands in the San Juan Archipelago (Washington State, United States) and adjacent northwest mainland to: (i) investigate the effects of island area, distance from mainland, and habitat on bat diversity; and (ii) evaluate whether differences in morphological (body mass, forearm length, wing loading) and ecological (dietary niche breadth, foraging guild) traits among species influence their prevalence across islands. We found that island size strongly influenced patterns of species richness, with larger islands having a greater number of bat species. However, neither island distance from mainland nor any measure of habitat availability was a significant predictor of species richness at the scale of this study. Additionally, we found that dietary niche breadth, as opposed to any morphological trait, best predicted the prevalence of species across the islands. Our results suggest that species with more specialized diets may be more vulnerable to habitat fragmentation, and provide insight into how geographic and ecological factors affect the diversity of insular bat communities, adding to growing knowledge about the role of species traits as mediators of their responses to large-scale landscape structure.

Vegetation canopy height shapes bats’ occupancy: a remote sensing approach

ABSTRACT Anthropogenic activities have significantly altered land cover on a global scale. These changes often have a negative effect on biodiversity limiting the distribution of species. The extent of the effect on species’ distribution depends on the landscape composition and configuration at a local and landscape level. To better understand this effect on a large scale, we evaluated how land cover and vegetation structure shape bat species’ occurrence while considering species’ imperfect detection. We hypothesize that intensification of anthropogenic activities in agriculture, for example, reduces heterogeneity of land cover and vegetation structure, and thereby, limits bat occurrence. To investigate this, we conducted acoustic bat sampling across 59 locations in southern Portugal, each with three spatial replicates. We derived fine-scale vegetation structural metrics by combining spaceborne LiDAR (GEDI) and synthetic aperture radar data (Sentinel-1 and ALOS/PALSAR-2). Additionally, we included land cover metrics and high-resolution climate data from CHELSA. Our findings revealed an important relationship between bat species’ occupancy and vegetation structure, particularly with vegetation canopy height. Moreover, forest and shrubland proportions were the main land cover types influencing bat species responses. All species’ best-ranking occupancy models included at least one climatic variable (temperature, humidity, or potential evapotranspiration), demonstrating the importance of climate when predicting bat distribution. Our acoustic surveys had a species’ detection probability varying from 0.19 to 0.86, and it was influenced by night conditions. These findings underscore the importance of modeling imperfect detection, especially for highly vagile and elusive organisms like bats. Our results demonstrate the effectiveness of using vegetation and landscape metrics derived from high-resolution remote sensing data to model species distribution in the context of biodiversity monitoring and conservation.

An Acoustical Survey of a Catholic Chapel

The modern-day celebration of the Roman Catholic Mass places an increased emphasis on the speech, music, and participation of the congregation in these ministries. Thus, there is an increasing emphasis on ensuring churches provide acceptable experiences for speech and music. This paper presents the results of an empirical study of a Roman Catholic chapel. The acoustic metrics were determined from impulse response data gathered via microphones recording balloon pops. Sound pressure level measurements were gathered with a sound pressure meter recording a boombox noise source. The balloon pops and boombox noise source were placed in various locations indicative of sources during the celebration of the Mass. Findings include that the RT30 measurements of the chapel were consistent with the EDT values suggesting that early reflections should not provide a significant impact on the perceived reverberation of the chapel. Based on the results, we conclude that the chapel does not have any significant deficiencies in its ability to transmit speech and music. This provides a basis for the chapel administration to invest in a professional assessment as they begin a project to install a new organ and choir loft.

Comparing the cost-effectiveness of drones, camera trapping and passive acoustic recorders in detecting changes in koala occupancy.

Quantifying the cost-effectiveness of alternative sampling methods is crucial for efficient biodiversity monitoring and detection of population trends. In this study, we compared the cost-effectiveness of three novel sampling methods for detecting changes in koala (Phascolarctos cinereus) occupancy: thermal drones, passive acoustic recorders and camera trapping. Specifically, we fitted single-season occupancy-detection models to data recorded from 46 sites in eight bioregions of New South Wales, Australia, between 2018 and 2022. We explored the effect of weather variables on daily detection probability for each method and, using these estimates, calculated the statistical power to detect 30%, 50% and 80% declines in koala occupancy. We calculated power for different combinations of sites (1-200) and repeat surveys (2-40) and developed a cost model that found the cheapest survey design that achieved 80% power to detect change. On average, detectability of koalas was highest with one 24-h period of acoustic surveys (0.32, 95% CI's: 0.26, 0.39) compared to a 25-ha flight of drone surveys (0.28, 95% 0.15, 0.48) or a 24-h period of camera trapping consisting of six cameras (0.019, 95% CI's: 0.014, 0.025). We found a negative quadratic relationship between detection probability and air temperature for all three methods. Our power and cost analysis suggested that 148 sites surveyed with acoustic recorders deployed for 14 days would be the cheapest method to sufficiently detect a 30% decline in occupancy with 80% power. We recommend passive acoustic recorders as the most efficient sampling method for monitoring koala occupancy compared to cameras or drones. Further comparative studies are needed to compare the relative effectiveness of these methods and others when the monitoring objective is to detect change in koala abundance over time.

Towards automated animal density estimation with acoustic spatial capture-recapture.

Passive acoustic monitoring can be an effective way of monitoring wildlife populations that are acoustically active but difficult to survey visually, but identifying target species calls in recordings is non-trivial. Machine learning (ML) techniques can do detection quickly but may miss calls and produce false positives, i.e., misidentify calls from other sources as being from the target species. While abundance estimation methods can address the former issue effectively, methods to deal with false positives are under-investigated. We propose an acoustic spatial capture-recapture (ASCR) method that deals with false positives by treating species identity as a latent variable. Individual-level outputs from ML techniques are treated as random variables whose distributions depend on the latent identity. This gives rise to a mixture model likelihood that we maximize to estimate call density. We compare our method to existing methods by applying it to an ASCR survey of frogs and simulated acoustic surveys of gibbons based on real gibbon acoustic data. Estimates from our method are closer to ASCR applied to the dataset without false positives than those from a widely used false positive "correction factor" method. Simulations show our method to have bias close to zero and accurate coverage probabilities and to perform substantially better than ASCR without accounting for false positives.

Mapping benthic sediment types and composition in a turbid Jamaican bay using hydroacoustic data and different spatially explicit interpolation techniques

Hydroacoustic sensors are better suited for mapping shallow turbid coastal waters, but surveys must take place along transects, necessitating the interpolation of processed, large, spatially dependent hydroacoustic datasets to produce continuous maps. The aim of this study was to evaluate the performance of interpolation techniques that can handle data obtained from a hydroacoustic survey (covering 23 km2) conducted in 2016 in a shallow bay in southwestern Jamaica. Independent data or a 5-fold cross validation were used to compare the accuracies of continuous maps of benthic sediment percentage composition and submerged aquatic vegetation cover (SAV), interpolated from the survey data using random forest (RF) regression (RFr), RFr kriging (RFrk) and Ensemble Machine Learning (EML) with rotated/oblique coordinates. Categorical/thematic maps of sediment types interpolated using RF (RFc) and EML classifiers and a spatial Markov chain model with indicator kriging predictor/simulation (Mc-IK) were also compared. RFr, RFc, RFrk and EML covariates included coordinates and auxiliary data. The following maps were more accurate than corresponding interpolated maps: EML - mud, sand, silt, clay and SAV; RFr - loss on ignition [LOI]; RFrk - calcium carbonate; Mc-IK – sediment types with 3 and 4 classes. The Mc-IK maps were also the most suitable maps for assessing spatial and quantitative changes. The predominant substrate type in the bay was muddy sand covering 15.7 km2. The highest percentage composition of mud, silt, clay, and organic matter [LOI] was found close to the mouth of the Black River where sandy mud was found, indicating a high input of terrigenous sediment. Calcium carbonate composition showed an opposite spatial trend and was positively associated with higher SAV percentage cover. The bay is now part of a protected area, as such, these maps will be used as a baseline to gauge the impact of increased protection/management on benthic dynamics.

The energy content and demographic composition of Antarctic krill (Euphausia superba) swarms in East Antarctica

Antarctic krill (Euphausia superba) aggregations are important in the Southern Ocean ecosystem; however, the amount of energy that krill swarms contain is unknown. We assessed the population structure of E. superba swarms in East Antarctica for two different sectors corresponding to the study areas of the ENRICH (2019) and TEMPO (2021) voyages. We also estimated swarm energy content based on calorific values for whole animals and key tissues (gonads, digestive gland, and muscle). A relationship between the energy content of muscle and gonads in gravid and adult females was detected: The higher the energy content in muscle tissue, the lower the energy content in the gonad tissue. This relationship was consistent for different female maturity stages, suggesting a trade-off in energy allocation between somatic growth and reproduction. Females will redirect energy to gonad development at the expense of somatic growth. Total swarm energy was calculated using the volumetric density of each swarm estimated from an acoustic survey. Four swarm types (A-D) were determined based on the body length and maturity stage of the individuals. The highest energy contents were found in Type B swarms, composed predominantly of adult males and females, and Type D swarms, comprising mostly gravid females. Trends in swarm demographic composition and energy content were consistent between surveys. We recommend swarm energy be incorporated into ecosystem models to represent energy transfer to top predators.

Natural disturbance allows multiple anuran taxa to persist in a dynamic wetland complex

AbstractThe maintenance of biological diversity is frequently enhanced in a heterogenous landscape by some level of disturbance. Thus, when a landscape becomes stabilized and homogenized through the spread of an invasive plant species, there may be consequences for native biodiversity, particularly for those biotas that depend on the pre‐existing, natural disturbance regime of that landscape. At Long Point, Ontario, Canada, a sandspit in Lake Erie, the community of anuran amphibians experienced the co‐incidence of 2 major disturbance events, 1 anthropogenic and 1 natural: an intervention to remove the invasive form of the common reed (Phragmites australis australis) and a spate of extensive dune washouts caused by high water levels and storms. As a result of the unequal distribution of disturbance in the landscape, different areas were affected by only the dune washouts, only the common reed treatment, both disturbances, or neither disturbance. From May to August 2021, we explored how the disturbances affected the resident frog and toad species habitat associations in the 4 different disturbance conditions in the landscape. By means of minnow traps, acoustic surveys, and visual surveys, we determined that the dune washouts had a greater influence on resident anurans than the common reed treatment. Ranid frogs (Ranidae) tended to inhabit non‐washout sites, whereas Fowler's toads (Anaxyrus fowleri) congregated in the newly formed, sand flats and shallow, de‐vegetated pools resulting from the washouts. Neither ranid frogs nor Fowler's toads demonstrated avoidance of the sites affected by the herbicide treatment and mechanical rollover of the common reed. This evidence of species sorting, which enables multiple species to persist in a dynamic and heterogeneous landscape, suggests that wildlife management in aid of threatened species recovery may find success by encouraging natural disturbance regimes in dynamic landscapes. Future management endeavors should consider refraining from complete dune stabilization in these landscapes to allow the natural washouts to occur at a higher frequency, while ensuring the continued effectiveness of the invasive common reed eradication.

Automated echolocation classifiers vary in accuracy for northeastern U.S. bat species.

Acoustic surveys of bat echolocation calls are an important management tool for determining presence and probable absence of threatened and endangered bat species. In the northeastern United States, software programs such as Bat Call Identification (BCID), Kaleidoscope Pro (KPro), and Sonobat can automatically classify ultrasonic detector sound files, yet the programs' accuracy in correctly classifying calls to species has not been independently assessed. We used 1,500 full-spectrum reference calls with known identities for nine northeastern United States bat species to test the accuracy of these programs using calculations of Positive Predictive Value (PPV), Negative Predictive Value (NPV), Sensitivity (SN), Specificity (SP), Overall Accuracy, and No Information Rate. We found that BCID performed less accurately than other programs, likely because it only operates on zero-crossing data and may be less accurate for recordings converted from full-spectrum to zero-crossing. NPV and SP values were high across all species categories for SonoBat and KPro, indicating these programs' success at avoiding false positives. However, PPV and SN values were relatively low, particularly for individual Myotis species, indicating these programs are prone to false negatives. SonoBat and KPro performed better when distinguishing Myotis species from non-Myotis species. We expect less accuracy from these programs for acoustic recordings collected under normal working conditions, and caution that a bat acoustic expert should verify automatically classified files when making species-specific regulatory or conservation decisions.

Pockmarks and associated fresh submarine groundwater discharge in the seafloor of Puck Bay, southern Baltic Sea

This is the first well-documented report on the occurrence of pockmarks in Puck Bay. Pockmarks in the seafloor of Puck Bay were discovered during a hydroacoustic survey carried out in 2020. They are located at a depth of 25–27 m in the southwestern part of the bay. Significant depletion of chloride (Cl−) concentrations in sediment pore water was found within the depressions. Most likely, the formation of pockmarks was due to groundwater flow through the Miocene–Pleistocene system of aquifers, which extends from land to the bay area. One-dimensional modeling of vertical Cl− concentration profiles in pore water revealed the upward flow of freshened groundwater within the pockmarks. The magnitude of submarine groundwater discharge (SGD) was estimated to vary from 1.53·10−2 to 18·10−2 L·m−2·h−1. The effect of groundwater seepage was also observed at 3 cm above the seafloor within the pockmarks, which was identified as a decrease in salinity of approximately 0.12 PSU compared to reference sites.Furthermore, due to the effect of water advection, SGD can be detected even several meters above the seafloor as a decrease in salinity values within the thermocline layer.